WLSRun Class Reference

#include <WLSRun.hh>

Inheritance diagram for WLSRun:

Public Member Functions
	WLSRun ()=default
	~WLSRun () override=default
void	AddTIR (G4int n)
void	AddExiting (G4int n)
void	AddEscapedEnd (G4int n)
void	AddEscapedMid (G4int n)
void	AddBounce (G4int n)
void	AddWLSBounce (G4int n)
void	AddClad1Bounce (G4int n)
void	AddClad2Bounce (G4int n)
void	AddReflected (G4int n)
void	AddEscaped (G4int n)
void	AddMirror (G4int n)
void	AddDetectorHits (G4int n)
void	EndOfRun ()
void	Merge (const G4Run *) override

Private Attributes
G4double	fNTIR = 0.
G4double	fNTIR2 = 0.
G4double	fNExiting = 0.
G4double	fNExiting2 = 0.
G4double	fEscapedEnd = 0.
G4double	fEscapedEnd2 = 0.
G4double	fEscapedMid = 0.
G4double	fEscapedMid2 = 0.
G4double	fBounce = 0.
G4double	fBounce2 = 0.
G4double	fWLSBounce = 0.
G4double	fWLSBounce2 = 0.
G4double	fClad1Bounce = 0.
G4double	fClad1Bounce2 = 0.
G4double	fClad2Bounce = 0.
G4double	fClad2Bounce2 = 0.
G4double	fReflected = 0.
G4double	fReflected2 = 0.
G4double	fEscaped = 0.
G4double	fEscaped2 = 0.
G4double	fMirror = 0.
G4double	fMirror2 = 0.
G4double	fDetectorHits = 0.
G4double	fDetectorHits2 = 0.

Detailed Description

Definition at line 39 of file WLSRun.hh.

Constructor & Destructor Documentation

WLSRun()

WLSRun::WLSRun ( )

default

~WLSRun()

WLSRun::~WLSRun ( )

overridedefault

Member Function Documentation

AddTIR()

void WLSRun::AddTIR ( G4int  n )

inline

Definition at line 45 of file WLSRun.hh.

  {
    G4double nd(n);
    fNTIR += nd;
    fNTIR2 += nd * nd;
  };

AddExiting()

void WLSRun::AddExiting ( G4int  n )

inline

Definition at line 51 of file WLSRun.hh.

  {
    G4double nd(n);
    fNExiting += nd;
    fNExiting2 += nd * nd;
  };

AddEscapedEnd()

void WLSRun::AddEscapedEnd ( G4int  n )

inline

Definition at line 57 of file WLSRun.hh.

  {
    G4double nd(n);
    fEscapedEnd += nd;
    fEscapedEnd2 += nd * nd;
  };

AddEscapedMid()

void WLSRun::AddEscapedMid ( G4int  n )

inline

Definition at line 63 of file WLSRun.hh.

  {
    G4double nd(n);
    fEscapedMid += nd;
    fEscapedMid2 += nd * nd;
  };

AddBounce()

void WLSRun::AddBounce ( G4int  n )

inline

Definition at line 69 of file WLSRun.hh.

  {
    G4double nd(n);
    fBounce += nd;
    fBounce2 += nd * nd;
  };

AddWLSBounce()

void WLSRun::AddWLSBounce ( G4int  n )

inline

Definition at line 75 of file WLSRun.hh.

  {
    G4double nd(n);
    fWLSBounce += nd;
    fWLSBounce2 += nd * nd;
  };

AddClad1Bounce()

void WLSRun::AddClad1Bounce ( G4int  n )

inline

Definition at line 81 of file WLSRun.hh.

  {
    G4double nd(n);
    fClad1Bounce += nd;
    fClad1Bounce2 += nd * nd;
  };

AddClad2Bounce()

void WLSRun::AddClad2Bounce ( G4int  n )

inline

Definition at line 87 of file WLSRun.hh.

  {
    G4double nd(n);
    fClad2Bounce += nd;
    fClad2Bounce2 += nd * nd;
  };

