ATLAS Offline Software
Loading...
Searching...
No Matches
HI::HIPileupTool Class Reference

#include <HIPileupTool.h>

Inheritance diagram for HI::HIPileupTool:
Collaboration diagram for HI::HIPileupTool:

Public Member Functions

 HIPileupTool (const std::string &myname="HIPileupTool")
 ~HIPileupTool ()
virtual StatusCode initialize () override
 Dummy implementation of the initialisation function.
virtual void print () const override
 Print the state of the tool.
void write (TFile *fOut)
double get_efficiency (const xAOD::HIEventShapeContainer &, double) const
double get_purity (const xAOD::HIEventShapeContainer &, double &)
bool is_pileup (const xAOD::HIEventShapeContainer &, const xAOD::ZdcModuleContainer &) const override
virtual const asg::AcceptInfogetAcceptInfo (const xAOD::HIEventShapeContainer &evShCont, const xAOD::ZdcModuleContainer &ZdcCont) const override
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 sysInitialize () override
 Perform system initialization for an algorithm.
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
Additional helper functions, not directly mimicking Athena
template<class T>
const T * getProperty (const std::string &name) const
 Get one of the tool's properties.
const std::string & msg_level_name () const __attribute__((deprecated))
 A deprecated function for getting the message level's name.
const std::string & getName (const void *ptr) const
 Get the name of an object that is / should be in the event store.
SG::sgkey_t getKey (const void *ptr) const
 Get the (hashed) key of an object that is in the event store.

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.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

double get_et (const xAOD::HIEventShapeContainer &evShCont) const
double get_nNeutrons (const xAOD::ZdcModuleContainer &ZdcCont) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

asg::AcceptInfo m_accept
std::vector< double > m_etCutAndPurity
std::string m_fname
double m_etThreshold = 0.0
double m_purityCut = 0.0
bool m_setup
const double m_npeak = 2510.
TH2D * m_hEvents = nullptr
TH2D * m_hPileUp = nullptr
TH1D * m_hCut = nullptr
TH1D * m_hEff = nullptr
TH1D * m_hPurity = nullptr
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

Detailed Description

Definition at line 26 of file HIPileupTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ HIPileupTool()

HI::HIPileupTool::HIPileupTool ( const std::string & myname = "HIPileupTool")

Definition at line 28 of file HIPileupTool.cxx.

28 : AsgTool(myname),
29 m_accept("HIPileup")
30 {
31 declareProperty("EtCutAndPurity", m_etCutAndPurity, "Et and purity for efficiency calculation");
32 declareProperty("inpFilePath", m_fname, "Calibration file path");
33
34 m_setup = 0;
35}
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
asg::AcceptInfo m_accept
std::vector< double > m_etCutAndPurity
std::string m_fname
AsgTool(const std::string &name)
Constructor specifying the tool instance's name.
Definition AsgTool.cxx:58

◆ ~HIPileupTool()

HI::HIPileupTool::~HIPileupTool ( )

Definition at line 37 of file HIPileupTool.cxx.

37 {
38
39 if (m_setup) {
40 delete m_hEvents; delete m_hPileUp;
41 delete m_hCut; delete m_hEff; delete m_hPurity;
42 }
43}

Member Function Documentation

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::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 {
160 hndl.value(),
161 hndl.documentation());
162
163 }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::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;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::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

◆ get_efficiency()

double HI::HIPileupTool::get_efficiency ( const xAOD::HIEventShapeContainer & evShCont,
double FCal_Et ) const

Definition at line 141 of file HIPileupTool.cxx.

141 {
142
143 FCal_Et = get_et(evShCont);
144 return m_hEff->GetBinContent(m_hEff->FindFixBin(FCal_Et));
145}
double get_et(const xAOD::HIEventShapeContainer &evShCont) const

◆ get_et()

double HI::HIPileupTool::get_et ( const xAOD::HIEventShapeContainer & evShCont) const
private

Definition at line 154 of file HIPileupTool.cxx.

154 {
155 double Fcal_Et=0;
156 Fcal_Et = evShCont.at(5)->et()*1e-6;
157 return Fcal_Et;
158}
const T * at(size_type n) const
Access an element, as an rvalue.

◆ get_nNeutrons()

double HI::HIPileupTool::get_nNeutrons ( const xAOD::ZdcModuleContainer & ZdcCont) const
private

Definition at line 160 of file HIPileupTool.cxx.

160 {
161
162 bool isCalib=1;
163 double ZdcE=0;
164 for (const auto *zdcModule : ZdcCont) {
165 if (zdcModule->zdcType()!=0) continue;
166
167 static const SG::ConstAccessor<float> CalibEnergyAcc("CalibEnergy");
168 if (!(CalibEnergyAcc.isAvailable(*zdcModule))) {isCalib = 0; continue;}
169 float modE = CalibEnergyAcc(*zdcModule);
170 ZdcE += modE;
171 }
172 if (!isCalib) throw std::invalid_argument("ZDC Module not Calibrated");
173
174 return ZdcE/m_npeak;
175}
bool isCalib()
const double m_npeak

◆ get_purity()

double HI::HIPileupTool::get_purity ( const xAOD::HIEventShapeContainer & evShCont,
double & FCal_Et )

Definition at line 147 of file HIPileupTool.cxx.

147 {
148
149 FCal_Et = get_et(evShCont);
150 return m_hPurity->GetBinContent(m_hPurity->FindFixBin(FCal_Et));
151}

◆ getAcceptInfo()

const asg::AcceptInfo & HI::HIPileupTool::getAcceptInfo ( const xAOD::HIEventShapeContainer & evShCont,
const xAOD::ZdcModuleContainer & ZdcCont ) const
overridevirtual

Implements HI::IHIPileupTool.

Definition at line 199 of file HIPileupTool.cxx.

200{
201 asg::AcceptData(&m_accept).clear();
202 bool pileup_dec = !is_pileup(evShCont, ZdcCont);
203 asg::AcceptData(&m_accept).setCutResult("pileupVeto", pileup_dec);
204 return m_accept;
205}
bool is_pileup(const xAOD::HIEventShapeContainer &, const xAOD::ZdcModuleContainer &) const override

◆ getKey()

SG::sgkey_t asg::AsgTool::getKey ( const void * ptr) const
inherited

Get the (hashed) key of an object that is in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the SG::sgkey_t key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getName
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The hashed key of the object in the store. If not found, an invalid (zero) key.

Definition at line 119 of file AsgTool.cxx.

119 {
120
121#ifdef XAOD_STANDALONE
122 // In case we use @c xAOD::TEvent, we have a direct function call
123 // for this.
124 return evtStore()->event()->getKey( ptr );
125#else
126 const SG::DataProxy* proxy = evtStore()->proxy( ptr );
127 return ( proxy == nullptr ? 0 : proxy->sgkey() );
128#endif // XAOD_STANDALONE
129 }
ServiceHandle< StoreGateSvc > & evtStore()

◆ getName()

const std::string & asg::AsgTool::getName ( const void * ptr) const
inherited

Get the name of an object that is / should be in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the std::string name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getKey
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The string name of the object in the store. If not found, an empty string.

Definition at line 106 of file AsgTool.cxx.

106 {
107
108#ifdef XAOD_STANDALONE
109 // In case we use @c xAOD::TEvent, we have a direct function call
110 // for this.
111 return evtStore()->event()->getName( ptr );
112#else
113 const SG::DataProxy* proxy = evtStore()->proxy( ptr );
114 static const std::string dummy = "";
115 return ( proxy == nullptr ? dummy : proxy->name() );
116#endif // XAOD_STANDALONE
117 }

◆ getProperty()

template<class T>
const T * asg::AsgTool::getProperty ( const std::string & name) const
inherited

Get one of the tool's properties.

◆ initialize()

StatusCode HI::HIPileupTool::initialize ( void )
overridevirtual

Dummy implementation of the initialisation function.

It's here to allow the dual-use tools to skip defining an initialisation function. Since many are doing so...

Reimplemented from asg::AsgTool.

Definition at line 45 of file HIPileupTool.cxx.

