Em10RunAction Class Reference

#include <Em10RunAction.hh>

Inheritance 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 G4Run *	GenerateRun ()
void	SetMaster (G4bool val=true)
G4bool	IsMaster () const

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

(

)

Definition at line 50 of file Em10RunAction.cc.

  :G4UserRunAction(),
   histName("histfile"),nbinStep(0),nbinEn(0),nbinTt(0),nbinTb(0),
   nbinTsec(0),nbinTh(0),nbinThback(0),nbinR(0),nbinGamma(0),
   nbinvertexz(0)
{
  runMessenger = new Em10RunMessenger(this);
  saveRndm = 1;  

  EnergySumAbs = EnergySquareSumAbs = tlSumAbs = tlsquareSumAbs = 
    nStepSumCharged = nStepSum2Charged= nStepSumNeutral = nStepSum2Neutral=
    TotNbofEvents = SumCharged= SumNeutral=Sum2Charged=Sum2Neutral=Selectron=
    Spositron=Transmitted=Reflected  =dStep=entryStep=underStep=overStep=
    dEn = entryEn= underEn=overEn=dTt = entryTt=underTt=overTt=Ttmean=
    Tt2mean=dTb = entryTb=underTb=overTb=Tbmean=Tb2mean=dTsec =entryTsec=
    underTsec=overTsec=dTh = entryTh=underTh=overTh=dThback =entryThback=
    underThback=overThback=dR  =entryR =underR =overR =Rmean=R2mean=dEGamma= 
    entryGamma = underGamma=overGamma=dz=entryvertexz=undervertexz=
    oververtexz=0.;
}

Em10RunAction::~Em10RunAction

(

)

Definition at line 73 of file Em10RunAction.cc.

{
  delete runMessenger;
}

Member Function Documentation

void Em10RunAction::AddEdeps

(

G4double

Eabs

)

Definition at line 716 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  EnergySumAbs += Eabs;
  EnergySquareSumAbs += Eabs*Eabs;
}

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

Definition at line 1051 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  Selectron += nele;
  Spositron += npos;
}

void Em10RunAction::AddnStepsCharged

(

G4double

ns

)

Definition at line 700 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  nStepSumCharged += nstp;
  nStepSum2Charged += nstp*nstp;
}

void Em10RunAction::AddnStepsNeutral

(

G4double

ns

)

Definition at line 708 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  nStepSumNeutral += nstp;
  nStepSum2Neutral += nstp*nstp;
}

void Em10RunAction::AddTrackLength

(

G4double

tlabs

)

Definition at line 724 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  tlSumAbs += tlabs;
  tlsquareSumAbs += tlabs*tlabs ;
}

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

Definition at line 732 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  Transmitted += tr ;
  Reflected   += ref;
}

void Em10RunAction::BeginOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

Definition at line 86 of file Em10RunAction.cc.

References G4UImanager::ApplyCommand(), G4cout, G4endl, G4VVisManager::GetConcreteInstance(), G4Run::GetRunID(), G4UImanager::GetUIpointer(), iz, CLHEP::HepRandom::saveEngineStatus(), and CLHEP::HepRandom::showEngineStatus().