AddReflected()

void WLSRun::AddReflected ( G4int  n )

inline

Definition at line 93 of file WLSRun.hh.

  {
    G4double nd(n);
    fReflected += nd;
    fReflected2 += nd * nd;
  };

AddEscaped()

void WLSRun::AddEscaped ( G4int  n )

inline

Definition at line 99 of file WLSRun.hh.

  {
    G4double nd(n);
    fEscaped += nd;
    fEscaped2 += nd * nd;
  };

AddMirror()

void WLSRun::AddMirror ( G4int  n )

inline

Definition at line 105 of file WLSRun.hh.

  {
    G4double nd(n);
    fMirror += nd;
    fMirror2 += nd * nd;
  };

AddDetectorHits()

void WLSRun::AddDetectorHits ( G4int  n )

inline

Definition at line 111 of file WLSRun.hh.

  {
    G4double nd(n);
    fDetectorHits += nd;
    fDetectorHits2 += nd * nd;
  };

EndOfRun()

void WLSRun::EndOfRun

(

)

Definition at line 70 of file WLSRun.cc.

{
  if(numberOfEvent == 0)
    return;
  auto TotNbofEvents = G4double(numberOfEvent);
 
  fNTIR           = fNTIR / TotNbofEvents;
  fNTIR2          = fNTIR2 / TotNbofEvents;
  G4double rmsTIR = fNTIR2 - fNTIR * fNTIR;
  if(rmsTIR > 0.)
    rmsTIR = std::sqrt(rmsTIR);
  else
    rmsTIR = 0.;
 
  fNExiting           = fNExiting / TotNbofEvents;
  fNExiting2          = fNExiting2 / TotNbofEvents;
  G4double rmsExiting = fNExiting2 - fNExiting * fNExiting;
  if(rmsExiting > 0.)
    rmsExiting = std::sqrt(rmsExiting);
  else
    rmsExiting = 0.;
 
  fEscapedEnd            = fEscapedEnd / TotNbofEvents;
  fEscapedEnd2           = fEscapedEnd2 / TotNbofEvents;
  G4double rmsEscapedEnd = fEscapedEnd2 - fEscapedEnd * fEscapedEnd;
  if(rmsEscapedEnd > 0.)
    rmsEscapedEnd = std::sqrt(rmsEscapedEnd);
  else
    rmsEscapedEnd = 0.;
 
  fEscapedMid            = fEscapedMid / TotNbofEvents;
  fEscapedMid2           = fEscapedMid2 / TotNbofEvents;
  G4double rmsEscapedMid = fEscapedMid2 - fEscapedMid * fEscapedMid;
  if(rmsEscapedMid > 0.)
    rmsEscapedMid = std::sqrt(rmsEscapedMid);
  else
    rmsEscapedMid = 0.;
 
  fBounce            = fBounce / TotNbofEvents;
  fBounce2           = fBounce2 / TotNbofEvents;
  G4double rmsBounce = fBounce2 - fBounce * fBounce;
  if(rmsBounce > 0.)
    rmsBounce = std::sqrt(rmsBounce);
  else
    rmsBounce = 0.;
 
  fWLSBounce            = fWLSBounce / TotNbofEvents;
  fWLSBounce2           = fWLSBounce2 / TotNbofEvents;
  G4double rmsWLSBounce = fWLSBounce2 - fWLSBounce * fWLSBounce;
  if(rmsWLSBounce > 0.)
    rmsWLSBounce = std::sqrt(rmsWLSBounce);
  else
    rmsWLSBounce = 0.;
 
  fClad1Bounce            = fClad1Bounce / TotNbofEvents;
  fClad1Bounce2           = fClad1Bounce2 / TotNbofEvents;
  G4double rmsClad1Bounce = fClad1Bounce2 - fClad1Bounce * fClad1Bounce;
  if(rmsClad1Bounce > 0.)
    rmsClad1Bounce = std::sqrt(rmsClad1Bounce);
  else
    rmsClad1Bounce = 0.;
 
  fClad2Bounce            = fClad2Bounce / TotNbofEvents;
  fClad2Bounce2           = fClad2Bounce2 / TotNbofEvents;
  G4double rmsClad2Bounce = fClad2Bounce2 - fClad2Bounce * fClad2Bounce;
  if(rmsClad2Bounce > 0.)
    rmsClad2Bounce = std::sqrt(rmsClad2Bounce);
  else
    rmsClad2Bounce = 0.;
 
  fReflected            = fReflected / TotNbofEvents;
  fReflected2           = fReflected2 / TotNbofEvents;
  G4double rmsReflected = fReflected2 - fReflected * fReflected;
  if(rmsReflected > 0.)
    rmsReflected = std::sqrt(rmsReflected);
  else
    rmsReflected = 0.;
 
  fEscaped            = fEscaped / TotNbofEvents;
  fEscaped2           = fEscaped2 / TotNbofEvents;
  G4double rmsEscaped = fEscaped2 - fEscaped * fEscaped;
  if(rmsEscaped > 0.)
    rmsEscaped = std::sqrt(rmsEscaped);
  else
    rmsEscaped = 0.;
 
  fMirror            = fMirror / TotNbofEvents;
  fMirror2           = fMirror2 / TotNbofEvents;
  G4double rmsMirror = fMirror2 - fMirror * fMirror;
  if(rmsMirror > 0.)
    rmsMirror = std::sqrt(rmsMirror);
  else
    rmsMirror = 0.;
 
  fDetectorHits            = fDetectorHits / TotNbofEvents;
  fDetectorHits2           = fDetectorHits2 / TotNbofEvents;
  G4double rmsDetectorHits = fDetectorHits2 - fDetectorHits * fDetectorHits;
  if(rmsDetectorHits > 0.)
    rmsDetectorHits = std::sqrt(rmsDetectorHits);
  else
    rmsDetectorHits = 0.;
 
  G4int prec = G4cout.precision(3);
  G4cout << "\n ======================== run summary =====================\n";
  G4cout << "Average number per event:" << G4endl;
  G4cout << " TIR:          " << fNTIR << " +- " << rmsTIR << G4endl
         << " Exiting:      " << fNExiting << " +- " << rmsExiting << G4endl
         << " Escaped Mid:  " << fEscapedMid << " +- " << rmsEscapedMid
         << G4endl << " Escaped End:  " << fEscapedEnd << " +- "
         << rmsEscapedEnd << G4endl << " Bounced:      " << fBounce << " +- "
         << rmsBounce << G4endl << " WLS Bounce:   " << fWLSBounce << " +- "
         << rmsWLSBounce << G4endl << " Clad1 Bounce: " << fClad1Bounce
         << " +- " << rmsClad1Bounce << G4endl
         << " Clad2 Bounce: " << fClad2Bounce << " +- " << rmsClad2Bounce
         << G4endl << " Reflected:    " << fReflected << " +- " << rmsReflected
         << G4endl << " Escaped:      " << fEscaped << " +- " << rmsEscaped
         << G4endl << " Mirror:       " << fMirror << " +- " << rmsMirror
         << G4endl << " Detector hit: " << fDetectorHits << " +- "
         << rmsDetectorHits << G4endl;
 
  G4cout << G4endl;
  G4cout.precision(prec);
}

