ATLAS Offline Software
AnalysisR3_Tier0.cxx
Go to the documentation of this file.
1 
11 #include <cmath>
12 
16 
18 
20  double /*pTCut*/,
21  double /*etaCut*/,
22  double /*d0Cut*/,
23  double /*z0Cut*/)
25  // m_pTCut(pTCut), /// left in commented for development purposes
26  // m_etaCut(etaCut),
27  // m_d0Cut(d0Cut),
28  // m_z0Cut(z0Cut),
29  m_debug(false),
30  m_eventid(0),
31  m_vtxanal(0),
32  m_monTool(0)
33 {}
34 
36  if ( m_monTool ) initialise_R3();
37 }
38 
40 
42 
43  ChainString cname = name();
44 
45  // std::cout << "AnalysisR3_Tier0::initialise() name " << name() << std::endl;
46 
47 #if 0
48  std::cout << "\nAnalysisR3_Tier0:: chain specification: " << cname << "\t" << cname.raw() << std::endl;
50  std::cout << "\tchain: " << cname.head() << std::endl;
51  std::cout << "\tkey: " << cname.tail() << std::endl;
52  std::cout << "\troi: " << cname.roi() << std::endl;
53  std::cout << "\tvtx: " << cname.vtx() << std::endl;
54  std::cout << "\tte: " << cname.element() << std::endl;
55 #endif
56 
57  m_debug = false;
58 
59  m_hchain = new TIDA::Histogram<float>( monTool(), "Chain" );
60 
61  m_hroieta = new TIDA::Histogram<float>( monTool(), "roi_eta" );
62 
65 
66  m_hntrk = new TIDA::Histogram<float>( monTool(), "reftrk_N" );
67 
69 
70  m_htrkpT = new TIDA::Histogram<float>( monTool(), "reftrk_pT" );
71  m_htrkphi = new TIDA::Histogram<float>( monTool(), "reftrk_phi" );
72  m_htrketa = new TIDA::Histogram<float>( monTool(), "reftrk_eta" );
73  if (name().find("LRT")!=std::string::npos || name().find("lrt")!=std::string::npos) {
74  m_htrkd0 = new TIDA::Histogram<float>( monTool(), "reftrk_d0" );
75  } else {
76  m_htrkd0 = new TIDA::Histogram<float>( monTool(), "reftrk_d0" );
77  }
78  m_htrkz0 = new TIDA::Histogram<float>( monTool(), "reftrk_z0" );
79 
81  m_htrkdd0 = new TIDA::Histogram<float>( monTool(), "reftrk_dd0" );
82  m_htrkdz0 = new TIDA::Histogram<float>( monTool(), "reftrk_dz0" );
83 
84  m_htrkd0sig = new TIDA::Histogram<float>( monTool(), "reftrk_d0sig" );
85 
86 
88 
91  // m_hntrk_rec = new TIDA::Histogram<float>( monTool(), "testtrk_N", "Test tracks", 100, vnbins );
92  m_hntrk_rec = new TIDA::Histogram<float>( monTool(), "testtrk_N" );
93 
94 
95  // m_htrkpT_rec = new TIDA::Histogram<float>( monTool(), "testtrk_pT" , "Test track pT", 25, 0., 100.);
96  m_htrkpT_rec = new TIDA::Histogram<float>( monTool(), "testtrk_pT" );
97  m_htrkphi_rec = new TIDA::Histogram<float>( monTool(), "testtrk_phi" );
98  m_htrketa_rec = new TIDA::Histogram<float>( monTool(), "testtrk_eta" );
99  if (name().find("LRT")!=std::string::npos || name().find("lrt")!=std::string::npos) {
100  m_htrkd0_rec = new TIDA::Histogram<float>( monTool(), "testtrk_d0" );
101  } else {
102  m_htrkd0_rec = new TIDA::Histogram<float>( monTool(), "testtrk_d0" );
103  }
104  m_htrkz0_rec = new TIDA::Histogram<float>( monTool(), "testtrk_z0" );
105 
106  m_htrkdd0_rec = new TIDA::Histogram<float>( monTool(), "testtrk_dd0" );
107  m_htrkdz0_rec = new TIDA::Histogram<float>( monTool(), "testtrk_dz0" );
108 
109  m_htrkd0sig_rec = new TIDA::Histogram<float>( monTool(), "testtrk_d0sig" );
110 
111 
112 
114 
115 
116  m_htotal_efficiency = new TIDA::Histogram<float>( monTool(), "Eff_overall" );
117 
118  m_hpTeff = new TIDA::Histogram<float>( monTool(), "Eff_pT" );
119  m_hetaeff = new TIDA::Histogram<float>( monTool(), "Eff_Eta" );
120  m_hphieff = new TIDA::Histogram<float>( monTool(), "Eff_Phi" );
121  if (name().find("LRT")!=std::string::npos || name().find("lrt")!=std::string::npos) {
122  m_hd0eff = new TIDA::Histogram<float>( monTool(), "Eff_d0" );
123  } else {
124  m_hd0eff = new TIDA::Histogram<float>( monTool(), "Eff_d0" );
125  }
126  m_hz0eff = new TIDA::Histogram<float>( monTool(), "Eff_z0" );
127  m_hnVtxeff = new TIDA::Histogram<float>( monTool(), "Eff_nVtx" );
128 
129 
130  m_hlbeff = new TIDA::Histogram<float>( monTool(), "Eff_lb" );
131 
132 
133  m_htrkvtx_x_lb = new TIDA::Histogram<float>( monTool(), "trkvtx_x_vs_lb" );
134  m_htrkvtx_y_lb = new TIDA::Histogram<float>( monTool(), "trkvtx_y_vs_lb" );
135  m_htrkvtx_z_lb = new TIDA::Histogram<float>( monTool(), "trkvtx_z_vs_lb" );
136 
139  // m_hotrkvtx_x_lb = new TIDA::Histogram<float>( monTool(), "otrkvtx_x_vs_lb" );
140  // m_hotrkvtx_y_lb = new TIDA::Histogram<float>( monTool(), "otrkvtx_y_vs_lb" );
141  // m_hotrkvtx_z_lb = new TIDA::Histogram<float>( monTool(), "otrkvtx_z_vs_lb" );
142 
143 
145  m_hnpixvseta = new TIDA::Histogram<float>( monTool(), "npix_vs_eta" );
146  m_hnpixvseta_rec = new TIDA::Histogram<float>( monTool(), "npix_vs_eta_rec" );
147 
148  m_hnsctvseta = new TIDA::Histogram<float>( monTool(), "nsct_vs_eta" );
149  m_hnsctvseta_rec = new TIDA::Histogram<float>( monTool(), "nsct_vs_eta_rec" );
150 
151  m_hntrtvseta = new TIDA::Histogram<float>( monTool(), "ntrt_vs_eta" );
152  m_hntrtvseta_rec = new TIDA::Histogram<float>( monTool(), "ntrt_vs_eta_rec" );
153 
154  m_hnpixvsphi = new TIDA::Histogram<float>( monTool(), "npix_vs_phi" );
155  m_hnpixvsphi_rec = new TIDA::Histogram<float>( monTool(), "npix_vs_phi_rec" );
156 
157  m_hnsctvsphi = new TIDA::Histogram<float>( monTool(), "nsct_vs_phi" );
158  m_hnsctvsphi_rec = new TIDA::Histogram<float>( monTool(), "nsct_vs_phi_rec" );
159 
160  m_hntrtvsphi = new TIDA::Histogram<float>( monTool(), "ntrt_vs_phi" );
161  m_hntrtvsphi_rec = new TIDA::Histogram<float>( monTool(), "ntrt_vs_phi_rec" );
162 
163  if (name().find("LRT")!=std::string::npos || name().find("lrt")!=std::string::npos) {
164  m_hnpixvsd0 = new TIDA::Histogram<float>( monTool(), "npix_vs_d0" );
165  m_hnpixvsd0_rec = new TIDA::Histogram<float>( monTool(), "npix_vs_d0_rec" );
166 
167  m_hnsctvsd0 = new TIDA::Histogram<float>( monTool(), "nsct_vs_d0" );
168  m_hnsctvsd0_rec = new TIDA::Histogram<float>( monTool(), "nsct_vs_d0_rec" );
169  } else {
170  m_hnpixvsd0 = new TIDA::Histogram<float>( monTool(), "npix_vs_d0" );
171  m_hnpixvsd0_rec = new TIDA::Histogram<float>( monTool(), "npix_vs_d0_rec" );
172 
173  m_hnsctvsd0 = new TIDA::Histogram<float>( monTool(), "nsct_vs_d0" );
174  m_hnsctvsd0_rec = new TIDA::Histogram<float>( monTool(), "nsct_vs_d0_rec" );
175  }
176 
177  m_hnpixvspT = new TIDA::Histogram<float>( monTool(), "npix_vs_pT" );
178  m_hnpixvspT_rec = new TIDA::Histogram<float>( monTool(), "npix_vs_pT_rec" );
179 
180  m_hnsctvspT = new TIDA::Histogram<float>( monTool(), "nsct_vs_pT" );
181  m_hnsctvspT_rec = new TIDA::Histogram<float>( monTool(), "nsct_vs_pT_rec" );
182 
183 
184  m_hnsihits_lb = new TIDA::Histogram<float>( monTool(), "nsihits_lb" );
185  m_hnsihits_lb_rec = new TIDA::Histogram<float>( monTool(), "nsihits_lb_rec" );
186 
187 
188  m_hlayer_rec = new TIDA::Histogram<float>( monTool(), "layer_rec" );
189  m_hlayer = new TIDA::Histogram<float>( monTool(), "layer" );
190 
192 
193 
194  m_hpTres = new TIDA::Histogram<float>( monTool(), "Res_pT" );
195  m_hipTres = new TIDA::Histogram<float>( monTool(), "Res_ipT" );
196  m_hetares = new TIDA::Histogram<float>( monTool(), "Res_eta" );
197  m_hphires = new TIDA::Histogram<float>( monTool(), "Res_phi" );
198  m_hd0res = new TIDA::Histogram<float>( monTool(), "Res_d0" );
199  m_hz0res = new TIDA::Histogram<float>( monTool(), "Res_z0" );
200 
201 
203 
204  m_htrkpT_residual = new TIDA::Histogram<float>( monTool(), "residual_pT" );
205  m_htrkipT_residual = new TIDA::Histogram<float>( monTool(), "residual_ipT" );
206  m_htrkphi_residual = new TIDA::Histogram<float>( monTool(), "residual_phi" );
207  m_htrketa_residual = new TIDA::Histogram<float>( monTool(), "residual_eta" );
208  m_htrkd0_residual = new TIDA::Histogram<float>( monTool(), "residual_d0" );
209  m_htrkz0_residual = new TIDA::Histogram<float>( monTool(), "residual_z0" );
210 
211  m_htrkdd0_residual = new TIDA::Histogram<float>( monTool(), "residual_dd0" );
212  m_htrkdz0_residual = new TIDA::Histogram<float>( monTool(), "residual_dz0" );
213 
214 
215 
216  m_hnpix = new TIDA::Histogram<float>( monTool(), "npix" );
217  m_hnpix_rec = new TIDA::Histogram<float>( monTool(), "npix_rec" );
218 
219  m_hnsct = new TIDA::Histogram<float>( monTool(), "nsct" );
220  m_hnsct_rec = new TIDA::Histogram<float>( monTool(), "nsct_rec" );
221 
222  m_hnsihits = new TIDA::Histogram<float>( monTool(), "nsiHits" );
223  m_hnsihits_rec = new TIDA::Histogram<float>( monTool(), "nsiHits_rec" );
224 
225  m_hntrt = new TIDA::Histogram<float>( monTool(), "ntrt" );
226  m_hntrt_rec = new TIDA::Histogram<float>( monTool(), "ntrt_rec" );
227 
228  m_chi2dof = new TIDA::Histogram<float>( monTool(), "chi2dof" );
229  m_chi2dof_rec = new TIDA::Histogram<float>( monTool(), "chi2dof_rec" );
230 
231 
233 
234  m_hd0vsphi = new TIDA::Histogram<float>( monTool(), "d0_vs_phi_prof" );
235  m_hd0vsphi_rec = new TIDA::Histogram<float>( monTool(), "d0_vs_phi_rec_prof" );
236 
240 
241  if ( m_vtxanal ) delete m_vtxanal;
242  m_vtxanal = 0;
243 
244 
245 #if 1
246 
248  // if ( name().find("vtx")!=std::string::npos || name().find("Vtx")!=std::string::npos ||
249  // name().find("vx")!=std::string::npos || name().find("Vx")!=std::string::npos ) {
250 
251  if ( ChainString(name()).vtx()!="" ) {
252 
253  m_vtxanal = new VtxR3Analysis("VTX");
254 
256 
258  store().insert( m_vtxanal, "VTX" );
259 
262 
263  }
264 
265 #endif
266 
267 }
268 
269 extern TIDARoiDescriptor* agroi;
270 
271 
272 void AnalysisR3_Tier0::execute(const std::vector<TIDA::Track*>& referenceTracks,
273  const std::vector<TIDA::Track*>& testTracks,
274  TrackAssociator* associator,
275  const TIDA::Event* tevt,
276  double* beamline,
277  TIDARoiDescriptor* roi)
278 {
279 
281  std::vector<TIDA::Track*>::const_iterator reference = referenceTracks.begin();
282  std::vector<TIDA::Track*>::const_iterator referenceEnd = referenceTracks.end();
283 
286  m_hchain->Fill( 0.5, 1 );
287 
288  if ( roi!=nullptr ) m_hroieta->Fill( roi->eta(), 1 );
289 
290  // if ( m_eventid != event()->event_number() ) {
291  if ( m_eventid != tevt->event_number() ) {
294  // m_eventid = event()->event_number();
295  m_eventid = tevt->event_number();
296  m_hchain->Fill( 1.5, 1 );
297  }
298 
299  m_hntrk->Fill( referenceTracks.size() );
300  m_hntrk_rec->Fill( testTracks.size() );
301 
303  m_hchain->Fill(4.5, testTracks.size() );
304 
305  for( ; reference!=referenceEnd ; ++reference ) {
306 
307  // Get reference parameters
308  double referenceEta = (*reference)->eta();
309  double referencePhi = phi((*reference)->phi());
310  double referenceZ0 = (*reference)->z0();
311  double referenceD0 = (*reference)->a0();
312  double referencePT = (*reference)->pT();
313 
314  double referenceDZ0 = (*reference)->dz0();
315  double referenceDD0 = (*reference)->da0();
316 
317  // Find matched tracks
318  const TIDA::Track* test = associator->matched(*reference);
319 
320  float eff_weight = 0;
321  if (test) eff_weight = 1;
322 
323  m_htotal_efficiency->Fill(0.5, eff_weight );
324 
325  m_hpTeff->Fill( std::fabs(referencePT)*0.001, eff_weight );
326  m_hz0eff->Fill( referenceZ0, eff_weight );
327  m_hd0eff->Fill( referenceD0, eff_weight );
328  m_hetaeff->Fill( referenceEta, eff_weight );
329  m_hphieff->Fill( referencePhi, eff_weight );
331  m_hnVtxeff->Fill( beamline[3], eff_weight );
332 
333  m_hlbeff->Fill( tevt->lumi_block(), eff_weight );
334 
335  m_htrkpT->Fill( std::fabs(referencePT)*0.001 );
336  m_htrketa->Fill( referenceEta );
337  m_htrkphi->Fill( referencePhi );
338  m_htrkd0->Fill( referenceD0 );
339  m_htrkz0->Fill( referenceZ0 );
340 
341  m_htrkdd0->Fill( referenceDD0 );
342  m_htrkdz0->Fill( referenceDZ0 );
343 
344  if ( referenceDD0!=0 ) m_htrkd0sig->Fill( referenceD0/referenceDD0 );
345 
346  m_hnpixvseta->Fill( referenceEta, int(((*reference)->pixelHits()+0.5)*0.5) );
347  m_hnsctvseta->Fill( referenceEta, (*reference)->sctHits() );
348  m_hntrtvseta->Fill( referenceEta, (*reference)->strawHits() );
349 
350  if ( (*reference)->dof()!=0 ) m_chi2dof->Fill( (*reference)->chi2()/(*reference)->dof() );
351 
352  m_hnpixvsphi->Fill( referencePhi, int(((*reference)->pixelHits()+0.5)*0.5) );
353  m_hnsctvsphi->Fill( referencePhi, (*reference)->sctHits() );
354  m_hntrtvsphi->Fill( referencePhi, (*reference)->strawHits() );
355 
356  m_hnpixvsd0->Fill( referenceD0, int(((*reference)->pixelHits()+0.5)*0.5) );
357  m_hnsctvsd0->Fill( referenceD0, (*reference)->sctHits() );
358 
359  m_hnpixvspT->Fill( std::fabs(referencePT)*0.001, int(((*reference)->pixelHits()+0.5)*0.5) );
360  m_hnsctvspT->Fill( std::fabs(referencePT)*0.001, (*reference)->sctHits() );
361 
362 
363  m_hnpix->Fill( int(((*reference)->pixelHits()+0.5)*0.5) );
364  m_hnsct->Fill( (*reference)->sctHits() );
365  m_hnsihits->Fill( (*reference)->siHits() );
366  m_hntrt->Fill( (*reference)->strawHits() );
367 
368 
369  // m_hnsihits_lb->Fill( event()->lumi_block(), (*reference)->siHits() );
370  m_hnsihits_lb->Fill( tevt->lumi_block(), (*reference)->siHits() );
371 
372  m_hd0vsphi->Fill(referencePhi, referenceD0 );
373 
375  m_hchain->Fill(2.5, 1);
376 
377 
378  for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
379  if ( (*reference)->hitPattern()&(1U<<ilayer) ) m_hlayer->Fill( ilayer );
380  }
381 
382  if (test) {
383 
384  m_hchain->Fill(3.5, 1);
385 
388 
390  // m_hpTres->Fill( referencePT*0.001, (test->pT() - referencePT)*0.001 );
391  // m_hipTres->Fill( 1000/referencePT, (1000/test->pT() - 1000/referencePT) );
392  // m_hetares->Fill( referenceEta, test->eta() - referenceEta );
393  // m_hphires->Fill( referencePhi, phi(test->phi() - referencePhi) );
394  // m_hd0res->Fill( referenceD0, test->a0() - referenceD0 );
395  // m_hz0res->Fill( referenceZ0, test->z0() - referenceZ0 );
396 
398  m_hpTres->Fill( referenceEta, (test->pT() - referencePT)*0.001 );
399  m_hipTres->Fill( referenceEta, (1000/test->pT() - 1000/referencePT) );
400  m_hetares->Fill( referenceEta, test->eta() - referenceEta );
401  m_hphires->Fill( referenceEta, phi(test->phi() - referencePhi) );
402  m_hd0res->Fill( referenceEta, test->a0() - referenceD0 );
403  m_hz0res->Fill( referenceEta, test->z0() - referenceZ0 );
404 
405  // m_htrkvtx_x_lb->Fill( event()->lumi_block(), beamTestx() );
406  // m_htrkvtx_y_lb->Fill( event()->lumi_block(), beamTesty() );
407  // m_htrkvtx_z_lb->Fill( event()->lumi_block(), beamTestz() );
408 
409  m_htrkvtx_x_lb->Fill( tevt->lumi_block(), beamline[0] );
410  m_htrkvtx_y_lb->Fill( tevt->lumi_block(), beamline[1] );
411  m_htrkvtx_z_lb->Fill( tevt->lumi_block(), beamline[2] );
412 
413  for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
414  if ( test->hitPattern()&(1U<<ilayer) ) m_hlayer_rec->Fill( ilayer );
415  }
416 
417  // std::cout << "SUTT beam x " << beamTestx() << " " << "\tx " << beamTesty() << " " << "\ty " << beamTestz() << std::endl;
418 
419 #if 0
420  m_htrkpT_rec->Fill( referencePT*0.001 );
422  m_htrketa_rec->Fill( referenceEta );
423  m_htrkphi_rec->Fill( referencePhi );
424  m_htrkd0_rec->Fill( referenceD0 );
425  m_htrkz0_rec->Fill( referenceZ0 );
426 
427 #endif
428 
430  m_htrkpT_rec->Fill( std::fabs(test->pT())*0.001 );
431  m_htrketa_rec->Fill( test->eta() );
432  m_htrkphi_rec->Fill( test->phi() );
433  m_htrkd0_rec->Fill( test->a0() );
434  m_htrkz0_rec->Fill( test->z0() );
435 
436  m_htrkdd0_rec->Fill( test->da0() );
437  m_htrkdz0_rec->Fill( test->dz0() );
438 
439  if ( test->da0()!=0 ) m_htrkd0sig_rec->Fill( test->a0()/test->da0() );
440 
441 
443  m_htrkpT_residual->Fill( (test->pT() - referencePT)*0.001 );
444  m_htrkipT_residual->Fill( (1000/test->pT() - 1000/referencePT) );
445  m_htrketa_residual->Fill( test->eta() - referenceEta );
446  m_htrkphi_residual->Fill( phi(test->phi() - referencePhi) );
447 
448  m_htrkd0_residual->Fill( test->a0() - referenceD0 );
449  m_htrkz0_residual->Fill( test->z0() - referenceZ0 );
450 
451  m_htrkdd0_residual->Fill( test->da0() - referenceDD0 );
452  m_htrkdz0_residual->Fill( test->dz0() - referenceDZ0 );
453 
454  m_hnpixvseta_rec->Fill( referenceEta, int((test->pixelHits()+0.5)*0.5) );
455  m_hnsctvseta_rec->Fill( referenceEta, test->sctHits() );
456 
457  m_hnpixvsphi_rec->Fill( referencePhi, int((test->pixelHits()+0.5)*0.5) );
458  m_hnsctvsphi_rec->Fill( referencePhi, test->sctHits() );
459 
460  m_hnpixvsd0_rec->Fill( referenceD0, int((test->pixelHits()+0.5)*0.5) );
461  m_hnsctvsd0_rec->Fill( referenceD0, test->sctHits() );
462 
463  m_hnpixvspT_rec->Fill( std::fabs(referencePT)*0.001, int((test->pixelHits()+0.5)*0.5) );
464  m_hnsctvspT_rec->Fill( std::fabs(referencePT)*0.001, test->sctHits() );
465 
466  m_hnpix_rec->Fill( int((test->pixelHits()+0.5)*0.5) );
467  m_hnsct_rec->Fill( test->sctHits() );
468  m_hnsihits_rec->Fill( test->siHits() );
469 
470  // m_hnsihits_lb_rec->Fill( event()->lumi_block(), test->siHits() );
471  m_hnsihits_lb_rec->Fill( tevt->lumi_block(), test->siHits() );
472 
473 
474  m_hntrt_rec->Fill( test->strawHits() );
475 
476  m_hntrtvseta_rec->Fill( referenceEta, test->strawHits() );
477  m_hntrtvsphi_rec->Fill( referencePhi, test->strawHits() );
478 
479  m_hd0vsphi_rec->Fill( test->phi(), test->a0() );
480 
481  if ( test->dof()!=0 ) m_chi2dof_rec->Fill( test->chi2()/test->dof() );
482 
483  }
484 
485  }
486 }
487 
488 
489 
490 
491 void AnalysisR3_Tier0::execute_vtx(const std::vector<TIDA::Vertex*>& vtx0,
492  const std::vector<TIDA::Vertex*>& vtx1,
493  const TIDA::Event* tevt ) {
494  if ( m_vtxanal ) m_vtxanal->execute( vtx0, vtx1, tevt );
495 }
496 
497 
499  if ( m_vtxanal ) {
500  m_vtxanal->finalise();
501  delete m_vtxanal;
502  m_vtxanal = 0;
503  }
504 }
505 
506 
507 
508 double AnalysisR3_Tier0::phi(double p) {
509  if(p < -M_PI) p += 2*M_PI;
510  if(p > M_PI) p -= 2*M_PI;
511  return p;
512 }
513 
514 
515 
516 
517 
AnalysisR3_Tier0::initialise_R3
virtual void initialise_R3()
Definition: AnalysisR3_Tier0.cxx:41
TIDA::Associator
Definition: TIDAAssociator.h:24
Track.h
AnalysisR3_Tier0::m_hntrtvseta
TIDA::Histogram< float > * m_hntrtvseta
Definition: AnalysisR3_Tier0.h:142
AnalysisR3_Tier0::m_hnpix
TIDA::Histogram< float > * m_hnpix
Definition: AnalysisR3_Tier0.h:182
AnalysisR3_Tier0::m_hnpixvseta_rec
TIDA::Histogram< float > * m_hnpixvseta_rec
Definition: AnalysisR3_Tier0.h:147
AnalysisR3_Tier0::m_htrkdd0_residual
TIDA::Histogram< float > * m_htrkdd0_residual
Definition: AnalysisR3_Tier0.h:223
ChainString::tail
std::string tail() const
Definition: ChainString.h:34
AnalysisR3_Tier0::m_hd0vsphi
TIDA::Histogram< float > * m_hd0vsphi
Definition: AnalysisR3_Tier0.h:162
AnalysisR3_Tier0::m_htrkz0_rec
TIDA::Histogram< float > * m_htrkz0_rec
Definition: AnalysisR3_Tier0.h:200
AnalysisR3_Tier0::m_hnsihits_lb_rec
TIDA::Histogram< float > * m_hnsihits_lb_rec
Definition: AnalysisR3_Tier0.h:160
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
AnalysisR3_Tier0::m_hz0eff
TIDA::Histogram< float > * m_hz0eff
Definition: AnalysisR3_Tier0.h:116
VtxR3Analysis
Definition: VtxR3Analysis.h:24
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
AnalysisR3_Tier0::m_htrketa_residual
TIDA::Histogram< float > * m_htrketa_residual
Definition: AnalysisR3_Tier0.h:218
AnalysisR3_Tier0::phi
double phi(double p)
Definition: AnalysisR3_Tier0.cxx:508
AnalysisR3_Tier0::m_hpTres
TIDA::Histogram< float > * m_hpTres
Definition: AnalysisR3_Tier0.h:120
AnalysisR3_Tier0::m_htrkphi_residual
TIDA::Histogram< float > * m_htrkphi_residual
Definition: AnalysisR3_Tier0.h:219
AnalysisR3_Tier0::m_htrkdz0_rec
TIDA::Histogram< float > * m_htrkdz0_rec
Definition: AnalysisR3_Tier0.h:203
AnalysisR3_Tier0::m_htrkd0sig_rec
TIDA::Histogram< float > * m_htrkd0sig_rec
Definition: AnalysisR3_Tier0.h:205
AnalysisR3_Tier0::m_htrkdz0_residual
TIDA::Histogram< float > * m_htrkdz0_residual
Definition: AnalysisR3_Tier0.h:224
AnalysisR3_Tier0::m_hnsctvseta_rec
TIDA::Histogram< float > * m_hnsctvseta_rec
Definition: AnalysisR3_Tier0.h:152
AnalysisR3_Tier0::m_hnpixvsd0
TIDA::Histogram< float > * m_hnpixvsd0
Definition: AnalysisR3_Tier0.h:134
AnalysisR3_Tier0::m_hnsihits
TIDA::Histogram< float > * m_hnsihits
Definition: AnalysisR3_Tier0.h:184
AnalysisR3_Tier0::m_hd0vsphi_rec
TIDA::Histogram< float > * m_hd0vsphi_rec
Definition: AnalysisR3_Tier0.h:163
AnalysisR3_Tier0::m_chi2dof_rec
TIDA::Histogram< float > * m_chi2dof_rec
Definition: AnalysisR3_Tier0.h:190
AnalysisR3_Tier0::m_hd0res
TIDA::Histogram< float > * m_hd0res
Definition: AnalysisR3_Tier0.h:124
VtxR3Analysis::set_monTool
void set_monTool(ToolHandle< GenericMonitoringTool > *m)
Definition: VtxR3Analysis.h:40
ChainString::head
std::string head() const
Definition: ChainString.h:33
AnalysisR3_Tier0::m_hntrt_rec
TIDA::Histogram< float > * m_hntrt_rec
Definition: AnalysisR3_Tier0.h:210
AnalysisR3_Tier0::m_hroieta
TIDA::Histogram< float > * m_hroieta
Definition: AnalysisR3_Tier0.h:167
M_PI
#define M_PI
Definition: ActiveFraction.h:11
AnalysisR3_Tier0::m_hnpixvspT
TIDA::Histogram< float > * m_hnpixvspT
Definition: AnalysisR3_Tier0.h:135
AnalysisR3_Tier0::m_htrkdd0_rec
TIDA::Histogram< float > * m_htrkdd0_rec
Definition: AnalysisR3_Tier0.h:202
TIDARoiDescriptor
Describes the Region of Ineterest geometry It has basically 8 parameters.
Definition: TIDARoiDescriptor.h:42
ChainString
Definition: ChainString.h:23
AnalysisR3_Tier0::m_hetaeff
TIDA::Histogram< float > * m_hetaeff
Definition: AnalysisR3_Tier0.h:113
reference
Definition: hcg.cxx:437
TrigInDetValidation_Base.test
test
Definition: TrigInDetValidation_Base.py:146
TrackAnalysis::roi
const TIDARoiDescriptor * roi() const
Definition: TrackAnalysis.h:134
AnalysisR3_Tier0::m_htrkz0
TIDA::Histogram< float > * m_htrkz0
Definition: AnalysisR3_Tier0.h:174
TIDA::Event
Definition: TIDAEvent.h:33
AnalysisR3_Tier0::m_htrkvtx_x_lb
TIDA::Histogram< float > * m_htrkvtx_x_lb
Definition: AnalysisR3_Tier0.h:127
TIDA::Event::lumi_block
void lumi_block(unsigned lb)
Definition: TIDAEvent.h:44
ChainString::roi
std::string roi() const
Definition: ChainString.h:35
AnalysisR3_Tier0::m_hnpix_rec
TIDA::Histogram< float > * m_hnpix_rec
Definition: AnalysisR3_Tier0.h:207
TIDA::Associator::matched
virtual const S * matched(T *t)
Definition: TIDAAssociator.h:45
AnalysisR3_Tier0::m_hchain
TIDA::Histogram< float > * m_hchain
Definition: AnalysisR3_Tier0.h:166
AnalysisR3_Tier0::m_htrkdz0
TIDA::Histogram< float > * m_htrkdz0
Definition: AnalysisR3_Tier0.h:178
TIDA::Histogram< float >
AnalysisR3_Tier0::m_hntrtvsphi_rec
TIDA::Histogram< float > * m_hntrtvsphi_rec
Definition: AnalysisR3_Tier0.h:158
AnalysisR3_Tier0::m_hz0res
TIDA::Histogram< float > * m_hz0res
Definition: AnalysisR3_Tier0.h:125
AnalysisR3_Tier0::finalise
virtual void finalise()
Definition: AnalysisR3_Tier0.cxx:498
TIDA::Histogram::Fill
void Fill(T d) const
Definition: TIDAHistogram.h:35
VtxR3Analysis::initialise
void initialise()
Definition: VtxR3Analysis.cxx:19
AnalysisR3_Tier0::m_htrkvtx_y_lb
TIDA::Histogram< float > * m_htrkvtx_y_lb
Definition: AnalysisR3_Tier0.h:128
AnalysisR3_Tier0::m_hipTres
TIDA::Histogram< float > * m_hipTres
Definition: AnalysisR3_Tier0.h:121
AnalysisR3_Tier0::m_hnpixvsphi_rec
TIDA::Histogram< float > * m_hnpixvsphi_rec
Definition: AnalysisR3_Tier0.h:148
AnalysisR3_Tier0::m_hlbeff
TIDA::Histogram< float > * m_hlbeff
Definition: AnalysisR3_Tier0.h:118
VtxR3Analysis::execute
void execute(const std::vector< TIDA::Vertex * > &vtx0, const std::vector< TIDA::Vertex * > &vtx1, const TIDA::Event *te=0)
Definition: VtxR3Analysis.cxx:63
AnalysisR3_Tier0::AnalysisR3_Tier0
AnalysisR3_Tier0(const std::string &name, double pTCut, double etaCut, double d0Cut, double z0Cut)
Definition: AnalysisR3_Tier0.cxx:19
AnalysisR3_Tier0::m_hphires
TIDA::Histogram< float > * m_hphires
Definition: AnalysisR3_Tier0.h:123
AnalysisR3_Tier0::monTool
ToolHandle< GenericMonitoringTool > * monTool()
Definition: AnalysisR3_Tier0.h:90
TIDA::Event::event_number
void event_number(unsigned long long e)
Definition: TIDAEvent.h:43
AnalysisR3_Tier0::m_hnsctvseta
TIDA::Histogram< float > * m_hnsctvseta
Definition: AnalysisR3_Tier0.h:137
AnalysisR3_Tier0::m_hnVtxeff
TIDA::Histogram< float > * m_hnVtxeff
Definition: AnalysisR3_Tier0.h:117
AnalysisR3_Tier0::m_hnsct
TIDA::Histogram< float > * m_hnsct
Definition: AnalysisR3_Tier0.h:183
AnalysisR3_Tier0::m_hnpixvseta
TIDA::Histogram< float > * m_hnpixvseta
Definition: AnalysisR3_Tier0.h:132
AnalysisR3_Tier0::m_hntrtvseta_rec
TIDA::Histogram< float > * m_hntrtvseta_rec
Definition: AnalysisR3_Tier0.h:157
TrackAnalysis
Definition: TrackAnalysis.h:32
TrackAnalysis::name
std::string name() const
return identifier
Definition: TrackAnalysis.h:48
AnalysisR3_Tier0::m_hpTeff
TIDA::Histogram< float > * m_hpTeff
Definition: AnalysisR3_Tier0.h:112
AnalysisR3_Tier0::m_htrkphi_rec
TIDA::Histogram< float > * m_htrkphi_rec
Definition: AnalysisR3_Tier0.h:198
AnalysisR3_Tier0::m_htrkvtx_z_lb
TIDA::Histogram< float > * m_htrkvtx_z_lb
Definition: AnalysisR3_Tier0.h:129
AnalysisR3_Tier0::m_hnsctvsphi
TIDA::Histogram< float > * m_hnsctvsphi
Definition: AnalysisR3_Tier0.h:138
TIDA::FeatureStore::insert
void insert(T *t, const std::string &key)
Definition: TIDAFeatureStore.h:44
AnalysisR3_Tier0::m_htrkd0
TIDA::Histogram< float > * m_htrkd0
Definition: AnalysisR3_Tier0.h:173
AnalysisR3_Tier0::m_htrkpT_rec
TIDA::Histogram< float > * m_htrkpT_rec
Definition: AnalysisR3_Tier0.h:196
AnalysisR3_Tier0::m_hntrtvsphi
TIDA::Histogram< float > * m_hntrtvsphi
Definition: AnalysisR3_Tier0.h:143
AnalysisR3_Tier0::m_vtxanal
VtxR3Analysis * m_vtxanal
Definition: AnalysisR3_Tier0.h:236
AnalysisR3_Tier0::m_hnpixvsd0_rec
TIDA::Histogram< float > * m_hnpixvsd0_rec
Definition: AnalysisR3_Tier0.h:149
VtxR3Analysis::finalise
void finalise()
Definition: VtxR3Analysis.cxx:151
ChainString::raw
std::string raw() const
Definition: ChainString.h:53
ChainString::element
std::string element() const
Definition: ChainString.h:37
AnalysisR3_Tier0::m_hd0eff
TIDA::Histogram< float > * m_hd0eff
Definition: AnalysisR3_Tier0.h:115
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:195
AnalysisR3_Tier0::m_htrkd0_residual
TIDA::Histogram< float > * m_htrkd0_residual
Definition: AnalysisR3_Tier0.h:220
AnalysisR3_Tier0::m_hnpixvsphi
TIDA::Histogram< float > * m_hnpixvsphi
Definition: AnalysisR3_Tier0.h:133
AnalysisR3_Tier0::m_hnsctvspT
TIDA::Histogram< float > * m_hnsctvspT
Definition: AnalysisR3_Tier0.h:140
ChainString.h
AnalysisR3_Tier0::m_htrkd0_rec
TIDA::Histogram< float > * m_htrkd0_rec
Definition: AnalysisR3_Tier0.h:199
AnalysisR3_Tier0::m_hnsctvsd0
TIDA::Histogram< float > * m_hnsctvsd0
Definition: AnalysisR3_Tier0.h:139
AnalysisR3_Tier0::m_hnsct_rec
TIDA::Histogram< float > * m_hnsct_rec
Definition: AnalysisR3_Tier0.h:208
TIDARoiDescriptor::eta
double eta() const
Definition: TIDARoiDescriptor.h:134
AnalysisR3_Tier0::m_htrkipT_residual
TIDA::Histogram< float > * m_htrkipT_residual
Definition: AnalysisR3_Tier0.h:217
AnalysisR3_Tier0::m_hnsihits_rec
TIDA::Histogram< float > * m_hnsihits_rec
Definition: AnalysisR3_Tier0.h:209
AnalysisR3_Tier0::m_htrkz0_residual
TIDA::Histogram< float > * m_htrkz0_residual
Definition: AnalysisR3_Tier0.h:221
AnalysisR3_Tier0::m_hnsctvsphi_rec
TIDA::Histogram< float > * m_hnsctvsphi_rec
Definition: AnalysisR3_Tier0.h:153
AnalysisR3_Tier0::execute_vtx
virtual void execute_vtx(const std::vector< TIDA::Vertex * > &vtx0, const std::vector< TIDA::Vertex * > &vtx1, const TIDA::Event *tevt=0)
Definition: AnalysisR3_Tier0.cxx:491
VtxR3Analysis.h
AnalysisR3_Tier0.h
AnalysisR3_Tier0::m_htrkphi
TIDA::Histogram< float > * m_htrkphi
Definition: AnalysisR3_Tier0.h:172
AnalysisR3_Tier0::m_monTool
ToolHandle< GenericMonitoringTool > * m_monTool
Definition: AnalysisR3_Tier0.h:238
ChainString::vtx
std::string vtx() const
Definition: ChainString.h:36
AnalysisR3_Tier0::m_htrkd0sig
TIDA::Histogram< float > * m_htrkd0sig
Definition: AnalysisR3_Tier0.h:180
agroi
TIDARoiDescriptor * agroi
AnalysisR3_Tier0::m_hntrk
TIDA::Histogram< float > * m_hntrk
Definition: AnalysisR3_Tier0.h:168
AnalysisR3_Tier0::m_hnsihits_lb
TIDA::Histogram< float > * m_hnsihits_lb
Definition: AnalysisR3_Tier0.h:145
dq_make_web_display.reference
reference
Definition: dq_make_web_display.py:44
AnalysisR3_Tier0::m_hntrk_rec
TIDA::Histogram< float > * m_hntrk_rec
Definition: AnalysisR3_Tier0.h:187
AnalysisR3_Tier0::m_htrkpT
TIDA::Histogram< float > * m_htrkpT
Definition: AnalysisR3_Tier0.h:169
TIDA::Track
Definition: Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/Track.h:26
AnalysisR3_Tier0::initialise
virtual void initialise()
standard operation interface
Definition: AnalysisR3_Tier0.cxx:35
AnalysisR3_Tier0::execute
virtual void execute(const std::vector< TIDA::Track * > &referenceTracks, const std::vector< TIDA::Track * > &testTracks, TrackAssociator *associator)
Definition: AnalysisR3_Tier0.h:45
AnalysisR3_Tier0::m_htrketa
TIDA::Histogram< float > * m_htrketa
Definition: AnalysisR3_Tier0.h:171
AnalysisR3_Tier0::m_hlayer
TIDA::Histogram< float > * m_hlayer
Definition: AnalysisR3_Tier0.h:193
AnalysisR3_Tier0::m_hphieff
TIDA::Histogram< float > * m_hphieff
Definition: AnalysisR3_Tier0.h:114
AnalysisR3_Tier0::m_eventid
unsigned long long m_eventid
Definition: AnalysisR3_Tier0.h:234
AnalysisR3_Tier0::m_hnpixvspT_rec
TIDA::Histogram< float > * m_hnpixvspT_rec
Definition: AnalysisR3_Tier0.h:150
AnalysisR3_Tier0::m_hlayer_rec
TIDA::Histogram< float > * m_hlayer_rec
Definition: AnalysisR3_Tier0.h:213
AnalysisR3_Tier0::m_hetares
TIDA::Histogram< float > * m_hetares
Definition: AnalysisR3_Tier0.h:122
AnalysisR3_Tier0::m_htrkdd0
TIDA::Histogram< float > * m_htrkdd0
Definition: AnalysisR3_Tier0.h:177
AnalysisR3_Tier0::initialise_R2
virtual void initialise_R2()
Definition: AnalysisR3_Tier0.cxx:39
TrackAnalysis::store
TIDA::FeatureStore & store()
Definition: TrackAnalysis.h:129
AnalysisR3_Tier0::m_hnsctvspT_rec
TIDA::Histogram< float > * m_hnsctvspT_rec
Definition: AnalysisR3_Tier0.h:155
AnalysisR3_Tier0::m_debug
bool m_debug
Definition: AnalysisR3_Tier0.h:232
AnalysisR3_Tier0::m_chi2dof
TIDA::Histogram< float > * m_chi2dof
Definition: AnalysisR3_Tier0.h:189
AnalysisR3_Tier0::m_htrkpT_residual
TIDA::Histogram< float > * m_htrkpT_residual
Definition: AnalysisR3_Tier0.h:216
AnalysisR3_Tier0::m_htotal_efficiency
TIDA::Histogram< float > * m_htotal_efficiency
Monitorwd::AScalar Histogram wrapper class.
Definition: AnalysisR3_Tier0.h:111
AnalysisR3_Tier0::m_htrketa_rec
TIDA::Histogram< float > * m_htrketa_rec
Definition: AnalysisR3_Tier0.h:197
AnalysisR3_Tier0::m_hntrt
TIDA::Histogram< float > * m_hntrt
Definition: AnalysisR3_Tier0.h:185
AnalysisR3_Tier0::m_hnsctvsd0_rec
TIDA::Histogram< float > * m_hnsctvsd0_rec
Definition: AnalysisR3_Tier0.h:154