{  
  G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
  
  // save Rndm status
  if (saveRndm > 0)
  { 
      CLHEP::HepRandom::showEngineStatus();
      CLHEP::HepRandom::saveEngineStatus("beginOfRun.rndm");
  }  
  G4UImanager* UI = G4UImanager::GetUIpointer();
   
  G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance();

  if(pVVisManager)    UI->ApplyCommand("/vis/scene/notifyHandlers");

      
  EnergySumAbs = 0. ;
  EnergySquareSumAbs = 0.;
  tlSumAbs = 0. ;
  tlsquareSumAbs = 0. ;
  nStepSumCharged = 0. ;
  nStepSum2Charged= 0. ;
  nStepSumNeutral = 0. ;
  nStepSum2Neutral= 0. ;
  TotNbofEvents = 0. ;
  SumCharged=0.;
  SumNeutral=0.;
  Sum2Charged=0.;
  Sum2Neutral=0.;
  Selectron=0.;
  Spositron=0.;

  Transmitted=0.;
  Reflected  =0.;

  //  plot definitions 
   
  if(nbinStep>0)
  {
    dStep=(Stephigh-Steplow)/nbinStep;
    entryStep=0.;
    underStep=0.;
    overStep=0.;
    for(G4int ist=0; ist<200; ist++)
    {
      distStep[ist]=0.;
    }
  }      
  if(nbinEn>0)
  {
    dEn = (Enhigh-Enlow)/nbinEn ;
    entryEn=0.;
    underEn=0.;
    overEn=0.;

    for (G4int ien=0; ien<200; ien++)   distEn[ien]=0.;
  }
  if(nbinTt>0)
  {
    dTt = (Tthigh-Ttlow)/nbinTt ;
    entryTt=0.;
    underTt=0.;
    overTt=0.;

    for (G4int itt=0; itt<200; itt++)  distTt[itt]=0.;

    Ttmean=0.;
    Tt2mean=0.;
  }
  if(nbinTb>0)
  {
    dTb = (Tbhigh-Tblow)/nbinTb ;
    entryTb=0.;
    underTb=0.;
    overTb=0.;
    for (G4int itt=0; itt<200; itt++)
    {
      distTb[itt]=0.;
    }
    Tbmean=0.;
    Tb2mean=0.;
  }
  if(nbinTsec>0)
  {
    dTsec = (Tsechigh-Tseclow)/nbinTsec ;
    entryTsec=0.;
    underTsec=0.;
    overTsec=0.;
    for (G4int its=0; its<200; its++)
    {
      distTsec[its]=0.;
    }
  }
  if(nbinTh>0)
  {
    dTh = (Thhigh-Thlow)/nbinTh ;
    entryTh=0.;
    underTh=0.;
    overTh=0.;
    for (G4int ith=0; ith<200; ith++)
    {
      distTh[ith]=0.;
    }
  }

  if(nbinThback>0)
  {
    dThback = (Thhighback-Thlowback)/nbinThback ;
    entryThback=0.;
    underThback=0.;
    overThback=0.;
    for (G4int ithback=0; ithback<200; ithback++)
    {
      distThback[ithback]=0.;
    }
  }


  if(nbinR >0)
  {
    dR  = (Rhigh-Rlow)/nbinR  ;
    entryR =0.;
    underR =0.;
    overR =0.;
    for (G4int ir =0; ir<200; ir++)
    {
      distR[ir]=0.;
    }
    Rmean=0.;
    R2mean=0.;
  }

  if(nbinGamma>0)
  {
    dEGamma = std::log(EhighGamma/ElowGamma)/nbinGamma ;
    entryGamma = 0.;
    underGamma=0.;
    overGamma=0.;
    for (G4int ig=0; ig<200; ig++)
    {
      distGamma[ig]=0.;
    }
  } 
  if(nbinvertexz>0)
  {
    dz=(zhigh-zlow)/nbinvertexz;
    entryvertexz=0.;
    undervertexz=0.;
    oververtexz=0.;
    for(G4int iz=0; iz<200; iz++)
    {
      distvertexz[iz]=0.;
    }
  }

  bookHisto();
}

void Em10RunAction::CountEvent

(

)

Definition at line 693 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  TotNbofEvents += 1. ;
}

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

Definition at line 1043 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
  SumCharged += nch ;
  SumNeutral += nne ;
  Sum2Charged += nch*nch ;
  Sum2Neutral += nne*nne ;
}

void Em10RunAction::EndOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 247 of file Em10RunAction.cc.

References G4UImanager::ApplyCommand(), G4cout, G4endl, G4VVisManager::GetConcreteInstance(), G4UImanager::GetUIpointer(), python.hepunit::keV, python.hepunit::MeV, python.hepunit::mm, python.hepunit::pi, CLHEP::HepRandom::saveEngineStatus(), CLHEP::HepRandom::showEngineStatus(), and z.

