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)

Private Member Functions
void	writeScintHitsToFile ()
void	writePmtHitsToFile (const DMXPmtHitsCollection *)
void	drawTracks (const G4Event *)

Private Attributes
G4int	event_id
const long *	seeds
G4double	energy_pri
G4int	scintillatorCollID
G4int	pmtCollID
G4int	S_hits
G4int	P_hits
G4double	aveTimePmtHits
G4double	totEnergy
G4double	totEnergyGammas
G4double	totEnergyNeutrons
G4double	hitEnergy
G4double	firstLXeHitTime
G4double	firstParticleE
G4double	particleEnergy
G4String	particleName
G4String	firstParticleName
G4bool	gamma_ev
G4bool	neutron_ev
G4bool	positron_ev
G4bool	electron_ev
G4bool	proton_ev
G4bool	other_ev
G4bool	start_gamma
G4bool	start_neutron
G4String	drawTrksFlag
G4String	drawColsFlag
G4int	drawHitsFlag
G4int	savePmtFlag
G4int	saveHitsFlag
G4int	printModulo
DMXEventActionMessenger *	eventMessenger
const DMXRunAction *	runAct
const DMXPrimaryGeneratorAction *	genAction
std::ofstream *	hitsfile
std::ofstream *	pmtfile

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::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::~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

BeginOfEventAction()

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:

drawTracks()

void DMXEventAction::drawTracks ( const G4Event *  evt )

private

Definition at line 484 of file DMXEventAction.cc.

                                                  {

  if(G4VVisManager::GetConcreteInstance()) {
    G4UImanager::GetUIpointer()->ApplyCommand("/vis/scene/notifyHandlers");    
    G4TrajectoryContainer* trajContainer = evt->GetTrajectoryContainer();
    G4int n_trajectories = 0;

    if(trajContainer) n_trajectories = trajContainer->entries();
    for (G4int i=0; i<n_trajectories; i++) {
      G4Trajectory* trj = (G4Trajectory*)(*trajContainer)[i];
      if (drawTrksFlag == "all") 
    trj->DrawTrajectory();
      else if ((drawTrksFlag == "charged") && (trj->GetCharge() != 0.))
    trj->DrawTrajectory();
      else if ((drawTrksFlag == "noscint") 
           && (trj->GetParticleName() != "opticalphoton"))
    trj->DrawTrajectory();
    }
    
    // G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/update");    
  } 

}

Here is the call graph for this function:

Here is the caller graph for this function:

EndOfEventAction()

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:

GetDrawColsFlag()

G4String DMXEventAction::GetDrawColsFlag ( ) const

inline

Definition at line 76 of file DMXEventAction.hh.

{return drawColsFlag;};

Here is the caller graph for this function:

GetDrawTrksFlag()

G4String DMXEventAction::GetDrawTrksFlag ( ) const

inline

Definition at line 73 of file DMXEventAction.hh.

{return drawTrksFlag;};

Here is the caller graph for this function:

SetDrawColsFlag()

void DMXEventAction::SetDrawColsFlag ( G4String  val )

inline

Definition at line 75 of file DMXEventAction.hh.

{drawColsFlag    = val;};

Here is the caller graph for this function:

SetDrawHitsFlag()

void DMXEventAction::SetDrawHitsFlag ( G4int  val )

inline

Definition at line 78 of file DMXEventAction.hh.

{drawHitsFlag    = val;};

Here is the caller graph for this function:

SetDrawTrksFlag()

void DMXEventAction::SetDrawTrksFlag ( G4String  val )

inline

Definition at line 72 of file DMXEventAction.hh.

{drawTrksFlag    = val;};

Here is the caller graph for this function:

SetPrintModulo()

void DMXEventAction::SetPrintModulo ( G4int  val )

inline

Definition at line 81 of file DMXEventAction.hh.

{printModulo     = val;};

Here is the caller graph for this function:

SetSaveHitsFlag()

