- Timestamp:
- Nov 5, 2010, 3:45:55 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/standard/src/G4ionIonisation.cc
r1337 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4ionIonisation.cc,v 1.7 0 2009/11/27 20:06:32vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: G4ionIonisation.cc,v 1.72 2010/10/26 10:42:04 vnivanch Exp $ 27 // GEANT4 tag $Name: emstand-V09-03-24 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 66 66 #include "G4Electron.hh" 67 67 #include "G4Proton.hh" 68 //#include "G4Alpha.hh"69 68 #include "G4GenericIon.hh" 70 69 #include "G4BraggModel.hh" … … 83 82 G4ionIonisation::G4ionIonisation(const G4String& name) 84 83 : G4VEnergyLossProcess(name), 85 corr(0),86 84 theParticle(0), 87 85 isInitialised(false), 88 stopDataActive(true), 89 nuclearStopping(true) 86 stopDataActive(true) 90 87 { 91 88 SetLinearLossLimit(0.02); … … 95 92 // SetVerboseLevel(1); 96 93 corr = G4LossTableManager::Instance()->EmCorrections(); 94 eth = 2*MeV; 97 95 } 98 96 … … 136 134 const G4ParticleDefinition* theBaseParticle = 0; 137 135 138 if(part == ion) theBaseParticle = 0;139 else if(bpart == 0) theBaseParticle = ion;140 else theBaseParticle = bpart;136 if(part == ion) { theBaseParticle = 0; } 137 else if(bpart == 0) { theBaseParticle = ion; } 138 else { theBaseParticle = bpart; } 141 139 142 140 SetBaseParticle(theBaseParticle); 143 141 SetSecondaryParticle(G4Electron::Electron()); 144 142 145 if (!EmModel(1)) SetEmModel(new G4BraggIonModel(), 1);143 if (!EmModel(1)) { SetEmModel(new G4BraggIonModel(), 1); } 146 144 EmModel(1)->SetLowEnergyLimit(MinKinEnergy()); 147 145 … … 150 148 EmModel(1)->SetHighEnergyLimit(eth); 151 149 152 if (!FluctModel()) SetFluctModel(new G4IonFluctuations());150 if (!FluctModel()) { SetFluctModel(new G4IonFluctuations()); } 153 151 AddEmModel(1, EmModel(1), FluctModel()); 154 152 155 if (!EmModel(2)) SetEmModel(new G4BetheBlochModel(),2);153 if (!EmModel(2)) { SetEmModel(new G4BetheBlochModel(),2); } 156 154 EmModel(2)->SetLowEnergyLimit(eth); 157 155 EmModel(2)->SetHighEnergyLimit(MaxKinEnergy()); … … 166 164 } 167 165 // reinitialisation of corrections for the new run 168 EmModel(1)->ActivateNuclearStopping(nuclearStopping); 169 EmModel(2)->ActivateNuclearStopping(nuclearStopping); 170 if(part == ion) corr->InitialiseForNewRun(); 166 if(part == ion) { corr->InitialiseForNewRun(); } 171 167 } 172 168 … … 178 174 G4cout << " Stopping Power data for " 179 175 << corr->GetNumberOfStoppingVectors() 180 << " ion/material pairs , nuclearStopping: " << nuclearStopping176 << " ion/material pairs " 181 177 << G4endl; 182 178 }
Note: See TracChangeset
for help on using the changeset viewer.