{
  G4double sAbs,sigAbs,sigstep,sigcharged,signeutral;
  if(0.0 >= TotNbofEvents) { return; }

  tlSumAbs /= TotNbofEvents ;
  sAbs = tlsquareSumAbs/TotNbofEvents-tlSumAbs*tlSumAbs ;
  if(sAbs>0.)
    sAbs = std::sqrt(sAbs/TotNbofEvents) ;
  else
    sAbs = 0. ;
  
  EnergySumAbs /= TotNbofEvents ;
  sigAbs = EnergySquareSumAbs/TotNbofEvents-EnergySumAbs*EnergySumAbs;
  if(sigAbs>0.)
    sigAbs = std::sqrt(sigAbs/TotNbofEvents);
  else
    sigAbs = 0.;

  nStepSumCharged /= TotNbofEvents ;
  sigstep = nStepSum2Charged/TotNbofEvents-nStepSumCharged*nStepSumCharged;
  if(sigstep>0.)
    sigstep = std::sqrt(sigstep/TotNbofEvents);
  else
    sigstep = 0.;
  G4double sigch=sigstep ;
  
  nStepSumNeutral /= TotNbofEvents ;
  sigstep = nStepSum2Neutral/TotNbofEvents-nStepSumNeutral*nStepSumNeutral;
  if(sigstep>0.)
    sigstep = std::sqrt(sigstep/TotNbofEvents);
  else
    sigstep = 0.;
  G4double signe=sigstep ;
  
  SumCharged /= TotNbofEvents;
  sigcharged = Sum2Charged/TotNbofEvents-SumCharged*SumCharged; 
  if(sigcharged>0.)
    sigcharged = std::sqrt(sigcharged/TotNbofEvents);
  else
    sigcharged = 0. ;
 
  SumNeutral /= TotNbofEvents;
  signeutral = Sum2Neutral/TotNbofEvents-SumNeutral*SumNeutral; 
  if(signeutral>0.)
    signeutral = std::sqrt(signeutral/TotNbofEvents);
  else
    signeutral = 0. ;
 
  Selectron /= TotNbofEvents ;
  Spositron /= TotNbofEvents ;

  Transmitted /=TotNbofEvents ;
  Reflected   /=TotNbofEvents ;
  G4cout << " ================== run summary =====================" << G4endl;
  G4int prec = G4cout.precision(6);
  G4cout << " end of Run TotNbofEvents = " <<  
           TotNbofEvents << G4endl ;
  G4cout << "    mean charged track length   in absorber=" <<
           tlSumAbs/mm      << " +- " << sAbs/mm    <<
          "  mm  " << G4endl; 
  G4cout << G4endl;
  G4cout << "            mean energy deposit in absorber=" <<
           EnergySumAbs/MeV << " +- " << sigAbs/MeV <<
          "  MeV " << G4endl ;
  G4cout << G4endl ;
  G4cout << " mean number of steps in absorber (charged) =" <<
           nStepSumCharged         << " +- " << sigch     <<
          "      " << G4endl ;
  G4cout << " mean number of steps in absorber (neutral) =" <<
           nStepSumNeutral         << " +- " << signe     <<
          "      " << G4endl ;
  G4cout << G4endl ;
  G4cout << "   mean number of charged secondaries = " <<
           SumCharged << " +- " << sigcharged << G4endl;  
  G4cout << G4endl ;
  G4cout << "   mean number of neutral secondaries = " <<
           SumNeutral << " +- " << signeutral << G4endl;  
  G4cout << G4endl ;
  
  G4cout << "   mean number of e-s =" << Selectron << 
            "  and e+s =" << Spositron << G4endl;
  G4cout << G4endl; 
  
  G4cout << "(number) transmission coeff=" << Transmitted <<
            "  reflection coeff=" << Reflected << G4endl;
  G4cout << G4endl; 

  if(nbinStep>0)
  {G4double E , dnorm, norm ;
   G4cout << "   step number/event distribution " << G4endl ;
   G4cout << "#entries=" << entryStep << "    #underflows=" << underStep <<
             "    #overflows=" << overStep << G4endl ;
   if( entryStep>0.)
   {
     E = Steplow - dStep ;
     norm = TotNbofEvents ;
     G4cout << " bin nb   nsteplow     entries     normalized " << G4endl ;
     for(G4int iss=0; iss<nbinStep; iss++)
     {
      E += dStep ;
      dnorm = distStep[iss]/norm;
      G4cout << std::setw(5) << iss << std::setw(10) << E << 
                std::setw(12) << distStep[iss] <<
                std::setw(12) << dnorm << G4endl ;
     }
     G4cout << G4endl;
   }     
  }
  if(nbinEn > 0)
  {
    std::ofstream fileOut("distribution.out", std::ios::out ) ;
    fileOut.setf( std::ios::scientific, std::ios::floatfield );

    std::ofstream normOut("normDist.out", std::ios::out ) ;
    normOut.setf( std::ios::scientific, std::ios::floatfield );

    G4double E , dnorm, norm,fmax,Emp ;
    Emp=-999.999 ;
    G4cout << " energy deposit distribution " << G4endl ;
    G4cout << "#entries=" << entryEn << "    #underflows=" << underEn <<
             "    #overflows=" << overEn << G4endl ;
    if( entryEn>0.)
    {
      E = Enlow - dEn ;
      norm = TotNbofEvents*1.0 ;   // *dEn ;
      G4cout << " bin nb      Elow      entries     normalized " << G4endl ;
      fmax = 0. ;

      for(G4int ien=0; ien<nbinEn; ien++)
      {
        E += dEn ;

        if(distEn[ien]>fmax)
        {
          fmax = distEn[ien] ;
          Emp = E ;                // most probable roughly
        }
        dnorm = distEn[ien]/norm;

        G4cout << std::setw(5) << ien << std::setw(10) << E/keV << 
                  std::setw(12) << distEn[ien] <<
                  std::setw(12) << dnorm << G4endl ;

        fileOut << E/keV << "\t"<< distEn[ien] << G4endl ;
        normOut << E/keV << "\t"<< dnorm << G4endl ;
      }
      G4cout << G4endl;
      G4int ii ;
      G4double E1,E2 ;
      E1=-1.e6 ;
      E2=+1.e6 ;
      E = Enlow -dEn ;
      ii = -1;

      for(G4int i1=0; i1<nbinEn; i1++)
      {
        E += dEn ;
        if(ii<0)
        {
          if(distEn[i1] >= 0.5*fmax)
          {
            E1=E ;
            ii=i1 ;
          }
        }
      }
      E = Enlow -dEn ;

      for(G4int i2=0; i2<nbinEn; i2++)
      {
        E += dEn ;

        if(distEn[i2] >= 0.5*fmax)   E2=E ;
      }
      G4cout << " Emp = " << std::setw(15) << Emp/MeV << "   width="
            << std::setw(15) << (E2-E1)/MeV <<   "  MeV " << G4endl;
      G4cout << G4endl ;
    }     
  }
  if(nbinTt>0)
  {
     G4double E , dnorm, norm ,sig;
   G4cout << " transmitted energy distribution " << G4endl ;
   G4cout << "#entries=" << entryTt << "    #underflows=" << underTt <<
             "    #overflows=" << overTt << G4endl ;
   if( entryTt>0.)
   {
     Ttmean /= entryTt;
     sig=Tt2mean/entryTt-Ttmean*Ttmean ;
     if(sig<=0.)
       sig=0.;
     else
       sig=std::sqrt(sig/entryTt) ;
     G4cout << " mean energy of transmitted particles=" << Ttmean/keV << 
               " +- " << sig/keV << "  keV." << G4endl;
     E = Ttlow - dTt ;
     norm = TotNbofEvents*dTt ;
     G4cout << " bin nb      Elow      entries     normalized " << G4endl ;
     for(G4int itt=0; itt<nbinTt; itt++)
     {
      E += dTt ;
      dnorm = distTt[itt]/norm;
      G4cout << std::setw(5) << itt << std::setw(10) << E << 
                std::setw(12) << distTt[itt] <<
                std::setw(12) << dnorm << G4endl ;
     }
     G4cout << G4endl;
   }     
  }
  if(nbinTb>0)
  {
     G4double E , dnorm, norm ,sig;
   G4cout << " backscattered energy distribution " << G4endl ;
   G4cout << "#entries=" << entryTb << "    #underflows=" << underTb <<
             "    #overflows=" << overTb << G4endl ;
   if( entryTb>0.)
   {
     Tbmean /= entryTb;
     sig=Tb2mean/entryTb-Tbmean*Tbmean ;
     if(sig<=0.)
       sig=0.;
     else
       sig=std::sqrt(sig/entryTb) ;
     G4cout << " mean energy of backscattered particles=" << Tbmean/keV << 
               " +- " << sig/keV << "  keV." << G4endl;
     E = Tblow - dTb ;
     norm = TotNbofEvents*dTb ;
     G4cout << " bin nb      Elow      entries     normalized " << G4endl ;
     for(G4int itt=0; itt<nbinTb; itt++)
     {
      E += dTb ;
      dnorm = distTb[itt]/norm;
      G4cout << std::setw(5) << itt << std::setw(10) << E << 
                std::setw(12) << distTb[itt] <<
                std::setw(12) << dnorm << G4endl ;
     }
     G4cout << G4endl;
   }     
  }
  if(nbinTsec>0)
  {G4double E , dnorm, norm ;
   G4cout << " energy distribution of charged secondaries " << G4endl ;
   G4cout << "#entries=" << entryTsec << "    #underflows=" << underTsec <<
             "    #overflows=" << overTsec << G4endl ;
   if( entryTsec>0.)
   {
     E = Tseclow - dTsec ;
     norm = TotNbofEvents*dTsec ;
     G4cout << " bin nb      Elow      entries     normalized " << G4endl ;
     for(G4int itt=0; itt<nbinTsec; itt++)
     {
      E += dTsec ;
      dnorm = distTsec[itt]/norm;
      G4cout << std::setw(5) << itt << std::setw(10) << E << 
                std::setw(12) << distTsec[itt] <<
                std::setw(12) << dnorm << G4endl ;
     }
     G4cout << G4endl;
   }     
  }

  if(nbinR >0)
  {G4double R , dnorm, norm,sig  ;
   G4cout << "  R  distribution " << G4endl ;
   G4cout << "#entries=" << entryR  << "    #underflows=" << underR  <<
             "    #overflows=" << overR  << G4endl ;
   if( entryR >0.)
   {
     Rmean /= entryR;
     sig = R2mean/entryR - Rmean*Rmean;
     if(sig<=0.) sig=0. ;
     else        sig = std::sqrt(sig/entryR) ;
     G4cout << " mean lateral displacement at exit=" << Rmean/mm << " +- "
            << sig/mm << "  mm." << G4endl ; 
     R = Rlow - dR  ;
     norm = TotNbofEvents*dR  ;
     G4cout << " bin nb      Rlow      entries     normalized " << G4endl ;
     for(G4int ier=0; ier<nbinR ; ier++)
     {
      R+= dR  ;
      dnorm = distR[ier]/norm;
      G4cout << std::setw(5) << ier << std::setw(10) << R  <<
                std::setw(12) << distR[ier] <<
                std::setw(12) << dnorm << G4endl ;
     }
     G4cout << G4endl;
   }
  }

  if(nbinTh>0)
  {G4double Th,Thdeg, dnorm, norm,fac0,fnorm,pere,Thpere,Thmean,sum;
   G4cout << "      angle   distribution " << G4endl ;
   G4cout << "#entries=" << entryTh << "    #underflows=" << underTh <<
             "    #overflows=" << overTh << G4endl ;
   if( entryTh>0.)
   {
     Th= Thlow - dTh ;
     norm = TotNbofEvents ;
     if(distTh[0] == 0.)
       fac0 = 1. ;
     else
       fac0 = 1./distTh[0] ;
     pere = 1./std::exp(1.) ;

     G4cout << " bin nb  Thlowdeg      Thlowrad      " <<
               " entries         normalized " << G4endl ;
     Thpere = 0. ;
     sum = 0. ;
     Thmean = 0. ;
     for(G4int ien=0; ien<nbinTh; ien++)
     {
      Th+= dTh ;
      Thdeg = Th*180./pi ;
      sum += distTh[ien] ;
      Thmean += distTh[ien]*(Th+0.5*dTh) ;
      dnorm = distTh[ien]/norm;
      fnorm = fac0*distTh[ien] ;
      if( fnorm > pere)
        Thpere = Th ; 
      G4cout << std::setw(5) << ien << std::setw(10) << Thdeg << "   " <<
                std::setw(10) << Th << "  " <<   
                std::setw(12) << distTh[ien] << "  " <<
                std::setw(12) << dnorm << "  " << std::setw(12) << fnorm <<G4endl ;
     }
     Thmean /= sum ;
     G4cout << G4endl;
     G4cout << " mean = " << Thmean << "  rad  or " << 180.*Thmean/pi <<
               " deg." << G4endl;
     G4cout << " theta(1/e)=" << Thpere << " - " << Thpere+dTh << " rad   "
            << " or " << 180.*Thpere/pi << " - " << 180.*(Thpere+dTh)/pi 
            << " deg." << G4endl;
     G4cout << G4endl;
   }
  }

  if(nbinThback>0)
  {G4double Thb,Thdegb, dnormb, normb,fac0b,fnormb,pereb,Thpereb,Thmeanb,sumb;
   G4cout << " backscattering angle   distribution " << G4endl ;
   G4cout << "#entries=" << entryThback << "    #underflows=" << underThback <<
             "    #overflows=" << overThback << G4endl ;
   if( entryThback>0.)
   {
     Thb= Thlowback - dThback ;
     normb = TotNbofEvents ;
     if(distThback[0] == 0.)
       fac0b = 1. ;
     else
       fac0b = 1./distThback[0] ;
     pereb = 1./std::exp(1.) ;

     G4cout << " bin nb  Thlowdeg      Thlowrad      " <<
               " entries         normalized " << G4endl ;
     Thpereb = 0. ;
     sumb = 0. ;
     Thmeanb = 0. ;
     for(G4int ien=0; ien<nbinThback; ien++)
     {
      Thb+= dThback ;
      Thdegb = Thb*180./pi ;
      sumb += distThback[ien] ;
      Thmeanb += distThback[ien]*(Thb+0.5*dThback) ;
      dnormb = distThback[ien]/normb;
      fnormb = fac0b*distThback[ien] ;
      if( fnormb > pereb)
        Thpereb = Thb ;
      G4cout << std::setw(5) << ien << std::setw(10) << Thdegb << "   " <<
                std::setw(10) << Thb << "  " <<
                std::setw(12) << distThback[ien] << "  " <<
                std::setw(12) << dnormb << "  " << std::setw(12) << fnormb <<G4endl ;
     }
     Thmeanb /= sumb ;
     G4cout << G4endl;
     G4cout << " mean = " << Thmeanb << "  rad  or " << 180.*Thmeanb/pi <<
               " deg." << G4endl;
     G4cout << " theta(1/e)=" << Thpereb << " - " << Thpereb+dThback << " rad   "
            << " or " << 180.*Thpereb/pi << " - " << 180.*(Thpereb+dThback)/pi
            << " deg." << G4endl;
     G4cout << G4endl;
   }
  }

  if(nbinGamma>0)
  {G4double E , fact,dnorm, norm  ;
   G4cout << " gamma energy distribution " << G4endl ;
   G4cout << "#entries=" << entryGamma << "    #underflows=" << underGamma <<
             "    #overflows=" << overGamma << G4endl ;
   if( entryGamma>0.)
   {
     fact=std::exp(dEGamma) ;
     E = ElowGamma/fact  ;
     norm = TotNbofEvents*dEGamma;
     G4cout << " bin nb         Elow      entries       normalized " << G4endl ;
     for(G4int itt=0; itt<nbinGamma; itt++)
     {
      E *= fact ;
      dnorm = distGamma[itt]/norm;
      G4cout << std::setw(5) << itt << std::setw(13) << E << 
                std::setw(12) << distGamma[itt] <<
                std::setw(15) << dnorm << G4endl ;
     }
     G4cout << G4endl;
   }     
  }

  if(nbinvertexz >0)
  {G4double z , dnorm, norm  ;
   G4cout << " vertex Z  distribution " << G4endl ;
   G4cout << "#entries=" << entryvertexz  << "    #underflows=" << undervertexz  <<
             "    #overflows=" << oververtexz  << G4endl ;
   if( entryvertexz >0.)
   {
     z =zlow - dz  ;
     norm = TotNbofEvents*dz  ;
     G4cout << " bin nb      zlow      entries     normalized " << G4endl ;
     for(G4int iez=0; iez<nbinvertexz ; iez++)
     {
      z+= dz  ;
      if(std::fabs(z)<1.e-12) z=0.;
      dnorm = distvertexz[iez]/norm;
      G4cout << std::setw(5) << iez << std::setw(10) << z  <<
                std::setw(12) << distvertexz[iez] <<
                std::setw(12) << dnorm << G4endl ;
     }
     G4cout << G4endl;
   }
  }
  
  G4cout.precision(prec);
  
  if (G4VVisManager::GetConcreteInstance())
  {
    G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/update");
  }

  // save Rndm status

  if (saveRndm == 1)
  { 
    CLHEP::HepRandom::showEngineStatus();
    CLHEP::HepRandom::saveEngineStatus("endOfRun.rndm");
  }     
}

