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
0051 class G4NuclearLevelData;
0052 class G4DeexPrecoParameters;
0053 class G4VCoulombBarrier;
0054 class G4InterfaceToXS;
0055
0056 class G4VPreCompoundFragment
0057 {
0058 public:
0059
0060   explicit G4VPreCompoundFragment(const G4ParticleDefinition*,
0061                   G4VCoulombBarrier*);
0062
0063   virtual ~G4VPreCompoundFragment();
0064
0065   friend std::ostream&
0066   operator<<(std::ostream&, const G4VPreCompoundFragment*);
0067   friend std::ostream&
0068   operator<<(std::ostream&, const G4VPreCompoundFragment&);
0069
0070   // =====================
0071   // Pure Virtual methods
0072   // =====================
0073
0074   // Run time initialization method
0075   G4bool Initialize(const G4Fragment& aFragment);
0076
0077   // Methods for calculating the emission probability
0078   // ------------------------------------------------
0079
0080   // Calculates the total (integrated over kinetic energy) emission
0081   // probability of a fragment
0082   virtual G4double CalcEmissionProbability(const G4Fragment&) = 0;
0083
0084   // sample kinetic energy of emitted fragment
0085   virtual G4double SampleKineticEnergy(const G4Fragment&) = 0;
0086
0087   inline G4ReactionProduct* GetReactionProduct() const;
0088
0089   G4int GetA() const { return theA; }
0090
0091   G4int GetZ() const { return theZ; }
0092
0093   G4int GetRestA() const { return theResA; }
0094
0095   G4int GetRestZ() const { return theResZ; }
0096
0097   G4double GetBindingEnergy() const { return theBindingEnergy; }
0098
0099   G4double GetEnergyThreshold() const
0100   {
0101     return theMaxKinEnergy - theCoulombBarrier;
0102   }
0103
0104   G4double GetEmissionProbability() const { return theEmissionProbability; }
0105
0106   G4double GetNuclearMass() const { return theMass; }
0107
0108   G4double GetRestNuclearMass() const { return theResMass; }
0109
0110   const G4LorentzVector& GetMomentum() const { return theMomentum; }
0111
0112   void SetMomentum(const G4LorentzVector& lv) { theMomentum = lv; }
0113
0114   //for inverse cross section choice
0115   void SetOPTxs(G4int) {}
0116   //for superimposed Coulomb Barrier for inverse cross sections
0117   void UseSICB(G4bool use) { useSICB = use; }
0118
0119   G4VPreCompoundFragment(const G4VPreCompoundFragment &right) = delete;
0120   const G4VPreCompoundFragment&
0121   operator= (const G4VPreCompoundFragment &right) = delete;
0122   G4bool operator==(const G4VPreCompoundFragment &right) const = delete;
0123   G4bool operator!=(const G4VPreCompoundFragment &right) const = delete;
0124
0125 protected:
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
0136   G4int theA;
0137   G4int theZ;
0138   G4int theResA{0};
0139   G4int theResZ{0};
0140   G4int theFragA{0};
0141   G4int theFragZ{0};
0142   //for inverse cross section choice
0143   G4int OPTxs;
0144   G4int index{0};
0145
0146   G4double theResA13{0.0};
0147   G4double theBindingEnergy{0.0};
0148   G4double theMinKinEnergy{0.0};
0149   G4double theMaxKinEnergy{0.0};
0150   G4double theResMass{0.0};
0151   G4double theReducedMass{0.0};
0152   G4double theMass;
0153
0154   G4double theEmissionProbability{0.0};
0155   G4double theCoulombBarrier{0.0};
0156
0157   //for superimposed Coulomb Barrier for inverse cross sections
0158   G4bool useSICB{true};
0159
0160 private:
0161
0162   const G4ParticleDefinition* particle;
0163   G4VCoulombBarrier* theCoulombBarrierPtr;
0164   G4LorentzVector theMomentum{0., 0., 0., 0.};
0165 };
0166
0167 inline G4ReactionProduct* G4VPreCompoundFragment::GetReactionProduct() const
0168 {
0169   G4ReactionProduct* theReactionProduct = new G4ReactionProduct(particle);
0170   theReactionProduct->SetMomentum(GetMomentum().vect());
0171   theReactionProduct->SetTotalEnergy(GetMomentum().e());
0172   return theReactionProduct;
0173 }
0174
0175 #endif