215 std::vector<TIDA::Track*>::const_iterator
reference = referenceTracks.begin();
216 std::vector<TIDA::Track*>::const_iterator referenceEnd = referenceTracks.end();
223 if ( roi!=
nullptr )
m_hroieta->Fill( roi->eta(), 1 );
225 if ( tevt!=
nullptr && m_eventid != tevt->event_number() ) {
229 m_eventid = tevt->event_number();
232 m_hmu->Fill( tevt->mu() );
236 m_hntrk->Fill( referenceTracks.size() );
240 m_hchain->Fill(4.5, testTracks.size() );
245 double referenceEta = (*reference)->eta();
246 double referencePhi =
phi((*reference)->phi());
247 double referenceZ0 = (*reference)->z0();
248 double referenceD0 = (*reference)->a0();
249 double referencePT = (*reference)->pT();
251 double referenceDZ0 = (*reference)->dz0();
252 double referenceDD0 = (*reference)->da0();
257 float eff_weight = 0;
258 if (test) eff_weight = 1;
262 m_hpTeff->Fill( std::fabs(referencePT)*0.001, eff_weight );
263 m_hz0eff->Fill( referenceZ0, eff_weight );
264 m_hd0eff->Fill( referenceD0, eff_weight );
265 m_hetaeff->Fill( referenceEta, eff_weight );
266 m_hphieff->Fill( referencePhi, eff_weight );
268 if (beamline)
m_hnVtxeff->Fill( beamline[3], eff_weight );
270 if (tevt)
m_hlbeff->Fill( tevt->lumi_block(), eff_weight );
272 m_htrkpT->Fill( std::fabs(referencePT)*0.001 );
281 if ( referenceDD0!=0 )
m_htrkd0sig->Fill( referenceD0/referenceDD0 );
283 m_hnpixvseta->Fill( referenceEta,
int(((*reference)->pixelHits()+0.5)*0.5) );
284 m_hnsctvseta->Fill( referenceEta, (*reference)->sctHits() );
285 m_hntrtvseta->Fill( referenceEta, (*reference)->strawHits() );
287 if ( (*reference)->dof()!=0 )
m_chi2dof->Fill( (*reference)->chi2()/(*reference)->dof() );
289 m_hnpixvsphi->Fill( referencePhi,
int(((*reference)->pixelHits()+0.5)*0.5) );
290 m_hnsctvsphi->Fill( referencePhi, (*reference)->sctHits() );
291 m_hntrtvsphi->Fill( referencePhi, (*reference)->strawHits() );
293 m_hnpixvsd0->Fill( referenceD0,
int(((*reference)->pixelHits()+0.5)*0.5) );
294 m_hnsctvsd0->Fill( referenceD0, (*reference)->sctHits() );
296 m_hnpixvspT->Fill( std::fabs(referencePT)*0.001,
int(((*reference)->pixelHits()+0.5)*0.5) );
297 m_hnsctvspT->Fill( std::fabs(referencePT)*0.001, (*reference)->sctHits() );
300 m_hnpix->Fill(
int(((*reference)->pixelHits()+0.5)*0.5) );
301 m_hnsct->Fill( (*reference)->sctHits() );
303 m_hntrt->Fill( (*reference)->strawHits() );
307 if (tevt)
m_hnsihits_lb->Fill( tevt->lumi_block(), (*reference)->siHits() );
315 for (
size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
316 if ( (*reference)->hitPattern()&(1U<<ilayer) )
m_hlayer->Fill( ilayer );
335 m_hpTres->Fill( referenceEta, (test->pT() - referencePT)*0.001 );
336 m_hipTres->Fill( referenceEta, (1000/test->pT() - 1000/referencePT) );
337 m_hetares->Fill( referenceEta, test->eta() - referenceEta );
338 m_hphires->Fill( referenceEta,
phi(test->phi() - referencePhi) );
339 m_hd0res->Fill( referenceEta, test->a0() - referenceD0 );
340 m_hz0res->Fill( referenceEta, test->z0() - referenceZ0 );
346 if (tevt && beamline) {
352 for (
size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
353 if ( test->hitPattern()&(1U<<ilayer) )
m_hlayer_rec->Fill( ilayer );
399 m_hnpixvsd0_rec->Fill( referenceD0,
int((test->pixelHits()+0.5)*0.5) );
402 m_hnpixvspT_rec->Fill( std::fabs(referencePT)*0.001,
int((test->pixelHits()+0.5)*0.5) );
403 m_hnsctvspT_rec->Fill( std::fabs(referencePT)*0.001, test->sctHits() );
405 m_hnpix_rec->Fill(
int((test->pixelHits()+0.5)*0.5) );
419 if ( test->dof()!=0 )
m_chi2dof_rec->Fill( test->chi2()/test->dof() );