void Em10RunAction::FillEn

(

G4double

En

)

Definition at line 746 of file Em10RunAction.cc.

References plottest35::bin.

Referenced by Em10EventAction::EndOfEventAction().

{
  G4double bin ;
  G4int ibin;

  entryEn += 1. ;
  if(Enlow < Enhigh) {

    if(En < Enlow)          underEn += 1. ;
    else if( En >= Enhigh)  overEn  += 1. ;
    else
    {
      bin = (En-Enlow)/dEn;
      ibin= (G4int)bin;
      if(ibin < 0) { ibin = 0; }
      if(ibin > 199) { ibin = 199; }
      distEn[ibin] += 1.;
    }
  }
}

void Em10RunAction::FillGammaSpectrum

(

G4double

E

)

Definition at line 787 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

void Em10RunAction::FillNbOfSteps

(

G4double

nstep

)

Definition at line 740 of file Em10RunAction.cc.

Referenced by Em10EventAction::EndOfEventAction().

{
}

void Em10RunAction::FillR

(

G4double

R

)

Definition at line 805 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

void Em10RunAction::FillTb

(

G4double

Tt

)

Definition at line 775 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

void Em10RunAction::FillTh

(

G4double

Th

)

Definition at line 793 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

void Em10RunAction::FillThBack