45 {
46
47 ATH_MSG_VERBOSE( "Initialising tool " << name() );
48
49 m_accept.addCut("pileupVeto","for pileup rejection");
50
51 const char* fullInpFilePath;
52 if (!m_fname.empty()) fullInpFilePath = gSystem->ExpandPathName (m_fname.c_str());
53 else fullInpFilePath = gSystem->ExpandPathName ("$ROOTCOREBIN/data/HIEventUtils/HIRun2PileUp_PbPb5p02_v2.root");
54 if (gSystem->AccessPathName(fullInpFilePath)) throw std::invalid_argument(("File not found at \"" + (std::string)fullInpFilePath + "\"").c_str());
55
56 ATH_MSG_INFO(("Read pileup cuts from " + (std::string)fullInpFilePath).c_str());
57 TFile* fIn = new TFile(fullInpFilePath);
58
59 std::string hname;
60 hname = "hEvents"; if (!(fIn->GetListOfKeys()->Contains(hname.c_str()))) {
61 fIn->Close(); throw std::invalid_argument(("Expected histogram \"" + hname + "\" not found in \"" + m_fname + "\"").c_str());}
62 hname = "hPileUp"; if (!(fIn->GetListOfKeys()->Contains(hname.c_str()))) {
63 fIn->Close(); throw std::invalid_argument(("Expected histogram \"" + hname + "\" not found in \"" + m_fname + "\"").c_str());}
64 hname = "hCut"; if (!(fIn->GetListOfKeys()->Contains(hname.c_str()))) {
65 fIn->Close(); throw std::invalid_argument(("Expected histogram \"" + hname + "\" not found in \"" + m_fname + "\"").c_str());}
66 hname = "hEff"; if (!(fIn->GetListOfKeys()->Contains(hname.c_str()))) {
67 fIn->Close(); throw std::invalid_argument(("Expected histogram \"" + hname + "\" not found in \"" + m_fname + "\"").c_str());}
68 hname = "hPurity"; if (!(fIn->GetListOfKeys()->Contains(hname.c_str()))) {
69 fIn->Close(); throw std::invalid_argument(("Expected histogram \"" + hname + "\" not found in \"" + m_fname + "\"").c_str());}
70
71 m_hEvents = (TH2D*)((TH2D*)fIn->Get("hEvents"))->Clone("hEvents_HIPileupTool"); m_hEvents->SetDirectory(0);
72 m_hPileUp = (TH2D*)((TH2D*)fIn->Get("hPileUp"))->Clone("hPileUp_HIPileupTool"); m_hPileUp->SetDirectory(0);
73 m_hCut = (TH1D*)((TH1D*)fIn->Get("hCut"))->Clone("hCut_HIPileupTool"); m_hCut->SetDirectory(0);
74 m_hEff = (TH1D*)((TH1D*)fIn->Get("hEff"))->Clone("hEff_HIPileupTool"); m_hEff->SetDirectory(0);
75 m_hPurity = (TH1D*)((TH1D*)fIn->Get("hPurity"))->Clone("hPurity_HIPileupTool"); m_hPurity->SetDirectory(0);
76 m_setup = 1;
77
78 fIn->Close();
79
80 m_etThreshold = 6.0;
81 m_purityCut = -1.0;
82 if (!m_etCutAndPurity.empty()) {
83 if (m_etCutAndPurity.size()<2) ATH_MSG_WARNING("Cuts vector should have two elements. Will use default settings");
84 else {
85 if (m_etCutAndPurity.at(1)<0 || m_etCutAndPurity.at(1)>1) ATH_MSG_WARNING("Purity cut should be between 0 and 1. Will use default settings");
86 else {
89 if (fabs(m_purityCut-1)<1e-6) m_purityCut -= 1e-3;
90 }
91 }
92 }
93
94 if (m_purityCut<0) {
95 ATH_MSG_INFO("ET threshold for purity cut = NONE");
96 ATH_MSG_INFO("Purity cut = DEFAULT");
97 }
98 else {
99 ATH_MSG_INFO("ET threshold for purity cut = " << m_etThreshold << " TeV");
100 ATH_MSG_INFO("Purity cut = " << m_purityCut );
101 }
102
103 for (int ibx=0;ibx<m_hPileUp->GetNbinsX();ibx++) {
104 if (m_hPileUp->GetXaxis()->GetBinCenter(ibx+1) < m_etThreshold) continue;
105 int NY = m_hPileUp->GetNbinsY();
106 if (m_hPileUp->Integral(ibx+1,ibx+1,1,NY) == 0) {m_hCut->SetBinContent(ibx+1,0); continue;}
107 //Cut on 5.0 because efficiency is calculated only for FCal < 5 TeV
108 if (m_hPileUp->GetXaxis()->GetBinCenter(ibx+1)>5.0 && fabs(m_purityCut-1)<1e-6) {m_hCut->SetBinContent(ibx+1,0); continue;}
109 for (int iby=0;iby<NY;iby++) {
110 if (m_hPileUp->Integral(ibx+1,ibx+1,NY-iby,NY)/m_hPileUp->Integral(ibx+1,ibx+1,1,NY) < m_purityCut) continue;
111 m_hCut->SetBinContent(ibx+1,m_hPileUp->GetYaxis()->GetBinCenter(NY-iby)); break;
112 }
113 }
114
115 for (int ibx=0;ibx<m_hEvents->GetNbinsX();ibx++) {
116 if (m_hEvents->GetXaxis()->GetBinCenter(ibx+1) < m_etThreshold) {m_hEff->SetBinContent(ibx+1,1-m_hEff->GetBinContent(ibx+1)); continue;}
117 double cut = m_hEvents->Integral(ibx+1,ibx+1,m_hEvents->GetYaxis()->FindBin(m_hCut->GetBinContent(ibx+1)), m_hEvents->GetNbinsY());
118 double all = m_hEvents->Integral(ibx+1,ibx+1,1,m_hEvents->GetYaxis()->GetNbins()) - m_hPileUp->Integral(ibx+1,ibx+1,1,m_hPileUp->GetYaxis()->GetNbins());
119 double pu_cut = m_hPileUp->Integral(ibx+1,ibx+1,m_hPileUp->GetYaxis()->FindBin(m_hCut->GetBinContent(ibx+1)), m_hPileUp->GetNbinsY());
120 double pu_all = m_hPileUp->Integral(ibx+1,ibx+1,1,m_hPileUp->GetYaxis()->GetNbins());
121 double eff = 1 - (cut - pu_cut)/all; if ((cut - pu_cut) < 0) eff=1;
122 double pty = pu_cut/pu_all;
123 if (all!=0) m_hEff->SetBinContent(ibx+1,eff);
124 if (pu_all!=0) m_hPurity->SetBinContent(ibx+1,pty);
125 }
126
127 return StatusCode::SUCCESS;
128}
#define ATH_MSG_INFO(x)
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
cut
This script demonstrates how to call a C++ class from Python Also how to use PyROOT is shown.

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

◆ is_pileup()

bool HI::HIPileupTool::is_pileup ( const xAOD::HIEventShapeContainer & evShCont,
const xAOD::ZdcModuleContainer & ZdcCont ) const
overridevirtual

Implements HI::IHIPileupTool.

Definition at line 130 of file HIPileupTool.cxx.

130 {
131
132 bool kPileup = false;
133 double FCal_Et = get_et(evShCont);
134 double nNeutrons;
135 nNeutrons = get_nNeutrons(ZdcCont);
136 if (nNeutrons > m_hCut->GetBinContent(m_hCut->FindFixBin(FCal_Et))) kPileup = true;
137
138 return kPileup;
139}
double get_nNeutrons(const xAOD::ZdcModuleContainer &ZdcCont) const

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg_level_name()

const std::string & asg::AsgTool::msg_level_name ( ) const
inherited

A deprecated function for getting the message level's name.

Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:

MSG::name( msg().level() )

This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.

Returns
The string name of the current minimum message level that's printed

Definition at line 101 of file AsgTool.cxx.

101 {
102
103 return MSG::name( msg().level() );
104 }
MsgStream & msg() const
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition MsgLevel.cxx:19

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

◆ print()

void HI::HIPileupTool::print ( ) const
overridevirtual

Print the state of the tool.

Reimplemented from asg::AsgTool.

Definition at line 177 of file HIPileupTool.cxx.

177 {
178 AsgTool::print();
179 if (m_purityCut<0) {
180 ATH_MSG_INFO("ET threshold for purity cut = NONE");
181 ATH_MSG_INFO("Purity cut = DEFAULT");
182 }
183 else {
184 ATH_MSG_INFO("ET threshold for purity cut = " << m_etThreshold << " TeV");
185 ATH_MSG_INFO("Purity cut = " << m_purityCut );
186 }
187}

◆ 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< AlgTool > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::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) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

◆ write()

void HI::HIPileupTool::write ( TFile * fOut)

Definition at line 189 of file HIPileupTool.cxx.

189 {
190
191 fOut->cd(0);
192 m_hEvents->Write();
193 m_hPileUp->Write();
194 m_hCut->Write();
195 m_hEff->Write();
196 m_hPurity->Write();
197}

Member Data Documentation

◆ m_accept

asg::AcceptInfo HI::HIPileupTool::m_accept
private

Definition at line 52 of file HIPileupTool.h.

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_etCutAndPurity

std::vector<double> HI::HIPileupTool::m_etCutAndPurity
private

Definition at line 54 of file HIPileupTool.h.

◆ m_etThreshold

double HI::HIPileupTool::m_etThreshold = 0.0
private

Definition at line 57 of file HIPileupTool.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_fname

std::string HI::HIPileupTool::m_fname
private

Definition at line 55 of file HIPileupTool.h.

◆ m_hCut

TH1D* HI::HIPileupTool::m_hCut = nullptr
private

Definition at line 65 of file HIPileupTool.h.

◆ m_hEff

TH1D* HI::HIPileupTool::m_hEff = nullptr
private

Definition at line 66 of file HIPileupTool.h.

◆ m_hEvents

TH2D* HI::HIPileupTool::m_hEvents = nullptr
private

Definition at line 63 of file HIPileupTool.h.

◆ m_hPileUp

TH2D* HI::HIPileupTool::m_hPileUp = nullptr
private

Definition at line 64 of file HIPileupTool.h.

◆ m_hPurity

TH1D* HI::HIPileupTool::m_hPurity = nullptr
private

Definition at line 67 of file HIPileupTool.h.

◆ m_npeak

const double HI::HIPileupTool::m_npeak = 2510.
private

Definition at line 61 of file HIPileupTool.h.

◆ m_purityCut

double HI::HIPileupTool::m_purityCut = 0.0
private

Definition at line 58 of file HIPileupTool.h.

◆ m_setup

bool HI::HIPileupTool::m_setup
private

Definition at line 59 of file HIPileupTool.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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