Geant4  10.02.p03
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::G4EmLEPTSPhysics ( const G4String name = "G4EmLEPTSPhysics")

Definition at line 68 of file G4EmLEPTSPhysics.cc.

69  : G4VPhysicsConstructor(name)
70 {
72 }
static G4EmParameters * Instance()
G4VPhysicsConstructor(const G4String &="")
Here is the call graph for this function:

◆ ~G4EmLEPTSPhysics()

virtual G4EmLEPTSPhysics::~G4EmLEPTSPhysics ( )
inlinevirtual

Definition at line 50 of file G4EmLEPTSPhysics.hh.

50 {};
Here is the call graph for this function:

Member Function Documentation

◆ ConstructParticle()

void G4EmLEPTSPhysics::ConstructParticle ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 75 of file G4EmLEPTSPhysics.cc.

76 {
80 
81 // baryons
83 
85 
86  G4DNAGenericIonsManager * genericIonsManager;
87  genericIonsManager=G4DNAGenericIonsManager::Instance();
88  genericIonsManager->GetIon("alpha++");
89  genericIonsManager->GetIon("alpha+");
90  genericIonsManager->GetIon("helium");
91  genericIonsManager->GetIon("hydrogen");
92  genericIonsManager->GetIon("carbon");
93  genericIonsManager->GetIon("nitrogen");
94  genericIonsManager->GetIon("oxygen");
95  genericIonsManager->GetIon("iron");
96 
97 }
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:
Here is the caller graph for this function:

◆ ConstructProcess()

void G4EmLEPTSPhysics::ConstructProcess ( void  )
virtual

Implements G4VPhysicsConstructor.

Definition at line 100 of file G4EmLEPTSPhysics.cc.

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

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