ATLAS Offline Software
PFO_v1.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 // EDM include(s):
6 #include "xAODCore/AuxStoreAccessorMacros.h"
7 #include "AthLinks/ElementLink.h"
8 
9 
10 // Local include(s):
11 #include "xAODPFlow/versions/PFO_v1.h"
12 #include "xAODPFlow/versions/PFOAttributesAccessor_v1.h"
13 
14 namespace xAOD {
15 
16  struct PFOParticleTypeMapper_temp {
17  static const SG::AuxElement::Accessor<std::vector<ElementLink<IParticleContainer > > >* getAccessor(PFODetails::PFOParticleType type) {
18  switch (type){
19  case PFODetails::CaloCluster:
20  const static SG::AuxElement::Accessor<std::vector<ElementLink<IParticleContainer > > > acc_cc("pfo_ClusterLinks");
21  return &acc_cc;
22  case PFODetails::Track:
23  const static SG::AuxElement::Accessor<std::vector<ElementLink<IParticleContainer > > > acc_tp("pfo_TrackLinks");
24  return &acc_tp;
25  // Note - clients should NEVER set both links to TauShot and HadroniCaloCluster - Tau WG agreed they NEVER will do such a thing (and they are only client)
26  case PFODetails::TauShot:
27  case PFODetails::HadronicCalo:
28  const static SG::AuxElement::Accessor<std::vector<ElementLink<IParticleContainer > > > acc_ts("pfo_TauShotLinks");
29  return &acc_ts;
30  case PFODetails::ChargedPFO:
31  const static SG::AuxElement::Accessor<std::vector<ElementLink<IParticleContainer > > > acc_cpfo("pfo_Charged");
32  return &acc_cpfo;
33  case PFODetails::NeutralPFO:
34  const static SG::AuxElement::Accessor<std::vector<ElementLink<IParticleContainer > > > acc_npfo("pfo_Neutral");
35  return &acc_npfo;
36  case PFODetails::TauTrack:
37  const static SG::AuxElement::Accessor<std::vector<ElementLink<IParticleContainer > > > acc_tt("pfo_TauTrack");
38  return &acc_tt;
39  }//switch
40  return nullptr;
41  }//getAccessor
42  };
43 
44  PFO_v1::PFO_v1() : IParticle() {
45  }
46 
47  PFO_v1::PFO_v1(const PFO_v1& other) : IParticle(other){
48  this->makePrivateStore(other);
49  }
50 
51 
52  double PFO_v1::pt() const {
53 
54  const static ConstAccessor<float> accPt("pt");
55  float pt = accPt(*this);
56 
57  return pt;
58  }
59 
60  double PFO_v1::eta() const {
61 
62  const static ConstAccessor<float> accEta("eta");
63  float eta = accEta(*this);
64  return eta;
65  }
66 
67  double PFO_v1::phi() const {
68 
69  const static ConstAccessor<float> accPhi("phi");
70  float phi = accPhi(*this);
71  return phi;
72  }
73 
74  double PFO_v1::m() const {
75 
76  const static ConstAccessor<float> accM("m");
77  float M = accM(*this);
78  return M;
79  }
80 
81  double PFO_v1::e() const {
82 
83  const static ConstAccessor<float> accPt("pt");
84  float pt = accPt(*this);
85 
86  if (pt < 0.0) return -genvecP4().E();
87  else return genvecP4().E();
88  }
89 
90  double PFO_v1::rapidity() const {
91 
92  return genvecP4().Rapidity();
93  }
94 
95  PFO_v1::FourMom_t PFO_v1::p4() const {
96  FourMom_t p4;
97  p4.SetPtEtaPhiM( pt(), eta(), phi(),m());
98  return p4;
99  }
100 
101  PFO_v1::GenVecFourMom_t PFO_v1::genvecP4() const {
102  return GenVecFourMom_t( pt(), eta(), phi(),m());
103  }
104 
105 
107  void PFO_v1::setP4(const FourMom_t& vec) {
108 
109  const static Accessor<float> accPt("pt");
110  accPt(*this) = vec.Pt();
111 
112  const static Accessor<float> accEta("eta");
113  accEta(*this) = vec.Eta();
114 
115  const static Accessor<float> accPhi("phi");
116  accPhi(*this) = vec.Phi();
117 
118  const static Accessor<float> accM("m");
119  accM(*this) = vec.M();
120 
121  }
122 
124  void PFO_v1::setP4(float pt, float eta, float phi, float m) {
125 
126  const static Accessor<float> accPt("pt");
127  accPt(*this) = pt;
128 
129  const static Accessor<float> accEta("eta");
130  accEta(*this) = eta;
131 
132  const static Accessor<float> accPhi("phi");
133  accPhi(*this) = phi;
134 
135  const static Accessor<float> accM("m");
136  accM(*this) = m;
137  }
138 
139 
140  Type::ObjectType PFO_v1::type() const {
141  return Type::ParticleFlow;
142  }
143 
144  PFO_v1::FourMom_t PFO_v1::p4EM() const {
145 
146  if (this->isCharged()) return p4();
147 
148  FourMom_t p4EM;
149 
150  //change to use pt, eta, phi ,e
151  const static ConstAccessor<float> accPt("ptEM");
152  const static ConstAccessor<float> accEta("eta");
153  const static ConstAccessor<float> accPhi("phi");
154  const static ConstAccessor<float> accM("mEM");
155 
156  p4EM.SetPtEtaPhiM(accPt(*this), accEta(*this), accPhi(*this), accM(*this));
157  return p4EM;
158  }
159 
160  PFO_v1::GenVecFourMom_t PFO_v1::genvecP4EM() const {
161 
162  if (this->isCharged()) return genvecP4();
163 
164  //change to use pt, eta, phi ,e
165  const static ConstAccessor<float> accPt("ptEM");
166  const static ConstAccessor<float> accEta("eta");
167  const static ConstAccessor<float> accPhi("phi");
168  const static ConstAccessor<float> accM("mEM");
169 
170  return GenVecFourMom_t(accPt(*this), accEta(*this), accPhi(*this), accM(*this));
171  }
172 
173  void PFO_v1::setP4EM(const FourMom_t& p4EM) {
174  const static Accessor<float> accPt("ptEM");
175  accPt(*this) = p4EM.Pt();
176 
177  const static Accessor<float> accEta("eta");
178  accEta(*this) = p4EM.Eta();
179 
180  const static Accessor<float> accPhi("phi");
181  accPhi(*this) = p4EM.Phi();
182 
183  const static Accessor<float> accM("mEM");
184  accM(*this) = p4EM.M();
185 
186  }
187 
188  void PFO_v1::setP4EM(float pt, float eta, float phi, float m) {
189 
190  const static Accessor<float> accPt("ptEM");
191  accPt(*this) = pt;
192 
193  const static Accessor<float> accEta("eta");
194  accEta(*this) = eta;
195 
196  const static Accessor<float> accPhi("phi");
197  accPhi(*this) = phi;
198 
199  const static Accessor<float> accM("mEM");
200  accM(*this) = m;
201 
202  }
203 
204  double PFO_v1::ptEM() const {
205 
206  if (this->isCharged()) return this->pt();
207 
208  const static ConstAccessor<float> accPt("ptEM");
209  float pt = accPt(*this);
210 
211  return pt;
212  }
213 
214  // JM: this seem to be idetical to eta(), right?
215  double PFO_v1::etaEM() const {
216  return eta();
217  }
218 
219  // JM: this seem to be idetical to phi(), right?
220  double PFO_v1::phiEM() const {
221  return phi();
222  }
223 
224  double PFO_v1::mEM() const {
225 
226  if (this->isCharged()) return this->m();
227  const static ConstAccessor<float> accM("mEM");
228  float M = accM(*this);
229 
230  return M;
231  }
232 
233  double PFO_v1::eEM() const {
234 
235  if (this->isCharged()) return this->e();
236 
237  const static ConstAccessor<float> accPt("ptEM");
238  float pt = accPt(*this);
239 
240  if (pt < 0.0) return -p4EM().E();
241  else return p4EM().E();
242  }
243 
244 
245  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER(PFO_v1, float, bdtPi0Score, setBDTPi0Score)
246  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER(PFO_v1, float, centerMag, setCenterMag)
247  AUXSTORE_PRIMITIVE_SETTER_AND_GETTER(PFO_v1, float, charge, setCharge)
248 
249  // Define function to avoid issues due to numeric precision when comparing pFO charge to zero
250  // Discussed in https://its.cern.ch/jira/browse/ATLJETMET-796
251  bool PFO_v1::isCharged() const{
252  return (fabs(this->charge())>1e-6);
253  }
254 
257  template<> void PFO_v1::setAttribute(const std::string& AttributeType, const xAOD::PFODetails::PFOLeptonType& anAttribute){
258  uint16_t uint16_variable = static_cast<uint16_t>(anAttribute);
259  this->setAttribute<uint16_t>(AttributeType, uint16_variable);
260  }
261 
262  template<> bool PFO_v1::attribute(const std::string& AttributeType, xAOD::PFODetails::PFOLeptonType& anAttribute) const {
263  bool isValid = false;
264  uint16_t internalAttribute;
265  isValid = attribute<uint16_t>(AttributeType,internalAttribute);
266  if (!isValid) return false;
267  else{
268  anAttribute = static_cast<xAOD::PFODetails::PFOLeptonType>(internalAttribute);
269  return true;
270  }
271 
272  }
273 
274  template<> void PFO_v1::setAttribute(PFODetails::PFOAttributes AttributeType, const float& anAttribute) {
275  if (this->isJetETMissFloatForCompression(AttributeType)){
276  float dummy = anAttribute*s_floatCompressionFactor;
277  const static int maxIntSize = 1000000000;
278  int internalAttribute = maxIntSize;
279  if (dummy < 0) internalAttribute *= -1;//if we had a large -ve energy, then we should set the max size to a -ve value
280  if (dummy < maxIntSize && dummy > -maxIntSize) internalAttribute = static_cast<int>(dummy);
281  setAttribute<int>(AttributeType, internalAttribute);
282  }
283  else {
284  float internalAttribute = anAttribute;
285  (*(PFOAttributesAccessor_v1<float>::constAccessor(AttributeType)))(*this) = internalAttribute;
286  }
287  }
288 
289  template<> bool PFO_v1::attribute(PFODetails::PFOAttributes AttributeType, float& anAttribute) const {
290  bool isValid = false;
291  if (this->isJetETMissFloatForCompression(AttributeType)){
292  int internalAttribute;
293  isValid = attribute<int>(AttributeType,internalAttribute);
294  if (isValid && 0 != internalAttribute) anAttribute = static_cast<float>(internalAttribute)/s_floatCompressionFactor;
295  else anAttribute = 0.0;
296  return isValid;
297  }
298  else{
299  const ConstAccessor<float>* p_acc = PFOAttributesAccessor_v1<float>::constAccessor(AttributeType);
300  //check if accessor pointer is NULL
301  if( ! p_acc ) { return false ;}
302  //check if variable is avaialable
303  if( ! p_acc->isAvailable( *this ) ) return false;
304  //set variable and return true
305  anAttribute =( *p_acc )( *this );
306  return true;
307 
308  }
309  return isValid;
310  }
311 
312  bool PFO_v1::isJetETMissFloatForCompression(PFODetails::PFOAttributes AttributeType) const{
313  switch (AttributeType) {
314  case PFODetails::PFOAttributes::eflowRec_LATERAL:
315  case PFODetails::PFOAttributes::eflowRec_LONGITUDINAL:
316  case PFODetails::PFOAttributes::eflowRec_SECOND_R:
317  case PFODetails::PFOAttributes::eflowRec_CENTER_LAMBDA:
318  case PFODetails::PFOAttributes::eflowRec_FIRST_ENG_DENS:
319  case PFODetails::PFOAttributes::eflowRec_ENG_FRAC_MAX:
320  case PFODetails::PFOAttributes::eflowRec_ISOLATION:
321  case PFODetails::PFOAttributes::eflowRec_ENG_BAD_CELLS:
322  case PFODetails::PFOAttributes::eflowRec_N_BAD_CELLS:
323  case PFODetails::PFOAttributes::eflowRec_BADLARQ_FRAC:
324  case PFODetails::PFOAttributes::eflowRec_ENG_POS:
325  case PFODetails::PFOAttributes::eflowRec_SIGNIFICANCE:
326  case PFODetails::PFOAttributes::eflowRec_AVG_LAR_Q:
327  case PFODetails::PFOAttributes::eflowRec_AVG_TILE_Q:
328  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EM:
329  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EM3:
330  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_HEC0:
331  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_HEC:
332  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_Tile0:
333  case PFODetails::PFOAttributes::eflowRec_TIMING:
334  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_PreSamplerB:
335  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EMB1:
336  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EMB2:
337  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EMB3:
338  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_PreSamplerE:
339  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EME1:
340  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EME2:
341  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_EME3:
342  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_HEC1:
343  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_HEC2:
344  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_HEC3:
345  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileBar0:
346  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileBar1:
347  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileBar2:
348  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileGap1:
349  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileGap2:
350  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileGap3:
351  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileExt0:
352  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileExt1:
353  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_TileExt2:
354  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_FCAL0:
355  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_FCAL1:
356  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_FCAL2:
357  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_MINIFCAL0:
358  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_MINIFCAL1:
359  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_MINIFCAL2:
360  case PFODetails::PFOAttributes::eflowRec_LAYERENERGY_MINIFCAL3:
361  return true;
362  default:
363  return false;
364  }
365  }
366 
367 
372  template<> void PFO_v1::setAttribute(PFODetails::PFOAttributes AttributeType, const double& anAttribute) {
373  float internalAttribute = static_cast<float>(anAttribute);
374  setAttribute<float>(AttributeType, internalAttribute);
375  }
376 
377  template<> bool PFO_v1::attribute(PFODetails::PFOAttributes AttributeType, double& anAttribute) const {
378 
379  float internalAttribute;
380  bool isValid = attribute<float>(AttributeType,internalAttribute);
381  if (isValid) anAttribute = static_cast<double>(internalAttribute);
382  return isValid;
383 
384  }
385 
386  template<> void PFO_v1::setAttribute(const std::string& AttributeType, const double& anAttribute) {
387 
388  float internalAttribute = static_cast<float>(anAttribute);
389  setAttribute<float>(AttributeType, internalAttribute);
390 
391  }
392 
393  template<> bool PFO_v1::attribute(const std::string& AttributeType, double& anAttribute) const {
394 
395  float internalAttribute;
396  bool isValid = attribute<float>(AttributeType,internalAttribute);
397  if (isValid) anAttribute = static_cast<float>(internalAttribute);
398  return isValid;
399 
400  }
401 
402  bool PFO_v1::getClusterMoment(float& theMoment, xAOD::CaloCluster::MomentType momentType) const {
403 
404  const xAOD::CaloCluster *theCluster = this->cluster(0);
405 
406  double myMoment = -99.0;
407 
408  if (theCluster){
409  bool gotMoment = theCluster->retrieveMoment( momentType, myMoment);
410  //if cluster does not have the moment, check if the PFO stored it - this can happen due to the way we redirect the ElementLinks to clusters in PFONeutralCreatorAlgorithm
411  if (!gotMoment) return getClusterMomentFromPFO(theMoment,momentType);
412  else if (gotMoment) {
413  theMoment = static_cast<float>(myMoment);
414  return true;
415  }
416  }
417  else return getClusterMomentFromPFO(theMoment,momentType);
418 
419  return false;
420 
421  }
422 
423  bool PFO_v1::getClusterMomentFromPFO(float& theMoment, xAOD::CaloCluster::MomentType momentType) const{
424  xAOD::PFODetails::PFOAttributes myAttribute;
425  bool haveGotAttributeInMap = this->getAttributeName_FromClusterMoment(myAttribute,momentType);
426  if (!haveGotAttributeInMap) return false;
427  else {
428  bool isRetrieved = this->attribute(myAttribute, theMoment);
429  return isRetrieved;
430  }
431  }
432 
433  bool PFO_v1::getAttributeName_FromClusterMoment(xAOD::PFODetails::PFOAttributes& myAttribute, xAOD::CaloCluster::MomentType& momentType) const{
434 
435  if (momentType == xAOD::CaloCluster::ENG_FRAC_CORE) {
436  myAttribute = xAOD::PFODetails::eflowRec_ENG_FRAC_CORE;
437  return true;
438  }
439  else if (momentType == xAOD::CaloCluster::FIRST_ENG_DENS) {
440  myAttribute = xAOD::PFODetails::eflowRec_FIRST_ENG_DENS;
441  return true;
442  }
443  else if (momentType == xAOD::CaloCluster::CENTER_LAMBDA) {
444  myAttribute = xAOD::PFODetails::eflowRec_CENTER_LAMBDA;
445  return true;
446  }
447  else if (momentType == xAOD::CaloCluster::SECOND_R) {
448  myAttribute = xAOD::PFODetails::eflowRec_SECOND_R;
449  return true;
450  }
451  else if (momentType == xAOD::CaloCluster::DELTA_ALPHA) {
452  myAttribute = xAOD::PFODetails::eflowRec_DELTA_ALPHA;
453  return true;
454  }
455  else if (momentType == xAOD::CaloCluster::LATERAL) {
456  myAttribute = xAOD::PFODetails::eflowRec_LATERAL;
457  return true;
458  }
459  else if (momentType == xAOD::CaloCluster::LONGITUDINAL) {
460  myAttribute = xAOD::PFODetails::eflowRec_LONGITUDINAL;
461  return true;
462  }
463  else if (momentType == xAOD::CaloCluster::SECOND_LAMBDA) {
464  myAttribute = xAOD::PFODetails::eflowRec_SECOND_LAMBDA;
465  return true;
466  }
467  else if (momentType == xAOD::CaloCluster::ISOLATION) {
468  myAttribute = xAOD::PFODetails::eflowRec_ISOLATION;
469  return true;
470  }
471  else if (momentType == xAOD::CaloCluster::ENG_FRAC_MAX) {
472  myAttribute = xAOD::PFODetails::eflowRec_ENG_FRAC_MAX;
473  return true;
474  }
475  else if (momentType == xAOD::CaloCluster::ENG_BAD_CELLS) {
476  myAttribute = xAOD::PFODetails::eflowRec_ENG_BAD_CELLS;
477  return true;
478  }
479  else if (momentType == xAOD::CaloCluster::N_BAD_CELLS) {
480  myAttribute = xAOD::PFODetails::eflowRec_N_BAD_CELLS;
481  return true;
482  }
483  else if (momentType == xAOD::CaloCluster::BADLARQ_FRAC) {
484  myAttribute = xAOD::PFODetails::eflowRec_BADLARQ_FRAC;
485  return true;
486  }
487  else if (momentType == xAOD::CaloCluster::ENG_POS) {
488  myAttribute = xAOD::PFODetails::eflowRec_ENG_POS;
489  return true;
490  }
491  else if (momentType == xAOD::CaloCluster::SIGNIFICANCE) {
492  myAttribute = xAOD::PFODetails::eflowRec_SIGNIFICANCE;
493  return true;
494  }
495  else if (momentType == xAOD::CaloCluster::CELL_SIGNIFICANCE) {
496  myAttribute = xAOD::PFODetails::eflowRec_CELL_SIGNIFICANCE;
497  return true;
498  }
499  else if (momentType == xAOD::CaloCluster::CELL_SIG_SAMPLING) {
500  myAttribute = xAOD::PFODetails::eflowRec_CELL_SIG_SAMPLING;
501  return true;
502  }
503  else if (momentType == xAOD::CaloCluster::AVG_LAR_Q) {
504  myAttribute = xAOD::PFODetails::eflowRec_AVG_LAR_Q;
505  return true;
506  }
507  else if (momentType == xAOD::CaloCluster::AVG_TILE_Q) {
508  myAttribute = xAOD::PFODetails::eflowRec_AVG_TILE_Q;
509  return true;
510  }
511  else if (momentType == xAOD::CaloCluster::EM_PROBABILITY){
512  myAttribute = xAOD::PFODetails::eflowRec_EM_PROBABILITY;
513  return true;
514  }
515  else if (momentType == xAOD::CaloCluster::ENG_CALIB_TOT){
516  myAttribute = xAOD::PFODetails::eflowRec_ENG_CALIB_TOT;
517  return true;
518  }
519  else if (momentType == xAOD::CaloCluster::ENG_CALIB_FRAC_EM){
520  myAttribute = xAOD::PFODetails::eflowRec_ENG_CALIB_FRAC_EM;
521  return true;
522  }
523  else if (momentType == xAOD::CaloCluster::ENG_CALIB_FRAC_HAD){
524  myAttribute = xAOD::PFODetails::eflowRec_ENG_CALIB_FRAC_HAD;
525  return true;
526  }
527  else if (momentType == xAOD::CaloCluster::ENG_CALIB_FRAC_REST){
528  myAttribute = xAOD::PFODetails::eflowRec_ENG_CALIB_FRAC_REST;
529  return true;
530  }
531 
532  else return false;
533 
534  }
535 
536  bool PFO_v1::setVertexLink(const ElementLink< xAOD::VertexContainer>& theVertexLink){
537  const static Accessor<ElementLink<xAOD::VertexContainer> > acc("pfo_vertex");
538  acc(*this) = theVertexLink;
539  acc(*this).toPersistent();
540  return true;
541  }
542 
543  bool PFO_v1::setTrackLink(const ElementLink<xAOD::TrackParticleContainer>& theTrack){
544  ElementLink< xAOD::IParticleContainer > myTrackLink;
545  this->convertLink(myTrackLink,theTrack);
546  return this->setAssociatedParticleLink(PFODetails::Track, myTrackLink);
547  }
548 
549  bool PFO_v1::setClusterLink(const ElementLink<xAOD::CaloClusterContainer>& theCluster){
550  ElementLink< xAOD::IParticleContainer > myClusterLink;
551  this->convertLink(myClusterLink,theCluster);
552  return this->setAssociatedParticleLink(PFODetails::CaloCluster, myClusterLink);
553  }
554 
555  bool PFO_v1::addClusterLink(const ElementLink<xAOD::CaloClusterContainer>& theCluster){
556  ElementLink< xAOD::IParticleContainer > myClusterLink;
557  this->convertLink(myClusterLink,theCluster);
558  return this->addAssociatedParticleLink(PFODetails::CaloCluster, myClusterLink);
559  }
560 
561  bool PFO_v1::setAssociatedParticleLink(PFODetails::PFOParticleType ParticleType, const ElementLink<IParticleContainer>& theParticle){
562  std::vector<ElementLink<xAOD::IParticleContainer> > theLinks(1,theParticle);
563  return this->setAssociatedParticleLinks(ParticleType,theLinks);
564  }
565 
566  void PFO_v1::setAssociatedParticleLink(const std::string& ParticleType, const ElementLink<IParticleContainer>& theParticle){
567  std::vector<ElementLink<xAOD::IParticleContainer> > theLinks(1,theParticle);
568  this->setAssociatedParticleLinks(ParticleType,theLinks);
569  }
570 
571  bool PFO_v1::addAssociatedParticleLink(PFODetails::PFOParticleType ParticleType, const ElementLink<IParticleContainer>& theParticle) {
572  const Accessor<std::vector<ElementLink<IParticleContainer > > >* p_acc = PFOParticleTypeMapper_temp::getAccessor(ParticleType);
573  if (!p_acc) return false;
574  else{
575  if (!p_acc->isAvailable(*this)) return false;
576  else{
577  std::vector<ElementLink<IParticleContainer> > storedContainer = (*p_acc)(*this);
578  storedContainer.push_back(theParticle);
579  (*p_acc)(*this) = storedContainer;
580  return true;
581  }
582  }
583  }
584 
585  void PFO_v1::addAssociatedParticleLink(const std::string& ParticleType, const ElementLink<IParticleContainer>& theParticle) {
586 
587  const static Accessor<std::vector<ElementLink<IParticleContainer > > > acc(ParticleType);
588  std::vector<ElementLink<IParticleContainer> > storedContainer = acc(*this);
589 
590  ElementLink<xAOD::IParticleContainer> newLink;
591  newLink.setElement(*theParticle);
592  newLink.setStorableObject(theParticle.getStorableObjectRef());
593  newLink.toPersistent();
594  storedContainer.push_back(newLink);
595 
596  acc(*this) = storedContainer;
597 
598  }
599 
600  bool PFO_v1::setAssociatedParticleLinks(PFODetails::PFOParticleType ParticleType, const std::vector<ElementLink<IParticleContainer> >& theParticles) {
601 
602  const Accessor<std::vector<ElementLink<IParticleContainer > > >* p_acc = PFOParticleTypeMapper_temp::getAccessor(ParticleType);
603  if (!p_acc) return false;
604  else{
605  (*p_acc)(*this) = theParticles;
606  return true;
607  }
608  }
609 
610  bool PFO_v1::associatedParticles(PFODetails::PFOParticleType ParticleType, std::vector<const IParticle*>& theParticles ) const{
611 
612  const ConstAccessor<std::vector<ElementLink<IParticleContainer > > >* p_acc = PFOParticleTypeMapper_temp::getAccessor(ParticleType);
613  if (!p_acc) return false;
614  else{
615  if (!p_acc->isAvailable(*this)) return false;
616  else{
617  const std::vector<ElementLink<IParticleContainer> >& theLinks = (*p_acc)(*this);
618  std::vector<ElementLink<IParticleContainer> >::const_iterator firstLink = theLinks.begin();
619  for (; firstLink != theLinks.end(); ++firstLink) {
620  if ( (*firstLink).isValid()) theParticles.push_back(**firstLink);
621  else theParticles.push_back(nullptr);
622  }
623  return true;
624  }
625  }
626 
627  }
628 
629  void PFO_v1::setAssociatedParticleLinks(const std::string& ParticleType, const std::vector<ElementLink<IParticleContainer> >& theParticles) {
630 
631  //Given we do not know in advance in the POOL convertors about these containers, we set toPersistent() internally here.
632 
633  std::vector<ElementLink<IParticleContainer> > storedContainer;
634  std::vector<ElementLink<IParticleContainer> >::const_iterator firstParticle = theParticles.begin();
635  for (; firstParticle != theParticles.end(); ++firstParticle){
636  ElementLink<xAOD::IParticleContainer> myLink = *firstParticle;
637  myLink.toPersistent();
638  storedContainer.push_back( myLink );
639  }
640 
641  const static Accessor<std::vector<ElementLink<IParticleContainer > > > acc(ParticleType);
642  acc(*this) = storedContainer;
643  }
644 
645  bool PFO_v1::associatedParticles(const std::string& ParticleType, std::vector<const IParticle*>& theParticles ) const{
646  const static ConstAccessor<std::vector<ElementLink<IParticleContainer > > > acc(ParticleType);
647  if (!acc.isAvailable(*this)) return false;
648  else{
649  const std::vector<ElementLink<IParticleContainer> >& theLinks = acc(*this);
650  std::vector<ElementLink<IParticleContainer> >::const_iterator firstLink = theLinks.begin();
651  for (; firstLink != theLinks.end(); ++firstLink) {
652  if ( (*firstLink).isValid()) theParticles.push_back(**firstLink);
653  else theParticles.push_back(nullptr);
654  }
655  return true;
656  }
657  }
658 
659  unsigned int PFO_v1::nCaloCluster() const {
660  const ConstAccessor<std::vector<ElementLink<IParticleContainer > > >* p_acc = PFOParticleTypeMapper_temp::getAccessor(PFODetails::CaloCluster);
661  if (!p_acc) return 0;
662  else if (!p_acc->isAvailable(*this)) {return 0;}
663  else{
664  const std::vector<ElementLink<IParticleContainer> >& theLinks = (*p_acc)(*this);
665  return theLinks.size();
666  }
667  }
668 
669  const CaloCluster* PFO_v1::cluster(unsigned int index) const {
670 
671  const ConstAccessor<std::vector<ElementLink<IParticleContainer > > >* p_acc = PFOParticleTypeMapper_temp::getAccessor(PFODetails::CaloCluster);
672  if (!p_acc) return nullptr;
673  else if (!p_acc->isAvailable(*this)) {return nullptr;}
674  else {
675  const std::vector<ElementLink<IParticleContainer> >& theLinks = (*p_acc)(*this);
676  if(index<theLinks.size()) {
677  ElementLink<IParticleContainer> theLink = theLinks[index];
678  if (theLink.isValid()){
679  const IParticle *theParticle = *theLink;
680  if (nullptr == theParticle) return nullptr;
681  else if (Type::CaloCluster == theParticle->type()) return static_cast<const CaloCluster*>(theParticle);
682  else return nullptr;
683  }
684  else return nullptr;
685  }
686  else return nullptr;
687  }
688 
689  }
690 
691  const TrackParticle* PFO_v1::track(unsigned int index) const {
692 
693  const ConstAccessor<std::vector<ElementLink<IParticleContainer > > >* p_acc = PFOParticleTypeMapper_temp::getAccessor(PFODetails::Track);
694  if (!p_acc) return nullptr;
695  else if (!p_acc->isAvailable(*this)) {return nullptr;}
696  else {
697  const std::vector<ElementLink<IParticleContainer> >& theLinks = (*p_acc)(*this);
698  if(index<theLinks.size()) {
699  ElementLink<IParticleContainer> theLink = theLinks[index];
700  if (theLink.isValid()){
701  const IParticle *theParticle = *theLinks[index];
702  if (nullptr == theParticle) return nullptr;
703  else if (Type::TrackParticle == theParticle->type()) return static_cast<const TrackParticle*>(theParticle);
704  else return nullptr;
705  }
706  else return nullptr;
707  }
708  else return nullptr;
709  }
710  }
711 
712  const xAOD::Vertex* PFO_v1::vertex() const{
713  const static ConstAccessor<ElementLink<xAOD::VertexContainer> > acc("pfo_vertex");
714  if (!acc.isAvailable(*this)) return nullptr;
715  else{
716  const ElementLink<xAOD::VertexContainer>& tempVertexLink = acc(*this);
717  if (tempVertexLink.isValid()) return *acc(*this);
718  else return nullptr;
719  }
720  }
721 
722  TLorentzVector PFO_v1::GetVertexCorrectedFourVec(const xAOD::Vertex& vertexToCorrectTo) const{
723  TVector3 theVertexVector(vertexToCorrectTo.x(), vertexToCorrectTo.y(), vertexToCorrectTo.z());
724  return GetVertexCorrectedFourVec(theVertexVector);
725  }
726 
727  TLorentzVector PFO_v1::GetVertexCorrectedFourVec(const TVector3& vertexToCorrectTo) const{
728 
729  TLorentzVector theNewVector(0.0,0.0,0.0,0.0);
730  theNewVector.SetPtEtaPhiM(this->pt(), this->eta(), this->phi(), this->m());
731 
732  this->VertexCorrectTheFourVector(vertexToCorrectTo,theNewVector);
733  return theNewVector;
734 
735  }
736 
737  TLorentzVector PFO_v1::GetVertexCorrectedEMFourVec(const xAOD::Vertex& vertexToCorrectTo) const{
738  TVector3 theVertexVector(vertexToCorrectTo.x(), vertexToCorrectTo.y(), vertexToCorrectTo.z());
739  return GetVertexCorrectedEMFourVec(theVertexVector);
740  }
741 
742  TLorentzVector PFO_v1::GetVertexCorrectedEMFourVec(const TVector3& vertexToCorrectTo) const{
743 
744  TLorentzVector theNewVector(0.0,0.0,0.0,0.0);
745  theNewVector.SetPtEtaPhiM(this->ptEM(), this->etaEM(), this->phiEM(), this->mEM());
746 
747  this->VertexCorrectTheFourVector(vertexToCorrectTo,theNewVector);
748  return theNewVector;
749 
750  }
751 
752  void PFO_v1::VertexCorrectTheFourVector(const TVector3& vertexToCorrectTo, TLorentzVector& theFourVector) const{
753 
754  float clusterEta = theFourVector.Eta();
755  float centerMag = this->centerMag();
756 
757  float radius = centerMag/cosh(clusterEta);
758 
759  float EtaVertexCorrection = 0.0, PhiVertexCorrection = 0.0;
760  float clusterPhi = theFourVector.Phi();
761 
762  if (radius > 1.0 && centerMag > 1e-3){
763  EtaVertexCorrection = (-vertexToCorrectTo.Z()/cosh(clusterEta) + (vertexToCorrectTo.X()*cos(clusterPhi) + vertexToCorrectTo.Y()*sin(clusterPhi))*tanh(clusterEta))/radius;
764  PhiVertexCorrection = (vertexToCorrectTo.X()*sin(clusterPhi) - vertexToCorrectTo.Y()*cos(clusterPhi))/radius;
765  }
766 
767  float etaVertexCorrected = clusterEta + EtaVertexCorrection;
768  float p = std::sqrt(theFourVector.E()*theFourVector.E()-theFourVector.M()*theFourVector.M());
769  float ptVertexCorrected = p/cosh(etaVertexCorrected);
770  theFourVector.SetPtEtaPhiM(ptVertexCorrected, etaVertexCorrected, clusterPhi + PhiVertexCorrection, theFourVector.M());
771 
772  }
773 
774  void PFO_v1::toPersistent(){
775 
776  //setup known links for persistification, if they exist
777  //if you added your own consituents the links will not be correctly persistified
778 
779  //clusters
780  const Accessor<std::vector<ElementLink<IParticleContainer > > >* p_accClusters = PFOParticleTypeMapper_temp::getAccessor(PFODetails::CaloCluster);
781  if (p_accClusters){
782  const Accessor<std::vector<ElementLink<IParticleContainer > > >& accClusters = *p_accClusters;
783  if ( accClusters.isAvailableWritable(*this) ){
784  std::vector<ElementLink<IParticleContainer> >& theClusterLinks = accClusters(*this);
785  std::vector< ElementLink< IParticleContainer > >::iterator firstClus = theClusterLinks.begin();
786  std::vector< ElementLink< IParticleContainer > >::iterator lastClus = theClusterLinks.end();
787  for (; firstClus != lastClus; ++firstClus) firstClus->toPersistent();
788  }
789  }
790 
791  //tracks
792  const Accessor<std::vector<ElementLink<IParticleContainer > > >* p_accTracks = PFOParticleTypeMapper_temp::getAccessor(PFODetails::Track);
793  if (p_accTracks){
794  const Accessor<std::vector<ElementLink<IParticleContainer > > >& accTracks = *p_accTracks;
795  if ( accTracks.isAvailableWritable(*this) ){
796  std::vector<ElementLink<IParticleContainer> >& theTrackLinks = accTracks(*this);
797  std::vector< ElementLink< IParticleContainer > >::iterator firstTrack = theTrackLinks.begin();
798  std::vector< ElementLink< IParticleContainer > >::iterator lastTrack = theTrackLinks.end();
799  for (; firstTrack != lastTrack; ++firstTrack) firstTrack->toPersistent();
800 
801  }
802  }
803 
804  //shots
805  const Accessor<std::vector<ElementLink<IParticleContainer > > >* p_accShots = PFOParticleTypeMapper_temp::getAccessor(PFODetails::TauShot);
806  if (p_accShots){
807  const Accessor<std::vector<ElementLink<IParticleContainer > > >& accShots = *p_accShots;
808  if ( accShots.isAvailableWritable(*this) ){
809  std::vector<ElementLink<IParticleContainer> >& theShotLinks = accShots(*this);
810  std::vector< ElementLink< IParticleContainer > >::iterator firstShot = theShotLinks.begin();
811  std::vector< ElementLink< IParticleContainer > >::iterator lastShot = theShotLinks.end();
812  for (; firstShot != lastShot; ++firstShot) firstShot->toPersistent();
813  }
814  }
815 
816  }
817 
818 
819 } // namespace xAOD
xAOD::CaloCluster_v1::SECOND_R
@ SECOND_R
Second Moment in .
Definition: CaloCluster_v1.h:123
xAOD::AUXSTORE_PRIMITIVE_SETTER_AND_GETTER
AUXSTORE_PRIMITIVE_SETTER_AND_GETTER(BTagging_v1, float, IP2D_pb, setIP2D_pb) AUXSTORE_PRIMITIVE_SETTER_AND_GETTER(BTagging_v1
xAOD::PFODetails::eflowRec_LAYERENERGY_EME1
@ eflowRec_LAYERENERGY_EME1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:67
xAOD::PFODetails::eflowRec_LONGITUDINAL
@ eflowRec_LONGITUDINAL
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:41
xAOD::Vertex_v1::x
float x() const
Returns the x position.
xAOD::PFODetails::eflowRec_LAYERENERGY_HEC1
@ eflowRec_LAYERENERGY_HEC1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:70
xAOD::PFODetails::PFOLeptonType
PFOLeptonType
This enum is used to label the lepton type of the neutral PFO.
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:179
xAOD::CaloCluster_v1::ENG_FRAC_CORE
@ ENG_FRAC_CORE
Energy fraction of the sum of the hottest cells in each sampling.
Definition: CaloCluster_v1.h:142
xAOD::PFODetails::eflowRec_LAYERENERGY_EMB3
@ eflowRec_LAYERENERGY_EMB3
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:65
xAOD::PFO_v1::getClusterMomentFromPFO
bool getClusterMomentFromPFO(float &theMoment, xAOD::CaloCluster::MomentType momentType) const
Gives access to cluster moments direct from the PFO - getClusterMoment calls this,...
Definition: PFO_v1.cxx:423
xAOD::PFO_v1::GenVecFourMom_t
ROOT::Math::LorentzVector< ROOT::Math::PtEtaPhiM4D< double > > GenVecFourMom_t
Base 4 Momentum type (GenVector version)
Definition: PFO_v1.h:71
python.SystemOfUnits.m
int m
Definition: SystemOfUnits.py:91
xAOD::PFO_v1::pt
virtual double pt() const
The transverse momentum ( ) of the particle.
Definition: PFO_v1.cxx:52
xAOD::PFO_v1::convertLink
void convertLink(ElementLink< IParticleContainer > &IParticleLink, const ElementLink< T > &templateLink)
Convert a generic link to an IParticleLink - perhaps this will go in a dedicated helper tool?
xAOD::PFODetails::eflowRec_SIGNIFICANCE
@ eflowRec_SIGNIFICANCE
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:49
PFO_v1.h
xAOD::CaloCluster_v1::EM_PROBABILITY
@ EM_PROBABILITY
Classification probability to be em-like.
Definition: CaloCluster_v1.h:173
xAOD::PFODetails::eflowRec_CELL_SIGNIFICANCE
@ eflowRec_CELL_SIGNIFICANCE
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:50
SG::Accessor
Helper class to provide type-safe access to aux data.
Definition: Control/AthContainers/AthContainers/Accessor.h:68
AuxStoreAccessorMacros.h
xAOD::PFOParticleTypeMapper_temp::getAccessor
static const SG::AuxElement::Accessor< std::vector< ElementLink< IParticleContainer > > > * getAccessor(PFODetails::PFOParticleType type)
Definition: PFO_v1.cxx:17
xAOD::PFODetails::TauShot
@ TauShot
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:171
xAOD::PFODetails::PFOAttributes
PFOAttributes
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:28
xAOD::PFO_v1::genvecP4
GenVecFourMom_t genvecP4() const
The full 4-momentum of the particle : GenVector version.
Definition: PFO_v1.cxx:101
xAOD::PFODetails::TauTrack
@ TauTrack
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:175
index
Definition: index.py:1
xAOD::PFODetails::eflowRec_LAYERENERGY_TileGap2
@ eflowRec_LAYERENERGY_TileGap2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:77
xAOD::PFOAttributesAccessor_v1::constAccessor
static const SG::AuxElement::Accessor< T > * constAccessor(xAOD::PFODetails::PFOAttributes variable)
Definition: PFOAttributesAccessor_v1.h:32
xAOD::PFODetails::eflowRec_LAYERENERGY_HEC0
@ eflowRec_LAYERENERGY_HEC0
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:55
xAOD::PFO_v1::track
const TrackParticle * track(unsigned int index) const
Retrieve a const pointer to a Rec::TrackParticle.
Definition: PFO_v1.cxx:691
xAOD::PFODetails::eflowRec_LAYERENERGY_MINIFCAL3
@ eflowRec_LAYERENERGY_MINIFCAL3
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:88
xAOD::PFO_v1::setP4EM
void setP4EM(const FourMom_t &p4EM)
set EM scale 4-vector
Definition: PFO_v1.cxx:173
xAOD::PFODetails::eflowRec_TIMING
@ eflowRec_TIMING
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:58
xAOD::PFO_v1::phiEM
virtual double phiEM() const
get EM scale phi
Definition: PFO_v1.cxx:220
xAOD::PFO_v1::attribute
bool attribute(PFODetails::PFOAttributes AttributeType, T &anAttribute) const
get a PFO Variable via enum
xAOD::CaloCluster_v1::AVG_LAR_Q
@ AVG_LAR_Q
Sum(E_cell_LAr^2 Q_cell_LAr)/Sum(E_cell_LAr^2)
Definition: CaloCluster_v1.h:163
xAOD::PFO_v1::eta
virtual double eta() const
The pseudorapidity ( ) of the particle.
Definition: PFO_v1.cxx:60
xAOD::PFODetails::eflowRec_AVG_TILE_Q
@ eflowRec_AVG_TILE_Q
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:53
xAOD::PFODetails::eflowRec_LAYERENERGY_EM3
@ eflowRec_LAYERENERGY_EM3
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:54
xAOD::PFODetails::eflowRec_BADLARQ_FRAC
@ eflowRec_BADLARQ_FRAC
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:47
xAOD::IParticle::type
virtual Type::ObjectType type() const =0
The type of the object as a simple enumeration.
xAOD::PFODetails::eflowRec_LAYERENERGY_FCAL2
@ eflowRec_LAYERENERGY_FCAL2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:84
xAOD::other
@ other
Definition: TrackingPrimitives.h:509
xAOD::CaloCluster_v1::CENTER_LAMBDA
@ CENTER_LAMBDA
Shower depth at Cluster Centroid.
Definition: CaloCluster_v1.h:136
xAOD::PFODetails::eflowRec_EM_PROBABILITY
@ eflowRec_EM_PROBABILITY
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:89
xAOD::PFODetails::PFOParticleType
PFOParticleType
This enum is used to label the associated particles to the PFO object.
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:168
xAOD
ICaloAffectedTool is abstract interface for tools checking if 4 mom is in calo affected region.
Definition: ICaloAffectedTool.h:24
xAOD::PFODetails::eflowRec_SECOND_LAMBDA
@ eflowRec_SECOND_LAMBDA
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:42
xAOD::PFODetails::eflowRec_LAYERENERGY_EM
@ eflowRec_LAYERENERGY_EM
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:61
xAOD::PFODetails::ChargedPFO
@ ChargedPFO
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:173
vec
std::vector< size_t > vec
Definition: CombinationsGeneratorTest.cxx:12
xAOD::PFODetails::eflowRec_LAYERENERGY_HEC2
@ eflowRec_LAYERENERGY_HEC2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:71
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:55
xAOD::PFO_v1::addClusterLink
bool addClusterLink(const ElementLink< xAOD::CaloClusterContainer > &theCluster)
Add a cluster constituent appends to existing container.
Definition: PFO_v1.cxx:555
xAOD::PFODetails::eflowRec_LAYERENERGY_FCAL1
@ eflowRec_LAYERENERGY_FCAL1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:83
xAOD::PFODetails::eflowRec_LAYERENERGY_TileBar0
@ eflowRec_LAYERENERGY_TileBar0
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:73
xAOD::pt
setRcore setEtHad setFside pt
Definition: TrigPhoton_v1.cxx:106
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
xAOD::PFODetails::eflowRec_LAYERENERGY_TileExt1
@ eflowRec_LAYERENERGY_TileExt1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:80
isValid
bool isValid(const T &p)
Definition: AtlasPID.h:225
xAOD::PFO_v1::associatedParticles
bool associatedParticles(PFODetails::PFOParticleType ParticleType, std::vector< const IParticle * > &theParticles) const
get a vector of PFO constituent particle types via enum
Definition: PFO_v1.cxx:610
xAOD::CaloCluster_v1::ENG_FRAC_MAX
@ ENG_FRAC_MAX
Energy fraction of hottest cell.
Definition: CaloCluster_v1.h:140
xAOD::PFO_v1::setAssociatedParticleLinks
bool setAssociatedParticleLinks(PFODetails::PFOParticleType ParticleType, const std::vector< ElementLink< IParticleContainer > > &theParticles)
Set a vector of PFO constituent particle types via enum - overwrite is allowed.
Definition: PFO_v1.cxx:600
xAOD::IParticle
Class providing the definition of the 4-vector interface.
Definition: Event/xAOD/xAODBase/xAODBase/IParticle.h:41
xAOD::PFODetails::eflowRec_SECOND_R
@ eflowRec_SECOND_R
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:36
xAOD::PFODetails::eflowRec_LAYERENERGY_EMB2
@ eflowRec_LAYERENERGY_EMB2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:64
xAOD::PFO_v1::GetVertexCorrectedFourVec
TLorentzVector GetVertexCorrectedFourVec(const xAOD::Vertex &vertexToCorrectTo) const
Correct 4-vector to point at a vertex.
Definition: PFO_v1.cxx:722
xAOD::PFODetails::eflowRec_ENG_FRAC_CORE
@ eflowRec_ENG_FRAC_CORE
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:33
xAOD::PFO_v1::s_floatCompressionFactor
static const int s_floatCompressionFactor
this defines the factor to compress floats by
Definition: PFO_v1.h:203
xAOD::CaloCluster
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
Definition: Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloCluster.h:19
xAOD::PFO_v1::ptEM
virtual double ptEM() const
get EM scale pt
Definition: PFO_v1.cxx:204
xAOD::PFODetails::eflowRec_LAYERENERGY_PreSamplerB
@ eflowRec_LAYERENERGY_PreSamplerB
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:62
xAOD::PFO_v1::p4EM
FourMom_t p4EM() const
get EM scale 4-vector
Definition: PFO_v1.cxx:144
xAOD::PFODetails::eflowRec_ISOLATION
@ eflowRec_ISOLATION
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:43
PFOAttributesAccessor_v1.h
xAOD::CaloCluster_v1::SECOND_LAMBDA
@ SECOND_LAMBDA
Second Moment in .
Definition: CaloCluster_v1.h:124
xAOD::PFODetails::eflowRec_CENTER_LAMBDA
@ eflowRec_CENTER_LAMBDA
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:35
xAOD::PFO_v1::getAttributeName_FromClusterMoment
bool getAttributeName_FromClusterMoment(xAOD::PFODetails::PFOAttributes &myAttribute, xAOD::CaloCluster::MomentType &momentType) const
Map from cluster moment name onto PFOAttribute name.
Definition: PFO_v1.cxx:433
xAOD::IParticle::FourMom_t
TLorentzVector FourMom_t
Definition of the 4-momentum type.
Definition: Event/xAOD/xAODBase/xAODBase/IParticle.h:69
xAOD::CaloCluster_v1::MomentType
MomentType
Enums to identify different moments.
Definition: CaloCluster_v1.h:120
xAOD::PFO_v1::setClusterLink
bool setClusterLink(const ElementLink< xAOD::CaloClusterContainer > &theCluster)
Set a cluster constituent - does NOT append to existing container
Definition: PFO_v1.cxx:549
xAOD::PFODetails::eflowRec_LAYERENERGY_Tile0
@ eflowRec_LAYERENERGY_Tile0
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:56
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition: Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::phi
setEt phi
Definition: TrigEMCluster_v1.cxx:29
xAOD::PFODetails::eflowRec_DELTA_ALPHA
@ eflowRec_DELTA_ALPHA
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:37
xAOD::PFO_v1::type
virtual Type::ObjectType type() const
The type of the object as a simple enumeration.
Definition: PFO_v1.cxx:140
SG::IAuxElement::index
size_t index() const
Return the index of this element within its container.
xAOD::PFO_v1::e
virtual double e() const
The total energy of the particle.
Definition: PFO_v1.cxx:81
xAOD::PFO_v1::setVertexLink
bool setVertexLink(const ElementLink< xAOD::VertexContainer > &theVertexLink)
Set a vertex link.
Definition: PFO_v1.cxx:536
xAOD::PFOParticleTypeMapper_temp
Definition: PFO_v1.cxx:16
xAOD::CaloCluster_v1
Description of a calorimeter cluster.
Definition: CaloCluster_v1.h:59
xAOD::PFODetails::eflowRec_LAYERENERGY_TileBar2
@ eflowRec_LAYERENERGY_TileBar2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:75
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
xAOD::PFO_v1::nCaloCluster
unsigned int nCaloCluster() const
Find out how many CaloCluster are linked.
Definition: PFO_v1.cxx:659
xAOD::uint16_t
setWord1 uint16_t
Definition: eFexEMRoI_v1.cxx:88
xAOD::PFODetails::eflowRec_LAYERENERGY_MINIFCAL1
@ eflowRec_LAYERENERGY_MINIFCAL1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:86
xAOD::PFODetails::eflowRec_LAYERENERGY_MINIFCAL2
@ eflowRec_LAYERENERGY_MINIFCAL2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:87
xAOD::PFO_v1::setTrackLink
bool setTrackLink(const ElementLink< xAOD::TrackParticleContainer > &theTrack)
Set a track constituent - does NOT append to existing container
Definition: PFO_v1.cxx:543
xAOD::PFODetails::eflowRec_LAYERENERGY_HEC
@ eflowRec_LAYERENERGY_HEC
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:57
xAOD::e
setPy e
Definition: CompositeParticle_v1.cxx:166
xAOD::CaloCluster_v1::CELL_SIG_SAMPLING
@ CELL_SIG_SAMPLING
CaloSample of the cell with the largest |E|/sig.
Definition: CaloCluster_v1.h:161
xAOD::charge
charge
Definition: TrigElectron_v1.cxx:85
xAOD::CaloCluster_v1::ENG_CALIB_TOT
@ ENG_CALIB_TOT
Calibration Hit energy inside the cluster.
Definition: CaloCluster_v1.h:195
xAOD::CaloCluster_v1::ENG_CALIB_FRAC_EM
@ ENG_CALIB_FRAC_EM
Calibration Hit energy inside the cluster caused by e/gamma/pi0.
Definition: CaloCluster_v1.h:248
xAOD::CaloCluster_v1::ISOLATION
@ ISOLATION
Energy weighted fraction of non-clustered perimeter cells.
Definition: CaloCluster_v1.h:146
xAOD::PFODetails::eflowRec_LAYERENERGY_TileGap3
@ eflowRec_LAYERENERGY_TileGap3
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:78
xAOD::CaloCluster_v1::ENG_CALIB_FRAC_REST
@ ENG_CALIB_FRAC_REST
Calibration Hit energy inside the cluster caused by other particles.
Definition: CaloCluster_v1.h:253
xAOD::PFODetails::eflowRec_LAYERENERGY_EME2
@ eflowRec_LAYERENERGY_EME2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:68
xAOD::PFO_v1::p4
virtual FourMom_t p4() const
The full 4-momentum of the particle.
Definition: PFO_v1.cxx:95
python.xAODType.dummy
dummy
Definition: xAODType.py:4
xAODType::ParticleFlow
@ ParticleFlow
The object is a particle-flow object.
Definition: ObjectType.h:41
xAOD::PFODetails::eflowRec_LAYERENERGY_TileGap1
@ eflowRec_LAYERENERGY_TileGap1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:76
xAOD::PFO_v1::phi
virtual double phi() const
The azimuthal angle ( ) of the particle.
Definition: PFO_v1.cxx:67
xAOD::PFODetails::eflowRec_LAYERENERGY_PreSamplerE
@ eflowRec_LAYERENERGY_PreSamplerE
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:66
xAOD::CaloCluster_v1::retrieveMoment
bool retrieveMoment(MomentType type, double &value) const
Retrieve individual moment.
Definition: CaloCluster_v1.cxx:738
xAOD::PFODetails::eflowRec_LAYERENERGY_TileBar1
@ eflowRec_LAYERENERGY_TileBar1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:74
xAOD::PFO_v1::genvecP4EM
GenVecFourMom_t genvecP4EM() const
get EM scale 4-vector
Definition: PFO_v1.cxx:160
xAOD::Vertex_v1::z
float z() const
Returns the z position.
xAOD::PFO_v1::mEM
virtual double mEM() const
get EM scale mass
Definition: PFO_v1.cxx:224
xAOD::CaloCluster_v1::FIRST_ENG_DENS
@ FIRST_ENG_DENS
First Moment in E/V.
Definition: CaloCluster_v1.h:143
xAOD::PFODetails::Track
@ Track
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:170
xAOD::CaloCluster_v1::ENG_BAD_CELLS
@ ENG_BAD_CELLS
Total em-scale energy of bad cells in this cluster.
Definition: CaloCluster_v1.h:148
xAOD::PFO_v1::isCharged
bool isCharged() const
is a charged PFO
Definition: PFO_v1.cxx:251
xAOD::PFO_v1::eEM
virtual double eEM() const
get EM scale energy
Definition: PFO_v1.cxx:233
xAOD::PFO_v1
Class describing a particle flow object.
Definition: PFO_v1.h:35
xAOD::PFODetails::eflowRec_LAYERENERGY_MINIFCAL0
@ eflowRec_LAYERENERGY_MINIFCAL0
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:85
xAOD::PFO_v1::PFO_v1
PFO_v1()
Default constructor.
Definition: PFO_v1.cxx:44
xAOD::PFODetails::eflowRec_FIRST_ENG_DENS
@ eflowRec_FIRST_ENG_DENS
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:34
xAOD::PFODetails::eflowRec_LAYERENERGY_EMB1
@ eflowRec_LAYERENERGY_EMB1
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:63
xAOD::PFO_v1::addAssociatedParticleLink
bool addAssociatedParticleLink(PFODetails::PFOParticleType ParticleType, const ElementLink< IParticleContainer > &theParticle)
Add an IParticle constituent via enum - appends to existing container.
Definition: PFO_v1.cxx:571
SG::AuxElement::makePrivateStore
void makePrivateStore()
Create a new (empty) private store for this object.
Definition: AuxElement.cxx:172
xAOD::CaloCluster_v1::AVG_TILE_Q
@ AVG_TILE_Q
Sum(E_cell_Tile^2 Q_cell_Tile)/Sum(E_cell_Tile^2)
Definition: CaloCluster_v1.h:165
xAOD::PFO_v1::etaEM
virtual double etaEM() const
get EM scale eta
Definition: PFO_v1.cxx:215
xAOD::PFODetails::eflowRec_N_BAD_CELLS
@ eflowRec_N_BAD_CELLS
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:46
ParticleGun_SamplingFraction.radius
radius
Definition: ParticleGun_SamplingFraction.py:96
xAOD::PFO_v1::m
virtual double m() const
The invariant mass of the particle.
Definition: PFO_v1.cxx:74
xAOD::PFO_v1::setAssociatedParticleLink
bool setAssociatedParticleLink(PFODetails::PFOParticleType ParticleType, const ElementLink< IParticleContainer > &theParticle)
Set an IParticle constituent via enum - does NOT append to existing container.
Definition: PFO_v1.cxx:561
xAOD::PFODetails::eflowRec_ENG_BAD_CELLS
@ eflowRec_ENG_BAD_CELLS
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:45
xAOD::PFO_v1::VertexCorrectTheFourVector
void VertexCorrectTheFourVector(const TVector3 &vertexToCorrectTo, TLorentzVector &theFourVector) const
This does the vertex correction of neutral PFO.
Definition: PFO_v1.cxx:752
xAOD::PFO_v1::setP4
void setP4(const FourMom_t &vec)
set the 4-vec
Definition: PFO_v1.cxx:107
xAOD::PFO_v1::vertex
const xAOD::Vertex * vertex() const
Retrieve a const pointer to the xAOD::Vertex a charged PFO is associated to.
Definition: PFO_v1.cxx:712
xAOD::PFODetails::eflowRec_ENG_CALIB_FRAC_HAD
@ eflowRec_ENG_CALIB_FRAC_HAD
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:96
xAOD::PFO_v1::getClusterMoment
bool getClusterMoment(float &theMoment, xAOD::CaloCluster::MomentType momentType) const
Accessor for cluster moments.
Definition: PFO_v1.cxx:402
xAOD::PFODetails::eflowRec_CELL_SIG_SAMPLING
@ eflowRec_CELL_SIG_SAMPLING
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:51
isCharged
bool isCharged(const T &p)
Definition: AtlasPID.h:545
xAOD::Vertex_v1
Class describing a Vertex.
Definition: Vertex_v1.h:42
xAOD::CaloCluster_v1::SIGNIFICANCE
@ SIGNIFICANCE
Cluster significance.
Definition: CaloCluster_v1.h:157
xAOD::PFODetails::eflowRec_LAYERENERGY_HEC3
@ eflowRec_LAYERENERGY_HEC3
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:72
xAOD::PFODetails::CaloCluster
@ CaloCluster
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:169
SG::Accessor::isAvailableWritable
bool isAvailableWritable(ELT &e) const
Test to see if this variable exists in the store and is writable.
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
xAOD::PFODetails::eflowRec_LAYERENERGY_EME3
@ eflowRec_LAYERENERGY_EME3
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:69
xAOD::Vertex_v1::y
float y() const
Returns the y position.
xAOD::PFODetails::eflowRec_LATERAL
@ eflowRec_LATERAL
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:40
xAOD::PFO_v1::isJetETMissFloatForCompression
bool isJetETMissFloatForCompression(xAOD::PFODetails::PFOAttributes AttributeType) const
Performs a check as to whether a variable should be compressed.
Definition: PFO_v1.cxx:312
xAOD::PFODetails::eflowRec_ENG_FRAC_MAX
@ eflowRec_ENG_FRAC_MAX
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:44
xAOD::PFODetails::eflowRec_ENG_CALIB_FRAC_REST
@ eflowRec_ENG_CALIB_FRAC_REST
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:97
xAOD::PFODetails::HadronicCalo
@ HadronicCalo
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:172
xAOD::CaloCluster_v1::LATERAL
@ LATERAL
Normalized lateral moment.
Definition: CaloCluster_v1.h:137
SG::ConstAccessor::isAvailable
bool isAvailable(const ELT &e) const
Test to see if this variable exists in the store.
xAOD::PFODetails::eflowRec_ENG_POS
@ eflowRec_ENG_POS
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:48
xAOD::PFO_v1::cluster
const CaloCluster * cluster(unsigned int index) const
Retrieve a const pointer to a CaloCluster.
Definition: PFO_v1.cxx:669
xAOD::JetConstituentVector::iterator
Definition: JetConstituentVector.h:121
xAOD::CaloCluster_v1::DELTA_ALPHA
@ DELTA_ALPHA
Angular shower axis deviation ( ) from IP-to-Center.
Definition: CaloCluster_v1.h:130
xAOD::CaloCluster_v1::ENG_CALIB_FRAC_HAD
@ ENG_CALIB_FRAC_HAD
Calibration Hit energy inside the cluster caused by charged pi+ and pi-.
Definition: CaloCluster_v1.h:251
xAOD::CaloCluster_v1::CELL_SIGNIFICANCE
@ CELL_SIGNIFICANCE
Cell significance = E/sig of the cell with the largest |E|/sig.
Definition: CaloCluster_v1.h:159
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
xAOD::PFODetails::eflowRec_ENG_CALIB_TOT
@ eflowRec_ENG_CALIB_TOT
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:94
xAOD::PFODetails::eflowRec_AVG_LAR_Q
@ eflowRec_AVG_LAR_Q
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:52
xAOD::PFODetails::eflowRec_LAYERENERGY_FCAL0
@ eflowRec_LAYERENERGY_FCAL0
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:82
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
xAOD::PFO_v1::rapidity
virtual double rapidity() const
The true rapidity (y) of the particle.
Definition: PFO_v1.cxx:90
xAODType::ObjectType
ObjectType
Type of objects that have a representation in the xAOD EDM.
Definition: ObjectType.h:32
xAOD::PFODetails::eflowRec_ENG_CALIB_FRAC_EM
@ eflowRec_ENG_CALIB_FRAC_EM
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:95
xAOD::PFO_v1::setAttribute
void setAttribute(PFODetails::PFOAttributes AttributeType, const T &anAttribute)
Set a PFO Variable via enum - overwrite is allowed.
xAOD::CaloCluster_v1::LONGITUDINAL
@ LONGITUDINAL
Normalized longitudinal moment.
Definition: CaloCluster_v1.h:138
xAOD::PFO_v1::GetVertexCorrectedEMFourVec
TLorentzVector GetVertexCorrectedEMFourVec(const xAOD::Vertex &vertexToCorrectTo) const
Correct EM scale 4-vector to point at a vertex.
Definition: PFO_v1.cxx:737
ParticleType
ParticleType
Definition: TruthClasses.h:8
xAOD::bool
setBGCode setTAP setLVL2ErrorBits bool
Definition: TrigDecision_v1.cxx:60
xAOD::PFO_v1::toPersistent
void toPersistent()
prepare all links for persistification
Definition: PFO_v1.cxx:774
xAOD::CaloCluster_v1::N_BAD_CELLS
@ N_BAD_CELLS
number of bad cells
Definition: CaloCluster_v1.h:149
xAOD::PFODetails::NeutralPFO
@ NeutralPFO
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:174
xAOD::PFO_v1::centerMag
float centerMag() const
get CenterMag moment needed for vertex correction
xAOD::PFODetails::eflowRec_LAYERENERGY_TileExt2
@ eflowRec_LAYERENERGY_TileExt2
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:81
xAOD::CaloCluster_v1::ENG_POS
@ ENG_POS
Total positive Energy of this cluster.
Definition: CaloCluster_v1.h:156
xAOD::PFODetails::eflowRec_LAYERENERGY_TileExt0
@ eflowRec_LAYERENERGY_TileExt0
Definition: Event/xAOD/xAODPFlow/xAODPFlow/PFODefs.h:79
xAOD::CaloCluster_v1::BADLARQ_FRAC
@ BADLARQ_FRAC
Energy fraction of LAr cells with quality larger than a given cut.
Definition: CaloCluster_v1.h:155
Generated on Thu Nov 7 2024 21:23:26 for ATLAS Offline Software by   doxygen 1.8.18