|
ATLAS Offline Software
|
#include <VoronoiWeightTool.h>
|
| VoronoiWeightTool (const std::string &name) |
|
| ~VoronoiWeightTool () |
|
StatusCode | initialize () |
| Dummy implementation of the initialisation function. More...
|
|
StatusCode | process (xAOD::IParticleContainer *cont) const override final |
|
virtual void | print () const |
| Print the state of the tool. More...
|
|
ServiceHandle< StoreGateSvc > & | evtStore () |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | evtStore () const |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | detStore () const |
| The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
virtual StatusCode | sysInitialize () override |
| Perform system initialization for an algorithm. More...
|
|
virtual StatusCode | sysStart () override |
| Handle START transition. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | inputHandles () const override |
| Return this algorithm's input handles. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | outputHandles () const override |
| Return this algorithm's output handles. More...
|
|
Gaudi::Details::PropertyBase & | declareProperty (Gaudi::Property< T > &t) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc="none") |
| Declare a new Gaudi property. More...
|
|
void | updateVHKA (Gaudi::Details::PropertyBase &) |
|
MsgStream & | msg () const |
|
MsgStream & | msg (const MSG::Level lvl) const |
|
bool | msgLvl (const MSG::Level lvl) const |
|
|
StatusCode | process_impl (xAOD::IParticleContainer *cont) const |
|
StatusCode | makeVoronoiParticles (std::vector< fastjet::PseudoJet > &particles, std::vector< std::pair< fastjet::PseudoJet, std::vector< float > > > &) const |
|
void | spreadPt (std::vector< std::pair< fastjet::PseudoJet, std::vector< float > > > &correctedptvec, float spreadr=0.4, float alpha=2) const |
|
StatusCode | setEtaPhi (xAOD::IParticle *obj, float eta, float phi) const |
|
StatusCode | setEnergyPt (xAOD::IParticle *obj, float e, float pt, const SG::AuxElement::Accessor< float > *weightAcc=nullptr) const |
|
StatusCode | setP4 (xAOD::IParticle *obj, const xAOD::JetFourMom_t &p4, const SG::AuxElement::Accessor< float > *weightAcc=nullptr) const |
|
void | renounceArray (SG::VarHandleKeyArray &handlesArray) |
| remove all handles from I/O resolution More...
|
|
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. More...
|
|
- Author
- Aviv Cukierman
- Date
- May 2016
This is a tool to reweight clusters by subtracting pileup based on Voronoi areas It creates a new cluster collection with reweighted clusters Documentation: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/VoronoiWeightTool
Definition at line 26 of file VoronoiWeightTool.h.
◆ StoreGateSvc_t
◆ VoronoiWeightTool()
VoronoiWeightTool::VoronoiWeightTool |
( |
const std::string & |
name | ) |
|
◆ ~VoronoiWeightTool()
VoronoiWeightTool::~VoronoiWeightTool |
( |
| ) |
|
|
default |
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
253 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleKey
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 221 of file AthCommonDataStore.h.
229 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This is the generic version, for types that do not derive from SG::VarHandleKey
. It just forwards to the base class version of declareProperty
.
Definition at line 333 of file AthCommonDataStore.h.
338 return PBASE::declareProperty(
name, property,
doc);
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This dispatches to either the generic declareProperty
or the one for VarHandle/Key/KeyArray.
Definition at line 352 of file AthCommonDataStore.h.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ extraDeps_update_handler()
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
◆ getKey()
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
-
ptr | The 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.
121 #ifdef XAOD_STANDALONE
127 return (
proxy ==
nullptr ? 0 :
proxy->sgkey() );
128 #endif // XAOD_STANDALONE
◆ 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
-
ptr | The 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.
108 #ifdef XAOD_STANDALONE
114 static const std::string
dummy =
"";
116 #endif // XAOD_STANDALONE
◆ getProperty()
template<class T >
const T* asg::AsgTool::getProperty |
( |
const std::string & |
name | ) |
const |
|
inherited |
Get one of the tool's properties.
◆ initialize()
StatusCode VoronoiWeightTool::initialize |
( |
| ) |
|
|
virtual |
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 89 of file VoronoiWeightTool.cxx.
93 ATH_MSG_ERROR(
"Incompatible configuration: setting both IgnoreChargedPFO and ApplyToChargedPFO to true"
94 <<
"will set all cPFOs to zero");
95 return StatusCode::FAILURE;
98 ATH_MSG_ERROR(
"Incompatible configuration: ApplyToNeutralPFO=False -- what kind of pileup do you wish to suppress?");
99 return StatusCode::FAILURE;
102 return StatusCode::SUCCESS;
◆ inputHandles()
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.
◆ makeVoronoiParticles()
StatusCode VoronoiWeightTool::makeVoronoiParticles |
( |
std::vector< fastjet::PseudoJet > & |
particles, |
|
|
std::vector< std::pair< fastjet::PseudoJet, std::vector< float > > > & |
correctedptvec |
|
) |
| const |
|
protected |
Definition at line 213 of file VoronoiWeightTool.cxx.
214 std::vector<fastjet::PseudoJet> & inputConst =
particles;
215 fastjet::Selector jselector = fastjet::SelectorAbsRapRange(0.0,2.1);
218 fastjet::JetDefinition
jetDef(algo, jetR,fastjet::E_scheme, fastjet::Best);
219 fastjet::AreaDefinition area_def(fastjet::voronoi_area, fastjet::VoronoiAreaSpec(0.9));
221 fastjet::ClusterSequenceArea clustSeq(inputConst,
jetDef, area_def);
222 fastjet::JetMedianBackgroundEstimator bge(jselector,
jetDef,area_def);
224 bge.set_particles(inputConst);
225 std::vector<fastjet::PseudoJet> inclusiveJets = sorted_by_pt(clustSeq.inclusive_jets(0));
228 float rho = bge.rho();
229 float sigma = bge.sigma();
230 for(
unsigned int iJet = 0 ; iJet < inclusiveJets.size() ; iJet++){
231 fastjet::PseudoJet
jet = inclusiveJets[iJet];
232 std::vector<fastjet::PseudoJet> constituents =
jet.constituents();
233 for(
const auto& cons : constituents){
234 float pt = cons.pt();
238 float voro0pt = subPt * (subPt > 0);
240 std::vector<float> algopts;
241 algopts.push_back(subPt);
242 algopts.push_back(voro0pt);
243 algopts.push_back(voro1pt);
244 algopts.push_back(0);
245 std::pair <fastjet::PseudoJet,std::vector<float> > pjcptpair (cons,algopts);
246 correctedptvec.push_back(pjcptpair);
252 return StatusCode::SUCCESS;
◆ msg() [1/2]
◆ msg() [2/2]
◆ 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.
◆ msgLvl()
◆ outputHandles()
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 asg::AsgTool::print |
( |
| ) |
const |
|
virtualinherited |
Print the state of the tool.
Implements asg::IAsgTool.
Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ process()
|
finaloverridevirtualinherited |
◆ process_impl()
Reimplemented from JetConstituentModifierBase.
Definition at line 108 of file VoronoiWeightTool.cxx.
141 std::vector< std::pair< fastjet::PseudoJet, std::vector<float> > > ptvec;
180 bool endContainer =
part->e()<FLT_MIN;
181 bool endVec =
i>=ptvec.size();
182 if(endVec && endContainer){
185 else if(endContainer || endVec){
186 ATH_MSG_ERROR(
"Filtered particle list doesn't have same number of elements as the list returned by FastJet.");
187 return StatusCode::FAILURE;
191 float Containerpt =
part->pt();
192 float PJpt = ptvec[
i].first.pt();
193 ATH_MSG_VERBOSE(
"Container pt: " << Containerpt <<
", fastjet initial pt: " << PJpt
194 <<
", subtracted pt: " << ptvec[
i].
second[
alg]);
195 if (fabs(Containerpt-PJpt) > 0.1){
198 return StatusCode::FAILURE;
200 newPt = ptvec[
i].second[
alg];
207 float w =
part->pt()>FLT_MIN ? newPt/
part->pt() : 0.;
210 return StatusCode::SUCCESS;
◆ renounce()
◆ renounceArray()
◆ setEnergyPt()
Definition at line 61 of file JetConstituentModifierBase.cxx.
70 if(weightAcc) (*weightAcc)(*clus) = clus->
calE() > 0. ?
e / clus->calE() : 0.;
79 if(weightAcc) (*weightAcc)(*pfo) = pfo->
pt() > 0. ?
pt / pfo->pt() : 0.;
82 pfo->setP4(
pt, pfo->eta(), pfo->phi());
91 if(weightAcc) (*weightAcc)(*pfo) = pfo->
pt() > 0. ?
pt / pfo->pt() : 0.;
92 pfo->setP4(
pt, pfo->eta(), pfo->phi(), 0.);
99 if( tcc->
taste() != 0) {
100 if(weightAcc) (*weightAcc)(*tcc) = tcc->
pt() > 0. ?
pt / tcc->pt() : 0.;
101 tcc->setParameters(
pt, tcc->eta(), tcc->phi(), tcc->m(),
xAOD::TrackCaloCluster::Taste(tcc->taste()), tcc->trackParticleLink(), tcc->caloClusterLinks());
108 return StatusCode::FAILURE;
110 return StatusCode::SUCCESS;
◆ setEtaPhi()
StatusCode JetConstituentModifierBase::setEtaPhi |
( |
xAOD::IParticle * |
obj, |
|
|
float |
eta, |
|
|
float |
phi |
|
) |
| const |
|
protectedinherited |
◆ setP4()
Definition at line 113 of file JetConstituentModifierBase.cxx.
120 if(weightAcc) (*weightAcc)(*clus) = clus->
calE() > 0. ? p4.e() / clus->calE() : 0.;
121 clus->setCalE(p4.e());
122 clus->setCalEta(p4.eta());
123 clus->setCalPhi(p4.phi());
132 if(weightAcc) (*weightAcc)(*pfo) = pfo->
pt() > 0. ? p4.pt() / pfo->pt() : 0.;
133 pfo->setP4(p4.pt(),p4.eta(),p4.phi(),p4.mass());
142 if(weightAcc) (*weightAcc)(*pfo) = pfo->
pt() > 0. ? p4.pt() / pfo->pt() : 0.;
143 pfo->setP4(p4.pt(),p4.eta(),p4.phi(),p4.mass());
150 if( tcc->
taste() != 0) {
151 if(weightAcc) (*weightAcc)(*tcc) = tcc->
pt() > 0. ? p4.pt() / tcc->pt() : 0.;
152 tcc->setParameters(p4.pt(), p4.eta(), p4.phi(), p4.mass(),
xAOD::TrackCaloCluster::Taste(tcc->taste()), tcc->trackParticleLink(), tcc->caloClusterLinks());
160 return StatusCode::FAILURE;
162 return StatusCode::SUCCESS;
◆ spreadPt()
void VoronoiWeightTool::spreadPt |
( |
std::vector< std::pair< fastjet::PseudoJet, std::vector< float > > > & |
correctedptvec, |
|
|
float |
spreadr = 0.4 , |
|
|
float |
alpha = 2 |
|
) |
| const |
|
protected |
Definition at line 255 of file VoronoiWeightTool.cxx.
258 size_t Nparticles = correctedptvec.size();
259 std::vector<float> spreadPT(Nparticles);
260 std::vector<bool> isPositive(Nparticles);
261 for(
size_t iPart = 0; iPart < Nparticles; iPart++){
262 spreadPT[iPart] = correctedptvec[iPart].second[0];
263 isPositive[iPart] = spreadPT[iPart]>0;
266 std::vector<std::vector<std::pair<size_t,float> > > particle_drs;
267 for(
size_t iPart = 0; iPart < Nparticles; iPart++){
268 fastjet::PseudoJet iparticle = correctedptvec[iPart].first;
269 std::vector<std::pair<size_t,float> > this_particle_drs;
270 for(
size_t jPart = 0; jPart < Nparticles; jPart++){
271 if(iPart == jPart)
continue;
272 if(!isPositive[jPart])
continue;
273 fastjet::PseudoJet jparticle = correctedptvec[jPart].first;
274 float dphi = iparticle.delta_phi_to(jparticle);
275 float deta = iparticle.eta() - jparticle.eta();
276 float dr2 =
pow(dphi,2) +
pow(deta,2);
277 if(dr2 >
pow(spreadr,2))
continue;
278 std::pair<size_t,float> jdrpair (jPart,dr2);
279 this_particle_drs.push_back(jdrpair);
281 particle_drs.push_back(this_particle_drs);
284 for(
size_t i = 0;
i < Nparticles;
i++){
285 if(!(spreadPT[
i]<0))
continue;
289 for(
size_t j=0; j<particle_drs[
i].size(); j++){
290 float dr = particle_drs[
i][j].second;
291 if(
dr>FLT_MIN) sumdR2 += 1./(
pow(
dr,alpha/2));
294 if(sumdR2 > FLT_MIN){
295 float spreadPT_orig = spreadPT[
i];
296 for(
size_t j=0; j<particle_drs[
i].size(); j++){
297 float dr = particle_drs[
i][j].second;
298 float realid = particle_drs[
i][j].first;
301 if(fabs(
weight*spreadPT_orig)>spreadPT[realid]){
302 spreadPT[
i]+=spreadPT[realid];
306 spreadPT[realid]+=
weight*spreadPT_orig;
307 spreadPT[
i]-=
weight*spreadPT_orig;
314 for(
size_t iPart = 0; iPart < Nparticles; iPart++){
315 correctedptvec[iPart].second[3] = spreadPT[iPart] * (spreadPT[iPart] > 0);
◆ sysInitialize()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_applyToChargedPFO
bool JetConstituentModifierBase::m_applyToChargedPFO =true |
|
protectedinherited |
◆ m_applyToNeutralPFO
bool JetConstituentModifierBase::m_applyToNeutralPFO =true |
|
protectedinherited |
◆ m_detStore
◆ m_doSpread
bool VoronoiWeightTool::m_doSpread = true |
|
protected |
◆ m_evtStore
◆ m_ignoreChargedPFOs
bool VoronoiWeightTool::m_ignoreChargedPFOs =true |
|
protected |
◆ m_inputType
unsigned int JetConstituentModifierBase::m_inputType =xAOD::Type::Other |
|
protectedinherited |
◆ m_maxArea
float VoronoiWeightTool::m_maxArea = 1e9 |
|
protected |
◆ m_nSigma
int VoronoiWeightTool::m_nSigma = 0 |
|
protected |
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
virtual double pt() const
The transverse momentum ( ) of the particle.
Scalar phi() const
phi method
JetConstituentModifierBase(const std::string &name)
Helper class to provide type-safe access to aux data.
Scalar eta() const
pseudorapidity method
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
flt_t calE() const
Geet Energy in signal state CALIBRATED.
Class describing a TrackCaloCluster.
StatusCode accept(const xAOD::Muon *mu)
virtual StatusCode process_impl(xAOD::IParticleContainer *) const
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
void setCalEta(flt_t)
Set for signal state CALIBRATED.
virtual double pt() const override
#define ATH_MSG_VERBOSE(x)
Class providing the definition of the 4-vector interface.
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
virtual void setOwner(IDataHandleHolder *o)=0
TrackCaloCluster_v1 TrackCaloCluster
Reference the current persistent version:
virtual double pt() const
The transverse momentum ( ) of the particle.
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Description of a calorimeter cluster.
void setCalPhi(flt_t)
Set for signal state CALIBRATED.
FlowElement_v1 FlowElement
Definition of the current "pfo version".
float charge() const
get charge of PFO
signal_t signalType() const
const T * front() const
Access the first element in the collection as an rvalue.
@ ParticleFlow
The object is a particle-flow object.
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
bool isCharged() const
is a charged PFO
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Class describing a particle flow object.
StatusCode setEnergyPt(xAOD::IParticle *obj, float e, float pt, const SG::AuxElement::Accessor< float > *weightAcc=nullptr) const
virtual int taste() const
The taste of the particle.
T sort_container_pt(T *inCont)
void sort(typename std::reverse_iterator< DataModel_detail::iterator< DVL > > beg, typename std::reverse_iterator< DataModel_detail::iterator< DVL > > end, const Compare &comp)
Specialization of sort for DataVector/List.
Taste
Type of TrackCaloCluster - Charged - Neutral - Combined.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
constexpr int pow(int base, int exp) noexcept
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
bool empty() const noexcept
Returns true if the collection is empty.
A detector object made of other lower level object(s)