LBNERunAction Class Reference

#include <LBNERunAction.hh>

Inheritance diagram for LBNERunAction:

Public Member Functions
	LBNERunAction ()
	~LBNERunAction ()
void	BeginOfRunAction (const G4Run *)
void	EndOfRunAction (const G4Run *)
void	CheckOKToRun ()
void	CheckMemoryUsage (double VMLimit) const
void	SetDebugLevel (bool t)
bool	GetDebugLevel () const
void	SetDoComputeEDepInGraphite (bool t) const
bool	GetDoComputeEDepInGraphite () const
void	storeEDepInGraphite (double edep) const
void	SetDoComputeEDepInArgonGas (bool t) const
bool	GetDoComputeEDepInArgonGas () const
void	storeEDepInArgonGas (double edepH1, double edepH2) const
void	SetDoComputeEDepInHorns (bool t) const
bool	GetDoComputeEDepInHorns () const
LBNFDeDxMap *	GetG4MARSDeDxMap () const

Private Attributes
int	fVerboseLevel
LBNERunActionMessenger *	runMessenger
bool	fDoComputeEDepInGraphite
bool	fDoComputeEDepInArgonGas
double	fMeanEDepInGraphite
double	fRMSEDepInGraphite
double	fMeanEDepInArgonGasHorn1
double	fRMSEDepInArgonGasHorn1
double	fMeanEDepInArgonGasHorn2
double	fRMSEDepInArgonGasHorn2
bool	fDoComputeEDepInHorns
LBNFDeDxMap *	fG4MARSDeDxMap

Detailed Description

Definition at line 16 of file LBNERunAction.hh.

Constructor & Destructor Documentation

LBNERunAction::LBNERunAction

(

)

Definition at line 28 of file LBNERunAction.cc.

                            :
fVerboseLevel(1),
fDoComputeEDepInGraphite(false),
fDoComputeEDepInArgonGas(false),
fMeanEDepInGraphite(0.),
fRMSEDepInGraphite(0.),
fMeanEDepInArgonGasHorn1(0.),
fRMSEDepInArgonGasHorn1(0.),
fMeanEDepInArgonGasHorn2(0.),
fRMSEDepInArgonGasHorn2(0.),
fDoComputeEDepInHorns(false),
fG4MARSDeDxMap(0)
{
   
   if(fVerboseLevel > 0)
   {
      std::cout << "LBNERunAction Constructor Called." << std::endl;
   }

   
   runMessenger = new LBNERunActionMessenger(this);
}

LBNERunAction::~LBNERunAction

(

)

Definition at line 52 of file LBNERunAction.cc.

{ 
   if(fVerboseLevel > 0)
   {
      std::cout << "LBNERunAction Destructor Called." << std::endl;
   }

  delete runMessenger;
}

Member Function Documentation

void LBNERunAction::BeginOfRunAction

(

const G4Run *

aRun

)

Definition at line 63 of file LBNERunAction.cc.

{
    if(fVerboseLevel > 0)
   {
      std::cout << "LBNERunAction::BeginOfRunAction() Called." << std::endl;
   }

   
   LBNERunManager* theRunManager = dynamic_cast<LBNERunManager*>(G4RunManager::GetRunManager());

   G4String spaces = "   ";
   //std::cout << std::endl;
   std::cout << "LBNERunAction::BeginOfRunAction() - Initializing Run " 
             << aRun->GetRunID() << "..." << std::endl;
   
   //std::cout << spaces << "Preparing \"" <<  fLBNEData -> GetSimulation() << "\" Simulation." << std::endl;

   //
   //Random number generator
   //
   G4String randomFile="rndm/beginOfRun_";
   char runN[4];
   snprintf(runN, 4, "%03d",aRun->GetRunID());
   randomFile.append(runN);
   randomFile.append(".rndm");
   CLHEP::HepRandom::saveEngineStatus(randomFile);
   std::cout << spaces << "Intializing Random Number generator named  " 
            << CLHEP::HepRandom::getTheEngine()->name() << " ... " <<  std::endl
             << spaces << "   Seed = " << CLHEP::HepRandom::getTheSeed() << std::endl
             << spaces << "   Saving Random engine status in "<< randomFile << std::endl;
   //

   fMeanEDepInGraphite = 0.;
   fRMSEDepInGraphite = 0.;
   fMeanEDepInArgonGasHorn1 = 0.;
   fRMSEDepInArgonGasHorn1 = 0.;
   fMeanEDepInArgonGasHorn2 = 0.;
   fRMSEDepInArgonGasHorn2 = 0.;
   //
   //Open input ntuples set Number of Events to loop over
   //
//   LBNEPrimaryGeneratorAction* primaryGeneratorAction = dynamic_cast<LBNEPrimaryGeneratorAction*>(pRunManager->userPrimaryGeneratorAction);
//   const LBNEPrimaryGeneratorAction* primaryGeneratorAction = 
//      dynamic_cast<const LBNEPrimaryGeneratorAction*>(pRunManager->GetUserPrimaryGeneratorAction());

   LBNEPrimaryGeneratorAction* primaryGeneratorAction = 
      (theRunManager->GetLBNEPrimaryGeneratorAction());

   
   if(!primaryGeneratorAction)
   {
      std::cout << spaces << "PROBLEM: INVALID LBNEPrimaryGeneratorAction POINTER" << std::endl;
//      fLBNEData -> SetOKToRun(false);
      return;
   }

//   if(fLBNEData->GetSimulation() == "Standard Neutrino Beam" ||
//      fLBNEData->GetSimulation() == "Target Tracking"        ||
//      fLBNEData->GetSimulation() == "Horn 1 Tracking"    ||
//      fLBNEData->GetSimulation() == "Horn 2 Tracking" )
   {

      if (theRunManager->GetUseFlukaInput() || theRunManager->GetUseMarsInput()) 
      {
         G4bool ntpOpen = primaryGeneratorAction->OpenNtuple(theRunManager->GetNptInputFileName());
//       
         if(!ntpOpen)
         {
            std::ostringstream mStrStr;
            mStrStr << "PROBLEM: FAILED TO OPEN INPUT NTUPLE. Fatal " << std::endl;
            G4String mStr(mStrStr.str());
            G4Exception("LBNERunAction::BeginOfRunAction", " ", RunMustBeAborted, mStr.c_str()); 
            
            return;
         } 
//       
         std::cout << spaces << "Successfully opened input ntuple \"" 
                   << theRunManager->GetNptInputFileName() << "\"" << std::endl;
//       
        theRunManager->SetNumberOfEventsLBNE(primaryGeneratorAction->GetNoOfPrimaries());
        // Note there are public accessors in geant4 v4.9.6 inside the G4RunManager... but not in v4.9.3... 
//       
      }
      else
      {
         primaryGeneratorAction->SetProtonBeam();

      }
      
   }
   //

   //
   //If still ok to run open output file
   //

   LBNEAnalysis* analysis = LBNEAnalysis::getInstance();
   if(theRunManager -> GetCreateOutput()) 
      if(!(analysis->CreateOutput())) {
            std::ostringstream mStrStr;
            mStrStr << "PROBLEM: FAILED TO OPEN OUTPUT NTUPLE. Fatal " << std::endl;
            G4String mStr(mStrStr.str());
            G4Exception("LBNERunAction::BeginOfRunAction", " ", RunMustBeAborted, mStr.c_str()); 
      } 
    
    if(theRunManager -> GetCreateDk2NuOutput()) { 
      if(!(analysis->CreateDk2NuOutput())) {
            std::ostringstream mStrStr;
            mStrStr << "PROBLEM: FAILED TO OPEN Dk2Nu OUTPUT NTUPLE. Fatal " << std::endl;
            G4String mStr(mStrStr.str());
            G4Exception("LBNERunAction::BeginOfRunAction", " ", RunMustBeAborted, mStr.c_str()); 
      }
      analysis->fillDkMeta();
   }
      
   std::cout << std::endl;
   const bool doLBNEQuickPiToNu = false; // default should be false!.  It will slow down execution otherise. 
   if (doLBNEQuickPiToNu) {
      LBNEQuickPiToNuVect::Instance()->YesDoThis();
      
      std::ostringstream fNameStr; 
      fNameStr << theRunManager->GetOutputNtpFileName() << "_";
      fNameStr.setf(std::ios::right); fNameStr.fill('0'); fNameStr.width(3);
      fNameStr << theRunManager->GetCurrentRun()->GetRunID();
      std::string fName(fNameStr.str());
      std::cerr << " QuickPitoNu file name " << fName << std::endl;
      LBNEQuickPiToNuVect::Instance()->open(fName);
   }                                     // if usefull 

   if (fDoComputeEDepInHorns) {
      std::string aNameHH("HornAConceptDesignJan2017");
      fG4MARSDeDxMap = new LBNFDeDxMap(aNameHH);
   }

   std::cout << "LBNERunAction::BeginOfRunAction() - ...completed run initialization. " << std::endl;
   std::cout << std::endl;

  
}

