DMXEventAction Class Reference

#include <DMXEventAction.hh>

Inheritance diagram for DMXEventAction:

Collaboration diagram for DMXEventAction:

Public Member Functions
	DMXEventAction ()
virtual	~DMXEventAction ()
virtual void	BeginOfEventAction (const G4Event *)
virtual void	EndOfEventAction (const G4Event *)
void	SetDrawTrksFlag (G4String val)
G4String	GetDrawTrksFlag () const
void	SetDrawColsFlag (G4String val)
G4String	GetDrawColsFlag () const
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 ()
virtual 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

(

)

Definition at line 76 of file DMXEventAction.cc.

  : runAct(0),genAction(0),hitsfile(0),pmtfile(0)
{

  // 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 103 of file DMXEventAction.cc.

                                {

  if (hitsfile)
    {
      hitsfile->close();
      delete hitsfile;
    }
  if (pmtfile)
    {
      pmtfile->close();
      delete pmtfile;
    }
  delete eventMessenger;
}

Member Function Documentation

void DMXEventAction::BeginOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 120 of file DMXEventAction.cc.

{

  //thread-local run action
  if (!runAct) 
    runAct = 
      dynamic_cast<const DMXRunAction*>
      (G4RunManager::GetRunManager()->GetUserRunAction());
  
  if (!genAction)
    genAction = dynamic_cast<const DMXPrimaryGeneratorAction*>
      (G4RunManager::GetRunManager()->GetUserPrimaryGeneratorAction());


  // 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");
  }

}

Here is the call graph for this function:

void DMXEventAction::EndOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 168 of file DMXEventAction.cc.

                                                        {

  // 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(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;  

}

Here is the call graph for this function:

G4String DMXEventAction::GetDrawColsFlag ( ) const

inline

Definition at line 76 of file DMXEventAction.hh.

{return drawColsFlag;};

Here is the caller graph for this function:

G4String DMXEventAction::GetDrawTrksFlag ( ) const

inline

Definition at line 73 of file DMXEventAction.hh.

{return drawTrksFlag;};

Here is the caller graph for this function:

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:

• source/geant4.10.03.p03/examples/advanced/underground_physics/include/DMXEventAction.hh
• source/geant4.10.03.p03/examples/advanced/underground_physics/src/DMXEventAction.cc