Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Classes | Public Member Functions | List of all members
RunAction Class Reference

Run action class. More...

#include <RunAction.hh>

Inherits G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, G4UserRunAction, and G4UserRunAction.

Collaboration diagram for RunAction:
Collaboration graph
[legend]

Public Member Functions

 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
void BeginOfRunAction (const G4Run *)
 
void EndOfRunAction (const G4Run *)
 
void GetCuts ()
 
void CriticalEnergy ()
 
 RunAction (DetectorConstruction *det, PrimaryGeneratorAction *prim=0)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *)
 
 ~RunAction ()
 
void BeginOfRunAction (const G4Run *)
 
void EndOfRunAction (const G4Run *)
 
void SetRndmFreq (G4int val)
 
G4int GetRndmFreq ()
 
void AddDoseN (G4double dose)
 
void SetDoseN (G4double dose)
 
G4double GetDoseN ()
 
void AddDoseC (G4double dose)
 
void SetDoseC (G4double dose)
 
G4double GetDoseC ()
 
G4int GetNumEvent ()
 
void SetNumEvent (G4int i)
 
G4int GetNbOfHitsGas ()
 
void AddNbOfHitsGas ()
 
void SetMassNucleus (G4double mN)
 
G4double GetMassNucleus ()
 
void SetMassCytoplasm (G4double mC)
 
G4double GetMassCytoplasm ()
 
void AddDoseBox (G4int i, G4double x)
 
G4double GetDoseBox (G4int i)
 
G4ThreeVector GetVectCell (G4int i)
 
 RunAction ()
 
 ~RunAction ()
 
void BeginOfRunAction (const G4Run *)
 
void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
void BeginOfRunAction (const G4Run *)
 
void EndOfRunAction (const G4Run *)
 
G4int GetRow ()
 
void AddRow ()
 
void AddToXVector (float v)
 
void AddToYVector (float v)
 
void AddToThetaVector (float v)
 
void AddToPhiVector (float v)
 
 RunAction ()
 
 ~RunAction ()
 
void BeginOfRunAction (const G4Run *)
 
void EndOfRunAction (const G4Run *)
 
 RunAction (HistoManager *)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void FillPerEvent (G4double, G4double, G4double, G4double)
 
 RunAction (HistoManager *)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void FillPerEvent (G4double, G4double, G4double, G4double)
 
 RunAction (HistoManager *)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void FillPerEvent (G4double, G4double, G4double, G4double)
 
 RunAction (HistoManager *)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void FillPerEvent (G4double, G4double, G4double, G4double)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void GetCuts ()
 
void CriticalEnergy ()
 
 RunAction (DetectorConstruction *det, PrimaryGeneratorAction *prim=0)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void SetRndmFreq (G4int val)
 
G4int GetRndmFreq ()
 
 RunAction (DetectorConstruction *, PhysicsList *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
virtual G4Run * GenerateRun ()
 
 RunAction (DetectorConstruction *, PhysicsList *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
virtual G4Run * GenerateRun ()
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *det, PrimaryGeneratorAction *prim=0)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void CountProcesses (G4String)
 
void SumPathLength (G4double truepl, G4double geompl)
 
void SumLateralDisplacement (G4double displa)
 
void SumPsi (G4double psi)
 
void SumTetaPlane (G4double teta)
 
void SumPhiCorrel (G4double correl)
 
G4double ComputeMscHighland (G4double pathLength)
 
 RunAction ()
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *, HistoManager *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void CountProcesses (G4String)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void AddEnergyDeposit (G4double edep)
 
void AddTrackLength (G4double step)
 
void AddChargedSecondary (G4double ekin)
 
void AddNeutralSecondary (G4double ekin)
 
G4double GetEnergyFromRestrictedRange (G4double, G4ParticleDefinition *, G4Material *, G4double)
 
G4double GetEnergyFromCSDARange (G4double, G4ParticleDefinition *, G4Material *, G4double)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
virtual ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void SetVerbose (G4int val)
 
void SetHistoName (G4String &val)
 
void SetEdepAndRMS (G4ThreeVector)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *prim=0)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void SetEdepAndRMS (G4int, G4double, G4double, G4double)
 
void SetApplyLimit (G4bool val)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *det, PrimaryGeneratorAction *prim=0)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void CountProcesses (G4String)
 
 RunAction (DetectorConstruction *, PhysicsList *, PrimaryGeneratorAction *)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void FillTallyEdep (G4int n, G4double e)
 
void FillEdep (G4double de, G4double eni)
 
void AddProjRange (G4double x)
 
void AddPrimaryStep ()
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void FillHisto (G4int id, G4double x, G4double weight=1.0)
 
void SetBinSize (G4double size)
 
G4double GetOffsetX ()
 
void SetVerbose (G4int verbose)
 
G4int GetVerbose ()
 
void AddProjRange (G4double x)
 
 RunAction (const DetectorConstruction *det=0)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void SetPrintFlag (G4bool)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
virtual G4Run * GenerateRun ()
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void CreateHistogram ()
 
void WriteHistogram ()
 
void SetHistoName (G4String &val)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *run)
 
virtual void EndOfRunAction (const G4Run *run)
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
virtual G4Run * GenerateRun ()
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
virtual G4Run * GenerateRun ()
 
 RunAction ()
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
virtual G4Run * GenerateRun ()
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *prim=0)
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
virtual G4Run * GenerateRun ()
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
virtual ~RunAction ()
 
void BeginOfRunAction (const G4Run *)
 
void EndOfRunAction (const G4Run *)
 
 RunAction ()
 
 ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *arun)
 
virtual void EndOfRunAction (const G4Run *arun)
 
 RunAction ()
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *arun)
 
virtual void EndOfRunAction (const G4Run *arun)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
virtual ~RunAction ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
void CountProcesses (G4String)
 
void FillData (const G4ParticleDefinition *particle, G4double kinEnergy, G4double costheta, G4double phi, G4double longitudinalPolarization)
 
