Changeset 1340 for trunk/source/processes/electromagnetic/highenergy/src
- Timestamp:
- Nov 5, 2010, 3:45:55 PM (14 years ago)
- Location:
- trunk/source/processes/electromagnetic/highenergy/src
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/highenergy/src/G4GammaConversionToMuons.cc
r1337 r1340 25 25 // 26 26 // 27 // $Id: G4GammaConversionToMuons.cc,v 1. 7 2008/10/16 14:29:48vnivanch Exp $28 // GEANT4 tag $Name: geant4-09-04-beta-01$27 // $Id: G4GammaConversionToMuons.cc,v 1.8 2010/10/26 14:15:40 vnivanch Exp $ 28 // GEANT4 tag $Name: emhighenergy-V09-03-02 $ 29 29 // 30 30 // ------------ G4GammaConversionToMuons physics process ------ … … 53 53 { 54 54 SetProcessSubType(15); 55 MeanFreePath = DBL_MAX; 55 56 } 56 57 -
trunk/source/processes/electromagnetic/highenergy/src/G4eeToHadronsModel.cc
r1337 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4eeToHadronsModel.cc,v 1. 9 2008/07/10 18:06:39vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: G4eeToHadronsModel.cc,v 1.10 2010/10/26 14:15:40 vnivanch Exp $ 27 // GEANT4 tag $Name: emhighenergy-V09-03-02 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 76 76 { 77 77 theGamma = G4Gamma::Gamma(); 78 highKinEnergy = HighEnergyLimit(); 79 lowKinEnergy = LowEnergyLimit(); 80 emin = lowKinEnergy; 81 emax = highKinEnergy; 82 peakKinEnergy = highKinEnergy; 83 epeak = emax; 78 84 } 79 85 … … 92 98 const G4DataVector&) 93 99 { 94 if(isInitialised) return;100 if(isInitialised) { return; } 95 101 isInitialised = true; 96 102 -
trunk/source/processes/electromagnetic/highenergy/src/G4eeToHadronsMultiModel.cc
r1337 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4eeToHadronsMultiModel.cc,v 1. 8 2009/04/12 17:48:21vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: G4eeToHadronsMultiModel.cc,v 1.9 2010/10/26 14:15:40 vnivanch Exp $ 27 // GEANT4 tag $Name: emhighenergy-V09-03-02 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 68 68 verbose(ver), 69 69 isInitialised(false) 70 {} 70 { 71 thKineticEnergy = DBL_MAX; 72 maxKineticEnergy = 1.2*GeV; 73 } 71 74 72 75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... … … 90 93 if(!isInitialised) { 91 94 isInitialised = true; 92 93 thKineticEnergy = DBL_MAX;94 maxKineticEnergy = 1.2*GeV;95 95 96 96 cross = new G4eeCrossSections(); -
trunk/source/processes/electromagnetic/highenergy/src/G4hhIonisation.cc
r1337 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4hhIonisation.cc,v 1.1 0 2010/06/04 10:23:31vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: G4hhIonisation.cc,v 1.11 2010/10/26 14:15:40 vnivanch Exp $ 27 // GEANT4 tag $Name: emhighenergy-V09-03-02 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 102 102 const G4ParticleDefinition* bpart) 103 103 { 104 if(isInitialised) return;104 if(isInitialised) { return; } 105 105 106 106 theParticle = part; 107 if(bpart) G4cout << "G4hhIonisation::InitialiseEnergyLossProcess WARNING: no " 108 << "base particle should be defined for the process " 109 << GetProcessName() << G4endl; 110 107 if(bpart) { 108 G4cout << "G4hhIonisation::InitialiseEnergyLossProcess WARNING: no " 109 << "base particle should be defined for the process " 110 << GetProcessName() << G4endl; 111 } 111 112 SetBaseParticle(0); 112 113 SetSecondaryParticle(G4Electron::Electron()); … … 114 115 mass = theParticle->GetPDGMass(); 115 116 ratio = electron_mass_c2/mass; 116 eth = 2.0*MeV*mass/proton_mass_c2;117 G4double eth = 2*MeV*mass/proton_mass_c2; 117 118 flucModel = new G4BohrFluctuations(); 118 119 119 120 G4int nm = 1; 120 121 121 minKinEnergy = MinKinEnergy();122 G4double minKinEnergy = MinKinEnergy(); 122 123 123 124 if(eth > minKinEnergy) { -
trunk/source/processes/electromagnetic/highenergy/src/G4mplIonisation.cc
r1337 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4mplIonisation.cc,v 1.1 0 2010/03/28 16:45:38vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: G4mplIonisation.cc,v 1.11 2010/10/26 15:40:03 vnivanch Exp $ 27 // GEANT4 tag $Name: emhighenergy-V09-03-02 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 49 49 #include "G4Electron.hh" 50 50 #include "G4mplIonisationModel.hh" 51 #include "G4mplIonisationWithDeltaModel.hh" 51 52 #include "G4BohrFluctuations.hh" 52 53 … … 85 86 const G4ParticleDefinition*) 86 87 { 87 if(isInitialised) return;88 if(isInitialised) { return; } 88 89 89 90 SetBaseParticle(0); 90 91 SetSecondaryParticle(G4Electron::Electron()); 91 92 92 G4mplIonisationModel* ion = new G4mplIonisationModel(magneticCharge,"PAI"); 93 G4mplIonisationWithDeltaModel* ion = 94 new G4mplIonisationWithDeltaModel(magneticCharge,"PAI"); 93 95 ion->SetLowEnergyLimit(MinKinEnergy()); 94 96 ion->SetHighEnergyLimit(MaxKinEnergy()); -
trunk/source/processes/electromagnetic/highenergy/src/G4mplIonisationModel.cc
r1337 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4mplIonisationModel.cc,v 1. 7 2009/04/12 17:35:41vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: G4mplIonisationModel.cc,v 1.8 2010/10/26 15:40:03 vnivanch Exp $ 27 // GEANT4 tag $Name: emhighenergy-V09-03-02 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 74 74 { 75 75 nmpl = G4int(abs(magCharge) * 2 * fine_structure_const + 0.5); 76 if(nmpl > 6) nmpl = 6;77 else if(nmpl < 1) nmpl = 1;76 if(nmpl > 6) { nmpl = 6; } 77 else if(nmpl < 1) { nmpl = 1; } 78 78 pi_hbarc2_over_mc2 = pi * hbarc * hbarc / electron_mass_c2; 79 79 chargeSquare = magCharge * magCharge; 80 80 dedxlim = 45.*nmpl*nmpl*GeV*cm2/g; 81 81 fParticleChange = 0; 82 mass = 0.0; 82 83 } 83 84 … … 94 95 monopole = p; 95 96 mass = monopole->GetPDGMass(); 96 if(!fParticleChange) fParticleChange = GetParticleChangeForLoss();97 if(!fParticleChange) { fParticleChange = GetParticleChangeForLoss(); } 97 98 } 98 99
Note: See TracChangeset
for help on using the changeset viewer.