ATLAS Offline Software
Classes | Public Member Functions | Private Member Functions | Private Attributes | Friends | List of all members
DerivationFramework::EGElectronAmbiguityTool Class Reference

#include <EGElectronAmbiguityTool.h>

Inheritance diagram for DerivationFramework::EGElectronAmbiguityTool:
Collaboration diagram for DerivationFramework::EGElectronAmbiguityTool:

Classes

struct  DecorHandles
 

Public Member Functions

 EGElectronAmbiguityTool (const std::string &t, const std::string &n, const IInterface *p)
 
virtual StatusCode initialize () override final
 
virtual StatusCode addBranches (const EventContext &ctx) const override final
 

Private Member Functions

StatusCode decorateSimple (DecorHandles &dh, std::unique_ptr< ConstDataVector< xAOD::TrackParticleContainer >> &tpC, const xAOD::Electron *ele, const xAOD::Vertex *pvtx) const
 

Private Attributes

SG::ReadHandleKey< xAOD::ElectronContainerm_containerName
 
SG::ReadHandleKey< xAOD::VertexContainerm_VtxContainerName
 
SG::ReadHandleKey< xAOD::TrackParticleContainerm_tpContainerName
 
SG::ReadHandleKey< xAOD::TrackParticleContainerm_tpCName
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_drv { this, "DFCommonSimpleConvRadius", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dphiv { this, "DFCommonSimpleConvPhi", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dmee { this, "DFCommonSimpleMee", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dmeeVtx { this, "DFCommonSimpleMeeAtVtx", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dsep { this, "DFCommonSimpleSeparation", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dambi { this, "DFCommonAddAmbiguity", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dtrv { this, "DFCommonProdTrueRadius", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dtpv { this, "DFCommonProdTruePhi", "", "" }
 
SG::WriteDecorHandleKey< xAOD::ElectronContainerm_dtzv { this, "DFCommonProdTrueZ", "", "" }
 
bool m_isMC {}
 
double m_elepTCut
 
std::string m_idCut
 
unsigned int m_nSiCut
 
double m_dctCut
 
double m_sepCut
 
double m_dzCut
 
double m_rvECCut
 
double m_meeAtVtxECCut
 
double m_meeICCut
 

Friends

struct DecorHandles
 

Detailed Description

Definition at line 33 of file EGElectronAmbiguityTool.h.

Constructor & Destructor Documentation

◆ EGElectronAmbiguityTool()

DerivationFramework::EGElectronAmbiguityTool::EGElectronAmbiguityTool ( const std::string &  t,
const std::string &  n,
const IInterface *  p 
)

Definition at line 52 of file EGElectronAmbiguityTool.cxx.

55  : base_class(t, n, p)
56 {
57 
58  declareProperty("isMC", m_isMC);
59 
60  declareProperty(
61  "pTCut", m_elepTCut = 9e3, "minimum pT for an electron to be studied");
62  declareProperty("idCut",
63  m_idCut = "DFCommonElectronsLHLoose",
64  "minimal quality for an electron to be studied");
65 
66  declareProperty(
67  "nSiCut", m_nSiCut = 7, "minimum number of Si hits in the other track");
68  declareProperty(
69  "dzsinTCut", m_dzCut = 0.5, "max dz sinTheta between ele and other tracks");
70  declareProperty("SeparationCut", m_sepCut = 1., "first separation cut");
71  declareProperty("DCTCut", m_dctCut = 0.02, "second separation cut");
72 
73  declareProperty("radiusCut",
74  m_rvECCut = 20,
75  "minimum radius to be classified as external conversion");
76  declareProperty(
77  "meeAtVtxCut",
78  m_meeAtVtxECCut = 100,
79  "maximal mass at vertex to be classified as external conversion");
80  declareProperty("meeCut",
81  m_meeICCut = 100,
82  "maximal mass at primary vertex to be classified as gamma*");
83 }

Member Function Documentation

◆ addBranches()

StatusCode DerivationFramework::EGElectronAmbiguityTool::addBranches ( const EventContext &  ctx) const
finaloverridevirtual

Definition at line 133 of file EGElectronAmbiguityTool.cxx.

134 {
135 
136  DecorHandles dh (*this, ctx);
137 
138  static const SG::AuxElement::ConstAccessor<char> aidCut(m_idCut);
139 
140  // retrieve primary vertex
141  const xAOD::Vertex* pvtx(nullptr);
143  if (inputVtx.ptr()){
144  const xAOD::VertexContainer* vtxC = inputVtx.ptr();
145  for (const auto* vertex : *vtxC) {
146  if (vertex->vertexType() == xAOD::VxType::VertexType::PriVtx) {
147  pvtx = vertex;
148  break;
149  }
150  }
151  }
152 
153  // retrieve electron container
155  ctx };
156  const xAOD::ElectronContainer* eleC = inputElectrons.ptr();
157 
158  if (!eleC) {
160  "Couldn't retrieve Electron container with key: " << m_containerName);
161  return StatusCode::FAILURE;
162  }
163  if (!pvtx) {
164  ATH_MSG_DEBUG("No primary vertex found. Setting default values.");
165  for (const xAOD::Electron* iele : *eleC) {
166  dh.drv(*iele) = -1;
167  dh.dphiv(*iele) = -1;
168  dh.dmee(*iele) = -1;
169  dh.dmeeVtx(*iele) = -1;
170  dh.dsep(*iele) = -1;
171  dh.dambi(*iele) = -1;
172  dh.dtrv(*iele) = -1;
173  dh.dtpv(*iele) = -1;
174  dh.dtzv(*iele) = -1;
175  }
176  return StatusCode::SUCCESS;
177  }
178  ATH_MSG_DEBUG("Pvx z = " << pvtx->z() << ", number of electrons "
179  << eleC->size());
180 
181  // Make a container of selected tracks : with Si hits, close to electron track
183  m_tpContainerName, ctx
184  };
185  const xAOD::TrackParticleContainer* idtpC = inputTrackParticles.ptr();
186 
187  std::set<const xAOD::TrackParticle*> alreadyStored;
188  std::set<const xAOD::TrackParticle*> eleIDtpStored, eleGSFtpStored;
189  auto closeByTracks =
190  std::make_unique<ConstDataVector<xAOD::TrackParticleContainer>>(
192 
193  for (const auto* ele : *eleC) {
194 
195  dh.dambi(*ele) = -1;
196 
197  // Electron preselection
198  if (ele->pt() < m_elepTCut ||
199  m_idCut.empty() || !aidCut.isAvailable(*ele) || !aidCut(*ele))
200  continue;
201 
202  // Just for debug
203  const xAOD::TrackParticle* eleGSFtp = ele->trackParticle();
204  eleGSFtpStored.insert(eleGSFtp);
205 
206  const xAOD::TrackParticle* eleIDtp =
208  eleIDtpStored.insert(eleIDtp);
209 
210  // The loop on track
211  for (const auto* tp : *idtpC) {
212 
213  // Keep the electron track (I build a container to run vertexing on it...)
214  if (tp == eleIDtp) {
215  closeByTracks->push_back(tp);
216  alreadyStored.insert(tp);
217  continue;
218  }
219 
220  // potential candidate to store if not already there
221  if (alreadyStored.find(tp) != alreadyStored.end())
222  continue;
223 
225  // if (tp->charge() * ele->charge() > 0)
226  // continue;
227 
228  // Close-by
229  double dR = eleIDtp->p4().DeltaR(tp->p4());
230  double dz = std::abs(eleIDtp->z0() - tp->z0()) * sin(eleIDtp->theta());
231  if (dR >= 0.3 || dz >= m_dzCut)
232  continue;
233 
234  // With minimum number of Si hits
236  continue;
237 
238  alreadyStored.insert(tp);
239 
240  closeByTracks->push_back(tp);
241  }
242  }
243 
244  if (closeByTracks->empty())
245  return StatusCode::SUCCESS;
246 
247  if (msgLvl(MSG::DEBUG)) {
249  ctx };
250  const xAOD::TrackParticleContainer* tpC = tpCReadHandle.ptr();
251 
252  ATH_MSG_DEBUG("Number of input tracks "
253  << idtpC->size() << " , number of selected close-by tracks "
254  << closeByTracks->size() << " , number of GSF tracks "
255  << tpC->size());
256  for (const auto* trk : eleIDtpStored)
257  ATH_MSG_DEBUG("ele ID trk "
258  << trk << " pt = " << trk->pt() * 1e-3
259  << " eta = " << trk->eta() << " phi = " << trk->phi()
260  << " nSi = " << xAOD::EgammaHelpers::numberOfSiHits(trk));
261  for (const auto* trk : eleGSFtpStored)
262  ATH_MSG_DEBUG("ele GSF trk "
263  << trk << " pt = " << trk->pt() * 1e-3
264  << " eta = " << trk->eta() << " phi = " << trk->phi()
265  << " nSi = " << xAOD::EgammaHelpers::numberOfSiHits(trk));
266  for (const xAOD::TrackParticle* trk : *closeByTracks)
267  ATH_MSG_DEBUG("closeby trk "
268  << trk << " pt = " << trk->pt() * 1e-3
269  << " eta = " << trk->eta() << " phi = " << trk->phi()
270  << " nSi = " << xAOD::EgammaHelpers::numberOfSiHits(trk));
271  }
272 
273  for (const auto* ele : *eleC) {
274 
275  // Electron preselection
276  if (ele->pt() < m_elepTCut ||
277  m_idCut.empty() || !aidCut.isAvailable(*ele) || !aidCut(*ele))
278  continue;
279 
280  // Henri's circles
281  if (decorateSimple(dh, closeByTracks, ele, pvtx).isFailure()) {
282  ATH_MSG_ERROR("Cannot decorate the electron with the simple info");
283  return StatusCode::FAILURE;
284  }
285 
286  } // loop on electrons to decorate
287 
288  return StatusCode::SUCCESS;
289 }

◆ decorateSimple()

StatusCode DerivationFramework::EGElectronAmbiguityTool::decorateSimple ( DecorHandles dh,
std::unique_ptr< ConstDataVector< xAOD::TrackParticleContainer >> &  tpC,
const xAOD::Electron ele,
const xAOD::Vertex pvtx 
) const
private

HelixToCircle Main Track Electron

HelixToCircle Other Electron Conv Track

Definition at line 293 of file EGElectronAmbiguityTool.cxx.

298 {
299  // This is the GSF electron track
300  const xAOD::TrackParticle* eleGSFtrkP = ele->trackParticle();
301 
302  // And the ID one
303  const xAOD::TrackParticle* eleIDtrkP =
305 
306  // For the time being, use the ID one, to be consistent when we only find a
307  // good ID to make a conversion and no GSF
308  bool useGSF = false; // hardcoded because it seems we will not use true. Kept
309  // for the time being, as not 100% sure
310  const xAOD::TrackParticle* eletrkP = useGSF ? eleGSFtrkP : eleIDtrkP;
311 
312  ATH_MSG_DEBUG("Electron pt = " << ele->pt() * 1e-3 << " eta = " << ele->eta()
313  << " phi = " << ele->phi() << " GSF trk ptr = "
314  << eleGSFtrkP << " ID trk ptr " << eleIDtrkP);
315 
316  if (m_isMC) {
317  const xAOD::TruthParticle* truthEl =
319  double tpvr = -1, tpvp = 9e9, tpvz = 9e9;
320  if (truthEl && MC::isElectron(truthEl) &&
321  truthEl->prodVtx() != nullptr) {
322  tpvr = truthEl->prodVtx()->perp();
323  tpvp = truthEl->prodVtx()->phi();
324  tpvz = truthEl->prodVtx()->z();
325  }
326  dh.dtrv(*ele) = tpvr;
327  dh.dtpv(*ele) = tpvp;
328  dh.dtzv(*ele) = tpvz;
329  }
330 
331  // Find the closest track particle with opposite charge and a minimum nb of Si
332  // hits
333  const xAOD::TrackParticle* otrkP(nullptr);
334  double detaMin = 9e9;
335  for (const xAOD::TrackParticle* tp : *tpC) {
336  // Keep only opposite charge
337  if (tp->charge() * eletrkP->charge() > 0)
338  continue;
339 
340  // Close-by
341  double dR = eletrkP->p4().DeltaR(tp->p4());
342  double dz = std::abs(eletrkP->z0() - tp->z0()) * sin(eletrkP->theta());
343  if (dR >= 0.3 || dz >= m_dzCut)
344  continue;
345 
346  double deta = std::abs(eletrkP->eta() - tp->eta());
347  if (deta < detaMin) {
348  otrkP = tp;
349  detaMin = deta;
350  }
351  }
352 
353  double rv = -9e9;
354  double pv = -9e9;
355  double mee = -1.;
356  double meeAtVtx = -1.;
357  double sep = -9e9;
358  bool goodConv = false;
359 
360  if (otrkP) {
361 
362  // To be consistent with the other, use the ID track.
363  TLorentzVector ep4;
364  ep4.SetPtEtaPhiM(eletrkP->pt(), eletrkP->eta(), eletrkP->phi(), ParticleConstants::electronMassInMeV);
365 
366  // Maybe could see if a GSF tp exists for this ID tp and use it if yes ?
367  TLorentzVector op4;
368  op4.SetPtEtaPhiM(otrkP->pt(), otrkP->eta(), otrkP->phi(), ParticleConstants::electronMassInMeV);
369 
370  // Simple masses
371  mee = (ep4 + op4).M();
372  op4.SetPhi(eletrkP->phi());
373  meeAtVtx = (ep4 + op4).M();
374 
375  // And the conversion point
376  std::vector<double> helix1, helix2;
377  helix1.resize(5);
378  helix2.resize(5);
379  helix(eletrkP, pvtx, helix1);
380  helix(otrkP, pvtx, helix2);
381 
382  double beta(0.);
383  if (helix1[4] < helix2[4])
384  beta = TMath::PiOver2() - helix1[4];
385  else
386  beta = TMath::PiOver2() - helix2[4];
387 
388  double phi1(helix1[4] + beta);
389  if (phi1 > TMath::TwoPi())
390  phi1 -= TMath::TwoPi();
391  if (phi1 < 0.)
392  phi1 += TMath::TwoPi();
393 
394  double phi2(helix2[4] + beta);
395  if (phi2 > TMath::TwoPi())
396  phi2 -= TMath::TwoPi();
397  if (phi2 < 0.)
398  phi2 += TMath::TwoPi();
399 
401  double r1 = 1 / (2. * std::abs(helix1[1]));
402 
403  double charge1(1.);
404  if (helix1[1] < 0.)
405  charge1 = -1.;
406  double rcenter1(helix1[3] / charge1 + r1);
407  double phicenter1(phi1 + TMath::PiOver2() * charge1);
408 
409  double x1 = rcenter1 * cos(phicenter1);
410  double y1 = rcenter1 * sin(phicenter1);
411 
413  double r2 = 1 / (2. * std::abs(helix2[1]));
414 
415  double charge2(1.);
416  if (helix2[1] < 0.)
417  charge2 = -1.;
418  double rcenter2(helix2[3] / charge2 + r2);
419  double phicenter2(phi2 + TMath::PiOver2() * charge2);
420 
421  double x2 = rcenter2 * cos(phicenter2);
422  double y2 = rcenter2 * sin(phicenter2);
424 
425  double dx(x1 - x2);
426  if (dx < 1e-9 && dx > 0.)
427  dx = 1e-9;
428  if (dx > -1e-9 && dx < 0.)
429  dx = -1e-9;
430  double slope((y1 - y2) / dx);
431  double b(y1 - slope * x1);
432  double alpha(atan(slope));
433  double d(sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2)));
434  // only keeping opposite sign option
435  double separation = d - r1 - r2;
436  double cpx1, cpx2;
437  if (x1 > x2) {
438  cpx1 = x1 - r1 * cos(alpha);
439  cpx2 = x2 + r2 * cos(alpha);
440  } else {
441  cpx1 = x1 + r1 * cos(alpha);
442  cpx2 = x2 - r2 * cos(alpha);
443  }
444 
445  double temp1 = (cpx1 + cpx2) / 2;
446  double temp2 = slope * temp1 + b;
447  double convX = cos(beta) * temp1 + sin(beta) * temp2;
448  double convY = -sin(beta) * temp1 + cos(beta) * temp2;
449 
450  double dct(helix1[0] - helix2[0]);
451 
453  if (std::abs(separation) < m_sepCut && std::abs(dct) < m_dctCut) {
454  goodConv = true;
455  sep = separation;
456  pv = std::atan2(convY, convX);
457  rv = sqrt(convX * convX + convY * convY);
458  if (convX * cos(eletrkP->phi()) + convY * sin(eletrkP->phi()) < 0)
459  rv *= -1.;
460  }
461  } else {
462  dh.dambi(*ele) = -1;
463  }
464  dh.drv(*ele) = rv;
465  dh.dphiv(*ele) = pv;
466  dh.dmee(*ele) = mee;
467  dh.dmeeVtx(*ele) = meeAtVtx;
468  dh.dsep(*ele) = sep;
469  if (goodConv && rv > m_rvECCut && meeAtVtx < m_meeAtVtxECCut)
470  dh.dambi(*ele) = 2;
471  else if (otrkP) {
472  if (mee < m_meeICCut)
473  dh.dambi(*ele) = 1;
474  else
475  dh.dambi(*ele) = 0;
476  }
477  return StatusCode::SUCCESS;
478 }

◆ initialize()

StatusCode DerivationFramework::EGElectronAmbiguityTool::initialize ( )
finaloverridevirtual

Definition at line 86 of file EGElectronAmbiguityTool.cxx.

87 {
88 
90  ATH_CHECK(m_VtxContainerName.initialize());
93 
94  const std::string baseName = m_containerName.key();
95  m_drv = baseName + ".DFCommonSimpleConvRadius";
96  m_dphiv = baseName + ".DFCommonSimpleConvPhi";
97  m_dmee = baseName + ".DFCommonSimpleMee";
98  m_dmeeVtx = baseName + ".DFCommonSimpleMeeAtVtx";
99  m_dsep = baseName + ".DFCommonSimpleSeparation";
100  m_dambi = baseName + ".DFCommonAddAmbiguity";
101  m_dtrv = baseName + ".DFCommonProdTrueRadius";
102  m_dtpv = baseName + ".DFCommonProdTruePhi";
103  m_dtzv = baseName + ".DFCommonProdTrueZ";
104 
114 
115  return StatusCode::SUCCESS;
116 }

Friends And Related Function Documentation

◆ DecorHandles

friend struct DecorHandles
friend

Definition at line 111 of file EGElectronAmbiguityTool.h.

Member Data Documentation

◆ m_containerName

SG::ReadHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_containerName
private
Initial value:
{
this,
"ContainerName",
"Electrons",
"SG key of electron container"
}

Definition at line 45 of file EGElectronAmbiguityTool.h.

◆ m_dambi

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dambi { this, "DFCommonAddAmbiguity", "", "" }
private

Definition at line 88 of file EGElectronAmbiguityTool.h.

◆ m_dctCut

double DerivationFramework::EGElectronAmbiguityTool::m_dctCut
private

Definition at line 127 of file EGElectronAmbiguityTool.h.

◆ m_dmee

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dmee { this, "DFCommonSimpleMee", "", "" }
private

Definition at line 79 of file EGElectronAmbiguityTool.h.

◆ m_dmeeVtx

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dmeeVtx { this, "DFCommonSimpleMeeAtVtx", "", "" }
private

Definition at line 82 of file EGElectronAmbiguityTool.h.

◆ m_dphiv

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dphiv { this, "DFCommonSimpleConvPhi", "", "" }
private

Definition at line 76 of file EGElectronAmbiguityTool.h.

◆ m_drv

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_drv { this, "DFCommonSimpleConvRadius", "", "" }
private

Definition at line 73 of file EGElectronAmbiguityTool.h.

◆ m_dsep

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dsep { this, "DFCommonSimpleSeparation", "", "" }
private

Definition at line 85 of file EGElectronAmbiguityTool.h.

◆ m_dtpv

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dtpv { this, "DFCommonProdTruePhi", "", "" }
private

Definition at line 94 of file EGElectronAmbiguityTool.h.

◆ m_dtrv

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dtrv { this, "DFCommonProdTrueRadius", "", "" }
private

Definition at line 91 of file EGElectronAmbiguityTool.h.

◆ m_dtzv

SG::WriteDecorHandleKey<xAOD::ElectronContainer> DerivationFramework::EGElectronAmbiguityTool::m_dtzv { this, "DFCommonProdTrueZ", "", "" }
private

Definition at line 97 of file EGElectronAmbiguityTool.h.

◆ m_dzCut

double DerivationFramework::EGElectronAmbiguityTool::m_dzCut
private

Definition at line 127 of file EGElectronAmbiguityTool.h.

◆ m_elepTCut

double DerivationFramework::EGElectronAmbiguityTool::m_elepTCut
private

Definition at line 122 of file EGElectronAmbiguityTool.h.

◆ m_idCut

std::string DerivationFramework::EGElectronAmbiguityTool::m_idCut
private

Definition at line 123 of file EGElectronAmbiguityTool.h.

◆ m_isMC

bool DerivationFramework::EGElectronAmbiguityTool::m_isMC {}
private

Definition at line 119 of file EGElectronAmbiguityTool.h.

◆ m_meeAtVtxECCut

double DerivationFramework::EGElectronAmbiguityTool::m_meeAtVtxECCut
private

Definition at line 132 of file EGElectronAmbiguityTool.h.

◆ m_meeICCut

double DerivationFramework::EGElectronAmbiguityTool::m_meeICCut
private

Definition at line 132 of file EGElectronAmbiguityTool.h.

◆ m_nSiCut

unsigned int DerivationFramework::EGElectronAmbiguityTool::m_nSiCut
private

Definition at line 126 of file EGElectronAmbiguityTool.h.

◆ m_rvECCut

double DerivationFramework::EGElectronAmbiguityTool::m_rvECCut
private

Definition at line 132 of file EGElectronAmbiguityTool.h.

◆ m_sepCut

double DerivationFramework::EGElectronAmbiguityTool::m_sepCut
private

Definition at line 127 of file EGElectronAmbiguityTool.h.

◆ m_tpCName

SG::ReadHandleKey<xAOD::TrackParticleContainer> DerivationFramework::EGElectronAmbiguityTool::m_tpCName
private
Initial value:
{
this,
"tpCName",
"GSFTrackParticles",
"SG key of TrackParticleInputContainer"
}

Definition at line 64 of file EGElectronAmbiguityTool.h.

◆ m_tpContainerName

SG::ReadHandleKey<xAOD::TrackParticleContainer> DerivationFramework::EGElectronAmbiguityTool::m_tpContainerName
private
Initial value:
{
this,
"tpContainerName",
"InDetTrackParticles",
"SG key of track particles container"
}

Definition at line 57 of file EGElectronAmbiguityTool.h.

◆ m_VtxContainerName

SG::ReadHandleKey<xAOD::VertexContainer> DerivationFramework::EGElectronAmbiguityTool::m_VtxContainerName
private
Initial value:
{
this,
"VtxContainerName",
"PrimaryVertices",
"SG key of vertex container"
}

Definition at line 51 of file EGElectronAmbiguityTool.h.


The documentation for this class was generated from the following files:
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition: TrackParticle_v1.cxx:74
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
plotBeamSpotCompare.x1
x1
Definition: plotBeamSpotCompare.py:215
add-xsec-uncert-quadrature-N.alpha
alpha
Definition: add-xsec-uncert-quadrature-N.py:110
DerivationFramework::EGElectronAmbiguityTool::m_dtrv
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dtrv
Definition: EGElectronAmbiguityTool.h:91
xAOD::TruthVertex_v1::phi
float phi() const
Vertex azimuthal angle.
Definition: TruthVertex_v1.cxx:177
DerivationFramework::EGElectronAmbiguityTool::m_dmee
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dmee
Definition: EGElectronAmbiguityTool.h:79
DerivationFramework::EGElectronAmbiguityTool::m_dambi
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dambi
Definition: EGElectronAmbiguityTool.h:88
DerivationFramework::EGElectronAmbiguityTool::m_isMC
bool m_isMC
Definition: EGElectronAmbiguityTool.h:119
DerivationFramework::EGElectronAmbiguityTool::m_tpContainerName
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_tpContainerName
Definition: EGElectronAmbiguityTool.h:57
SG::VIEW_ELEMENTS
@ VIEW_ELEMENTS
this data object is a view, it does not own its elmts
Definition: OwnershipPolicy.h:18
xAOD::TrackParticle_v1::charge
float charge() const
Returns the charge.
Definition: TrackParticle_v1.cxx:151
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:67
hist_file_dump.d
d
Definition: hist_file_dump.py:142
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: TrackParticle_v1.cxx:78
plotBeamSpotCompare.x2
x2
Definition: plotBeamSpotCompare.py:217
DerivationFramework::EGElectronAmbiguityTool::m_drv
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_drv
Definition: EGElectronAmbiguityTool.h:73
DerivationFramework::EGElectronAmbiguityTool::m_dctCut
double m_dctCut
Definition: EGElectronAmbiguityTool.h:127
xAOD::TrackParticle_v1::z0
float z0() const
Returns the parameter.
DerivationFramework::EGElectronAmbiguityTool::m_rvECCut
double m_rvECCut
Definition: EGElectronAmbiguityTool.h:132
ParticleTest.tp
tp
Definition: ParticleTest.py:25
DerivationFramework::EGElectronAmbiguityTool::m_tpCName
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_tpCName
Definition: EGElectronAmbiguityTool.h:64
xAOD::Electron_v1::trackParticle
const xAOD::TrackParticle * trackParticle(size_t index=0) const
Pointer to the xAOD::TrackParticle/s that match the electron candidate.
Definition: Electron_v1.cxx:55
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:55
PyPoolBrowser.dh
dh
Definition: PyPoolBrowser.py:102
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
SG::VarHandleKey::key
const std::string & key() const
Return the StoreGate ID for the referenced object.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:141
MCP::ScaleSmearParam::r2
@ r2
DerivationFramework::EGElectronAmbiguityTool::m_sepCut
double m_sepCut
Definition: EGElectronAmbiguityTool.h:127
DerivationFramework::EGElectronAmbiguityTool::DecorHandles
friend struct DecorHandles
Definition: EGElectronAmbiguityTool.h:111
DerivationFramework::EGElectronAmbiguityTool::m_elepTCut
double m_elepTCut
Definition: EGElectronAmbiguityTool.h:122
makeTRTBarrelCans.y1
tuple y1
Definition: makeTRTBarrelCans.py:15
drawFromPickle.atan
atan
Definition: drawFromPickle.py:36
xAOD::TrackParticle_v1::p4
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
Definition: TrackParticle_v1.cxx:130
DerivationFramework::EGElectronAmbiguityTool::m_dsep
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dsep
Definition: EGElectronAmbiguityTool.h:85
DerivationFramework::EGElectronAmbiguityTool::m_VtxContainerName
SG::ReadHandleKey< xAOD::VertexContainer > m_VtxContainerName
Definition: EGElectronAmbiguityTool.h:51
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
CheckAppliedSFs.e3
e3
Definition: CheckAppliedSFs.py:264
DerivationFramework::EGElectronAmbiguityTool::m_nSiCut
unsigned int m_nSiCut
Definition: EGElectronAmbiguityTool.h:126
DerivationFramework::EGElectronAmbiguityTool::m_dzCut
double m_dzCut
Definition: EGElectronAmbiguityTool.h:127
beamspotman.n
n
Definition: beamspotman.py:727
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
xAOD::EgammaHelpers::isElectron
bool isElectron(const xAOD::Egamma *eg)
is the object an electron (not Fwd)
Definition: EgammaxAODHelpers.cxx:12
makeTRTBarrelCans.y2
tuple y2
Definition: makeTRTBarrelCans.py:18
xAOD::TruthParticle_v1
Class describing a truth particle in the MC record.
Definition: TruthParticle_v1.h:37
xAOD::TruthVertex_v1::perp
float perp() const
Vertex transverse distance from the beam line.
Definition: TruthVertex_v1.cxx:164
xAOD::VxType::PriVtx
@ PriVtx
Primary vertex.
Definition: TrackingPrimitives.h:572
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
DerivationFramework::EGElectronAmbiguityTool::m_idCut
std::string m_idCut
Definition: EGElectronAmbiguityTool.h:123
xAOD::Egamma_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
Definition: Egamma_v1.cxx:75
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
DataVector
Derived DataVector<T>.
Definition: DataVector.h:795
xAOD::TruthParticle_v1::prodVtx
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
Definition: TruthParticle_v1.cxx:75
DerivationFramework::EGElectronAmbiguityTool::m_dmeeVtx
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dmeeVtx
Definition: EGElectronAmbiguityTool.h:82
DerivationFramework::EGElectronAmbiguityTool::m_dphiv
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dphiv
Definition: EGElectronAmbiguityTool.h:76
grepfile.sep
sep
Definition: grepfile.py:38
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:76
xAOD::TruthHelpers::getTruthParticle
const xAOD::TruthParticle * getTruthParticle(const xAOD::IParticle &p)
Return the truthParticle associated to the given IParticle (if any)
Definition: xAODTruthHelpers.cxx:25
xAOD::EgammaHelpers::numberOfSiHits
std::size_t numberOfSiHits(const xAOD::TrackParticle *tp)
return the number of Si hits in the track particle
Definition: ElectronxAODHelpers.cxx:66
xAOD::Electron_v1
Definition: Electron_v1.h:34
dq_make_web_display.rv
def rv
Definition: dq_make_web_display.py:218
DerivationFramework::EGElectronAmbiguityTool::m_dtpv
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dtpv
Definition: EGElectronAmbiguityTool.h:94
Trk::vertex
@ vertex
Definition: MeasurementType.h:21
SG::WriteDecorHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
xAOD::Vertex_v1
Class describing a Vertex.
Definition: Vertex_v1.h:42
xAOD::TruthVertex_v1::z
float z() const
Vertex longitudinal distance along the beam line form the origin.
DEBUG
#define DEBUG
Definition: page_access.h:11
makeTRTBarrelCans.dx
tuple dx
Definition: makeTRTBarrelCans.py:20
python.changerun.pv
pv
Definition: changerun.py:79
xAOD::Egamma_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition: Egamma_v1.cxx:65
DerivationFramework::EGElectronAmbiguityTool::m_containerName
SG::ReadHandleKey< xAOD::ElectronContainer > m_containerName
Definition: EGElectronAmbiguityTool.h:45
xAOD::EgammaHelpers::getOriginalTrackParticle
const xAOD::TrackParticle * getOriginalTrackParticle(const xAOD::Electron *el)
Helper function for getting the "Original" Track Particle (i.e before GSF) via the electron.
Definition: ElectronxAODHelpers.cxx:11
xAOD::Egamma_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: Egamma_v1.cxx:70
DerivationFramework::EGElectronAmbiguityTool::decorateSimple
StatusCode decorateSimple(DecorHandles &dh, std::unique_ptr< ConstDataVector< xAOD::TrackParticleContainer >> &tpC, const xAOD::Electron *ele, const xAOD::Vertex *pvtx) const
Definition: EGElectronAmbiguityTool.cxx:293
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
DerivationFramework::EGElectronAmbiguityTool::m_dtzv
SG::WriteDecorHandleKey< xAOD::ElectronContainer > m_dtzv
Definition: EGElectronAmbiguityTool.h:97
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
MCP::ScaleSmearParam::r1
@ r1
MuonParameters::beta
@ beta
Definition: MuonParamDefs.h:144
xAOD::TrackParticle_v1::theta
float theta() const
Returns the parameter, which has range 0 to .
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
DerivationFramework::EGElectronAmbiguityTool::m_meeAtVtxECCut
double m_meeAtVtxECCut
Definition: EGElectronAmbiguityTool.h:132
DerivationFramework::EGElectronAmbiguityTool::m_meeICCut
double m_meeICCut
Definition: EGElectronAmbiguityTool.h:132
xAOD::TrackParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .)
ParticleConstants::PDG2011::electronMassInMeV
constexpr double electronMassInMeV
the mass of the electron (in MeV)
Definition: ParticleConstants.h:26