void EventFinished ()
 
 RunAction (PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
 RunAction (DetectorConstruction *, PrimaryGeneratorAction *)
 
 ~RunAction ()
 
virtual G4Run * GenerateRun ()
 
virtual void BeginOfRunAction (const G4Run *)
 
virtual void EndOfRunAction (const G4Run *)
 
- Public Member Functions inherited from G4UserRunAction
 G4UserRunAction ()
 
virtual ~G4UserRunAction ()
 
virtual void SetMaster (G4bool val=true)
 
G4bool IsMaster () const
 

Additional Inherited Members

- Protected Attributes inherited from G4UserRunAction
G4bool isMaster
 

Detailed Description

Run action class.

Definition at line 45 of file RunAction.hh.

Constructor & Destructor Documentation

RunAction::RunAction ( DetectorConstruction *  det,
PrimaryGeneratorAction *  kin 
)

Definition at line 46 of file RunAction.cc.

47 :detector(det), primary(kin)
48 { }
RunAction::~RunAction ( )

Definition at line 52 of file RunAction.cc.

53 { }
RunAction::RunAction ( DetectorConstruction *  det,
PrimaryGeneratorAction *  prim = 0 
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  det)

Definition at line 43 of file RunAction.cc.

44 :fDetector(det)
45 {
46  fSaveRndm = 0;
47 }
RunAction::~RunAction ( )
RunAction::RunAction ( )

Definition at line 42 of file RunAction.cc.

43 {
44  fpTrackingAction = 0;
45  fInitialized = 0;
46  fDebug = false;
47 }
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( HistoManager *  histo)

Definition at line 44 of file RunAction.cc.

45 : G4UserRunAction(),
46  fHistoManager(histo),
47  fSumEAbs(0.), fSum2EAbs(0.),
48  fSumEGap(0.), fSum2EGap(0.),
49  fSumLAbs(0.), fSum2LAbs(0.),
50  fSumLGap(0.), fSum2LGap(0.)
51 {}
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( HistoManager *  )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( HistoManager *  )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( HistoManager *  )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  det,
PrimaryGeneratorAction *  prim = 0 
)
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  det,
PhysicsList *  phys,
PrimaryGeneratorAction *  kin 
)

Definition at line 52 of file RunAction.cc.

54 :G4UserRunAction(),fDetector(det),fPhysics(phys),fPrimary(kin),fRun(0),
55  fHistoManager(0)
56 {
57  // Book predefined histograms
58  fHistoManager = new HistoManager();
59 }
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PhysicsList *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  det,
PrimaryGeneratorAction *  prim = 0 
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  det,
PrimaryGeneratorAction *  prim,
HistoManager *  HistM 
)

Definition at line 51 of file RunAction.cc.

53  : G4UserRunAction(),
54  fDetector(det), fPrimary(prim), fProcCounter(0), fHistoManager(HistM)
55 {}
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *  prim = 0 
)
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  det,
PrimaryGeneratorAction *  prim = 0 
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  )
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PhysicsList *  ,
PrimaryGeneratorAction *   
)
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( const DetectorConstruction *  det = 0)

Definition at line 51 of file RunAction.cc.

52  : fDetector(det)
53 {
54  fHisto = new HistoManager();
55 }
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *  prim = 0 
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )
RunAction::RunAction ( )
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( )
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
virtual RunAction::~RunAction ( )
virtual
RunAction::RunAction ( PrimaryGeneratorAction *  kin)

Definition at line 48 of file RunAction.cc.

49 :G4UserRunAction(),
50  fPrimary(kin), fRun(0), fHistoManager(0)
51 {
52  fHistoManager = new HistoManager();
53 }
RunAction::~RunAction ( )
RunAction::RunAction ( DetectorConstruction *  ,
PrimaryGeneratorAction *   
)
RunAction::~RunAction ( )

Member Function Documentation

void RunAction::AddChargedSecondary ( G4double  ekin)
inline

Definition at line 69 of file RunAction.hh.

70  {fEnergyCharged += ekin; fNbCharged++;
71  if (ekin<fEmin[0]) fEmin[0] = ekin;
72  if (ekin>fEmax[0]) fEmax[0] = ekin;
73  };
void RunAction::AddDoseBox ( G4int  i,
G4double  x 
)
inline

Definition at line 80 of file RunAction.hh.

80 { fDose3DDose[i] +=x;}
test.x
tuple x
Definition: test.py:50
void RunAction::AddDoseC ( G4double  dose)
inline

Definition at line 64 of file RunAction.hh.

64 { fDoseC += dose;}
void RunAction::AddDoseN ( G4double  dose)
inline

Definition at line 60 of file RunAction.hh.

60 { fDoseN += dose;}
void RunAction::AddEnergyDeposit ( G4double  edep)
inline

Definition at line 63 of file RunAction.hh.

64  {fEnergyDeposit += edep;};
void RunAction::AddNbOfHitsGas ( )
inline

Definition at line 72 of file RunAction.hh.

72 {fNbOfHitsGas = fNbOfHitsGas+1;}
void RunAction::AddNeutralSecondary ( G4double  ekin)
inline

Definition at line 75 of file RunAction.hh.

76  {fEnergyNeutral += ekin; fNbNeutral++;
77  if (ekin<fEmin[1]) fEmin[1] = ekin;
78  if (ekin>fEmax[1]) fEmax[1] = ekin;
79  };
void RunAction::AddPrimaryStep ( )
inline

Definition at line 106 of file RunAction.hh.

107 {
108  ++fNbPrimarySteps;
109 }
void RunAction::AddProjRange ( G4double  x)
inline

Definition at line 97 of file RunAction.hh.

98 {
99  fProjRange += x;
100  fProjRange2 += x*x;
101  ++fRange;
102 }
test.x
tuple x
Definition: test.py:50
void RunAction::AddProjRange ( G4double  x)
inline

Definition at line 67 of file RunAction.hh.

67 { fProjRange += x; fProjRange2 += x*x; };
test.x
tuple x
Definition: test.py:50
void RunAction::AddRow ( )
inline

