MCMonitorBeam.cxx

Go to the documentation of this file.

//
// MCMonitorBeam
//    Modified version of Kregg Arm's MCMonitor class to generate
//    histograms specific to comparison of ptsim to gminos beam 
//    events
#include "TClonesArray.h"

#include "MCMonitorBeam.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro
#include "HistMan/HistMan.h"
#include "MCNtuple/NtpMCStdHep.h"
#include "MCNtuple/NtpMCTruth.h"

JOBMODULE(MCMonitorBeam, "MCMonitorBeam",
          "Produce MC plots for beam ray validation.");
CVSID("$Id: MCMonitorBeam.cxx,v 1.1 2008/02/05 20:30:05 schubert Exp $");
ClassImp(MCMonitorBeam)

//___________________________________________________________________________
MCMonitorBeam::MCMonitorBeam() : 
  fHists(new MCMonitorBeamHistograms("MCMonBeam")),fOutRootFileName("") {
  // Default constructor

  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam constructor" << endl;

}
//___________________________________________________________________________
MCMonitorBeam::~MCMonitorBeam()
{
  // Destructor
  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam destructor" << endl;

}
//___________________________________________________________________________
void MCMonitorBeam::BeginJob() {
  
  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam::BeginJob" << endl;

}
//___________________________________________________________________________
void MCMonitorBeam::EndJob() {
  //
  //
  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam::EndJob" << endl;

  if (fNumSnarlsInCurrFile !=  0) EndFile();

  unsigned int numFiles = fNumSnarlsInFile.size();
  unsigned int numSnarls = 0;
  for ( unsigned int ifile = 0; ifile < numFiles; ifile++ ) {
    numSnarls += fNumSnarlsInFile[ifile];
  }
  
  fHists->SetCanRebin();
  fHists->AllSumw2();
  
  // Send all of the histograms to output root file
  MSG("MCMon",Msg::kInfo)
     << "MCMonitor::EndJob storing histograms to file "
     << fOutRootFileName.c_str() << endl;
  HistMan hm((fHists->GetName()).c_str());
  hm.WriteOut(fOutRootFileName.c_str());

  // Print a summary of the module
  cout << std::string(75,'-') << endl;
  cout << "Processed " << numFiles << " files" << endl
       << "Processed " << numSnarls << " snarls" << endl
       << "Produced  " << fHists->NumHistograms()
       << " histograms written to file: " << fOutRootFileName.c_str()
       << endl;
  cout << std::string(75,'-') << endl;

}
//___________________________________________________________________________
void MCMonitorBeam::BeginFile() {
  //

  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam::BeginFile" << endl;

}
//___________________________________________________________________________
void MCMonitorBeam::EndFile() {

  MSG("MCMon",Msg::kDebug)
    << "In MCMonitorBeam::EndFile \n" 
    << "Num Snarls in Current File           : "
    << fNumSnarlsInCurrFile << endl;
  // Send 'number of X per file' info to vectors; increment fFileNum
  fNumSnarlsInFile.push_back(fNumSnarlsInCurrFile);

  fNumSnarlsInCurrFile = 0;

}

//___________________________________________________________________________
JobCResult MCMonitorBeam::Ana(const MomNavigator* mom)
{

  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam::Ana" << endl;

  // get NtpSt from the file
  NtpStRecord* strec =
    dynamic_cast<NtpStRecord*>( mom->GetFragment("NtpStRecord") );
  // (just in case strec isn't present . . .)
  NtpMCRecord* mcrec = 0;
  
  if (!strec) {
    mcrec = dynamic_cast<NtpMCRecord*>( mom->GetFragment("NtpMCRecord") );
    if (!mcrec) {
      MSG("MCMon",Msg::kWarning)
        << "No NtpStRecord or NtpMCRecord in Mom! Aborting . . ."
        << endl;
      return JobCResult::kFailed;
    }
  } 

  fHists -> FillHistograms(strec, mcrec);

  // Increment #snarls/file
  fNumSnarlsInCurrFile++;

  MSG("MCMon",Msg::kDebug) << "\n-----> Processed " << fNumSnarlsInCurrFile 
                           << " snarls in this file." << endl;

  return JobCResult::kPassed; 

}
//___________________________________________________________________________
const Registry& MCMonitorBeam::DefaultConfig() const
{
  //
  // Supply the default configuration for the module
  //
  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam::DefaultConfig" << endl;

  static Registry r; // Default configuration for module

  // Set name of config
  std::string name = this->GetName();
  name += ".config.default";
  r.SetName(name.c_str());

  // Set values in configuration
  r.UnLockValues();
  r.Set("OutFile","mcmonitor.root");    // Output root file for histograms
  r.LockValues();

  return r;
}
//___________________________________________________________________________
void MCMonitorBeam::Config(const Registry& r)
{
  //
  // Configure the module given the Registry r
  //

  MSG("MCMon",Msg::kDebug) << "In MCMonitorBeam::Config" << endl;

  fOutRootFileName   = r.GetCharString("OutFile");
}
//___________________________________________________________________________

---