LXeEventAction Class Reference

#include <Doxymodules_optical.h>

Inheritance diagram for LXeEventAction:

Public Member Functions
	LXeEventAction (const LXeDetectorConstruction *)
	~LXeEventAction () override
void	BeginOfEventAction (const G4Event *) override
void	EndOfEventAction (const G4Event *) override
void	SetEventVerbose (G4int v)
void	SetPMTThreshold (G4int t)
void	SetForceDrawPhotons (G4bool b)
void	SetForceDrawNoPhotons (G4bool b)
void	IncPhotonCount_Scint ()
void	IncPhotonCount_Ceren ()
void	IncEDep (G4double dep)
void	IncAbsorption ()
void	IncBoundaryAbsorption ()
void	IncHitCount (G4int i=1)
void	SetEWeightPos (const G4ThreeVector &p)
void	SetReconPos (const G4ThreeVector &p)
void	SetConvPos (const G4ThreeVector &p)
void	SetPosMax (const G4ThreeVector &p, G4double edep)
G4int	GetPhotonCount_Scint () const
G4int	GetPhotonCount_Ceren () const
G4int	GetHitCount () const
G4double	GetEDep () const
G4int	GetAbsorptionCount () const
G4int	GetBoundaryAbsorptionCount () const
G4ThreeVector	GetEWeightPos ()
G4ThreeVector	GetReconPos ()
G4ThreeVector	GetConvPos ()
G4ThreeVector	GetPosMax ()
G4double	GetEDepMax ()
G4double	IsConvPosSet ()
G4int	GetPhotonCount ()
void	IncPMTSAboveThreshold ()
G4int	GetPMTSAboveThreshold ()

Private Attributes
LXeEventMessenger *	fEventMessenger = nullptr
const LXeDetectorConstruction *	fDetector = nullptr
G4int	fScintCollID = -1
G4int	fPMTCollID = -1
G4int	fVerbose = 0
G4int	fPMTThreshold = 1
G4bool	fForcedrawphotons = false
G4bool	fForcenophotons = false
G4int	fHitCount = 0
G4int	fPhotonCount_Scint = 0
G4int	fPhotonCount_Ceren = 0
G4int	fAbsorptionCount = 0
G4int	fBoundaryAbsorptionCount = 0
G4double	fTotE = 0.
G4ThreeVector	fEWeightPos
G4ThreeVector	fReconPos
G4ThreeVector	fConvPos
G4bool	fConvPosSet = false
G4ThreeVector	fPosMax
G4double	fEdepMax = 0.
G4int	fPMTsAboveThreshold = 0

Detailed Description

Definition at line 51 of file Doxymodules_optical.h.

Constructor & Destructor Documentation

LXeEventAction()

LXeEventAction::LXeEventAction

(

const LXeDetectorConstruction *

det

)

Definition at line 53 of file LXeEventAction.cc.

  : fDetector(det)
{
  fEventMessenger = new LXeEventMessenger(this);
}

~LXeEventAction()

LXeEventAction::~LXeEventAction ( )

override

Definition at line 61 of file LXeEventAction.cc.

{ delete fEventMessenger; }

Member Function Documentation

BeginOfEventAction()

void LXeEventAction::BeginOfEventAction ( const G4Event *  )

override

Definition at line 65 of file LXeEventAction.cc.

{
  fHitCount                = 0;
  fPhotonCount_Scint       = 0;
  fPhotonCount_Ceren       = 0;
  fAbsorptionCount         = 0;
  fBoundaryAbsorptionCount = 0;
  fTotE                    = 0.0;
 
  fConvPosSet = false;
  fEdepMax    = 0.0;
 
  fPMTsAboveThreshold = 0;
 
  G4SDManager* SDman = G4SDManager::GetSDMpointer();
  if(fScintCollID < 0)
    fScintCollID = SDman->GetCollectionID("scintCollection");
  if(fPMTCollID < 0)
    fPMTCollID = SDman->GetCollectionID("pmtHitCollection");
}

