DMXEventAction Class Reference

#include <DMXEventAction.hh>

Inheritance diagram for DMXEventAction:

Public Member Functions
	DMXEventAction (DMXRunAction *, DMXPrimaryGeneratorAction *)
virtual	~DMXEventAction ()
virtual void	BeginOfEventAction (const G4Event *)
virtual void	EndOfEventAction (const G4Event *)
void	SetDrawTrksFlag (G4String val)
G4String	GetDrawTrksFlag ()
void	SetDrawColsFlag (G4String val)
G4String	GetDrawColsFlag ()
void	SetDrawHitsFlag (G4int val)
void	SetSavePmtFlag (G4int val)
void	SetSaveHitsFlag (G4int val)
void	SetPrintModulo (G4int val)
Public Member Functions inherited from G4UserEventAction
	G4UserEventAction ()
virtual	~G4UserEventAction ()
void	SetEventManager (G4EventManager *value)

Additional Inherited Members
Protected Attributes inherited from G4UserEventAction
G4EventManager *	fpEventManager

Detailed Description

Definition at line 58 of file DMXEventAction.hh.

Constructor & Destructor Documentation

DMXEventAction::DMXEventAction	(	DMXRunAction *	DMXRun,
		DMXPrimaryGeneratorAction *	DMXGenerator
	)

Definition at line 74 of file DMXEventAction.cc.

  : runAct(DMXRun),genAction(DMXGenerator)
{

  // create messenger
  eventMessenger = new DMXEventActionMessenger(this);

  // defaults for messenger
  drawColsFlag = "standard";
  drawTrksFlag = "all";
  drawHitsFlag = 1;
  savePmtFlag  = 0;
  saveHitsFlag = 1;
  
  printModulo = 1;

  // hits collections
  scintillatorCollID = -1;
  pmtCollID = -1;

  energy_pri=0;
  seeds=NULL;

}

DMXEventAction::~DMXEventAction ( )

virtual

Definition at line 102 of file DMXEventAction.cc.

                                {

  delete eventMessenger;
}

Member Function Documentation

void DMXEventAction::BeginOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 109 of file DMXEventAction.cc.

References G4BestUnit, G4cout, G4endl, G4SDManager::GetCollectionID(), DMXPrimaryGeneratorAction::GetEnergyPrimary(), G4Event::GetEventID(), DMXPrimaryGeneratorAction::GetEventSeeds(), and G4SDManager::GetSDMpointer().

                                                          {


  // grab seeds
  seeds = genAction->GetEventSeeds();

  // grab energy of primary
  energy_pri = genAction->GetEnergyPrimary();

  event_id = evt->GetEventID();
 
  // print this information event by event (modulo n)   
  if (event_id%printModulo == 0)
    {
      G4cout << "\n---> Begin of event: " << event_id << G4endl;
      G4cout << "\n     Primary Energy: " << G4BestUnit(energy_pri,"Energy") 
         << G4endl;
      //      HepRandom::showEngineStatus(); 
    }


  // get ID for scintillator hits collection
  if (scintillatorCollID==-1) {
    G4SDManager *SDman = G4SDManager::GetSDMpointer();
    scintillatorCollID = SDman->GetCollectionID("scintillatorCollection");
  }

  // get ID for pmt hits collection
  if (pmtCollID==-1) {
    G4SDManager *SDman = G4SDManager::GetSDMpointer();
    pmtCollID = SDman->GetCollectionID("pmtCollection");
  }

}

void DMXEventAction::EndOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 146 of file DMXEventAction.cc.

