ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | Protected Attributes | Private Attributes | List of all members
xAODPrivate::Loader Class Reference
Inheritance diagram for xAODPrivate::Loader:
Collaboration diagram for xAODPrivate::Loader:

Public Member Functions

 Loader (xAOD::Event &event, const std::string &name, const std::type_info &ti)
void setProxy (SG::DataProxy &proxy)
virtual StatusCode createObj (IOpaqueAddress *, DataObject *&obj) override
virtual long repSvcType () const override
void addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > & tracks () const
void clear ()
void selectTracks (const TrigInDetTrackCollection *trigtracks)
void selectTracks (const Rec::TrackParticleContainer *trigtracks)
void selectTracks (const Analysis::MuonContainer *muontracks)

Protected Member Functions

double phiCorr (double phi)
void ipCorr (double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)

Protected Attributes

double m_beamX
double m_beamY
double m_beamZ
std::vector< TIDA::Track * > m_tracks

Private Attributes

xAOD::Event & m_event
std::string m_name
const std::type_info & m_ti
SG::DataProxy * m_proxy

Detailed Description

Definition at line 142 of file EventIProxyDict.cxx.

Constructor & Destructor Documentation

◆ Loader()

xAODPrivate::Loader::Loader ( xAOD::Event & event,
const std::string & name,
const std::type_info & ti )
inline

Definition at line 144 of file EventIProxyDict.cxx.

145 : Converter(0, CLID_NULL),
146 m_event(event),
147 m_name(name),
148 m_ti(ti),
149 m_proxy(nullptr) {}
xAODPrivate::Loader::m_proxy
SG::DataProxy * m_proxy
Definition EventIProxyDict.cxx:177
xAODPrivate::Loader::m_ti
const std::type_info & m_ti
Definition EventIProxyDict.cxx:176

Member Function Documentation

◆ addTrack()

void Converter::addTrack ( TIDA::Track * t)
inlineinherited

Definition at line 45 of file Converter.h.

45 {
46 m_tracks.push_back(t);
47 }
Converter::m_tracks
std::vector< TIDA::Track * > m_tracks
Definition Converter.h:321

◆ clear()

void Converter::clear ( )
inlineinherited

Definition at line 53 of file Converter.h.

53{m_tracks.clear();}

◆ createObj()

virtual StatusCode xAODPrivate::Loader::createObj ( IOpaqueAddress * ,
DataObject *& obj )
inlineoverridevirtual

Definition at line 153 of file EventIProxyDict.cxx.

153 {
154 static const bool SILENT = true;
155 static const bool METADATA = false;
156
157 // Try to find the object amongst the output objects first:
158 if (m_event.getOutputObject(m_name, m_ti, METADATA)) {
159 obj = new xAODPrivate::HolderBucket(m_name, m_ti, m_event);
160 return StatusCode::SUCCESS;
161 }
162 // If it's not on the output, try the input:
163 else if (m_event.getInputObject(m_name, m_ti, SILENT, METADATA)) {
164 obj = new xAODPrivate::HolderBucket(m_name, m_ti, m_event);
165 m_proxy->setConst();
166 return StatusCode::SUCCESS;
167 }
168
169 return StatusCode::FAILURE;
170 }
CP::OutOfValidityAction::SILENT
@ SILENT
don't print anything and return success
Definition OutOfValidityHelper.h:59

◆ ipCorr()

void Converter::ipCorr ( double d0,
double z0,
double & d0c,
double & z0c,
double phi0,
double eta,
double pt )
inlineprotectedinherited

Definition at line 257 of file Converter.h.