EndOfEventAction()

void LXeEventAction::EndOfEventAction ( const G4Event *  anEvent )

override

Definition at line 88 of file LXeEventAction.cc.

{
  G4TrajectoryContainer* trajectoryContainer =
    anEvent->GetTrajectoryContainer();
 
  G4int n_trajectories = 0;
  if(trajectoryContainer)
    n_trajectories = trajectoryContainer->entries();
 
  // extract the trajectories and draw them
  if(G4VVisManager::GetConcreteInstance())
  {
    for(G4int i = 0; i < n_trajectories; ++i)
    {
      auto trj = (LXeTrajectory*) ((*(anEvent->GetTrajectoryContainer()))[i]);
      if(trj->GetParticleName() == "opticalphoton")
      {
        trj->SetForceDrawTrajectory(fForcedrawphotons);
        trj->SetForceNoDrawTrajectory(fForcenophotons);
      }
      trj->DrawTrajectory();
    }
  }
 
  LXeScintHitsCollection* scintHC = nullptr;
  LXePMTHitsCollection* pmtHC     = nullptr;
  G4HCofThisEvent* hitsCE         = anEvent->GetHCofThisEvent();
 
  // Get the hit collections
  if(hitsCE)
  {
    if(fScintCollID >= 0)
    {
      scintHC = (LXeScintHitsCollection*) (hitsCE->GetHC(fScintCollID));
    }
    if(fPMTCollID >= 0)
    {
      pmtHC = (LXePMTHitsCollection*) (hitsCE->GetHC(fPMTCollID));
    }
  }
 
  // Hits in scintillator
  if(scintHC)
  {
    size_t n_hit = scintHC->entries();
    G4ThreeVector eWeightPos(0.);
    G4double edep;
    G4double edepMax = 0;
 
    for(size_t i = 0; i < n_hit; ++i)
    {  // gather info on hits in scintillator
      edep = (*scintHC)[i]->GetEdep();
      fTotE += edep;
      eWeightPos +=
        (*scintHC)[i]->GetPos() * edep;  // calculate energy weighted pos
      if(edep > edepMax)
      {
        edepMax = edep;  // store max energy deposit
        G4ThreeVector posMax = (*scintHC)[i]->GetPos();
        fPosMax              = posMax;
        fEdepMax             = edep;
      }
    }
 
    G4AnalysisManager::Instance()->FillH1(7, fTotE);
 
    if(fTotE == 0.)
    {
      if(fVerbose > 0)
        G4cout << "No hits in the scintillator this event." << G4endl;
    }
    else
    {
      // Finish calculation of energy weighted position
      eWeightPos /= fTotE;
      fEWeightPos = eWeightPos;
      if(fVerbose > 0)
      {
        G4cout << "\tEnergy weighted position of hits in LXe : "
               << eWeightPos / mm << G4endl;
      }
    }
    if(fVerbose > 0)
    {
      G4cout << "\tTotal energy deposition in scintillator : " << fTotE / keV
             << " (keV)" << G4endl;
    }
  }
 
  if(pmtHC)
  {
    G4ThreeVector reconPos(0., 0., 0.);
    size_t pmts = pmtHC->entries();
    // Gather info from all PMTs
    for(size_t i = 0; i < pmts; ++i)
    {
      fHitCount += (*pmtHC)[i]->GetPhotonCount();
      reconPos += (*pmtHC)[i]->GetPMTPos() * (*pmtHC)[i]->GetPhotonCount();
      if((*pmtHC)[i]->GetPhotonCount() >= fPMTThreshold)
      {
        ++fPMTsAboveThreshold;
      }
      else
      {  // wasn't above the threshold, turn it back off
        (*pmtHC)[i]->SetDrawit(false);
      }
    }
 
    G4AnalysisManager::Instance()->FillH1(1, fHitCount);
    G4AnalysisManager::Instance()->FillH1(2, fPMTsAboveThreshold);
 
    if(fHitCount > 0)
    {  // don't bother unless there were hits
      reconPos /= fHitCount;
      if(fVerbose > 0)
      {
        G4cout << "\tReconstructed position of hits in LXe : " << reconPos / mm
               << G4endl;
      }
      fReconPos = reconPos;
    }
    pmtHC->DrawAllHits();
  }
 
  G4AnalysisManager::Instance()->FillH1(3, fPhotonCount_Scint);
  G4AnalysisManager::Instance()->FillH1(4, fPhotonCount_Ceren);
  G4AnalysisManager::Instance()->FillH1(5, fAbsorptionCount);
  G4AnalysisManager::Instance()->FillH1(6, fBoundaryAbsorptionCount);
 
  if(fVerbose > 0)
  {
    // End of event output. later to be controlled by a verbose level
    G4cout << "\tNumber of photons that hit PMTs in this event : " << fHitCount
           << G4endl;
    G4cout << "\tNumber of PMTs above threshold(" << fPMTThreshold
           << ") : " << fPMTsAboveThreshold << G4endl;
    G4cout << "\tNumber of photons produced by scintillation in this event : "
           << fPhotonCount_Scint << G4endl;
    G4cout << "\tNumber of photons produced by cerenkov in this event : "
           << fPhotonCount_Ceren << G4endl;
    G4cout << "\tNumber of photons absorbed (OpAbsorption) in this event : "
           << fAbsorptionCount << G4endl;
    G4cout << "\tNumber of photons absorbed at boundaries (OpBoundary) in "
           << "this event : " << fBoundaryAbsorptionCount << G4endl;
    G4cout << "Unaccounted for photons in this event : "
           << (fPhotonCount_Scint + fPhotonCount_Ceren - fAbsorptionCount -
               fHitCount - fBoundaryAbsorptionCount)
           << G4endl;
  }
 
  // update the run statistics
  auto run = static_cast<LXeRun*>(
    G4RunManager::GetRunManager()->GetNonConstCurrentRun());
 
  run->IncHitCount(fHitCount);
  run->IncPhotonCount_Scint(fPhotonCount_Scint);
  run->IncPhotonCount_Ceren(fPhotonCount_Ceren);
  run->IncEDep(fTotE);
  run->IncAbsorption(fAbsorptionCount);
  run->IncBoundaryAbsorption(fBoundaryAbsorptionCount);
  run->IncHitsAboveThreshold(fPMTsAboveThreshold);
 
  // If we have set the flag to save 'special' events, save here
  if(fPhotonCount_Scint + fPhotonCount_Ceren < fDetector->GetSaveThreshold())
  {
    G4RunManager::GetRunManager()->rndmSaveThisEvent();
  }
}