void DMXEventAction::SetSaveHitsFlag ( G4int  val )

inline

Definition at line 80 of file DMXEventAction.hh.

{saveHitsFlag    = val;};

Here is the caller graph for this function:

SetSavePmtFlag()

void DMXEventAction::SetSavePmtFlag ( G4int  val )

inline

Definition at line 79 of file DMXEventAction.hh.

{savePmtFlag     = val;};

Here is the caller graph for this function:

writePmtHitsToFile()

void DMXEventAction::writePmtHitsToFile ( const DMXPmtHitsCollection *  hits )

private

Definition at line 418 of file DMXEventAction.cc.

{ 
  G4String filename="pmt.out";
  if (runAct)
    filename=runAct->GetsavepmtFile();
  

  //first time it is invoked: create it
  if (!pmtfile)
    {
      //check that we are in a worker: returns -1 in a master and -2 in sequential
      //one file per thread is produced ending with ".N", with N= thread number
      if (G4Threading::G4GetThreadId() >= 0)
    {
      std::stringstream sss;
      sss << filename.c_str() << "." << G4Threading::G4GetThreadId();    
      filename = sss.str();
      //G4cout << "Filename is: " << filename << G4endl;
    }
      pmtfile = new std::ofstream;
      pmtfile->open(filename);
    }


  G4double x; G4double y; G4double z;
  G4AnalysisManager* man = G4AnalysisManager::Instance();

  if(pmtfile->is_open()) {
    (*pmtfile) << "Hit#    X, mm   Y, mm   Z, mm" << G4endl;       
    (*pmtfile) << std::setiosflags(std::ios::fixed)
           << std::setprecision(3)
           << std::setiosflags(std::ios::left)
           << std::setw(6);
    for (G4int i=0; i<P_hits; i++)
      {
    x = ((*hits)[i]->GetPos()).x()/mm;
    y = ((*hits)[i]->GetPos()).y()/mm;
    z = ((*hits)[i]->GetPos()).z()/mm;
    (*pmtfile) << i << "\t"
           << x << "\t" 
           << y << "\t"
           << z << G4endl;
    
    man->FillH2(1,x/mm, y/mm);  // fill(x,y)
    if (event_id == 0 ) {
      man->FillH2(2,x,y);
    }

    //Fill ntuple #3
    man->FillNtupleDColumn(3,0,event_id);
    man->FillNtupleDColumn(3,1,(G4double) i);
    man->FillNtupleDColumn(3,2,x);
    man->FillNtupleDColumn(3,3,y);
    man->FillNtupleDColumn(3,4,z);
    man->AddNtupleRow(3);

      }
    if (event_id%printModulo == 0 && P_hits > 0) 
      G4cout << "     " << P_hits << " PMT hits in " << filename << G4endl;  
  }
  
}

Here is the call graph for this function:

Here is the caller graph for this function:

writeScintHitsToFile()

void DMXEventAction::writeScintHitsToFile ( )

private

Definition at line 307 of file DMXEventAction.cc.

