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 "G4XMesonBaryonElastic.hh" 31 #include "G4KineticTrack.hh" 32 #include "G4Gamma.hh" 33 #include "G4PionPlus.hh" 34 #include "G4Proton.hh" 35 #include "G4XAqmElastic.hh" 36 #include "G4XPDGElastic.hh" 37 38 G4XMesonBaryonElastic::G4XMesonBaryonElastic() 39 { 40 // As a first approximation the model is ass 41 // the entire energy range 42 lowLimit = 0.; 43 highLimit = DBL_MAX; 44 } 45 46 47 G4XMesonBaryonElastic::~G4XMesonBaryonElastic( 48 { } 49 50 51 G4bool G4XMesonBaryonElastic::operator==(const 52 { 53 return (this == (G4XMesonBaryonElastic *) &r 54 } 55 56 57 G4bool G4XMesonBaryonElastic::operator!=(const 58 { 59 return (this != (G4XMesonBaryonElastic *) &r 60 } 61 62 63 G4double G4XMesonBaryonElastic::CrossSection(c 64 { 65 G4double sigma; 66 67 // No gamma-baryon elastic scattering 68 const G4ParticleDefinition* defLight = FindL 69 if (defLight == G4Gamma::GammaDefinition()) 70 { 71 sigma = 0; 72 } 73 else 74 { 75 76 G4LorentzVector p41 = trk1.Get4Momentum( 77 G4ThreeVector p3 = trk1.GetPosition(); 78 G4ParticleDefinition* def = G4PionPlus:: 79 80 G4KineticTrack piTrk(def, 81 trk1.GetFormationTime(), 82 p3, 83 (G4LorentzVector&)p41); 84 85 G4LorentzVector p42 = trk2.Get4Momentum( 86 G4KineticTrack pTrk(((G4ParticleDefiniti 87 trk2.GetFormationTime(), 88 (G4ThreeVector)trk2.GetPosition(), 89 (G4LorentzVector&)p42); 90 91 G4XAqmElastic aqm; 92 G4double xAqmDummy = aqm.CrossSection(pi 93 G4double xAqm = aqm.CrossSection(trk1,tr 94 G4double factor = 1.; 95 if (xAqmDummy != 0.0) 96 { 97 factor = xAqm / xAqmDummy; 98 } 99 G4XPDGElastic pdg; 100 101 sigma = pdg.CrossSection(piTrk,pTrk); 102 sigma = sigma * factor; 103 } 104 105 return sigma; 106 } 107 108 109 G4String G4XMesonBaryonElastic::Name() const 110 { 111 G4String name("MesonBaryonElasticCrossSectio 112 return name; 113 } 114 115 116 117 G4bool G4XMesonBaryonElastic::IsValid(G4double 118 { 119 G4bool answer = InLimits(e,lowLimit,highLimi 120 121 return answer; 122 } 123