Definition at line 53 of file RunAction.hh.

53 {fRow=fRow+1;}
void RunAction::AddToPhiVector ( float  v)
inline

Definition at line 58 of file RunAction.hh.

58 {fPhiVector(fRow)=v;}
test.v
tuple v
Definition: test.py:18
void RunAction::AddToThetaVector ( float  v)
inline

Definition at line 57 of file RunAction.hh.

57 {fThetaVector(fRow)=v;}
test.v
tuple v
Definition: test.py:18
void RunAction::AddToXVector ( float  v)
inline

Definition at line 55 of file RunAction.hh.

55 {fXVector(fRow)=v;}
test.v
tuple v
Definition: test.py:18
void RunAction::AddToYVector ( float  v)
inline

Definition at line 56 of file RunAction.hh.

56 {fYVector(fRow)=v;}
test.v
tuple v
Definition: test.py:18
void RunAction::AddTrackLength ( G4double  step)
inline

Definition at line 66 of file RunAction.hh.

67  {fTrackLength += step; fNbSteps++;};
virtual void RunAction::BeginOfRunAction ( const G4Run *  arun)
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  arun)
virtual

Reimplemented from G4UserRunAction.

void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  run)
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::BeginOfRunAction ( const G4Run *  run)
virtual

G4cout << "\n---> stepMax from histos 1 and 8 = " << G4BestUnit(stepMax,"Length") << G4endl;

Reimplemented from G4UserRunAction.

Definition at line 57 of file RunAction.cc.

