MCAnaProcessor.cxx

Go to the documentation of this file.

 
#include "MCAnaProcessor.h"
#include <iostream>
 
MCAnaProcessor::MCAnaProcessor(const std::string& name, Process& process) : Processor(name,process){}
//TODO CHECK THIS DESTRUCTOR
MCAnaProcessor::~MCAnaProcessor(){}
 
 
void MCAnaProcessor::configure(const ParameterSet& parameters) {
    std::cout << "Configuring MCAnaProcessor" << std::endl;
    try
    {
        debug_           = parameters.getInteger("debug");
        anaName_         = parameters.getString("anaName");
        partColl_        = parameters.getString("partColl");
        trkrHitColl_     = parameters.getString("trkrHitColl");
        ecalHitColl_     = parameters.getString("ecalHitColl");
        histCfgFilename_ = parameters.getString("histCfg");
        analysis_        = parameters.getString("analysis");
    }
    catch (std::runtime_error& error)
    {
        std::cout << error.what() << std::endl;
    }
}
void MCAnaProcessor::configure(const ParameterSet& parameters) {…}
 
void MCAnaProcessor::initialize(TTree* tree) {
    tree_= tree;
    // init histos
    histos = new MCAnaHistos(anaName_);
    histos->loadHistoConfig(histCfgFilename_);
    histos->DefineHistos();
    histos->Define2DHistos();
 
    // init TTree
    tree_->SetBranchAddress(partColl_.c_str(), &mcParts_, &bmcParts_);
    if (tree_->FindBranch(trkrHitColl_.c_str()))
        tree_->SetBranchAddress(trkrHitColl_.c_str(), &mcTrkrHits_, &bmcTrkrHits_);
    else
        std::cout<<"WARNING: No tracker hit collection, will skip FillMCTrackerHits! "<<std::endl;
 
    if ( tree_->FindBranch(ecalHitColl_.c_str()))
        tree_->SetBranchAddress(ecalHitColl_.c_str(), &mcEcalHits_, &bmcEcalHits_);
    else
        std::cout<<"WARNING: No Ecal hit collection, will skip FillMCEcalHits"<<std::endl;
 
}
void MCAnaProcessor::initialize(TTree* tree) {…}
 
bool MCAnaProcessor::process(IEvent* ievent) {
 
  histos->FillMCParticles(mcParts_, analysis_);
    if(mcTrkrHits_)
        histos->FillMCTrackerHits(mcTrkrHits_);
    if(mcEcalHits_)
        histos->FillMCEcalHits(mcEcalHits_);
 
    return true;
}
bool MCAnaProcessor::process(IEvent* ievent) {…}
 
void MCAnaProcessor::finalize() {
 
    histos->saveHistos(outF_, anaName_);
    delete histos;
    histos = nullptr;
}
void MCAnaProcessor::finalize() {…}
 
DECLARE_PROCESSOR(MCAnaProcessor);