SetEventVerbose()

void LXeEventAction::SetEventVerbose ( G4int  v )

inline

Definition at line 52 of file LXeEventAction.hh.

{ fVerbose = v; }

SetPMTThreshold()

void LXeEventAction::SetPMTThreshold ( G4int  t )

inline

Definition at line 54 of file LXeEventAction.hh.

{ fPMTThreshold = t; }

SetForceDrawPhotons()

void LXeEventAction::SetForceDrawPhotons ( G4bool  b )

inline

Definition at line 56 of file LXeEventAction.hh.

{ fForcedrawphotons = b; }

SetForceDrawNoPhotons()

void LXeEventAction::SetForceDrawNoPhotons ( G4bool  b )

inline

Definition at line 57 of file LXeEventAction.hh.

{ fForcenophotons = b; }

IncPhotonCount_Scint()

void LXeEventAction::IncPhotonCount_Scint ( )

inline

Definition at line 59 of file LXeEventAction.hh.

{ ++fPhotonCount_Scint; }

IncPhotonCount_Ceren()

void LXeEventAction::IncPhotonCount_Ceren ( )

inline

Definition at line 60 of file LXeEventAction.hh.

{ ++fPhotonCount_Ceren; }

IncEDep()

void LXeEventAction::IncEDep ( G4double  dep )

