IDTPM2::Efficiencies Class Reference

#include <Efficiencies.h>

Inheritance diagram for IDTPM2::Efficiencies:

Collaboration diagram for IDTPM2::Efficiencies:

Public Member Functions
	Efficiencies (const std::string &name, ToolHandle< GenericMonitoringTool > *m=0)
	Constructor.
virtual	~Efficiencies ()=default
	Destructor.
void	initialize ()
	Dedicated fill method (for tracks and/or truth particles).
void	execute (const TrackAnalysisStore &collections) const
void	finalize ()
	Print out final stats on histograms.
const std::string &	name () const
	generic utilities, just in case
ToolHandle< GenericMonitoringTool > *	monTool ()
	mon tool
const ToolHandle< GenericMonitoringTool > *	monTool () const
void	set_monTool (ToolHandle< GenericMonitoringTool > *m)
	do we need / should we have this specific setter ??
std::string	realPath (const std::string &tag="")
	This is a horrific function, has to get ALL the histograms, just to be able to work out how to get the actual path to the histograms for this tool - there simply MUST be a better way ...
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

Static Public Member Functions
template<typename Derived>
static std::unique_ptr< IPlotBase >	create (const std::string &name, ToolHandle< GenericMonitoringTool > *m)
	factory creation mathod

Protected Member Functions
void	copyHist (TH1 *hsrc, TH1 *hdest, bool copyoverflow=true) const
	resolution porocessing stuff these don't really even need to be class functions, they could be simply utility functions
void	getHist (const std::string &hname, TH1 *hsrc)
void	getHist (ServiceHandle< ITHistSvc > &histsvc, const std::string &hname, TH1 *hsrc) const
void	processResolution (const std::string &path, bool writeflag=false)
void	processResolution (ServiceHandle< ITHistSvc > &histsvc, const std::string &path, bool writeflag=false) const
virtual ServiceHandle< ITHistSvc > &	histSvc ()
	we may well need these at some point ...

Protected Attributes
const std::string &	m_name
ToolHandle< GenericMonitoringTool > *	m_monTool

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
TIDA::Histogram2D< float, bool >	m_eff_vs_inclusive
	NB: the egregious root convention for a TEfficiency is Fill( ismatched, variable ), but that is awful, and inconsistent, so we swap that round.
TIDA::Histogram2D< float, bool >	m_eff_vs_ntracks
TIDA::Histogram2D< float, bool >	m_eff_vs_pt
TIDA::Histogram2D< float, bool >	m_eff_vs_eta
TIDA::Histogram2D< float, bool >	m_eff_vs_phi
TIDA::Histogram2D< float, bool >	m_eff_vs_d0
TIDA::Histogram2D< float, bool >	m_eff_vs_z0
TIDA::Histogram2D< float, bool >	m_eff_vs_Et
TIDA::Histogram2D< float, bool >	m_eff_vs_EtovPt
TIDA::Histogram2D< float, bool >	m_eff_vs_actualMu
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels).
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging).

Detailed Description

Definition at line 18 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

Constructor & Destructor Documentation

Efficiencies()

IDTPM2::Efficiencies::Efficiencies	(	const std::string &	name,
		ToolHandle< GenericMonitoringTool > *	m = 0 )

Constructor.

---

--— Constructor --—

Definition at line 16 of file Efficiencies.cxx.

  : IPlotBase( name, m )
{ }

~Efficiencies()

virtual IDTPM2::Efficiencies::~Efficiencies ( )

virtualdefault

Destructor.

Member Function Documentation

copyHist()

void IDTPM2::IPlotBase::copyHist ( TH1 * hsrc, TH1 * hdest, bool copyoverflow = true ) const

protectedinherited

resolution porocessing stuff these don't really even need to be class functions, they could be simply utility functions

Definition at line 20 of file IPlotBase.cxx.

                                                                               {
  
  if ( hsrc->GetNbinsX()!=hdest->GetNbinsX() ) return; // do we want to give an error message ?? throw and exception ??
 
  int imin = 0;
  int imax = hsrc->GetNbinsX()+1; 
 
  if ( !copyoverflow ) {
    imin = 1;
    imax = hsrc->GetNbinsX(); 
  }
  
  for  ( int i=imin ; i<=imax ; i++ ) {
    hdest->SetBinContent( i, hsrc->GetBinContent(i) );
    hdest->SetBinError( i, hsrc->GetBinError(i) );
  }
 
  hdest->SetEntries( hsrc->GetEntries() );
}