References G4THitsCollection< T >::DrawAllHits(), G4THitsCollection< T >::entries(), G4VAnalysisManager::FillH1(), G4BestUnit, G4cout, G4endl, G4HCofThisEvent::GetHC(), and G4Event::GetHCofThisEvent().

                                                        {

  // check that both hits collections have been defined
  if(scintillatorCollID<0||pmtCollID<0) return;

  G4AnalysisManager* man = G4AnalysisManager::Instance();

  // address hits collections
  DMXScintHitsCollection* SHC = NULL;
  DMXPmtHitsCollection*   PHC = NULL;
  G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
  if(HCE) {
    SHC = (DMXScintHitsCollection*)(HCE->GetHC(scintillatorCollID));
    PHC = (DMXPmtHitsCollection*)(HCE->GetHC(pmtCollID));
  }

  // event summary
  totEnergy         = 0.;
  totEnergyGammas   = 0.;
  totEnergyNeutrons = 0.;
  firstParticleE    = 0.;
  particleEnergy    = 0.;
  firstLXeHitTime   = 0.;
  aveTimePmtHits    = 0.;

  firstParticleName = "";
  particleName      = "";


  // particle flags
  gamma_ev          = false;
  neutron_ev        = false;
  positron_ev       = false;
  electron_ev       = false;
  proton_ev         = false;
  other_ev          = false;
  start_gamma       = false;
  start_neutron     = false;


  // scintillator hits
  if(SHC) {
    S_hits = SHC->entries();
    
    for (G4int i=0; i<S_hits; i++) {
      if(i==0) {
    firstParticleName = (*SHC)[0]->GetParticle();
    firstLXeHitTime   = (*SHC)[0]->GetTime();
    firstParticleE = (*SHC)[0]->GetParticleEnergy();
    if (event_id%printModulo == 0 && S_hits > 0) {
      G4cout << "     First hit in LXe: " << firstParticleName << G4endl;
      G4cout << "     Number of hits in LXe: " << S_hits << G4endl; 
    }
      }
      hitEnergy         = (*SHC)[i]->GetEdep();
      totEnergy        += hitEnergy;
      
      particleName      = (*SHC)[i]->GetParticle();
      particleEnergy    = (*SHC)[i]->GetParticleEnergy();

      if(particleName == "gamma") {
    gamma_ev = true;
    start_gamma = true;
    start_neutron = false;
      }
      else if(particleName == "neutron") 
    neutron_ev = true;
      else if(particleName == "e+") 
    positron_ev = true;
      else if(particleName == "e-") 
    electron_ev = true;
      else if(particleName == "proton") 
    proton_ev = true;
      else {
    other_ev = true;
    start_gamma = false;
    start_neutron = true;
      }

      if(start_gamma && !start_neutron) 
    totEnergyGammas += hitEnergy;
      if(start_neutron && !start_gamma) 
    totEnergyNeutrons += hitEnergy;
    }
    
    if (event_id%printModulo == 0)
      G4cout << "     Total energy in LXe: " 
         << G4BestUnit(totEnergy,"Energy") << G4endl;  
    
  }
  
  
  // PMT hits
  if(PHC) {
    P_hits = PHC->entries();
    
    // average time of PMT hits
    for (G4int i=0; i<P_hits; i++) {
      G4double time = ( (*PHC)[i]->GetTime() - firstLXeHitTime );
      aveTimePmtHits += time / (G4double)P_hits;
      ////      if (event_id == 0) {
      if(P_hits >= 0) {
    man->FillH1(7,time);
      }
    }
  
    if (event_id%printModulo == 0 && P_hits > 0) {
      G4cout << "     Average light collection time: "
         << G4BestUnit(aveTimePmtHits,"Time") << G4endl;
      G4cout << "     Number of PMT hits (photoelectron equivalent): " 
         << P_hits << G4endl;     
    }
    // write out (x,y,z) of PMT hits
    if (savePmtFlag)
      writePmtHitsToFile(PHC);
  }
  

  // write out event summary
  if(saveHitsFlag) 
    writeScintHitsToFile();
  
  // draw trajectories
  if(drawColsFlag=="standard" && drawTrksFlag!="none")
    drawTracks(evt);

  // hits in PMT
  if(drawHitsFlag && PHC) 
    PHC->DrawAllHits();

  // print this event by event (modulo n)   
  if (event_id%printModulo == 0) 
    G4cout << "\n---> End of event: " << event_id << G4endl << G4endl;  

}

G4String DMXEventAction::GetDrawColsFlag ( )

inline

Definition at line 76 of file DMXEventAction.hh.

Referenced by DMXSteppingAction::UserSteppingAction().

{return drawColsFlag;};

G4String DMXEventAction::GetDrawTrksFlag ( )

inline

Definition at line 73 of file DMXEventAction.hh.

Referenced by DMXSteppingAction::UserSteppingAction().

{return drawTrksFlag;};

void DMXEventAction::SetDrawColsFlag ( G4String  val )

inline

Definition at line 75 of file DMXEventAction.hh.

{drawColsFlag    = val;};

void DMXEventAction::SetDrawHitsFlag ( G4int  val )

inline

Definition at line 78 of file DMXEventAction.hh.

{drawHitsFlag    = val;};

void DMXEventAction::SetDrawTrksFlag ( G4String  val )

inline

Definition at line 72 of file DMXEventAction.hh.

{drawTrksFlag    = val;};

void DMXEventAction::SetPrintModulo ( G4int  val )

inline

Definition at line 81 of file DMXEventAction.hh.

{printModulo     = val;};

void DMXEventAction::SetSaveHitsFlag ( G4int  val )

inline

Definition at line 80 of file DMXEventAction.hh.

{saveHitsFlag    = val;};

void DMXEventAction::SetSavePmtFlag ( G4int  val )

inline

Definition at line 79 of file DMXEventAction.hh.

{savePmtFlag     = val;};

---

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

• DMXEventAction.hh
• DMXEventAction.cc