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