Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4EmLEPTSPhysics Class Reference

#include <G4EmLEPTSPhysics.hh>

Inheritance diagram for G4EmLEPTSPhysics:
Collaboration diagram for G4EmLEPTSPhysics:

Public Member Functions

 G4EmLEPTSPhysics (const G4String &name="G4EmLEPTSPhysics")
 
virtual ~G4EmLEPTSPhysics ()
 
virtual void ConstructParticle ()
 
virtual void ConstructProcess ()
 
- Public Member Functions inherited from G4VPhysicsConstructor
 G4VPhysicsConstructor (const G4String &="")
 
 G4VPhysicsConstructor (const G4String &name, G4int physics_type)
 
virtual ~G4VPhysicsConstructor ()
 
void SetPhysicsName (const G4String &="")
 
const G4StringGetPhysicsName () const
 
void SetPhysicsType (G4int)
 
G4int GetPhysicsType () const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
G4int GetInstanceID () const
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManagerGetSubInstanceManager ()
 
- Protected Member Functions inherited from G4VPhysicsConstructor
G4bool RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
 
G4ParticleTable::G4PTblDicIteratorGetParticleIterator () const
 
- Protected Attributes inherited from G4VPhysicsConstructor
G4int verboseLevel
 
G4String namePhysics
 
G4int typePhysics
 
G4ParticleTabletheParticleTable
 
G4int g4vpcInstanceID
 
- Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager subInstanceManager
 

Detailed Description

Definition at line 46 of file G4EmLEPTSPhysics.hh.

Constructor & Destructor Documentation

G4EmLEPTSPhysics::G4EmLEPTSPhysics ( const G4String name = "G4EmLEPTSPhysics")
explicit

Definition at line 68 of file G4EmLEPTSPhysics.cc.

69  : G4VPhysicsConstructor(name)
70 {
73 }
static G4EmParameters * Instance()
G4VPhysicsConstructor(const G4String &="")

Here is the call graph for this function:

virtual G4EmLEPTSPhysics::~G4EmLEPTSPhysics ( )
inlinevirtual

Definition at line 51 of file G4EmLEPTSPhysics.hh.

51 {};

Member Function Documentation

void G4EmLEPTSPhysics::ConstructParticle ( )
virtual

Implements G4VPhysicsConstructor.

Definition at line 76 of file G4EmLEPTSPhysics.cc.

77 {
81 
82 // baryons
84 
86 
87  G4DNAGenericIonsManager * genericIonsManager;
88  genericIonsManager=G4DNAGenericIonsManager::Instance();
89  genericIonsManager->GetIon("alpha++");
90  genericIonsManager->GetIon("alpha+");
91  genericIonsManager->GetIon("helium");
92  genericIonsManager->GetIon("hydrogen");
93  genericIonsManager->GetIon("carbon");
94  genericIonsManager->GetIon("nitrogen");
95  genericIonsManager->GetIon("oxygen");
96  genericIonsManager->GetIon("iron");
97 
98 }
static G4GenericIon * GenericIonDefinition()
Definition: G4GenericIon.cc:88
static G4Proton * Proton()
Definition: G4Proton.cc:93
static G4DNAGenericIonsManager * Instance(void)
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
static G4Positron * Positron()
Definition: G4Positron.cc:94
static G4Electron * Electron()
Definition: G4Electron.cc:94
G4ParticleDefinition * GetIon(const G4String &name)

Here is the call graph for this function:

void G4EmLEPTSPhysics::ConstructProcess ( )
virtual

Implements G4VPhysicsConstructor.

Definition at line 101 of file G4EmLEPTSPhysics.cc.