inline

Definition at line 61 of file LXeEventAction.hh.

{ fTotE += dep; }

IncAbsorption()

void LXeEventAction::IncAbsorption ( )

inline

Definition at line 62 of file LXeEventAction.hh.

{ ++fAbsorptionCount; }

IncBoundaryAbsorption()

void LXeEventAction::IncBoundaryAbsorption ( )

inline

Definition at line 63 of file LXeEventAction.hh.

{ ++fBoundaryAbsorptionCount; }

IncHitCount()

void LXeEventAction::IncHitCount ( G4int  i = 1 )

inline

Definition at line 64 of file LXeEventAction.hh.

{ fHitCount += i; }

SetEWeightPos()

void LXeEventAction::SetEWeightPos ( const G4ThreeVector &  p )

inline

Definition at line 66 of file LXeEventAction.hh.

{ fEWeightPos = p; }

SetReconPos()

void LXeEventAction::SetReconPos ( const G4ThreeVector &  p )

inline

Definition at line 67 of file LXeEventAction.hh.

{ fReconPos = p; }

SetConvPos()

void LXeEventAction::SetConvPos ( const G4ThreeVector &  p )

inline

Definition at line 68 of file LXeEventAction.hh.

  {
    fConvPos    = p;
    fConvPosSet = true;
  }

SetPosMax()

void LXeEventAction::SetPosMax ( const G4ThreeVector &  p, G4double  edep  )

inline

Definition at line 73 of file LXeEventAction.hh.

  {
    fPosMax  = p;
    fEdepMax = edep;
  }

GetPhotonCount_Scint()

G4int LXeEventAction::GetPhotonCount_Scint ( ) const

inline

Definition at line 79 of file LXeEventAction.hh.

{ return fPhotonCount_Scint; }

GetPhotonCount_Ceren()

G4int LXeEventAction::GetPhotonCount_Ceren ( ) const

inline

Definition at line 80 of file LXeEventAction.hh.

{ return fPhotonCount_Ceren; }

GetHitCount()

G4int LXeEventAction::GetHitCount ( ) const

inline

Definition at line 81 of file LXeEventAction.hh.

{ return fHitCount; }

GetEDep()

G4double LXeEventAction::GetEDep ( ) const

inline

Definition at line 82 of file LXeEventAction.hh.

{ return fTotE; }

GetAbsorptionCount()

G4int LXeEventAction::GetAbsorptionCount ( ) const

inline

Definition at line 83 of file LXeEventAction.hh.

{ return fAbsorptionCount; }

GetBoundaryAbsorptionCount()

G4int LXeEventAction::GetBoundaryAbsorptionCount ( ) const

inline

Definition at line 84 of file LXeEventAction.hh.

{ return fBoundaryAbsorptionCount; }

GetEWeightPos()

G4ThreeVector LXeEventAction::GetEWeightPos ( )

inline

Definition at line 86 of file LXeEventAction.hh.

{ return fEWeightPos; }

GetReconPos()

G4ThreeVector LXeEventAction::GetReconPos ( )

inline

Definition at line 87 of file LXeEventAction.hh.

{ return fReconPos; }

GetConvPos()

G4ThreeVector LXeEventAction::GetConvPos ( )

inline

Definition at line 88 of file LXeEventAction.hh.

{ return fConvPos; }

GetPosMax()

G4ThreeVector LXeEventAction::GetPosMax ( )

inline

Definition at line 89 of file LXeEventAction.hh.

{ return fPosMax; }

GetEDepMax()

G4double LXeEventAction::GetEDepMax ( )

inline

Definition at line 90 of file LXeEventAction.hh.

{ return fEdepMax; }

IsConvPosSet()

G4double LXeEventAction::IsConvPosSet ( )

inline

