LBNEEventAction.cc

Go to the documentation of this file.

//
//LBNEEventAction.cc
//

//C++
#include <string>
#include <stdlib.h>

#include "LBNEEventAction.hh"
#include "LBNERunManager.hh"
#include "LBNEAnalysis.hh"
#include "LBNETrajectory.hh"

#include "G4RunManager.hh"
#include "G4Run.hh"
#include "G4Event.hh"
#include "G4Track.hh"
#include "G4EventManager.hh"
#include "G4TrajectoryContainer.hh"
#include "G4Trajectory.hh"
#include "G4ios.hh"
#include "LBNEQuickPiToNu.hh"

//------------------------------------------------------------------------------------- 
LBNEEventAction::LBNEEventAction()
{
   pRunManager=(LBNERunManager*)LBNERunManager::GetRunManager();

   if(pRunManager->GetVerboseLevel() > 0)
   {
      std::cout << "LBNEEventAction Constructor Called." << std::endl;
   }
   
}

//-------------------------------------------------------------------------------------
LBNEEventAction::~LBNEEventAction()
{

   if(pRunManager->GetVerboseLevel() > 0)
   {
      std::cout << "LBNEEventAction Destructor Called." << std::endl;
   }
}

//-------------------------------------------------------------------------------------
void LBNEEventAction::BeginOfEventAction(const G4Event* evt)
{
   if(pRunManager->GetVerboseLevel() > 1) 
   { 
      G4cout << "LBNEEventAction::BeginOfEventAction called...Beginning Event #" << evt->GetEventID() << G4endl;
   }
   const LBNERunAction *runUserAction =  reinterpret_cast<const LBNERunAction *>(G4RunManager::GetRunManager()->GetUserRunAction());
   if (runUserAction->GetDoComputeEDepInGraphite()) {
       LBNESteppingAction *stepAction = (LBNESteppingAction *) G4RunManager::GetRunManager()->GetUserSteppingAction();
       stepAction->ResetEnergyDepInGraphite();
//       std::cerr << " Ready to tally energy deposition in the graphite ... " << std::endl;
   }
   if (runUserAction->GetDoComputeEDepInArgonGas()) {
       LBNESteppingAction *stepAction = (LBNESteppingAction *) G4RunManager::GetRunManager()->GetUserSteppingAction();
       stepAction->ResetEnergyDepInArgonGas();
//       std::cerr << " Ready to tally energy deposition in the graphite ... " << std::endl;
   }
   // This code does not work.. Need to be studied.. 
//   if (evt->GetEventID() == -2222) {
//    
//     CLHEP::HepRandom::restoreEngineStatus("./rndm/beginRun_103Evt_32135.rndm");
//     std::cerr << " Reloaded random from 103Evt_32135 " << std::endl;
//   }
// 
//   if (evt->GetEventID() > 32130) {
//     const G4Run* theRun = pRunManager->GetCurrentRun();
//     std::ostringstream fOutRndmStrStr; 
//     fOutRndmStrStr << "./rndm/beginRun_" << theRun->GetRunID() << "Evt_" << evt->GetEventID() << ".rndm";
//     std::string fOutRndmStr(fOutRndmStrStr.str());
//     CLHEP::HepRandom::saveEngineStatus(fOutRndmStr.c_str());
     // delete the previous previous file 
//     std::ostringstream fOutRndmPStrStr; 
//     fOutRndmPStrStr << "rm -f ./rndm/beginRun_" << theRun->GetRunID() << "Evt_" << evt->GetEventID()-2 << ".rndm";
//     std::string fOutRndmPStr(fOutRndmPStrStr.str());
//     system(fOutRndmPStr.c_str());
     // Also get the seeds May this what we have to save... 
//     std::cerr << " The seed is now " << CLHEP::HepRandom::getTheSeed() << std::endl;
//   }

  if (pRunManager->GetCreateTrkPlaneOutput()||pRunManager->GetCreateAlcoveTrackingOutput()){
    //--for tracking planes
    LBNEAnalysis* analysis = LBNEAnalysis::getInstance();
    analysis->ResetEvent();
  }

}


//-------------------------------------------------------------------------------------
void LBNEEventAction::EndOfEventAction(const G4Event* evt)
{
  if(pRunManager->GetVerboseLevel() > 1) 
    { 
      G4cout << "LBNEEventAction::EndOfEventAction called...Ending Event # " << evt-> GetEventID() << G4endl;
    }
   if( LBNEQuickPiToNuVect::Instance()->doIt()) LBNEQuickPiToNuVect::Instance()->evtOut();

   const LBNERunAction *runUserAction =  reinterpret_cast<const LBNERunAction *>(G4RunManager::GetRunManager()->GetUserRunAction());
   if (runUserAction->GetDoComputeEDepInGraphite()) {
     const LBNESteppingAction *stepAction = 
        reinterpret_cast<const LBNESteppingAction *> 
                (G4RunManager::GetRunManager()->GetUserSteppingAction());
     runUserAction->storeEDepInGraphite(stepAction->GetEnergyDepInGraphite());
   }
   if (runUserAction->GetDoComputeEDepInArgonGas()) {
     const LBNESteppingAction *stepAction = 
        reinterpret_cast<const LBNESteppingAction *> 
                (G4RunManager::GetRunManager()->GetUserSteppingAction());
     const double eDepH1 = stepAction->GetEnergyDepInArgonGasH1();
     const double eDepH2 = stepAction->GetEnergyDepInArgonGasH2();
     runUserAction->storeEDepInArgonGas(eDepH1, eDepH2);
   }

   //--for tracking planes
   if (pRunManager->GetCreateTrkPlaneOutput()||pRunManager->GetCreateAlcoveTrackingOutput()){
    LBNEAnalysis* analysis = LBNEAnalysis::getInstance();
    analysis->FillEvent();
   }

}