58 {
59  //set precision for printing
60  G4int prec = G4cout.precision(6);
61 
62  // get particle
63  G4ParticleDefinition* particle = primary->GetParticleGun()
64  ->GetParticleDefinition();
65  G4String partName = particle->GetParticleName();
66  G4double charge = particle->GetPDGCharge();
67  G4double energy = primary->GetParticleGun()->GetParticleEnergy();
68 
69  // get material
70  G4Material* material = detector->GetMaterial();
71  G4String matName = material->GetName();
72  G4double density = material->GetDensity();
73  G4double radl = material->GetRadlen();
74 
75  G4cout << "\n " << partName << " ("
76  << G4BestUnit(energy,"Energy") << ") in "
77  << material->GetName() << " (density: "
78  << G4BestUnit(density,"Volumic Mass") << "; radiation length: "
79  << G4BestUnit(radl, "Length") << ")" << G4endl;
80 
81  // get cuts
82  GetCuts();
83  if (charge != 0.) {
84  G4cout << "\n Range cuts : \t gamma "
85  << std::setw(8) << G4BestUnit(rangeCut[0],"Length")
86  << "\t e- " << std::setw(8) << G4BestUnit(rangeCut[1],"Length");
87  G4cout << "\n Energy cuts : \t gamma "
88  << std::setw(8) << G4BestUnit(energyCut[0],"Energy")
89  << "\t e- " << std::setw(8) << G4BestUnit(energyCut[1],"Energy")
90  << G4endl;
91  }
92 
93  // get processList and extract EM processes (but not MultipleScattering)
94  G4ProcessVector* plist = particle->GetProcessManager()->GetProcessList();
95  G4String procName;
96  G4double cut;
97  std::vector<G4String> emName;
98  std::vector<G4double> enerCut;
99  size_t length = plist->size();
100  for (size_t j=0; j<length; j++) {
101  procName = (*plist)[j]->GetProcessName();
102  cut = energyCut[1];
103  if ((procName == "eBrem")||(procName == "muBrems")) cut = energyCut[0];
104  if (((*plist)[j]->GetProcessType() == fElectromagnetic) &&
105  (procName != "msc")) {
106  emName.push_back(procName);
107  enerCut.push_back(cut);
108  }
109  }
110 
111  // print list of processes
112  G4cout << "\n processes : ";
113  for (size_t j=0; j<emName.size();j++)
114  G4cout << "\t" << std::setw(13) << emName[j] << "\t";
115  G4cout << "\t" << std::setw(13) <<"total";
116 
117  //instanciate EmCalculator
118  G4EmCalculator emCal;
119  // emCal.SetVerbose(2);
120 
121  //compute cross section per atom (only for single material)
122  if (material->GetNumberOfElements() == 1) {
123  G4double Z = material->GetZ();
124  G4double A = material->GetA();
125 
126  std::vector<G4double> sigma0;
127  G4double sig, sigtot = 0.;
128 
129  for (size_t j=0; j<emName.size();j++) {
130  sig = emCal.ComputeCrossSectionPerAtom
131  (energy,particle,emName[j],Z,A,enerCut[j]);
132  sigtot += sig;
133  sigma0.push_back(sig);
134  }
135  sigma0.push_back(sigtot);
136 
137  G4cout << "\n \n cross section per atom : ";
138  for (size_t j=0; j<sigma0.size();j++) {
139  G4cout << "\t" << std::setw(13) << G4BestUnit(sigma0[j], "Surface");
140  }
141  G4cout << G4endl;
142  }
143 
144  //get cross section per volume
145  std::vector<G4double> sigma1;
146  std::vector<G4double> sigma2;
147  G4double Sig, Sigtot = 0.;
148 
149  for (size_t j=0; j<emName.size();j++) {
150  Sig = emCal.GetCrossSectionPerVolume(energy,particle,emName[j],material);
151  if (Sig == 0.) Sig = emCal.ComputeCrossSectionPerVolume
152  (energy,particle,emName[j],material,enerCut[j]);
153  Sigtot += Sig;
154  sigma1.push_back(Sig);
155  sigma2.push_back(Sig/density);
156  }
157  sigma1.push_back(Sigtot);
158  sigma2.push_back(Sigtot/density);
159 
160  //print cross sections
161  G4cout << "\n \n cross section per volume : ";
162  for (size_t j=0; j<sigma1.size();j++) {
163  G4cout << "\t" << std::setw(13) << sigma1[j]*cm << " cm^-1";
164  }
165 
166  G4cout << "\n cross section per mass : ";
167  for (size_t j=0; j<sigma2.size();j++) {
168  G4cout << "\t" << std::setw(13) << G4BestUnit(sigma2[j], "Surface/Mass");
169  }
170 
171  //print mean free path
172 
173  G4double lambda;
174 
175  G4cout << "\n \n mean free path : ";
176  for (size_t j=0; j<sigma1.size();j++) {
177  lambda = DBL_MAX;
178  if (sigma1[j] > 0.) lambda = 1/sigma1[j];
179  G4cout << "\t" << std::setw(13) << G4BestUnit( lambda, "Length");
180  }
181 
182  //mean free path (g/cm2)
183  G4cout << "\n (g/cm2) : ";
184  for (size_t j=0; j<sigma2.size();j++) {
185  lambda = DBL_MAX;
186  if (sigma2[j] > 0.) lambda = 1/sigma2[j];
187  G4cout << "\t" << std::setw(13) << G4BestUnit( lambda, "Mass/Surface");
188  }
189  G4cout << G4endl;
190 
191  if (charge == 0.) {
192  G4cout.precision(prec);
193  G4cout << "\n-------------------------------------------------------------\n"
194  << G4endl;
195  return;
196  }
197 
198  //get stopping power
199  std::vector<G4double> dedx1;
200  std::vector<G4double> dedx2;
201  G4double dedx, dedxtot = 0.;
202 
203  for (size_t j=0; j<emName.size();j++) {
204  dedx = emCal.ComputeDEDX(energy,particle,emName[j],material,enerCut[j]);
205  dedx1.push_back(dedx);
206  dedx2.push_back(dedx/density);
207  }
208  dedxtot = emCal.GetDEDX(energy,particle,material);
209  dedx1.push_back(dedxtot);
210  dedx2.push_back(dedxtot/density);
211 
212  //print stopping power
213  G4cout << "\n \n restricted dE/dx : ";
214  for (size_t j=0; j<sigma1.size();j++) {
215  G4cout << "\t" << std::setw(13) << G4BestUnit(dedx1[j],"Energy/Length");
216  }
217 
218  G4cout << "\n (MeV/g/cm2) : ";
219  for (size_t j=0; j<sigma2.size();j++) {
220  G4cout << "\t" << std::setw(13) << G4BestUnit(dedx2[j],"Energy*Surface/Mass");
221  }
222 
223  //get range from restricted dedx
224  G4double range1 = emCal.GetRangeFromRestricteDEDX(energy,particle,material);
225  G4double range2 = range1*density;
226 
227  //get range from full dedx
228  G4double Range1 = emCal.GetCSDARange(energy,particle,material);
229  G4double Range2 = Range1*density;
230 
231  //print range
232  G4cout << "\n \n range from restrict dE/dx: "
233  << "\t" << std::setw(8) << G4BestUnit(range1,"Length")
234  << " (" << std::setw(8) << G4BestUnit(range2,"Mass/Surface") << ")";
235 
236  G4cout << "\n range from full dE/dx : "
237  << "\t" << std::setw(8) << G4BestUnit(Range1,"Length")
238  << " (" << std::setw(8) << G4BestUnit(Range2,"Mass/Surface") << ")";
239 
240  //get transport mean free path (for multiple scattering)
241  G4double MSmfp1 = emCal.GetMeanFreePath(energy,particle,"msc",material);
242  G4double MSmfp2 = MSmfp1*density;
243 
244  //print transport mean free path
245  G4cout << "\n \n transport mean free path : "
246  << "\t" << std::setw(8) << G4BestUnit(MSmfp1,"Length")
247  << " (" << std::setw(8) << G4BestUnit(MSmfp2,"Mass/Surface") << ")";
248 
249  if (particle == G4Electron::Electron()) CriticalEnergy();
250 
251  G4cout << "\n-------------------------------------------------------------\n";
252  G4cout << G4endl;
253 
254  // reset default precision
255  G4cout.precision(prec);
256 }
G4Material::GetZ
G4double GetZ() const
Definition: G4Material.cc:623
fElectromagnetic
Definition: G4ProcessType.hh:47
G4Material::GetName
const G4String & GetName() const
Definition: G4Material.hh:178
G4Material::GetDensity
G4double GetDensity() const
Definition: G4Material.hh:180
G4EmCalculator::GetCSDARange
G4double GetCSDARange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:204
G4BestUnit
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
Definition: G4SteppingVerbose.cc:53
G4int
int G4int
Definition: G4Types.hh:78
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
eplot.material
string material
Definition: eplot.py:19
CLHEP::prec
static const double prec
Definition: RanecuEngine.cc:58
G4cout
G4GLOB_DLL std::ostream G4cout
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:60
G4EmCalculator::ComputeDEDX
G4double ComputeDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=DBL_MAX)
Definition: G4EmCalculator.cc:381
G4EmCalculator::GetRangeFromRestricteDEDX
G4double GetRangeFromRestricteDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:181
G4EmCalculator::GetDEDX
G4double GetDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:140
cm
static constexpr double cm
Definition: G4SIunits.hh:119
RunAction::CriticalEnergy
void CriticalEnergy()
Definition: RunAction.cc:300
G4Material::GetRadlen
G4double GetRadlen() const
Definition: G4Material.hh:220
G4Material::GetA
G4double GetA() const
Definition: G4Material.cc:636
G4EmCalculator::GetCrossSectionPerVolume
G4double GetCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:273
G4ProcessVector::size
G4int size() const
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
Definition: G4ParticleDefinition.cc:247
PrimaryGeneratorAction::GetParticleGun
G4ParticleGun * GetParticleGun()
Definition: PrimaryGeneratorAction.hh:53
G4INCL::KinematicsUtils::energy
G4double energy(const ThreeVector &p, const G4double m)
Definition: G4INCLKinematicsUtils.cc:157
G4EmCalculator::ComputeCrossSectionPerVolume
G4double ComputeCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
Definition: G4EmCalculator.cc:613
G4ParticleGun::GetParticleDefinition
G4ParticleDefinition * GetParticleDefinition() const
Definition: G4ParticleGun.hh:106
G4Electron::Electron
static G4Electron * Electron()
Definition: G4Electron.cc:94
G4endl
#define G4endl
Definition: G4ios.hh:61
G4EmCalculator::ComputeCrossSectionPerAtom
G4double ComputeCrossSectionPerAtom(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4double Z, G4double A, G4double cut=0.0)
Definition: G4EmCalculator.cc:648
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:186
G4double
double G4double
Definition: G4Types.hh:76
RunAction::GetCuts
void GetCuts()
Definition: RunAction.cc:269
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:124
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4EmCalculator::GetMeanFreePath
G4double GetMeanFreePath(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:332
G4ParticleGun::GetParticleEnergy
G4double GetParticleEnergy() const
Definition: G4ParticleGun.hh:110
G4ProcessManager::GetProcessList
G4ProcessVector * GetProcessList() const
Z
G4int Z
Definition: G4ParticleHPJENDLHEData.cc:262

Here is the call graph for this function:

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::BeginOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

G4double RunAction::ComputeMscHighland ( G4double  pathLength)

Definition at line 230 of file RunAction.cc.

231 {
232  //compute the width of the Gaussian central part of the MultipleScattering
233  //projected angular distribution.
234  //Eur. Phys. Jour. C15 (2000) page 166, formule 23.9
235 
236  G4double t = pathLength/(fDetector->GetMaterial()->GetRadlen());
237  if (t < DBL_MIN) return 0.;
238 
239  G4ParticleGun* particle = fPrimary->GetParticleGun();
240  G4double T = particle->GetParticleEnergy();
241  G4double M = particle->GetParticleDefinition()->GetPDGMass();
242  G4double z = std::abs(particle->GetParticleDefinition()->GetPDGCharge()/eplus);
243 
244  G4double bpc = T*(T+2*M)/(T+M);
245  G4double teta0 = 13.6*MeV*z*std::sqrt(t)*(1.+0.038*std::log(t))/bpc;
246  return teta0;
247 }
eplus
static constexpr double eplus
Definition: G4SIunits.hh:199
G4Material::GetRadlen
G4double GetRadlen() const
Definition: G4Material.hh:220
G4ParticleDefinition::GetPDGMass
G4double GetPDGMass() const
Definition: G4ParticleDefinition.hh:122
PrimaryGeneratorAction::GetParticleGun
G4ParticleGun * GetParticleGun()
Definition: PrimaryGeneratorAction.hh:53
DBL_MIN
#define DBL_MIN
Definition: templates.hh:75
G4ParticleGun::GetParticleDefinition
G4ParticleDefinition * GetParticleDefinition() const
Definition: G4ParticleGun.hh:106
test.z
tuple z
Definition: test.py:28
MeV
static constexpr double MeV
Definition: G4SIunits.hh:214
G4double
double G4double
Definition: G4Types.hh:76
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:124
G4ParticleGun::GetParticleEnergy
G4double GetParticleEnergy() const
Definition: G4ParticleGun.hh:110

Here is the call graph for this function:

void RunAction::CountProcesses ( G4String  procName)

Definition at line 91 of file RunAction.cc.

92 {
93  //does the process already encounted ?
94  size_t nbProc = fProcCounter->size();
95  size_t i = 0;
96  while ((i<nbProc)&&((*fProcCounter)[i]->GetName()!=procName)) i++;
97  if (i == nbProc) fProcCounter->push_back( new OneProcessCount(procName));
98 
99  (*fProcCounter)[i]->Count();
100 }
void RunAction::CountProcesses ( G4String  )
void RunAction::CountProcesses ( G4String  )
void RunAction::CountProcesses ( G4String  )
void RunAction::CreateHistogram ( )
void RunAction::CriticalEnergy ( )

Definition at line 300 of file RunAction.cc.

301 {
302  // compute e- critical energy (Rossi definition) and Moliere radius.
303  // Review of Particle Physics - Eur. Phys. J. C3 (1998) page 147
304  //
305  G4EmCalculator emCal;
306 
307  const G4Material* material = detector->GetMaterial();
308  const G4double radl = material->GetRadlen();
309  G4double ekin = 5*MeV;
310  G4double deioni;
311  G4double err = 1., errmax = 0.001;
312  G4int iter = 0 , itermax = 10;
313  while (err > errmax && iter < itermax) {
314  iter++;
315  deioni = radl*
316  emCal.ComputeDEDX(ekin,G4Electron::Electron(),"eIoni",material);
317  err = std::abs(deioni - ekin)/ekin;
318  ekin = deioni;
319  }
320  G4cout << "\n \n critical energy (Rossi) : "
321  << "\t" << std::setw(8) << G4BestUnit(ekin,"Energy");
322 
323  //Pdg formula (only for single material)
324  G4double pdga[2] = { 610*MeV, 710*MeV };
325  G4double pdgb[2] = { 1.24, 0.92 };
326  G4double EcPdg = 0.;
327 
328  if (material->GetNumberOfElements() == 1) {
329  G4int istat = 0;
330  if (material->GetState() == kStateGas) istat = 1;
331  G4double Zeff = material->GetZ() + pdgb[istat];
332  EcPdg = pdga[istat]/Zeff;
333  G4cout << "\t\t\t (from Pdg formula : "
334  << std::setw(8) << G4BestUnit(EcPdg,"Energy") << ")";
335  }
336 
337  const G4double Es = 21.2052*MeV;
338  G4double rMolier1 = Es/ekin, rMolier2 = rMolier1*radl;
339  G4cout << "\n Moliere radius : "
340  << "\t" << std::setw(8) << rMolier1 << " X0 "
341  << "= " << std::setw(8) << G4BestUnit(rMolier2,"Length");
342 
343  if (material->GetNumberOfElements() == 1) {
344  G4double rMPdg = radl*Es/EcPdg;
345  G4cout << "\t (from Pdg formula : "
346  << std::setw(8) << G4BestUnit(rMPdg,"Length") << ")";
347  }
348 }
G4Material::GetZ
G4double GetZ() const
Definition: G4Material.cc:623
G4BestUnit
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
Definition: G4SteppingVerbose.cc:53
kStateGas
Definition: G4Material.hh:114
G4int
int G4int
Definition: G4Types.hh:78
eplot.material
string material
Definition: eplot.py:19
G4cout
G4GLOB_DLL std::ostream G4cout
G4EmCalculator::ComputeDEDX
G4double ComputeDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=DBL_MAX)
Definition: G4EmCalculator.cc:381
G4Material::GetRadlen
G4double GetRadlen() const
Definition: G4Material.hh:220
G4Electron::Electron
static G4Electron * Electron()
Definition: G4Electron.cc:94
MeV
static constexpr double MeV
Definition: G4SIunits.hh:214
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:186
G4Material::GetState
G4State GetState() const
Definition: G4Material.hh:181
G4double
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