create()

template<typename Derived>

std::unique_ptr< IPlotBase > IDTPM2::IPlotBase::create ( const std::string & name, ToolHandle< GenericMonitoringTool > * m )

inlinestaticinherited

factory creation mathod

Definition at line 50 of file IPlotBase.h.

                                                               {
    return std::make_unique<Derived>( name, m );
  } 

execute()

void IDTPM2::Efficiencies::execute ( const TrackAnalysisStore & collections ) const

virtual

---

— Dedicated fill method —

don't need the actual test tracks

is there a matched track ??

Compute track parameters

Fill the histograms do we really need the weighted histograms ???

do we really even want these ???

truth mu ??

Implements IDTPM2::IPlotBase.

Definition at line 61 of file Efficiencies.cxx.

                                                                              {  
 
  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 );
  }
}

finalize()

void IDTPM2::Efficiencies::finalize ( )

virtual

Print out final stats on histograms.

---

--— finalizePlots --—

Implements IDTPM2::IPlotBase.

Definition at line 130 of file Efficiencies.cxx.

{  MY_MSG_DEBUG <<  "Finalising efficiency plots"  << endmsg; }

getHist() [1/2]

void IDTPM2::IPlotBase::getHist ( const std::string & hname, TH1 * hsrc )

inlineprotectedinherited

Definition at line 88 of file IPlotBase.h.

{ return getHist( histSvc(), hname, hsrc ); }

getHist() [2/2]

void IDTPM2::IPlotBase::getHist ( ServiceHandle< ITHistSvc > & histsvc, const std::string & hname, TH1 * hsrc ) const

protectedinherited

Definition at line 42 of file IPlotBase.cxx.

                                                                                                          {
  TH1* h = 0;
  if ( histsvc->getHist( hname, h ).isFailure() ) MY_MSG_WARNING << "Could not retieve histogram: " << hname << endmsg;
  else copyHist( hsrc, h );
}

histSvc()

virtual ServiceHandle< ITHistSvc > & IDTPM2::IPlotBase::histSvc ( )

inlineprotectedvirtualinherited

we may well need these at some point ...

this method is dangerous - should not be used anywhere except finalise or stop

Definition at line 105 of file IPlotBase.h.

                                              {
    if ( monTool()==0 ) throw std::runtime_error("IPlotBase: monTool not set");
    return (*monTool())->histogramService();
  }

initialize()

void IDTPM2::Efficiencies::initialize ( )

virtual

Dedicated fill method (for tracks and/or truth particles).

---

— Book the histograms —

why would we ever not want the efficiency vs mu ?

for later ...

m_eff_vs_eta_vs_pt = TIDA::Histogram<float,bool>( monTool(), "eff_vs_eta_vs_pt" ); m_eff_vs_eta_vs_phi = TIDA::Histogram<float,bool>( monTool(), "eff_vs_eta_vs_phi" ); m_eff_vs_z0_vs_d0 = TIDA::Histogram<float,bool>( monTool(), "eff_vs_z0_vs_d0" ); m_eff_vs_z0sin_vs_d0 = TIDA::Histogram<float,bool>( monTool(), "eff_vs_z0sin_vs_d0" );

Implements IDTPM2::IPlotBase.

Definition at line 23 of file Efficiencies.cxx.

                                    {
  
  //   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" );
  // }

 
 
}

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

monTool() [1/2]

ToolHandle< GenericMonitoringTool > * IDTPM2::IPlotBase::monTool ( )

inlineinherited

mon tool

Definition at line 71 of file IPlotBase.h.

{ return m_monTool; }

monTool() [2/2]

const ToolHandle< GenericMonitoringTool > * IDTPM2::IPlotBase::monTool ( ) const

inlineinherited

Definition at line 73 of file IPlotBase.h.

{ return m_monTool; }

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 167 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 182 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  // If user did not set explicit message level we have to initialize
  // the messaging and retrieve the default via the MessageSvc.
  if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
 
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

name()

const std::string & IDTPM2::IPlotBase::name ( ) const

inlineinherited

generic utilities, just in case

Definition at line 68 of file IPlotBase.h.

{ return m_name; }

processResolution() [1/2]

void IDTPM2::IPlotBase::processResolution ( const std::string & path, bool writeflag = false )

inlineprotectedinherited

Definition at line 89 of file IPlotBase.h.

{ return processResolution( histSvc(), path, writeflag ); }

processResolution() [2/2]

void IDTPM2::IPlotBase::processResolution ( ServiceHandle< ITHistSvc > & histsvc, const std::string & path, bool writeflag = false ) const

protectedinherited

process the histogramss NB: might be better if instead of adding the 2D histogram and then copying the 1D hjistograms back, we passed in the 1D histogram pointers and filled them directly in the Resplot

use the provides 1D, maan, sigme histograms etc would like to do this ...

use the provides 1D, maan, sigme histograms etc

Definition at line 49 of file IPlotBase.cxx.

                                                                                                                        {
 
  TH2* hist2d = 0;
  
  if ( histsvc->getHist( path+"2D", hist2d ).isFailure() ) MY_MSG_WARNING << "Could not retieve 2D: " << (path+"2D") << endmsg;
  else {
    
    
    Resplot::Resplot rp( path, hist2d );
 
    rp.Finalise( Resplot::FitNull95 );
 
    if ( !writeflag ) {
      //  if ( histsvc->regHist( path+"1D",    rp.H1D()   ).isFailure() ) MY_MSG_WARNING << "Unable to register histogram: 1D:    " << name() << endmsg; 
      //  if ( histsvc->regHist( path+"mean",  rp.Mean()  ).isFailure() ) MY_MSG_WARNING << "Unable to register histogram: mean:  " << name() << endmsg; 
      //  if ( histsvc->regHist( path+"sigma", rp.Sigma() ).isFailure() ) MY_MSG_WARNING << "Unable to register histogram: sigma: " << name() << endmsg; 
      //  }
      //  else {
      getHist( histsvc, path+"1D",    rp.H1D()   );
      getHist( histsvc, path+"mean",  rp.Mean()  );
      getHist( histsvc, path+"sigma", rp.Sigma() );
    }
  }
 
}

realPath()

std::string IDTPM2::IPlotBase::realPath ( const std::string & tag = "" )

inherited

This is a horrific function, has to get ALL the histograms, just to be able to work out how to get the actual path to the histograms for this tool - there simply MUST be a better way ...

can add on our "fixed" histogram "Vertex" domain

Definition at line 86 of file IPlotBase.cxx.

                                                          {
 
  std::string dir  = (*monTool())->getPath();

  std::string stub = dir.substr( dir.find('/', 1) + 1 ) + tag; 
  
  const std::vector<std::string>& allhists = histSvc()->getHists();
 
  std::string basestr = "";
  
  for ( size_t i=0 ; i<allhists.size() ; i++ ) {
    if ( allhists[i].find(stub)!=std::string::npos ) {
      std::string s = allhists[i];
      basestr = s.substr(0,s.find("/"+stub)+1);
      break;
    }
  }
 
  return basestr + stub;
  
}

set_monTool()

void IDTPM2::IPlotBase::set_monTool ( ToolHandle< GenericMonitoringTool > * m )

inlineinherited

do we need / should we have this specific setter ??

Definition at line 76 of file IPlotBase.h.

{ m_monTool=m; }

setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging).

Definition at line 141 of file AthMessaging.h.

m_eff_vs_actualMu

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_actualMu

private

Definition at line 58 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_d0

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_d0

private

Definition at line 52 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_Et

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_Et

private

Definition at line 55 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_eta

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_eta

private

Definition at line 50 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_EtovPt

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_EtovPt

private

Definition at line 56 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_inclusive

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_inclusive

private

NB: the egregious root convention for a TEfficiency is Fill( ismatched, variable ), but that is awful, and inconsistent, so we swap that round.

Definition at line 43 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_ntracks

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_ntracks

private

Definition at line 45 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_phi

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_phi

private

Definition at line 51 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_pt

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_pt

private

Definition at line 47 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_eff_vs_z0

TIDA::Histogram2D<float,bool> IDTPM2::Efficiencies::m_eff_vs_z0

private

Definition at line 53 of file TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

m_monTool

ToolHandle<GenericMonitoringTool>* IDTPM2::IPlotBase::m_monTool

protectedinherited

Definition at line 116 of file IPlotBase.h.

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels).

Definition at line 132 of file AthMessaging.h.

m_name

const std::string& IDTPM2::IPlotBase::m_name

protectedinherited

Definition at line 114 of file IPlotBase.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

---

The documentation for this class was generated from the following files:

• TrigMonitoring/TrigIDR4Monitoring/src/Efficiencies.h
• Efficiencies.cxx