(

G4double

Th

)

Definition at line 799 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

void Em10RunAction::FillTsec

(

G4double

T

)

Definition at line 781 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

void Em10RunAction::FillTt

(

G4double

Tt

)

Definition at line 769 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

void Em10RunAction::Fillvertexz

(

G4double

z

)

Definition at line 811 of file Em10RunAction.cc.

Referenced by Em10SteppingAction::UserSteppingAction().

{
}

G4int Em10RunAction::GetRndmFreq ( )

inline

Definition at line 125 of file Em10RunAction.hh.

Referenced by Em10EventAction::EndOfEventAction().

{return saveRndm;}

void Em10RunAction::SetEhighGamma

(

G4double

Ehigh

)

Definition at line 890 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(Ehigh > 0.0) {
    EhighGamma = Ehigh;
    G4cout << " Ehigh in the gamma spectrum plot = " << EhighGamma << G4endl ;
  }
}

void Em10RunAction::SetElowGamma

(

G4double

Elow

)

Definition at line 882 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(Elow >= 0.0) {
    ElowGamma = Elow;
    G4cout << " Elow  in the gamma spectrum plot = " << ElowGamma << G4endl ;
  }
}

void Em10RunAction::SetEnhigh

(

G4double

Enhigh

)

Definition at line 864 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(Ehigh > 0.0) {
    Enhigh = Ehigh ;
    G4cout << " Ehigh in the  Edep plot = " << Enhigh << G4endl ;
  }
}