257 {
258
259 double sn = sin(phi0);
260 double cs = cos(phi0);
261 double sd0 = (d0 != 0 ? d0/fabs(d0) : 1);
262 double spt = (pt != 0 ? pt/fabs(pt) : 1);
263
264 if (fabs(pt) >= 1*CLHEP::TeV) {
265
266 d0c = d0 + m_beamX*sn - m_beamY*cs;
267
268 } else {
269
270 double rc = fabs(pt)*15.0/(9.0*1.042);
271
272 double xc = (fabs(d0)-spt*sd0*rc)*cos(phi0+M_PI/2*sd0) - m_beamX;
273 double yc = (fabs(d0)-spt*sd0*rc)*sin(phi0+M_PI/2*sd0) - m_beamY;
274
275 double newphi;
276 double xd01,yd01,xd02,yd02;
277
278 if (xc == 0) {
279 xd01 = 0; yd01 = rc+yc;
280 xd02 = 0; yd02 = -rc+yc;
281 } else {
282 xd01 = xc+yc/xc*yc+sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd01 = yc/xc*xd01;
283 xd02 = xc+yc/xc*yc-sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd02 = yc/xc*xd02;
284 }
285
286 double r1 = sqrt(xd01*xd01+yd01*yd01);
287 double r2 = sqrt(xd02*xd02+yd02*yd02);
288
289 double phiV;
290
291 if (r1 < r2)
292 phiV = atan2(yd01,xd01);
293 else
294 phiV = atan2(yd02,xd02);
295
296 double phi1 = phiCorr(phiV+M_PI/2);
297 double phi2 = phiCorr(phiV-M_PI/2);
298
299 if (fabs(phiCorr(phi1-phi0))<=fabs(phiCorr(phi2-phi0)))
300 newphi = phi1;
301 else
302 newphi = phi2;
303
304 d0c = fabs(sqrt(xc*xc+yc*yc)-rc)*sin(phiV-newphi);
305
306 double theta=2*atan2(exp(-eta),1);
307 double theta_save=theta;
308 theta = theta - (1+spt)/2*M_PI;
309 if (theta>0) theta = theta_save;
310
311 double deltaz0= -spt*rc/tan(theta)*phiCorr(phiCorr(newphi)-phiCorr(phi0));
312 z0c = z0 + deltaz0;
313
314 }
315 }
M_PI
#define M_PI
Definition ActiveFraction.h:11
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
theta
Scalar theta() const
theta method
Definition AmgMatrixBasePlugin.h:75
rc
static Double_t rc
Definition LArPhysWaveHECTool.cxx:37
pow
constexpr int pow(int base, int exp) noexcept
Definition ap_fixedTest.cxx:15
Converter::m_beamX
double m_beamX
Definition Converter.h:320
Converter::phiCorr
double phiCorr(double phi)
Definition Converter.h:248
Converter::m_beamY
double m_beamY
Definition Converter.h:320
test_pyathena.pt
pt
Definition test_pyathena.py:11

◆ phiCorr()

double Converter::phiCorr ( double phi)
inlineprotectedinherited

Definition at line 248 of file Converter.h.

248 {
249
250 if (phi < -TMath::Pi()) phi += 2*TMath::Pi();
251 if (phi > TMath::Pi()) phi -= 2*TMath::Pi();
252
253 return phi;
254 }
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67

◆ repSvcType()

virtual long xAODPrivate::Loader::repSvcType ( ) const
inlineoverridevirtual

Definition at line 171 of file EventIProxyDict.cxx.

171{ return 0; }

◆ selectTracks() [1/3]

void Converter::selectTracks ( const Analysis::MuonContainer * muontracks)
inlineinherited

Definition at line 177 of file Converter.h.

