include/Geant4/G4VPreCompoundFragment.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 // J. M. Quesada (August 2008).
0028 // Based  on previous work by V. Lara
0029 //
0030 // Modif (03 September 2008) by J. M. Quesada for external choice of inverse
0031 // cross section option
0032 // JMQ (06 September 2008) Also external choice has been added for:
0033 //                      - superimposed Coulomb barrier (if useSICB=true)
0034 // 20.08.2010 V.Ivanchenko added int Z and A and cleanup; added
0035 //                        G4ParticleDefinition to constructor,
0036 //                        inline method to build G4ReactionProduct;
0037 //                        remove string name
0038 //
0039
0040 #ifndef G4VPreCompoundFragment_h
0041 #define G4VPreCompoundFragment_h 1
0042
0043 #include "G4ios.hh"
0044 #include <iomanip>
0045 #include "G4ParticleDefinition.hh"
0046 #include "G4IonTable.hh"
0047 #include "G4Fragment.hh"
0048 #include "G4ReactionProduct.hh"
0049 #include "G4Pow.hh"
0050 #include "G4VSIntegration.hh"
0051
0052 class G4NuclearLevelData;
0053 class G4DeexPrecoParameters;
0054 class G4VCoulombBarrier;
0055 class G4InterfaceToXS;
0056
0057 class G4VPreCompoundFragment : G4VSIntegration
0058 {
0059 public:
0060
0061   explicit G4VPreCompoundFragment(const G4ParticleDefinition*,
0062                   G4VCoulombBarrier*);
0063
0064   ~G4VPreCompoundFragment() override;
0065
0066   // Run time initialization method
0067   G4bool Initialize(const G4Fragment& aFragment);
0068
0069   // Calculates the total (integrated over kinetic energy) emission
0070   // probability of a fragment
0071   virtual G4double CalcEmissionProbability(const G4Fragment&);
0072
0073   // sample kinetic energy of emitted fragment
0074   virtual G4double SampleKineticEnergy(const G4Fragment&);
0075
0076   G4double ProbabilityDensityFunction(G4double energy) override;
0077
0078   inline G4ReactionProduct* GetReactionProduct() const;
0079
0080   G4int GetA() const { return theA; }
0081
0082   G4int GetZ() const { return theZ; }
0083
0084   G4int GetRestA() const { return theResA; }
0085
0086   G4int GetRestZ() const { return theResZ; }
0087
0088   G4double GetBindingEnergy() const { return theBindingEnergy; }
0089
0090   G4double GetEnergyThreshold() const
0091   {
0092     return theMaxKinEnergy - theCoulombBarrier;
0093   }
0094
0095   G4double GetEmissionProbability() const { return theEmissionProbability; }
0096
0097   G4double GetNuclearMass() const { return theMass; }
0098
0099   G4double GetRestNuclearMass() const { return theResMass; }
0100
0101   const G4LorentzVector& GetMomentum() const { return theMomentum; }
0102
0103   void SetMomentum(const G4LorentzVector& lv) { theMomentum = lv; }
0104
0105   //for inverse cross section choice
0106   void SetOPTxs(G4int) {}
0107   //for superimposed Coulomb Barrier for inverse cross sections
0108   void UseSICB(G4bool use) { useSICB = use; }
0109
0110   friend std::ostream&
0111   operator<<(std::ostream&, const G4VPreCompoundFragment*);
0112   friend std::ostream&
0113   operator<<(std::ostream&, const G4VPreCompoundFragment&);
0114
0115   G4VPreCompoundFragment(const G4VPreCompoundFragment &right) = delete;
0116   const G4VPreCompoundFragment&
0117   operator= (const G4VPreCompoundFragment &right) = delete;
0118   G4bool operator==(const G4VPreCompoundFragment &right) const = delete;
0119   G4bool operator!=(const G4VPreCompoundFragment &right) const = delete;
0120
0121 protected:
0122
0123   virtual G4double
0124   ProbabilityDistributionFunction(G4double, const G4Fragment&)
0125   { return 0.0; };
0126
0127   virtual G4double GetAlpha() const = 0;
0128
0129   virtual G4double GetBeta() const { return -theCoulombBarrier; }
0130
0131   G4NuclearLevelData* fNucData;
0132   G4DeexPrecoParameters* theParameters;
0133   G4Pow* g4calc;
0134   G4InterfaceToXS* fXSection{nullptr};
0135   const G4Fragment* pFragment{nullptr};
0136
0137   G4int theA;
0138   G4int theZ;
0139   G4int theResA{0};
0140   G4int theResZ{0};
0141   G4int theFragA{0};
0142   G4int theFragZ{0};
0143   //for inverse cross section choice
0144   G4int OPTxs;
0145   G4int index{0};
0146
0147   G4double theResA13{0.0};
0148   G4double theBindingEnergy{0.0};
0149   G4double theMinKinEnergy{0.0};
0150   G4double theMaxKinEnergy{0.0};
0151   G4double theResMass{0.0};
0152   G4double theReducedMass{0.0};
0153   G4double theMass;
0154
0155   G4double theEmissionProbability{0.0};
0156   G4double theCoulombBarrier{0.0};
0157
0158   //for superimposed Coulomb Barrier for inverse cross sections
0159   G4bool useSICB{true};
0160
0161 private:
0162
0163   const G4ParticleDefinition* particle;
0164   G4VCoulombBarrier* theCoulombBarrierPtr;
0165   G4LorentzVector theMomentum{0., 0., 0., 0.};
0166 };
0167
0168 inline G4ReactionProduct* G4VPreCompoundFragment::GetReactionProduct() const
0169 {
0170   G4ReactionProduct* theReactionProduct = new G4ReactionProduct(particle);
0171   theReactionProduct->SetMomentum(GetMomentum().vect());
0172   theReactionProduct->SetTotalEnergy(GetMomentum().e());
0173   return theReactionProduct;
0174 }
0175
0176 #endif