ATLAS Offline Software
TrigTauPrecisionDiKaonHypoTool.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 /********************************************************************
6  *
7  * NAME: TrigTauPrecisionDiKaonHypoTool.cxx
8  * PACKAGE: Trigger/TrigHypothesis/TrigTauHypo
9  *
10  * AUTHOR: J. Silva based on TrigTauPrecisionIDHypoTool
11  * CREATED: Feb. 17, 2021
12  *
13  *********************************************************************/
14 
15 #include "AthenaMonitoringKernel/Monitored.h"
16 #include "TrigCompositeUtils/TrigCompositeUtils.h"
17 #include "xAODTau/TauJetContainer.h"
18 #include "xAODTracking/TrackParticle.h"
19 #include "GaudiKernel/SystemOfUnits.h"
20 #include "TruthUtils/ParticleConstants.h"
21 
22 #include "TrigTauPrecisionDiKaonHypoTool.h"
23 
24 
25 using namespace TrigCompositeUtils;
26 
27 TrigTauPrecisionDiKaonHypoTool::TrigTauPrecisionDiKaonHypoTool(const std::string& type,
28  const std::string& name,
29  const IInterface* parent)
30  : base_class(type, name, parent), m_decisionId(HLT::Identifier::fromToolName(name))
31 {
32 
33 }
34 
35 
36 StatusCode TrigTauPrecisionDiKaonHypoTool::initialize()
37 {
38  ATH_MSG_DEBUG(name() << ": in initialize()");
39 
40  ATH_MSG_DEBUG("TrigTauPrecisionDiKaonHypoTool will cut on:");
41  ATH_MSG_DEBUG(" - PtMin: " << m_ptMin.value());
42  ATH_MSG_DEBUG(" - NTracksMin: " << m_numTrackMin.value());
43  ATH_MSG_DEBUG(" - NTracksMax: " << m_numTrackMax.value());
44  ATH_MSG_DEBUG(" - NIsoTracksMax: " << m_numIsoTrackMax.value());
45 
46  ATH_MSG_DEBUG(" - massTrkSysMin: " << m_massTrkSysMin.value());
47  ATH_MSG_DEBUG(" - massTrkSysMax: " << m_massTrkSysMax.value());
48  ATH_MSG_DEBUG(" - massTrkSysKaonMin: " << m_massTrkSysKaonMin.value());
49  ATH_MSG_DEBUG(" - massTrkSysKaonMax: " << m_massTrkSysKaonMax.value());
50  ATH_MSG_DEBUG(" - massTrkSysKaonPiMin: " << m_massTrkSysKaonPiMin.value());
51  ATH_MSG_DEBUG(" - massTrkSysKaonPiMax: " << m_massTrkSysKaonPiMax.value());
52  ATH_MSG_DEBUG(" - targetMassTrkSysKaonPi: " << m_targetMassTrkSysKaonPi.value());
53  ATH_MSG_DEBUG(" - leadTrkPtMin: " << m_leadTrkPtMin.value());
54  ATH_MSG_DEBUG(" - dRmaxMax: " << m_dRmaxMax.value());
55  ATH_MSG_DEBUG(" - etOverPtLeadTrkMin: " << m_etOverPtLeadTrkMin.value());
56  ATH_MSG_DEBUG(" - etOverPtLeadTrkMax: " << m_etOverPtLeadTrkMax.value());
57  ATH_MSG_DEBUG(" - EMPOverTrkSysPMax: " << m_EMPOverTrkSysPMax.value());
58 
59  if(m_massTrkSysMin > m_massTrkSysMax || m_numTrackMin > m_numTrackMax
60  || m_massTrkSysKaonPiMin > m_massTrkSysKaonPiMax || m_massTrkSysKaonMin > m_massTrkSysKaonMax
61  || m_etOverPtLeadTrkMin > m_etOverPtLeadTrkMax) {
62  ATH_MSG_ERROR("Invalid tool configuration!");
63  return StatusCode::FAILURE;
64  }
65 
66  return StatusCode::SUCCESS;
67 }
68 
69 
70 bool TrigTauPrecisionDiKaonHypoTool::decide(const ITrigTauPrecisionHypoTool::ToolInfo& input) const
71 {
72  ATH_MSG_DEBUG(name() << ": in execute()");
73 
74  auto NInputTaus = Monitored::Scalar<int>("NInputTaus", -1);
75  auto passedCuts = Monitored::Scalar<int>("CutCounter", 0);
76  auto PtAccepted = Monitored::Scalar<float>("PtAccepted", -1);
77  auto NTracksAccepted = Monitored::Scalar<int>("NTracksAccepted", -1);
78  auto NIsoTracksAccepted = Monitored::Scalar<int>("NIsoTracksAccepted", -1);
79  auto massTrkSysAccepted = Monitored::Scalar<float>("massTrkSysAccepted", -10);
80  auto massTrkSysKaonAccepted = Monitored::Scalar<float>("massTrkSysKaonAccepted", -10);
81  auto massTrkSysKaonPiAccepted = Monitored::Scalar<float>("massTrkSysKaonPiAccepted", -10);
82  auto leadTrkPtAccepted = Monitored::Scalar<float>("leadTrkPtAccepted", -10);
83  auto dRAccepted = Monitored::Scalar<float>("dRAccepted", -1);
84  auto etOverPtLeadTrkAccepted = Monitored::Scalar<float>("etOverPtLeadTrkAccepted", -10);
85  auto EMOverTrkSysPAccepted = Monitored::Scalar<float>("EMOverTrkSysPAccepted", -10);
86 
87  auto monitorIt = Monitored::Group(m_monTool,
88  NInputTaus, passedCuts,
89  PtAccepted, NTracksAccepted, NIsoTracksAccepted,
90  massTrkSysAccepted, massTrkSysKaonAccepted, massTrkSysKaonPiAccepted,
91  leadTrkPtAccepted, dRAccepted, etOverPtLeadTrkAccepted, EMOverTrkSysPAccepted);
92 
93 
94  // Tau pass flag
95  bool pass = false;
96 
97  if(m_acceptAll) {
98  pass = true;
99  ATH_MSG_DEBUG("AcceptAll property is set: taking all events");
100  }
101 
102  // Debugging location of the TauJet RoI
103  ATH_MSG_DEBUG("Input RoI eta: " << input.roi->eta() << ", phi: " << input.roi->phi() << ", z: " << input.roi->zed());
104 
105  const xAOD::TauJetContainer* TauContainer = input.tauContainer;
106  NInputTaus = TauContainer->size();
107  // There should only be a single TauJet in the TauJetContainer; just in case we still run the loop
108  for(const xAOD::TauJet* Tau : *TauContainer) {
109  ATH_MSG_DEBUG(" New HLT TauJet candidate:");
110  passedCuts++;
111 
112 
113  //---------------------------------------------------------------
114  // Calibrated tau pT cut ('idperf' step)
115  //---------------------------------------------------------------
116  float pT = Tau->pt();
117  ATH_MSG_DEBUG(" pT: " << pT / Gaudi::Units::GeV);
118 
119  if(!(pT > m_ptMin)) continue;
120  passedCuts++;
121  PtAccepted = pT / Gaudi::Units::GeV;
122 
123 
124  //---------------------------------------------------------------
125  // Track counting ('perf' step)
126  //---------------------------------------------------------------
127  int numTrack = Tau->nTracks();
128  int numIsoTrack = Tau->nTracksIsolation();
129 
130  ATH_MSG_DEBUG(" N Tracks: " << numTrack);
131  ATH_MSG_DEBUG(" N Iso Tracks: " << numIsoTrack);
132 
133  // Apply NTrackMin/Max and NIsoTrackMax cuts:
134  if(!(numTrack >= m_numTrackMin && numTrack <= m_numTrackMax)) continue;
135  if(!(numIsoTrack <= m_numIsoTrackMax)) continue;
136  passedCuts++;
137  NTracksAccepted = numTrack;
138  NIsoTracksAccepted = numIsoTrack;
139 
140 
141  //---------------------------------------------------------------
142  // Leading track pT cut
143  //---------------------------------------------------------------
144  float leadTrkPt = -1;
145  Tau->detail(xAOD::TauJetParameters::leadTrkPt, leadTrkPt);
146  ATH_MSG_DEBUG(" leadTrkPt: " << leadTrkPt / Gaudi::Units::GeV);
147  if(!(leadTrkPt > m_leadTrkPtMin)) continue;
148  passedCuts++;
149  leadTrkPtAccepted = leadTrkPt / Gaudi::Units::GeV;
150 
151 
152  //---------------------------------------------------------------
153  // Cuts on the track-system mass, for different meson hypothesis
154  //---------------------------------------------------------------
155  float massTrkSys = -1.;
156  Tau->detail(xAOD::TauJetParameters::massTrkSys, massTrkSys);
157  ATH_MSG_DEBUG(" massTrkSys: " << massTrkSys / Gaudi::Units::GeV);
158 
159  // For the dikaon mass hypothesis, compute the invariant mass with the kaon mass
160  // Also cache the p4 of tracks, to use for other mass hypothesis
161  TLorentzVector my_kaons;
162  std::vector<TLorentzVector> my_trks;
163  for(unsigned int i = 0; i < Tau->nTracks(); ++i) {
164  const xAOD::TrackParticle* trk = nullptr;
165  TLorentzVector tmpKaon;;
166 
167  try {
168  trk = Tau->track(i)->track();
169  } catch(const std::exception& e) {
170  ATH_MSG_WARNING(" Failed to get tau track link!");
171  }
172 
173  if(trk) {
174  tmpKaon.SetPtEtaPhiM(trk->pt(), trk->eta(), trk->phi(), ParticleConstants::chargedKaonMassInMeV);
175  my_trks.push_back(trk->p4());
176  }
177 
178  my_kaons = my_kaons + tmpKaon;
179  }
180 
181  float massTrkSysKaon = my_kaons.M();
182  ATH_MSG_DEBUG(" massTrkSys with kaon mass hypo: " << massTrkSysKaon / Gaudi::Units::GeV);
183 
184  // kaon+pi mass hypo
185  float finalKPiMass = 0;
186  if(my_trks.size() == 2) {
187  TLorentzVector tmpKaon;
188 
189  tmpKaon.SetPtEtaPhiM(my_trks.at(0).Pt(), my_trks.at(0).Eta(), my_trks.at(0).Phi(), ParticleConstants::chargedKaonMassInMeV);
190  TLorentzVector tmpPion = my_trks.at(1);
191  float kPiMass1 = (tmpKaon+tmpPion).M();
192 
193  tmpKaon.SetPtEtaPhiM(my_trks.at(1).Pt(), my_trks.at(1).Eta(), my_trks.at(1).Phi(), ParticleConstants::chargedKaonMassInMeV);
194  tmpPion = my_trks.at(0);
195  float kPiMass2 = (tmpKaon+tmpPion).M();
196 
197  if(std::abs(kPiMass1 - m_targetMassTrkSysKaonPi) < std::abs(kPiMass2 - m_targetMassTrkSysKaonPi)) {
198  finalKPiMass = kPiMass1;
199  } else {
200  finalKPiMass = kPiMass2;
201  }
202  }
203  float massTrkSysKaonPi = finalKPiMass;
204  ATH_MSG_DEBUG(" massTrkSys with kaon+pi mass hypo: " << massTrkSysKaonPi / Gaudi::Units::GeV);
205 
206  if(!(massTrkSys > m_massTrkSysMin && massTrkSys < m_massTrkSysMax)) continue;
207  passedCuts++;
208  massTrkSysAccepted = massTrkSys / Gaudi::Units::GeV;
209 
210  if(!(massTrkSysKaon > m_massTrkSysKaonMin && massTrkSysKaon < m_massTrkSysKaonMax)) continue;
211  passedCuts++;
212  massTrkSysKaonAccepted = massTrkSysKaon / Gaudi::Units::GeV;
213 
214  // Use '>=' here, otherwise the singlepion chain would fail!
215  if(!(massTrkSysKaonPi >= m_massTrkSysKaonPiMin && massTrkSysKaonPi < m_massTrkSysKaonPiMax)) continue;
216  passedCuts++;
217  massTrkSysKaonPiAccepted = massTrkSysKaonPi / Gaudi::Units::GeV;
218 
219 
220  //---------------------------------------------------------------
221  // Cut on E_{T}^{EM} / p_T^{trk sys}
222  //---------------------------------------------------------------
223  float EMPOverTrkSysP = -1;
224  Tau->detail(xAOD::TauJetParameters::EMPOverTrkSysP, EMPOverTrkSysP);
225  ATH_MSG_DEBUG(" EMPOverTrkSysP: " << EMPOverTrkSysP);
226  if(!(EMPOverTrkSysP < m_EMPOverTrkSysPMax)) continue;
227  passedCuts++;
228  EMOverTrkSysPAccepted = EMPOverTrkSysP;
229 
230 
231  //---------------------------------------------------------------
232  // Cut on E_{T} / p_T^{lead trk}
233  //---------------------------------------------------------------
234  float etOverPtLeadTrk = -1;
235  Tau->detail(xAOD::TauJetParameters::etOverPtLeadTrk, etOverPtLeadTrk);
236  ATH_MSG_DEBUG(" etOverPtLeadTrk: " << etOverPtLeadTrk);
237  if(!(etOverPtLeadTrk > m_etOverPtLeadTrkMin && etOverPtLeadTrk < m_etOverPtLeadTrkMax)) continue;
238  passedCuts++;
239  etOverPtLeadTrkAccepted = etOverPtLeadTrk;
240 
241 
242  //---------------------------------------------------------------
243  // Cut on Max dR(track, tau)
244  //---------------------------------------------------------------
245  float dRmax = -1;
246  Tau->detail(xAOD::TauJetParameters::dRmax, dRmax);
247  ATH_MSG_DEBUG(" dRmax: " << dRmax);
248  if(!(dRmax < m_dRmaxMax)) continue;
249  passedCuts++;
250  dRAccepted = dRmax;
251 
252 
253  //---------------------------------------------------------
254  // At least one Tau passed all the cuts. Accept the event!
255  //---------------------------------------------------------
256  pass = true;
257 
258  ATH_MSG_DEBUG(" Pass hypo tool: " << pass);
259  }
260 
261  return pass;
262 }
263 
264 StatusCode TrigTauPrecisionDiKaonHypoTool::decide(std::vector<ITrigTauPrecisionHypoTool::ToolInfo>& input) const
265 {
266  for(auto& i : input) {
267  if(passed(m_decisionId.numeric(), i.previousDecisionIDs)) {
268  if(decide(i)) {
269  addDecisionID(m_decisionId, i.decision);
270  }
271  }
272  }
273 
274  return StatusCode::SUCCESS;
275 }
276 
TrigTauPrecisionDiKaonHypoTool::initialize
virtual StatusCode initialize() override
Definition: TrigTauPrecisionDiKaonHypoTool.cxx:36
TrigTauPrecisionDiKaonHypoTool::m_massTrkSysMax
Gaudi::Property< float > m_massTrkSysMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:49
TrigTauPrecisionDiKaonHypoTool::m_numTrackMin
Gaudi::Property< int > m_numTrackMin
Definition: TrigTauPrecisionDiKaonHypoTool.h:44
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
TrigTauPrecisionDiKaonHypoTool::m_numIsoTrackMax
Gaudi::Property< int > m_numIsoTrackMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:45
CalculateHighPtTerm.pT
pT
Definition: ICHEP2016/CalculateHighPtTerm.py:57
GeV
#define GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:18
TrigTauPrecisionDiKaonHypoTool::m_EMPOverTrkSysPMax
Gaudi::Property< float > m_EMPOverTrkSysPMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:62
TrigDefs::Group
Group
Properties of a chain group.
Definition: GroupProperties.h:13
TrigCompositeUtils::passed
bool passed(DecisionID id, const DecisionIDContainer &idSet)
checks if required decision ID is in the set of IDs in the container
Definition: TrigCompositeUtilsRoot.cxx:117
HLT::Identifier::numeric
TrigCompositeUtils::DecisionID numeric() const
numeric ID
Definition: TrigCompositeUtils/TrigCompositeUtils/HLTIdentifier.h:41
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: TrackParticle_v1.cxx:78
TrigTauPrecisionDiKaonHypoTool::decide
virtual StatusCode decide(std::vector< ITrigTauPrecisionHypoTool::ToolInfo > &input) const override
Definition: TrigTauPrecisionDiKaonHypoTool.cxx:264
TrigCompositeUtils::addDecisionID
void addDecisionID(DecisionID id, Decision *d)
Appends the decision (given as ID) to the decision object.
Definition: TrigCompositeUtilsRoot.cxx:58
TrigTauPrecisionDiKaonHypoTool::m_decisionId
HLT::Identifier m_decisionId
Definition: TrigTauPrecisionDiKaonHypoTool.h:40
TrigTauPrecisionDiKaonHypoTool::m_ptMin
Gaudi::Property< float > m_ptMin
Definition: TrigTauPrecisionDiKaonHypoTool.h:42
ParticleConstants::PDG2011::chargedKaonMassInMeV
constexpr double chargedKaonMassInMeV
the mass of the charged kaon (in MeV)
Definition: ParticleConstants.h:50
Tau
Definition: EfficiencyPtPlots.cxx:10
ITrigTauPrecisionHypoTool::ToolInfo
Definition: ITrigTauPrecisionHypoTool.h:25
TrigTauPrecisionDiKaonHypoTool.h
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
TrigTauPrecisionDiKaonHypoTool::m_etOverPtLeadTrkMin
Gaudi::Property< float > m_etOverPtLeadTrkMin
Definition: TrigTauPrecisionDiKaonHypoTool.h:61
xAOD::TauJetParameters::dRmax
@ dRmax
Get maximal dR of tracks associated to calo-seeded tau.
Definition: TauDefs.h:226
xAOD::TrackParticle_v1::p4
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
Definition: TrackParticle_v1.cxx:130
TrigCompositeUtils.h
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
HLT
It used to be useful piece of code for replacing actual SG with other store of similar functionality ...
Definition: HLTResultReader.h:26
ParticleConstants.h
lumiFormat.i
int i
Definition: lumiFormat.py:85
TrigTauPrecisionDiKaonHypoTool::m_massTrkSysKaonMax
Gaudi::Property< float > m_massTrkSysKaonMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:52
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
xAOD::TauJet_v3
Class describing a tau jet.
Definition: TauJet_v3.h:41
TrigTauPrecisionDiKaonHypoTool::TrigTauPrecisionDiKaonHypoTool
TrigTauPrecisionDiKaonHypoTool(const std::string &type, const std::string &name, const IInterface *parent)
Definition: TrigTauPrecisionDiKaonHypoTool.cxx:27
PlotPulseshapeFromCool.input
input
Definition: PlotPulseshapeFromCool.py:106
calibdata.exception
exception
Definition: calibdata.py:495
test_pyathena.parent
parent
Definition: test_pyathena.py:15
xAOD::TauJetParameters::etOverPtLeadTrk
@ etOverPtLeadTrk
Definition: TauDefs.h:158
TrigTauPrecisionDiKaonHypoTool::m_massTrkSysKaonPiMax
Gaudi::Property< float > m_massTrkSysKaonPiMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:55
TauJetContainer.h
TrigTauPrecisionDiKaonHypoTool::m_monTool
ToolHandle< GenericMonitoringTool > m_monTool
Definition: TrigTauPrecisionDiKaonHypoTool.h:66
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
Monitored.h
Header file to be included by clients of the Monitored infrastructure.
xAOD::TauJetParameters::massTrkSys
@ massTrkSys
Definition: TauDefs.h:161
TrigTauPrecisionDiKaonHypoTool::m_etOverPtLeadTrkMax
Gaudi::Property< float > m_etOverPtLeadTrkMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:60
TrigTauPrecisionDiKaonHypoTool::m_targetMassTrkSysKaonPi
Gaudi::Property< float > m_targetMassTrkSysKaonPi
Definition: TrigTauPrecisionDiKaonHypoTool.h:56
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
TrigTauPrecisionDiKaonHypoTool::m_massTrkSysKaonMin
Gaudi::Property< float > m_massTrkSysKaonMin
Definition: TrigTauPrecisionDiKaonHypoTool.h:51
TrigTauPrecisionDiKaonHypoTool::m_dRmaxMax
Gaudi::Property< float > m_dRmaxMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:59
TrigTauPrecisionDiKaonHypoTool::m_massTrkSysMin
Gaudi::Property< float > m_massTrkSysMin
Definition: TrigTauPrecisionDiKaonHypoTool.h:48
TrackParticle.h
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
TrigTauPrecisionDiKaonHypoTool::m_numTrackMax
Gaudi::Property< int > m_numTrackMax
Definition: TrigTauPrecisionDiKaonHypoTool.h:43
TrigCompositeUtils
Definition: Event/xAOD/xAODTrigger/xAODTrigger/TrigComposite.h:19
TrigTauPrecisionDiKaonHypoTool::m_leadTrkPtMin
Gaudi::Property< float > m_leadTrkPtMin
Definition: TrigTauPrecisionDiKaonHypoTool.h:58
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
Monitored::Scalar
Declare a monitored scalar variable.
Definition: MonitoredScalar.h:34
TrigTauPrecisionDiKaonHypoTool::m_acceptAll
Gaudi::Property< bool > m_acceptAll
Definition: TrigTauPrecisionDiKaonHypoTool.h:64
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
TrigTauPrecisionDiKaonHypoTool::m_massTrkSysKaonPiMin
Gaudi::Property< float > m_massTrkSysKaonPiMin
Definition: TrigTauPrecisionDiKaonHypoTool.h:54
xAOD::TauJetParameters::EMPOverTrkSysP
@ EMPOverTrkSysP
Definition: TauDefs.h:279
xAOD::TrackParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .)
xAOD::TauJetParameters::leadTrkPt
@ leadTrkPt
Definition: TauDefs.h:159
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Tue Sep 2 2025 21:24:48 for ATLAS Offline Software by   doxygen 1.8.18