177 {
178
179 static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
180
181 Analysis::MuonContainer::const_iterator muonitr = muontracks->begin();
182 Analysis::MuonContainer::const_iterator muonend = muontracks->end();
183
184 while(muonitr!=muonend) {
185
186 bool hasInDet = (*muonitr)->hasInDetTrackParticle();
187 const Rec::TrackParticle *trackitr=NULL;
188 if(hasInDet) trackitr=(*muonitr)->inDetTrackParticle();
189 muonitr++;
190 if(!hasInDet) continue;
191
192 const Trk::MeasuredPerigee* measPer = trackitr->measuredPerigee();
193
194 CLHEP::HepVector perigeeParams = measPer->parameters();
195 double pT = measPer->pT();
196 double eta = measPer->eta();
197 double phi = perigeeParams[Trk::phi0];
198 double z0 = perigeeParams[Trk::z0];
199 double d0 = perigeeParams[Trk::d0];
200 //double theta = perigeeParams[Trk::theta];
201
202 double deta = 1;
203 double dphi = 1;
204 double dz0 = 1;
205 double dd0 = 1;
206 double dpT = 1;
207
208 // Check number of hits
209 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
210 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
211 // stereo clusters making a spacepoint are two "hits"
212 const Trk::TrackSummary *summary = trackitr->trackSummary();
213 int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
214 int nPixelHits = 2*summary->get(Trk::numberOfPixelHits);
215 int nSctHits = summary->get(Trk::numberOfSCTHits);
216 int nStrawHits = summary->get(Trk::numberOfTRTHits);
217 int nTrHits = summary->get(Trk::numberOfTRTHighThresholdHits);
218
219 int nSiHits = nPixelHits + nSctHits;
220
221 const Trk::FitQuality *quality = trackitr->fitQuality();
222 double chi2 = quality->chiSquared();
223 double dof = quality->quality->numberDoF();
224
225 unsigned bitmap = 0;
226
227 for ( int ih=0 ; ih<20 ; ih++ ) {
228 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
229 }
230
231 // Shift d0 and z0 according to beam position
232 ipCorr(d0, z0, d0, z0, phi, eta, pT);
233
234 // Create and save Track
235 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
236 deta, dphi, dz0, dd0, dpT,
237 nBlayerHits, nPixelHits, nSctHits, nSiHits,
238 nStrawHits, nTrHits, bitmap, 0,
239 -1) ;
240
241 addTrack(t);
242 }
243
244 }
Converter::addTrack
void addTrack(TIDA::Track *t)
Definition Converter.h:45
Converter::ipCorr
void ipCorr(double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)
Definition Converter.h:257
DataVector< Analysis::Muon >::const_iterator
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
DataVector::end
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
DataVector::begin
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
Rec::TrackParticle::measuredPerigee
const Trk::Perigee * measuredPerigee() const
Accessor method for Perigee.
Definition Reconstruction/Particle/Particle/TrackParticle.h:180
Trk::FitQuality::numberDoF
int numberDoF() const
returns the number of degrees of freedom of the overall track or vertex fit as integer
Definition FitQuality.h:60
Trk::FitQuality::chiSquared
double chiSquared() const
returns the of the overall track fit
Definition FitQuality.h:56
Trk::ParametersBase::eta
double eta() const
Access method for pseudorapidity - from momentum.
Trk::ParametersBase::pT
double pT() const
Access method for transverse momentum.
Trk::TrackParticleBase::fitQuality
const FitQuality * fitQuality() const
accessor function for FitQuality.
Definition TrackParticleBase.h:259
Trk::TrackParticleBase::trackSummary
const TrackSummary * trackSummary() const
accessor function for TrackSummary.
Definition TrackParticleBase.h:247
chi2
double chi2(TH1 *h0, TH1 *h1)
Definition comparitor.cxx:525
IDTPM::nSiHits
float nSiHits(const U &p)
Definition TrackParametersHelper.h:427
IDTPM::nPixelHits
float nPixelHits(const U &p)
Definition TrackParametersHelper.h:354
SCT_Monitoring::summary
@ summary
Definition SCT_MonitoringNumbers.h:65
Trk::phi0
@ phi0
Definition ParamDefs.h:65
Trk::d0
@ d0
Definition ParamDefs.h:63
Trk::z0
@ z0
Definition ParamDefs.h:64
Trk::DetectorType
DetectorType
enumerates the various detector types currently accessible from the isHit() method.
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:231
Trk::numberOfPixelHits
@ numberOfPixelHits
number of pixel layers on track with absence of hits
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:57
Trk::numberOfBLayerHits
@ numberOfBLayerHits
these are the hits in the 0th pixel layer?
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:51
Trk::numberOfTRTHighThresholdHits
@ numberOfTRTHighThresholdHits
total number of TRT hits which pass the high threshold
Definition Tracking/TrkEvent/TrkTrackSummary/TrkTrackSummary/TrackSummary.h:87
Trk::MeasuredPerigee
Perigee MeasuredPerigee
Definition MeasuredPerigeeCnv_p1.h:23
read_hist_ntuple.t
t
Definition read_hist_ntuple.py:5

◆ selectTracks() [2/3]

void Converter::selectTracks ( const Rec::TrackParticleContainer * trigtracks)
inlineinherited

Definition at line 111 of file Converter.h.

111 {
112
113 Rec::TrackParticleContainer::const_iterator trackitr = trigtracks->begin();
114 Rec::TrackParticleContainer::const_iterator trackend = trigtracks->end();
115
116 static int hpmap[20] = { 0, 1, 2, 7, 8, 9, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
117
118 while ( trackitr!=trackend ) {
119
120 const Trk::MeasuredPerigee* measPer = (*trackitr)->measuredPerigee();
121
122 CLHEP::HepVector perigeeParams = measPer->parameters();
123 double pT = measPer->pT();
124 double eta = measPer->eta();
125 double phi = perigeeParams[Trk::phi0];
126 double z0 = perigeeParams[Trk::z0];
127 double d0 = perigeeParams[Trk::d0];
128 //double theta = perigeeParams[Trk::theta];
129
130 double deta = 1;
131 double dphi = 1;
132 double dz0 = 1;
133 double dd0 = 1;
134 double dpT = 1;
135
136 // Check number of hits
137 // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really
138 // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
139 // stereo clusters making a spacepoint are two "hits"
140 const Trk::TrackSummary *summary = (*trackitr)->trackSummary();
141 int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits);
142 int nPixelHits = 2*summary->get(Trk::numberOfPixelHits);
143 int nSctHits = summary->get(Trk::numberOfSCTHits);
144 int nStrawHits = summary->get(Trk::numberOfTRTHits);
145 int nTrHits = summary->get(Trk::numberOfTRTHighThresholdHits);
146 int nSiHits = nPixelHits + nSctHits;
147 bool expectBL = false; // Not stored for Rec::TrackParticle
148
149 const Trk::FitQuality *quality = (*trackitr)->fitQuality();
150 double chi2 = quality->chiSquared();
151 double dof = quality->numberDoF();
152
153 unsigned bitmap = 0;
154
155 for ( int ih=0 ; ih<20 ; ih++ ) {
156 if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );
157 }
158
159 // Shift d0 and z0 according to beam position
160 ipCorr(d0, z0, d0, z0, phi, eta, pT);
161
162 // Create and save Track
163 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
164 deta, dphi, dz0, dd0, dpT,
165 nBlayerHits, nPixelHits, nSctHits, nSiHits,
166 nStrawHits, nTrHits, bitmap, 0,
167 -1,
168 expectBL) ;
169
170 addTrack(t);
171 trackitr++;
172 }
173
174 }

