ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
EFTrackingSmearingAlg Class Reference

#include <EFTrackingSmearingAlg.h>

Inheritance diagram for EFTrackingSmearingAlg:
Collaboration diagram for EFTrackingSmearingAlg:

Public Member Functions

 EFTrackingSmearingAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~EFTrackingSmearingAlg ()
virtual StatusCode initialize () override
virtual StatusCode execute () override
virtual StatusCode sysInitialize ()
 Initialization method invoked by the framework.
const ServiceHandle< ITHistSvc > & histSvc () const
 The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.
StatusCode configAthHistogramming (const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
 To be called by the derived classes to fill the internal configuration.
TH1 * bookGetPointer (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TH1 * bookGetPointer (TH1 &histRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
TTree * bookGetPointer (const TTree &treeRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
TGraph * bookGetPointer (const TGraph &graphRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
TEfficiency * bookGetPointer (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency * bookGetPointer (TEfficiency &effRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (TH1 &histRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms.
StatusCode book (const TTree &treeRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TTrees.
StatusCode book (const TGraph &graphRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs.
StatusCode book (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode book (TEfficiency &effRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency.
TH1 * hist (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered histograms of any type.
TH2 * hist2d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 2-d histograms.
TH3 * hist3d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 3-d histograms.
TTree * tree (const std::string &treeName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TTrees.
TGraph * graph (const std::string &graphName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TGraphs.
TEfficiency * efficiency (const std::string &effName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TEfficiency.

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
typedef uint32_t hash_t
 typedef for the internal hash
typedef std::map< const hash_t, TH1 * > HistMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TEfficiency * > EffMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TTree * > TreeMap_t
 Typedef for convenience.
typedef std::map< const hash_t, TGraph * > GraphMap_t
 Typedef for convenience.

Private Member Functions

StatusCode book_histograms ()
StatusCode smearTruthParticles (const EventContext &ctx)
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>
void buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
 Method to build individual booking string.
void myReplace (std::string &str, const std::string &oldStr, const std::string &newStr)
 Helper method to replace sub-string.
hash_t hash (const std::string &histName) const
 Method to calculate a 32-bit hash from a string.

Private Attributes

SG::ReadHandleKey< xAOD::TrackParticleContainer > m_inputTrackParticleKey
SG::WriteHandleKey< xAOD::TrackParticleContainer > m_outputTrackParticleKey
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_inputTruthParticleKey
SG::WriteHandleKey< xAOD::TruthParticleContainer > m_outputTruthParticleKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_d0DecoratorKey {this, "d0", "TruthParticles.d0", "Particle d0 decoration, set at initialisation"}
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_z0DecoratorKey {this, "z0", "TruthParticles.z0", "Particle z0 decoration, set at initialisation"}
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_ptDecoratorKey {this, "pt", "TruthParticles.pt", "Particle pt decoration, set at initialisation"}
DoubleProperty m_SigmaScaleFactor {this, "SmearingScaleFactor", 1, "Set the smearing SF value for the sigmas"}
BooleanProperty m_UseResolutionPtCutOff {this, "UseResolutionPtCutOff", false, "Apply ptCutoff on smearing"}
DoubleProperty m_SetResolutionPtCutOff {this, "SetResolutionPtCutOff", 0., "Set ptCutoff off for semaring"}
DoubleProperty m_inputTracksPtCut {this, "InputTracksPtCutGeV", 0.0, "Set input track pT cut"}
DoubleProperty m_outputTracksPtCut {this, "OutputTracksPtCutGeV", 1.0, "Set the output track pT cut"}
DoubleProperty m_smearedTrackEfficiency {this, "SmearedTrackEfficiency", 1.0, "Set track efficiency for smearing"}
BooleanProperty m_parameterizedTrackEfficiency {this, "ParameterizedTrackEfficiency", false, "Enable parameterized efficiency for smearing"}
BooleanProperty m_enableMonitoring {this, "EnableMonitoring", false, "Enable debugging monitoring of the algorithm"}
BooleanProperty m_smearTruthParticle {this, "SmearTruthParticle", false, "Enable smearing on truth particles, disabling the track smearing"}
BooleanProperty m_parameterizedTrackEfficiency_LRT {this, "ParameterizedTrackEfficiency_LRT", false, "Enable parameterized efficiency for LRT smearing"}
DoubleProperty m_smearedTrackEfficiency_d0high_LRT {this, "SmearedTrackEfficiency_d0high_LRT", 1.0, "Set high d0 value for smeared track efficiency LRT"}
DoubleProperty m_smearedTrackEfficiency_d0low_LRT {this, "SmearedTrackEfficiency_d0low_LRT", 1.0, "Set low d0 value for smeared track efficiency LRT"}
BooleanProperty m_EnableFakes {this,"IncludeDuplicatesAndFakes",false,"Enable fake track production"}
BooleanProperty m_FakeKillerEnable {this,"FakeKillerEnable",false, "disable inclusion of broad fakes"}
BooleanProperty m_IncludeFakesInResolutionCalculation
BooleanProperty m_UseCoinToss
LongLongProperty m_RandomSeed {this,"RandomSeed",0}
std::unique_ptr< FakeTrackSmearer > m_mySmearer
ServiceHandle< ITHistSvc > m_histSvc
 Default constructor: AthHistogramAlgorithm();.
std::string m_prefix
 Name of the ROOT output stream (file).
std::string m_rootDir
 Name of the ROOT directory.
std::string m_histNamePrefix
 The prefix for the histogram THx name.
std::string m_histNamePostfix
 The postfix for the histogram THx name.
std::string m_histTitlePrefix
 The prefix for the histogram THx title.
std::string m_histTitlePostfix
 The postfix for the histogram THx title.
DataObjIDColl m_extendedExtraObjects
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default).
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared
HistMap_t m_histMap
 The map of histogram names to their pointers.
EffMap_t m_effMap
 The map of histogram names to their pointers.
TreeMap_t m_treeMap
 The map of TTree names to their pointers.
GraphMap_t m_graphMap
 The map of TGraph names to their pointers.
std::string m_streamName
 Name of the ROOT output stream (file).
std::string m_name
 Instance name.
MsgStream m_msg
 Cached Message Stream.

Detailed Description

Definition at line 19 of file EFTrackingSmearingAlg.h.

Member Typedef Documentation

◆ EffMap_t

typedef std::map< const hash_t, TEfficiency* > AthHistogramming::EffMap_t
privateinherited

Typedef for convenience.

Definition at line 207 of file AthHistogramming.h.

◆ GraphMap_t

typedef std::map< const hash_t, TGraph* > AthHistogramming::GraphMap_t
privateinherited

Typedef for convenience.

Definition at line 221 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 171 of file AthHistogramming.h.

◆ HistMap_t

typedef std::map< const hash_t, TH1* > AthHistogramming::HistMap_t
privateinherited

Typedef for convenience.

Definition at line 200 of file AthHistogramming.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ TreeMap_t

typedef std::map< const hash_t, TTree* > AthHistogramming::TreeMap_t
privateinherited

Typedef for convenience.

Definition at line 214 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ EFTrackingSmearingAlg()

EFTrackingSmearingAlg::EFTrackingSmearingAlg ( const std::string & name,
ISvcLocator * pSvcLocator )

Definition at line 94 of file EFTrackingSmearingAlg.cxx.

95: AthHistogramAlgorithm( name, pSvcLocator ){}
AthHistogramAlgorithm::AthHistogramAlgorithm
AthHistogramAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor with parameters:
Definition AthHistogramAlgorithm.cxx:31

◆ ~EFTrackingSmearingAlg()

EFTrackingSmearingAlg::~EFTrackingSmearingAlg ( )
virtual

Definition at line 98 of file EFTrackingSmearingAlg.cxx.

98{}

Member Function Documentation

◆ book() [1/8]

StatusCode AthHistogramming::book ( const TEfficiency & eff,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 337 of file AthHistogramming.h.

338{
339 // We need to create a non-const clone
340 TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
341 if ( !effClone ) {
342 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
343 return StatusCode::FAILURE;
344 }
345 return this->book( *effClone, tDir, stream );
346}
endmsg
#define endmsg
Definition AnalysisConfig_Ntuple.cxx:63
AthHistogramming::book
StatusCode book(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.
Definition AthHistogramming.h:305
AthHistogramming::m_msg
MsgStream m_msg
Cached Message Stream.
Definition AthHistogramming.h:250
dqt_zlumi_alleff_HIST.eff
int eff
Definition dqt_zlumi_alleff_HIST.py:113

◆ book() [2/8]

StatusCode AthHistogramming::book ( const TGraph & graphRef,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

◆ book() [3/8]

StatusCode AthHistogramming::book ( const TH1 & hist,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 305 of file AthHistogramming.h.

306{
307 // We need to create a non-const clone
308 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
309 if ( !histClone ) {
310 m_msg << MSG::ERROR << "Couldn't create a TH1 clone" << endmsg;
311 return StatusCode::FAILURE;
312 }
313 return this->book( *histClone, tDir, stream );
314}
AthHistogramming::hist
TH1 * hist(const std::string &histName, const std::string &tDir="", const std::string &stream="")
Simplify the retrieval of registered histograms of any type.
Definition AthHistogramming.cxx:183

◆ book() [4/8]

StatusCode AthHistogramming::book ( const TTree & treeRef,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 405 of file AthHistogramming.h.

406{
407 // Call the other Book method and see if it returns a valid pointer
408 TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
409 if ( treePointer )
410 {
411 return StatusCode::SUCCESS;
412 }
413 else
414 {
415 return StatusCode::FAILURE;
416 }
417}
AthHistogramming::bookGetPointer
TH1 * bookGetPointer(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.
Definition AthHistogramming.h:262

◆ book() [5/8]

StatusCode AthHistogramming::book ( TEfficiency & effRef,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 358 of file AthHistogramming.h.

359{
360 // Call the other Book method and see if it returns a valid pointer
361 TEfficiency* effPointer = this->bookGetPointer( effRef, tDir, stream );
362 if ( !effPointer ) {
363 m_msg << MSG::ERROR << "Couldn't book a TEfficiency" << endmsg;
364 return StatusCode::FAILURE;
365 }
366 return StatusCode::SUCCESS;
367}

◆ book() [6/8]

StatusCode AthHistogramming::book ( TEfficiency * eff,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 348 of file AthHistogramming.h.

349{
350 if ( !eff ) {
351 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency" << endmsg;
352 return StatusCode::FAILURE;
353 }
354 return this->book( *eff, tDir, stream );
355}

◆ book() [7/8]

StatusCode AthHistogramming::book ( TH1 & histRef,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 326 of file AthHistogramming.h.

327{
328 // Call the other Book method and see if it returns a valid pointer
329 TH1* histPointer = this->bookGetPointer( histRef, tDir, stream );
330 if ( !histPointer ) {
331 m_msg << MSG::ERROR << "Couldn't book a TH1" << endmsg;
332 return StatusCode::FAILURE;
333 }
334 return StatusCode::SUCCESS;
335}

◆ book() [8/8]

StatusCode AthHistogramming::book ( TH1 * hist,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 316 of file AthHistogramming.h.

317{
318 if ( !hist ) {
319 m_msg << MSG::ERROR << "Got a zero pointer to a TH1" << endmsg;
320 return StatusCode::FAILURE;
321 }
322 return this->book( *hist, tDir, stream );
323}

◆ book_histograms()

StatusCode EFTrackingSmearingAlg::book_histograms ( )
private

Definition at line 545 of file EFTrackingSmearingAlg.cxx.

546{
547 double b_eta[3]={50,-5,5};
548 CHECK(book(new TH1F("track_input_eta","#eta of input tracks",b_eta[0],b_eta[1],b_eta[2])));
549 CHECK(book(new TH1F("track_input_theta","#Theta of input tracks",50,0.0,4.0)));
550 CHECK(book(new TH1F("track_input_pt" ,"p_T of input tracks",50.0,-10.0,10.0)));
551 CHECK(book(new TH1F("track_input_phi","#phi of input tracks",50,-6.28,6.28)));
552 CHECK(book(new TH1F("track_input_z0" ,"z_0 of input tracks",100,-50,50)));
553 CHECK(book(new TH1F("track_input_d0" ,"d_0 of input tracks",50,-5,5)));
554
555 CHECK(book(new TH1F("track_input_sigma_theta","#sigma_{#Theta} of input tracks" ,50,0.0,0.001)));
556 CHECK(book(new TH1F("track_input_sigma_qOverP" ,"#sigma_{q/P} of input tracks" ,50.0,0.0,0.0001)));
557 CHECK(book(new TH1F("track_input_sigma_phi" ,"#sigma_{#phi} of input tracks" ,50,0.0,0.002)));
558 CHECK(book(new TH1F("track_input_sigma_z0" ,"#sigma_{z_0} of input tracks" ,100,0.0,5)));
559 CHECK(book(new TH1F("track_input_sigma_d0" ,"#sigma_{d_0} of input tracks" ,50,0.0,5)));
560
561 int maxtracks=500;
562 CHECK(book(new TH1F("n_input_tracks" ,"Number of input tracks" ,maxtracks,0,maxtracks)));
563 CHECK(book(new TH1F("n_output_tracks" ,"Number of output tracks" ,maxtracks,0,maxtracks)));
564 CHECK(book(new TH1F("n_output_narrow_tracks","Number of output tracks (narrow)" ,maxtracks,0,maxtracks)));
565 CHECK(book(new TH1F("n_output_broad_tracks" ,"Number of output tracks (broad)" ,maxtracks,0,maxtracks)));
566
567 CHECK(book(new TH1F("track_output_eta","#eta of output tracks",50,-5,5)));
568 CHECK(book(new TH1F("track_output_theta","#Theta of output tracks",50,0.0,4.0)));
569 CHECK(book(new TH1F("track_output_pt" ,"p_T of output tracks [GeV]",50,-10,10)));
570 CHECK(book(new TH1F("track_output_phi","#phi of output tracks",50,-6.28,6.28)));
571 CHECK(book(new TH1F("track_output_z0" ,"z_0 of output tracks",100,-50,50)));
572 CHECK(book(new TH1F("track_output_d0" ,"d_0 of output tracks",50,-5,5)));
573
574 CHECK(book(new TH1F("track_outputcoll_eta","#eta of output tracks collection",50,-5,5)));
575 CHECK(book(new TH1F("track_outputcoll_theta","#Theta of output tracks collection",50,0.0,4.0)));
576 CHECK(book(new TH1F("track_outputcoll_pt" ,"p_T of output tracks collection [GeV]",50,-10,10)));
577 CHECK(book(new TH1F("track_outputcoll_phi","#phi of output tracks collection",50,-6.28,6.28)));
578 CHECK(book(new TH1F("track_outputcoll_z0" ,"z_0 of output tracks collection",100,-50,50)));
579 CHECK(book(new TH1F("track_outputcoll_d0" ,"d_0 of output tracks collection",50,-5,5)));
580
581 CHECK(book(new TH1F("track_outputcoll_sigma_theta","#sigma_{#Theta} of output tracks collection" ,50,0.0,0.001)));
582 CHECK(book(new TH1F("track_outputcoll_sigma_qOverP" ,"#sigma_{q/P} of output tracks collection",50.0,0.0,0.0001)));
583 CHECK(book(new TH1F("track_outputcoll_sigma_phi" ,"#sigma_{#phi} of output tracks collection" ,50,0.0,0.002)));
584 CHECK(book(new TH1F("track_outputcoll_sigma_z0" ,"#sigma_{z_0} of output tracks collection" ,100,0.0,5)));
585 CHECK(book(new TH1F("track_outputcoll_sigma_d0" ,"#sigma_{d_0} of output tracks collection" ,50,0.0,5)));
586
587 CHECK(book(new TH1F("track_delta_sigma_theta", "#sigma_{#Theta} of output tracks collection",100,-0.001,0.001)));
588 CHECK(book(new TH1F("track_delta_sigma_qOverP", "#sigma_{q/P} of output tracks collection",100,-0.0001,0.0001)));
589 CHECK(book(new TH1F("track_delta_sigma_phi", "#sigma_{#phi} of output tracks collection" ,100,-0.002,0.002)));
590 CHECK(book(new TH1F("track_delta_sigma_z0", "#sigma_{z_0} of output tracks collection" ,100,-1.,1.)));
591 CHECK(book(new TH1F("track_delta_sigma_d0", "#sigma_{d_0} of output tracks collection" ,100,-0.5,0.5)));
592
593 CHECK(book(new TH1F("track_delta_eta","tracks #Delta #eta",50,-1,1)));
594 CHECK(book(new TH1F("track_delta_pt","tracks #Delta #pt [GeV]",50,-2.,2.)));
595 CHECK(book(new TH1F("track_delta_crv" ,"tracks #Delta crv",50,-1,1)));
596 CHECK(book(new TH1F("track_delta_phi","tracks #Delta #phi",50,-0.5,0.5)));
597 CHECK(book(new TH1F("track_delta_z0" ,"tracks #Delta z_0 ",100,-10,10)));
598 CHECK(book(new TH1F("track_delta_d0" ,"tracks #Delta d_0 ",50,-2,2)));
599 return StatusCode::SUCCESS;
600}
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition Control/AthenaKernel/AthenaKernel/errorcheck.h:422
python.TrigEgammaMonitorHelper.TH1F
TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
Definition TrigEgammaMonitorHelper.py:24

◆ bookGetPointer() [1/8]

TEfficiency * AthHistogramming::bookGetPointer ( const TEfficiency & eff,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 283 of file AthHistogramming.h.

284{
285 // We need to create a non-const clone
286 TEfficiency* histClone = dynamic_cast< TEfficiency* >( hist.Clone() );
287 if ( !histClone ) {
288 m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone in bookGetPointer" << endmsg;
289 return 0;
290 }
291 return this->bookGetPointer( *histClone, tDir, stream );
292
293}

◆ bookGetPointer() [2/8]

TGraph * AthHistogramming::bookGetPointer ( const TGraph & graphRef,
std::string tDir = "",
std::string stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

Definition at line 412 of file AthHistogramming.cxx.

413{
414 // Get a pointer
415 const TGraph* graphPointer = &graphRef;
416
417 // Check that we got a valid pointer
418 if ( !graphPointer )
419 {
420 m_msg << MSG::WARNING
421 << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
422 << "!" << endmsg;
423 return NULL;
424 }
425
426 // Modify the name and title according to the prefixes of this classes instance
427 std::string graphName = graphPointer->GetName();
428 const std::string graphTitle = graphPointer->GetTitle();
429
430 // Check if the hash for this graphName already exists, i.e., if we have a hash collision
431 const hash_t graphHash = this->hash(graphName);
432 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
433 if ( it != m_graphMap.end() ) // It does exist!
434 {
435 m_msg << MSG::WARNING
436 << "Detected a hash collision. The hash for the TGraph with name=" << graphName
437 << " already exists and points to a TGraph with name=" << it->second->GetName()
438 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
439 return NULL;
440 }
441
442 // Create a clone that has the new name
443 TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
444 if( !graphClone )
445 {
446 m_msg << MSG::WARNING
447 << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
448 << "!" << endmsg;
449 return NULL;
450 }
451 graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
452
453 // Massage the final string to book things
454 std::string bookingString("");
455 this->buildBookingString( bookingString, graphName, tDir, stream );
456
457 // Register the TGraph into the THistSvc
458 if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
459 {
460 m_msg << MSG::WARNING
461 << "Problem registering TGraph with name " << graphName
462 << ", title " << graphTitle
463 << " in " << m_name << "!" << endmsg;
464 return NULL;
465 }
466
467 // Also register it in the local map of string to pointer
468 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
469
470 return graphClone;
471}
AthHistogramming::histSvc
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
Definition AthHistogramming.h:423
AthHistogramming::m_histNamePostfix
std::string m_histNamePostfix
The postfix for the histogram THx name.
Definition AthHistogramming.h:237
AthHistogramming::hash
hash_t hash(const std::string &histName) const
Method to calculate a 32-bit hash from a string.
Definition AthHistogramming.h:430
AthHistogramming::hash_t
uint32_t hash_t
typedef for the internal hash
Definition AthHistogramming.h:171
AthHistogramming::m_histTitlePostfix
std::string m_histTitlePostfix
The postfix for the histogram THx title.
Definition AthHistogramming.h:243
AthHistogramming::m_histTitlePrefix
std::string m_histTitlePrefix
The prefix for the histogram THx title.
Definition AthHistogramming.h:240
AthHistogramming::m_name
std::string m_name
Instance name.
Definition AthHistogramming.h:247
AthHistogramming::buildBookingString
void buildBookingString(std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
Method to build individual booking string.
Definition AthHistogramming.cxx:545
AthHistogramming::m_histNamePrefix
std::string m_histNamePrefix
The prefix for the histogram THx name.
Definition AthHistogramming.h:234
AthHistogramming::m_graphMap
GraphMap_t m_graphMap
The map of TGraph names to their pointers.
Definition AthHistogramming.h:224

◆ bookGetPointer() [3/8]

TH1 * AthHistogramming::bookGetPointer ( const TH1 & hist,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 262 of file AthHistogramming.h.

263{
264 // We need to create a non-const clone
265 TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
266 if ( !histClone ) {
267 m_msg << MSG::ERROR << "Couldn't create a TH1 clone in bookGetPointer" << endmsg;
268 return 0;
269 }
270 return this->bookGetPointer( *histClone, tDir, stream );
271
272}

◆ bookGetPointer() [4/8]

TTree * AthHistogramming::bookGetPointer ( const TTree & treeRef,
std::string tDir = "",
std::string stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 297 of file AthHistogramming.cxx.

298{
299 // Get a pointer
300 const TTree* treePointer = &treeRef;
301
302 // Check that we got a valid pointer
303 if ( !treePointer )
304 {
305 m_msg << MSG::WARNING
306 << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
307 << "!" << endmsg;
308 return NULL;
309 }
310
311 // Modify the name and title according to the prefixes of this classes instance
312 std::string treeName = treePointer->GetName();
313 const std::string treeTitle = treePointer->GetTitle();
314
315 // Check if the hash for this treeName already exists, i.e., if we have a hash collision
316 const hash_t treeHash = this->hash(treeName);
317 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
318 if ( it != m_treeMap.end() ) // It does exist!
319 {
320 m_msg << MSG::WARNING
321 << "Detected a hash collision. The hash for the TTree with name=" << treeName
322 << " already exists and points to a TTree with name=" << it->second->GetName()
323 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
324 return NULL;
325 }
326
327 // Create a clone that has the new name
328 TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
329 if( !treeClone )
330 {
331 m_msg << MSG::WARNING
332 << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
333 << "!" << endmsg;
334 return NULL;
335 }
336 treeClone->SetTitle (treeTitle.c_str());
337
338 // Massage the final string to book things
339 std::string bookingString("");
340 this->buildBookingString( bookingString, treeName, tDir, stream );
341
342 // Register the TTree into the THistSvc
343 if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
344 {
345 m_msg << MSG::WARNING
346 << "Problem registering TTree with name " << treeName
347 << ", title " << treeTitle
348 << " in " << m_name << "!" << endmsg;
349 return NULL;
350 }
351
352 // Also register it in the local map of string to pointer
353 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
354
355 return treeClone;
356}
AthHistogramming::m_treeMap
TreeMap_t m_treeMap
The map of TTree names to their pointers.
Definition AthHistogramming.h:217

◆ bookGetPointer() [5/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency & effRef,
std::string tDir = "",
std::string stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 131 of file AthHistogramming.cxx.

132{
133 // Modify the name and title according to the prefixes of this classes instance
134 std::string effName(effRef.GetName());
135 const std::string effTitle(effRef.GetTitle());
136 std::string bookingString("");
137
138 this->buildBookingString( bookingString, effName, tDir, stream );
139 effRef.SetTitle((m_histTitlePrefix+effTitle+m_histTitlePostfix).c_str() );
140 effRef.SetName(effName.c_str());
141
142 // Check if the hash for this effName already exists, i.e., if we have a hash collision
143 const hash_t effHash = this->hash(effName);
144 EffMap_t::const_iterator it = m_effMap.find( effHash );
145 if ( it != m_effMap.end() ) // It does exist!
146 {
147 m_msg << MSG::WARNING
148 << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
149 << " already exists and points to a TEfficiency with name=" << it->second->GetName()
150 << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
151 return NULL;
152 }
153
154 // Set the new name and title for the TEfficiency, based on the prefixes that the user set for this class instance
155 // Create a clone that has the new name
156
157 // Massage the final string to book things
158
159 // Register the TEfficiency into the THistSvc
160 if ( !((histSvc()->regEfficiency(bookingString, &effRef)).isSuccess()) )
161 {
162 m_msg << MSG::WARNING
163 << "Problem registering TEfficiency with name " << effName
164 << ", name prefix " << m_histNamePrefix
165 << ", title " << effTitle
166 << ", tile prefix " << m_histTitlePrefix
167 << ", and tile postfix " << m_histTitlePostfix
168 << " in " << m_name << "!" << endmsg;
169 return NULL;
170 }
171
172 // Also register it in the local map of string to pointer
173 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, &effRef ) );
174
175 return &effRef;
176}
AthHistogramming::m_effMap
EffMap_t m_effMap
The map of histogram names to their pointers.
Definition AthHistogramming.h:210

◆ bookGetPointer() [6/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency * eff,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 295 of file AthHistogramming.h.

296{
297 if ( !hist ) {
298 m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency in bookGetPointer" << endmsg;
299 return 0;
300 }
301 return this->bookGetPointer( *hist, tDir, stream );
302}

◆ bookGetPointer() [7/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 & histRef,
std::string tDir = "",
std::string stream = "" )
protectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 83 of file AthHistogramming.cxx.

84{
85 // Modify the name and title according to the prefixes of this classes instance
86 std::string histName(histRef.GetName());
87 const std::string histTitle(histRef.GetTitle());
88 std::string bookingString("");
89
90 this->buildBookingString( bookingString, histName, tDir, stream );
91 histRef.SetTitle((m_histTitlePrefix+histTitle+m_histTitlePostfix).c_str() );
92 histRef.SetName(histName.c_str());
93
94 // Check if the hash for this histName already exists, i.e., if we have a hash collision
95 const hash_t histHash = this->hash(histName);
96 HistMap_t::const_iterator it = m_histMap.find( histHash );
97 if ( it != m_histMap.end() ) // It does exist!
98 {
99 m_msg << MSG::WARNING
100 << "Detected a hash collision. The hash for the histogram with name=" << histName
101 << " already exists and points to a histogram with name=" << it->second->GetName()
102 << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
103 return NULL;
104 }
105
106 // Set the new name and title for the histogram, based on the prefixes that the user set for this class instance
107 // Create a clone that has the new name
108
109 // Massage the final string to book things
110
111 // Register the histogram into the THistSvc
112 if ( !((histSvc()->regHist(bookingString, &histRef)).isSuccess()) )
113 {
114 m_msg << MSG::WARNING
115 << "Problem registering histogram with name " << histName
116 << ", name prefix " << m_histNamePrefix
117 << ", title " << histTitle
118 << ", tile prefix " << m_histTitlePrefix
119 << ", and tile postfix " << m_histTitlePostfix
120 << " in " << m_name << "!" << endmsg;
121 return NULL;
122 }
123
124 // Also register it in the local map of string to pointer
125 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, &histRef ) );
126
127 return &histRef;
128}
AthHistogramming::m_histMap
HistMap_t m_histMap
The map of histogram names to their pointers.
Definition AthHistogramming.h:203

◆ bookGetPointer() [8/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 * hist,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 274 of file AthHistogramming.h.

275{
276 if ( !hist ) {
277 m_msg << MSG::ERROR << "Got a zero pointer to a TH1 in bookGetPointer" << endmsg;
278 return 0;
279 }
280 return this->bookGetPointer( *hist, tDir, stream );
281}

◆ buildBookingString()

void AthHistogramming::buildBookingString ( std::string & bookingString,
std::string & histName,
std::string & tDir,
std::string & stream,
bool usePrefixPostfix = false )
privateinherited

Method to build individual booking string.

Definition at line 545 of file AthHistogramming.cxx.

550{
551 // Massage the final string to book things
552 if(tDir.empty()) tDir = m_rootDir;
553 size_t pos = histName.rfind('/');
554 if(pos != std::string::npos){
555 tDir+='/';
556 tDir.append(histName, 0,pos);
557 histName.erase(0,pos+1);
558 };
559 if(stream.empty()) stream = m_streamName;
560
561 if(usePrefixPostfix){
562 bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
563 } else {
564 bookingString = "/"+stream+"/"+tDir+"/"+histName;
565 }
566 while(bookingString.find("//") != std::string::npos){
567 this->myReplace(bookingString,"//","/");
568 }
569
570 return;
571}
AthHistogramming::m_rootDir
std::string m_rootDir
Name of the ROOT directory.
Definition AthHistogramming.h:231
AthHistogramming::m_streamName
std::string m_streamName
Name of the ROOT output stream (file).
Definition AthHistogramming.h:228
AthHistogramming::myReplace
void myReplace(std::string &str, const std::string &oldStr, const std::string &newStr)
Helper method to replace sub-string.
Definition AthHistogramming.cxx:575

◆ configAthHistogramming()

StatusCode AthHistogramming::configAthHistogramming ( const ServiceHandle< ITHistSvc > & histSvc,
const std::string & prefix,
const std::string & rootDir,
const std::string & histNamePrefix,
const std::string & histNamePostfix,
const std::string & histTitlePrefix,
const std::string & histTitlePostfix )
protectedinherited

To be called by the derived classes to fill the internal configuration.

Definition at line 51 of file AthHistogramming.cxx.

55{
56 m_histSvc = histSvc;
57 m_streamName = prefix;
58 m_rootDir = rootDir;
59 m_histNamePrefix = histNamePrefix;
60 m_histNamePostfix = histNamePostfix;
61 m_histTitlePrefix = histTitlePrefix;
62 m_histTitlePostfix = histTitlePostfix;
63
64 return StatusCode::SUCCESS;
65}
AthHistogramming::m_histSvc
ServiceHandle< ITHistSvc > m_histSvc
Pointer to the THistSvc (event store by default).
Definition AthHistogramming.h:196

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl,
const SG::VarHandleKeyType &  )
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158 {
159 return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160 hndl.value(),
161 hndl.documentation());
162
163 }
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition AthCommonDataStore.h:145

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
147 return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148 }
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition AthCommonDataStore.h:156

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

◆ efficiency()

TEfficiency * AthHistogramming::efficiency ( const std::string & effName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered TEfficiency.

Definition at line 235 of file AthHistogramming.cxx.

236{
237 // Build a 32 bit hash out of the name
238 const hash_t effHash = this->hash(effName);
239
240 // See if this entry exists in the map
241 EffMap_t::const_iterator it = m_effMap.find( effHash );
242 if ( it == m_effMap.end() ) // It doesn't exist!
243 { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
244
245 // Need to copy the strings as we will massage them from here on
246 std::string effNameCopy = effName;
247 std::string tDirCopy = tDir;
248 std::string streamCopy = stream;
249
250 // Massage the final string to book things
251 std::string bookingString("");
252 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
253
254 TEfficiency* effPointer(NULL);
255 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
256 {
257 // Massage the final string to book things
258 std::string bookingString("");
259 this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
260
261 if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
262 {
263 m_msg << MSG::WARNING
264 << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
265 << m_histNamePrefix + effNameCopy + m_histNamePostfix
266 << " or with name " << effNameCopy
267 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
268 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
269 return NULL;
270 }
271 // If we get to here, we actually found the TEfficiency in the THistSvc.
272 // So let's add it to the local cache map and return its pointer
273 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
274 return effPointer;
275 }
276 // If we get to here, we actually found the TEfficiency in the THistSvc.
277 // So let's add it to the local cache map and return its pointer
278 m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
279 return effPointer;
280 }
281
282 // Return the pointer to the TEfficiency that we got from the local cache map
283 return it->second;
284}
AthenaInternal::getEfficiency
std::pair< StatusCode, TEfficiency * > getEfficiency(ITHistSvc &svc, const std::string &name)
Definition AthenaPythonDict.h:59

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85{ return m_evtStore; }

◆ execute()

StatusCode EFTrackingSmearingAlg::execute ( )
overridevirtual

Definition at line 342 of file EFTrackingSmearingAlg.cxx.

342 {
343 ATH_MSG_DEBUG ("Executing " << name() << "...");
344
345 auto ctx = getContext() ;
346 if (m_smearTruthParticle)
347 return smearTruthParticles(ctx);
348
349 SG::ReadHandle<xAOD::TrackParticleContainer> inputTracks_handle( m_inputTrackParticleKey, ctx );
350 const xAOD::TrackParticleContainer* inputTracks = inputTracks_handle.cptr();
351 if (not inputTracks) {
352 ATH_MSG_FATAL("Unable to retrieve input ID tacks");
353 return StatusCode::FAILURE;
354 }
355
356 SG::WriteHandle<xAOD::TrackParticleContainer> outputTracks_handle( m_outputTrackParticleKey, ctx );
357 ATH_CHECK( outputTracks_handle.record( std::make_unique<xAOD::TrackParticleContainer>(), std::make_unique<xAOD::TrackParticleAuxContainer>() ) );
358 auto outputTracks = outputTracks_handle.ptr();
359
360 // clear the smearear
361 m_mySmearer->Clear();
362
363 //int trackno=0;
364 int n_input_tracks=0;
365 int n_output_tracks=0;
366 int n_output_broad_tracks=0;
367 int n_output_narrow_tracks=0;
368 ATH_MSG_DEBUG ("Found "<<inputTracks->size()<< " input tracks");
369 for ( const auto* trk : *inputTracks )
370 {
371 // get Cov matrix of input track
372 xAOD::ParametersCovMatrix_t trkcov = trk->definingParametersCovMatrix();
373 auto trkcovvec = trk->definingParametersCovMatrixVec();
374 double theta=trk->theta();
375 double pt = std::sin(theta)/std::abs(trk->qOverP());// MeV
376
377 ATH_MSG_DEBUG ("===> New Track: "
378 <<" curv=" << 1./trk->pt()
379 <<" phi=" << trk->phi0()
380 <<" eta=" << trk->eta()
381 <<" d0=" << trk->d0()
382 <<" z0=" << trk->z0()
383 <<" pT=" << trk->pt()
384 <<" cov_d0=" << trkcov(Trk::d0,Trk::d0)
385 <<" cov_z0=" << trkcov(Trk::z0,Trk::z0)
386 <<" sigma_d0=" << std::sqrt(std::abs(trkcov(Trk::d0,Trk::d0)))
387 <<" sigma_z0=" << std::sqrt(std::abs(trkcov(Trk::z0,Trk::z0))) );
388
389
390
391 if (std::abs(pt) > m_inputTracksPtCut)
392 {
393 n_input_tracks++;
394 if (m_enableMonitoring) {
395 hist("track_input_eta")->Fill(trk->eta());
396 hist("track_input_theta")->Fill(trk->theta());
397 hist("track_input_pt" )->Fill(pt/1000.);
398 hist("track_input_phi")->Fill(trk->phi0());
399 hist("track_input_z0" )->Fill(trk->z0());
400 hist("track_input_d0" )->Fill(trk->d0());
401
402 hist("track_input_sigma_theta") ->Fill(std::sqrt(trkcov(Trk::theta,Trk::theta)));
403 hist("track_input_sigma_qOverP")->Fill(std::sqrt(trkcov(Trk::qOverP,Trk::qOverP)));
404 hist("track_input_sigma_phi") ->Fill(std::sqrt(trkcov(Trk::phi0,Trk::phi0)));
405 hist("track_input_sigma_z0") ->Fill(std::sqrt(trkcov(Trk::z0,Trk::z0)));
406 hist("track_input_sigma_d0") ->Fill(std::sqrt(trkcov(Trk::d0,Trk::d0)));
407 }
408
409 // get Cov matrix of input track
410 auto trkcovvec = trk->definingParametersCovMatrixVec();
411 double qoverPt = trk->charge()*1000./pt; //this must be in GeV
412 m_mySmearer->AddTrack(trk->d0(),trk->z0(),qoverPt,trk->eta(),trk->phi0());
413 n_output_tracks += m_mySmearer->GetNTracks();
414
415 ATH_MSG_DEBUG ("Looping on output tracks #"<< m_mySmearer->GetNTracks());
416 for (const auto& otrack : m_mySmearer->Tracks)
417 {
418 xAOD::TrackParticle * newtrk = new xAOD::TrackParticle(*trk);
419 outputTracks->push_back(newtrk);
420 *newtrk = *trk;
421
422 double sintheta=std::sin(otrack.theta());
423 trkcov = trk->definingParametersCovMatrix();
424 auto newtrkcov = trkcov;
425 if (m_SigmaScaleFactor !=0) {
426 // modify the cov matrix
427 ATH_MSG_DEBUG ("Setting parameters in covariance vector");
428 for (int ii=0;ii<5;ii++) for (int jj=0;jj<5;jj++) {
429 // correct the CM to be consistent with the smearing ofthe parameters
430 newtrkcov(ii,jj)=trkcov(ii,jj) * (std::pow(m_SigmaScaleFactor,2) + 1);
431 }
432 trkcovvec.clear();
433 EigenHelpers::eigenMatrixToVector(trkcovvec,newtrkcov,"");
434 newtrk->setDefiningParametersCovMatrixVec(trkcovvec);
435 ATH_MSG_DEBUG ("Setting parameters covariance");
436 newtrk->setDefiningParametersCovMatrix(newtrkcov);
437
438 /* old method stays for reference, since it's not yet validated
439 for (int ii=0;ii<5;ii++) for (int jj=0;jj<5;jj++) if (ii!=jj) trkcov(ii,jj)=0.0;
440 trkcov(Trk::d0 ,Trk::d0 ) = otrack.sigma_d0() * otrack.sigma_d0() ;
441 trkcov(Trk::z0 ,Trk::z0 ) = otrack.sigma_z0() * otrack.sigma_z0() ;
442 trkcov(Trk::theta ,Trk::theta ) = otrack.sigma_theta()* otrack.sigma_theta();
443 trkcov(Trk::phi0 ,Trk::phi0 ) = otrack.sigma_phi() * otrack.sigma_phi() ;
444 trkcov(Trk::qOverP,Trk::qOverP) = otrack.sigma_curv()/1000. * otrack.sigma_curv()/1000. *sintheta*sintheta;
445
446 trkcovvec.clear();
447 EigenHelpers::eigenMatrixToVector(trkcovvec,trkcov,"");
448 newtrk->setDefiningParametersCovMatrixVec(trkcovvec);
449 ATH_MSG_DEBUG ("Setting parameters covariance");
450 newtrk->setDefiningParametersCovMatrix(trkcov);
451 */
452
453 //(float d0, float z0, float phi0, float theta, float qOverP)
454 newtrk->setDefiningParameters(
455 otrack.d0(),
456 otrack.z0(),
457 otrack.phi(),
458 otrack.theta(),
459 (otrack.curv()*sintheta/1000.) //oneOverp in MeV
460 );
461 }
462 xAOD::ParametersCovMatrix_t trkcov_out = newtrk->definingParametersCovMatrix();
463
464
465 ATH_MSG_DEBUG ("Smeared Track: "
466 <<" curv=" << 1./newtrk->pt()
467 <<" phi=" << newtrk->phi()
468 <<" eta=" << newtrk->eta()
469 <<" d0=" << newtrk->d0()
470 <<" z0=" << newtrk->z0()
471 <<" pT=" << newtrk->pt()
472 <<" cov_d0=" << trkcov_out(Trk::d0,Trk::d0)
473 <<" cov_z0=" << trkcov_out(Trk::z0,Trk::z0)
474 <<" sigma_d0=" << std::sqrt(std::abs(trkcov_out(Trk::d0,Trk::d0)))
475 <<" sigma_z0=" << std::sqrt(std::abs(trkcov_out(Trk::z0,Trk::z0))) );
476
477
478 if (m_enableMonitoring) {
479 hist("track_output_eta")->Fill(otrack.eta());
480 hist("track_output_theta")->Fill(otrack.theta());
481 hist("track_output_pt" )->Fill(trk->charge()*otrack.pt() );
482 hist("track_output_phi")->Fill(otrack.phi());
483 hist("track_output_z0" )->Fill(otrack.z0() );
484 hist("track_output_d0" )->Fill(otrack.d0() );
485
486 hist("track_outputcoll_eta")->Fill(newtrk->eta());
487 hist("track_outputcoll_theta")->Fill(newtrk->theta());
488 hist("track_outputcoll_pt" )->Fill(trk->charge()*std::sin(newtrk->theta())/(1000.0*newtrk->qOverP()));
489 hist("track_outputcoll_phi")->Fill(newtrk->phi0());
490 hist("track_outputcoll_z0" )->Fill(newtrk->z0());
491 hist("track_outputcoll_d0" )->Fill(newtrk->d0());
492
493 hist("track_delta_eta")->Fill(newtrk->eta() - trk->eta());
494 hist("track_delta_pt")->Fill((newtrk->pt() - trk->pt())/1000.);
495 hist("track_delta_crv")->Fill((1000.0*newtrk->qOverP()/std::sin(theta))-((1000.0*trk->qOverP())/std::sin(theta)));
496 hist("track_delta_phi")->Fill(newtrk->phi() - trk->phi0());
497 hist("track_delta_z0" )->Fill(newtrk->z0() - trk->z0());
498 hist("track_delta_d0" )->Fill(newtrk->d0() - trk->d0());
499
500 auto trkcov_original = trk->definingParametersCovMatrix();
501 hist("track_delta_sigma_theta") ->Fill(std::sqrt(trkcov_out(Trk::theta,Trk::theta)) - std::sqrt(trkcov_original(Trk::theta,Trk::theta)));
502 hist("track_delta_sigma_qOverP")->Fill(std::sqrt(trkcov_out(Trk::qOverP,Trk::qOverP)) - std::sqrt(trkcov_original(Trk::qOverP,Trk::qOverP)));
503 hist("track_delta_sigma_phi") ->Fill(std::sqrt(trkcov_out(Trk::phi0,Trk::phi0)) - std::sqrt(trkcov_original(Trk::phi0,Trk::phi0)));
504 hist("track_delta_sigma_z0") ->Fill(std::sqrt(trkcov_out(Trk::z0,Trk::z0)) - std::sqrt(trkcov_original(Trk::z0,Trk::z0)));
505 hist("track_delta_sigma_d0") ->Fill(std::sqrt(trkcov_out(Trk::d0,Trk::d0)) - std::sqrt(trkcov_original(Trk::d0,Trk::d0)));
506
507
508 hist("track_outputcoll_sigma_theta") ->Fill(std::sqrt(trkcov(Trk::theta,Trk::theta)));
509 hist("track_outputcoll_sigma_qOverP")->Fill(std::sqrt(trkcov(Trk::qOverP,Trk::qOverP)));
510 hist("track_outputcoll_sigma_phi") ->Fill(std::sqrt(trkcov(Trk::phi0,Trk::phi0)));
511 hist("track_outputcoll_sigma_z0") ->Fill(std::sqrt(trkcov(Trk::z0,Trk::z0)));
512 hist("track_outputcoll_sigma_d0") ->Fill(std::sqrt(trkcov(Trk::d0,Trk::d0)));
513
514 }
515 // do we need this hack?
516 // hack!!!
517 // Chi2: 1 -> "Broad Fake"
518 // 0 -> "core tracks"
519 //
520 // NDF: index of parent track in input track collection
521 //newtrk->setFitQuality(otrack.FakeFlags(),trackno);
522
523
524 }
525
526 }
527 m_mySmearer->Clear(); // clear teh smearer after each input track
528 //trackno++;
529 }
530
531 ATH_MSG_DEBUG ("End of loop track #"<<n_input_tracks<<" ---> "<<" "<< n_output_tracks
532 <<" "<<n_output_narrow_tracks<<" "<<n_output_broad_tracks);
533 if (m_enableMonitoring) {
534 hist("n_input_tracks")->Fill(n_input_tracks);
535 hist("n_output_tracks")->Fill(n_output_tracks);
536 hist("n_output_narrow_tracks")->Fill(n_output_narrow_tracks);
537 hist("n_output_broad_tracks")->Fill(n_output_broad_tracks);
538 }
539 return StatusCode::SUCCESS;
540}
theta
Scalar theta() const
theta method
Definition AmgMatrixBasePlugin.h:75
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition AthMsgStreamMacros.h:34
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
EFTrackingSmearingAlg::m_enableMonitoring
BooleanProperty m_enableMonitoring
Definition EFTrackingSmearingAlg.h:55
EFTrackingSmearingAlg::m_inputTrackParticleKey
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_inputTrackParticleKey
Definition EFTrackingSmearingAlg.h:29
EFTrackingSmearingAlg::smearTruthParticles
StatusCode smearTruthParticles(const EventContext &ctx)
Definition EFTrackingSmearingAlg.cxx:204
EFTrackingSmearingAlg::m_mySmearer
std::unique_ptr< FakeTrackSmearer > m_mySmearer
Definition EFTrackingSmearingAlg.h:72
EFTrackingSmearingAlg::m_smearTruthParticle
BooleanProperty m_smearTruthParticle
Definition EFTrackingSmearingAlg.h:56
EFTrackingSmearingAlg::m_outputTrackParticleKey
SG::WriteHandleKey< xAOD::TrackParticleContainer > m_outputTrackParticleKey
Definition EFTrackingSmearingAlg.h:32
xAOD::TrackParticle_v1::setDefiningParametersCovMatrix
void setDefiningParametersCovMatrix(const ParametersCovMatrix_t &cov)
Set the defining parameters covariance matrix.
Definition TrackParticle_v1.cxx:227
xAOD::TrackParticle_v1::z0
float z0() const
Returns the parameter.
xAOD::TrackParticle_v1::theta
float theta() const
Returns the parameter, which has range 0 to .
xAOD::TrackParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .).
xAOD::TrackParticle_v1::definingParametersCovMatrix
const ParametersCovMatrix_t definingParametersCovMatrix() const
Returns the 5x5 symmetric matrix containing the defining parameters covariance matrix.
Definition TrackParticle_v1.cxx:260
xAOD::TrackParticle_v1::d0
float d0() const
Returns the parameter.
xAOD::TrackParticle_v1::setDefiningParameters
void setDefiningParameters(float d0, float z0, float phi0, float theta, float qOverP)
Set the defining parameters.
Definition TrackParticle_v1.cxx:181
xAOD::TrackParticle_v1::qOverP
float qOverP() const
Returns the parameter.
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition TrackParticle_v1.cxx:75
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TrackParticle_v1.cxx:79
xAOD::TrackParticle_v1::phi0
float phi0() const
Returns the parameter, which has range to .
Definition TrackParticle_v1.cxx:150
xAOD::TrackParticle_v1::setDefiningParametersCovMatrixVec
void setDefiningParametersCovMatrixVec(const std::vector< float > &cov)
Definition TrackParticle_v1.cxx:412
EigenHelpers::eigenMatrixToVector
static void eigenMatrixToVector(VECTOR &vec, COVARIANCE &cov, const char *)
Definition EFTrackingSmearingAlg.cxx:87
Trk::phi0
@ phi0
Definition ParamDefs.h:65
Trk::theta
@ theta
Definition ParamDefs.h:66
Trk::qOverP
@ qOverP
perigee
Definition ParamDefs.h:67
Trk::d0
@ d0
Definition ParamDefs.h:63
Trk::z0
@ z0
Definition ParamDefs.h:64
test_pyathena.pt
pt
Definition test_pyathena.py:11
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::TrackParticleContainer
TrackParticleContainer_v1 TrackParticleContainer
Definition of the current "TrackParticle container version".
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticleContainer.h:14

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

51{
52 // If we didn't find any symlinks to add, just return the collection
53 // from the base class. Otherwise, return the extended collection.
54 if (!m_extendedExtraObjects.empty()) {
55 return m_extendedExtraObjects;
56 }
57 return Algorithm::extraOutputDeps();
58}
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition AthAlgorithm.h:79

◆ graph()

TGraph * AthHistogramming::graph ( const std::string & graphName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered TGraphs.

Definition at line 477 of file AthHistogramming.cxx.

478{
479 // Build a 32 bit hash out of the name
480 const hash_t graphHash = this->hash(graphName);
481
482 // See if this entry exists in the map
483 GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
484 if ( it == m_graphMap.end() ) // It doesn't exist!
485 { // Let's see into the THistSvc if somebody else has registered the TGraph...
486
487 // Need to copy the strings as we will massage them from here on
488 std::string graphNameCopy = graphName;
489 std::string tDirCopy = tDir;
490 std::string streamCopy = stream;
491
492 // Massage the final string to book things
493 std::string bookingString("");
494 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true);
495
496 TGraph* graphPointer(NULL);
497 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
498 {
499 // Massage the final string to book things
500 std::string bookingString("");
501 this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false );
502
503 if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
504 {
505 m_msg << MSG::WARNING
506 << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
507 << m_histNamePrefix + graphNameCopy + m_histNamePostfix
508 << " or with name " << graphNameCopy
509 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
510 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
511 return NULL;
512 }
513 // If we get to here, we actually found the TGraph in the THistSvc.
514 // So let's add it to the local cache map and return its pointer
515 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
516 return graphPointer;
517 }
518 // If we get to here, we actually found the TGraph in the THistSvc.
519 // So let's add it to the local cache map and return its pointer
520 m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
521 return graphPointer;
522 }
523
524
525 // Return the pointer to the TGraph that we got from the local cache map
526 return it->second;
527}
AthenaInternal::getGraph
std::pair< StatusCode, TGraph * > getGraph(ITHistSvc &svc, const std::string &name)
Definition AthenaPythonDict.h:51

◆ hash()

AthHistogramming::hash_t AthHistogramming::hash ( const std::string & histName) const
inlineprivateinherited

Method to calculate a 32-bit hash from a string.

Definition at line 430 of file AthHistogramming.h.

431{
432 const uint64_t hash64 = CxxUtils::crc64( histName );
433 return (hash_t)(hash64 & 0xFFFFFFFF);
434}
CxxUtils::crc64
uint64_t crc64(const CRCTable &table, const char *data, size_t data_len)
Find the CRC-64 of a string,.
Definition crc64.cxx:696
xxh3::hash64
std::uint64_t hash64(const void *data, std::size_t size)
Passthrough to XXH3_64bits.
Definition XXH.cxx:9

◆ hist()

TH1 * AthHistogramming::hist ( const std::string & histName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered histograms of any type.

Definition at line 183 of file AthHistogramming.cxx.

184{
185 // Build a 32 bit hash out of the name
186 const hash_t histHash = this->hash(histName);
187
188 // See if this entry exists in the map
189 HistMap_t::const_iterator it = m_histMap.find( histHash );
190 if ( it == m_histMap.end() ) // It doesn't exist!
191 { // Let's see into the THistSvc if somebody else has registered the histogram...
192
193 // Need to copy the strings as we will massage them from here on
194 std::string histNameCopy = histName;
195 std::string tDirCopy = tDir;
196 std::string streamCopy = stream;
197
198 // Massage the final string to book things
199 std::string bookingString("");
200 this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy ,false);
201
202 TH1* histPointer(NULL);
203 if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
204 {
205 // Massage the final string to book things
206 std::string bookingString("");
207 this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy, true );
208
209 if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
210 {
211 m_msg << MSG::WARNING
212 << "Problem retrieving the histogram with name (including pre- and post-fixes) "
213 << m_histNamePrefix + histNameCopy + m_histNamePostfix
214 << " or with name " << histNameCopy
215 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
216 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
217 return NULL;
218 }
219 // If we get to here, we actually found the histogram in the THistSvc.
220 // So let's add it to the local cache map and return its pointer
221 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
222 return histPointer;
223 }
224 // If we get to here, we actually found the histogram in the THistSvc.
225 // So let's add it to the local cache map and return its pointer
226 m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
227 return histPointer;
228 }
229
230
231 // Return the pointer to the histogram that we got from the local cache map
232 return it->second;
233}
AthenaInternal::getHist
std::pair< StatusCode, TH1 * > getHist(ITHistSvc &svc, const std::string &name, size_t index=0)
Definition AthenaPythonDict.h:42

◆ hist2d()

TH2 * AthHistogramming::hist2d ( const std::string & histName,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 2-d histograms.

Definition at line 371 of file AthHistogramming.h.

372{
373 // Get the TH1 pointer
374 TH1* th1Pointer = this->hist(histName, tDir, stream);
375 if ( !th1Pointer )
376 {
377 m_msg << MSG::ERROR
378 << "Cannot get a 2-d histogram with name " << histName
379 << "... will probably seg-fault!" << endmsg;
380 return NULL;
381 }
382 // If the TH1 pointer is valid, simply return the dynamic_cast
383 return dynamic_cast<TH2*>( th1Pointer );
384}

◆ hist3d()

TH3 * AthHistogramming::hist3d ( const std::string & histName,
const std::string & tDir = "",
const std::string & stream = "" )
inlineprotectedinherited

Simplify the retrieval of registered 3-d histograms.

Definition at line 388 of file AthHistogramming.h.

389{
390 // Get the TH1 pointer
391 TH1* th1Pointer = this->hist(histName, tDir, stream);
392 if ( !th1Pointer )
393 {
394 m_msg << MSG::ERROR
395 << "Cannot get a 3-d histogram with name " << histName
396 << "... will probably seg-fault!" << endmsg;
397 return NULL;
398 }
399 // If the TH1 pointer is valid, simply return the dynamic_cast
400 return dynamic_cast<TH3*>( th1Pointer );
401}

◆ histSvc()

const ServiceHandle< ITHistSvc > & AthHistogramAlgorithm::histSvc ( ) const
inlineinherited

The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.

Definition at line 113 of file AthHistogramAlgorithm.h.

114{
115 return m_histSvc;
116}
AthHistogramAlgorithm::m_histSvc
ServiceHandle< ITHistSvc > m_histSvc
Default constructor: AthHistogramAlgorithm();.
Definition AthHistogramAlgorithm.h:83

◆ initialize()

StatusCode EFTrackingSmearingAlg::initialize ( )
overridevirtual

Definition at line 101 of file EFTrackingSmearingAlg.cxx.

101 {
102 ATH_MSG_INFO ("Initializing " << name() << "...");
103
104 ATH_CHECK( m_inputTrackParticleKey.initialize() );
105 ATH_CHECK( m_outputTrackParticleKey.initialize() );
106
107 ATH_CHECK( m_inputTruthParticleKey.initialize() );
108 ATH_CHECK( m_outputTruthParticleKey.initialize() );
109
110 // Decoration keys
111 m_d0DecoratorKey = m_outputTruthParticleKey.key() + ".d0";
112 m_z0DecoratorKey = m_outputTruthParticleKey.key() + ".z0";
113 m_ptDecoratorKey = m_outputTruthParticleKey.key() + ".pt";
114
115 ATH_CHECK(m_d0DecoratorKey.initialize());
116 ATH_CHECK(m_z0DecoratorKey.initialize());
117 ATH_CHECK(m_ptDecoratorKey.initialize());
118
119 ATH_MSG_INFO("########## EFTrackingSmearingAlg Configurations are ########## ");
120 ATH_MSG_INFO("------- InputTrackParticleKey: "<<m_inputTrackParticleKey.key());
121 ATH_MSG_INFO("------- OutputTrackParticleKey: "<<m_outputTrackParticleKey.key());
122 ATH_MSG_INFO("------- InputTruthParticleKey: " <<m_inputTruthParticleKey.key());
123 ATH_MSG_INFO("------- OutputTruthParticleKey: "<<m_outputTruthParticleKey.key());
124 ATH_MSG_INFO("------- inputTracksPtCut [GeV]: "<<m_inputTracksPtCut);
125 ATH_MSG_INFO("------- outputTracksPtCut [GeV]: "<<m_outputTracksPtCut);
126 ATH_MSG_INFO("------- SmearingSigma: "<<m_SigmaScaleFactor);
127 ATH_MSG_INFO("------- trackEfficiency: "<<m_smearedTrackEfficiency);
128 ATH_MSG_INFO("------- parameterizedEfficiency: "<<m_parameterizedTrackEfficiency);
129 ATH_MSG_INFO("------- parameterizedEfficiency LRT: "<<m_parameterizedTrackEfficiency_LRT);
130 ATH_MSG_INFO("------- parameterizedEfficiency LRT high d0 cut: "<<m_smearedTrackEfficiency_d0high_LRT );
131 ATH_MSG_INFO("------- parameterizedEfficiency LRT low d0 cut: "<<m_smearedTrackEfficiency_d0low_LRT );
132 ATH_MSG_INFO("------- UseResolutionPtCutOff: "<<m_UseResolutionPtCutOff);
133 ATH_MSG_INFO("------- SetResolutionPtCutOff: "<<m_SetResolutionPtCutOff);
134 ATH_MSG_INFO("------- EnableMonitoring:" <<m_enableMonitoring);
135 ATH_MSG_INFO("------- SmearTruthParticle:"<< m_smearTruthParticle);
136
137 ATH_MSG_INFO("------- IncludeDuplicatesAndFakes: "<<m_EnableFakes);
138 ATH_MSG_INFO("------- RandomSeed: " <<m_RandomSeed);
139 ATH_MSG_INFO("------- UseCoinToss: " <<m_UseCoinToss);
140 ATH_MSG_INFO("------- FakeKillerEnable: "<<m_FakeKillerEnable);
141 ATH_MSG_INFO("------- IncludeFakesInResolutionCalculation: "<<m_IncludeFakesInResolutionCalculation);
142 ATH_MSG_INFO("########## EFTrackingSmearingAlg Configurations: That's it. ########## ");
143
144 std::string smearerName;
145 if (m_smearTruthParticle){
146 ATH_MSG_INFO("Will output new truth container with name "<<m_outputTruthParticleKey.key());
147 smearerName = m_outputTruthParticleKey.key()+"_smearer";
148 }
149 else {
150 ATH_MSG_INFO("Will output new track container with name "<<m_outputTrackParticleKey.key());
151 smearerName = m_outputTrackParticleKey.key()+"_smearer";
152 }
153
154 // configure the Smearer
155 m_mySmearer = std::make_unique<FakeTrackSmearer>(smearerName, m_RandomSeed, msgLvl (MSG::DEBUG));
156 m_mySmearer->SetInputTracksPtCut(m_inputTracksPtCut);
157 m_mySmearer->SetOutputTracksPtCut(m_outputTracksPtCut);
158 m_mySmearer->SetTrackingEfficiency(m_smearedTrackEfficiency);
159 m_mySmearer->SetParameterizedEfficiency(m_parameterizedTrackEfficiency);
160 m_mySmearer->SetParameterizedEfficiency_LRT(m_parameterizedTrackEfficiency_LRT);
161 m_mySmearer->SetParameterizedEfficiency_highd0_LRT(m_smearedTrackEfficiency_d0high_LRT);
162 m_mySmearer->SetParameterizedEfficiency_lowd0_LRT(m_smearedTrackEfficiency_d0low_LRT);
163
164 m_mySmearer->SetSigmaScaleFactor(m_SigmaScaleFactor.value());
165 m_mySmearer->UseResolutionPtCutOff(m_UseResolutionPtCutOff.value());
166 m_mySmearer->SetResolutionPtCutOff(m_SetResolutionPtCutOff.value());
167
168 m_mySmearer->EnableFakes(m_EnableFakes.value());
169 m_mySmearer->UseCoinToss(m_UseCoinToss.value());
170 m_mySmearer->FakeKillerEnable(m_FakeKillerEnable.value());
171 m_mySmearer->IncludeFakesInResolutionCalculation(m_IncludeFakesInResolutionCalculation.value());
172
173 if (m_enableMonitoring) {
174 // store the smearing functions
175 TF1 *d0res_eta = m_mySmearer->d0res_eta;
176 TF1 *z0res_eta = m_mySmearer->z0res_eta;
177 TF1 *curvres_eta = m_mySmearer->curvres_eta;
178 TF1 *d0res_pt = m_mySmearer->d0res_pt;
179 TF1 *z0res_pt = m_mySmearer->z0res_pt;
180 TF1 *curvres_pt = m_mySmearer->curvres_pt;
181 TF1 *effLRT_d0 = m_mySmearer->effLRT_d0;
182
183 CHECK(book(new TH1F("d0res_function_vs_eta","#eta of track (p_{T}=10GeV);#eta",100, 0.0,4.0)));
184 CHECK(book(new TH1F("z0res_function_vs_eta","#eta of track (p_{T}=10GeV);#eta",100, 0.0,4.0)));
185 CHECK(book(new TH1F("curvres_function_vs_eta","#eta of track (p_{T}=10GeV);#eta",100, 0.0,4.0)));
186 CHECK(book(new TH1F("d0res_function_vs_pt","p_{T} of track (#eta=1);p_{T} [GeV]",100, 1.0,200.0)));
187 CHECK(book(new TH1F("z0res_function_vs_pt","p_{T} of track (#eta=1);p_{T} [GeV]",100, 1.0,200.0)));
188 CHECK(book(new TH1F("curvres_function_vs_pt","p_{T} of track (#eta=1);p_{T} [GeV]",100, 1.0,200.0)));
189 CHECK(book(new TH1F("effLRT_function_vs_d0","d_{0} of track;d_{0} [mm]",100, 0.001,600.0)));
190 hist("d0res_function_vs_eta")->Add(d0res_eta);
191 hist("z0res_function_vs_eta")->Add(z0res_eta);
192 hist("curvres_function_vs_eta")->Add(curvres_eta);
193 hist("d0res_function_vs_pt")->Add(d0res_pt);
194 hist("z0res_function_vs_pt")->Add(z0res_pt);
195 hist("curvres_function_vs_pt")->Add(curvres_pt);
196 hist("effLRT_function_vs_d0")->Add(effLRT_d0);
197 // book historgams
198 CHECK(book_histograms());
199 }
200 return StatusCode::SUCCESS;
201}
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
AthCommonMsg< Algorithm >::msgLvl
bool msgLvl(const MSG::Level lvl) const
Definition AthCommonMsg.h:30
EFTrackingSmearingAlg::m_IncludeFakesInResolutionCalculation
BooleanProperty m_IncludeFakesInResolutionCalculation
Definition EFTrackingSmearingAlg.h:65
EFTrackingSmearingAlg::m_SetResolutionPtCutOff
DoubleProperty m_SetResolutionPtCutOff
Definition EFTrackingSmearingAlg.h:50
EFTrackingSmearingAlg::m_smearedTrackEfficiency
DoubleProperty m_smearedTrackEfficiency
Definition EFTrackingSmearingAlg.h:53
EFTrackingSmearingAlg::m_ptDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_ptDecoratorKey
Definition EFTrackingSmearingAlg.h:44
EFTrackingSmearingAlg::m_z0DecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_z0DecoratorKey
Definition EFTrackingSmearingAlg.h:43
EFTrackingSmearingAlg::m_parameterizedTrackEfficiency
BooleanProperty m_parameterizedTrackEfficiency
Definition EFTrackingSmearingAlg.h:54
EFTrackingSmearingAlg::m_UseResolutionPtCutOff
BooleanProperty m_UseResolutionPtCutOff
Definition EFTrackingSmearingAlg.h:49
EFTrackingSmearingAlg::m_d0DecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_d0DecoratorKey
Definition EFTrackingSmearingAlg.h:42
EFTrackingSmearingAlg::m_smearedTrackEfficiency_d0low_LRT
DoubleProperty m_smearedTrackEfficiency_d0low_LRT
Definition EFTrackingSmearingAlg.h:60
EFTrackingSmearingAlg::m_FakeKillerEnable
BooleanProperty m_FakeKillerEnable
Definition EFTrackingSmearingAlg.h:64
EFTrackingSmearingAlg::m_inputTruthParticleKey
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_inputTruthParticleKey
Definition EFTrackingSmearingAlg.h:35
EFTrackingSmearingAlg::m_parameterizedTrackEfficiency_LRT
BooleanProperty m_parameterizedTrackEfficiency_LRT
Definition EFTrackingSmearingAlg.h:58
EFTrackingSmearingAlg::m_smearedTrackEfficiency_d0high_LRT
DoubleProperty m_smearedTrackEfficiency_d0high_LRT
Definition EFTrackingSmearingAlg.h:59
EFTrackingSmearingAlg::m_outputTruthParticleKey
SG::WriteHandleKey< xAOD::TruthParticleContainer > m_outputTruthParticleKey
Definition EFTrackingSmearingAlg.h:38

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ msg()

MsgStream & AthCommonMsg< Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

bool AthCommonMsg< Algorithm >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ myReplace()

void AthHistogramming::myReplace ( std::string & str,
const std::string & oldStr,
const std::string & newStr )
privateinherited

Helper method to replace sub-string.

Definition at line 575 of file AthHistogramming.cxx.

578{
579 size_t pos = 0;
580 while((pos = str.find(oldStr, pos)) != std::string::npos)
581 {
582 str.replace(pos, oldStr.length(), newStr);
583 pos += newStr.length();
584 }
585}

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
383 PBASE::renounce (h);
384 }
AthCommonDataStore::renounce
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)
Definition AthCommonDataStore.h:380

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
365 handlesArray.renounce();
366 }

◆ smearTruthParticles()

StatusCode EFTrackingSmearingAlg::smearTruthParticles ( const EventContext & ctx)
private

Definition at line 204 of file EFTrackingSmearingAlg.cxx.

204 {
205
206 SG::ReadHandle<xAOD::TruthParticleContainer> inputTruth_handle( m_inputTruthParticleKey, ctx );
207 const xAOD::TruthParticleContainer* inputTruth = inputTruth_handle.cptr();
208 if (not inputTruth) {
209 ATH_MSG_FATAL("Unable to retrieve input truth particle");
210 return StatusCode::FAILURE;
211 }
212 SG::WriteHandle<xAOD::TruthParticleContainer> outputTruth_handle( m_outputTruthParticleKey, ctx );
213 ATH_CHECK( outputTruth_handle.record( std::make_unique<xAOD::TruthParticleContainer>(), std::make_unique<xAOD::TruthParticleAuxContainer>() ) );
214 auto outputTruth = outputTruth_handle.ptr();
215
216 // create decorators
217 SG::WriteDecorHandle<xAOD::TruthParticleContainer, float > d0Decorator(m_d0DecoratorKey, ctx);
218 SG::WriteDecorHandle<xAOD::TruthParticleContainer, float > z0Decorator(m_z0DecoratorKey, ctx);
219 SG::WriteDecorHandle<xAOD::TruthParticleContainer, float > ptDecorator(m_ptDecoratorKey, ctx);
220
221
222 // clear the smearear
223 m_mySmearer->Clear();
224
225 static const SG::ConstAccessor<float> ptAcc("pt");
226 static const SG::ConstAccessor<float> d0Acc("d0");
227 static const SG::ConstAccessor<float> z0Acc("z0");
228 static const SG::ConstAccessor<float> thetaAcc("theta");
229
230 int n_input_tracks=0;
231 int n_output_tracks=0;
232 int n_output_broad_tracks=0;
233 int n_output_narrow_tracks=0;
234 ATH_MSG_DEBUG ("Found "<<inputTruth->size()<< " input truth particles");
235 for ( const auto* part : *inputTruth )
236 {
237 double pt = part->pt();// MeV
238 float theta = thetaAcc(*part);
239 float z0 = z0Acc(*part);
240 float d0 = d0Acc(*part);
241 float eta = part->eta();
242 float phi = part->phi();
243 if (part->isNeutral()) continue;
244 if (pt <=0.) continue;
245
246 ATH_MSG_DEBUG ("===> New Truth: "
247 <<" curv=" << 1./part->pt()
248 <<" phi=" << part->phi()
249 <<" eta=" << part->eta()
250 <<" d0=" << d0Acc(*part)
251 <<" z0=" << z0Acc(*part)
252 <<" pT=" << part->pt()
253 <<" PDGID=" << part->pdgId()
254 <<" status=" << part->status()
255 );
256 if (part->parent(0)) ATH_MSG_DEBUG (" parent status=" << part->parent(0)->pdgId());
257
258 if (std::abs(pt)/1000. > m_inputTracksPtCut) //GeV cut
259 {
260 n_input_tracks++;
261 if (m_enableMonitoring) {
262 hist("track_input_eta")->Fill(eta);
263 hist("track_input_theta")->Fill(theta);
264 hist("track_input_pt" )->Fill(pt/1000.);
265 hist("track_input_phi")->Fill(phi);
266 hist("track_input_z0" )->Fill(z0);
267 hist("track_input_d0" )->Fill(d0);
268 }
269 double qoverPt = part->charge()*1000./pt; //this must be in GeV
270 m_mySmearer->AddTrack(d0,z0,qoverPt,eta,phi); // smearing here
271 n_output_tracks += m_mySmearer->GetNTracks();
272
273 ATH_MSG_DEBUG ("Looping on output tracks #"<< m_mySmearer->GetNTracks());
274 for (const auto& otrack : m_mySmearer->Tracks)
275 {
276 xAOD::TruthParticle * newtrk = new xAOD::TruthParticle(*part);
277 outputTruth->push_back(newtrk);
278 *newtrk = *part;
279 auto newpt = part->pt();
280 if (m_SigmaScaleFactor !=0) {
281 // set the decorators
282 d0Decorator(*newtrk) = otrack.d0();
283 z0Decorator(*newtrk) = otrack.z0();
284 ptDecorator(*newtrk) = otrack.pt()*1000.; //MeV
285 //TrackParticle has already ::pt(), so the smeared value is in the decorator
286 // and can be accessed by
287 newpt = ptAcc(*newtrk);
288 }
289 if (newpt==0.) continue;
290 ATH_MSG_DEBUG ("Smeared Truth: "
291 <<" curv=" << 1./newpt
292 <<" phi=" << newtrk->phi()
293 <<" eta=" << newtrk->eta()
294 <<" d0=" << d0Acc(*newtrk)
295 <<" z0=" << z0Acc(*newtrk)
296 <<" pT=" << newpt
297 <<" PDGID=" << newtrk->pdgId()
298 <<" status=" << newtrk->status()
299 );
300 if (newtrk->parent(0)) ATH_MSG_DEBUG (" parent status=" << newtrk->parent(0)->pdgId());
301
302 if (m_enableMonitoring) {
303 hist("track_output_eta")->Fill(otrack.eta());
304 hist("track_output_theta")->Fill(otrack.theta());
305 hist("track_output_pt" )->Fill(part->charge()*otrack.pt() );
306 hist("track_output_phi")->Fill(otrack.phi());
307 hist("track_output_z0" )->Fill(otrack.z0() );
308 hist("track_output_d0" )->Fill(otrack.d0() );
309
310 hist("track_outputcoll_eta")->Fill(newtrk->eta());
311 hist("track_outputcoll_theta")->Fill(thetaAcc(*newtrk));
312 hist("track_outputcoll_pt" )->Fill(part->charge()* newpt/1000.);
313 hist("track_outputcoll_phi")->Fill(newtrk->phi());
314 hist("track_outputcoll_z0" )->Fill(z0Acc(*newtrk));
315 hist("track_outputcoll_d0" )->Fill(d0Acc(*newtrk));
316
317 hist("track_delta_eta")->Fill(newtrk->eta() - part->eta());
318 hist("track_delta_pt") ->Fill((newpt - part->pt())/1000.);
319 hist("track_delta_crv")->Fill(newtrk->charge()*1000./newpt - ((part->charge()*1000./part->pt())));
320 hist("track_delta_phi")->Fill(newtrk->phi() - part->phi());
321 hist("track_delta_z0" )->Fill(z0Acc(*newtrk) - z0Acc(*part));
322 hist("track_delta_d0" )->Fill(d0Acc(*newtrk) - d0Acc(*part));
323 }
324 } // end of loop
325 }
326 m_mySmearer->Clear(); // clear the smearer after each input track
327 }
328
329 ATH_MSG_DEBUG ("End of loop track #"<<n_input_tracks<<" ---> "<<" "<< n_output_tracks
330 <<" "<<n_output_narrow_tracks<<" "<<n_output_broad_tracks);
331 if (m_enableMonitoring) {
332 hist("n_input_tracks")->Fill(n_input_tracks);
333 hist("n_output_tracks")->Fill(n_output_tracks);
334 hist("n_output_narrow_tracks")->Fill(n_output_narrow_tracks);
335 hist("n_output_broad_tracks")->Fill(n_output_broad_tracks);
336 }
337 return StatusCode::SUCCESS;
338}
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
xAOD::TruthParticle_v1::status
int status() const
Status code.
xAOD::TruthParticle_v1::pdgId
int pdgId() const
PDG ID code.
xAOD::TruthParticle_v1::parent
const TruthParticle_v1 * parent(size_t i) const
Retrieve the i-th mother (TruthParticle) of this TruthParticle.
Definition TruthParticle_v1.cxx:126
xAOD::TruthParticle_v1::charge
double charge() const
Physical charge.
xAOD::TruthParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TruthParticle_v1.cxx:169
xAOD::TruthParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
Definition TruthParticle_v1.cxx:176
xAOD::TruthParticle
TruthParticle_v1 TruthParticle
Typedef to implementation.
Definition Event/xAOD/xAODTruth/xAODTruth/TruthParticle.h:15
xAOD::TruthParticleContainer
TruthParticleContainer_v1 TruthParticleContainer
Declare the latest version of the truth particle container.
Definition Event/xAOD/xAODTruth/xAODTruth/TruthParticleContainer.h:17

◆ sysInitialize()

StatusCode AthHistogramAlgorithm::sysInitialize ( )
virtualinherited

Initialization method invoked by the framework.

This method is responsible for any bookkeeping of initialization required by the framework itself. It will in turn invoke the initialize() method of the derived algorithm, and of any sub-algorithms which it creates.

Reimplemented from AthAlgorithm.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

76{
77 // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
78 // Bypass the initialization if the algorithm
79 // has already been initialized.
80 if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
81
82 // Set the Algorithm's properties
83 bindPropertiesTo( serviceLocator()->getOptsSvc() );
84
85 // Bypass the initialization if the algorithm is disabled.
86 // Need to do this after bindPropertiesTo.
87 if ( !isEnabled( ) ) return StatusCode::SUCCESS;
88
89 // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
90
91
92 // Get the THistSvc
93 ATH_CHECK ( histSvc().retrieve() );
94
95 // Configure the underlying AthHistogramming helper
96 ATH_CHECK ( AthHistogramming::configAthHistogramming( m_histSvc,
97 m_prefix, m_rootDir,
98 m_histNamePrefix, m_histNamePostfix,
99 m_histTitlePrefix, m_histTitlePostfix ) );
100
101 // Print some setup information into the log file
102 ATH_MSG_DEBUG ("Initializing " << name() << "...");
103 ATH_MSG_DEBUG (" using THistService = " << m_histSvc );
104 ATH_MSG_DEBUG (" using RootStreamName = " << m_prefix );
105 ATH_MSG_DEBUG (" using RootDirName = " << m_rootDir );
106 ATH_MSG_DEBUG (" using HistNamePrefix = " << m_histNamePrefix );
107 ATH_MSG_DEBUG (" using HistNamePostfix = " << m_histNamePostfix );
108 ATH_MSG_DEBUG (" using HistTitlePrefix = " << m_histTitlePrefix );
109 ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
110
111
112 // re-direct to base class...
113 return AthAlgorithm::sysInitialize();
114}
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition AthAlgorithm.cxx:66
AthHistogramAlgorithm::histSvc
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
Definition AthHistogramAlgorithm.h:113
AthHistogramAlgorithm::m_histNamePrefix
std::string m_histNamePrefix
The prefix for the histogram THx name.
Definition AthHistogramAlgorithm.h:94
AthHistogramAlgorithm::m_prefix
std::string m_prefix
Name of the ROOT output stream (file).
Definition AthHistogramAlgorithm.h:88
AthHistogramAlgorithm::m_histNamePostfix
std::string m_histNamePostfix
The postfix for the histogram THx name.
Definition AthHistogramAlgorithm.h:97
AthHistogramAlgorithm::m_histTitlePostfix
std::string m_histTitlePostfix
The postfix for the histogram THx title.
Definition AthHistogramAlgorithm.h:103
AthHistogramAlgorithm::m_rootDir
std::string m_rootDir
Name of the ROOT directory.
Definition AthHistogramAlgorithm.h:91
AthHistogramAlgorithm::m_histTitlePrefix
std::string m_histTitlePrefix
The prefix for the histogram THx title.
Definition AthHistogramAlgorithm.h:100
AthHistogramming::configAthHistogramming
StatusCode configAthHistogramming(const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
To be called by the derived classes to fill the internal configuration.
Definition AthHistogramming.cxx:51

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ tree()

TTree * AthHistogramming::tree ( const std::string & treeName,
const std::string & tDir = "",
const std::string & stream = "" )
protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 363 of file AthHistogramming.cxx.

364{
365 // Build a 32 bit hash out of the name
366 const hash_t treeHash = this->hash(treeName);
367
368 // See if this entry exists in the map
369 TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
370 if ( it == m_treeMap.end() ) // It doesn't exist!
371 { // Let's see into the THistSvc if somebody else has registered the TTree...
372
373 // Need to copy the strings as we will massage them from here on
374 std::string treeNameCopy = treeName;
375 std::string tDirCopy = tDir;
376 std::string streamCopy = stream;
377
378 // Massage the final string to book things
379 std::string bookingString("");
380 this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
381
382 TTree* treePointer(NULL);
383 if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
384 {
385 m_msg << MSG::WARNING
386 << "Problem retrieving the TTree with name " << treeNameCopy
387 << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
388 << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
389 return NULL;
390 }
391 // If we get to here, we actually found the TTree in the THistSvc.
392 // So let's add it to the local cache map and return its pointer
393 m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
394 return treePointer;
395 }
396
397 // Return the pointer to the TTree that we got from the local cache map
398 return it->second;
399}
AthenaInternal::getTree
std::pair< StatusCode, TTree * > getTree(ITHistSvc &svc, const std::string &name)
Definition AthenaPythonDict.h:67

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
312 std::vector<SG::VarHandleKey*> keys = a->keys();
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition AthCommonDataStore.h:398

Member Data Documentation

◆ m_d0DecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> EFTrackingSmearingAlg::m_d0DecoratorKey {this, "d0", "TruthParticles.d0", "Particle d0 decoration, set at initialisation"}
private

Definition at line 42 of file EFTrackingSmearingAlg.h.

42{this, "d0", "TruthParticles.d0", "Particle d0 decoration, set at initialisation"};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_effMap

EffMap_t AthHistogramming::m_effMap
privateinherited

The map of histogram names to their pointers.

Definition at line 210 of file AthHistogramming.h.

◆ m_EnableFakes

BooleanProperty EFTrackingSmearingAlg::m_EnableFakes {this,"IncludeDuplicatesAndFakes",false,"Enable fake track production"}
private

Definition at line 63 of file EFTrackingSmearingAlg.h.

63{this,"IncludeDuplicatesAndFakes",false,"Enable fake track production"};

◆ m_enableMonitoring

BooleanProperty EFTrackingSmearingAlg::m_enableMonitoring {this, "EnableMonitoring", false, "Enable debugging monitoring of the algorithm"}
private

Definition at line 55 of file EFTrackingSmearingAlg.h.

55{this, "EnableMonitoring", false, "Enable debugging monitoring of the algorithm"};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_FakeKillerEnable

BooleanProperty EFTrackingSmearingAlg::m_FakeKillerEnable {this,"FakeKillerEnable",false, "disable inclusion of broad fakes"}
private

Definition at line 64 of file EFTrackingSmearingAlg.h.

64{this,"FakeKillerEnable",false, "disable inclusion of broad fakes"};

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 224 of file AthHistogramming.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 203 of file AthHistogramming.h.

◆ m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix
privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

◆ m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix
privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc
privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

◆ m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix
privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

◆ m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix
privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

◆ m_IncludeFakesInResolutionCalculation

BooleanProperty EFTrackingSmearingAlg::m_IncludeFakesInResolutionCalculation
private
Initial value:
{this,"IncludeFakesInResolutionCalculation",false,
"modify calculation of covariance parameters to include weighted combination of core and broad RMS"}

Definition at line 65 of file EFTrackingSmearingAlg.h.

65 {this,"IncludeFakesInResolutionCalculation",false,
66 "modify calculation of covariance parameters to include weighted combination of core and broad RMS"};

◆ m_inputTrackParticleKey

SG::ReadHandleKey<xAOD::TrackParticleContainer> EFTrackingSmearingAlg::m_inputTrackParticleKey
private
Initial value:
{ this, "InputTrackParticleContainer", "InDetTrackParticles_tosmear",
"key for retrieval of input TrackParticles" }

Definition at line 29 of file EFTrackingSmearingAlg.h.

29 { this, "InputTrackParticleContainer", "InDetTrackParticles_tosmear",
30 "key for retrieval of input TrackParticles" };

◆ m_inputTracksPtCut

DoubleProperty EFTrackingSmearingAlg::m_inputTracksPtCut {this, "InputTracksPtCutGeV", 0.0, "Set input track pT cut"}
private

Definition at line 51 of file EFTrackingSmearingAlg.h.

51{this, "InputTracksPtCutGeV", 0.0, "Set input track pT cut"};

◆ m_inputTruthParticleKey

SG::ReadHandleKey<xAOD::TruthParticleContainer> EFTrackingSmearingAlg::m_inputTruthParticleKey
private
Initial value:
{this,"InputTruthParticleContainer","TruthParticles_tosmear",
"key for retrieval of input Truth particle"}

Definition at line 35 of file EFTrackingSmearingAlg.h.

35 {this,"InputTruthParticleContainer","TruthParticles_tosmear",
36 "key for retrieval of input Truth particle"};

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 250 of file AthHistogramming.h.

◆ m_mySmearer

std::unique_ptr<FakeTrackSmearer> EFTrackingSmearingAlg::m_mySmearer
private

Definition at line 72 of file EFTrackingSmearingAlg.h.

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 247 of file AthHistogramming.h.

◆ m_outputTrackParticleKey

SG::WriteHandleKey<xAOD::TrackParticleContainer> EFTrackingSmearingAlg::m_outputTrackParticleKey
private
Initial value:
{ this, "OutputTrackParticleContainer", "InDetTrackParticles_smeared",
"key for retrieval of output TrackParticles" }

Definition at line 32 of file EFTrackingSmearingAlg.h.

32 { this, "OutputTrackParticleContainer", "InDetTrackParticles_smeared",
33 "key for retrieval of output TrackParticles" };

◆ m_outputTracksPtCut

DoubleProperty EFTrackingSmearingAlg::m_outputTracksPtCut {this, "OutputTracksPtCutGeV", 1.0, "Set the output track pT cut"}
private

Definition at line 52 of file EFTrackingSmearingAlg.h.

52{this, "OutputTracksPtCutGeV", 1.0, "Set the output track pT cut"};

◆ m_outputTruthParticleKey

SG::WriteHandleKey<xAOD::TruthParticleContainer> EFTrackingSmearingAlg::m_outputTruthParticleKey
private
Initial value:
{this,"OutputTruthParticleContainer","TruthParticles_smeared",
"key for retrieval of output Truth particle"}

Definition at line 38 of file EFTrackingSmearingAlg.h.

38 {this,"OutputTruthParticleContainer","TruthParticles_smeared",
39 "key for retrieval of output Truth particle"};

◆ m_parameterizedTrackEfficiency

BooleanProperty EFTrackingSmearingAlg::m_parameterizedTrackEfficiency {this, "ParameterizedTrackEfficiency", false, "Enable parameterized efficiency for smearing"}
private

Definition at line 54 of file EFTrackingSmearingAlg.h.

54{this, "ParameterizedTrackEfficiency", false, "Enable parameterized efficiency for smearing"};

◆ m_parameterizedTrackEfficiency_LRT

BooleanProperty EFTrackingSmearingAlg::m_parameterizedTrackEfficiency_LRT {this, "ParameterizedTrackEfficiency_LRT", false, "Enable parameterized efficiency for LRT smearing"}
private

Definition at line 58 of file EFTrackingSmearingAlg.h.

58{this, "ParameterizedTrackEfficiency_LRT", false, "Enable parameterized efficiency for LRT smearing"};

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_ptDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> EFTrackingSmearingAlg::m_ptDecoratorKey {this, "pt", "TruthParticles.pt", "Particle pt decoration, set at initialisation"}
private

Definition at line 44 of file EFTrackingSmearingAlg.h.

44{this, "pt", "TruthParticles.pt", "Particle pt decoration, set at initialisation"};

◆ m_RandomSeed

LongLongProperty EFTrackingSmearingAlg::m_RandomSeed {this,"RandomSeed",0}
private

Definition at line 70 of file EFTrackingSmearingAlg.h.

70{this,"RandomSeed",0};

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_SetResolutionPtCutOff

DoubleProperty EFTrackingSmearingAlg::m_SetResolutionPtCutOff {this, "SetResolutionPtCutOff", 0., "Set ptCutoff off for semaring"}
private

Definition at line 50 of file EFTrackingSmearingAlg.h.

50{this, "SetResolutionPtCutOff", 0., "Set ptCutoff off for semaring"};

◆ m_SigmaScaleFactor

DoubleProperty EFTrackingSmearingAlg::m_SigmaScaleFactor {this, "SmearingScaleFactor", 1, "Set the smearing SF value for the sigmas"}
private

Definition at line 48 of file EFTrackingSmearingAlg.h.

48{this, "SmearingScaleFactor", 1, "Set the smearing SF value for the sigmas"};

◆ m_smearedTrackEfficiency

DoubleProperty EFTrackingSmearingAlg::m_smearedTrackEfficiency {this, "SmearedTrackEfficiency", 1.0, "Set track efficiency for smearing"}
private

Definition at line 53 of file EFTrackingSmearingAlg.h.

53{this, "SmearedTrackEfficiency", 1.0, "Set track efficiency for smearing"};

◆ m_smearedTrackEfficiency_d0high_LRT

DoubleProperty EFTrackingSmearingAlg::m_smearedTrackEfficiency_d0high_LRT {this, "SmearedTrackEfficiency_d0high_LRT", 1.0, "Set high d0 value for smeared track efficiency LRT"}
private

Definition at line 59 of file EFTrackingSmearingAlg.h.

59{this, "SmearedTrackEfficiency_d0high_LRT", 1.0, "Set high d0 value for smeared track efficiency LRT"};

◆ m_smearedTrackEfficiency_d0low_LRT

DoubleProperty EFTrackingSmearingAlg::m_smearedTrackEfficiency_d0low_LRT {this, "SmearedTrackEfficiency_d0low_LRT", 1.0, "Set low d0 value for smeared track efficiency LRT"}
private

Definition at line 60 of file EFTrackingSmearingAlg.h.

60{this, "SmearedTrackEfficiency_d0low_LRT", 1.0, "Set low d0 value for smeared track efficiency LRT"};

◆ m_smearTruthParticle

BooleanProperty EFTrackingSmearingAlg::m_smearTruthParticle {this, "SmearTruthParticle", false, "Enable smearing on truth particles, disabling the track smearing"}
private

Definition at line 56 of file EFTrackingSmearingAlg.h.

56{this, "SmearTruthParticle", false, "Enable smearing on truth particles, disabling the track smearing"};

◆ m_streamName

std::string AthHistogramming::m_streamName
privateinherited

Name of the ROOT output stream (file).

Definition at line 228 of file AthHistogramming.h.

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap
privateinherited

The map of TTree names to their pointers.

Definition at line 217 of file AthHistogramming.h.

◆ m_UseCoinToss

BooleanProperty EFTrackingSmearingAlg::m_UseCoinToss
private
Initial value:
{this,"UseCoinToss",false,
"if True, fakes are generated with flat inefficiency, else use poissonian criteria"}

Definition at line 67 of file EFTrackingSmearingAlg.h.

67 {this,"UseCoinToss",false,
68 "if True, fakes are generated with flat inefficiency, else use poissonian criteria"};

◆ m_UseResolutionPtCutOff

BooleanProperty EFTrackingSmearingAlg::m_UseResolutionPtCutOff {this, "UseResolutionPtCutOff", false, "Apply ptCutoff on smearing"}
private

Definition at line 49 of file EFTrackingSmearingAlg.h.

49{this, "UseResolutionPtCutOff", false, "Apply ptCutoff on smearing"};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_z0DecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> EFTrackingSmearingAlg::m_z0DecoratorKey {this, "z0", "TruthParticles.z0", "Particle z0 decoration, set at initialisation"}
private

Definition at line 43 of file EFTrackingSmearingAlg.h.

43{this, "z0", "TruthParticles.z0", "Particle z0 decoration, set at initialisation"};
The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.16.1