FCALTBEventAction Class Reference

#include <FCALTBEventAction.hh>

Inheritance diagram for FCALTBEventAction:

Public Member Functions
	FCALTBEventAction (FCALSteppingAction *)
virtual	~FCALTBEventAction ()
virtual void	BeginOfEventAction (const G4Event *)
virtual void	EndOfEventAction (const G4Event *)
void	SetDrawFlag (G4String 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 45 of file FCALTBEventAction.hh.

Constructor & Destructor Documentation

FCALTBEventAction::FCALTBEventAction

(

FCALSteppingAction *

SA

)

Definition at line 60 of file FCALTBEventAction.cc.

  :drawFlag("all"),printModulo(10), StepAction(SA), eventMessenger(0)
{
  eventMessenger = new FCALTBEventActionMessenger(this);
}  

FCALTBEventAction::~FCALTBEventAction ( )

virtual

Definition at line 68 of file FCALTBEventAction.cc.

{
  delete eventMessenger;
  eventMessenger = 0;
}

Member Function Documentation

void FCALTBEventAction::BeginOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 76 of file FCALTBEventAction.cc.

References G4cout, G4endl, G4Event::GetEventID(), and FCALSteppingAction::initialize().

{  
  G4int evtNb = evt->GetEventID();
    if (evtNb%printModulo == 0)
    { 
      G4cout << "\n---> Begin of event: " << evtNb+1 << G4endl;
//      HepRandom::showEngineStatus();
    }
    
  NTracksOutOfWorld = 0;
  NSecondaries = 0;
 
  StepAction->initialize(evtNb+1);
}

void FCALTBEventAction::EndOfEventAction ( const G4Event *  )

virtual

Reimplemented from G4UserEventAction.

Definition at line 93 of file FCALTBEventAction.cc.

References testem0::app, G4VAnalysisManager::FillH1(), G4cout, G4endl, FCALSteppingAction::GetEdepFCAL(), FCALSteppingAction::GetOutOfWorldTracks(), FCALSteppingAction::GetSecondaries(), and python.hepunit::MeV.

{
  // Fill histograms
  G4AnalysisManager* man = G4AnalysisManager::Instance();

  G4int i,j;
  NTracksOutOfWorld = StepAction->GetOutOfWorldTracks(0, 0); 
  G4cout << "N Tracks out of world " << NTracksOutOfWorld << G4endl;

  // Write Leaving Particles in File
  //--------------------------------
  G4String FileName1 = "OutTracks_802_1mm.dat";
  std::ios::openmode iostemp1;
  if(Init1 == 0) {
    iostemp1 = std::ios::out;
    Init1++;
  } else {
    iostemp1 = std::ios::out|std::ios::app; // std::ios::app;  
  };
  std::ofstream OutTracks(FileName1, iostemp1);

  OutTracks << NTracksOutOfWorld << G4endl;

  G4double OutOfWorld;
  for(i=1; i<= NTracksOutOfWorld ; i++){
    for(j=1; j<11 ; j++) {
      //      G4double OutOfWorld = StepAction->GetOutOfWorldTracks(i,j);
      OutOfWorld = StepAction->GetOutOfWorldTracks(i,j);
      OutTracks << OutOfWorld << " " ; 
    }
    OutTracks << std::endl;
    // G4double OutOfWorld2 = StepAction->GetOutOfWorldTracks(i,j);
  } 
  OutTracks.close();

  man->FillH1(1,NTracksOutOfWorld);

  NSecondaries = StepAction->GetSecondaries(0,0);
  G4cout << "N Scondaries " << NSecondaries << G4endl;   
  
  // Write Secondary Particles in File
  //----------------------------------
  G4String FileName2 = "SecndTracks_802_1mm.dat";
  std::ios::openmode iostemp2;
  if(Init2 == 0) {
    iostemp2 = std::ios::out;
    Init2++;
  } else {
    iostemp2 = std::ios::out|std::ios::app; // std::ios::app;  
  };
  
  std::ofstream SecndTracks(FileName2, iostemp2);
  
  SecndTracks << NSecondaries << std::endl;

  G4double Secondary;  
  for(i=1; i<= NSecondaries ; i++){
    for(j=1; j<11 ; j++) {
      Secondary = StepAction->GetSecondaries(i,j);
      SecndTracks << Secondary  << " " ; 
    }
    SecndTracks << std::endl;
    // G4double Secondary2 = StepAction->GetSecondaries(i,j);
  }
  SecndTracks.close();
  
  man->FillH1(2,NSecondaries);

  // Write Edep in FCAL1 and FCAL2 
  G4String FileName3 = "EdepFCAL_802_1mm.dat";
  std::ios::openmode iostemp3;
  if(Init3 == 0) {
    iostemp3 = std::ios::out;
    Init3++;
  } else {
    iostemp3 = std::ios::out|std::ios::app; // std::ios::app;  
  };
  
  std::ofstream EdepFCAL(FileName3, iostemp3);
  
  G4double EmEdep  = StepAction->GetEdepFCAL("FCALEm");
  G4double HadEdep = StepAction->GetEdepFCAL("FCALHad");

  EdepFCAL << EmEdep << " ";
  EdepFCAL << HadEdep; 
  EdepFCAL << std::endl;
  EdepFCAL.close();

  man->FillH1(3,EmEdep/MeV);
  man->FillH1(4,HadEdep/MeV);


  G4cout<<"EmEdep is="<<EmEdep/MeV << " MeV " << G4endl;
  G4cout<<"HadEdep is="<<HadEdep << " MeV" << G4endl;

  G4cout << "Edep in FCAL1 FCAl2 : " << StepAction->GetEdepFCAL("FCALEm") << " ";
  G4cout << StepAction->GetEdepFCAL("FCALHad") << G4endl;
}

void FCALTBEventAction::SetDrawFlag ( G4String  val )

inline

Definition at line 55 of file FCALTBEventAction.hh.

Referenced by FCALTBEventActionMessenger::SetNewValue().

{drawFlag = val;};

void FCALTBEventAction::SetPrintModulo ( G4int  val )

inline

Definition at line 56 of file FCALTBEventAction.hh.

Referenced by FCALTBEventActionMessenger::SetNewValue().

{printModulo = val;};

---

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

• FCALTBEventAction.hh
• FCALTBEventAction.cc