{

  G4String filename="hits.out";
  if (runAct)
    filename=runAct->GetsavehitsFile();

 
  
  //First time it is inkoved
  if (!hitsfile)
    {
      //check that we are in a worker: returns -1 in a master and -2 in sequential
      //one file per thread is produced ending with ".N", with N= thread number
      if (G4Threading::G4GetThreadId() >= 0)
    {
      std::stringstream sss;
      sss << filename.c_str() << "." << G4Threading::G4GetThreadId();    
      filename = sss.str();
      //G4cout << "Filename is: " << filename << G4endl;
    }
      
      hitsfile = new std::ofstream;
      hitsfile->open(filename);
      (*hitsfile) <<"Evt     Eprim   Etot    LXe     LXeTime PMT     PMTTime Seed1           Seed2           First   Flags" 
           << G4endl;
      (*hitsfile) <<"#       MeV     MeV     hits    ns      hits    ns                                      hit"
           << G4endl
           << G4endl;
    }

  if(S_hits) {

    if(hitsfile->is_open()) {


      (*hitsfile) << std::setiosflags(std::ios::fixed)
          << std::setprecision(4)
          << std::setiosflags(std::ios::left)
          << std::setw(6)
          << event_id << "\t"
          << energy_pri/MeV << "\t" 
          << totEnergy/MeV << "\t"
          << S_hits  << "\t"
          << std::setiosflags(std::ios::scientific) 
          << std::setprecision(2)
          << firstLXeHitTime/nanosecond << "\t"
          << P_hits << "\t"
          << std::setiosflags(std::ios::fixed) 
          << std::setprecision(4)
          << aveTimePmtHits/nanosecond << "\t"
          << *seeds     << "\t"
          << *(seeds+1) << "\t"
          << firstParticleName << "\t"
          << (gamma_ev    ? "gamma " : "") 
          << (neutron_ev  ? "neutron " : "") 
          << (positron_ev ? "positron " : "") 
          << (electron_ev ? "electron " : "") 
          << (other_ev    ? "other " : "") 
          << G4endl;

      if (event_id%printModulo == 0)
    G4cout << "     Event summary in file " << filename << G4endl;  
    }

    G4AnalysisManager* man = G4AnalysisManager::Instance();
    G4int firstparticleIndex = 0;
    if(firstParticleName == "gamma") firstparticleIndex = 1;
    else if (firstParticleName == "neutron") firstparticleIndex = 2;
    else if(firstParticleName == "electron") firstparticleIndex = 3;
    else if(firstParticleName == "positron") firstparticleIndex = 4;
    else{
      firstparticleIndex = 5;
      man->FillH1(3,totEnergy);
    }

    man->FillH1(4,P_hits,10); //weight
    man->FillH1(5,P_hits);

    man->FillH1(1,energy_pri/keV);
    man->FillH1(2,totEnergy/keV);
    man->FillH1(6,aveTimePmtHits/ns);

    long seed1 = *seeds;
    long seed2 = *(seeds+1);    
    
    //Fill ntuple #2
    man->FillNtupleDColumn(2,0,event_id);
    man->FillNtupleDColumn(2,1,energy_pri/keV);
    man->FillNtupleDColumn(2,2,totEnergy);
    man->FillNtupleDColumn(2,3,S_hits);
    man->FillNtupleDColumn(2,4,firstLXeHitTime);
    man->FillNtupleDColumn(2,5,P_hits);
    man->FillNtupleDColumn(2,6,aveTimePmtHits);
    man->FillNtupleDColumn(2,7,firstparticleIndex);
    man->FillNtupleDColumn(2,8,firstParticleE);
    man->FillNtupleDColumn(2,9,gamma_ev);
    man->FillNtupleDColumn(2,10,neutron_ev);
    man->FillNtupleDColumn(2,11,positron_ev);
    man->FillNtupleDColumn(2,12,electron_ev);
    man->FillNtupleDColumn(2,13,other_ev);
    man->FillNtupleDColumn(2,14,seed1);
    man->FillNtupleDColumn(2,15,seed2);
    man->AddNtupleRow(2);
    
  }

}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

aveTimePmtHits

G4double DMXEventAction::aveTimePmtHits

private

Definition at line 96 of file DMXEventAction.hh.

drawColsFlag

G4String DMXEventAction::drawColsFlag

private

Definition at line 118 of file DMXEventAction.hh.

drawHitsFlag

G4int DMXEventAction::drawHitsFlag

private

Definition at line 119 of file DMXEventAction.hh.

drawTrksFlag

G4String DMXEventAction::drawTrksFlag

private

Definition at line 117 of file DMXEventAction.hh.

electron_ev

G4bool DMXEventAction::electron_ev

private

Definition at line 110 of file DMXEventAction.hh.

energy_pri