Merge()

void WLSRun::Merge ( const G4Run *  run )

override

Definition at line 36 of file WLSRun.cc.

{
  const auto localRun = static_cast<const WLSRun*>(run);
 
  fNTIR += localRun->fNTIR;
  fNTIR2 += localRun->fNTIR2;
  fNExiting += localRun->fNExiting;
  fNExiting2 += localRun->fNExiting2;
  fEscapedEnd += localRun->fEscapedEnd;
  fEscapedEnd2 += localRun->fEscapedEnd2;
  fEscapedMid += localRun->fEscapedMid;
  fEscapedMid2 += localRun->fEscapedMid2;
  fBounce += localRun->fBounce;
  fBounce2 += localRun->fBounce2;
  fWLSBounce += localRun->fWLSBounce;
  fWLSBounce2 += localRun->fWLSBounce2;
  fClad1Bounce += localRun->fClad1Bounce;
  fClad1Bounce2 += localRun->fClad1Bounce2;
  fClad2Bounce += localRun->fClad2Bounce;
  fClad2Bounce2 += localRun->fClad2Bounce2;
  fReflected += localRun->fReflected;
  fReflected2 += localRun->fReflected2;
  fEscaped += localRun->fEscaped;
  fEscaped2 += localRun->fEscaped2;
  fMirror += localRun->fMirror;
  fMirror2 += localRun->fMirror2;
  fDetectorHits += localRun->fDetectorHits;
  fDetectorHits2 += localRun->fDetectorHits2;
 
  G4Run::Merge(run);
}

