226 {
227
228#if 0
229
231 std::vector<TIDA::Track*>::const_iterator
reference = referenceTracks.begin();
232 std::vector<TIDA::Track*>::const_iterator referenceEnd = referenceTracks.end();
233
237
238
239 if ( roi!=
nullptr )
m_hroieta->Fill( roi->eta(), 1 );
240
241 if ( tevt!=nullptr && m_eventid != tevt->event_number() ) {
244
245 m_eventid = tevt->event_number();
247
248 m_hmu->Fill( tevt->mu() );
249
250 }
251
252 m_hntrk->Fill( referenceTracks.size() );
254
256 m_hchain->Fill(4.5, testTracks.size() );
257
259
260
261 double referenceEta = (*reference)->eta();
262 double referencePhi =
phi((*reference)->phi());
263 double referenceZ0 = (*reference)->z0();
264 double referenceD0 = (*reference)->a0();
265 double referencePT = (*reference)->pT();
266
267 double referenceDZ0 = (*reference)->dz0();
268 double referenceDD0 = (*reference)->da0();
269
270
271 const TIDA::Track*
test = associator->matched(*reference);
272
273 float eff_weight = 0;
274 if (test) eff_weight = 1;
275
277
278 m_hpTeff->Fill( std::fabs(referencePT)*0.001, eff_weight );
279 m_hz0eff->Fill( referenceZ0, eff_weight );
280 m_hd0eff->Fill( referenceD0, eff_weight );
281 m_hetaeff->Fill( referenceEta, eff_weight );
282 m_hphieff->Fill( referencePhi, eff_weight );
284 if (beamline)
m_hnVtxeff->Fill( beamline[3], eff_weight );
285
286 if (tevt)
m_hlbeff->Fill( tevt->lumi_block(), eff_weight );
287
288 m_htrkpT->Fill( std::fabs(referencePT)*0.001 );
293
296
297 if ( referenceDD0!=0 )
m_htrkd0sig->Fill( referenceD0/referenceDD0 );
298
299 m_hnpixvseta->Fill( referenceEta,
int(((*reference)->pixelHits()+0.5)*0.5) );
300 m_hnsctvseta->Fill( referenceEta, (*reference)->sctHits() );
301 m_hntrtvseta->Fill( referenceEta, (*reference)->strawHits() );
302
303 if ( (*reference)->dof()!=0 )
m_chi2dof->Fill( (*reference)->chi2()/(*reference)->dof() );
304
305 m_hnpixvsphi->Fill( referencePhi,
int(((*reference)->pixelHits()+0.5)*0.5) );
306 m_hnsctvsphi->Fill( referencePhi, (*reference)->sctHits() );
307 m_hntrtvsphi->Fill( referencePhi, (*reference)->strawHits() );
308
309 m_hnpixvsd0->Fill( referenceD0,
int(((*reference)->pixelHits()+0.5)*0.5) );
310 m_hnsctvsd0->Fill( referenceD0, (*reference)->sctHits() );
311
312 m_hnpixvspT->Fill( std::fabs(referencePT)*0.001,
int(((*reference)->pixelHits()+0.5)*0.5) );
313 m_hnsctvspT->Fill( std::fabs(referencePT)*0.001, (*reference)->sctHits() );
314
315
316 m_hnpix->Fill(
int(((*reference)->pixelHits()+0.5)*0.5) );
317 m_hnsct->Fill( (*reference)->sctHits() );
319 m_hntrt->Fill( (*reference)->strawHits() );
320
321
322
323 if (tevt)
m_hnsihits_lb->Fill( tevt->lumi_block(), (*reference)->siHits() );
324
326
329
330
331 for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
332 if ( (*reference)->hitPattern()&(1U<<ilayer) )
m_hlayer->Fill( ilayer );
333 }
334
335 if (test) {
336
338
341
343
344
345
346
347
348
349
351 m_hpTres->Fill( referenceEta, (
test->pT() - referencePT)*0.001 );
352 m_hipTres->Fill( referenceEta, (1000/
test->pT() - 1000/referencePT) );
355 m_hd0res->Fill( referenceEta,
test->a0() - referenceD0 );
356 m_hz0res->Fill( referenceEta,
test->z0() - referenceZ0 );
357
358
359
360
361
362 if (tevt && beamline) {
366 }
367
368 for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
370 }
371
372
373
374#if 0
381
382#endif
383
390
393
395
396
402
405
408
411
414
417
420
424
426
427
429
432
434
436
437 }
438
439 }
440
441#endif
442
443}
Scalar phi() const
phi method
TIDA::Histogram< float > m_htrkpT_rec
TIDA::Histogram< float > m_hnsihits
TIDA::Histogram< float > m_hnVtxeff
TIDA::Histogram< float > m_hntrtvsphi
TIDA::Histogram< float > m_htotal_efficiency
Monitorwd::AScalar Histogram wrapper class.
TIDA::Histogram< float > m_hnsctvseta
TIDA::Histogram< float > m_htrkdd0_rec
TIDA::Histogram< float > m_hnsctvsphi_rec
TIDA::Histogram< float > m_htrkpT_residual
TIDA::Histogram< float > m_hlayer_rec
TIDA::Histogram< float > m_hnpixvspT
TIDA::Histogram< float > m_hnsihits_lb_rec
TIDA::Histogram< float > m_hnsctvsd0
TIDA::Histogram< float > m_hnsctvsphi
TIDA::Histogram< float > m_hnpixvsd0_rec
TIDA::Histogram< float > m_hd0vsphi_rec
TIDA::Histogram< float > m_hnpixvsd0
TIDA::Histogram< float > m_hetares
TIDA::Histogram< float > m_hntrt_rec
TIDA::Histogram< float > m_htrketa_residual
TIDA::Histogram< float > m_htrkd0
TIDA::Histogram< float > m_hnsihits_rec
TIDA::Histogram< float > m_hntrtvseta_rec
TIDA::Histogram< float > m_htrkvtx_y_lb
TIDA::Histogram< float > m_hnsct_rec
TIDA::Histogram< float > m_chi2dof_rec
TIDA::Histogram< float > m_hnsct
TIDA::Histogram< float > m_hnpix
TIDA::Histogram< float > m_hmu
TIDA::Histogram< float > m_htrketa_rec
TIDA::Histogram< float > m_hnsctvspT
TIDA::Histogram< float > m_hpTeff
TIDA::Histogram< float > m_hlayer
TIDA::Histogram< float > m_htrkphi
TIDA::Histogram< float > m_hd0eff
TIDA::Histogram< float > m_htrkd0sig
TIDA::Histogram< float > m_htrkdz0_residual
TIDA::Histogram< float > m_hnpixvsphi
TIDA::Histogram< float > m_hnsihits_lb
TIDA::Histogram< float > m_htrkdz0
TIDA::Histogram< float > m_hipTres
TIDA::Histogram< float > m_hetaeff
TIDA::Histogram< float > m_hnsctvspT_rec
TIDA::Histogram< float > m_hd0res
TIDA::Histogram< float > m_htrkvtx_z_lb
TIDA::Histogram< float > m_hchain
TIDA::Histogram< float > m_htrkphi_residual
TIDA::Histogram< float > m_hz0eff
TIDA::Histogram< float > m_htrkd0_rec
TIDA::Histogram< float > m_htrkpT
TIDA::Histogram< float > m_htrkdd0_residual
TIDA::Histogram< float > m_htrkz0_rec
TIDA::Histogram< float > m_htrkz0
TIDA::Histogram< float > m_htrkd0_residual
TIDA::Histogram< float > m_hd0vsphi
TIDA::Histogram< float > m_htrkipT_residual
TIDA::Histogram< float > m_htrkd0sig_rec
TIDA::Histogram< float > m_hnpixvseta_rec
TIDA::Histogram< float > m_hpTres
TIDA::Histogram< float > m_hnpixvspT_rec
TIDA::Histogram< float > m_hnpixvsphi_rec
TIDA::Histogram< float > m_chi2dof
TIDA::Histogram< float > m_htrkdd0
TIDA::Histogram< float > m_hntrt
TIDA::Histogram< float > m_hphires
TIDA::Histogram< float > m_hntrtvseta
TIDA::Histogram< float > m_htrkvtx_x_lb
TIDA::Histogram< float > m_hntrk
TIDA::Histogram< float > m_hz0res
TIDA::Histogram< float > m_hntrk_rec
TIDA::Histogram< float > m_hnsctvseta_rec
TIDA::Histogram< float > m_htrketa
TIDA::Histogram< float > m_hnpixvseta
TIDA::Histogram< float > m_hntrtvsphi_rec
TIDA::Histogram< float > m_hnpix_rec
TIDA::Histogram< float > m_htrkz0_residual
TIDA::Histogram< float > m_hphieff
TIDA::Histogram< float > m_htrkdz0_rec
TIDA::Histogram< float > m_hlbeff
TIDA::Histogram< float > m_htrkphi_rec
TIDA::Histogram< float > m_hnsctvsd0_rec
TIDA::Histogram< float > m_hroieta