Here is the caller graph for this function:

void RunAction::CriticalEnergy ( )
virtual void RunAction::EndOfRunAction ( const G4Run *  arun)
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  arun)
virtual

Reimplemented from G4UserRunAction.

void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  run)
virtual

Reimplemented from G4UserRunAction.

void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::EndOfRunAction ( const G4Run *  run)
virtual

Reimplemented from G4UserRunAction.

Definition at line 260 of file RunAction.cc.

261 { }
virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

virtual void RunAction::EndOfRunAction ( const G4Run *  )
virtual

Reimplemented from G4UserRunAction.

void RunAction::EventFinished ( )

Definition at line 246 of file RunAction.cc.

247 {
248  ++fTotalEventCount;
249  fPhotonStats.EventFinished();
250  fElectronStats.EventFinished();
251  fPositronStats.EventFinished();
252 }
void RunAction::FillData ( const G4ParticleDefinition *  particle,
G4double  kinEnergy,
G4double  costheta,
G4double  phi,
G4double  longitudinalPolarization 
)

Definition at line 94 of file RunAction.cc.

98 {
99  G4int id = -1;
100  if (particle == fGamma) {
101  fPhotonStats.FillData(kinEnergy, costheta, longitudinalPolarization);
102  if(fAnalysisManager) { id = 1; }
103  } else if (particle == fElectron) {
104  fElectronStats.FillData(kinEnergy, costheta, longitudinalPolarization);
105  if(fAnalysisManager) { id = 5; }
106  } else if (particle == fPositron) {
107  fPositronStats.FillData(kinEnergy, costheta, longitudinalPolarization);
108  if(fAnalysisManager) { id = 9; }
109  }
110  if(id > 0) {
111  fAnalysisManager->FillH1(id,kinEnergy,1.0);
112  fAnalysisManager->FillH1(id+1,costheta,1.0);
113  fAnalysisManager->FillH1(id+2,phi,1.0);
114  fAnalysisManager->FillH1(id+3,longitudinalPolarization,1.0);
115  }
116 }
G4int
int G4int
Definition: G4Types.hh:78
void RunAction::FillEdep ( G4double  de,
G4double  eni 
)
inline