void LBNERunAction::CheckMemoryUsage

(

double

VMLimit

)

const

Definition at line 367 of file LBNERunAction.cc.

                                                         { // in kb (I think..) 

// http://stackoverflow.com/questions/669438/how-to-get-memory-usage-at-run-time-in-c
   using std::ios_base;
   using std::ifstream;
   using std::string;

   double vm_usage     = 0.0;
//   double resident_set = 0.0;

   // 'file' stat seems to give the most reliable results
   //
   ifstream stat_stream("/proc/self/stat",ios_base::in);

   // dummy vars for leading entries in stat that we don't care about
   //
   string pid, comm, state, ppid, pgrp, session, tty_nr;
   string tpgid, flags, minflt, cminflt, majflt, cmajflt;
   string utime, stime, cutime, cstime, priority, nice;
   string O, itrealvalue, starttime;

   // the two fields we want
   //
   unsigned long vsize;
   long rss;

   stat_stream >> pid >> comm >> state >> ppid >> pgrp >> session >> tty_nr
               >> tpgid >> flags >> minflt >> cminflt >> majflt >> cmajflt
               >> utime >> stime >> cutime >> cstime >> priority >> nice
               >> O >> itrealvalue >> starttime >> vsize >> rss; // don't care about the rest

   stat_stream.close();

//   long page_size_kb = sysconf(_SC_PAGE_SIZE) / 1024; // in case x86-64 is configured to use 2MB pages
// Don't care about resident memory.. 
//
   vm_usage     = vsize / 1024.0;
//   std::cerr << " Current memory usage " << vm_usage << std::endl;
   if (vm_usage > VMlimit) {
      std::cerr << " Current memory usage " << vm_usage << " Too much above limit " << VMlimit << std::endl;
      sleep(10);
     exit(2);
   }

}

void LBNERunAction::CheckOKToRun

(

)

Definition at line 335 of file LBNERunAction.cc.

{
   //check the simulation
   // Obsolte..
   /*
   G4bool knownSim = false;
   const vstring_t simvec = fLBNEData -> GetListOfSimulations();
   const G4String simulation = fLBNEData -> GetSimulation();
   
   for(vstring_t::const_iterator sit = simvec.begin(); sit != simvec.end(); ++sit)
   {
      if (simulation == *sit) { knownSim = true; break;}
   }

   if(knownSim) { fLBNEData -> SetOKToRun(true);}
   else
   { 
      std::cout << "LBNERunAction::CheckOKToRun() - PROBLEM: Simulation "
                << "\"" << simulation << "\"" << " unknown. Currently "
                << "implemented simulations are " << std::endl;

      for(vstring_t::const_iterator sit = simvec.begin(); sit != simvec.end(); ++sit)
      {
         std::cout << "     " << *sit << std::endl;
      }
      std::cout << std::endl;

      fLBNEData -> SetOKToRun(false);
   }
   */
   //
}

void LBNERunAction::EndOfRunAction

(

const G4Run *

aRun

)

Definition at line 204 of file LBNERunAction.cc.

{
//   const bool doLBNEQuickPiToNu = true; // default should be false!.  Historical.. 
   if(LBNEQuickPiToNuVect::Instance()->doIt()) LBNEQuickPiToNuVect::Instance()->close();

   LBNERunManager* theRunManager = dynamic_cast<LBNERunManager*>(G4RunManager::GetRunManager());
    if(theRunManager->GetVerboseLevel() > 0)
   {
      std::cout << "LBNERunAction::EndOfRunAction() Called." << std::endl;
   }

   const LBNESteppingAction *pStep = dynamic_cast<const LBNESteppingAction*>(theRunManager->GetUserSteppingAction());
   std::cout << " End of Run.  Number of tracks killed because stuck " 
           <<  pStep->GetNumTracksKilledAsStuck() << std::endl;

   

   G4String spaces = "   ";
   std::cout << std::endl;
   std::cout << "LBNERunAction::EndOfRunAction() - Terminating Run " 
             << aRun->GetRunID() << "..." << std::endl;

   //
   //Random number generator
   //
   G4String randomFile="rndm/endOfRun_";
   char runN[4];
   snprintf(runN, 4, "%03d",aRun->GetRunID());
   randomFile.append(runN);
   randomFile.append(".rndm");
   CLHEP::HepRandom::saveEngineStatus(randomFile);
   std::cout << spaces << "Closing Random Number generator... " << std::endl
             << spaces << "   Seed = " << CLHEP::HepRandom::getTheSeed() << std::endl
             << spaces << "   Random engine status saved in "<< randomFile << std::endl;
   
   //
   //Close Input and Output Ntuple
   //
      
   std::cout << spaces << "Closing Input File... " << std::endl;
   LBNEPrimaryGeneratorAction* primaryGeneratorAction = (theRunManager->GetLBNEPrimaryGeneratorAction());
   
   if(primaryGeneratorAction)
   {
      primaryGeneratorAction->CloseNtuple();
   }
   
   if(theRunManager -> GetCreateOutput()) 
   {
      std::cout << spaces << "Closing Output File... " << std::endl;
      LBNEAnalysis* analysis = LBNEAnalysis::getInstance();
      analysis->CloseOutput();
   }
   if(theRunManager -> GetCreateDk2NuOutput()) 
   {
      std::cout << spaces << "Closing Dk2Nu Output File... " << std::endl;
      LBNEAnalysis* analysis = LBNEAnalysis::getInstance();
      analysis->CloseDk2NuOutput();
   }
   //
   if (fDoComputeEDepInGraphite) { 
     int numEvts = aRun->GetNumberOfEvent();
     fMeanEDepInGraphite /= numEvts;
     double aa = fRMSEDepInGraphite;
     fRMSEDepInGraphite = (std::sqrt(aa - fMeanEDepInGraphite*fMeanEDepInGraphite*numEvts)/(numEvts - 1));
     std::cout << "LBNERunAction::EndOfRunAction .. Average energy deposition in target [GeV]  " 
               << fMeanEDepInGraphite/CLHEP::GeV << " +- " << fRMSEDepInGraphite/CLHEP::GeV << std::endl;
   }
   if (fDoComputeEDepInArgonGas) { 
     int numEvts = aRun->GetNumberOfEvent();
     fMeanEDepInArgonGasHorn1 /= numEvts;
     double aa = fRMSEDepInArgonGasHorn1;
     fRMSEDepInArgonGasHorn1 = (std::sqrt(aa - fMeanEDepInArgonGasHorn1*fMeanEDepInArgonGasHorn1*numEvts)/(numEvts - 1));
     std::cout << "LBNERunAction::EndOfRunAction .. Average energy deposition in Argon Gas Horn1 [MeV]  "
               << fMeanEDepInArgonGasHorn1/CLHEP::MeV <<  " +- " << fRMSEDepInArgonGasHorn1/CLHEP::MeV << std::endl;
     fMeanEDepInArgonGasHorn2 /= numEvts;
     aa = fRMSEDepInArgonGasHorn2;
     fRMSEDepInArgonGasHorn2 = (std::sqrt(aa - fMeanEDepInArgonGasHorn2*fMeanEDepInArgonGasHorn2*numEvts)/(numEvts - 1));
     std::cout << "LBNERunAction::EndOfRunAction .. Average energy deposition in Argon Gas Horn2 [MeV]  " 
              << fMeanEDepInArgonGasHorn2/CLHEP::MeV << " +- " << fRMSEDepInArgonGasHorn2/CLHEP::MeV << std::endl;
   }

   if (fDoComputeEDepInHorns) {
     std::string aDir("./");
     const char* aDirGG = getenv("_CONDOR_SCRATCH_DIR");
     if (aDirGG != 0) aDir = std::string(aDirGG);
     fG4MARSDeDxMap->dumpASCII(aDir, aRun->GetNumberOfEvent());
     std::cerr << " Before deleting fG4MARSDeDxMap " << std::endl;
//          delete(fG4MARSDeDxMap);
//     std::cerr << " After deleting fG4MARSDeDxMap " << std::endl; 
   }

   std::cout << "LBNERunAction::EndOfRunAction() - ...completed run termination. " << std::endl;
   std::cout << std::endl;

/*

   G4cout << G4endl;
  G4cout << G4endl;
  G4cout << "********************************************************************" << G4endl;
  G4cout << "********************************************************************" << G4endl;
  G4cout << "LBNERunAction::EndOfRunAction..." << G4endl;
  G4cout << "********************************************************************" << G4endl;

  G4String randomFile="rndm/endOfRun_";
  char runN[4];
  sprintf(runN,"%04d",aRun->GetRunID());
  randomFile.append(runN);
  randomFile.append(".rndm");
  CLHEP::HepRandom::saveEngineStatus(randomFile);
  G4cout << "  Random engine status at the end of the run saved in "<<randomFile<<G4endl;
  LBNEAnalysis* analysis = LBNEAnalysis::getInstance();
  analysis->finish();

  G4cout << "********************************************************************" << G4endl;
  G4cout << "LBNERunAction::EndOfRunAction - Completed." << G4endl;
  G4cout << "********************************************************************" << G4endl;
  G4cout << "********************************************************************" << G4endl;
  G4cout << G4endl;
  G4cout << G4endl;
*/
}

