Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 28 #include "globals.hh" 29 #include "G4ios.hh" 30 #include "G4KineticTrack.hh" 31 #include "G4XResonance.hh" 32 #include "Randomize.hh" 33 #include "G4Proton.hh" 34 35 G4XResonance::G4XResonance(const G4ParticleDef 36 const G4ParticleDefinition* in2, 37 G4int iIsospinOut1, G4double spinOut1 38 G4int iIsospinOut2, G4double spinOut2 39 G4String subType1, G4String subType2, 40 const G4VXResonanceTable& sigmaTable) 41 isoOut1(iIsospinOut1), iSpinOut1(spinOut1), 42 isoOut2(iIsospinOut2), iSpinOut2(spinOut2), 43 44 { 45 table = sigmaTable.CrossSectionTable(); 46 // Check if there is a valid cross section t 47 if (table == 0) 48 throw G4HadronicException(__FILE__, __LINE 49 50 name = in1->GetParticleName() + in2->GetPart 51 } 52 53 54 G4XResonance::~G4XResonance() 55 { 56 delete table; 57 table = 0; 58 } 59 60 61 G4bool G4XResonance::operator==(const G4XReson 62 { 63 return (this == (G4XResonance *) &right); 64 } 65 66 67 G4bool G4XResonance::operator!=(const G4XReson 68 { 69 return (this != (G4XResonance *) &right); 70 } 71 72 73 G4String G4XResonance::Name() const 74 { 75 return name; 76 } 77 78 G4double G4XResonance::CrossSection(const G4Ki 79 { 80 G4bool dummy = false; 81 G4double sqrtS = (trk1.Get4Momentum() + trk2 82 83 // pp -> trk1 + trk2 cross section 84 G4double sigma = table->GetValue(sqrtS,dummy 85 86 // Isospin correction 87 // from G4VXResonance; 88 // G4cout << "SigmaCheck "<<sigma; 89 sigma *= IsospinCorrection(trk1,trk2,isoOut1 90 // G4cout << " "<<sigma<<G4endl; 91 92 // Detailed balance 93 if (trk1.GetDefinition()->IsShortLived() || 94 { 95 sigma *= DetailedBalance(trk1,trk2, isoOut1 96 } 97 98 return sigma; 99 100 } 101 102