Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // ------------------------------------------------------------------- 27 // GEANT4 Class file 28 // 29 // 30 // File name: G4LightIonQMDReaction.hh 31 // 32 // Author: Koi, Tatsumi (tkoi@slac.stanford.edu) 33 // 34 // Creation date: 02 April 2007 35 // ----------------------------------------------------------------------------- 36 // 37 // 081107 Add UnUseGEM (then use the default channel of G4Evaporation) 38 // UseFrag (chage criterion of a inelastic reaction) 39 // 090331 Change member shenXS and genspaXS object to pointer 40 // 41 42 #ifndef G4LightIonQMDReaction_hh 43 #define G4LightIonQMDReaction_hh 44 45 #include "G4QMDSystem.hh" 46 #include "G4LightIonQMDCollision.hh" 47 #include "G4LightIonQMDMeanField.hh" 48 #include "G4QMDParticipant.hh" 49 50 #include "G4HadronicInteraction.hh" 51 52 #include "G4Evaporation.hh" 53 #include "G4ExcitationHandler.hh" 54 55 class G4VCrossSectionDataSet; 56 class G4BGGPionElasticXS; 57 class G4BGGPionInelasticXS; 58 59 class G4LightIonQMDReaction : public G4HadronicInteraction 60 { 61 public: 62 G4LightIonQMDReaction(); 63 ~G4LightIonQMDReaction() override; 64 65 std::vector< G4QMDSystem* > GetFinalStates(); 66 67 G4HadFinalState* ApplyYourself(const G4HadProjectile &aTrack, 68 G4Nucleus & targetNucleus) override; 69 70 G4ExcitationHandler* GetExcitationHandler() {return excitationHandler;}; 71 72 void UnUseGEM() {gem = false; setEvaporationCh();}; 73 void UseFRAG() {frag = true;}; 74 75 void SetTMAX(G4int i) { maxTime = i; }; 76 void SetDT(G4double t) { deltaT = t; }; 77 void SetEF(G4double x) { envelopF = x; }; 78 79 void ModelDescription(std::ostream& outFile) const override; 80 81 //copy is unexpected 82 G4LightIonQMDReaction(const G4LightIonQMDReaction& right) = delete; 83 const G4LightIonQMDReaction& operator = (const G4LightIonQMDReaction& right) = delete; 84 85 private: 86 87 void setEvaporationCh(); 88 void setHighEnergyModel(); 89 90 G4LightIonQMDMeanField* meanField; 91 92 G4LightIonQMDCollision* collision; 93 94 void doCollision(); 95 std::vector< G4QMDSystem* > doClusterJudgment(); 96 97 G4QMDSystem* system; 98 G4double deltaT; 99 G4int maxTime; 100 G4double envelopF; 101 102 G4ExcitationHandler* excitationHandler; 103 104 // b pd_proj pd_targ z_p a_p z_t a_t plab elab 105 // G4double offSetOfCollision( G4double , G4ParticleDefinition* , G4ParticleDefinition* , G4int , G4int , G4int , G4int , G4double , G4double ); 106 // b pd_proj pd_targ plab elab bmax boostToCM 107 108 void calcOffSetOfCollision(G4double, const G4ParticleDefinition*, 109 const G4ParticleDefinition*, G4double, 110 G4double, G4double, G4ThreeVector); 111 112 G4double coulomb_collision_gamma_proj; 113 G4double coulomb_collision_rx_proj; 114 G4double coulomb_collision_rz_proj; 115 G4double coulomb_collision_px_proj; 116 G4double coulomb_collision_pz_proj; 117 118 G4double coulomb_collision_gamma_targ; 119 G4double coulomb_collision_rx_targ; 120 G4double coulomb_collision_rz_targ; 121 G4double coulomb_collision_px_targ; 122 G4double coulomb_collision_pz_targ; 123 124 G4VCrossSectionDataSet* theXS; 125 126 G4BGGPionElasticXS* pipElNucXS; 127 G4BGGPionElasticXS* pimElNucXS; 128 G4BGGPionInelasticXS* pipInelNucXS; 129 G4BGGPionInelasticXS* pimInelNucXS; 130 131 G4bool gem; 132 G4bool frag; 133 G4int secID; // Creator model ID for the secondaries 134 }; 135 136 #endif 137