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