ClusterOnTrackAnaProcessor.cxx

Go to the documentation of this file.

#include "ClusterOnTrackAnaProcessor.h"
#include "TBranch.h"
 
ClusterOnTrackAnaProcessor::ClusterOnTrackAnaProcessor(const std::string& name, Process& process) : Processor(name,process){}
//TODO CHECK THIS DESTRUCTOR
ClusterOnTrackAnaProcessor::~ClusterOnTrackAnaProcessor(){}
 
 
void ClusterOnTrackAnaProcessor::configure(const ParameterSet& parameters) {
 
    std::cout << "Configuring ClusterOnTrackAnaProcessor" << std::endl;
    debug_        = parameters.getInteger("debug");
    anaName_      = parameters.getString("anaName");
    trkColl_      = parameters.getString("trkColl");
    baselineFits_ = parameters.getString("BaselineFits");
    baselineRun_  = parameters.getString("BaselineRun");
    if(debug_ > 0) std::cout << "Configured: " << baselineFits_ << " " << baselineRun_ << std::endl;
 
}
void ClusterOnTrackAnaProcessor::configure(const ParameterSet& parameters) {…}
 
void ClusterOnTrackAnaProcessor::initialize(TTree* tree) {
    tree_= tree;
    clusterHistos = new ClusterHistos(anaName_);
 
    if(debug_ > 0) std::cout << "Loading: " << baselineFits_ << " " << baselineRun_ << std::endl;
    if (!baselineFits_.empty() && !baselineRun_.empty()) {
        /*
        clusterHistos->setBaselineFitsDir(baselineFits_);
        if (!clusterHistos->LoadBaselineHistos(baselineRun_))
            std::cout<<"WARNING: baselines not loaded in Cluster on Track histos."<<std::endl;
            */
        clusterHistos->LoadOfflineBaselines(baselineFits_);
    }
    if(debug_ > 0) std::cout << "Loaded: " << baselineFits_ << " " << baselineRun_ << std::endl;
 
    clusterHistos->Define1DHistos();
    clusterHistos->Define2DHistos();
    if(debug_ > 0) std::cout << "Histos Defined" << std::endl;
 
    //TODO Change this.
    tree_->SetBranchAddress(trkColl_.c_str(),&tracks_,&btracks_);
    if(debug_ > 0) std::cout << "Branch changed to " << trkColl_ << std::endl;
 
}
void ClusterOnTrackAnaProcessor::initialize(TTree* tree) {…}
 
bool ClusterOnTrackAnaProcessor::process(IEvent* ievent) {
    
 
    for (int itrack = 0; itrack<tracks_->size();itrack++) {
        Track *track = tracks_->at(itrack);
        //Loop on hits
        if (track->getSvtHits().GetEntries()==0) {
            std::cout<<"WARNING::track doesn't have hits associated to it"<<std::endl;
            return false;
        }
 
        for (int ihit = 0; ihit<track->getSvtHits().GetEntries(); ++ihit) {
            TrackerHit* hit3d = (TrackerHit*) track->getSvtHits().At(ihit);
            clusterHistos->FillHistograms(hit3d, 1.);
        }
    }
    return true;
}
bool ClusterOnTrackAnaProcessor::process(IEvent* ievent) {…}
 
void ClusterOnTrackAnaProcessor::finalize() {
 
    clusterHistos->saveHistos(outF_, anaName_);
    delete clusterHistos;
    clusterHistos = nullptr;
}
void ClusterOnTrackAnaProcessor::finalize() {…}
 
DECLARE_PROCESSOR(ClusterOnTrackAnaProcessor);