6060
6161#include " G4SystemOfUnits.hh"
6262
63- CMSEmStandardPhysicsEMZ::CMSEmStandardPhysicsEMZ (G4int ver)
64- : G4VPhysicsConstructor(" CMSEmStandard_emz" ) {
63+ CMSEmStandardPhysicsEMZ::CMSEmStandardPhysicsEMZ (G4int ver) : G4VPhysicsConstructor(" CMSEmStandard_emz" ) {
6564 SetVerboseLevel (ver);
6665 G4EmParameters* param = G4EmParameters::Instance ();
6766 param->SetDefaults ();
6867 param->SetVerbose (ver);
69- param->SetMinEnergy (100 * CLHEP::eV);
70- param->SetLowestElectronEnergy (100 * CLHEP::eV);
68+ param->SetMinEnergy (100 * CLHEP::eV);
69+ param->SetLowestElectronEnergy (100 * CLHEP::eV);
7170 param->SetNumberOfBinsPerDecade (20 );
7271 param->ActivateAngularGeneratorForIonisation (true );
73- param->SetStepFunction (0.2 , 10 * CLHEP::um);
74- param->SetStepFunctionMuHad (0.1 , 50 * CLHEP::um);
75- param->SetStepFunctionLightIons (0.1 , 20 * CLHEP::um);
76- param->SetStepFunctionIons (0.1 , 1 * CLHEP::um);
77- param->SetUseMottCorrection (true ); // use Mott-correction for e-/e+ msc gs
78- param->SetMscStepLimitType (fUseSafetyPlus ); // for e-/e+ msc gs
79- param->SetMscSkin (3 ); // error-free stepping for e-/e+ msc gs
80- param->SetMscRangeFactor (0.08 ); // error-free stepping for e-/e+ msc gs
72+ param->SetStepFunction (0.2 , 10 * CLHEP::um);
73+ param->SetStepFunctionMuHad (0.1 , 50 * CLHEP::um);
74+ param->SetStepFunctionLightIons (0.1 , 20 * CLHEP::um);
75+ param->SetStepFunctionIons (0.1 , 1 * CLHEP::um);
76+ param->SetUseMottCorrection (true ); // use Mott-correction for e-/e+ msc gs
77+ param->SetMscStepLimitType (fUseSafetyPlus ); // for e-/e+ msc gs
78+ param->SetMscSkin (3 ); // error-free stepping for e-/e+ msc gs
79+ param->SetMscRangeFactor (0.08 ); // error-free stepping for e-/e+ msc gs
8180 param->SetMuHadLateralDisplacement (true );
8281 param->SetFluo (true );
8382 param->SetUseICRU90Data (true );
84- param->SetMaxNIELEnergy (1 * CLHEP::MeV);
83+ param->SetMaxNIELEnergy (1 * CLHEP::MeV);
8584 SetPhysicsType (bElectromagnetic);
86-
8785}
8886
8987CMSEmStandardPhysicsEMZ::~CMSEmStandardPhysicsEMZ () {}
@@ -110,7 +108,7 @@ void CMSEmStandardPhysicsEMZ::ConstructProcess() {
110108 // G4NuclearStopping* pnuc(nullptr);
111109 G4double nielEnergyLimit = G4EmParameters::Instance ()->MaxNIELEnergy ();
112110 G4NuclearStopping* pnuc = nullptr ;
113- if (nielEnergyLimit > 0.0 ) {
111+ if (nielEnergyLimit > 0.0 ) {
114112 pnuc = new G4NuclearStopping ();
115113 pnuc->SetMaxKinEnergy (nielEnergyLimit);
116114 }
@@ -130,7 +128,7 @@ void CMSEmStandardPhysicsEMZ::ConstructProcess() {
130128 G4ComptonScattering* cs = new G4ComptonScattering;
131129 cs->SetEmModel (new G4KleinNishinaModel ());
132130 G4VEmModel* theLowEPComptonModel = new G4LowEPComptonModel ();
133- theLowEPComptonModel->SetHighEnergyLimit (20 * CLHEP::MeV);
131+ theLowEPComptonModel->SetHighEnergyLimit (20 * CLHEP::MeV);
134132 cs->AddEmModel (0 , theLowEPComptonModel);
135133
136134 // Gamma conversion
@@ -167,9 +165,9 @@ void CMSEmStandardPhysicsEMZ::ConstructProcess() {
167165 msc->SetEmModel (msc1);
168166 msc->SetEmModel (msc2);
169167
170- G4eCoulombScatteringModel* ssm = new G4eCoulombScatteringModel ();
168+ G4eCoulombScatteringModel* ssm = new G4eCoulombScatteringModel ();
171169 G4CoulombScattering* ss = new G4CoulombScattering ();
172- ss->SetEmModel (ssm);
170+ ss->SetEmModel (ssm);
173171 ss->SetMinKinEnergy (highEnergyLimit);
174172 ssm->SetLowEnergyLimit (highEnergyLimit);
175173 ssm->SetActivationLowEnergyLimit (highEnergyLimit);
@@ -185,8 +183,8 @@ void CMSEmStandardPhysicsEMZ::ConstructProcess() {
185183 // ionisation
186184 G4eIonisation* eioni = new G4eIonisation ();
187185 G4VEmModel* theIoniLiv = new G4LivermoreIonisationModel ();
188- theIoniLiv->SetHighEnergyLimit (0.1 * CLHEP::MeV);
189- eioni->AddEmModel (0 , theIoniLiv, new G4UniversalFluctuation () );
186+ theIoniLiv->SetHighEnergyLimit (0.1 * CLHEP::MeV);
187+ eioni->AddEmModel (0 , theIoniLiv, new G4UniversalFluctuation ());
190188
191189 // bremsstrahlung
192190 G4eBremsstrahlung* brem = new G4eBremsstrahlung ();
@@ -260,7 +258,8 @@ void CMSEmStandardPhysicsEMZ::ConstructProcess() {
260258 ionIoni->SetEmModel (new G4IonParametrisedLossModel ());
261259 ph->RegisterProcess (hmsc, particle);
262260 ph->RegisterProcess (ionIoni, particle);
263- if (nullptr != pnuc) ph->RegisterProcess (pnuc, particle);
261+ if (nullptr != pnuc)
262+ ph->RegisterProcess (pnuc, particle);
264263
265264 // muons, hadrons, ions
266265 G4EmBuilder::ConstructCharged (hmsc, pnuc);
0 commit comments