Definition at line 89 of file RunAction.hh.

90 {
91  fEdeptot += de;
92  fEniel += eni;
93 }
void RunAction::FillHisto ( G4int  id,
G4double  x,
G4double  weight = 1.0 
)

Definition at line 186 of file RunAction.cc.

187 {
188  if(GetVerbose() > 1) {
189  G4cout << "FillHisto " << ih << " x=" << x << " weight= " << weight
190  << G4endl;
191  }
192  fHisto->Fill(ih, x, weight);
193 }
test.x
tuple x
Definition: test.py:50
Histo::Fill
void Fill(G4int, G4double, G4double)
Definition: Histo.cc:221
G4cout
G4GLOB_DLL std::ostream G4cout
RunAction::GetVerbose
G4int GetVerbose()
Definition: RunAction.hh:65
G4endl
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

void RunAction::FillPerEvent ( G4double  ,
G4double  ,
G4double  ,
G4double   
)
void RunAction::FillPerEvent ( G4double  ,
G4double  ,
G4double  ,
G4double   
)
void RunAction::FillPerEvent ( G4double  ,
G4double  ,
G4double  ,
G4double   
)
void RunAction::FillPerEvent ( G4double  EAbs,
G4double  EGap,
G4double  LAbs,
G4double  LGap 
)

Definition at line 76 of file RunAction.cc.

78 {
79  //accumulate statistic
80  //
81  fSumEAbs += EAbs; fSum2EAbs += EAbs*EAbs;
82  fSumEGap += EGap; fSum2EGap += EGap*EGap;
83 
84  fSumLAbs += LAbs; fSum2LAbs += LAbs*LAbs;
85  fSumLGap += LGap; fSum2LGap += LGap*LGap;
86 }
void RunAction::FillTallyEdep ( G4int  n,
G4double  e 
)
inline

Definition at line 82 of file RunAction.hh.

83 {
84  fTallyEdep[n] += e;
85 }
n
const G4int n
Definition: G4UrQMD1_3Interface.hh:144
virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

G4Run * RunAction::GenerateRun ( void  )
virtual

Reimplemented from G4UserRunAction.

Definition at line 66 of file RunAction.cc.

67 {
68  fRun = new Run(fDetector);
69  return fRun;
70 }
Run
Definition: Run.hh:46
virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

virtual G4Run* RunAction::GenerateRun ( )
virtual

Reimplemented from G4UserRunAction.

void RunAction::GetCuts ( )

Definition at line 269 of file RunAction.cc.

