114{
115 if(!isInitialised)
116 {
117 theParticle = part;
118 theBaseParticle = bpart;
119
120 mass = theParticle->GetPDGMass();
121 SetSecondaryParticle(G4Electron::Electron());
122
123 flucModel = new G4UniversalFluctuation();
124
125 G4VEmModel* em = new G4BraggModel();
126 em->SetLowEnergyLimit(0.1*keV);
127 em->SetHighEnergyLimit(0.2*MeV);
128 AddEmModel(1, em, flucModel);
129 G4VEmModel* em1 = new G4BetheBlochModel();
130 em1->SetLowEnergyLimit(0.2*MeV);
131 em1->SetHighEnergyLimit(1.0*GeV);
132 AddEmModel(2, em1, flucModel);
133 G4VEmModel* em2 = new G4MuBetheBlochModel();
134 em2->SetLowEnergyLimit(1.0*GeV);
135 em2->SetHighEnergyLimit(100.0*TeV);
136 AddEmModel(3, em2, flucModel);
137 SetStepFunction(0.2, 1.0*mm);
138 ratio = electron_mass_c2/mass;
139 isInitialised = true;
140 }
141}