102 {
103 
104  auto myParticleIterator=GetParticleIterator();
105  myParticleIterator->reset();
106 
107  while( (*myParticleIterator)() ){
108  G4ParticleDefinition * particle = myParticleIterator->value();
109  G4String particleName = particle->GetParticleName();
110  G4ProcessManager * manager = particle->GetProcessManager();
111  G4cout << " particle " << particle->GetParticleName() << " manager " << manager << G4endl; //GDEB
112  if(particleName=="e-") {
113  G4DNAElastic* elastic = new G4DNAElastic("e-_G4LEPTSElastic");
114  elastic->SetEmModel( new G4LEPTSElasticModel );
115  manager->AddDiscreteProcess(elastic );
116 
117  G4DNAIonisation* ioni = new G4DNAIonisation("e-_G4LEPTSIonisation");
118  ioni->SetEmModel( new G4LEPTSIonisationModel );
119  manager->AddDiscreteProcess( ioni );
120 
121  G4DNAExcitation* excit = new G4DNAExcitation("e-_G4LEPTSExcitation");
122  excit->SetEmModel( new G4LEPTSExcitationModel );
123  manager->AddDiscreteProcess( excit );
124 
125  manager->AddDiscreteProcess(new G4DNADissociation("e-_G4LEPTSDissocNeutr") );
126 
127  G4DNAVibExcitation* vibExcit = new G4DNAVibExcitation("e-_G4LEPTSExcitVibrat");
128  vibExcit->SetEmModel( new G4LEPTSVibExcitationModel );
129  manager->AddDiscreteProcess( vibExcit );
130 
131  manager->AddDiscreteProcess( new G4DNARotExcitation("e-_G4LEPTSExcitRotat") );
132 
133  G4DNAAttachment* attach = new G4DNAAttachment("e-_G4LEPTSAttachment");
134  attach->SetEmModel( new G4LEPTSAttachmentModel );
135  manager->AddDiscreteProcess( attach );
136 
137  } else if(particleName=="e+") {
138  G4DNAElastic* elastic = new G4DNAElastic("e+_G4LEPTSElastic");
139  elastic->SetEmModel( new G4LEPTSElasticModel );
140  manager->AddDiscreteProcess(elastic );
141 
142  G4DNAIonisation* ioni = new G4DNAIonisation("e+_G4LEPTSIonisation");
143  ioni->SetEmModel( new G4LEPTSIonisationModel );
144  manager->AddDiscreteProcess( ioni );
145 
146  G4DNAExcitation* excit = new G4DNAExcitation("e+_G4LEPTSExcitation");
147  excit->SetEmModel( new G4LEPTSExcitationModel );
148  manager->AddDiscreteProcess( excit );
149 
150  manager->AddDiscreteProcess(new G4DNADissociation("e+_G4LEPTSDissocNeutr") );
151 
152  G4DNAVibExcitation* vibExcit = new G4DNAVibExcitation("e+_G4LEPTSExcitVibrat");
153  vibExcit->SetEmModel( new G4LEPTSVibExcitationModel );
154  manager->AddDiscreteProcess( vibExcit );
155 
156  manager->AddDiscreteProcess( new G4DNARotExcitation("e+_G4LEPTSExcitRotat") );
157 
158  G4DNAAttachment* attach = new G4DNAAttachment("e+_G4LEPTSAttachment");
159  attach->SetEmModel( new G4LEPTSAttachmentModel );
160  manager->AddDiscreteProcess( attach );
161 
162  manager->AddDiscreteProcess(new G4DNAPositronium("e+_G4LEPTSPositronium") );
163 
164  /* } else if ( particleName == "proton" ) {
165  manager->AddDiscreteProcess(new G4DNAExcitation("proton_G4DNAExcitation"));
166  manager->AddDiscreteProcess(new G4DNAIonisation("proton_G4DNAIonisation"));
167  manager->AddDiscreteProcess(new G4DNAChargeDecrease("proton_G4DNAChargeDecrease"));
168 
169  } else if ( particleName == "hydrogen" ) {
170  manager->AddDiscreteProcess(new G4DNAExcitation("hydrogen_G4DNAExcitation"));
171  manager->AddDiscreteProcess(new G4DNAIonisation("hydrogen_G4DNAIonisation"));
172  manager->AddDiscreteProcess(new G4DNAChargeIncrease("hydrogen_G4DNAChargeIncrease"));
173 
174  } else if ( particleName == "alpha" ) {
175  manager->AddDiscreteProcess(new G4DNAExcitation("alpha_G4DNAExcitation"));
176  manager->AddDiscreteProcess(new G4DNAIonisation("alpha_G4DNAIonisation"));
177  manager->AddDiscreteProcess(new G4DNAChargeDecrease("alpha_G4DNAChargeDecrease"));
178 
179  } else if ( particleName == "alpha+" ) {
180  manager->AddDiscreteProcess(new G4DNAExcitation("alpha+_G4DNAExcitation"));
181  manager->AddDiscreteProcess(new G4DNAIonisation("alpha+_G4DNAIonisation"));
182  manager->AddDiscreteProcess(new G4DNAChargeDecrease("alpha+_G4DNAChargeDecrease"));
183  manager->AddDiscreteProcess(new G4DNAChargeIncrease("alpha+_G4DNAChargeIncrease"));
184 
185  } else if ( particleName == "helium" ) {
186  manager->AddDiscreteProcess(new G4DNAExcitation("helium_G4DNAExcitation"));
187  manager->AddDiscreteProcess(new G4DNAIonisation("helium_G4DNAIonisation"));
188  manager->AddDiscreteProcess(new G4DNAChargeIncrease("helium_G4DNAChargeIncrease"));
189 
190  // Extension to HZE proposed by Z. Francis
191 
192  } else if ( particleName == "carbon" ) {
193  manager->AddDiscreteProcess(new G4DNAIonisation("carbon_G4DNAIonisation"));
194 
195  } else if ( particleName == "nitrogen" ) {
196  manager->AddDiscreteProcess(new G4DNAIonisation("nitrogen_G4DNAIonisation"));
197 
198  } else if ( particleName == "oxygen" ) {
199  manager->AddDiscreteProcess(new G4DNAIonisation("oxygen_G4DNAIonisation"));
200 
201  } else if ( particleName == "iron" ) {
202  manager->AddDiscreteProcess(new G4DNAIonisation("iron_G4DNAIonisation"));
203  */
204  } else if (particleName == "gamma") {
205 
206  G4double LivermoreHighEnergyLimit = CLHEP::GeV;
207 
208  G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect();
209  G4LivermorePhotoElectricModel* theLivermorePhotoElectricModel =
211  theLivermorePhotoElectricModel->SetHighEnergyLimit(LivermoreHighEnergyLimit);
212  thePhotoElectricEffect->AddEmModel(0, theLivermorePhotoElectricModel);
213  manager->AddDiscreteProcess(thePhotoElectricEffect);
214 
215  G4ComptonScattering* theComptonScattering = new G4ComptonScattering();
216  G4LivermoreComptonModel* theLivermoreComptonModel =
218  theLivermoreComptonModel->SetHighEnergyLimit(LivermoreHighEnergyLimit);
219  theComptonScattering->AddEmModel(0, theLivermoreComptonModel);
220  manager->AddDiscreteProcess(theComptonScattering);
221 
222  G4GammaConversion* theGammaConversion = new G4GammaConversion();
223  G4LivermoreGammaConversionModel* theLivermoreGammaConversionModel =
225  theLivermoreGammaConversionModel->SetHighEnergyLimit(LivermoreHighEnergyLimit);
226  theGammaConversion->AddEmModel(0, theLivermoreGammaConversionModel);
227  manager->AddDiscreteProcess(theGammaConversion);
228 
229  G4RayleighScattering* theRayleigh = new G4RayleighScattering();
230  G4LivermoreRayleighModel* theRayleighModel = new G4LivermoreRayleighModel();
231  theRayleighModel->SetHighEnergyLimit(LivermoreHighEnergyLimit);
232  theRayleigh->AddEmModel(0, theRayleighModel);
233  manager->AddDiscreteProcess(theRayleigh);
234  }
235 
236  }
237 
238 }
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:724
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
void SetEmModel(G4VEmModel *, G4int index=1)
G4GLOB_DLL std::ostream G4cout
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=nullptr)
static constexpr double GeV
G4ProcessManager * GetProcessManager() const
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76
G4double elastic(Particle const *const p1, Particle const *const p2)

Here is the call graph for this function:


The documentation for this class was generated from the following files: