ATLAS Offline Software
Public Member Functions | Private Member Functions | Private Attributes | List of all members
eflowCaloObject Class Reference

An internal EDM object which stores information about systems of associated tracks and calorimeter clusters. More...

#include <eflowCaloObject.h>

Collaboration diagram for eflowCaloObject:

Public Member Functions

 eflowCaloObject ()=default
 
 ~eflowCaloObject ()
 
void addTrack (eflowRecTrack *track)
 
void addCluster (eflowRecCluster *cluster)
 
void addTrackClusterLinks (const std::vector< eflowTrackClusterLink * > &trackClusterLink)
 
void addTracks (const std::vector< eflowRecTrack * > &tracks)
 
void addClusters (const std::vector< eflowRecCluster * > &clusters)
 
void setTrackClusterLinkSubtractionStatus (unsigned int index, std::pair< float, float > energyRatio_energyValPair)
 
const eflowRecTrack * efRecTrack (int i) const
 
eflowRecTrack * efRecTrack (int i)
 
unsigned nTracks () const
 
void clearTracks ()
 
const eflowRecCluster * efRecCluster (int i) const
 
eflowRecCluster * efRecCluster (int i)
 
unsigned nClusters () const
 
void clearClusters ()
 
const std::vector< std::pair< eflowTrackClusterLink *, std::pair< float, float > > > & efRecLink () const
 
void clearLinks ()
 
double getExpectedEnergy () const
 
double getExpectedVariance () const
 
double getClusterEnergy () const
 
void simulateShower (eflowLayerIntegrator *integrator, const eflowEEtaBinnedParameters *binnedParameters, const PFEnergyPredictorTool *energyP, bool useLegacyEnergyBinIndexing)
 

Private Member Functions

void addTrackClusterLink (eflowTrackClusterLink *trackClusterLink)
 

Private Attributes

std::vector< eflowRecCluster * > m_eflowRecClusters
 
std::vector< std::pair< eflowTrackClusterLink *, std::pair< float, float > > > m_trackClusterLinks
 
std::vector< eflowRecTrack * > m_eflowRecTracks
 

Detailed Description

An internal EDM object which stores information about systems of associated tracks and calorimeter clusters.

Specifically it stores vectors of pointers to eflowRecTracks, eflowRecClusters and eflowTrackClusterLinks. In addition it stores links to an xAOD::CaloClusterContainer and its associated aux container. This class also calculates the expected energy deposit in the calorimeter from a track in the system, and stores that information so that clients can retrieve it. It also calculates the calorimeter cell ordering to be used in the subtraction. Both of these things are done in the simulateShower method which uses the data stored in an eflowEEtaBinnedParameters object, which is filled by e.g the eflowCellEOverP_mc12_JetETMiss tool.

Definition at line 34 of file eflowCaloObject.h.

Constructor & Destructor Documentation

◆ eflowCaloObject()

eflowCaloObject::eflowCaloObject ( )
default

◆ ~eflowCaloObject()

eflowCaloObject::~eflowCaloObject ( )
default

Member Function Documentation

◆ addCluster()

void eflowCaloObject::addCluster ( eflowRecCluster *  cluster)
inline

Definition at line 41 of file eflowCaloObject.h.

41 { m_eflowRecClusters.push_back(cluster); }

◆ addClusters()

void eflowCaloObject::addClusters ( const std::vector< eflowRecCluster * > &  clusters)

Definition at line 36 of file eflowCaloObject.cxx.

36  {
37  m_eflowRecClusters.insert(m_eflowRecClusters.end(), clusters.begin(), clusters.end());
38  }

◆ addTrack()

void eflowCaloObject::addTrack ( eflowRecTrack *  track)
inline

Definition at line 40 of file eflowCaloObject.h.

40 { m_eflowRecTracks.push_back(track); }

◆ addTrackClusterLink()

void eflowCaloObject::addTrackClusterLink ( eflowTrackClusterLink *  trackClusterLink)
inlineprivate

Definition at line 78 of file eflowCaloObject.h.

78 { m_trackClusterLinks.push_back(std::pair(trackClusterLink,std::pair(NAN,NAN))); }

◆ addTrackClusterLinks()

void eflowCaloObject::addTrackClusterLinks ( const std::vector< eflowTrackClusterLink * > &  trackClusterLink)

Definition at line 26 of file eflowCaloObject.cxx.

26  {
27  for (auto *ptr : trackClusterLink) {
28  addTrackClusterLink(ptr);
29  }
30 }

◆ addTracks()

void eflowCaloObject::addTracks ( const std::vector< eflowRecTrack * > &  tracks)

Definition at line 32 of file eflowCaloObject.cxx.

32  {
33  m_eflowRecTracks.insert(m_eflowRecTracks.end(), tracks.begin(), tracks.end());
34  }

◆ clearClusters()

void eflowCaloObject::clearClusters ( )
inline

Definition at line 61 of file eflowCaloObject.h.

61 { m_eflowRecClusters.clear(); }

◆ clearLinks()

void eflowCaloObject::clearLinks ( )
inline

Definition at line 66 of file eflowCaloObject.h.

66 { m_trackClusterLinks.clear(); }

◆ clearTracks()

void eflowCaloObject::clearTracks ( )
inline

Definition at line 55 of file eflowCaloObject.h.

55 { m_eflowRecTracks.clear(); }

◆ efRecCluster() [1/2]

eflowRecCluster* eflowCaloObject::efRecCluster ( int  i)
inline

Definition at line 59 of file eflowCaloObject.h.

59 { return m_eflowRecClusters[i]; }

◆ efRecCluster() [2/2]

const eflowRecCluster* eflowCaloObject::efRecCluster ( int  i) const
inline

Definition at line 58 of file eflowCaloObject.h.

58 { return m_eflowRecClusters[i]; }

◆ efRecLink()

const std::vector<std::pair<eflowTrackClusterLink*,std::pair<float,float> > >& eflowCaloObject::efRecLink ( ) const
inline

Definition at line 65 of file eflowCaloObject.h.

65 { return m_trackClusterLinks; }

◆ efRecTrack() [1/2]

eflowRecTrack* eflowCaloObject::efRecTrack ( int  i)
inline

Definition at line 53 of file eflowCaloObject.h.

53 { return m_eflowRecTracks[i]; }

◆ efRecTrack() [2/2]

const eflowRecTrack* eflowCaloObject::efRecTrack ( int  i) const
inline

Definition at line 52 of file eflowCaloObject.h.

52 { return m_eflowRecTracks[i]; }

◆ getClusterEnergy()

double eflowCaloObject::getClusterEnergy ( ) const

Definition at line 56 of file eflowCaloObject.cxx.

56  {
57  double clusterEnergy(0.0);
58  for(unsigned iCluster=0; iCluster < nClusters(); ++iCluster) {
59  clusterEnergy += efRecCluster(iCluster)->getCluster()->e();
60  }
61  return clusterEnergy;
62 }

◆ getExpectedEnergy()

double eflowCaloObject::getExpectedEnergy ( ) const

Definition at line 40 of file eflowCaloObject.cxx.

40  {
41  double expectedEnergy(0.0);
42  for(unsigned iTrack=0; iTrack < nTracks(); ++iTrack) {
43  if (!efRecTrack(iTrack)->isInDenseEnvironment()) expectedEnergy += efRecTrack(iTrack)->getEExpect();
44  }
45  return expectedEnergy;
46 }

◆ getExpectedVariance()

double eflowCaloObject::getExpectedVariance ( ) const

Definition at line 48 of file eflowCaloObject.cxx.

48  {
49  double expectedVariance(0.0);
50  for(unsigned iTrack=0; iTrack < nTracks(); ++iTrack) {
51  if (!efRecTrack(iTrack)->isInDenseEnvironment()) expectedVariance += efRecTrack(iTrack)->getVarEExpect();
52  }
53  return expectedVariance;
54 }

◆ nClusters()

unsigned eflowCaloObject::nClusters ( ) const
inline

Definition at line 60 of file eflowCaloObject.h.

60 { return m_eflowRecClusters.size(); }

◆ nTracks()

unsigned eflowCaloObject::nTracks ( ) const
inline

Definition at line 54 of file eflowCaloObject.h.

54 { return m_eflowRecTracks.size(); }

◆ setTrackClusterLinkSubtractionStatus()

void eflowCaloObject::setTrackClusterLinkSubtractionStatus ( unsigned int  index,
std::pair< float, float >  energyRatio_energyValPair 
)
inline

Definition at line 49 of file eflowCaloObject.h.

49 { m_trackClusterLinks[index].second = energyRatio_energyValPair; }

◆ simulateShower()

void eflowCaloObject::simulateShower ( eflowLayerIntegrator *  integrator,
const eflowEEtaBinnedParameters *  binnedParameters,
const PFEnergyPredictorTool *  energyP,
bool  useLegacyEnergyBinIndexing 
)

Definition at line 64 of file eflowCaloObject.cxx.

65  {
66 
67  for (auto *thisEfRecTrack : m_eflowRecTracks) {
68 
69  std::vector<eflowRecCluster*> matchedClusters;
70  matchedClusters.clear();
71  std::vector<eflowTrackClusterLink*> links = thisEfRecTrack->getClusterMatches();
72  std::vector<eflowTrackClusterLink*>::iterator itLink = links.begin();
73  std::vector<eflowTrackClusterLink*>::iterator endLink = links.end();
74  for (; itLink != endLink; ++itLink) {
75  matchedClusters.push_back((*itLink)->getCluster());
76  }
77 
78  double trackEM1eta = thisEfRecTrack->getTrackCaloPoints().getEM1eta();
79  /* if a track is in the forward EM (2.5 < eta < 3.2) then there is no EM1 -> need to use EM2 */
80  if(trackEM1eta<-998.) trackEM1eta = thisEfRecTrack->getTrackCaloPoints().getEM2eta();
81  /* if a track is not in the EM region (3.2 < eta < 4.0) then should use FCAL0 */
82  if(trackEM1eta<-998.) trackEM1eta = thisEfRecTrack->getTrackCaloPoints().getFCAL0eta();
83 
84  double trackE = thisEfRecTrack->getTrack()->e();
85 
86  if (!binnedParameters->binExists(trackE, trackEM1eta)) {
87  thisEfRecTrack->setHasBin(false);
88  return;
89  }
90 
91  /* Determine the LFI */
92  integrator->measureNewClus(matchedClusters, thisEfRecTrack);
93  eflowFirstIntENUM j1st = integrator->getFirstIntLayer();
94 
95  /*Save j1st info */
96  thisEfRecTrack->setLayerHED(j1st);
97 
98  /* Get parameters for j1st */
99  eflowRingSubtractionManager& cellSubtractionManager = thisEfRecTrack->getCellSubtractionManager();
100  cellSubtractionManager.getOrdering(binnedParameters, trackE, trackEM1eta, j1st,useLegacyEnergyBinIndexing);
101 
102  /* Set expected energy in the eflowRecTrack object */
103  const double expectedEnergy = energyP ? energyP->nnEnergyPrediction(thisEfRecTrack) : cellSubtractionManager.fudgeMean() * thisEfRecTrack->getTrack()->e();
104  const double expectedEnergySigma = fabs(cellSubtractionManager.fudgeStdDev() * thisEfRecTrack->getTrack()->e());
105 
106  const std::vector<eflowTrackClusterLink*>* bestClusters_015 = thisEfRecTrack->getAlternativeClusterMatches("cone_015");
107  const std::vector<eflowTrackClusterLink*>* bestClusters_02 = thisEfRecTrack->getAlternativeClusterMatches("cone_02");
108 
109  //First we calculate how much energy is found in a cone of 0.15
110  float totalE_015 = 0.0;
111 
112  //This pointer can be a nullptr if no clusters were matched to a track in dR < 0.15
113  if (bestClusters_015){
114  for (eflowTrackClusterLink* thisLink : *bestClusters_015){
115  eflowRecCluster* thisBestCluster = thisLink->getCluster();
116  if (thisBestCluster){
117  const xAOD::CaloCluster* theCluster = thisBestCluster->getCluster();
118  if (theCluster) {
119  if (theCluster->e()>0.0){
120  totalE_015 += theCluster->e();
121  }
122  }
123  }
124  }
125  }//if vector of 0.15 clusters exists
126 
127  double pull_015 = NAN;
128  if (expectedEnergySigma > 1e-6 ) pull_015 = (totalE_015-expectedEnergy)/expectedEnergySigma;
129  thisEfRecTrack->setpull15(pull_015);
130 
131  double trackPt = thisEfRecTrack->getTrack()->pt();
132  //If the looked up expected energy deposit was 0.0, then pull_015 is NAN. In that case we should not try to apply the 2D cut described below.
133  if (!std::isnan(pull_015)){
134  //We use a 2D cut in the pull_015 and log10 of track pt plane to define a dense environment - if too dense then we disable the charged shower subtraction
135  //The cut values are such that above a track pt of 40 GeV we are very agrressive about disabling the cell by cell charged shower subtraction.
136  //The specific cut values were found by optimising the jet resolutions in the run 2 data taking environment.
137  //Further details can be found in these slides from C. Young:
138  //https://indico.cern.ch/event/396923/contributions/949333/attachments/795878/1090871/Presentation.pdf
139  if (pull_015 > 0.0 + (log10(40000)-log10(trackPt))*33.2 && 0.0 != expectedEnergySigma && bestClusters_015 && bestClusters_02){
140  thisEfRecTrack->setSubtracted(); //this tricks eflowRec into thinking this track was subtracted, and hence no further subtraction will be done
141  thisEfRecTrack->setIsInDenseEnvironment();
142  //recalculate the LHED and the ordering and find the new expected E + sigma of expected E (the new LHED can change the latter two values we find in the look up tables)
143  //we use a larger cone of 0.2 for this
144  std::vector<eflowRecCluster*> theBestEfRecClusters_02;
145  for (eflowTrackClusterLink* thisLink : *bestClusters_02) if (thisLink->getCluster()->getCluster()->e() > 0.0) theBestEfRecClusters_02.push_back(thisLink->getCluster());
146  integrator->measureNewClus(theBestEfRecClusters_02, thisEfRecTrack);
147  j1st = integrator->getFirstIntLayer();
148  cellSubtractionManager.getOrdering(binnedParameters, trackE, trackEM1eta, j1st,useLegacyEnergyBinIndexing);
149  thisEfRecTrack->setEExpect(cellSubtractionManager.fudgeMean() * trackE, fabs(cellSubtractionManager.fudgeStdDev()*trackE)*fabs(cellSubtractionManager.fudgeStdDev()*trackE));
150  }
151  else {
152  thisEfRecTrack->setEExpect(expectedEnergy, expectedEnergySigma*expectedEnergySigma);
153  }//ok to do subtraction, and so we just set the usual expected E + sigma of expected E needed for subtraction
154  }
155  }
156 }

Member Data Documentation

◆ m_eflowRecClusters

std::vector<eflowRecCluster*> eflowCaloObject::m_eflowRecClusters
private

Definition at line 82 of file eflowCaloObject.h.

◆ m_eflowRecTracks

std::vector<eflowRecTrack*> eflowCaloObject::m_eflowRecTracks
private

Definition at line 92 of file eflowCaloObject.h.

◆ m_trackClusterLinks

std::vector<std::pair<eflowTrackClusterLink*,std::pair<float,float> > > eflowCaloObject::m_trackClusterLinks
private

Definition at line 89 of file eflowCaloObject.h.

The documentation for this class was generated from the following files:
xAOD::iterator
JetConstituentVector::iterator iterator
Definition: JetConstituentVector.cxx:68
eflowRecCluster
This class extends the information about a xAOD::CaloCluster.
Definition: eflowRecCluster.h:40
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
eflowCaloObject::nClusters
unsigned nClusters() const
Definition: eflowCaloObject.h:60
eflowRecTrack::isInDenseEnvironment
bool isInDenseEnvironment() const
Definition: eflowRecTrack.h:103
eflowRingSubtractionManager::fudgeStdDev
double fudgeStdDev() const
Definition: eflowRingSubtractionManager.h:38
dbg::ptr
void * ptr(T *p)
Definition: SGImplSvc.cxx:74
eflowLayerIntegrator::measureNewClus
void measureNewClus(const xAOD::CaloCluster *clus, const eflowTrackCaloPoints &trackCalo)
Definition: eflowLayerIntegrator.cxx:154
eflowCaloObject::addTrackClusterLink
void addTrackClusterLink(eflowTrackClusterLink *trackClusterLink)
Definition: eflowCaloObject.h:78
xAOD::CaloCluster_v1
Description of a calorimeter cluster.
Definition: CaloCluster_v1.h:59
eflowCaloObject::efRecTrack
const eflowRecTrack * efRecTrack(int i) const
Definition: eflowCaloObject.h:52
lumiFormat.i
int i
Definition: lumiFormat.py:85
DMTest::links
links
Definition: CLinks_v1.cxx:22
eflowRecTrack::getEExpect
double getEExpect() const
Definition: eflowRecTrack.h:84
eflowRingSubtractionManager::fudgeMean
double fudgeMean() const
Definition: eflowRingSubtractionManager.h:37
eflowRecCluster::getCluster
xAOD::CaloCluster * getCluster()
Definition: eflowRecCluster.h:49
eflowCaloObject::nTracks
unsigned nTracks() const
Definition: eflowCaloObject.h:54
eflowCaloObject::m_eflowRecClusters
std::vector< eflowRecCluster * > m_eflowRecClusters
Definition: eflowCaloObject.h:82
TauGNNUtils::Variables::Track::trackPt
bool trackPt(const xAOD::TauJet &, const xAOD::TauTrack &track, double &out)
Definition: TauGNNUtils.cxx:472
eflowRingSubtractionManager
This stores information, a rank and ring thickness, about cell rings in an ordered way.
Definition: eflowRingSubtractionManager.h:31
eflowRingSubtractionManager::getOrdering
bool getOrdering(const eflowEEtaBinnedParameters *binnedParameters, double e, double eta, eflowFirstIntENUM j1st, bool useLegacyEnergyBinIndexing)
Definition: eflowRingSubtractionManager.cxx:103
eflowCaloObject::m_trackClusterLinks
std::vector< std::pair< eflowTrackClusterLink *, std::pair< float, float > > > m_trackClusterLinks
Definition: eflowCaloObject.h:89
PFEnergyPredictorTool::nnEnergyPrediction
float nnEnergyPrediction(const eflowRecTrack *ptr) const
Definition: PFEnergyPredictorTool.cxx:134
DeMoScan.index
string index
Definition: DeMoScan.py:364
eflowCaloObject::efRecCluster
const eflowRecCluster * efRecCluster(int i) const
Definition: eflowCaloObject.h:58
eflowRecTrack::getVarEExpect
double getVarEExpect() const
Definition: eflowRecTrack.h:85
RunTileMonitoring.clusters
clusters
Definition: RunTileMonitoring.py:133
eflowLayerIntegrator::getFirstIntLayer
eflowFirstIntENUM getFirstIntLayer() const
Definition: eflowLayerIntegrator.cxx:230
eflowEEtaBinBase::binExists
bool binExists(double e, double eta) const
Definition: eflowEEtaBinBase.h:32
eflowCaloObject::m_eflowRecTracks
std::vector< eflowRecTrack * > m_eflowRecTracks
Definition: eflowCaloObject.h:92
xAOD::track
@ track
Definition: TrackingPrimitives.h:512
xAOD::CaloCluster_v1::e
virtual double e() const
The total energy of the particle.
Definition: CaloCluster_v1.cxx:265
eflowFirstIntRegions::J1STLAYER
J1STLAYER
Definition: eflowCaloRegions.h:58
Generated on Sun Dec 22 2024 21:31:38 for ATLAS Offline Software by   doxygen 1.8.18