Geant4  10.02.p03
Em10RunAction Class Reference

#include <Em10RunAction.hh>

Inheritance diagram for Em10RunAction:
Collaboration diagram for Em10RunAction:

Public Member Functions

 Em10RunAction ()
 
 ~Em10RunAction ()
 
void BeginOfRunAction (const G4Run *)
 
void EndOfRunAction (const G4Run *)
 
void CountEvent ()
 
void CountParticles (G4double, G4double)
 
void AddEP (G4double, G4double)
 
void AddEdeps (G4double Eabs)
 
void AddTrackLength (G4double tlabs)
 
void AddnStepsCharged (G4double ns)
 
void AddnStepsNeutral (G4double ns)
 
void AddTrRef (G4double tr, G4double ref)
 
void FillEn (G4double En)
 
void FillTh (G4double Th)
 
void FillThBack (G4double Th)
 
void FillR (G4double R)
 
void FillTt (G4double Tt)
 
void FillTb (G4double Tt)
 
void FillTsec (G4double T)
 
void FillGammaSpectrum (G4double E)
 
void FillNbOfSteps (G4double nstep)
 
void Fillvertexz (G4double z)
 
void SethistName (G4String name)
 
void SetnbinStep (G4int nbin)
 
void SetSteplow (G4double Slow)
 
void SetStephigh (G4double Shigh)
 
void SetnbinEn (G4int nbin)
 
void SetEnlow (G4double Elow)
 
void SetEnhigh (G4double Enhigh)
 
void SetnbinTt (G4int nbin)
 
void SetTtlow (G4double Ttlow)
 
void SetTthigh (G4double Tthigh)
 
void SetnbinTb (G4int nbin)
 
void SetTblow (G4double Tblow)
 
void SetTbhigh (G4double Tbhigh)
 
void SetnbinTsec (G4int nbin)
 
void SetTseclow (G4double Tlow)
 
void SetTsechigh (G4double Thigh)
 
void SetnbinTh (G4int nbin)
 
void SetThlow (G4double Thlow)
 
void SetThhigh (G4double Thhigh)
 
void SetnbinThBack (G4int nbin)
 
void SetThlowBack (G4double Thlow)
 
void SetThhighBack (G4double Thhigh)
 
void SetnbinR (G4int nbin)
 
void SetRlow (G4double Rlow)
 
void SetRhigh (G4double Rhigh)
 
void SetnbinGamma (G4int nbin)
 
void SetElowGamma (G4double Elow)
 
void SetEhighGamma (G4double Ehigh)
 
void Setnbinzvertex (G4int nbin)
 
void Setzlow (G4double z)
 
void Setzhigh (G4double z)
 
void SetRndmFreq (G4int val)
 
G4int GetRndmFreq ()
 
- Public Member Functions inherited from G4UserRunAction
 G4UserRunAction ()
 
virtual ~G4UserRunAction ()
 
virtual G4RunGenerateRun ()
 
void SetMaster (G4bool val=true)
 
G4bool IsMaster () const
 

Private Member Functions

void bookHisto ()
 

Private Attributes

G4String histName
 
G4double EnergySumAbs
 
G4double EnergySquareSumAbs
 
G4double tlSumAbs
 
G4double tlsquareSumAbs
 
G4double nStepSumCharged
 
G4double nStepSum2Charged
 
G4double nStepSumNeutral
 
G4double nStepSum2Neutral
 
G4double TotNbofEvents
 
G4double SumCharged
 
G4double Sum2Charged
 
G4double SumNeutral
 
G4double Sum2Neutral
 
G4double Selectron
 
G4double Spositron
 
G4double Transmitted
 
G4double Reflected
 
G4double entryStep
 
G4double underStep
 
G4double overStep
 
G4double distStep [200]
 
G4double Steplow
 
G4double Stephigh
 
G4double dStep
 
G4int nbinStep
 
G4double entryEn
 
G4double underEn
 
G4double overEn
 
G4double distEn [200]
 
G4double Enlow
 
G4double Enhigh
 
G4double dEn
 
G4int nbinEn
 
G4double entryTt
 
G4double underTt
 
G4double overTt
 
G4double distTt [200]
 
G4double Ttlow
 
G4double Tthigh
 
G4double dTt
 
G4int nbinTt
 
G4double Ttmean
 
G4double Tt2mean
 
G4double entryTb
 
G4double underTb
 
G4double overTb
 
G4double distTb [200]
 
G4double Tblow
 
G4double Tbhigh
 
G4double dTb
 
G4int nbinTb
 
G4double Tbmean
 
G4double Tb2mean
 
G4double entryTsec
 
G4double underTsec
 
G4double overTsec
 
G4double distTsec [200]
 
G4double Tseclow
 
G4double Tsechigh
 
G4double dTsec
 
G4int nbinTsec
 
G4double entryTh
 
G4double underTh
 
G4double overTh
 
G4double distTh [200]
 
G4double Thlow
 
G4double Thhigh
 
G4double dTh
 
G4int nbinTh
 
G4double entryThback
 
G4double underThback
 
G4double overThback
 
G4double distThback [200]
 
G4double Thlowback
 
G4double Thhighback
 
G4double dThback
 
G4int nbinThback
 
G4double entryR
 
G4double underR
 
G4double overR
 
G4double distR [200]
 
G4double Rlow
 
G4double Rhigh
 
G4double dR
 
G4int nbinR
 
G4double Rmean
 
G4double R2mean
 
G4double entryGamma
 
G4double underGamma
 
G4double overGamma
 
G4double distGamma [200]
 
G4double ElowGamma
 
G4double EhighGamma
 
G4double dEGamma
 
G4int nbinGamma
 
G4double entryvertexz
 
G4double undervertexz
 
G4double oververtexz
 
G4double distvertexz [200]
 
G4double zlow
 
G4double zhigh
 
G4double dz
 
G4int nbinvertexz
 
Em10RunMessengerrunMessenger
 
G4int saveRndm
 

Additional Inherited Members

- Protected Attributes inherited from G4UserRunAction
G4bool isMaster
 

Detailed Description

Definition at line 51 of file Em10RunAction.hh.

Constructor & Destructor Documentation

◆ Em10RunAction()

Em10RunAction::Em10RunAction ( )

Definition at line 50 of file Em10RunAction.cc.

51  :G4UserRunAction(),
52  histName("histfile"),nbinStep(0),nbinEn(0),nbinTt(0),nbinTb(0),
53  nbinTsec(0),nbinTh(0),nbinThback(0),nbinR(0),nbinGamma(0),
54  nbinvertexz(0)
55 {
56  runMessenger = new Em10RunMessenger(this);
57  saveRndm = 1;
58 
68  oververtexz=0.;
69 }
G4double nStepSumCharged
G4double Sum2Charged
G4double entryThback
G4double underTh
G4double TotNbofEvents
G4double Transmitted
G4double dEGamma
G4double entryGamma
G4double underEn
G4double entryTb
G4double overTsec
G4double dThback
G4double overThback
G4double SumCharged
G4double entryEn
G4double overGamma
G4double Selectron
G4double nStepSum2Neutral
G4double underStep
G4double undervertexz
G4double tlSumAbs
G4double Sum2Neutral
G4double nStepSum2Charged
G4String histName
G4double oververtexz
G4double tlsquareSumAbs
G4double underTsec
G4double entryvertexz
G4double nStepSumNeutral
G4double entryTsec
G4double underThback
G4double entryTh
G4double entryTt
G4double Reflected
G4double SumNeutral
G4double Tt2mean
G4double Spositron
G4double EnergySumAbs
G4double underTt
Em10RunMessenger * runMessenger
G4double entryStep
G4double underGamma
G4double Tb2mean
G4double overStep
G4double EnergySquareSumAbs
G4double underTb

◆ ~Em10RunAction()

Em10RunAction::~Em10RunAction ( )

Definition at line 73 of file Em10RunAction.cc.

74 {
75  delete runMessenger;
76 }
Em10RunMessenger * runMessenger

Member Function Documentation

◆ AddEdeps()

void Em10RunAction::AddEdeps ( G4double  Eabs)

Definition at line 716 of file Em10RunAction.cc.

717 {
718  EnergySumAbs += Eabs;
719  EnergySquareSumAbs += Eabs*Eabs;
720 }
G4double EnergySumAbs
G4double EnergySquareSumAbs
Here is the caller graph for this function:

◆ AddEP()

void Em10RunAction::AddEP ( G4double  nele,
G4double  npos 
)

Definition at line 1051 of file Em10RunAction.cc.

1052 {
1053  Selectron += nele;
1054  Spositron += npos;
1055 }
G4double Selectron
G4double Spositron
Here is the caller graph for this function:

◆ AddnStepsCharged()

void Em10RunAction::AddnStepsCharged ( G4double  ns)

Definition at line 700 of file Em10RunAction.cc.

701 {
702  nStepSumCharged += nstp;
703  nStepSum2Charged += nstp*nstp;
704 }
G4double nStepSumCharged
G4double nStepSum2Charged
Here is the caller graph for this function:

◆ AddnStepsNeutral()

void Em10RunAction::AddnStepsNeutral ( G4double  ns)

Definition at line 708 of file Em10RunAction.cc.

709 {
710  nStepSumNeutral += nstp;
711  nStepSum2Neutral += nstp*nstp;
712 }
G4double nStepSum2Neutral
G4double nStepSumNeutral
Here is the caller graph for this function:

◆ AddTrackLength()

void Em10RunAction::AddTrackLength ( G4double  tlabs)

Definition at line 724 of file Em10RunAction.cc.

725 {
726  tlSumAbs += tlabs;
727  tlsquareSumAbs += tlabs*tlabs ;
728 }
G4double tlSumAbs
G4double tlsquareSumAbs
Here is the caller graph for this function:

◆ AddTrRef()

void Em10RunAction::AddTrRef ( G4double  tr,
G4double  ref 
)

Definition at line 732 of file Em10RunAction.cc.

733 {
734  Transmitted += tr ;
735  Reflected += ref;
736 }
G4double Transmitted
G4double Reflected
Here is the caller graph for this function:

◆ BeginOfRunAction()

void Em10RunAction::BeginOfRunAction ( const G4Run aRun)
virtual

Reimplemented from G4UserRunAction.

Definition at line 86 of file Em10RunAction.cc.

87 {
88  G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
89 
90  // save Rndm status
91  if (saveRndm > 0)
92  {
94  CLHEP::HepRandom::saveEngineStatus("beginOfRun.rndm");
95  }
97 
99 
100  if(pVVisManager) UI->ApplyCommand("/vis/scene/notifyHandlers");
101 
102 
103  EnergySumAbs = 0. ;
104  EnergySquareSumAbs = 0.;
105  tlSumAbs = 0. ;
106  tlsquareSumAbs = 0. ;
107  nStepSumCharged = 0. ;
108  nStepSum2Charged= 0. ;
109  nStepSumNeutral = 0. ;
110  nStepSum2Neutral= 0. ;
111  TotNbofEvents = 0. ;
112  SumCharged=0.;
113  SumNeutral=0.;
114  Sum2Charged=0.;
115  Sum2Neutral=0.;
116  Selectron=0.;
117  Spositron=0.;
118 
119  Transmitted=0.;
120  Reflected =0.;
121 
122  // plot definitions
123 
124  if(nbinStep>0)
125  {
127  entryStep=0.;
128  underStep=0.;
129  overStep=0.;
130  for(G4int ist=0; ist<200; ist++)
131  {
132  distStep[ist]=0.;
133  }
134  }
135  if(nbinEn>0)
136  {
137  dEn = (Enhigh-Enlow)/nbinEn ;
138  entryEn=0.;
139  underEn=0.;
140  overEn=0.;
141 
142  for (G4int ien=0; ien<200; ien++) distEn[ien]=0.;
143  }
144  if(nbinTt>0)
145  {
146  dTt = (Tthigh-Ttlow)/nbinTt ;
147  entryTt=0.;
148  underTt=0.;
149  overTt=0.;
150 
151  for (G4int itt=0; itt<200; itt++) distTt[itt]=0.;
152 
153  Ttmean=0.;
154  Tt2mean=0.;
155  }
156  if(nbinTb>0)
157  {
158  dTb = (Tbhigh-Tblow)/nbinTb ;
159  entryTb=0.;
160  underTb=0.;
161  overTb=0.;
162  for (G4int itt=0; itt<200; itt++)
163  {
164  distTb[itt]=0.;
165  }
166  Tbmean=0.;
167  Tb2mean=0.;
168  }
169  if(nbinTsec>0)
170  {
172  entryTsec=0.;
173  underTsec=0.;
174  overTsec=0.;
175  for (G4int its=0; its<200; its++)
176  {
177  distTsec[its]=0.;
178  }
179  }
180  if(nbinTh>0)
181  {
182  dTh = (Thhigh-Thlow)/nbinTh ;
183  entryTh=0.;
184  underTh=0.;
185  overTh=0.;
186  for (G4int ith=0; ith<200; ith++)
187  {
188  distTh[ith]=0.;
189  }
190  }
191 
192  if(nbinThback>0)
193  {
195  entryThback=0.;
196  underThback=0.;
197  overThback=0.;
198  for (G4int ithback=0; ithback<200; ithback++)
199  {
200  distThback[ithback]=0.;
201  }
202  }
203 
204 
205  if(nbinR >0)
206  {
207  dR = (Rhigh-Rlow)/nbinR ;
208  entryR =0.;
209  underR =0.;
210  overR =0.;
211  for (G4int ir =0; ir<200; ir++)
212  {
213  distR[ir]=0.;
214  }
215  Rmean=0.;
216  R2mean=0.;
217  }
218 
219  if(nbinGamma>0)
220  {
221  dEGamma = std::log(EhighGamma/ElowGamma)/nbinGamma ;
222  entryGamma = 0.;
223  underGamma=0.;
224  overGamma=0.;
225  for (G4int ig=0; ig<200; ig++)
226  {
227  distGamma[ig]=0.;
228  }
229  }
230  if(nbinvertexz>0)
231  {
233  entryvertexz=0.;
234  undervertexz=0.;
235  oververtexz=0.;
236  for(G4int iz=0; iz<200; iz++)
237  {
238  distvertexz[iz]=0.;
239  }
240  }
241 
242  bookHisto();
243 }
G4double nStepSumCharged
G4double distvertexz[200]
G4double Sum2Charged
G4double entryThback
G4double underTh
G4double TotNbofEvents
G4double Transmitted
G4double dEGamma
G4double entryGamma
G4double distTsec[200]
G4double distTh[200]
static G4VVisManager * GetConcreteInstance()
G4double underEn
G4double entryTb
G4double overTsec
G4double distTt[200]
G4double Thlowback
G4double dThback
G4double distEn[200]
G4double distStep[200]
G4double overThback
G4double Stephigh
G4double SumCharged
int G4int
Definition: G4Types.hh:78
G4double entryEn
static G4UImanager * GetUIpointer()
Definition: G4UImanager.cc:58
G4double overGamma
G4double distGamma[200]
G4GLOB_DLL std::ostream G4cout
G4double Selectron
G4double nStepSum2Neutral
G4double distThback[200]
G4double underStep
G4double EhighGamma
G4double undervertexz
G4double tlSumAbs
G4double Sum2Neutral
G4double nStepSum2Charged
G4double distTb[200]
G4double ElowGamma
G4double iz
Definition: TRTMaterials.hh:39
G4double distR[200]
G4double oververtexz
G4double tlsquareSumAbs
G4double underTsec
G4double entryvertexz
G4int GetRunID() const
Definition: G4Run.hh:76
static void showEngineStatus()
Definition: Random.cc:204
G4double Steplow
G4double Thhighback
G4double nStepSumNeutral
G4double entryTsec
G4double underThback
G4double entryTh
G4double entryTt
G4double Reflected
G4double SumNeutral
G4double Tseclow
G4double Tt2mean
G4double Spositron
G4double EnergySumAbs
G4double underTt
static void saveEngineStatus(const char filename[]="Config.conf")
Definition: Random.cc:176
G4double entryStep
#define G4endl
Definition: G4ios.hh:61
G4double underGamma
G4double Tsechigh
G4double Tb2mean
G4double overStep
G4double EnergySquareSumAbs
G4int ApplyCommand(const char *aCommand)
Definition: G4UImanager.cc:446
G4double underTb
Here is the call graph for this function:

◆ bookHisto()

void Em10RunAction::bookHisto ( )
private

Definition at line 80 of file Em10RunAction.cc.

81 {
82 }
Here is the caller graph for this function:

◆ CountEvent()

void Em10RunAction::CountEvent ( )

Definition at line 693 of file Em10RunAction.cc.

694 {
695  TotNbofEvents += 1. ;
696 }
G4double TotNbofEvents
Here is the caller graph for this function:

◆ CountParticles()

void Em10RunAction::CountParticles ( G4double  nch,
G4double  nne 
)

Definition at line 1043 of file Em10RunAction.cc.

1044 {
1045  SumCharged += nch ;
1046  SumNeutral += nne ;
1047  Sum2Charged += nch*nch ;
1048  Sum2Neutral += nne*nne ;
1049 }
G4double Sum2Charged
G4double SumCharged
G4double Sum2Neutral
G4double SumNeutral
Here is the caller graph for this function:

◆ EndOfRunAction()

void Em10RunAction::EndOfRunAction ( const G4Run )
virtual

Reimplemented from G4UserRunAction.

Definition at line 247 of file Em10RunAction.cc.

248 {
249  G4double sAbs,sigAbs,sigstep,sigcharged,signeutral;
250  if(0.0 >= TotNbofEvents) { return; }
251 
254  if(sAbs>0.)
255  sAbs = std::sqrt(sAbs/TotNbofEvents) ;
256  else
257  sAbs = 0. ;
258 
261  if(sigAbs>0.)
262  sigAbs = std::sqrt(sigAbs/TotNbofEvents);
263  else
264  sigAbs = 0.;
265 
268  if(sigstep>0.)
269  sigstep = std::sqrt(sigstep/TotNbofEvents);
270  else
271  sigstep = 0.;
272  G4double sigch=sigstep ;
273 
276  if(sigstep>0.)
277  sigstep = std::sqrt(sigstep/TotNbofEvents);
278  else
279  sigstep = 0.;
280  G4double signe=sigstep ;
281 
284  if(sigcharged>0.)
285  sigcharged = std::sqrt(sigcharged/TotNbofEvents);
286  else
287  sigcharged = 0. ;
288 
291  if(signeutral>0.)
292  signeutral = std::sqrt(signeutral/TotNbofEvents);
293  else
294  signeutral = 0. ;
295 
298 
301  G4cout << " ================== run summary =====================" << G4endl;
302  G4int prec = G4cout.precision(6);
303  G4cout << " end of Run TotNbofEvents = " <<
304  TotNbofEvents << G4endl ;
305  G4cout << " mean charged track length in absorber=" <<
306  tlSumAbs/mm << " +- " << sAbs/mm <<
307  " mm " << G4endl;
308  G4cout << G4endl;
309  G4cout << " mean energy deposit in absorber=" <<
310  EnergySumAbs/MeV << " +- " << sigAbs/MeV <<
311  " MeV " << G4endl ;
312  G4cout << G4endl ;
313  G4cout << " mean number of steps in absorber (charged) =" <<
314  nStepSumCharged << " +- " << sigch <<
315  " " << G4endl ;
316  G4cout << " mean number of steps in absorber (neutral) =" <<
317  nStepSumNeutral << " +- " << signe <<
318  " " << G4endl ;
319  G4cout << G4endl ;
320  G4cout << " mean number of charged secondaries = " <<
321  SumCharged << " +- " << sigcharged << G4endl;
322  G4cout << G4endl ;
323  G4cout << " mean number of neutral secondaries = " <<
324  SumNeutral << " +- " << signeutral << G4endl;
325  G4cout << G4endl ;
326 
327  G4cout << " mean number of e-s =" << Selectron <<
328  " and e+s =" << Spositron << G4endl;
329  G4cout << G4endl;
330 
331  G4cout << "(number) transmission coeff=" << Transmitted <<
332  " reflection coeff=" << Reflected << G4endl;
333  G4cout << G4endl;
334 
335  if(nbinStep>0)
336  {G4double E , dnorm, norm ;
337  G4cout << " step number/event distribution " << G4endl ;
338  G4cout << "#entries=" << entryStep << " #underflows=" << underStep <<
339  " #overflows=" << overStep << G4endl ;
340  if( entryStep>0.)
341  {
342  E = Steplow - dStep ;
343  norm = TotNbofEvents ;
344  G4cout << " bin nb nsteplow entries normalized " << G4endl ;
345  for(G4int iss=0; iss<nbinStep; iss++)
346  {
347  E += dStep ;
348  dnorm = distStep[iss]/norm;
349  G4cout << std::setw(5) << iss << std::setw(10) << E <<
350  std::setw(12) << distStep[iss] <<
351  std::setw(12) << dnorm << G4endl ;
352  }
353  G4cout << G4endl;
354  }
355  }
356  if(nbinEn > 0)
357  {
358  std::ofstream fileOut("distribution.out", std::ios::out ) ;
359  fileOut.setf( std::ios::scientific, std::ios::floatfield );
360 
361  std::ofstream normOut("normDist.out", std::ios::out ) ;
362  normOut.setf( std::ios::scientific, std::ios::floatfield );
363 
364  G4double E , dnorm, norm,fmax,Emp ;
365  Emp=-999.999 ;
366  G4cout << " energy deposit distribution " << G4endl ;
367  G4cout << "#entries=" << entryEn << " #underflows=" << underEn <<
368  " #overflows=" << overEn << G4endl ;
369  if( entryEn>0.)
370  {
371  E = Enlow - dEn ;
372  norm = TotNbofEvents*1.0 ; // *dEn ;
373  G4cout << " bin nb Elow entries normalized " << G4endl ;
374  fmax = 0. ;
375 
376  for(G4int ien=0; ien<nbinEn; ien++)
377  {
378  E += dEn ;
379 
380  if(distEn[ien]>fmax)
381  {
382  fmax = distEn[ien] ;
383  Emp = E ; // most probable roughly
384  }
385  dnorm = distEn[ien]/norm;
386 
387  G4cout << std::setw(5) << ien << std::setw(10) << E/keV <<
388  std::setw(12) << distEn[ien] <<
389  std::setw(12) << dnorm << G4endl ;
390 
391  fileOut << E/keV << "\t"<< distEn[ien] << G4endl ;
392  normOut << E/keV << "\t"<< dnorm << G4endl ;
393  }
394  G4cout << G4endl;
395  G4int ii ;
396  G4double E1,E2 ;
397  E1=-1.e6 ;
398  E2=+1.e6 ;
399  E = Enlow -dEn ;
400  ii = -1;
401 
402  for(G4int i1=0; i1<nbinEn; i1++)
403  {
404  E += dEn ;
405  if(ii<0)
406  {
407  if(distEn[i1] >= 0.5*fmax)
408  {
409  E1=E ;
410  ii=i1 ;
411  }
412  }
413  }
414  E = Enlow -dEn ;
415 
416  for(G4int i2=0; i2<nbinEn; i2++)
417  {
418  E += dEn ;
419 
420  if(distEn[i2] >= 0.5*fmax) E2=E ;
421  }
422  G4cout << " Emp = " << std::setw(15) << Emp/MeV << " width="
423  << std::setw(15) << (E2-E1)/MeV << " MeV " << G4endl;
424  G4cout << G4endl ;
425  }
426  }
427  if(nbinTt>0)
428  {
429  G4double E , dnorm, norm ,sig;
430  G4cout << " transmitted energy distribution " << G4endl ;
431  G4cout << "#entries=" << entryTt << " #underflows=" << underTt <<
432  " #overflows=" << overTt << G4endl ;
433  if( entryTt>0.)
434  {
435  Ttmean /= entryTt;
437  if(sig<=0.)
438  sig=0.;
439  else
440  sig=std::sqrt(sig/entryTt) ;
441  G4cout << " mean energy of transmitted particles=" << Ttmean/keV <<
442  " +- " << sig/keV << " keV." << G4endl;
443  E = Ttlow - dTt ;
444  norm = TotNbofEvents*dTt ;
445  G4cout << " bin nb Elow entries normalized " << G4endl ;
446  for(G4int itt=0; itt<nbinTt; itt++)
447  {
448  E += dTt ;
449  dnorm = distTt[itt]/norm;
450  G4cout << std::setw(5) << itt << std::setw(10) << E <<
451  std::setw(12) << distTt[itt] <<
452  std::setw(12) << dnorm << G4endl ;
453  }
454  G4cout << G4endl;
455  }
456  }
457  if(nbinTb>0)
458  {
459  G4double E , dnorm, norm ,sig;
460  G4cout << " backscattered energy distribution " << G4endl ;
461  G4cout << "#entries=" << entryTb << " #underflows=" << underTb <<
462  " #overflows=" << overTb << G4endl ;
463  if( entryTb>0.)
464  {
465  Tbmean /= entryTb;
467  if(sig<=0.)
468  sig=0.;
469  else
470  sig=std::sqrt(sig/entryTb) ;
471  G4cout << " mean energy of backscattered particles=" << Tbmean/keV <<
472  " +- " << sig/keV << " keV." << G4endl;
473  E = Tblow - dTb ;
474  norm = TotNbofEvents*dTb ;
475  G4cout << " bin nb Elow entries normalized " << G4endl ;
476  for(G4int itt=0; itt<nbinTb; itt++)
477  {
478  E += dTb ;
479  dnorm = distTb[itt]/norm;
480  G4cout << std::setw(5) << itt << std::setw(10) << E <<
481  std::setw(12) << distTb[itt] <<
482  std::setw(12) << dnorm << G4endl ;
483  }
484  G4cout << G4endl;
485  }
486  }
487  if(nbinTsec>0)
488  {G4double E , dnorm, norm ;
489  G4cout << " energy distribution of charged secondaries " << G4endl ;
490  G4cout << "#entries=" << entryTsec << " #underflows=" << underTsec <<
491  " #overflows=" << overTsec << G4endl ;
492  if( entryTsec>0.)
493  {
494  E = Tseclow - dTsec ;
495  norm = TotNbofEvents*dTsec ;
496  G4cout << " bin nb Elow entries normalized " << G4endl ;
497  for(G4int itt=0; itt<nbinTsec; itt++)
498  {
499  E += dTsec ;
500  dnorm = distTsec[itt]/norm;
501  G4cout << std::setw(5) << itt << std::setw(10) << E <<
502  std::setw(12) << distTsec[itt] <<
503  std::setw(12) << dnorm << G4endl ;
504  }
505  G4cout << G4endl;
506  }
507  }
508 
509  if(nbinR >0)
510  {G4double R , dnorm, norm,sig ;
511  G4cout << " R distribution " << G4endl ;
512  G4cout << "#entries=" << entryR << " #underflows=" << underR <<
513  " #overflows=" << overR << G4endl ;
514  if( entryR >0.)
515  {
516  Rmean /= entryR;
517  sig = R2mean/entryR - Rmean*Rmean;
518  if(sig<=0.) sig=0. ;
519  else sig = std::sqrt(sig/entryR) ;
520  G4cout << " mean lateral displacement at exit=" << Rmean/mm << " +- "
521  << sig/mm << " mm." << G4endl ;
522  R = Rlow - dR ;
523  norm = TotNbofEvents*dR ;
524  G4cout << " bin nb Rlow entries normalized " << G4endl ;
525  for(G4int ier=0; ier<nbinR ; ier++)
526  {
527  R+= dR ;
528  dnorm = distR[ier]/norm;
529  G4cout << std::setw(5) << ier << std::setw(10) << R <<
530  std::setw(12) << distR[ier] <<
531  std::setw(12) << dnorm << G4endl ;
532  }
533  G4cout << G4endl;
534  }
535  }
536 
537  if(nbinTh>0)
538  {G4double Th,Thdeg, dnorm, norm,fac0,fnorm,pere,Thpere,Thmean,sum;
539  G4cout << " angle distribution " << G4endl ;
540  G4cout << "#entries=" << entryTh << " #underflows=" << underTh <<
541  " #overflows=" << overTh << G4endl ;
542  if( entryTh>0.)
543  {
544  Th= Thlow - dTh ;
545  norm = TotNbofEvents ;
546  if(distTh[0] == 0.)
547  fac0 = 1. ;
548  else
549  fac0 = 1./distTh[0] ;
550  pere = 1./std::exp(1.) ;
551 
552  G4cout << " bin nb Thlowdeg Thlowrad " <<
553  " entries normalized " << G4endl ;
554  Thpere = 0. ;
555  sum = 0. ;
556  Thmean = 0. ;
557  for(G4int ien=0; ien<nbinTh; ien++)
558  {
559  Th+= dTh ;
560  Thdeg = Th*180./pi ;
561  sum += distTh[ien] ;
562  Thmean += distTh[ien]*(Th+0.5*dTh) ;
563  dnorm = distTh[ien]/norm;
564  fnorm = fac0*distTh[ien] ;
565  if( fnorm > pere)
566  Thpere = Th ;
567  G4cout << std::setw(5) << ien << std::setw(10) << Thdeg << " " <<
568  std::setw(10) << Th << " " <<
569  std::setw(12) << distTh[ien] << " " <<
570  std::setw(12) << dnorm << " " << std::setw(12) << fnorm <<G4endl ;
571  }
572  Thmean /= sum ;
573  G4cout << G4endl;
574  G4cout << " mean = " << Thmean << " rad or " << 180.*Thmean/pi <<
575  " deg." << G4endl;
576  G4cout << " theta(1/e)=" << Thpere << " - " << Thpere+dTh << " rad "
577  << " or " << 180.*Thpere/pi << " - " << 180.*(Thpere+dTh)/pi
578  << " deg." << G4endl;
579  G4cout << G4endl;
580  }
581  }
582 
583  if(nbinThback>0)
584  {G4double Thb,Thdegb, dnormb, normb,fac0b,fnormb,pereb,Thpereb,Thmeanb,sumb;
585  G4cout << " backscattering angle distribution " << G4endl ;
586  G4cout << "#entries=" << entryThback << " #underflows=" << underThback <<
587  " #overflows=" << overThback << G4endl ;
588  if( entryThback>0.)
589  {
590  Thb= Thlowback - dThback ;
591  normb = TotNbofEvents ;
592  if(distThback[0] == 0.)
593  fac0b = 1. ;
594  else
595  fac0b = 1./distThback[0] ;
596  pereb = 1./std::exp(1.) ;
597 
598  G4cout << " bin nb Thlowdeg Thlowrad " <<
599  " entries normalized " << G4endl ;
600  Thpereb = 0. ;
601  sumb = 0. ;
602  Thmeanb = 0. ;
603  for(G4int ien=0; ien<nbinThback; ien++)
604  {
605  Thb+= dThback ;
606  Thdegb = Thb*180./pi ;
607  sumb += distThback[ien] ;
608  Thmeanb += distThback[ien]*(Thb+0.5*dThback) ;
609  dnormb = distThback[ien]/normb;
610  fnormb = fac0b*distThback[ien] ;
611  if( fnormb > pereb)
612  Thpereb = Thb ;
613  G4cout << std::setw(5) << ien << std::setw(10) << Thdegb << " " <<
614  std::setw(10) << Thb << " " <<
615  std::setw(12) << distThback[ien] << " " <<
616  std::setw(12) << dnormb << " " << std::setw(12) << fnormb <<G4endl ;
617  }
618  Thmeanb /= sumb ;
619  G4cout << G4endl;
620  G4cout << " mean = " << Thmeanb << " rad or " << 180.*Thmeanb/pi <<
621  " deg." << G4endl;
622  G4cout << " theta(1/e)=" << Thpereb << " - " << Thpereb+dThback << " rad "
623  << " or " << 180.*Thpereb/pi << " - " << 180.*(Thpereb+dThback)/pi
624  << " deg." << G4endl;
625  G4cout << G4endl;
626  }
627  }
628 
629  if(nbinGamma>0)
630  {G4double E , fact,dnorm, norm ;
631  G4cout << " gamma energy distribution " << G4endl ;
632  G4cout << "#entries=" << entryGamma << " #underflows=" << underGamma <<
633  " #overflows=" << overGamma << G4endl ;
634  if( entryGamma>0.)
635  {
636  fact=std::exp(dEGamma) ;
637  E = ElowGamma/fact ;
638  norm = TotNbofEvents*dEGamma;
639  G4cout << " bin nb Elow entries normalized " << G4endl ;
640  for(G4int itt=0; itt<nbinGamma; itt++)
641  {
642  E *= fact ;
643  dnorm = distGamma[itt]/norm;
644  G4cout << std::setw(5) << itt << std::setw(13) << E <<
645  std::setw(12) << distGamma[itt] <<
646  std::setw(15) << dnorm << G4endl ;
647  }
648  G4cout << G4endl;
649  }
650  }
651 
652  if(nbinvertexz >0)
653  {G4double z , dnorm, norm ;
654  G4cout << " vertex Z distribution " << G4endl ;
655  G4cout << "#entries=" << entryvertexz << " #underflows=" << undervertexz <<
656  " #overflows=" << oververtexz << G4endl ;
657  if( entryvertexz >0.)
658  {
659  z =zlow - dz ;
660  norm = TotNbofEvents*dz ;
661  G4cout << " bin nb zlow entries normalized " << G4endl ;
662  for(G4int iez=0; iez<nbinvertexz ; iez++)
663  {
664  z+= dz ;
665  if(std::fabs(z)<1.e-12) z=0.;
666  dnorm = distvertexz[iez]/norm;
667  G4cout << std::setw(5) << iez << std::setw(10) << z <<
668  std::setw(12) << distvertexz[iez] <<
669  std::setw(12) << dnorm << G4endl ;
670  }
671  G4cout << G4endl;
672  }
673  }
674 
675  G4cout.precision(prec);
676 
678  {
679  G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/update");
680  }
681 
682  // save Rndm status
683 
684  if (saveRndm == 1)
685  {
687  CLHEP::HepRandom::saveEngineStatus("endOfRun.rndm");
688  }
689 }
G4double nStepSumCharged
G4double distvertexz[200]
G4double Sum2Charged
static const double MeV
Definition: G4SIunits.hh:211
G4double entryThback
G4double underTh
G4double TotNbofEvents
G4double Transmitted
G4double dEGamma
const G4double fact
G4double entryGamma
G4double distTsec[200]
G4double distTh[200]
static G4VVisManager * GetConcreteInstance()
G4double underEn
G4double entryTb
G4double overTsec
G4double distTt[200]
G4double Thlowback
Float_t norm
G4double dThback
G4double distEn[200]
G4double distStep[200]
G4double overThback
G4double SumCharged
int G4int
Definition: G4Types.hh:78
G4double entryEn
static G4UImanager * GetUIpointer()
Definition: G4UImanager.cc:58
G4double overGamma
static const double prec
Definition: RanecuEngine.cc:58
G4double distGamma[200]
G4GLOB_DLL std::ostream G4cout
G4double Selectron
G4double nStepSum2Neutral
G4double distThback[200]
G4double underStep
G4double undervertexz
G4double tlSumAbs
G4double Sum2Neutral
G4double nStepSum2Charged
G4double distTb[200]
G4double ElowGamma
G4double distR[200]
G4double oververtexz
G4double tlsquareSumAbs
G4double underTsec
G4double entryvertexz
static void showEngineStatus()
Definition: Random.cc:204
G4double Steplow
G4double nStepSumNeutral
G4double entryTsec
G4double underThback
G4double entryTh
G4double entryTt
G4double Reflected
G4double SumNeutral
static const double pi
Definition: G4SIunits.hh:74
G4double Tseclow
G4double Tt2mean
G4double Spositron
G4double EnergySumAbs
G4double underTt
static void saveEngineStatus(const char filename[]="Config.conf")
Definition: Random.cc:176
G4double entryStep
#define G4endl
Definition: G4ios.hh:61
static const double keV
Definition: G4SIunits.hh:213
G4double underGamma
double G4double
Definition: G4Types.hh:76
static const double mm
Definition: G4SIunits.hh:114
G4double Tb2mean
G4double overStep
G4double EnergySquareSumAbs
G4int ApplyCommand(const char *aCommand)
Definition: G4UImanager.cc:446
G4double underTb
Here is the call graph for this function:

◆ FillEn()

void Em10RunAction::FillEn ( G4double  En)

Definition at line 746 of file Em10RunAction.cc.

747 {
748  G4double bin ;
749  G4int ibin;
750 
751  entryEn += 1. ;
752  if(Enlow < Enhigh) {
753 
754  if(En < Enlow) underEn += 1. ;
755  else if( En >= Enhigh) overEn += 1. ;
756  else
757  {
758  bin = (En-Enlow)/dEn;
759  ibin= (G4int)bin;
760  if(ibin < 0) { ibin = 0; }
761  if(ibin > 199) { ibin = 199; }
762  distEn[ibin] += 1.;
763  }
764  }
765 }
G4double underEn
float bin[41]
Definition: plottest35.C:14
G4double distEn[200]
int G4int
Definition: G4Types.hh:78
G4double entryEn
double G4double
Definition: G4Types.hh:76
Here is the caller graph for this function:

◆ FillGammaSpectrum()

void Em10RunAction::FillGammaSpectrum ( G4double  E)

Definition at line 787 of file Em10RunAction.cc.

788 {
789 }
Here is the caller graph for this function:

◆ FillNbOfSteps()

void Em10RunAction::FillNbOfSteps ( G4double  nstep)

Definition at line 740 of file Em10RunAction.cc.

741 {
742 }
Here is the caller graph for this function:

◆ FillR()

void Em10RunAction::FillR ( G4double  R)

Definition at line 805 of file Em10RunAction.cc.

806 {
807 }
Here is the caller graph for this function:

◆ FillTb()

void Em10RunAction::FillTb ( G4double  Tt)

Definition at line 775 of file Em10RunAction.cc.

776 {
777 }
Here is the caller graph for this function:

◆ FillTh()

void Em10RunAction::FillTh ( G4double  Th)

Definition at line 793 of file Em10RunAction.cc.

794 {
795 }
Here is the caller graph for this function:

◆ FillThBack()

void Em10RunAction::FillThBack ( G4double  Th)

Definition at line 799 of file Em10RunAction.cc.

800 {
801 }
Here is the caller graph for this function:

◆ FillTsec()

void Em10RunAction::FillTsec ( G4double  T)

Definition at line 781 of file Em10RunAction.cc.

782 {
783 }
Here is the caller graph for this function:

◆ FillTt()

void Em10RunAction::FillTt ( G4double  Tt)

Definition at line 769 of file Em10RunAction.cc.

770 {
771 }
Here is the caller graph for this function:

◆ Fillvertexz()

void Em10RunAction::Fillvertexz ( G4double  z)

Definition at line 811 of file Em10RunAction.cc.

812 {
813 }
Here is the caller graph for this function:

◆ GetRndmFreq()

G4int Em10RunAction::GetRndmFreq ( )
inline

Definition at line 125 of file Em10RunAction.hh.

125 {return saveRndm;}
Here is the call graph for this function:
Here is the caller graph for this function:

◆ SetEhighGamma()

void Em10RunAction::SetEhighGamma ( G4double  Ehigh)

Definition at line 890 of file Em10RunAction.cc.

891 {
892  if(Ehigh > 0.0) {
893  EhighGamma = Ehigh;
894  G4cout << " Ehigh in the gamma spectrum plot = " << EhighGamma << G4endl ;
895  }
896 }
G4GLOB_DLL std::ostream G4cout
G4double EhighGamma
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetElowGamma()

void Em10RunAction::SetElowGamma ( G4double  Elow)

Definition at line 882 of file Em10RunAction.cc.

883 {
884  if(Elow >= 0.0) {
885  ElowGamma = Elow;
886  G4cout << " Elow in the gamma spectrum plot = " << ElowGamma << G4endl ;
887  }
888 }
G4GLOB_DLL std::ostream G4cout
G4double ElowGamma
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetEnhigh()

void Em10RunAction::SetEnhigh ( G4double  Enhigh)

Definition at line 864 of file Em10RunAction.cc.

865 {
866  if(Ehigh > 0.0) {
867  Enhigh = Ehigh ;
868  G4cout << " Ehigh in the Edep plot = " << Enhigh << G4endl ;
869  }
870 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetEnlow()

void Em10RunAction::SetEnlow ( G4double  Elow)

Definition at line 856 of file Em10RunAction.cc.

857 {
858  if(Elow >= 0.0) {
859  Enlow = Elow ;
860  G4cout << " Elow in the Edep plot = " << Enlow << G4endl ;
861  }
862 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SethistName()

void Em10RunAction::SethistName ( G4String  name)

Definition at line 817 of file Em10RunAction.cc.

818 {
819  histName = name ;
820  G4cout << " hist file = " << histName << G4endl;
821 }
G4String name
Definition: TRTMaterials.hh:40
G4GLOB_DLL std::ostream G4cout
G4String histName
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinEn()

void Em10RunAction::SetnbinEn ( G4int  nbin)

Definition at line 848 of file Em10RunAction.cc.

849 {
850  if(nbin > 0 && nbin <= 200) {
851  nbinEn = nbin;
852  G4cout << " Nb of bins in Edep plot = " << nbinEn << G4endl ;
853  }
854 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinGamma()

void Em10RunAction::SetnbinGamma ( G4int  nbin)

Definition at line 874 of file Em10RunAction.cc.

875 {
876  if(nbin > 0 && nbin <= 200) {
877  nbinGamma = nbin;
878  G4cout << " Nb of bins in gamma spectrum plot = " << nbinGamma << G4endl ;
879  }
880 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinR()

void Em10RunAction::SetnbinR ( G4int  nbin)

Definition at line 963 of file Em10RunAction.cc.

964 {
965  if(nbin > 0 && nbin <= 200) {
966  nbinR = nbin;
967  G4cout << " Nb of bins in R plot = " << nbinR << G4endl ;
968  }
969 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinStep()

void Em10RunAction::SetnbinStep ( G4int  nbin)

Definition at line 823 of file Em10RunAction.cc.

824 {
825  if(nbin> 0 && nbin<= 200) {
826  nbinStep = nbin;
827  G4cout << " Nb of bins in #step plot = " << nbinStep << G4endl ;
828  }
829 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinTb()

void Em10RunAction::SetnbinTb ( G4int  nbin)

Definition at line 922 of file Em10RunAction.cc.

923 {
924  if(nbin > 0 && nbin <= 200) {
925  nbinTb = nbin;
926  G4cout << " Nb of bins in Ebackscattered plot = " << nbinTb << G4endl ;
927  }
928 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinTh()

void Em10RunAction::SetnbinTh ( G4int  nbin)

Definition at line 1003 of file Em10RunAction.cc.

1004 {
1005  if(nbin > 0 && nbin <= 200) {
1006  nbinTh = nbin;
1007  G4cout << " Nb of bins in Theta plot = " << nbinTh << G4endl ;
1008  }
1009 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinThBack()

void Em10RunAction::SetnbinThBack ( G4int  nbin)

Definition at line 1023 of file Em10RunAction.cc.

1024 {
1025  if(nbin > 0 && nbin <= 200) {
1026  nbinThback = nbin;
1027  G4cout << " Nb of bins in Theta plot = " << nbinThback << G4endl ;
1028  }
1029 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinTsec()

void Em10RunAction::SetnbinTsec ( G4int  nbin)

Definition at line 943 of file Em10RunAction.cc.

944 {
945  if(nbin > 0 && nbin <= 200) {
946  nbinTsec = nbin;
947  G4cout << " Nb of bins in Tsecondary plot = " << nbinTsec << G4endl ;
948  }
949 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetnbinTt()

void Em10RunAction::SetnbinTt ( G4int  nbin)

Definition at line 898 of file Em10RunAction.cc.

899 {
900  if(nbin > 0 && nbin <= 200) {
901  nbinTt = nbin;
902  G4cout << " Nb of bins in Etransmisssion plot = " << nbinTt << G4endl ;
903  }
904 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ Setnbinzvertex()

void Em10RunAction::Setnbinzvertex ( G4int  nbin)

Definition at line 983 of file Em10RunAction.cc.

984 {
985  if(nbin > 0 && nbin <= 200) {
986  nbinvertexz = nbin;
987  G4cout << " Nb of bins in Z plot = " << nbinvertexz << G4endl ;
988  }
989 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetRhigh()

void Em10RunAction::SetRhigh ( G4double  Rhigh)

Definition at line 977 of file Em10RunAction.cc.

978 {
979  Rhigh = rhigh ;
980  G4cout << " Rhigh in the R plot = " << Rhigh << G4endl ;
981 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetRlow()

void Em10RunAction::SetRlow ( G4double  Rlow)

Definition at line 971 of file Em10RunAction.cc.

972 {
973  Rlow = rlow ;
974  G4cout << " Rlow in the R plot = " << Rlow << G4endl ;
975 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetRndmFreq()

void Em10RunAction::SetRndmFreq ( G4int  val)
inline

Definition at line 124 of file Em10RunAction.hh.

124 {saveRndm = val;}
Here is the caller graph for this function:

◆ SetStephigh()

void Em10RunAction::SetStephigh ( G4double  Shigh)

Definition at line 838 of file Em10RunAction.cc.

839 {
840  if(high > 0.0) {
841  Stephigh = high;
842  G4cout << " high in the #step plot = " << Stephigh << G4endl ;
843  }
844 }
G4double Stephigh
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetSteplow()

void Em10RunAction::SetSteplow ( G4double  Slow)

Definition at line 831 of file Em10RunAction.cc.

832 {
833  if(low >= 0.0){
834  Steplow = low ;
835  G4cout << " low in the #step plot = " << Steplow << G4endl ;
836  }
837 }
G4GLOB_DLL std::ostream G4cout
G4double Steplow
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetTbhigh()

void Em10RunAction::SetTbhigh ( G4double  Tbhigh)

Definition at line 937 of file Em10RunAction.cc.

938 {
939  Tbhigh = Ehigh ;
940  G4cout << " Ehigh in the Ebackscattered plot = " << Tbhigh << G4endl ;
941 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetTblow()

void Em10RunAction::SetTblow ( G4double  Tblow)

Definition at line 930 of file Em10RunAction.cc.

931 {
932 
933  Tblow = Elow ;
934  G4cout << " Elow in the Ebackscattered plot = " << Tblow << G4endl ;
935 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetThhigh()

void Em10RunAction::SetThhigh ( G4double  Thhigh)

Definition at line 1017 of file Em10RunAction.cc.

1018 {
1019  Thhigh = Thigh ;
1020  G4cout << " Thigh in the Theta plot = " << Thhigh << G4endl ;
1021 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetThhighBack()

void Em10RunAction::SetThhighBack ( G4double  Thhigh)

Definition at line 1037 of file Em10RunAction.cc.

1038 {
1039  Thhighback = Thigh ;
1040  G4cout << " Thigh in the Theta plot = " << Thhighback << G4endl ;
1041 }
G4GLOB_DLL std::ostream G4cout
G4double Thhighback
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetThlow()

void Em10RunAction::SetThlow ( G4double  Thlow)

Definition at line 1011 of file Em10RunAction.cc.

1012 {
1013  Thlow = Tlow ;
1014  G4cout << " Tlow in the Theta plot = " << Thlow << G4endl ;
1015 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetThlowBack()

void Em10RunAction::SetThlowBack ( G4double  Thlow)

Definition at line 1031 of file Em10RunAction.cc.

1032 {
1033  Thlowback = Tlow ;
1034  G4cout << " Tlow in the Theta plot = " << Thlowback << G4endl ;
1035 }
G4double Thlowback
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetTsechigh()

void Em10RunAction::SetTsechigh ( G4double  Thigh)

Definition at line 957 of file Em10RunAction.cc.

958 {
959  Tsechigh = Ehigh ;
960  G4cout << " Ehigh in the Tsecondary plot = " << Tsechigh << G4endl ;
961 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
G4double Tsechigh
Here is the caller graph for this function:

◆ SetTseclow()

void Em10RunAction::SetTseclow ( G4double  Tlow)

Definition at line 951 of file Em10RunAction.cc.

952 {
953  Tseclow = Elow ;
954  G4cout << " Elow in the Tsecondary plot = " << Tseclow << G4endl ;
955 }
G4GLOB_DLL std::ostream G4cout
G4double Tseclow
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetTthigh()

void Em10RunAction::SetTthigh ( G4double  Tthigh)

Definition at line 914 of file Em10RunAction.cc.

915 {
916  if(Ehigh > 0.0) {
917  Tthigh = Ehigh;
918  G4cout << " Ehigh in the Etransmission plot = " << Tthigh << G4endl ;
919  }
920 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ SetTtlow()

void Em10RunAction::SetTtlow ( G4double  Ttlow)

Definition at line 906 of file Em10RunAction.cc.

907 {
908  if(Elow >= 0.0) {
909  Ttlow = Elow;
910  G4cout << " Elow in the Etransmission plot = " << Ttlow << G4endl ;
911  }
912 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ Setzhigh()

void Em10RunAction::Setzhigh ( G4double  z)

Definition at line 997 of file Em10RunAction.cc.

998 {
999  zhigh = z ;
1000  G4cout << " zhigh in the Z plot = " << zhigh << G4endl ;
1001 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

◆ Setzlow()

void Em10RunAction::Setzlow ( G4double  z)

Definition at line 991 of file Em10RunAction.cc.

992 {
993  zlow = z ;
994  G4cout << " zlow in the Z plot = " << zlow << G4endl ;
995 }
G4GLOB_DLL std::ostream G4cout
#define G4endl
Definition: G4ios.hh:61
Here is the caller graph for this function:

Member Data Documentation

◆ dEGamma

G4double Em10RunAction::dEGamma
private

Definition at line 173 of file Em10RunAction.hh.

◆ dEn

G4double Em10RunAction::dEn
private

Definition at line 149 of file Em10RunAction.hh.

◆ distEn

G4double Em10RunAction::distEn[200]
private

Definition at line 148 of file Em10RunAction.hh.

◆ distGamma

G4double Em10RunAction::distGamma[200]
private

Definition at line 172 of file Em10RunAction.hh.

◆ distR

G4double Em10RunAction::distR[200]
private

Definition at line 168 of file Em10RunAction.hh.

◆ distStep

G4double Em10RunAction::distStep[200]
private

Definition at line 145 of file Em10RunAction.hh.

◆ distTb

G4double Em10RunAction::distTb[200]
private

Definition at line 155 of file Em10RunAction.hh.

◆ distTh

G4double Em10RunAction::distTh[200]
private

Definition at line 162 of file Em10RunAction.hh.

◆ distThback

G4double Em10RunAction::distThback[200]
private

Definition at line 165 of file Em10RunAction.hh.

◆ distTsec

G4double Em10RunAction::distTsec[200]
private

Definition at line 159 of file Em10RunAction.hh.

◆ distTt

G4double Em10RunAction::distTt[200]
private

Definition at line 151 of file Em10RunAction.hh.

◆ distvertexz

G4double Em10RunAction::distvertexz[200]
private

Definition at line 175 of file Em10RunAction.hh.

◆ dR

G4double Em10RunAction::dR
private

Definition at line 169 of file Em10RunAction.hh.

◆ dStep

G4double Em10RunAction::dStep
private

Definition at line 146 of file Em10RunAction.hh.

◆ dTb

G4double Em10RunAction::dTb
private

Definition at line 156 of file Em10RunAction.hh.

◆ dTh

G4double Em10RunAction::dTh
private

Definition at line 163 of file Em10RunAction.hh.

◆ dThback

G4double Em10RunAction::dThback
private

Definition at line 166 of file Em10RunAction.hh.

◆ dTsec

G4double Em10RunAction::dTsec
private

Definition at line 160 of file Em10RunAction.hh.

◆ dTt

G4double Em10RunAction::dTt
private

Definition at line 152 of file Em10RunAction.hh.

◆ dz

G4double Em10RunAction::dz
private

Definition at line 176 of file Em10RunAction.hh.

◆ EhighGamma

G4double Em10RunAction::EhighGamma
private

Definition at line 173 of file Em10RunAction.hh.

◆ ElowGamma

G4double Em10RunAction::ElowGamma
private

Definition at line 173 of file Em10RunAction.hh.

◆ EnergySquareSumAbs

G4double Em10RunAction::EnergySquareSumAbs
private

Definition at line 135 of file Em10RunAction.hh.

◆ EnergySumAbs

G4double Em10RunAction::EnergySumAbs
private

Definition at line 135 of file Em10RunAction.hh.

◆ Enhigh

G4double Em10RunAction::Enhigh
private

Definition at line 149 of file Em10RunAction.hh.

◆ Enlow

G4double Em10RunAction::Enlow
private

Definition at line 149 of file Em10RunAction.hh.

◆ entryEn

G4double Em10RunAction::entryEn
private

Definition at line 148 of file Em10RunAction.hh.

◆ entryGamma

G4double Em10RunAction::entryGamma
private

Definition at line 172 of file Em10RunAction.hh.

◆ entryR

G4double Em10RunAction::entryR
private

Definition at line 168 of file Em10RunAction.hh.

◆ entryStep

G4double Em10RunAction::entryStep
private

Definition at line 145 of file Em10RunAction.hh.

◆ entryTb

G4double Em10RunAction::entryTb
private

Definition at line 155 of file Em10RunAction.hh.

◆ entryTh

G4double Em10RunAction::entryTh
private

Definition at line 162 of file Em10RunAction.hh.

◆ entryThback

G4double Em10RunAction::entryThback
private

Definition at line 165 of file Em10RunAction.hh.

◆ entryTsec

G4double Em10RunAction::entryTsec
private

Definition at line 159 of file Em10RunAction.hh.

◆ entryTt

G4double Em10RunAction::entryTt
private

Definition at line 151 of file Em10RunAction.hh.

◆ entryvertexz

G4double Em10RunAction::entryvertexz
private

Definition at line 175 of file Em10RunAction.hh.

◆ histName

G4String Em10RunAction::histName
private

Definition at line 133 of file Em10RunAction.hh.

◆ nbinEn

G4int Em10RunAction::nbinEn
private

Definition at line 150 of file Em10RunAction.hh.

◆ nbinGamma

G4int Em10RunAction::nbinGamma
private

Definition at line 174 of file Em10RunAction.hh.

◆ nbinR

G4int Em10RunAction::nbinR
private

Definition at line 170 of file Em10RunAction.hh.

◆ nbinStep

G4int Em10RunAction::nbinStep
private

Definition at line 147 of file Em10RunAction.hh.

◆ nbinTb

G4int Em10RunAction::nbinTb
private

Definition at line 157 of file Em10RunAction.hh.

◆ nbinTh

G4int Em10RunAction::nbinTh
private

Definition at line 164 of file Em10RunAction.hh.

◆ nbinThback

G4int Em10RunAction::nbinThback
private

Definition at line 167 of file Em10RunAction.hh.

◆ nbinTsec

G4int Em10RunAction::nbinTsec
private

Definition at line 161 of file Em10RunAction.hh.

◆ nbinTt

G4int Em10RunAction::nbinTt
private

Definition at line 153 of file Em10RunAction.hh.

◆ nbinvertexz

G4int Em10RunAction::nbinvertexz
private

Definition at line 177 of file Em10RunAction.hh.

◆ nStepSum2Charged

G4double Em10RunAction::nStepSum2Charged
private

Definition at line 137 of file Em10RunAction.hh.

◆ nStepSum2Neutral

G4double Em10RunAction::nStepSum2Neutral
private

Definition at line 138 of file Em10RunAction.hh.

◆ nStepSumCharged

G4double Em10RunAction::nStepSumCharged
private

Definition at line 137 of file Em10RunAction.hh.

◆ nStepSumNeutral

G4double Em10RunAction::nStepSumNeutral
private

Definition at line 138 of file Em10RunAction.hh.

◆ overEn

G4double Em10RunAction::overEn
private

Definition at line 148 of file Em10RunAction.hh.

◆ overGamma

G4double Em10RunAction::overGamma
private

Definition at line 172 of file Em10RunAction.hh.

◆ overR

G4double Em10RunAction::overR
private

Definition at line 168 of file Em10RunAction.hh.

◆ overStep

G4double Em10RunAction::overStep
private

Definition at line 145 of file Em10RunAction.hh.

◆ overTb

G4double Em10RunAction::overTb
private

Definition at line 155 of file Em10RunAction.hh.

◆ overTh

G4double Em10RunAction::overTh
private

Definition at line 162 of file Em10RunAction.hh.

◆ overThback

G4double Em10RunAction::overThback
private

Definition at line 165 of file Em10RunAction.hh.

◆ overTsec

G4double Em10RunAction::overTsec
private

Definition at line 159 of file Em10RunAction.hh.

◆ overTt

G4double Em10RunAction::overTt
private

Definition at line 151 of file Em10RunAction.hh.

◆ oververtexz

G4double Em10RunAction::oververtexz
private

Definition at line 175 of file Em10RunAction.hh.

◆ R2mean

G4double Em10RunAction::R2mean
private

Definition at line 171 of file Em10RunAction.hh.

◆ Reflected

G4double Em10RunAction::Reflected
private

Definition at line 143 of file Em10RunAction.hh.

◆ Rhigh

G4double Em10RunAction::Rhigh
private

Definition at line 169 of file Em10RunAction.hh.

◆ Rlow

G4double Em10RunAction::Rlow
private

Definition at line 169 of file Em10RunAction.hh.

◆ Rmean

G4double Em10RunAction::Rmean
private

Definition at line 171 of file Em10RunAction.hh.

◆ runMessenger

Em10RunMessenger* Em10RunAction::runMessenger
private

Definition at line 179 of file Em10RunAction.hh.

◆ saveRndm

G4int Em10RunAction::saveRndm
private

Definition at line 180 of file Em10RunAction.hh.

◆ Selectron

G4double Em10RunAction::Selectron
private

Definition at line 141 of file Em10RunAction.hh.

◆ Spositron

G4double Em10RunAction::Spositron
private

Definition at line 141 of file Em10RunAction.hh.

◆ Stephigh

G4double Em10RunAction::Stephigh
private

Definition at line 146 of file Em10RunAction.hh.

◆ Steplow

G4double Em10RunAction::Steplow
private

Definition at line 146 of file Em10RunAction.hh.

◆ Sum2Charged

G4double Em10RunAction::Sum2Charged
private

Definition at line 140 of file Em10RunAction.hh.

◆ Sum2Neutral

G4double Em10RunAction::Sum2Neutral
private

Definition at line 140 of file Em10RunAction.hh.

◆ SumCharged

G4double Em10RunAction::SumCharged
private

Definition at line 140 of file Em10RunAction.hh.

◆ SumNeutral

G4double Em10RunAction::SumNeutral
private

Definition at line 140 of file Em10RunAction.hh.

◆ Tb2mean

G4double Em10RunAction::Tb2mean
private

Definition at line 158 of file Em10RunAction.hh.

◆ Tbhigh

G4double Em10RunAction::Tbhigh
private

Definition at line 156 of file Em10RunAction.hh.

◆ Tblow

G4double Em10RunAction::Tblow
private

Definition at line 156 of file Em10RunAction.hh.

◆ Tbmean

G4double Em10RunAction::Tbmean
private

Definition at line 158 of file Em10RunAction.hh.

◆ Thhigh

G4double Em10RunAction::Thhigh
private

Definition at line 163 of file Em10RunAction.hh.

◆ Thhighback

G4double Em10RunAction::Thhighback
private

Definition at line 166 of file Em10RunAction.hh.

◆ Thlow

G4double Em10RunAction::Thlow
private

Definition at line 163 of file Em10RunAction.hh.

◆ Thlowback

G4double Em10RunAction::Thlowback
private

Definition at line 166 of file Em10RunAction.hh.

◆ tlsquareSumAbs

G4double Em10RunAction::tlsquareSumAbs
private

Definition at line 136 of file Em10RunAction.hh.

◆ tlSumAbs

G4double Em10RunAction::tlSumAbs
private

Definition at line 136 of file Em10RunAction.hh.

◆ TotNbofEvents

G4double Em10RunAction::TotNbofEvents
private

Definition at line 139 of file Em10RunAction.hh.

◆ Transmitted

G4double Em10RunAction::Transmitted
private

Definition at line 143 of file Em10RunAction.hh.

◆ Tsechigh

G4double Em10RunAction::Tsechigh
private

Definition at line 160 of file Em10RunAction.hh.

◆ Tseclow

G4double Em10RunAction::Tseclow
private

Definition at line 160 of file Em10RunAction.hh.

◆ Tt2mean

G4double Em10RunAction::Tt2mean
private

Definition at line 154 of file Em10RunAction.hh.

◆ Tthigh

G4double Em10RunAction::Tthigh
private

Definition at line 152 of file Em10RunAction.hh.

◆ Ttlow

G4double Em10RunAction::Ttlow
private

Definition at line 152 of file Em10RunAction.hh.

◆ Ttmean

G4double Em10RunAction::Ttmean
private

Definition at line 154 of file Em10RunAction.hh.

◆ underEn

G4double Em10RunAction::underEn
private

Definition at line 148 of file Em10RunAction.hh.

◆ underGamma

G4double Em10RunAction::underGamma
private

Definition at line 172 of file Em10RunAction.hh.

◆ underR

G4double Em10RunAction::underR
private

Definition at line 168 of file Em10RunAction.hh.

◆ underStep

G4double Em10RunAction::underStep
private

Definition at line 145 of file Em10RunAction.hh.

◆ underTb

G4double Em10RunAction::underTb
private

Definition at line 155 of file Em10RunAction.hh.

◆ underTh

G4double Em10RunAction::underTh
private

Definition at line 162 of file Em10RunAction.hh.

◆ underThback

G4double Em10RunAction::underThback
private

Definition at line 165 of file Em10RunAction.hh.

◆ underTsec

G4double Em10RunAction::underTsec
private

Definition at line 159 of file Em10RunAction.hh.

◆ underTt

G4double Em10RunAction::underTt
private

Definition at line 151 of file Em10RunAction.hh.

◆ undervertexz

G4double Em10RunAction::undervertexz
private

Definition at line 175 of file Em10RunAction.hh.

◆ zhigh

G4double Em10RunAction::zhigh
private

Definition at line 176 of file Em10RunAction.hh.

◆ zlow

G4double Em10RunAction::zlow
private

Definition at line 176 of file Em10RunAction.hh.


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