210 {
211
212#if 0
213
215 std::vector<TIDA::Track*>::const_iterator
reference = referenceTracks.begin();
216 std::vector<TIDA::Track*>::const_iterator referenceEnd = referenceTracks.end();
217
221
222
223 if ( roi!=
nullptr )
m_hroieta->Fill( roi->eta(), 1 );
224
225 if ( tevt!=nullptr && m_eventid != tevt->event_number() ) {
228
229 m_eventid = tevt->event_number();
231
232 m_hmu->Fill( tevt->mu() );
233
234 }
235
236 m_hntrk->Fill( referenceTracks.size() );
238
240 m_hchain->Fill(4.5, testTracks.size() );
241
243
244
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();
250
251 double referenceDZ0 = (*reference)->dz0();
252 double referenceDD0 = (*reference)->da0();
253
254
255 const TIDA::Track*
test = associator->matched(*reference);
256
257 float eff_weight = 0;
258 if (test) eff_weight = 1;
259
261
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 );
269
270 if (tevt)
m_hlbeff->Fill( tevt->lumi_block(), eff_weight );
271
272 m_htrkpT->Fill( std::fabs(referencePT)*0.001 );
277
280
281 if ( referenceDD0!=0 )
m_htrkd0sig->Fill( referenceD0/referenceDD0 );
282
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() );
286
287 if ( (*reference)->dof()!=0 )
m_chi2dof->Fill( (*reference)->chi2()/(*reference)->dof() );
288
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() );
292
293 m_hnpixvsd0->Fill( referenceD0,
int(((*reference)->pixelHits()+0.5)*0.5) );
294 m_hnsctvsd0->Fill( referenceD0, (*reference)->sctHits() );
295
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() );
298
299
300 m_hnpix->Fill(
int(((*reference)->pixelHits()+0.5)*0.5) );
301 m_hnsct->Fill( (*reference)->sctHits() );
303 m_hntrt->Fill( (*reference)->strawHits() );
304
305
306
307 if (tevt)
m_hnsihits_lb->Fill( tevt->lumi_block(), (*reference)->siHits() );
308
310
313
314
315 for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
316 if ( (*reference)->hitPattern()&(1U<<ilayer) )
m_hlayer->Fill( ilayer );
317 }
318
319 if (test) {
320
322
325
327
328
329
330
331
332
333
335 m_hpTres->Fill( referenceEta, (
test->pT() - referencePT)*0.001 );
336 m_hipTres->Fill( referenceEta, (1000/
test->pT() - 1000/referencePT) );
339 m_hd0res->Fill( referenceEta,
test->a0() - referenceD0 );
340 m_hz0res->Fill( referenceEta,
test->z0() - referenceZ0 );
341
342
343
344
345
346 if (tevt && beamline) {
350 }
351
352 for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) {
354 }
355
356
357
358#if 0
365
366#endif
367
374
377
379
380
386
389
392
395
398
401
404
408
410
411
413
416
418
420
421 }
422
423 }
424
425#endif
426
427}
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