ATLAS Offline Software
G4mplAtlasIonisation.cxx
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4mplAtlasIonisation.cxx 729653 2016-03-14 15:55:47Z jchapman $
27 // GEANT4 tag $Name: not supported by cvs2svn $
28 //
29 // -------------------------------------------------------------------
30 //
31 // GEANT4 Class file
32 //
33 //
34 // File name: G4mplAtlasIonisation
35 //
36 // Author: Vladimir Ivanchenko
37 //
38 // Creation date: 25.08.2005
39 //
40 // Modifications:
41 //
42 //
43 // -------------------------------------------------------------------
44 //
45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
47 
48 // class header
49 #include "G4mplAtlasIonisation.hh"
50 // package headers
51 #include "G4mplAtlasIonisationWithDeltaModel.hh"
52 // Geant4 headers
53 #include "G4Version.hh"
54 
55 #if G4VERSION_NUMBER < 1100
56 #include "G4BohrFluctuations.hh"
57 #endif
58 
59 #include "G4Electron.hh"
60 #include "G4EmProcessSubType.hh"
61 #include "G4Version.hh"
62 // CLHEP headers
63 #include "CLHEP/Units/SystemOfUnits.h"
64 #include "CLHEP/Units/PhysicalConstants.h"
65 
66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
67 
68 G4mplAtlasIonisation::G4mplAtlasIonisation(G4double mCharge, const G4String& name)
69  : G4VEnergyLossProcess(name),
70  magneticCharge(mCharge),
71  isInitialised(false)
72 {
73 
74  G4cout <<"!!! G4mplAtlasIonisation constructor"<<G4endl;
75 
76  SetProcessSubType( fIonisation );
77 
78  // By default classical magnetic charge is used
79  if(magneticCharge == 0.0) magneticCharge = CLHEP::eplus*0.5/CLHEP::fine_structure_const;
80 
81  SetDEDXBinning(120);
82 #if G4VERSION_NUMBER < 1010
83  SetLambdaBinning(120); //G4 9.6 only
84 #endif
85  SetMinKinEnergy(0.1*CLHEP::keV);
86  SetMaxKinEnergy(100.0*CLHEP::TeV);
87  SetVerboseLevel(0);
88 }
89 
90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
91 
92 G4mplAtlasIonisation::~G4mplAtlasIonisation()
93 {
94  G4cout <<"!!! G4mplAtlasIonisation destructor"<<G4endl;
95 }
96 
97 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
98 
99 void G4mplAtlasIonisation::InitialiseEnergyLossProcess(const G4ParticleDefinition*,
100  const G4ParticleDefinition*)
101 {
102  if(isInitialised) return;
103 
104  G4cout <<"!!! G4mplAtlasIonisation::InitialiseEnergyLossProcess"<<G4endl;
105 
106  SetBaseParticle(0);
107  SetSecondaryParticle(G4Electron::Electron());
108 
109  // G4mplIonisationModel* ion = new G4mplIonisationModel(magneticCharge);
110  // G4mplIonisationModel* ion = new G4mplIonisationModel(magneticCharge,"PAI");
111  // G4mplAtlasIonisationModel* ion = new G4mplAtlasIonisationModel(magneticCharge,"PAI");
112  G4mplAtlasIonisationWithDeltaModel* ion = new G4mplAtlasIonisationWithDeltaModel(magneticCharge,"PAI");
113  ion->SetLowEnergyLimit(0.1*CLHEP::keV);
114  // ion->SetLowEnergyLimit(1000.0*keV);
115  // ion->SetHighEnergyLimit(100.*TeV);
116  ion->SetHighEnergyLimit(7.*CLHEP::TeV);
117  AddEmModel(0,ion,ion);
118 
119  SetStepFunction(0.2, 1*CLHEP::mm);
120 
121  isInitialised = true;
122 
123  G4cout <<"!!! G4mplAtlasIonisation::InitialiseEnergyLossProcess done"<<G4endl;
124 
125 }
126 
127 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
128 
129 void G4mplAtlasIonisation::PrintInfo()
130 {
131  // G4cout << "SB: My No delta-electron production, only dE/dx"
132  G4cout << "!!! G4mplAtlasIonisation: Delta-electron production"
133  << G4endl;
134 }
135 
136 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
xAOD::Electron
Electron_v1 Electron
Definition of the current "egamma version".
Definition: Event/xAOD/xAODEgamma/xAODEgamma/Electron.h:17
python.PhysicalConstants.fine_structure_const
float fine_structure_const
Definition: PhysicalConstants.py:103
python.SystemOfUnits.TeV
int TeV
Definition: SystemOfUnits.py:158
python.SystemOfUnits.keV
int keV
Definition: SystemOfUnits.py:156
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
python.SystemOfUnits.mm
int mm
Definition: SystemOfUnits.py:83
python.SystemOfUnits.eplus
int eplus
Definition: SystemOfUnits.py:137