Definition at line 91 of file LXeEventAction.hh.

{ return fConvPosSet; }

GetPhotonCount()

G4int LXeEventAction::GetPhotonCount ( )

inline

Definition at line 94 of file LXeEventAction.hh.

{ return fPhotonCount_Scint + fPhotonCount_Ceren; }

IncPMTSAboveThreshold()

void LXeEventAction::IncPMTSAboveThreshold ( )

inline

Definition at line 96 of file LXeEventAction.hh.

{ ++fPMTsAboveThreshold; }

GetPMTSAboveThreshold()

G4int LXeEventAction::GetPMTSAboveThreshold ( )

inline

Definition at line 97 of file LXeEventAction.hh.

{ return fPMTsAboveThreshold; }

Member Data Documentation

fEventMessenger

LXeEventMessenger* LXeEventAction::fEventMessenger = nullptr

private

Definition at line 100 of file LXeEventAction.hh.

fDetector

const LXeDetectorConstruction* LXeEventAction::fDetector = nullptr

private

Definition at line 101 of file LXeEventAction.hh.

fScintCollID

G4int LXeEventAction::fScintCollID = -1

private

Definition at line 103 of file LXeEventAction.hh.

fPMTCollID

G4int LXeEventAction::fPMTCollID = -1

private

Definition at line 104 of file LXeEventAction.hh.

fVerbose

G4int LXeEventAction::fVerbose = 0

private

Definition at line 106 of file LXeEventAction.hh.

fPMTThreshold

G4int LXeEventAction::fPMTThreshold = 1

private

Definition at line 108 of file LXeEventAction.hh.

fForcedrawphotons

G4bool LXeEventAction::fForcedrawphotons = false

private

Definition at line 110 of file LXeEventAction.hh.

fForcenophotons

G4bool LXeEventAction::fForcenophotons = false

private

Definition at line 111 of file LXeEventAction.hh.

fHitCount

G4int LXeEventAction::fHitCount = 0

private

Definition at line 113 of file LXeEventAction.hh.

fPhotonCount_Scint

G4int LXeEventAction::fPhotonCount_Scint = 0

private

Definition at line 114 of file LXeEventAction.hh.

fPhotonCount_Ceren

G4int LXeEventAction::fPhotonCount_Ceren = 0

private

Definition at line 115 of file LXeEventAction.hh.

fAbsorptionCount

G4int LXeEventAction::fAbsorptionCount = 0

private

Definition at line 116 of file LXeEventAction.hh.

fBoundaryAbsorptionCount

G4int LXeEventAction::fBoundaryAbsorptionCount = 0

private

Definition at line 117 of file LXeEventAction.hh.

fTotE

G4double LXeEventAction::fTotE = 0.

private

Definition at line 119 of file LXeEventAction.hh.

fEWeightPos

G4ThreeVector LXeEventAction::fEWeightPos

private

Definition at line 123 of file LXeEventAction.hh.

fReconPos

G4ThreeVector LXeEventAction::fReconPos

private

Definition at line 124 of file LXeEventAction.hh.

fConvPos

G4ThreeVector LXeEventAction::fConvPos

private

Definition at line 125 of file LXeEventAction.hh.

fConvPosSet

G4bool LXeEventAction::fConvPosSet = false

private

Definition at line 126 of file LXeEventAction.hh.

fPosMax

G4ThreeVector LXeEventAction::fPosMax

private

Definition at line 127 of file LXeEventAction.hh.

fEdepMax

G4double LXeEventAction::fEdepMax = 0.

private

Definition at line 128 of file LXeEventAction.hh.

fPMTsAboveThreshold

G4int LXeEventAction::fPMTsAboveThreshold = 0

private

Definition at line 130 of file LXeEventAction.hh.

---

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

• .doxygen/Doxymodules_optical.h
• extended/optical/LXe/include/LXeEventAction.hh
• extended/optical/LXe/src/LXeEventAction.cc