G4double DMXEventAction::energy_pri

private

Definition at line 87 of file DMXEventAction.hh.

event_id

G4int DMXEventAction::event_id

private

Definition at line 81 of file DMXEventAction.hh.

eventMessenger

DMXEventActionMessenger* DMXEventAction::eventMessenger

private

Definition at line 123 of file DMXEventAction.hh.

firstLXeHitTime

G4double DMXEventAction::firstLXeHitTime

private

Definition at line 101 of file DMXEventAction.hh.

firstParticleE

G4double DMXEventAction::firstParticleE

private

Definition at line 102 of file DMXEventAction.hh.

firstParticleName

G4String DMXEventAction::firstParticleName

private

Definition at line 105 of file DMXEventAction.hh.

gamma_ev

G4bool DMXEventAction::gamma_ev

private

Definition at line 107 of file DMXEventAction.hh.

genAction

const DMXPrimaryGeneratorAction* DMXEventAction::genAction

private

Definition at line 126 of file DMXEventAction.hh.

hitEnergy

G4double DMXEventAction::hitEnergy

private

Definition at line 100 of file DMXEventAction.hh.

hitsfile

std::ofstream* DMXEventAction::hitsfile

private

Definition at line 127 of file DMXEventAction.hh.

neutron_ev

G4bool DMXEventAction::neutron_ev

private

Definition at line 108 of file DMXEventAction.hh.

other_ev

G4bool DMXEventAction::other_ev

private

Definition at line 112 of file DMXEventAction.hh.

P_hits

G4int DMXEventAction::P_hits

private

Definition at line 93 of file DMXEventAction.hh.

particleEnergy

G4double DMXEventAction::particleEnergy

private

Definition at line 103 of file DMXEventAction.hh.

particleName

G4String DMXEventAction::particleName

private

Definition at line 104 of file DMXEventAction.hh.

pmtCollID

G4int DMXEventAction::pmtCollID

private

Definition at line 91 of file DMXEventAction.hh.

pmtfile

std::ofstream* DMXEventAction::pmtfile

private

Definition at line 128 of file DMXEventAction.hh.

positron_ev

G4bool DMXEventAction::positron_ev

private

Definition at line 109 of file DMXEventAction.hh.

printModulo

G4int DMXEventAction::printModulo

private

Definition at line 122 of file DMXEventAction.hh.

proton_ev

G4bool DMXEventAction::proton_ev

private

Definition at line 111 of file DMXEventAction.hh.

runAct

const DMXRunAction* DMXEventAction::runAct

private

Definition at line 125 of file DMXEventAction.hh.

S_hits

G4int DMXEventAction::S_hits

private

Definition at line 92 of file DMXEventAction.hh.

saveHitsFlag

G4int DMXEventAction::saveHitsFlag

private

Definition at line 121 of file DMXEventAction.hh.

savePmtFlag

G4int DMXEventAction::savePmtFlag

private

Definition at line 120 of file DMXEventAction.hh.

scintillatorCollID

G4int DMXEventAction::scintillatorCollID

private

Definition at line 90 of file DMXEventAction.hh.

seeds

const long* DMXEventAction::seeds

private

Definition at line 86 of file DMXEventAction.hh.

start_gamma

G4bool DMXEventAction::start_gamma

private

Definition at line 113 of file DMXEventAction.hh.

start_neutron

G4bool DMXEventAction::start_neutron

private

Definition at line 114 of file DMXEventAction.hh.

totEnergy

G4double DMXEventAction::totEnergy

private

Definition at line 97 of file DMXEventAction.hh.

totEnergyGammas

G4double DMXEventAction::totEnergyGammas

private

Definition at line 98 of file DMXEventAction.hh.

totEnergyNeutrons

G4double DMXEventAction::totEnergyNeutrons

private

Definition at line 99 of file DMXEventAction.hh.

---

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

• DMXEventAction.hh
• DMXEventAction.cc