void Em10RunAction::SetEnlow

(

G4double

Elow

)

Definition at line 856 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(Elow >= 0.0) {
    Enlow = Elow ;
    G4cout << " Elow  in the  Edep plot = " << Enlow << G4endl ;
  }
}

void Em10RunAction::SethistName

(

G4String

name

)

Definition at line 817 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  histName = name ;
  G4cout << " hist file = " << histName << G4endl;
}

void Em10RunAction::SetnbinEn

(

G4int

nbin

)

Definition at line 848 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinEn = nbin;
    G4cout << " Nb of bins in Edep plot = " << nbinEn << G4endl ;
  }
}

void Em10RunAction::SetnbinGamma

(

G4int

nbin

)

Definition at line 874 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinGamma = nbin;
    G4cout << " Nb of bins in gamma spectrum plot = " << nbinGamma << G4endl ;
  }
}

void Em10RunAction::SetnbinR

(

G4int

nbin

)

Definition at line 963 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinR  = nbin;
    G4cout << " Nb of bins in R plot = " << nbinR << G4endl ;
  }
}

void Em10RunAction::SetnbinStep

(

G4int

nbin

)

Definition at line 823 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin> 0 && nbin<= 200) {
    nbinStep = nbin;
    G4cout << " Nb of bins in #step plot = " << nbinStep << G4endl ;
  }
}