Member Data Documentation

fNTIR

G4double WLSRun::fNTIR = 0.

private

Definition at line 122 of file WLSRun.hh.

fNTIR2

G4double WLSRun::fNTIR2 = 0.

private

Definition at line 123 of file WLSRun.hh.

fNExiting

G4double WLSRun::fNExiting = 0.

private

Definition at line 124 of file WLSRun.hh.

fNExiting2

G4double WLSRun::fNExiting2 = 0.

private

Definition at line 125 of file WLSRun.hh.

fEscapedEnd

G4double WLSRun::fEscapedEnd = 0.

private

Definition at line 126 of file WLSRun.hh.

fEscapedEnd2

G4double WLSRun::fEscapedEnd2 = 0.

private

Definition at line 127 of file WLSRun.hh.

fEscapedMid

G4double WLSRun::fEscapedMid = 0.

private

Definition at line 128 of file WLSRun.hh.

fEscapedMid2

G4double WLSRun::fEscapedMid2 = 0.

private

Definition at line 129 of file WLSRun.hh.

fBounce

G4double WLSRun::fBounce = 0.

private

Definition at line 130 of file WLSRun.hh.

fBounce2

G4double WLSRun::fBounce2 = 0.

private

Definition at line 131 of file WLSRun.hh.

fWLSBounce

G4double WLSRun::fWLSBounce = 0.

private

Definition at line 132 of file WLSRun.hh.

fWLSBounce2

G4double WLSRun::fWLSBounce2 = 0.

private

Definition at line 133 of file WLSRun.hh.

fClad1Bounce

G4double WLSRun::fClad1Bounce = 0.

private

Definition at line 134 of file WLSRun.hh.

fClad1Bounce2

G4double WLSRun::fClad1Bounce2 = 0.

private

Definition at line 135 of file WLSRun.hh.

fClad2Bounce

G4double WLSRun::fClad2Bounce = 0.

private

Definition at line 136 of file WLSRun.hh.

fClad2Bounce2

G4double WLSRun::fClad2Bounce2 = 0.

private

Definition at line 137 of file WLSRun.hh.

fReflected

G4double WLSRun::fReflected = 0.

private

Definition at line 138 of file WLSRun.hh.

fReflected2

G4double WLSRun::fReflected2 = 0.

private

Definition at line 139 of file WLSRun.hh.

fEscaped

G4double WLSRun::fEscaped = 0.

private

Definition at line 140 of file WLSRun.hh.

fEscaped2

G4double WLSRun::fEscaped2 = 0.

private

Definition at line 141 of file WLSRun.hh.

fMirror

G4double WLSRun::fMirror = 0.

private

Definition at line 142 of file WLSRun.hh.

fMirror2

G4double WLSRun::fMirror2 = 0.

private

Definition at line 143 of file WLSRun.hh.

fDetectorHits

G4double WLSRun::fDetectorHits = 0.

private

Definition at line 144 of file WLSRun.hh.

fDetectorHits2

G4double WLSRun::fDetectorHits2 = 0.

private

Definition at line 145 of file WLSRun.hh.

---

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

• extended/optical/wls/include/WLSRun.hh
• extended/optical/wls/src/WLSRun.cc