Efficiencies.cxx

Go to the documentation of this file.

 
 
#include "Efficiencies.h"
 

IDTPM2::Efficiencies::Efficiencies( const std::string& name, ToolHandle<GenericMonitoringTool>* m ) 
  : IPlotBase( name, m )
{ }

void IDTPM2::Efficiencies::initialize() {
  
  //   MY_MSG_DEBUG <<  "Booking efficiency plots in " << getDirectory()  << endmsg; 
 
  m_eff_vs_inclusive     = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_inclusive" );
  m_eff_vs_ntracks   = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_ntracks" );
  m_eff_vs_pt            = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_pt" );
  // may still want these in the future ...
  //  m_eff_vs_logPt     = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_logPt" );
  //  m_eff_vs_lowPt     = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_lowPt" );
  m_eff_vs_eta           = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_eta" );
  m_eff_vs_phi           = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_phi" );
  m_eff_vs_d0            = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_d0" );
  m_eff_vs_z0            = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_z0" );
 
  m_eff_vs_Et            = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_Et" );
  m_eff_vs_EtovPt    = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_EtovPt" );
 
  
  m_eff_vs_actualMu      = TIDA::Histogram2D<float,bool>( monTool(), "eff_vs_actualMu" );
 
  // if( m_doTruthMuPlots ) m_eff_vs_truthMu = TIDA::Histogram<float,bool>( monTool(), "eff_vs_truthMu" );
  // }

 
 
}
 

void IDTPM2::Efficiencies::execute( const TrackAnalysisStore& collections ) const {  
 
  const std::vector<const TrackView*>& ref  = collections.reference_tracks();

  //  const std::vector<const TrackView*>& test = collections.test_tracks();
 
  double mu = collections.eventInfo()->actualInteractionsPerCrossing();
    
  const TIDA::Associator<const TrackView>* matcher = collections.matcher();
 
  for ( size_t i=0 ; i<ref.size(); i++ ) {
 
    const TrackView& p = *ref[i];
 
    bool matched = false;

    if ( matcher->matched(ref[i]) ) matched = true;
  
    float ppt    = p.pt() / Gaudi::Units::GeV;
    float peta   = p.eta();
    float pphi   = p.phi();
    float pd0    = p.d0();
    float pz0    = p.z0();
 
    float Et     = p.ET() / Gaudi::Units::GeV;
 
    if ( p.hasParent() ) Et = p.parent()->ET()*0.001;
    
    float EtovPt = Et/ppt; // pt should never be 0 so don't bother to check that
    
    m_eff_vs_inclusive.Fill( 1,    matched );
 
    m_eff_vs_ntracks.Fill( ref.size(),  matched );
 
    m_eff_vs_pt.Fill(        ppt,  matched );
    //    m_eff_vs_logPt.Fill(     ppt,  matched );
    //    m_eff_vs_lowPt.Fill(     ppt,  matched );
    m_eff_vs_eta.Fill(       peta, matched );
    m_eff_vs_phi.Fill(       pphi, matched );
    m_eff_vs_d0.Fill(        pd0,  matched );
    m_eff_vs_z0.Fill(        pz0,  matched );
 
    m_eff_vs_Et.Fill(         Et,  matched );
    m_eff_vs_EtovPt.Fill( EtovPt,  matched );
 
    m_eff_vs_actualMu.Fill(   mu,  matched );

    //  m_eff_vs_eta_vs_pt.Fill(    ppt, matched, ppt );
    //  m_eff_vs_eta_vs_phi.Fill(   pphi,  pet, matched );
    //  m_eff_vs_z0_vs_d0.Fill(     pd0,   pz0, matched );
    //  float ptheta = p.theta(); /// no point setting it unless it is needed
    //  m_eff_vs_z0sin_vs_d0.Fill(  pz0*std::sin(ptheta), pd0, matched );
    
    //    if( m_doTruthMuPlots ) m_eff_vs_truthMu.Fill( matched, matched );
  }
}
 
 
                                      
void IDTPM2::Efficiencies::finalize() {  MY_MSG_DEBUG <<  "Finalising efficiency plots"  << endmsg; }