◆ selectTracks() [3/3]

void Converter::selectTracks ( const TrigInDetTrackCollection * trigtracks)
inlineinherited

Definition at line 56 of file Converter.h.

56 {
57
58 TrigInDetTrackCollection::const_iterator trackitr = trigtracks->begin();
59 TrigInDetTrackCollection::const_iterator trackend = trigtracks->end();
60
61 while ( trackitr!=trackend ) {
62
63 double eta = (*trackitr)->param()->eta();
64 double phi = (*trackitr)->param()->phi0();
65 double z0 = (*trackitr)->param()->z0();
66 double pT = (*trackitr)->param()->pT();
67 double d0 = (*trackitr)->param()->a0();
68 //double theta = 2*atan2(exp(-(*trackitr)->param()->eta()),1);
69
70 double deta = (*trackitr)->param()->eeta();
71 double dphi = (*trackitr)->param()->ephi0();
72 double dz0 = (*trackitr)->param()->ez0();
73 double dpT = (*trackitr)->param()->epT();
74 double dd0 = (*trackitr)->param()->ea0();
75
76 int algoid = (*trackitr)->algorithmId();
77
78 int nBlayerHits = ((*trackitr)->HitPattern() & 0x1);
79 int nPixelHits = 2*(*trackitr)->NPixelSpacePoints(); // NB: for comparison with offline
80 int nSctHits = 2*(*trackitr)->NSCT_SpacePoints(); // a spacepoint is 2 "hits"
81 int nStrawHits = (*trackitr)->NStrawHits();
82 int nTrHits = (*trackitr)->NTRHits();
83
84 int nSiHits = nPixelHits + nSctHits;
85 bool expectBL = false; //not available with TrigInDetTrack
86
87 unsigned hitPattern = (*trackitr)->HitPattern();
88 unsigned multiPattern = 0;
89
90 double chi2 = (*trackitr)->chi2();
91 double dof = 0;
92
93 // Shift d0 and z0 according to beam position
94 ipCorr(d0, z0, d0, z0, phi, eta, pT);
95
96 // Create and save Track
97 TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
98 deta, dphi, dz0, dd0, dpT,
99 nBlayerHits, nPixelHits, nSctHits, nSiHits,
100 nStrawHits, nTrHits, hitPattern, multiPattern,
101 algoid,
102 expectBL ) ;
103
104 addTrack(t);
105 trackitr++;
106 }
107 }

◆ setProxy()

void xAODPrivate::Loader::setProxy ( SG::DataProxy & proxy)
inline

Definition at line 151 of file EventIProxyDict.cxx.

151{ m_proxy = &proxy; }
StateLessPT_NewConfig.proxy
proxy
Definition StateLessPT_NewConfig.py:407

◆ tracks()

const std::vector< TIDA::Track * > & Converter::tracks ( ) const
inlineinherited

Definition at line 50 of file Converter.h.

50{return m_tracks;}

Member Data Documentation

◆ m_beamX

double Converter::m_beamX
protectedinherited

Definition at line 320 of file Converter.h.

◆ m_beamY

double Converter::m_beamY
protectedinherited

Definition at line 320 of file Converter.h.

◆ m_beamZ

double Converter::m_beamZ
protectedinherited

Definition at line 320 of file Converter.h.

◆ m_event

xAOD::Event& xAODPrivate::Loader::m_event
private

Definition at line 174 of file EventIProxyDict.cxx.

◆ m_name

std::string xAODPrivate::Loader::m_name
private

Definition at line 175 of file EventIProxyDict.cxx.

◆ m_proxy

SG::DataProxy* xAODPrivate::Loader::m_proxy
private

Definition at line 177 of file EventIProxyDict.cxx.

◆ m_ti

const std::type_info& xAODPrivate::Loader::m_ti
private

Definition at line 176 of file EventIProxyDict.cxx.

◆ m_tracks

std::vector<TIDA::Track*> Converter::m_tracks
protectedinherited

Definition at line 321 of file Converter.h.

The documentation for this class was generated from the following file:
Generated on for ATLAS Offline Software by doxygen 1.14.0