270 {
271  G4ProductionCutsTable* theCoupleTable =
272  G4ProductionCutsTable::GetProductionCutsTable();
273 
274  size_t numOfCouples = theCoupleTable->GetTableSize();
275  const G4MaterialCutsCouple* couple = 0;
276  G4int index = 0;
277  for (size_t i=0; i<numOfCouples; i++) {
278  couple = theCoupleTable->GetMaterialCutsCouple(i);
279  if (couple->GetMaterial() == detector->GetMaterial()) {index = i; break;}
280  }
281 
282  rangeCut[0] =
283  (*(theCoupleTable->GetRangeCutsVector(idxG4GammaCut)))[index];
284  rangeCut[1] =
285  (*(theCoupleTable->GetRangeCutsVector(idxG4ElectronCut)))[index];
286  rangeCut[2] =
287  (*(theCoupleTable->GetRangeCutsVector(idxG4PositronCut)))[index];
288 
289  energyCut[0] =
290  (*(theCoupleTable->GetEnergyCutsVector(idxG4GammaCut)))[index];
291  energyCut[1] =
292  (*(theCoupleTable->GetEnergyCutsVector(idxG4ElectronCut)))[index];
293  energyCut[2] =
294  (*(theCoupleTable->GetEnergyCutsVector(idxG4PositronCut)))[index];
295 
296 }
G4ProductionCutsTable::GetEnergyCutsVector
const std::vector< G4double > * GetEnergyCutsVector(size_t pcIdx) const
Definition: G4ProductionCutsTable.hh:250
idxG4GammaCut
Definition: G4ProductionCuts.hh:53
G4int
int G4int
Definition: G4Types.hh:78
idxG4ElectronCut
Definition: G4ProductionCuts.hh:54
idxG4PositronCut
Definition: G4ProductionCuts.hh:55
G4ProductionCutsTable::GetRangeCutsVector
const std::vector< G4double > * GetRangeCutsVector(size_t pcIdx) const
Definition: G4ProductionCutsTable.hh:244
G4ProductionCutsTable::GetProductionCutsTable
static G4ProductionCutsTable * GetProductionCutsTable()
Definition: G4ProductionCutsTable.cc:63
G4ProductionCutsTable::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
Definition: G4ProductionCutsTable.hh:262
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:150

Here is the call graph for this function:

Here is the caller graph for this function:

void RunAction::GetCuts ( )
G4double RunAction::GetDoseBox ( G4int  i)
inline

Definition at line 81 of file RunAction.hh.

81 { return fDose3DDose[i];}
G4double RunAction::GetDoseC ( )
inline

Definition at line 66 of file RunAction.hh.

66 {return fDoseC;}
G4double RunAction::GetDoseN ( )
inline

Definition at line 62 of file RunAction.hh.

62 {return fDoseN;}
G4double RunAction::GetEnergyFromCSDARange ( G4double  range,
G4ParticleDefinition *  particle,
G4Material *  material,
G4double  Etry 
)

Definition at line 258 of file RunAction.cc.

260 {
261  G4EmCalculator emCal;
262 
263  G4double Energy = Etry, dE = 0., dEdx;
264  G4double r, dr;
265  G4double err = 1., errmax = 0.00001;
266  G4int iter = 0 , itermax = 10;
267  while (err > errmax && iter < itermax) {
268  iter++;
269  Energy -= dE;
270  r = emCal.GetCSDARange(Energy,particle,material);
271  dr = r - range;
272  dEdx = emCal.ComputeTotalDEDX(Energy,particle,material);
273  dE = dEdx*dr;
274  err = std::abs(dE)/Energy;
275  }
276  if (iter == itermax) {
277  G4cout
278  << "\n ---> warning: RunAction::GetEnergyFromCSDARange() did not converge"
279  << " Etry = " << G4BestUnit(Etry,"Energy")
280  << " Energy = " << G4BestUnit(Energy,"Energy")
281  << " err = " << err
282  << " iter = " << iter << G4endl;
283  }
284 
285  return Energy;
286 }
G4EmCalculator::ComputeTotalDEDX
G4double ComputeTotalDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, G4double cut=DBL_MAX)
Definition: G4EmCalculator.cc:575
G4EmCalculator::GetCSDARange
G4double GetCSDARange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:204
G4BestUnit
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
Definition: G4SteppingVerbose.cc:53
G4int
int G4int
Definition: G4Types.hh:78
dE
static const G4double dE
Definition: G4PhotoNuclearCrossSection.cc:47
G4cout
G4GLOB_DLL std::ostream G4cout
pyG4EmCalculator::range
const G4ParticleDefinition const G4Material *G4double range
Definition: pyG4EmCalculator.cc:185
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4double RunAction::GetEnergyFromRestrictedRange ( G4double  range,
G4ParticleDefinition *  particle,
G4Material *  material,
G4double  Etry 
)

Definition at line 226 of file RunAction.cc.

228 {
229  G4EmCalculator emCal;
230 
231  G4double Energy = Etry, dE = 0., dEdx;
232  G4double r, dr;
233  G4double err = 1., errmax = 0.00001;
234  G4int iter = 0 , itermax = 10;
235  while (err > errmax && iter < itermax) {
236  iter++;
237  Energy -= dE;
238  r = emCal.GetRangeFromRestricteDEDX(Energy,particle,material);
239  dr = r - range;
240  dEdx = emCal.GetDEDX(Energy,particle,material);
241  dE = dEdx*dr;
242  err = std::abs(dE)/Energy;
243  }
244  if (iter == itermax) {
245  G4cout
246  << "\n ---> warning: RunAction::GetEnergyFromRestRange() did not converge"
247  << " Etry = " << G4BestUnit(Etry,"Energy")
248  << " Energy = " << G4BestUnit(Energy,"Energy")
249  << " err = " << err
250  << " iter = " << iter << G4endl;
251  }
252 
253  return Energy;
254 }
G4BestUnit
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
Definition: G4SteppingVerbose.cc:53
G4int
int G4int
Definition: G4Types.hh:78
dE
static const G4double dE
Definition: G4PhotoNuclearCrossSection.cc:47
G4cout
G4GLOB_DLL std::ostream G4cout
pyG4EmCalculator::range
const G4ParticleDefinition const G4Material *G4double range
Definition: pyG4EmCalculator.cc:185
G4EmCalculator::GetRangeFromRestricteDEDX
G4double GetRangeFromRestricteDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:181
G4EmCalculator::GetDEDX
G4double GetDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
Definition: G4EmCalculator.cc:140
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4double RunAction::GetMassCytoplasm ( )
inline