bool LBNERunAction::GetDebugLevel ( ) const

inline

Definition at line 46 of file LBNERunAction.hh.

{ return fVerboseLevel; }  

bool LBNERunAction::GetDoComputeEDepInArgonGas ( ) const

inline

Definition at line 54 of file LBNERunAction.hh.

{return fDoComputeEDepInArgonGas;}

bool LBNERunAction::GetDoComputeEDepInGraphite ( ) const

inline

Definition at line 48 of file LBNERunAction.hh.

{return fDoComputeEDepInGraphite;}

bool LBNERunAction::GetDoComputeEDepInHorns ( ) const

inline

Definition at line 62 of file LBNERunAction.hh.

{ return fDoComputeEDepInHorns;}

LBNFDeDxMap* LBNERunAction::GetG4MARSDeDxMap ( ) const

inline

Definition at line 63 of file LBNERunAction.hh.

{return fG4MARSDeDxMap;}

void LBNERunAction::SetDebugLevel ( bool  t )

inline

Definition at line 45 of file LBNERunAction.hh.

{ fVerboseLevel = t; }  

void LBNERunAction::SetDoComputeEDepInArgonGas ( bool  t ) const

inline

Definition at line 53 of file LBNERunAction.hh.

{ fDoComputeEDepInArgonGas = t;}

void LBNERunAction::SetDoComputeEDepInGraphite ( bool  t ) const

inline

Definition at line 47 of file LBNERunAction.hh.

{ fDoComputeEDepInGraphite = t;}

void LBNERunAction::SetDoComputeEDepInHorns ( bool  t ) const

inline

Definition at line 61 of file LBNERunAction.hh.

{ fDoComputeEDepInHorns = t;}

void LBNERunAction::storeEDepInArgonGas ( double  edepH1, double  edepH2  ) const

inline

Definition at line 55 of file LBNERunAction.hh.

                                                                        {  /* not really... */
      fMeanEDepInArgonGasHorn1 += edepH1;
      fRMSEDepInArgonGasHorn1 += edepH1*edepH1;
      fMeanEDepInArgonGasHorn2 += edepH2;
      fRMSEDepInArgonGasHorn2 += edepH2*edepH2;
   }

void LBNERunAction::storeEDepInGraphite ( double  edep ) const

inline

Definition at line 49 of file LBNERunAction.hh.

                                                       {  /* not really... */
      fMeanEDepInGraphite += edep;
      fRMSEDepInGraphite += edep*edep;
   }

Member Data Documentation

bool LBNERunAction::fDoComputeEDepInArgonGas

mutableprivate

Definition at line 34 of file LBNERunAction.hh.

bool LBNERunAction::fDoComputeEDepInGraphite

mutableprivate

Definition at line 33 of file LBNERunAction.hh.

bool LBNERunAction::fDoComputeEDepInHorns

mutableprivate

Definition at line 41 of file LBNERunAction.hh.

LBNFDeDxMap* LBNERunAction::fG4MARSDeDxMap

mutableprivate

Definition at line 42 of file LBNERunAction.hh.

double LBNERunAction::fMeanEDepInArgonGasHorn1

mutableprivate

Definition at line 37 of file LBNERunAction.hh.

double LBNERunAction::fMeanEDepInArgonGasHorn2

mutableprivate

Definition at line 39 of file LBNERunAction.hh.

double LBNERunAction::fMeanEDepInGraphite

mutableprivate

Definition at line 35 of file LBNERunAction.hh.

double LBNERunAction::fRMSEDepInArgonGasHorn1

mutableprivate

Definition at line 38 of file LBNERunAction.hh.

double LBNERunAction::fRMSEDepInArgonGasHorn2

mutableprivate

Definition at line 40 of file LBNERunAction.hh.

double LBNERunAction::fRMSEDepInGraphite

mutableprivate

Definition at line 36 of file LBNERunAction.hh.

int LBNERunAction::fVerboseLevel

private

Definition at line 30 of file LBNERunAction.hh.

LBNERunActionMessenger* LBNERunAction::runMessenger

private

Definition at line 31 of file LBNERunAction.hh.

---

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

• g4lbne/include/LBNERunAction.hh
• g4lbne/src/LBNERunAction.cc