void Em10RunAction::SetnbinTb

(

G4int

nbin

)

Definition at line 922 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinTb = nbin;
    G4cout << " Nb of bins in Ebackscattered plot = " << nbinTb << G4endl ;
  }
}

void Em10RunAction::SetnbinTh

(

G4int

nbin

)

Definition at line 1003 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinTh = nbin;
    G4cout << " Nb of bins in Theta plot = " << nbinTh << G4endl ;
  }
}

void Em10RunAction::SetnbinThBack

(

G4int

nbin

)

Definition at line 1023 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinThback = nbin;
    G4cout << " Nb of bins in Theta plot = " << nbinThback << G4endl ;
  }
}

void Em10RunAction::SetnbinTsec

(

G4int

nbin

)

Definition at line 943 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinTsec = nbin;
    G4cout << " Nb of bins in Tsecondary  plot = " << nbinTsec << G4endl ;
  }
}

void Em10RunAction::SetnbinTt

(

G4int

nbin

)

Definition at line 898 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinTt = nbin;
    G4cout << " Nb of bins in Etransmisssion plot = " << nbinTt << G4endl ;
  }
}

void Em10RunAction::Setnbinzvertex

(

G4int

nbin

)

Definition at line 983 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(nbin > 0 && nbin <= 200) {
    nbinvertexz  = nbin;
    G4cout << " Nb of bins in Z plot = " << nbinvertexz << G4endl ;
  }
}