Definition at line 78 of file RunAction.hh.

78 {return fMassCytoplasm;}
G4double RunAction::GetMassNucleus ( )
inline

Definition at line 75 of file RunAction.hh.

75 {return fMassNucleus;}
G4int RunAction::GetNbOfHitsGas ( )
inline

Definition at line 71 of file RunAction.hh.

71 {return fNbOfHitsGas;}
G4int RunAction::GetNumEvent ( )
inline

Definition at line 68 of file RunAction.hh.

68 {return fNumEvent;}
G4double RunAction::GetOffsetX ( )
inline

Definition at line 62 of file RunAction.hh.

62 { return fOffsetX;}
G4int RunAction::GetRndmFreq ( )
inline

Definition at line 58 of file RunAction.hh.

58 {return fSaveRndm;}
G4int RunAction::GetRndmFreq ( )
inline

Definition at line 60 of file RunAction.hh.

60 { return fRndmFreq; }
G4int RunAction::GetRow ( )
inline

Definition at line 52 of file RunAction.hh.

52 {return fRow;}
G4ThreeVector RunAction::GetVectCell ( G4int  i)
inline

Definition at line 83 of file RunAction.hh.

83 {return fMapVoxels[i];}
G4int RunAction::GetVerbose ( )
inline

Definition at line 65 of file RunAction.hh.

65 { return fVerboseLevel;}

Here is the caller graph for this function:

void RunAction::SetApplyLimit ( G4bool  val)

Definition at line 132 of file RunAction.cc.

133 {
134  if (fRun) fRun->SetApplyLimit(val);
135 }
Run::SetApplyLimit
void SetApplyLimit(G4bool)
Definition: Run.cc:357

Here is the call graph for this function:

void RunAction::SetBinSize ( G4double  size)
inline

Definition at line 61 of file RunAction.hh.

61 { fBinLength = size; }
void RunAction::SetDoseC ( G4double  dose)
inline

Definition at line 65 of file RunAction.hh.

65 { fDoseC = dose;}
void RunAction::SetDoseN ( G4double  dose)
inline

Definition at line 61 of file RunAction.hh.

61 { fDoseN = dose;}
void RunAction::SetEdepAndRMS ( G4int  i,
G4double  edep,
G4double  rms,
G4double  lim 
)

Definition at line 125 of file RunAction.cc.

126 {
127  if (fRun) fRun->SetEdepAndRMS(i, edep, rms, lim);
128 }
Run::SetEdepAndRMS
void SetEdepAndRMS(G4int, G4double, G4double, G4double)
Definition: Run.cc:346

Here is the call graph for this function:

void RunAction::SetEdepAndRMS ( G4ThreeVector  Value)

Definition at line 185 of file RunAction.cc.

186 {
187  fEdeptrue = Value(0);
188  fRmstrue = Value(1);
189  fLimittrue= Value(2);
190 }

Here is the caller graph for this function:

void RunAction::SetHistoName ( G4String &  val)
inline

Definition at line 56 of file RunAction.hh.

57  {
58  fFileName = val;
59  };
void RunAction::SetHistoName ( G4String &  val)

Definition at line 141 of file RunAction.cc.

142 {
143  fHistoName[0] = val;
144 }

Here is the caller graph for this function:

void RunAction::SetMassCytoplasm ( G4double  mC)
inline

Definition at line 77 of file RunAction.hh.

77 { fMassCytoplasm = mC;}
void RunAction::SetMassNucleus ( G4double  mN)
inline

Definition at line 74 of file RunAction.hh.

74 { fMassNucleus = mN;}
void RunAction::SetNumEvent ( G4int  i)
inline

Definition at line 69 of file RunAction.hh.

69 {fNumEvent = i;}
void RunAction::SetPrintFlag ( G4bool  flag)

Definition at line 119 of file RunAction.cc.

120 {
121  fPrint = flag;
122 }
void RunAction::SetRndmFreq ( G4int  val)
inline

Definition at line 57 of file RunAction.hh.

57 {fSaveRndm = val;}

Here is the caller graph for this function:

void RunAction::SetRndmFreq ( G4int  val)
inline

Definition at line 59 of file RunAction.hh.

59 { fRndmFreq = val; }
void RunAction::SetVerbose ( G4int  val)

Definition at line 193 of file RunAction.cc.

194 {
195  fVerbose = val;
196  if (fRun) fRun->SetVerbose(val);
197 }
Run::SetVerbose
void SetVerbose(G4int val)
Definition: Run.hh:70

Here is the call graph for this function:

Here is the caller graph for this function:

void RunAction::SetVerbose ( G4int  verbose)
inline

Definition at line 64 of file RunAction.hh.

64 { fVerboseLevel = verbose;}
void RunAction::SumLateralDisplacement ( G4double  displa)
inline

Definition at line 67 of file RunAction.hh.

68  {fLDispl += displa; fLDispl2 += displa*displa;}
void RunAction::SumPathLength ( G4double  truepl,
G4double  geompl 
)
inline

Definition at line 61 of file RunAction.hh.

62  {fTotalCount++;
63  fTruePL += truepl; fTruePL2 += truepl*truepl;
64  fGeomPL += geompl; fGeomPL2 += geompl*geompl;
65  };
void RunAction::SumPhiCorrel ( G4double  correl)
inline

Definition at line 76 of file RunAction.hh.

77  {fPhiCor += correl; fPhiCor2 += correl*correl;}
void RunAction::SumPsi ( G4double  psi)
inline

Definition at line 70 of file RunAction.hh.

71  {fPsiSpa += psi; fPsiSpa2 += psi*psi;}
void RunAction::SumTetaPlane ( G4double  teta)
inline

Definition at line 73 of file RunAction.hh.

74  {fTetPrj += teta; fTetPrj2 += teta*teta;}
void RunAction::WriteHistogram ( )
The documentation for this class was generated from the following files: