 |
ATLAS Offline Software
|
#include <EGElectronAmbiguityTool.h>
|
| SG::ReadHandleKey< xAOD::ElectronContainer > | m_containerName |
| |
| SG::ReadHandleKey< xAOD::VertexContainer > | m_VtxContainerName |
| |
| SG::ReadHandleKey< xAOD::TrackParticleContainer > | m_tpContainerName |
| |
| SG::ReadHandleKey< xAOD::TrackParticleContainer > | m_tpCName |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_drv { this, "DFCommonSimpleConvRadius", m_containerName, "DFCommonSimpleConvRadius", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dphiv { this, "DFCommonSimpleConvPhi", m_containerName, "DFCommonSimpleConvPhi", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dmee { this, "DFCommonSimpleMee", m_containerName, "DFCommonSimpleMee", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dmeeVtx { this, "DFCommonSimpleMeeAtVtx", m_containerName, "DFCommonSimpleMeeAtVtx", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dsep { this, "DFCommonSimpleSeparation", m_containerName, "DFCommonSimpleSeparation", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dambi { this, "DFCommonAddAmbiguity", m_containerName, "DFCommonAddAmbiguity", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dtrv { this, "DFCommonProdTrueRadius", m_containerName, "DFCommonProdTrueRadius", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dtpv { this, "DFCommonProdTruePhi", m_containerName, "DFCommonProdTruePhi", "" } |
| |
| SG::WriteDecorHandleKey< xAOD::ElectronContainer > | m_dtzv { this, "DFCommonProdTrueZ", m_containerName, "DFCommonProdTrueZ", "" } |
| |
| Gaudi::Property< bool > | m_isMC {this, "isMC", false} |
| |
| Gaudi::Property< double > | m_elepTCut |
| |
| Gaudi::Property< std::string > | m_idCut |
| |
| Gaudi::Property< unsigned int > | m_nSiCut |
| |
| Gaudi::Property< double > | m_dctCut |
| |
| Gaudi::Property< double > | m_sepCut |
| |
| Gaudi::Property< double > | m_dzCut |
| |
| Gaudi::Property< double > | m_rvECCut |
| |
| Gaudi::Property< double > | m_meeAtVtxECCut |
| |
| Gaudi::Property< double > | m_meeICCut |
| |
Definition at line 29 of file EGElectronAmbiguityTool.h.
◆ addBranches()
| StatusCode DerivationFramework::EGElectronAmbiguityTool::addBranches |
( |
const EventContext & |
ctx | ) |
const |
|
finaloverridevirtual |
Definition at line 83 of file EGElectronAmbiguityTool.cxx.
95 for (
const auto*
vertex : *vtxC) {
111 return StatusCode::FAILURE;
114 ATH_MSG_DEBUG(
"No primary vertex found. Setting default values.");
117 dh.dphiv(*iele) = -1;
119 dh.dmeeVtx(*iele) = -1;
121 dh.dambi(*iele) = -1;
126 return StatusCode::SUCCESS;
128 ATH_MSG_DEBUG(
"Pvx z = " << pvtx->z() <<
", number of electrons "
137 std::set<const xAOD::TrackParticle*> alreadyStored;
138 std::set<const xAOD::TrackParticle*> eleIDtpStored, eleGSFtpStored;
140 std::make_unique<ConstDataVector<xAOD::TrackParticleContainer>>(
143 for (
const auto* ele : *eleC) {
149 m_idCut.empty() || !aidCut.isAvailable(*ele) || !aidCut(*ele))
154 eleGSFtpStored.insert(eleGSFtp);
158 eleIDtpStored.insert(eleIDtp);
161 for (
const auto*
tp : *idtpC) {
165 closeByTracks->push_back(
tp);
166 alreadyStored.insert(
tp);
171 if (alreadyStored.find(
tp) != alreadyStored.end())
179 double dR = eleIDtp->
p4().DeltaR(
tp->p4());
180 double dz = std::abs(eleIDtp->
z0() -
tp->z0()) *
sin(eleIDtp->
theta());
181 if (dR >= 0.3 || dz >=
m_dzCut)
188 alreadyStored.insert(
tp);
190 closeByTracks->push_back(
tp);
194 if (closeByTracks->empty())
195 return StatusCode::SUCCESS;
203 << idtpC->size() <<
" , number of selected close-by tracks "
204 << closeByTracks->size() <<
" , number of GSF tracks "
206 for (
const auto* trk : eleIDtpStored)
208 << trk <<
" pt = " << trk->pt() * 1
e-3
209 <<
" eta = " << trk->eta() <<
" phi = " << trk->phi()
211 for (
const auto* trk : eleGSFtpStored)
213 << trk <<
" pt = " << trk->pt() * 1
e-3
214 <<
" eta = " << trk->eta() <<
" phi = " << trk->phi()
218 << trk <<
" pt = " << trk->pt() * 1
e-3
219 <<
" eta = " << trk->eta() <<
" phi = " << trk->phi()
223 for (
const auto* ele : *eleC) {
227 m_idCut.empty() || !aidCut.isAvailable(*ele) || !aidCut(*ele))
232 ATH_MSG_ERROR(
"Cannot decorate the electron with the simple info");
233 return StatusCode::FAILURE;
238 return StatusCode::SUCCESS;
◆ decorateSimple()
HelixToCircle Main Track Electron
HelixToCircle Other Electron Conv Track
Definition at line 243 of file EGElectronAmbiguityTool.cxx.
263 <<
" phi = " << ele->
phi() <<
" GSF trk ptr = "
264 << eleGSFtrkP <<
" ID trk ptr " << eleIDtrkP);
269 double tpvr = -1, tpvp = 9e9, tpvz = 9e9;
271 truthEl->
prodVtx() !=
nullptr) {
276 dh.dtrv(*ele) = tpvr;
277 dh.dtpv(*ele) = tpvp;
278 dh.dtzv(*ele) = tpvz;
284 double detaMin = 9e9;
287 if (
tp->charge() * eletrkP->
charge() > 0)
291 double dR = eletrkP->
p4().DeltaR(
tp->p4());
292 double dz = std::abs(eletrkP->
z0() -
tp->z0()) *
sin(eletrkP->
theta());
293 if (dR >= 0.3 || dz >=
m_dzCut)
296 double deta = std::abs(eletrkP->
eta() -
tp->eta());
297 if (deta < detaMin) {
306 double meeAtVtx = -1.;
308 bool goodConv =
false;
321 mee = (ep4 + op4).M();
322 op4.SetPhi(eletrkP->
phi());
323 meeAtVtx = (ep4 + op4).M();
326 std::vector<double> helix1, helix2;
329 helix(eletrkP, pvtx, helix1);
330 helix(otrkP, pvtx, helix2);
333 if (helix1[4] < helix2[4])
334 beta = TMath::PiOver2() - helix1[4];
336 beta = TMath::PiOver2() - helix2[4];
338 double phi1(helix1[4] +
beta);
339 if (phi1 > TMath::TwoPi())
340 phi1 -= TMath::TwoPi();
342 phi1 += TMath::TwoPi();
344 double phi2(helix2[4] +
beta);
345 if (phi2 > TMath::TwoPi())
346 phi2 -= TMath::TwoPi();
348 phi2 += TMath::TwoPi();
351 double r1 = 1 / (2. * std::abs(helix1[1]));
356 double rcenter1(helix1[3] / charge1 + r1);
357 double phicenter1(phi1 + TMath::PiOver2() * charge1);
359 double x1 = rcenter1 *
cos(phicenter1);
360 double y1 = rcenter1 *
sin(phicenter1);
363 double r2 = 1 / (2. * std::abs(helix2[1]));
368 double rcenter2(helix2[3] / charge2 + r2);
369 double phicenter2(phi2 + TMath::PiOver2() * charge2);
371 double x2 = rcenter2 *
cos(phicenter2);
372 double y2 = rcenter2 *
sin(phicenter2);
376 if (dx < 1e-9 && dx > 0.)
378 if (
dx > -1
e-9 &&
dx < 0.)
380 double slope((
y1 -
y2) /
dx);
381 double b(
y1 - slope *
x1);
385 double separation =
d -
r1 -
r2;
395 double temp1 = (cpx1 + cpx2) / 2;
396 double temp2 = slope * temp1 +
b;
400 double dct(helix1[0] - helix2[0]);
406 pv = std::atan2(convY, convX);
407 rv = sqrt(convX * convX + convY * convY);
408 if (convX *
cos(eletrkP->
phi()) + convY *
sin(eletrkP->
phi()) < 0)
417 dh.dmeeVtx(*ele) = meeAtVtx;
427 return StatusCode::SUCCESS;
◆ initialize()
| StatusCode DerivationFramework::EGElectronAmbiguityTool::initialize |
( |
| ) |
|
|
finaloverridevirtual |
◆ DecorHandles
◆ m_containerName
Initial value:{
this,
"ContainerName",
"Electrons",
"SG key of electron container"
}
Definition at line 40 of file EGElectronAmbiguityTool.h.
◆ m_dambi
◆ m_dctCut
| Gaudi::Property<double> DerivationFramework::EGElectronAmbiguityTool::m_dctCut |
|
private |
◆ m_dmee
◆ m_dmeeVtx
◆ m_dphiv
◆ m_drv
◆ m_dsep
◆ m_dtpv
◆ m_dtrv
◆ m_dtzv
◆ m_dzCut
| Gaudi::Property<double> DerivationFramework::EGElectronAmbiguityTool::m_dzCut |
|
private |
Initial value:{
this, "dzsinTCut", 0.5, "max dz sinTheta between ele and other tracks"}
Definition at line 129 of file EGElectronAmbiguityTool.h.
◆ m_elepTCut
| Gaudi::Property<double> DerivationFramework::EGElectronAmbiguityTool::m_elepTCut |
|
private |
◆ m_idCut
| Gaudi::Property<std::string> DerivationFramework::EGElectronAmbiguityTool::m_idCut |
|
private |
Initial value:{
this, "idCut", "DFCommonElectronsLHLoose", "minimal quality for an electron to be studied"}
Definition at line 119 of file EGElectronAmbiguityTool.h.
◆ m_isMC
| Gaudi::Property<bool> DerivationFramework::EGElectronAmbiguityTool::m_isMC {this, "isMC", false} |
|
private |
◆ m_meeAtVtxECCut
| Gaudi::Property<double> DerivationFramework::EGElectronAmbiguityTool::m_meeAtVtxECCut |
|
private |
Initial value:{
this, "meeAtVtxCut", 100, "maximal mass at vertex to be classified as external conversion"}
Definition at line 137 of file EGElectronAmbiguityTool.h.
◆ m_meeICCut
| Gaudi::Property<double> DerivationFramework::EGElectronAmbiguityTool::m_meeICCut |
|
private |
Initial value:{
this, "meeCut", 100, "maximal mass at primary vertex to be classified as gamma*"}
Definition at line 139 of file EGElectronAmbiguityTool.h.
◆ m_nSiCut
| Gaudi::Property<unsigned int> DerivationFramework::EGElectronAmbiguityTool::m_nSiCut |
|
private |
Initial value:{
this, "nSiCut", 7, "minimum number of Si hits in the other track"}
Definition at line 123 of file EGElectronAmbiguityTool.h.
◆ m_rvECCut
| Gaudi::Property<double> DerivationFramework::EGElectronAmbiguityTool::m_rvECCut |
|
private |
Initial value:{
this, "radiusCut", 20, "minimum radius to be classified as external conversion"}
Definition at line 135 of file EGElectronAmbiguityTool.h.
◆ m_sepCut
| Gaudi::Property<double> DerivationFramework::EGElectronAmbiguityTool::m_sepCut |
|
private |
◆ m_tpCName
Initial value:{
this,
"tpCName",
"GSFTrackParticles",
"SG key of TrackParticleInputContainer"
}
Definition at line 59 of file EGElectronAmbiguityTool.h.
◆ m_tpContainerName
Initial value:{
this,
"tpContainerName",
"InDetTrackParticles",
"SG key of track particles container"
}
Definition at line 52 of file EGElectronAmbiguityTool.h.
◆ m_VtxContainerName
Initial value:{
this,
"VtxContainerName",
"PrimaryVertices",
"SG key of vertex container"
}
Definition at line 46 of file EGElectronAmbiguityTool.h.
The documentation for this class was generated from the following files:
virtual double pt() const override final
The transverse momentum ( ) of the particle.
float phi() const
Vertex azimuthal angle.
@ VIEW_ELEMENTS
this data object is a view, it does not own its elmts
float charge() const
Returns the charge.
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
float z0() const
Returns the parameter.
const xAOD::TrackParticle * trackParticle(size_t index=0) const
Pointer to the xAOD::TrackParticle/s that match the electron candidate.
Helper class to provide constant type-safe access to aux data.
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
bool isElectron(const xAOD::Egamma *eg)
is the object an electron (not Fwd)
Class describing a truth particle in the MC record.
float perp() const
Vertex transverse distance from the beam line.
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
const xAOD::TruthParticle * getTruthParticle(const xAOD::IParticle &p)
Return the truthParticle associated to the given IParticle (if any)
std::size_t numberOfSiHits(const xAOD::TrackParticle *tp)
return the number of Si hits in the track particle
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Class describing a Vertex.
float z() const
Vertex longitudinal distance along the beam line form the origin.
virtual double pt() const override final
The transverse momentum ( ) of the particle.
const xAOD::TrackParticle * getOriginalTrackParticle(const xAOD::Electron *el)
Helper function for getting the "Original" Track Particle (i.e before GSF) via the electron.
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Class describing a TrackParticle.
float theta() const
Returns the parameter, which has range 0 to .
size_type size() const noexcept
Returns the number of elements in the collection.
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .)
constexpr double electronMassInMeV
the mass of the electron (in MeV)