ATLAS Offline Software
xAODTauJetAuxContainerCnv_v2.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 // System include(s):
6 #include <stdexcept>
7 
8 // Gaudi/Athena include(s):
9 #include "GaudiKernel/MsgStream.h"
10 #include "GaudiKernel/ThreadLocalContext.h"
11 
12 // EDM include(s):
13 #include "xAODTau/versions/TauJetContainer_v2.h"
14 #include "xAODTau/TauJetContainer.h"
15 #include "xAODTau/TauTrackContainer.h"
16 #include "xAODTau/TauTrackAuxContainer.h"
17 #include "xAODTracking/TrackParticleContainer.h"
18 
19 // Local include(s):
20 #include "xAODTauJetAuxContainerCnv_v2.h"
21 
22 // StoreGateSvc
23 #include "StoreGate/StoreGateSvc.h"
24 #include "GaudiKernel/ServiceHandle.h"
25 
26 
27 xAODTauJetAuxContainerCnv_v2::xAODTauJetAuxContainerCnv_v2()
28 {
29 }
30 
31 
32 void xAODTauJetAuxContainerCnv_v2::
33 persToTransWithKey( const xAOD::TauJetAuxContainer_v2* oldObj,
34  xAOD::TauJetAuxContainer* newObj,
35  const std::string& key,
36  MsgStream& log ) const {
37 
38  ServiceHandle<StoreGateSvc> evtStore ("StoreGateSvc", "tauJetCnv_v2");
39  if (evtStore.retrieve().isFailure()) {
40  log << MSG::ERROR << "Cannot get StoreGateHandle" << endmsg;
41  return;
42  }
43 
44  std::string tauTrackContName=key;
45  if (tauTrackContName.ends_with ("Aux.")) {
46  tauTrackContName.resize (tauTrackContName.size()-4);
47  }
48 
49  //example names:
50  //TauJets : Jets --> Tracks
51  //HLT_xAOD__TauJetContainer_TrigTauRecMerged Jet --> Track; +=Tracks
52  //HLT_xAOD__TauJetContainer_TrigTauRecPreselection ""
53  if(tauTrackContName.find("Jet") != std::string::npos){
54  tauTrackContName.replace( tauTrackContName.find("Jet"), 3, "Track" );
55  if(tauTrackContName.find("HLT") != std::string::npos) tauTrackContName+="Tracks";
56  }
57  else if (key.length() > 0) {
58  log << MSG::ERROR << "Cannot decipher name TauTrackContainer should have" << endmsg;
59  return;
60  }
61 
62  std::string tauTrackAuxContName=tauTrackContName+"Aux.";
63 
64  xAOD::TauTrackContainer* pTracks = nullptr;
65  xAOD::TauTrackAuxContainer* pAuxTracks = nullptr;
66  if(key.length()){
67  //if reading data, then trigger calls T/P converter directly
68  //and the key is not set. In this case, forget about TauTracks
69  pTracks = new xAOD::TauTrackContainer();
70  pAuxTracks = new xAOD::TauTrackAuxContainer();
71  pTracks->setStore(pAuxTracks);
72 
73  if(evtStore->record(pTracks, tauTrackContName).isFailure() ||
74  evtStore->record(pAuxTracks, tauTrackAuxContName).isFailure()){
75  log << MSG::DEBUG << "Couldn't Record TauTracks" << endmsg;
76  return;
77  }
78  }
79 
80  // Clear the transient object:
81  newObj->resize( 0 );
82 
83  // Set up interface containers on top of them:
84  xAOD::TauJetContainer_v2 oldInt;
85  for( size_t i = 0; i < oldObj->size(); ++i ) {
86  oldInt.push_back( new xAOD::TauJet_v2() );
87  }
88  oldInt.setStore( oldObj );
89  xAOD::TauJetContainer newInt;
90  newInt.setStore( newObj );
91 
92  // Loop over the interface objects, and do the conversion with their help:
93  for( const xAOD::TauJet_v2* oldTau : oldInt ) {
94 
95  // Add an object to the output container:
96  xAOD::TauJet* newTau = new xAOD::TauJet;
97  newInt.push_back( newTau );
98 
99  //
100  // Copy the 4-momentum variables:
101  //
102  newTau->setP4( oldTau->pt(), oldTau->eta(), oldTau->phi(), oldTau->m() );
103  newTau->setP4( xAOD::TauJetParameters::JetSeed, oldTau->ptJetSeed(), oldTau->etaJetSeed(), oldTau->phiJetSeed(), oldTau->mJetSeed() );
104  newTau->setP4( xAOD::TauJetParameters::DetectorAxis, oldTau->ptDetectorAxis(), oldTau->etaDetectorAxis(), oldTau->phiDetectorAxis(), oldTau->mDetectorAxis() );
105  newTau->setP4( xAOD::TauJetParameters::IntermediateAxis, oldTau->ptIntermediateAxis(), oldTau->etaIntermediateAxis(), oldTau->phiIntermediateAxis(), oldTau->mIntermediateAxis() );
106  newTau->setP4( xAOD::TauJetParameters::TauEnergyScale, oldTau->ptTauEnergyScale(), oldTau->etaTauEnergyScale(), oldTau->phiTauEnergyScale(), oldTau->mTauEnergyScale() );
107  newTau->setP4( xAOD::TauJetParameters::TauEtaCalib, oldTau->ptTauEtaCalib(), oldTau->etaTauEtaCalib(), oldTau->phiTauEtaCalib(), oldTau->mTauEtaCalib() );
108  newTau->setP4( xAOD::TauJetParameters::PanTauCellBased, oldTau->ptPanTauCellBased(), oldTau->etaPanTauCellBased(), oldTau->phiPanTauCellBased(), oldTau->mPanTauCellBased() );
109  newTau->setP4( xAOD::TauJetParameters::PanTauCellBasedProto, oldTau->ptPanTauCellBasedProto(), oldTau->etaPanTauCellBasedProto(), oldTau->phiPanTauCellBasedProto(), oldTau->mPanTauCellBasedProto() );
110 
111  // ROI and charge
112  newTau->setROIWord( oldTau->ROIWord() );
113  newTau->setCharge( oldTau->charge() );
114  newTau->setTrackFilterProngs( oldTau->trackFilterProngs() );
115  newTau->setTrackFilterQuality( oldTau->trackFilterQuality() );
116  newTau->setPi0ConeDR( oldTau->pi0ConeDR() );
117 
118  //
119  //copy PID variables
120  //
121  newTau->setDiscriminant(xAOD::TauJetParameters::BDTJetScoreSigTrans , oldTau->discriminant(xAOD::TauJetParameters::BDTJetScoreSigTrans) );
122  newTau->setDiscriminant(xAOD::TauJetParameters::BDTJetScore , oldTau->discriminant(xAOD::TauJetParameters::BDTJetScore) );
123  newTau->setDiscriminant(xAOD::TauJetParameters::BDTEleScore , oldTau->discriminant(xAOD::TauJetParameters::BDTEleScore) );
124  newTau->setIsTau(xAOD::TauJetParameters::JetBDTSigLoose , oldTau->isTau(xAOD::TauJetParameters::JetBDTSigLoose) );
125  newTau->setIsTau(xAOD::TauJetParameters::JetBDTSigMedium , oldTau->isTau(xAOD::TauJetParameters::JetBDTSigMedium) );
126  newTau->setIsTau(xAOD::TauJetParameters::JetBDTSigTight , oldTau->isTau(xAOD::TauJetParameters::JetBDTSigTight) );
127  newTau->setIsTau(xAOD::TauJetParameters::EleBDTLoose , oldTau->isTau(xAOD::TauJetParameters::EleBDTLoose) );
128  newTau->setIsTau(xAOD::TauJetParameters::EleBDTMedium , oldTau->isTau(xAOD::TauJetParameters::EleBDTMedium) );
129  newTau->setIsTau(xAOD::TauJetParameters::EleBDTTight , oldTau->isTau(xAOD::TauJetParameters::EleBDTTight) );
130 
131  //
132  //set individual int type details variables
133  //
134  newTau->setDetail(xAOD::TauJetParameters::nCharged , oldTau->detail<int>(xAOD::TauJetParameters::nCharged) );
135  newTau->setDetail(xAOD::TauJetParameters::numCells , oldTau->detail<int>(xAOD::TauJetParameters::numCells) );
136  newTau->setDetail(xAOD::TauJetParameters::numTopoClusters , oldTau->detail<int>(xAOD::TauJetParameters::numTopoClusters) );
137  newTau->setDetail(xAOD::TauJetParameters::nStrip , oldTau->detail<int>(xAOD::TauJetParameters::nStrip) );
138 
139  //
140  //set individual float type details variables
141  //
142  newTau->setDetail(xAOD::TauJetParameters::ipZ0SinThetaSigLeadTrk , oldTau->detail<float>(xAOD::TauJetParameters::ipZ0SinThetaSigLeadTrk) );
143  newTau->setDetail(xAOD::TauJetParameters::etOverPtLeadTrk , oldTau->detail<float>(xAOD::TauJetParameters::etOverPtLeadTrk) );
144  newTau->setDetail(xAOD::TauJetParameters::leadTrkPt , oldTau->detail<float>(xAOD::TauJetParameters::leadTrkPt) );
145  newTau->setDetail(xAOD::TauJetParameters::ipSigLeadTrk , oldTau->detail<float>(xAOD::TauJetParameters::ipSigLeadTrk) );
146  newTau->setDetail(xAOD::TauJetParameters::massTrkSys , oldTau->detail<float>(xAOD::TauJetParameters::massTrkSys) );
147  newTau->setDetail(xAOD::TauJetParameters::trkWidth2 , oldTau->detail<float>(xAOD::TauJetParameters::trkWidth2) );
148  newTau->setDetail(xAOD::TauJetParameters::trFlightPathSig , oldTau->detail<float>(xAOD::TauJetParameters::trFlightPathSig) );
149  newTau->setDetail(xAOD::TauJetParameters::numEffTopoClusters , oldTau->detail<float>(xAOD::TauJetParameters::numEffTopoClusters) );
150  newTau->setDetail(xAOD::TauJetParameters::topoInvMass , oldTau->detail<float>(xAOD::TauJetParameters::topoInvMass) );
151  newTau->setDetail(xAOD::TauJetParameters::effTopoInvMass , oldTau->detail<float>(xAOD::TauJetParameters::effTopoInvMass) );
152  newTau->setDetail(xAOD::TauJetParameters::topoMeanDeltaR , oldTau->detail<float>(xAOD::TauJetParameters::topoMeanDeltaR) );
153  newTau->setDetail(xAOD::TauJetParameters::effTopoMeanDeltaR , oldTau->detail<float>(xAOD::TauJetParameters::effTopoMeanDeltaR) );
154  newTau->setDetail(xAOD::TauJetParameters::EMRadius , oldTau->detail<float>(xAOD::TauJetParameters::EMRadius) );
155  newTau->setDetail(xAOD::TauJetParameters::hadRadius , oldTau->detail<float>(xAOD::TauJetParameters::hadRadius) );
156  newTau->setDetail(xAOD::TauJetParameters::etEMAtEMScale , oldTau->detail<float>(xAOD::TauJetParameters::etEMAtEMScale) );
157  newTau->setDetail(xAOD::TauJetParameters::etHadAtEMScale , oldTau->detail<float>(xAOD::TauJetParameters::etHadAtEMScale) );
158  newTau->setDetail(xAOD::TauJetParameters::isolFrac , oldTau->detail<float>(xAOD::TauJetParameters::isolFrac) );
159  newTau->setDetail(xAOD::TauJetParameters::centFrac , oldTau->detail<float>(xAOD::TauJetParameters::centFrac) );
160  newTau->setDetail(xAOD::TauJetParameters::stripWidth2 , oldTau->detail<float>(xAOD::TauJetParameters::stripWidth2) );
161  newTau->setDetail(xAOD::TauJetParameters::trkAvgDist , oldTau->detail<float>(xAOD::TauJetParameters::trkAvgDist) );
162  newTau->setDetail(xAOD::TauJetParameters::trkRmsDist , oldTau->detail<float>(xAOD::TauJetParameters::trkRmsDist) );
163  newTau->setDetail(xAOD::TauJetParameters::lead2ClusterEOverAllClusterE , oldTau->detail<float>(xAOD::TauJetParameters::lead2ClusterEOverAllClusterE) );
164  newTau->setDetail(xAOD::TauJetParameters::lead3ClusterEOverAllClusterE , oldTau->detail<float>(xAOD::TauJetParameters::lead3ClusterEOverAllClusterE) );
165  newTau->setDetail(xAOD::TauJetParameters::caloIso , oldTau->detail<float>(xAOD::TauJetParameters::caloIso) );
166  newTau->setDetail(xAOD::TauJetParameters::caloIsoCorrected , oldTau->detail<float>(xAOD::TauJetParameters::caloIsoCorrected) );
167  newTau->setDetail(xAOD::TauJetParameters::dRmax , oldTau->detail<float>(xAOD::TauJetParameters::dRmax) );
168  newTau->setDetail(xAOD::TauJetParameters::secMaxStripEt , oldTau->detail<float>(xAOD::TauJetParameters::secMaxStripEt) );
169  newTau->setDetail(xAOD::TauJetParameters::sumEMCellEtOverLeadTrkPt , oldTau->detail<float>(xAOD::TauJetParameters::sumEMCellEtOverLeadTrkPt) );
170  newTau->setDetail(xAOD::TauJetParameters::hadLeakEt , oldTau->detail<float>(xAOD::TauJetParameters::hadLeakEt) );
171  newTau->setDetail(xAOD::TauJetParameters::cellBasedEnergyRing1 , oldTau->detail<float>(xAOD::TauJetParameters::cellBasedEnergyRing1) );
172  newTau->setDetail(xAOD::TauJetParameters::cellBasedEnergyRing2 , oldTau->detail<float>(xAOD::TauJetParameters::cellBasedEnergyRing2) );
173  newTau->setDetail(xAOD::TauJetParameters::cellBasedEnergyRing3 , oldTau->detail<float>(xAOD::TauJetParameters::cellBasedEnergyRing3) );
174  newTau->setDetail(xAOD::TauJetParameters::cellBasedEnergyRing4 , oldTau->detail<float>(xAOD::TauJetParameters::cellBasedEnergyRing4) );
175  newTau->setDetail(xAOD::TauJetParameters::cellBasedEnergyRing5 , oldTau->detail<float>(xAOD::TauJetParameters::cellBasedEnergyRing5) );
176  newTau->setDetail(xAOD::TauJetParameters::cellBasedEnergyRing6 , oldTau->detail<float>(xAOD::TauJetParameters::cellBasedEnergyRing6) );
177  newTau->setDetail(xAOD::TauJetParameters::cellBasedEnergyRing7 , oldTau->detail<float>(xAOD::TauJetParameters::cellBasedEnergyRing7) );
178  newTau->setDetail(xAOD::TauJetParameters::TRT_NHT_OVER_NLT , oldTau->detail<float>(xAOD::TauJetParameters::TRT_NHT_OVER_NLT) );
179  newTau->setDetail(xAOD::TauJetParameters::TauJetVtxFraction , oldTau->detail<float>(xAOD::TauJetParameters::TauJetVtxFraction) );
180  newTau->setDetail(xAOD::TauJetParameters::TauJetVtxFraction , oldTau->detail<float>(xAOD::TauJetParameters::TauJetVtxFraction) );
181  newTau->setDetail(xAOD::TauJetParameters::PSSFraction , oldTau->detail<float>(xAOD::TauJetParameters::PSSFraction) );
182  newTau->setDetail(xAOD::TauJetParameters::ChPiEMEOverCaloEME , oldTau->detail<float>(xAOD::TauJetParameters::ChPiEMEOverCaloEME) );
183  newTau->setDetail(xAOD::TauJetParameters::EMPOverTrkSysP , oldTau->detail<float>(xAOD::TauJetParameters::EMPOverTrkSysP) );
184  newTau->setDetail(xAOD::TauJetParameters::TESOffset , oldTau->detail<float>(xAOD::TauJetParameters::TESOffset) );
185  newTau->setDetail(xAOD::TauJetParameters::TESCalibConstant , oldTau->detail<float>(xAOD::TauJetParameters::TESCalibConstant) );
186  newTau->setDetail(xAOD::TauJetParameters::centFracCorrected , oldTau->detail<float>(xAOD::TauJetParameters::centFracCorrected) );
187  newTau->setDetail(xAOD::TauJetParameters::etOverPtLeadTrkCorrected , oldTau->detail<float>(xAOD::TauJetParameters::etOverPtLeadTrkCorrected) );
188  newTau->setDetail(xAOD::TauJetParameters::innerTrkAvgDist , oldTau->detail<float>(xAOD::TauJetParameters::innerTrkAvgDist) );
189  newTau->setDetail(xAOD::TauJetParameters::innerTrkAvgDistCorrected , oldTau->detail<float>(xAOD::TauJetParameters::innerTrkAvgDistCorrected) );
190  newTau->setDetail(xAOD::TauJetParameters::SumPtTrkFrac , oldTau->detail<float>(xAOD::TauJetParameters::SumPtTrkFrac) );
191  newTau->setDetail(xAOD::TauJetParameters::SumPtTrkFracCorrected , oldTau->detail<float>(xAOD::TauJetParameters::SumPtTrkFracCorrected) );
192 
193 
194  //
195  //set pantau details
196  //
197  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_isPanTauCandidate, oldTau->panTauDetail<int>(xAOD::TauJetParameters::PanTau_isPanTauCandidate ));
198  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_DecayModeProto, oldTau->panTauDetail<int>(xAOD::TauJetParameters::PanTau_DecayModeProto ));
199  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_DecayMode, oldTau->panTauDetail<int>(xAOD::TauJetParameters::PanTau_DecayMode ));
200  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTValue_1p0n_vs_1p1n, oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTValue_1p0n_vs_1p1n ));
201  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTValue_1p1n_vs_1pXn, oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTValue_1p1n_vs_1pXn ));
202  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTValue_3p0n_vs_3pXn, oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTValue_3p0n_vs_3pXn ));
203  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts, oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts ));
204  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1, oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1 ));
205  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2, oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2 ));
206  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts,oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts ));
207  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts, oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts ));
208  newTau->setPanTauDetail(xAOD::TauJetParameters::PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged,oldTau->panTauDetail<float>(xAOD::TauJetParameters::PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged ));
209 
210  //copy element links
211  newTau->setJetLink( oldTau->jetLink() );
212  newTau->setVertexLink( oldTau->vertexLink() );
213  newTau->setSecondaryVertexLink( oldTau->secondaryVertexLink() );
214 
215  newTau->setNeutralPFOLinks( oldTau->neutralPFOLinks() );
216  newTau->setChargedPFOLinks( oldTau->chargedPFOLinks() );
217  newTau->setPi0PFOLinks( oldTau->pi0PFOLinks() );
218  newTau->setShotPFOLinks( oldTau->shotPFOLinks() );
219 
220  //v2 doesn't have pfo element link with specific type name, so copy cellbased ones into proto
221  newTau->setProtoNeutralPFOLinks( oldTau->protoNeutralPFOLinks() );
222  newTau->setProtoChargedPFOLinks( oldTau->protoChargedPFOLinks() );
223  newTau->setProtoPi0PFOLinks( oldTau->protoPi0PFOLinks() );
224 
225  if(key.length()==0) continue;
226 
227  // Get context and look up hashed track container key.
228  const EventContext& ctx = Gaudi::Hive::currentContext();
229  ElementLink<xAOD::TauTrackContainer> dum (tauTrackContName, 0, ctx);
230  SG::sgkey_t track_sgkey = dum.key();
231  IProxyDict* sg = dum.source();
232 
233  for(unsigned int i = 0; i < oldTau->nTracks(); ++i){
234  ElementLink< xAOD::TrackParticleContainer > linkToTrackParticle = oldTau->trackLinks()[i];
235  linkToTrackParticle.toTransient (sg);
236  if(!linkToTrackParticle.isValid()) continue;
237  newTau->addTauTrackLink(ElementLink<xAOD::TauTrackContainer>(track_sgkey, pTracks->size(), ctx));
238  pTracks->push_back(std::make_unique<xAOD::TauTrack>());
239  xAOD::TauTrack& track = *pTracks->back();
240  const xAOD::TrackParticle* trackParticle=*linkToTrackParticle;
241  track.addTrackLink(linkToTrackParticle);
242  track.setP4(trackParticle->pt(), trackParticle->eta(), trackParticle->phi(), trackParticle->m());
243  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::coreTrack, true);
244  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::passTrkSelector, true);
245  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::classifiedCharged, true);
246  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::unclassified, true);
247  //ElementLink<xAOD::TauTrackContainer> linkToTauTrack;
248  //linkToTauTrack.resetWithKeyAndIndex(track_sgkey, pTracks->size()-1, ctx);
249  //newTau->addTauTrackLink(linkToTauTrack);
250  }
251 
252  for(unsigned int i = 0; i < oldTau->nWideTracks(); ++i){
253  ElementLink< xAOD::TrackParticleContainer > linkToTrackParticle = oldTau->wideTrackLinks()[i];
254  linkToTrackParticle.toTransient (sg);
255  if(!linkToTrackParticle.isValid()) continue;
256  newTau->addTauTrackLink(ElementLink<xAOD::TauTrackContainer>(track_sgkey, pTracks->size(), ctx));
257  pTracks->push_back(std::make_unique<xAOD::TauTrack>());
258  xAOD::TauTrack& track = *pTracks->back();
259  const xAOD::TrackParticle* trackParticle=*linkToTrackParticle;
260  track.addTrackLink(linkToTrackParticle);
261  track.setP4(trackParticle->pt(), trackParticle->eta(), trackParticle->phi(), trackParticle->m());
262  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::wideTrack, true);
263  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::passTrkSelector, true);
264  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::classifiedIsolation, true);
265  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::modifiedIsolationTrack, true);
266  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::unclassified, true);
267  }
268 
269  for(unsigned int i = 0; i < oldTau->nOtherTracks(); ++i){
270  ElementLink< xAOD::TrackParticleContainer > linkToTrackParticle = oldTau->otherTrackLinks()[i];
271  linkToTrackParticle.toTransient (sg);
272  if(!linkToTrackParticle.isValid()) continue;
273  newTau->addTauTrackLink(ElementLink<xAOD::TauTrackContainer>(track_sgkey, pTracks->size(), ctx));
274  pTracks->push_back(std::make_unique<xAOD::TauTrack>());
275  xAOD::TauTrack& track = *pTracks->back();
276  const xAOD::TrackParticle* trackParticle=*linkToTrackParticle;
277  track.addTrackLink(linkToTrackParticle);
278  track.setP4(trackParticle->pt(), trackParticle->eta(), trackParticle->phi(), trackParticle->m());
279  float dR=oldTau->p4(xAOD::TauJetParameters::IntermediateAxis).DeltaR(trackParticle->p4());
280  if(dR<=0.2) track.setFlag(xAOD::TauJetParameters::TauTrackFlag::coreTrack, true);
281  else track.setFlag(xAOD::TauJetParameters::TauTrackFlag::wideTrack, true);
282  track.setFlag(xAOD::TauJetParameters::TauTrackFlag::unclassified, true);
283  }
284 
285 
286  newTau->setDetail(xAOD::TauJetParameters::nChargedTracks, (int) newTau->nTracks());
287  newTau->setDetail(xAOD::TauJetParameters::nIsolatedTracks, (int) newTau->nTracks(xAOD::TauJetParameters::classifiedIsolation));
288 
289  }
290 
291  return;
292 }
293 
297 void xAODTauJetAuxContainerCnv_v2::
298 transToPersWithKey( const xAOD::TauJetAuxContainer*,
299  xAOD::TauJetAuxContainer_v2*,
300  const std::string& /*key*/,
301  MsgStream& log ) const {
302 
303  log << MSG::ERROR
304  << "Somebody called xAODTauJetAuxContainerCnv_v2::transToPers"
305  << endmsg;
306  throw std::runtime_error( "Somebody called xAODTauJetAuxContainerCnv_v2::"
307  "transToPers" );
308 
309  return;
310 }
311 
xAOD::TauJetParameters::cellBasedEnergyRing4
@ cellBasedEnergyRing4
Ring 4: 0.10 < R < 0.125.
Definition: TauDefs.h:252
xAOD::TauJetParameters::trkRmsDist
@ trkRmsDist
Get the RMS of track distance to calorimeter seed.
Definition: TauDefs.h:216
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition: TrackParticle_v1.cxx:74
xAOD::TauJetParameters::PanTau_DecayMode
@ PanTau_DecayMode
Definition: TauDefs.h:360
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts
@ PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts
Definition: TauDefs.h:372
xAOD::TauJetParameters::PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged
@ PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged
Definition: TauDefs.h:375
xAOD::TauJetParameters::nIsolatedTracks
@ nIsolatedTracks
Definition: TauDefs.h:323
xAOD::TrackParticle_v1::m
virtual double m() const override final
The invariant mass of the particle..
Definition: TrackParticle_v1.cxx:84
xAOD::TauJetParameters::trkWidth2
@ trkWidth2
Definition: TauDefs.h:162
xAOD::TauJetParameters::topoInvMass
@ topoInvMass
get invariant mass from topoclusters of jet associated to tau candidate
Definition: TauDefs.h:177
xAOD::TauJetAuxContainer_v3
Temporary container used until we have I/O for AuxStoreInternal.
Definition: TauJetAuxContainer_v3.h:39
xAOD::TauJetParameters::cellBasedEnergyRing7
@ cellBasedEnergyRing7
Ring 7: 0.2 < R < 0.4.
Definition: TauDefs.h:258
xAOD::TauJetParameters::hadLeakEt
@ hadLeakEt
Definition: TauDefs.h:232
xAOD::TauJetParameters::PanTauCellBasedProto
@ PanTauCellBasedProto
Definition: TauDefs.h:343
xAOD::TauJet_v3::setProtoNeutralPFOLinks
void setProtoNeutralPFOLinks(const PFOLinks_t &protoNeutralPFOs)
xAOD::TauJetParameters::IntermediateAxis
@ IntermediateAxis
Definition: TauDefs.h:338
xAOD::TauJetParameters::trFlightPathSig
@ trFlightPathSig
Definition: TauDefs.h:163
xAOD::TauTrackContainer
TauTrackContainer_v1 TauTrackContainer
Definition of the current TauTrack container version.
Definition: TauTrackContainer.h:16
xAOD::TauJetParameters::lead2ClusterEOverAllClusterE
@ lead2ClusterEOverAllClusterE
Get sum of 2 leading cluster energy over all cluster energy.
Definition: TauDefs.h:218
xAOD::TauJetParameters::nCharged
@ nCharged
Definition: TauDefs.h:266
xAOD::TauJet_v3::setJetLink
void setJetLink(const JetLink_t &jetLink)
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: TrackParticle_v1.cxx:78
xAOD::TauJet_v3::setPanTauDetail
void setPanTauDetail(TauJetParameters::PanTauDetails panTauDetail, int value)
Definition: TauJet_v3.cxx:367
xAOD::TauJetParameters::PanTau_BDTValue_3p0n_vs_3pXn
@ PanTau_BDTValue_3p0n_vs_3pXn
Definition: TauDefs.h:364
xAOD::TauJet_v3::setCharge
void setCharge(float)
TauTrackAuxContainer.h
xAOD::TauJet_v3::nTracks
size_t nTracks(TauJetParameters::TauTrackFlag flag=TauJetParameters::TauTrackFlag::classifiedCharged) const
Definition: TauJet_v3.cxx:488
xAOD::TauJet_v3::setProtoPi0PFOLinks
void setProtoPi0PFOLinks(const PFOLinks_t &protoPi0PFOs)
xAOD::TauJetParameters::classifiedCharged
@ classifiedCharged
Definition: TauDefs.h:406
xAOD::TauJet_v3::setNeutralPFOLinks
void setNeutralPFOLinks(const PFOLinks_t &neutralPFOs)
xAOD::TauJetParameters::TauEtaCalib
@ TauEtaCalib
Definition: TauDefs.h:340
xAOD::TauJetParameters::secMaxStripEt
@ secMaxStripEt
migrate only seedTrk_ variables which are used in reco and ID and without prefix
Definition: TauDefs.h:230
xAOD::TauJetParameters::lead3ClusterEOverAllClusterE
@ lead3ClusterEOverAllClusterE
Get sum of 3 leading cluster energy over all cluster energy.
Definition: TauDefs.h:220
xAOD::TauJetParameters::JetBDTSigMedium
@ JetBDTSigMedium
Definition: TauDefs.h:137
xAODTauJetAuxContainerCnv_v2::xAODTauJetAuxContainerCnv_v2
xAODTauJetAuxContainerCnv_v2()
Default constructor.
Definition: xAODTauJetAuxContainerCnv_v2.cxx:27
TauTrackContainer.h
xAOD::TauJetParameters::EleBDTLoose
@ EleBDTLoose
Definition: TauDefs.h:139
xAOD::TauJetParameters::cellBasedEnergyRing1
@ cellBasedEnergyRing1
EM+TES final scale.
Definition: TauDefs.h:246
xAOD::TauJetParameters::EleBDTMedium
@ EleBDTMedium
Definition: TauDefs.h:140
xAODTauJetAuxContainerCnv_v2::transToPersWithKey
virtual void transToPersWithKey(const xAOD::TauJetAuxContainer *, xAOD::TauJetAuxContainer_v2 *, const std::string &key, MsgStream &log) const override
Dummy function inherited from the base class.
Definition: xAODTauJetAuxContainerCnv_v2.cxx:298
IProxyDict
A proxy dictionary.
Definition: AthenaKernel/AthenaKernel/IProxyDict.h:47
xAOD::TauJetParameters::cellBasedEnergyRing6
@ cellBasedEnergyRing6
Ring 6: 0.15 < R < 0.2.
Definition: TauDefs.h:256
xAOD::AuxContainerBase::resize
virtual bool resize(size_t size) override
Resize the arrays to a given size.
Definition: AuxContainerBase.cxx:507
xAOD::TauJetAuxContainer_v2
Temporary container used until we have I/O for AuxStoreInternal.
Definition: TauJetAuxContainer_v2.h:39
xAOD::TauJet_v3::setIsTau
void setIsTau(TauJetParameters::IsTauFlag flag, bool value)
Set Flag for tau acceptance based on predefined arbitrary criteria.
Definition: TauJet_v3.cxx:252
xAOD::TauJetParameters::hadRadius
@ hadRadius
Get hadron calorimeter radius.
Definition: TauDefs.h:192
xAOD::TauJet_v3::setP4
void setP4(double pt, double eta, double phi, double m)
Set methods for IParticle values.
Definition: TauJet_v3.cxx:171
xAOD::TauJetParameters::ipSigLeadTrk
@ ipSigLeadTrk
Definition: TauDefs.h:160
xAOD::TauJetParameters::JetSeed
@ JetSeed
Definition: TauDefs.h:336
xAOD::TauJetParameters::stripWidth2
@ stripWidth2
Get strip width ^2.
Definition: TauDefs.h:202
xAOD::TauJetParameters::centFrac
@ centFrac
Get centrality fraction.
Definition: TauDefs.h:200
xAOD::TauJetParameters::dRmax
@ dRmax
Get maximal dR of tracks associated to calo-seeded tau.
Definition: TauDefs.h:226
xAOD::TauJetParameters::JetBDTSigLoose
@ JetBDTSigLoose
Definition: TauDefs.h:136
xAOD::AuxContainerBase::size
virtual size_t size() const override
Get the size of the container.
Definition: AuxContainerBase.cxx:432
xAOD::TrackParticle_v1::p4
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
Definition: TrackParticle_v1.cxx:130
xAOD::TauJetParameters::passTrkSelector
@ passTrkSelector
Definition: TauDefs.h:405
xAOD::TauJetParameters::etHadAtEMScale
@ etHadAtEMScale
Get Hadronic energy at EM scale.
Definition: TauDefs.h:196
xAOD::TauTrackAuxContainer
TauTrackAuxContainer_v1 TauTrackAuxContainer
Definition: TauTrackAuxContainer.h:14
xAOD::TauJetParameters::PanTau_BDTValue_1p1n_vs_1pXn
@ PanTau_BDTValue_1p1n_vs_1pXn
Definition: TauDefs.h:363
xAOD::TauJetParameters::topoMeanDeltaR
@ topoMeanDeltaR
get mean distance to leading topocluster for topoclusters of jet associated to tau candidate
Definition: TauDefs.h:181
xAOD::TauJetParameters::effTopoMeanDeltaR
@ effTopoMeanDeltaR
get mean distance to leading topocluster for effective topoclusters of jet associated to tau candidat...
Definition: TauDefs.h:183
lumiFormat.i
int i
Definition: lumiFormat.py:85
xAOD::TauJetParameters::caloIso
@ caloIso
Get sum of transvers energy of clusters around jet seed within 0.2 < dR < 0.4
Definition: TauDefs.h:222
xAOD::TauJetParameters::centFracCorrected
@ centFracCorrected
Definition: TauDefs.h:285
xAODTauJetAuxContainerCnv_v2.h
xAOD::TauTrackAuxContainer_v1
Definition: TauTrackAuxContainer_v1.h:21
xAOD::TauJetParameters::cellBasedEnergyRing2
@ cellBasedEnergyRing2
Ring 2: 0.05 < R < 0.075.
Definition: TauDefs.h:248
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
xAOD::TauJet_v3::setSecondaryVertexLink
void setSecondaryVertexLink(const VertexLink_t &vertexLink)
xAOD::TauJetParameters::trkAvgDist
@ trkAvgDist
Get calibrated EM transverse energy (DEPRECATED since r19)
Definition: TauDefs.h:214
xAOD::TauJet_v3
Class describing a tau jet.
Definition: TauJet_v3.h:41
xAOD::TauJetParameters::nChargedTracks
@ nChargedTracks
Definition: TauDefs.h:322
xAOD::TauJet_v3::setShotPFOLinks
void setShotPFOLinks(const PFOLinks_t &shotPFOs)
xAOD::TauJet_v3::setTrackFilterProngs
void setTrackFilterProngs(int)
xAOD::TauJetParameters::TESCalibConstant
@ TESCalibConstant
Definition: TauDefs.h:283
xAOD::TauJet_v3::addTauTrackLink
void addTauTrackLink(const ElementLink< TauTrackContainer > &tr)
add a TauTrack to the tau
Definition: TauJet_v3.cxx:523
xAOD::TauJet_v3::setTrackFilterQuality
void setTrackFilterQuality(int)
xAOD::TauJetParameters::etOverPtLeadTrk
@ etOverPtLeadTrk
Definition: TauDefs.h:158
xAOD::TauJetParameters::caloIsoCorrected
@ caloIsoCorrected
Get sum of transvers energy of clusters around jet seed within 0.2 < dR < 0.4 (pile up corrected)
Definition: TauDefs.h:224
xAOD::TauJet_v2
Class describing a tau jet.
Definition: TauJet_v2.h:38
xAOD::TauJetParameters::PanTau_BDTValue_1p0n_vs_1p1n
@ PanTau_BDTValue_1p0n_vs_1p1n
Definition: TauDefs.h:362
xAOD::TauJetParameters::coreTrack
@ coreTrack
Definition: TauDefs.h:403
xAOD::TauJetParameters::PanTauCellBased
@ PanTauCellBased
Definition: TauDefs.h:344
xAOD::TauJet
TauJet_v3 TauJet
Definition of the current "tau version".
Definition: Event/xAOD/xAODTau/xAODTau/TauJet.h:17
xAOD::TauJetParameters::SumPtTrkFracCorrected
@ SumPtTrkFracCorrected
Definition: TauDefs.h:290
TauJetContainer.h
xAOD::TauJetParameters::etOverPtLeadTrkCorrected
@ etOverPtLeadTrkCorrected
Definition: TauDefs.h:286
xAOD::TauJetParameters::numTopoClusters
@ numTopoClusters
get number of topocluster constituents of jet associated to tau candidate
Definition: TauDefs.h:173
DataVector::back
const T * back() const
Access the last element in the collection as an rvalue.
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
xAOD::TauJetParameters::EleBDTTight
@ EleBDTTight
Definition: TauDefs.h:141
xAOD::TauJetParameters::cellBasedEnergyRing5
@ cellBasedEnergyRing5
Ring 5: 0.125 < R < 0.15.
Definition: TauDefs.h:254
TauJetContainer_v2.h
xAOD::TauJetParameters::massTrkSys
@ massTrkSys
Definition: TauDefs.h:161
xAOD::TauJetParameters::TauEnergyScale
@ TauEnergyScale
Definition: TauDefs.h:339
xAOD::TauJetParameters::classifiedIsolation
@ classifiedIsolation
Definition: TauDefs.h:407
xAOD::TauJetParameters::JetBDTSigTight
@ JetBDTSigTight
Definition: TauDefs.h:138
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2
@ PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2
Definition: TauDefs.h:371
xAOD::TauJet_v3::setProtoChargedPFOLinks
void setProtoChargedPFOLinks(const PFOLinks_t &protoChargedPFOs)
xAOD::TauJetParameters::BDTJetScore
@ BDTJetScore
Boosted Decision Tree socre for Jet rejection (not transformed)
Definition: TauDefs.h:78
DataVector::push_back
value_type push_back(value_type pElem)
Add an element to the end of the collection.
xAOD::TauJetParameters::TauJetVtxFraction
@ TauJetVtxFraction
@Tau Jet Vertex Fraction
Definition: TauDefs.h:262
xAOD::TauJetParameters::nStrip
@ nStrip
Get number of strips.
Definition: TauDefs.h:204
SG::sgkey_t
uint32_t sgkey_t
Type used for hashed StoreGate key+CLID pairs.
Definition: CxxUtils/CxxUtils/sgkey_t.h:32
xAOD::TauJetParameters::TRT_NHT_OVER_NLT
@ TRT_NHT_OVER_NLT
TRT hits high threshold over low threshold.
Definition: TauDefs.h:260
xAOD::TauJet_v3::setChargedPFOLinks
void setChargedPFOLinks(const PFOLinks_t &chargedPFOs)
xAOD::TauJetParameters::cellBasedEnergyRing3
@ cellBasedEnergyRing3
Ring 3: 0.075 < R < 0.10.
Definition: TauDefs.h:250
xAOD::TauJetParameters::numEffTopoClusters
@ numEffTopoClusters
get number of effective topocluster constituents of jet associated to tau candidate !...
Definition: TauDefs.h:175
xAOD::TauJet_v3::setPi0ConeDR
void setPi0ConeDR(float)
xAOD::TauJetParameters::PanTau_DecayModeProto
@ PanTau_DecayModeProto
Definition: TauDefs.h:359
xAOD::TauJetParameters::BDTEleScore
@ BDTEleScore
Boosted Decision Tree score for electron rejection.
Definition: TauDefs.h:80
xAOD::TauJet_v3::setPi0PFOLinks
void setPi0PFOLinks(const PFOLinks_t &pi0PFOs)
xAOD::TauJetParameters::ipZ0SinThetaSigLeadTrk
@ ipZ0SinThetaSigLeadTrk
Definition: TauDefs.h:157
xAOD::TauJetParameters::numCells
@ numCells
Definition: TauDefs.h:171
xAOD::TauJetParameters::innerTrkAvgDist
@ innerTrkAvgDist
Definition: TauDefs.h:287
xAOD::TauJet_v3::setROIWord
void setROIWord(unsigned int)
xAOD::TauJet_v3::setVertexLink
void setVertexLink(const VertexLink_t &vertexLink)
xAOD::TauJetParameters::TESOffset
@ TESOffset
Definition: TauDefs.h:282
xAOD::TauJetParameters::PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts
@ PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts
Definition: TauDefs.h:368
xAOD::TauJetParameters::ChPiEMEOverCaloEME
@ ChPiEMEOverCaloEME
Definition: TauDefs.h:278
xAOD::TauJetParameters::sumEMCellEtOverLeadTrkPt
@ sumEMCellEtOverLeadTrkPt
Definition: TauDefs.h:231
xAOD::TauTrack_v1
Definition: TauTrack_v1.h:27
xAOD::TauJetParameters::modifiedIsolationTrack
@ modifiedIsolationTrack
Definition: TauDefs.h:412
DEBUG
#define DEBUG
Definition: page_access.h:11
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
xAOD::TauJetParameters::effTopoInvMass
@ effTopoInvMass
get invariant mass from effective topoclusters of jet associated to tau candidate
Definition: TauDefs.h:179
xAOD::TauJet_v3::setDetail
void setDetail(TauJetParameters::Detail detail, int value)
Definition: TauJet_v3.cxx:309
xAOD::TauJetParameters::innerTrkAvgDistCorrected
@ innerTrkAvgDistCorrected
Definition: TauDefs.h:288
xAOD::TauJetParameters::isolFrac
@ isolFrac
Get isolation fraction.
Definition: TauDefs.h:198
xAOD::TauJetParameters::wideTrack
@ wideTrack
Definition: TauDefs.h:404
xAOD::TauJetParameters::PanTau_isPanTauCandidate
@ PanTau_isPanTauCandidate
Definition: TauDefs.h:357
xAOD::track
@ track
Definition: TrackingPrimitives.h:513
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
xAOD::TauJetParameters::SumPtTrkFrac
@ SumPtTrkFrac
Definition: TauDefs.h:289
xAOD::TauJetParameters::DetectorAxis
@ DetectorAxis
Definition: TauDefs.h:337
xAODTauJetAuxContainerCnv_v2::persToTransWithKey
virtual void persToTransWithKey(const xAOD::TauJetAuxContainer_v2 *oldObj, xAOD::TauJetAuxContainer *newObj, const std::string &key, MsgStream &log) const override
Function converting from the old type to the current one.
Definition: xAODTauJetAuxContainerCnv_v2.cxx:33
xAOD::TauJetParameters::BDTJetScoreSigTrans
@ BDTJetScoreSigTrans
BDT score which is signal transformed/flattened.
Definition: TauDefs.h:84
xAOD::TauJetParameters::PSSFraction
@ PSSFraction
Definition: TauDefs.h:277
xAOD::TauJetParameters::EMRadius
@ EMRadius
Get E_T radius.
Definition: TauDefs.h:190
xAOD::TauJet_v3::setDiscriminant
void setDiscriminant(TauJetParameters::TauID discID, double disc)
Set value of discriminant.
Definition: TauJet_v3.cxx:237
StoreGateSvc.h
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
TrackParticleContainer.h
xAOD::TauJetParameters::EMPOverTrkSysP
@ EMPOverTrkSysP
Definition: TauDefs.h:279
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts
@ PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts
Definition: TauDefs.h:373
xAOD::TauJetParameters::etEMAtEMScale
@ etEMAtEMScale
Get EM energy at EM scale.
Definition: TauDefs.h:194
xAOD::TrackParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .)
xAOD::TauJetParameters::unclassified
@ unclassified
Definition: TauDefs.h:410
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1
@ PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1
Definition: TauDefs.h:370
ServiceHandle< StoreGateSvc >
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
xAOD::TauJetParameters::leadTrkPt
@ leadTrkPt
Definition: TauDefs.h:159
Generated on Mon Sep 29 2025 21:23:55 for ATLAS Offline Software by   doxygen 1.8.18