276 const std::vector<TIDA::Track*>& testTracks,
284 std::vector<TIDA::Track*>::const_iterator
reference = referenceTracks.begin();
285 std::vector<TIDA::Track*>::const_iterator referenceEnd = referenceTracks.end();
305 m_hntrk->Fill( referenceTracks.size() );
309 m_hchain->Fill(4.5, testTracks.size() );
314 double referenceEta = (*reference)->eta();
315 double referencePhi =
phi((*reference)->phi());
316 double referenceZ0 = (*reference)->z0();
317 double referenceD0 = (*reference)->a0();
318 double referencePT = (*reference)->pT();
320 double referenceDZ0 = (*reference)->dz0();
321 double referenceDD0 = (*reference)->da0();
326 float eff_weight = 0;
327 if (test) eff_weight = 1;
331 m_hpTeff->Fill( std::fabs(referencePT)*0.001, eff_weight );
332 m_hz0eff->Fill( referenceZ0, eff_weight );
333 m_hd0eff->Fill( referenceD0, eff_weight );
334 m_hetaeff->Fill( referenceEta, eff_weight );
335 m_hphieff->Fill( referencePhi, eff_weight );
337 if (beamline)
m_hnVtxeff->Fill( beamline[3], eff_weight );
341 m_hntraxeff->Fill( referenceTracks.size(), eff_weight );
343 m_htrkpT->Fill( std::fabs(referencePT)*0.001 );
352 if ( referenceDD0!=0 )
m_htrkd0sig->Fill( referenceD0/referenceDD0 );
354 m_hnpixvseta->Fill( referenceEta,
int(((*reference)->pixelHits()+0.5)*0.5) );
355 m_hnsctvseta->Fill( referenceEta, (*reference)->sctHits() );
356 m_hntrtvseta->Fill( referenceEta, (*reference)->strawHits() );
358 if ( (*reference)->dof()!=0 )
m_chi2dof->Fill( (*reference)->chi2()/(*reference)->dof() );
360 m_hnpixvsphi->Fill( referencePhi,
int(((*reference)->pixelHits()+0.5)*0.5) );
361 m_hnsctvsphi->Fill( referencePhi, (*reference)->sctHits() );
362 m_hntrtvsphi->Fill( referencePhi, (*reference)->strawHits() );
364 m_hnpixvsd0->Fill( referenceD0,
int(((*reference)->pixelHits()+0.5)*0.5) );
365 m_hnsctvsd0->Fill( referenceD0, (*reference)->sctHits() );
367 m_hnpixvspT->Fill( std::fabs(referencePT)*0.001,
int(((*reference)->pixelHits()+0.5)*0.5) );
368 m_hnsctvspT->Fill( std::fabs(referencePT)*0.001, (*reference)->sctHits() );
371 m_hnpix->Fill(
int(((*reference)->pixelHits()+0.5)*0.5) );
372 m_hnsct->Fill( (*reference)->sctHits() );
374 m_hntrt->Fill( (*reference)->strawHits() );
386 for (
size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
387 if ( (*reference)->hitPattern()&(1U<<ilayer) )
m_hlayer->Fill( ilayer );
406 m_hpTres->Fill( referenceEta, (test->pT() - referencePT)*0.001 );
407 m_hipTres->Fill( referenceEta, (1000/test->pT() - 1000/referencePT) );
408 m_hetares->Fill( referenceEta, test->eta() - referenceEta );
409 m_hphires->Fill( referenceEta,
phi(test->phi() - referencePhi) );
410 m_hd0res->Fill( referenceEta, test->a0() - referenceD0 );
411 m_hz0res->Fill( referenceEta, test->z0() - referenceZ0 );
417 if (tevt && beamline) {
423 for (
size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
424 if ( test->hitPattern()&(1U<<ilayer) )
m_hlayer_rec->Fill( ilayer );
470 m_hnpixvsd0_rec->Fill( referenceD0,
int((test->pixelHits()+0.5)*0.5) );
473 m_hnpixvspT_rec->Fill( std::fabs(referencePT)*0.001,
int((test->pixelHits()+0.5)*0.5) );
474 m_hnsctvspT_rec->Fill( std::fabs(referencePT)*0.001, test->sctHits() );
476 m_hnpix_rec->Fill(
int((test->pixelHits()+0.5)*0.5) );
490 if ( test->dof()!=0 )
m_chi2dof_rec->Fill( test->chi2()/test->dof() );