ATLAS Offline Software
Loading...
Searching...
No Matches
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
30
31
32void xAODTauJetAuxContainerCnv_v2::
33persToTransWithKey( 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
297void xAODTauJetAuxContainerCnv_v2::
298transToPersWithKey( 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
endmsg
#define endmsg
Definition AnalysisConfig_Ntuple.cxx:63
TauJetContainer_v2.h
DataVector::back
const T * back() const
Access the last element in the collection as an rvalue.
DataVector::push_back
value_type push_back(value_type pElem)
Add an element to the end of the collection.
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
xAODTauJetAuxContainerCnv_v2::xAODTauJetAuxContainerCnv_v2
xAODTauJetAuxContainerCnv_v2()
Default constructor.
Definition xAODTauJetAuxContainerCnv_v2.cxx:27
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
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::AuxContainerBase::resize
virtual bool resize(size_t size) override
Resize the arrays to a given size.
Definition AuxContainerBase.cxx:507
xAOD::AuxContainerBase::size
virtual size_t size() const override
Get the size of the container.
Definition AuxContainerBase.cxx:432
xAOD::TauJetAuxContainer_v2
Temporary container used until we have I/O for AuxStoreInternal.
Definition TauJetAuxContainer_v2.h:39
xAOD::TauJet_v2
Class describing a tau jet.
Definition TauJet_v2.h:38
xAOD::TauJet_v3::setShotPFOLinks
void setShotPFOLinks(const PFOLinks_t &shotPFOs)
xAOD::TauJet_v3::setProtoPi0PFOLinks
void setProtoPi0PFOLinks(const PFOLinks_t &protoPi0PFOs)
xAOD::TauJet_v3::setJetLink
void setJetLink(const JetLink_t &jetLink)
xAOD::TauJet_v3::setPi0PFOLinks
void setPi0PFOLinks(const PFOLinks_t &pi0PFOs)
xAOD::TauJet_v3::setROIWord
void setROIWord(unsigned int)
xAOD::TauJet_v3::setPanTauDetail
void setPanTauDetail(TauJetParameters::PanTauDetails panTauDetail, int value)
Definition TauJet_v3.cxx:367
xAOD::TauJet_v3::setCharge
void setCharge(float)
xAOD::TauJet_v3::setVertexLink
void setVertexLink(const VertexLink_t &vertexLink)
xAOD::TauJet_v3::setPi0ConeDR
void setPi0ConeDR(float)
xAOD::TauJet_v3::setProtoChargedPFOLinks
void setProtoChargedPFOLinks(const PFOLinks_t &protoChargedPFOs)
xAOD::TauJet_v3::setSecondaryVertexLink
void setSecondaryVertexLink(const VertexLink_t &vertexLink)
xAOD::TauJet_v3::setDiscriminant
void setDiscriminant(TauJetParameters::TauID discID, double disc)
Set value of discriminant.
Definition TauJet_v3.cxx:237
xAOD::TauJet_v3::setDetail
void setDetail(TauJetParameters::Detail detail, int value)
Definition TauJet_v3.cxx:309
xAOD::TauJet_v3::setTrackFilterQuality
void setTrackFilterQuality(int)
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::TauJet_v3::setTrackFilterProngs
void setTrackFilterProngs(int)
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::TauJet_v3::setProtoNeutralPFOLinks
void setProtoNeutralPFOLinks(const PFOLinks_t &protoNeutralPFOs)
xAOD::TauJet_v3::setNeutralPFOLinks
void setNeutralPFOLinks(const PFOLinks_t &neutralPFOs)
xAOD::TauJet_v3::addTauTrackLink
void addTauTrackLink(const ElementLink< TauTrackContainer > &tr)
add a TauTrack to the tau
Definition TauJet_v3.cxx:523
xAOD::TauJet_v3::setChargedPFOLinks
void setChargedPFOLinks(const PFOLinks_t &chargedPFOs)
xAOD::TauJet_v3::nTracks
size_t nTracks(TauJetParameters::TauTrackFlag flag=TauJetParameters::TauTrackFlag::classifiedCharged) const
Definition TauJet_v3.cxx:488
xAOD::TrackParticle_v1::m
virtual double m() const override final
The invariant mass of the particle..
Definition TrackParticle_v1.cxx:85
xAOD::TrackParticle_v1::p4
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
Definition TrackParticle_v1.cxx:122
xAOD::TrackParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .)
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition TrackParticle_v1.cxx:75
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TrackParticle_v1.cxx:79
SG::sgkey_t
uint32_t sgkey_t
Type used for hashed StoreGate key+CLID pairs.
Definition sgkey_t.h:32
xAOD::TauJetParameters::hadRadius
@ hadRadius
Get hadron calorimeter radius.
Definition TauDefs.h:192
xAOD::TauJetParameters::EMRadius
@ EMRadius
Get E_T radius.
Definition TauDefs.h:190
xAOD::TauJetParameters::secMaxStripEt
@ secMaxStripEt
migrate only seedTrk_ variables which are used in reco and ID and without prefix
Definition TauDefs.h:230
xAOD::TauJetParameters::TauJetVtxFraction
@ TauJetVtxFraction
@Tau Jet Vertex Fraction
Definition TauDefs.h:262
xAOD::TauJetParameters::stripWidth2
@ stripWidth2
Get strip width ^2.
Definition TauDefs.h:202
xAOD::TauJetParameters::topoInvMass
@ topoInvMass
get invariant mass from topoclusters of jet associated to tau candidate
Definition TauDefs.h:177
xAOD::TauJetParameters::etHadAtEMScale
@ etHadAtEMScale
Get Hadronic energy at EM scale.
Definition TauDefs.h:196
xAOD::TauJetParameters::isolFrac
@ isolFrac
Get isolation fraction.
Definition TauDefs.h:198
xAOD::TauJetParameters::lead2ClusterEOverAllClusterE
@ lead2ClusterEOverAllClusterE
Get sum of 2 leading cluster energy over all cluster energy.
Definition TauDefs.h:218
xAOD::TauJetParameters::topoMeanDeltaR
@ topoMeanDeltaR
get mean distance to leading topocluster for topoclusters of jet associated to tau candidate
Definition TauDefs.h:181
xAOD::TauJetParameters::trkRmsDist
@ trkRmsDist
Get the RMS of track distance to calorimeter seed.
Definition TauDefs.h:216
xAOD::TauJetParameters::nStrip
@ nStrip
Get number of strips.
Definition TauDefs.h:204
xAOD::TauJetParameters::cellBasedEnergyRing4
@ cellBasedEnergyRing4
Ring 4: 0.10 < R < 0.125.
Definition TauDefs.h:252
xAOD::TauJetParameters::effTopoMeanDeltaR
@ effTopoMeanDeltaR
get mean distance to leading topocluster for effective topoclusters of jet associated to tau candidat...
Definition TauDefs.h:183
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::trkAvgDist
@ trkAvgDist
Get calibrated EM transverse energy (DEPRECATED since r19)
Definition TauDefs.h:214
xAOD::TauJetParameters::cellBasedEnergyRing6
@ cellBasedEnergyRing6
Ring 6: 0.15 < R < 0.2.
Definition TauDefs.h:256
xAOD::TauJetParameters::cellBasedEnergyRing5
@ cellBasedEnergyRing5
Ring 5: 0.125 < R < 0.15.
Definition TauDefs.h:254
xAOD::TauJetParameters::effTopoInvMass
@ effTopoInvMass
get invariant mass from effective topoclusters of jet associated to tau candidate
Definition TauDefs.h:179
xAOD::TauJetParameters::cellBasedEnergyRing7
@ cellBasedEnergyRing7
Ring 7: 0.2 < R < 0.4.
Definition TauDefs.h:258
xAOD::TauJetParameters::etEMAtEMScale
@ etEMAtEMScale
Get EM energy at EM scale.
Definition TauDefs.h:194
xAOD::TauJetParameters::TRT_NHT_OVER_NLT
@ TRT_NHT_OVER_NLT
TRT hits high threshold over low threshold.
Definition TauDefs.h:260
xAOD::TauJetParameters::centFrac
@ centFrac
Get centrality fraction.
Definition TauDefs.h:200
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::TauJetParameters::cellBasedEnergyRing2
@ cellBasedEnergyRing2
Ring 2: 0.05 < R < 0.075.
Definition TauDefs.h:248
xAOD::TauJetParameters::numTopoClusters
@ numTopoClusters
get number of topocluster constituents of jet associated to tau candidate
Definition TauDefs.h:173
xAOD::TauJetParameters::numEffTopoClusters
@ numEffTopoClusters
get number of effective topocluster constituents of jet associated to tau candidate !...
Definition TauDefs.h:175
xAOD::TauJetParameters::cellBasedEnergyRing3
@ cellBasedEnergyRing3
Ring 3: 0.075 < R < 0.10.
Definition TauDefs.h:250
xAOD::TauJetParameters::cellBasedEnergyRing1
@ cellBasedEnergyRing1
EM+TES final scale.
Definition TauDefs.h:246
xAOD::TauJetParameters::dRmax
@ dRmax
Get maximal dR of tracks associated to calo-seeded tau.
Definition TauDefs.h:226
xAOD::TauJetParameters::lead3ClusterEOverAllClusterE
@ lead3ClusterEOverAllClusterE
Get sum of 3 leading cluster energy over all cluster energy.
Definition TauDefs.h:220
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2
@ PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2
Definition TauDefs.h:371
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts
@ PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts
Definition TauDefs.h:372
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1
@ PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1
Definition TauDefs.h:370
xAOD::TauJetParameters::PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged
@ PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged
Definition TauDefs.h:375
xAOD::TauJetParameters::PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts
@ PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts
Definition TauDefs.h:373
xAOD::TauJetParameters::PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts
@ PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts
Definition TauDefs.h:368
xAOD::TauJetParameters::BDTJetScore
@ BDTJetScore
Boosted Decision Tree socre for Jet rejection (not transformed)
Definition TauDefs.h:78
xAOD::TauJetParameters::BDTEleScore
@ BDTEleScore
Boosted Decision Tree score for electron rejection.
Definition TauDefs.h:80
xAOD::TauJetParameters::BDTJetScoreSigTrans
@ BDTJetScoreSigTrans
BDT score which is signal transformed/flattened.
Definition TauDefs.h:84
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::TauTrack
TauTrack_v1 TauTrack
Definition of the current version.
Definition TauTrack.h:16
xAOD::TauJet
TauJet_v3 TauJet
Definition of the current "tau version".
Definition Event/xAOD/xAODTau/xAODTau/TauJet.h:17
xAOD::TauJetAuxContainer
TauJetAuxContainer_v3 TauJetAuxContainer
Definition of the current taujet auxiliary container.
Definition TauJetAuxContainer.h:22
xAOD::TauTrackContainer
TauTrackContainer_v1 TauTrackContainer
Definition of the current TauTrack container version.
Definition TauTrackContainer.h:16
xAOD::TauJetContainer_v2
DataVector< xAOD::TauJet_v2 > TauJetContainer_v2
The container is a simple typedef for now.
Definition TauJetContainer_v2.h:19
xAOD::TauJetContainer
TauJetContainer_v3 TauJetContainer
Definition of the current "taujet container version".
Definition Event/xAOD/xAODTau/xAODTau/TauJetContainer.h:17
xAOD::TauTrackAuxContainer
TauTrackAuxContainer_v1 TauTrackAuxContainer
Definition TauTrackAuxContainer.h:14
xAODTauJetAuxContainerCnv_v2.h
Generated on for ATLAS Offline Software by doxygen 1.14.0