MCMonitorCosmic.cxx

Go to the documentation of this file.

//
// MCMonitorCosmic
//    Modified version of Kregg Arm's MCMonitor class to generate
//    histograms specific to comparison of ptsim to gminos cosmic ray
//    muons
#include "TClonesArray.h"

#include "MCMonitorCosmic.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"
#include "MCNtuple/NtpMCDigiScintHit.h"

JOBMODULE(MCMonitorCosmic, "MCMonitorCosmic",
          "Produce MC plots for cosmic ray muon validation.");
CVSID("$Id: MCMonitorCosmic.cxx,v 1.1 2007/10/10 03:11:22 schubert Exp $");
ClassImp(MCMonitorCosmic)

//___________________________________________________________________________
MCMonitorCosmic::MCMonitorCosmic() : 
  fHists(new MCMonitorCosmicHistograms("MCMonCosmic")),fOutRootFileName("") {
  // Default constructor

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

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

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

}
//___________________________________________________________________________
void MCMonitorCosmic::EndJob() {
  //
  //
  MSG("MCMon",Msg::kDebug) << "In MCMonitorCosmic::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 MCMonitorCosmic::BeginFile() {
  //

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

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

  MSG("MCMon",Msg::kDebug)
    << "In MCMonitorCosmic::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 MCMonitorCosmic::Ana(const MomNavigator* mom)
{

  MSG("MCMon",Msg::kDebug) << "In MCMonitorCosmic::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;
  NtpSRRecord* srrec = 0;
  NtpTHRecord* threc = 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;
    }
    srrec = dynamic_cast<NtpSRRecord*>(mom->GetFragment("NtpSRRecord"));
    threc = dynamic_cast<NtpTHRecord*>(mom->GetFragment("NtpTHRecord"));
  } 

  fHists -> FillHistograms(strec, mcrec, srrec, threc);

  // Increment #snarls/file
  fNumSnarlsInCurrFile++;

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

  return JobCResult::kPassed; 

}
//___________________________________________________________________________
const Registry& MCMonitorCosmic::DefaultConfig() const
{
  //
  // Supply the default configuration for the module
  //
  MSG("MCMon",Msg::kDebug) << "In MCMonitorCosmic::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 MCMonitorCosmic::Config(const Registry& r)
{
  //
  // Configure the module given the Registry r
  //

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

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

---