void Em10RunAction::SetRhigh

(

G4double

Rhigh

)

Definition at line 977 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Rhigh = rhigh ;
  G4cout << " Rhigh in the R plot = " << Rhigh << G4endl ;
}

void Em10RunAction::SetRlow

(

G4double

Rlow

)

Definition at line 971 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Rlow = rlow ;
  G4cout << " Rlow  in the  R plot = " << Rlow << G4endl ;
}

void Em10RunAction::SetRndmFreq ( G4int  val )

inline

Definition at line 124 of file Em10RunAction.hh.

Referenced by Em10RunMessenger::SetNewValue().

{saveRndm = val;}

void Em10RunAction::SetStephigh

(

G4double

Shigh

)

Definition at line 838 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(high > 0.0) {
    Stephigh = high;
    G4cout << " high in the #step plot = " << Stephigh << G4endl ;
  }
}

void Em10RunAction::SetSteplow

(

G4double

Slow

)

Definition at line 831 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(low >= 0.0){ 
    Steplow = low ;
    G4cout << " low  in the #step plot = " << Steplow << G4endl ;
  }
}

void Em10RunAction::SetTbhigh

(

G4double

Tbhigh

)

Definition at line 937 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Tbhigh = Ehigh ;
  G4cout << " Ehigh in the  Ebackscattered plot = " << Tbhigh << G4endl ;
}

void Em10RunAction::SetTblow

(

G4double

Tblow

)

Definition at line 930 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{

  Tblow = Elow ;
  G4cout << " Elow  in the  Ebackscattered plot = " << Tblow << G4endl ;
}

void Em10RunAction::SetThhigh

(

G4double

Thhigh

)

Definition at line 1017 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Thhigh = Thigh ;
  G4cout << " Thigh in the Theta plot = " << Thhigh << G4endl ;
}

void Em10RunAction::SetThhighBack

(

G4double

Thhigh

)

Definition at line 1037 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Thhighback = Thigh ;
  G4cout << " Thigh in the Theta plot = " << Thhighback << G4endl ;
}

void Em10RunAction::SetThlow

(

G4double

Thlow

)

Definition at line 1011 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Thlow = Tlow ;
  G4cout << " Tlow  in the  Theta plot = " << Thlow << G4endl ;
}

void Em10RunAction::SetThlowBack

(

G4double

Thlow

)

Definition at line 1031 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Thlowback = Tlow ;
  G4cout << " Tlow  in the  Theta plot = " << Thlowback << G4endl ;
}

void Em10RunAction::SetTsechigh

(

G4double

Thigh

)

Definition at line 957 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Tsechigh = Ehigh ;
  G4cout << " Ehigh in the  Tsecondary plot = " << Tsechigh << G4endl ;
}

void Em10RunAction::SetTseclow

(

G4double

Tlow

)

Definition at line 951 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  Tseclow = Elow ;
  G4cout << " Elow  in the  Tsecondary plot = " << Tseclow << G4endl ;
}

void Em10RunAction::SetTthigh

(

G4double

Tthigh

)

Definition at line 914 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(Ehigh > 0.0) {
    Tthigh = Ehigh;
    G4cout << " Ehigh in the  Etransmission plot = " << Tthigh << G4endl ;
  }
}

void Em10RunAction::SetTtlow

(

G4double

Ttlow

)

Definition at line 906 of file Em10RunAction.cc.

References G4cout, and G4endl.

Referenced by Em10RunMessenger::SetNewValue().

{
  if(Elow >= 0.0) {
    Ttlow = Elow;
    G4cout << " Elow  in the  Etransmission plot = " << Ttlow << G4endl ;
  }
}

void Em10RunAction::Setzhigh

(

G4double

z

)

Definition at line 997 of file Em10RunAction.cc.

References G4cout, G4endl, and z.

Referenced by Em10RunMessenger::SetNewValue().

{
  zhigh = z ;
  G4cout << " zhigh in the Z plot = " << zhigh << G4endl ;
}

void Em10RunAction::Setzlow

(

G4double

z

)

Definition at line 991 of file Em10RunAction.cc.

References G4cout, G4endl, and z.

Referenced by Em10RunMessenger::SetNewValue().

{
  zlow = z ;
  G4cout << " zlow  in the  Z plot = " << zlow << G4endl ;
}

---

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

• Em10RunAction.hh
• Em10RunAction.cc