keep the old ranges commented, since we need the new ranges to be consistent double tmp_maxEta = 3.;
keep the old eta bins commented, since we need the new binning to be consisten with the old int etaBins = int(30*binConfig.eta_NScale);
what is goping on here? the bins for the residuals should depend on the residual itself, not the x variable, how come these are all different?
87 double tmp_absResPt = 0.5;
95 double tmp_maxEta = 5.;
96 double tmp_absResEta = 0.04;
99 double tmp_maxPhi = 3.142;
100 double tmp_absResPhi = 0.02;
114 const double zMax = binConfig.
z0Max;
116 const int zresBins = 100;
117 const double zresMax = 10;
120 const double d0Max = binConfig.
d0Max;
122 const int d0resBins = 100;
123 const double d0resMax = 5;
127 const double a0Max = binConfig.
a0Max;
129 const int a0resBins = 100;
130 const double a0resMax = 5;
155 const int ptnbins = Npt;
156 std::vector<double> ptbinlimsv(ptnbins+1);
157 double* ptbinlims = &ptbinlimsv[0];
160 for (
int i=0 ;
i<=ptnbins ;
i++ ) ptbinlims[
i] =
std::pow(10, pt_a*
i/ptnbins+pt_b)/1000;
165 const int ptnbins2 = (2*ptnbins);
167 std::vector<double> ptbinlims2v(ptnbins2 + 1);
168 double* ptbinlims2 = &ptbinlims2v[0];
170 int ptnbin_counter = 0;
171 for (
int i=ptnbins;
i>0 ;
i-- ) {
172 ptbinlims2[ptnbin_counter] =
std::pow(10, pt_a*
i/ptnbins+pt_b)/(-2000);
177 for (
int i=0 ;
i<ptnbins+1 ;
i++ ) {
178 ptbinlims2[ptnbin_counter] =
std::pow(10, pt_a*
i/ptnbins+pt_b)/2000;
184 const int iptnbins = 20;
185 const double minmaxipt=0.5;
186 std::vector<double> iptbinlimsv(iptnbins+1);
187 double* iptbinlims = &iptbinlimsv[0];
188 for (
int i=0 ;
i<=iptnbins ;
i++ ) {
189 iptbinlims[
i] = -minmaxipt+
i*minmaxipt*2./iptnbins;
205 if (
name() !=
gDirectory->GetName() ) std::cerr <<
"ConfAnalysis::initialize() Directory: problem with directory " <<
gDirectory->GetName() <<
" " <<
name() << std::endl;
275 double d0bins[40] = { -5.0, -4.0, -3.0, -2.5,
276 -2.0, -1.8, -1.6, -1.4, -1.2,
277 -1.05, -0.95, -0.85, -0.75, -0.65, -0.55, -0.45, -0.35, -0.25, -0.15, -0.05,
278 0.05, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 0.95, 1.05,
279 1.2, 1.4, 1.6, 1.8, 2.0,
280 2.5, 3.0, 4.0, 5.0 };
341 m_invmass =
new TH1F(
"invmass",
"invariant mass;mass [GeV]", 320, 0, 200 );
342 m_invmassObj =
new TH1F(
"invmassObj",
"invariant mass;mass [GeV]", 320, 0, 200 );
348 TH2F* heta_vs_pt =
new TH2F(
"eta_vs_pt",
"p_{T} [GeV],#eta", ptnbins, ptbinlims, 10, -3, 3 );
349 TH2F* hd0_vs_pt =
new TH2F(
"d0_vs_pt",
"p_{T} [GeV],d_{0} [mm]", ptnbins, ptbinlims, 30, -300, 300 );
362 m_eff_pt->
Hist()->GetYaxis()->SetTitle(
"Efficiency [%]");
374 m_eff_z0->
Hist()->GetYaxis()->SetTitle(
"Efficiency [%]");
378 m_eff_d0->
Hist()->GetYaxis()->SetTitle(
"Efficiency [%]");
382 m_eff_a0->
Hist()->GetYaxis()->SetTitle(
"Efficiency [%]");
385 m_eff_ptm->
Hist()->GetXaxis()->SetTitle(
"Negative P_{T} [GeV]");
389 m_eff_ptp->
Hist()->GetXaxis()->SetTitle(
"Positive P_{T} [GeV]");
417 addHistogram(
new TH1F(
"phi_rec",
"phi_rec", phiBins, -tmp_maxPhi, tmp_maxPhi ) );
422 addHistogram2D(
new TH2F(
"eta_phi_rec",
"eta_phi_rec", (tmp_maxEta+1)*30 , -tmp_maxEta-1, tmp_maxEta+1, (tmp_maxPhi+1)*30, -tmp_maxPhi-1, tmp_maxPhi+1 ) );
423 addHistogram2D(
new TH2F(
"phi_d0_rec",
"phi_d0_rec", (2*tmp_maxPhi+2)*15, -tmp_maxPhi-1, tmp_maxPhi+1 ,d0Bins+20, -d0Max+7, d0Max-7 ));
445 addHistogram(
new TH1F(
"eta_pull_simple",
"eta_pull_simple", 100, -10, 10) );
446 addHistogram(
new TH1F(
"phi_pull_simple",
"phi_pull_simple", 100, -10, 10) );
453 TH1F* pT_res =
new TH1F(
"pT_res",
"pT_res", 4*pTResBins, -0.1, 0.1 );
454 pT_res->GetXaxis()->SetTitle(
"#Delta P_{T} [GeV]");
455 pT_res->GetYaxis()->SetTitle(
"Entries");
459 TH1F* spT_res =
new TH1F(
"spT_res",
"spT_res", 4*pTResBins, -0.1, 0.1 );
460 spT_res->GetXaxis()->SetTitle(
"#Delta sP_{T} [GeV]");
461 spT_res->GetYaxis()->SetTitle(
"Entries");
465 TH1F* ipT_res =
new TH1F(
"ipT_res",
"ipT_res", 4*pTResBins, -0.4, 0.4 );
466 ipT_res->GetXaxis()->SetTitle(
"#Delta 1/P_{T} [GeV^{-1}]");
467 ipT_res->GetYaxis()->SetTitle(
"Entries");
470 TH1F* eta_res =
new TH1F(
"eta_res",
"eta_res", etaResBins, -2*tmp_absResEta, 2*tmp_absResEta );
471 eta_res->GetXaxis()->SetTitle(
"#Delta #eta");
472 eta_res->GetYaxis()->SetTitle(
"Entries");
481 addHistogram(
new TH1F(
"phi_res",
"phi_res;#Delta #phi;Entries", 2*phiResBins, -2*tmp_absResPhi, 2*tmp_absResPhi ) );
482 addHistogram(
new TH1F(
"z0_res",
"z0_res;#Deltaz_{0};Entries", 16*zresBins, -8*zresMax, 8*zresMax ) );
483 addHistogram(
new TH1F(
"d0_res",
"d0_res;#Deltad_{0};Entries", 4*d0resBins, -0.2*d0resMax, 0.2*d0resMax ) );
484 addHistogram(
new TH1F(
"a0_res",
"a0_res;#Deltaa_{0};Entries", 4*a0resBins, -0.2*a0resMax, 0.2*a0resMax ) );
487 addHistogram(
new TH1F(
"dphi_res",
"dphi_res;#Delta #phi;Entries", 2*phiResBins, -0.2*tmp_absResPhi, 0.2*tmp_absResPhi ) );
488 addHistogram(
new TH1F(
"dz0_res",
"dz0_res;#Deltaz_{0};Entries", 8*zresBins, -0.8*zresMax, 0.8*zresMax ) );
489 addHistogram(
new TH1F(
"dd0_res",
"dd0_res;#Deltad_{0};Entries", 4*d0resBins, -0.05*d0resMax, 0.05*d0resMax ) );
490 addHistogram(
new TH1F(
"da0_res",
"da0_res;#Deltaa_{0};Entries", 4*a0resBins, -0.05*a0resMax, 0.05*a0resMax ) );
501 m_rDd0res.push_back(
new Resplot(
"rDd0_vs_pt", ptnbins, ptbinlims, 1200, -0.1, 0.1 ) );
504 m_rDd0res.push_back(
new Resplot(
"rDd0_vs_d0", 20, -1.5, 1.5, 1200, -0.1, 0.1 ) );
509 m_rDa0res.push_back(
new Resplot(
"rDa0_vs_pt", ptnbins, ptbinlims, 100, -0.1, 0.1 ) );
512 m_rDa0res.push_back(
new Resplot(
"rDa0_vs_da0", 20, -1.5, 1.5, 100, -0.1, 0.1 ) );
523 m_rDz0res.push_back(
new Resplot(
"rDz0_vs_pt", ptnbins, ptbinlims, 501, -1, 1 ) );
532 m_retares.push_back(
new Resplot(
"reta_vs_ipt", iptnbins, iptbinlims, 2*etaResBins, -wfactor*tmp_absResEta, wfactor*tmp_absResEta ) );
533 m_rphires.push_back(
new Resplot(
"rphi_vs_ipt", iptnbins, iptbinlims, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
534 m_rzedres.push_back(
new Resplot(
"rzed_vs_ipt", iptnbins, iptbinlims, 8*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
535 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_ipt", iptnbins, iptbinlims, 24*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
536 m_riptres.push_back(
new Resplot(
"ript_vs_ipt", iptnbins, iptbinlims, 16*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
537 m_rptres.push_back(
new Resplot(
"rpt_vs_ipt", iptnbins, iptbinlims, 8*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
538 m_rd0res.push_back(
new Resplot(
"rd0_vs_ipt", iptnbins, iptbinlims, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
541 m_retares.push_back(
new Resplot(
"reta_vs_pt", ptnbins, ptbinlims, 8*etaResBins, -wfactor*tmp_absResEta, wfactor*tmp_absResEta ) );
542 m_rphires.push_back(
new Resplot(
"rphi_vs_pt", ptnbins, ptbinlims, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
543 m_rzedres.push_back(
new Resplot(
"rzed_vs_pt", ptnbins, ptbinlims, 24*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
544 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_pt", ptnbins, ptbinlims, 24*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
545 m_riptres.push_back(
new Resplot(
"ript_vs_pt", ptnbins, ptbinlims, 16*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
546 m_rptres.push_back(
new Resplot(
"rpt_vs_pt", ptnbins, ptbinlims, 8*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
547 m_rd0res.push_back(
new Resplot(
"rd0_vs_pt", ptnbins, ptbinlims, factor*24*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
549 m_retares.push_back(
new Resplot(
"reta_vs_ET", ptnbins, ptbinlims, 8*etaResBins, -wfactor*tmp_absResEta, wfactor*tmp_absResEta ) );
550 m_rphires.push_back(
new Resplot(
"rphi_vs_ET", ptnbins, ptbinlims, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
551 m_rzedres.push_back(
new Resplot(
"rzed_vs_ET", ptnbins, ptbinlims, 24*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
552 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_ET", ptnbins, ptbinlims, 24*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
553 m_riptres.push_back(
new Resplot(
"ript_vs_ET", ptnbins, ptbinlims, 16*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
554 m_rptres.push_back(
new Resplot(
"rpt_vs_ET", ptnbins, ptbinlims, 8*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
555 m_rd0res.push_back(
new Resplot(
"rd0_vs_ET", ptnbins, ptbinlims, factor*24*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
559 m_retares.push_back(
new Resplot(
"reta_vs_eta",
etaBins, -tmp_maxEta, tmp_maxEta, 4*etaResBins, -wfactor*tmp_absResEta, wfactor*tmp_absResEta ) );
560 m_rphires.push_back(
new Resplot(
"rphi_vs_eta",
etaBins, -tmp_maxEta, tmp_maxEta, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
561 m_rzedres.push_back(
new Resplot(
"rzed_vs_eta",
etaBins, -tmp_maxEta, tmp_maxEta, 12*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
563 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_eta",
etaBins, -tmp_maxEta, tmp_maxEta, 24*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
565 m_riptres.push_back(
new Resplot(
"ript_vs_eta",
etaBins, -tmp_maxEta, tmp_maxEta, 16*pTResBins, -tmp_absResPt, tmp_absResPt ) );
566 m_rptres.push_back(
new Resplot(
"rpt_vs_eta",
etaBins, -tmp_maxEta, tmp_maxEta, 8*pTResBins, -tmp_absResPt, tmp_absResPt ) );
567 m_rd0res.push_back(
new Resplot(
"rd0_vs_eta",
etaBins, -tmp_maxEta, tmp_maxEta, factor*24*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
580 for (
unsigned ih=0 ; ih<3 ; ih++ ) {
585 m_retares.push_back(
new Resplot(
"reta_vs_zed", 0.2*
zBins, -zMax, zMax, 2*etaResBins, -tmp_absResEta, tmp_absResEta ) );
586 m_rphires.push_back(
new Resplot(
"rphi_vs_zed", 0.2*
zBins, -zMax, zMax, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
587 m_rzedres.push_back(
new Resplot(
"rzed_vs_zed", 0.2*
zBins, -zMax, zMax, 8*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
588 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_zed", 0.2*
zBins, -zMax, zMax, 24*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
591 m_riptres.push_back(
new Resplot(
"ript_vs_zed", 0.2*
zBins, -zMax, zMax, 2*pTResBins, -2*tmp_absResPt, 2*tmp_absResPt ) );
592 m_rptres.push_back(
new Resplot(
"rpt_vs_zed", 0.2*
zBins, -zMax, zMax, 8*pTResBins, -tmp_absResPt, tmp_absResPt ) );
593 m_rd0res.push_back(
new Resplot(
"rd0_vs_zed", 0.2*
zBins, -zMax, zMax, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
596 m_retares.push_back(
new Resplot(
"reta_vs_nvtx", 24, 0, 72, 4*etaResBins, -tmp_absResEta, tmp_absResEta ) );
597 m_rphires.push_back(
new Resplot(
"rphi_vs_nvtx", 24, 0, 72, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
598 m_rzedres.push_back(
new Resplot(
"rzed_vs_nvtx", 24, 0, 72, 4*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
599 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_nvtx", 24, 0, 72, 24*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
600 m_riptres.push_back(
new Resplot(
"ript_vs_nvtx", 24, 0, 72, 4*pTResBins, -tmp_absResPt, tmp_absResPt ) );
601 m_rptres.push_back(
new Resplot(
"rpt_vs_nvtx", 24, 0, 72, 8*pTResBins, -tmp_absResPt, tmp_absResPt ) );
602 m_rd0res.push_back(
new Resplot(
"rd0_vs_nvtx", 24, 0, 72, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
605 m_retares.push_back(
new Resplot(
"reta_vs_ntracks", 60, 0, 600, 4*etaResBins, -tmp_absResEta, tmp_absResEta ) );
606 m_rphires.push_back(
new Resplot(
"rphi_vs_ntracks", 60, 0, 600, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
607 m_rzedres.push_back(
new Resplot(
"rzed_vs_ntracks", 60, 0, 600, 4*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
608 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_ntracks", 60, 0, 600, 24*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
610 m_riptres.push_back(
new Resplot(
"ript_vs_ntracks", 60, 0, 600, 4*pTResBins, -tmp_absResPt, tmp_absResPt ) );
611 m_rptres.push_back(
new Resplot(
"rpt_vs_ntracks", 60, 0, 600, 8*pTResBins, -tmp_absResPt, tmp_absResPt ) );
612 m_rd0res.push_back(
new Resplot(
"rd0_vs_ntracks", 60, 0, 600, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
615 m_retares.push_back(
new Resplot(
"reta_vs_phi", 128, -
M_PI,
M_PI, 2*etaResBins, -wfactor*tmp_absResEta, wfactor*tmp_absResEta ) );
616 m_rphires.push_back(
new Resplot(
"rphi_vs_phi", 128, -
M_PI,
M_PI, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
617 m_rzedres.push_back(
new Resplot(
"rzed_vs_phi", 128, -
M_PI,
M_PI, 8*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
619 m_riptres.push_back(
new Resplot(
"ript_vs_phi", 128, -
M_PI,
M_PI, 16*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
620 m_rptres.push_back(
new Resplot(
"rpt_vs_phi", 128, -
M_PI,
M_PI, 8*pTResBins, -wfactor*tmp_absResPt, wfactor*tmp_absResPt ) );
621 m_rd0res.push_back(
new Resplot(
"rd0_vs_phi", 128, -
M_PI,
M_PI, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
625 m_retares.push_back(
new Resplot(
"reta_vs_mu", 24, 0, 72, 4*etaResBins, -tmp_absResEta, tmp_absResEta ) );
626 m_rphires.push_back(
new Resplot(
"rphi_vs_mu", 24, 0, 72, 8*phiResBins, -wfactor*tmp_absResPhi, wfactor*tmp_absResPhi ) );
627 m_rzedres.push_back(
new Resplot(
"rzed_vs_mu", 24, 0, 72, 4*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
628 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_mu", 24, 0, 72, 24*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
629 m_riptres.push_back(
new Resplot(
"ript_vs_mu", 24, 0, 72, 4*pTResBins, -tmp_absResPt, tmp_absResPt ) );
630 m_rptres.push_back(
new Resplot(
"rpt_vs_mu", 24, 0, 72, 8*pTResBins, -tmp_absResPt, tmp_absResPt ) );
631 m_rd0res.push_back(
new Resplot(
"rd0_vs_mu", 24, 0, 72, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
636 m_rzedres.push_back(
new Resplot(
"rzed_vs_signed_pt", ptnbins2, ptbinlims2, 4*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
637 m_rzedres.push_back(
new Resplot(
"rzed_vs_ABS_pt", ptnbins, ptbinlims, 4*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
638 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_signed_pt", ptnbins2, ptbinlims2, 24*zfactor*zresBins, -zfactor*zresMax, zfactor*zresMax ) );
639 m_rzedthetares.push_back(
new Resplot(
"rzedtheta_vs_ABS_pt", ptnbins, ptbinlims, 24*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax ) );
642 m_rd0res.push_back(
new Resplot(
"rd0_vs_signed_pt", ptnbins2, ptbinlims2, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
643 m_rd0res.push_back(
new Resplot(
"rd0_vs_ABS_pt", ptnbins, ptbinlims, factor*8*a0resBins, -wfactor*a0resMax, wfactor*a0resMax ) );
649 m_retaresPull.push_back(
new Resplot(
"retaPull_vs_ipt", iptnbins, iptbinlims, 2*etaResBins, -5,5));
650 m_rphiresPull.push_back(
new Resplot(
"rphiPull_vs_ipt", iptnbins, iptbinlims, 8*phiResBins, -5,5));
651 m_rzedresPull.push_back(
new Resplot(
"rzedPull_vs_ipt", iptnbins, iptbinlims, 4*zfactor*zresBins, -5,5));
653 m_rptresPull.push_back(
new Resplot(
"rptPull_vs_ipt", iptnbins, iptbinlims, 8*pTResBins, -5,5));
654 m_rd0resPull.push_back(
new Resplot(
"rd0Pull_vs_ipt", iptnbins, iptbinlims, factor*8*a0resBins, -5,5));
656 m_retaresPull.push_back(
new Resplot(
"retaPull_vs_pt", ptnbins, ptbinlims, 2*etaResBins, -5,5));
657 m_rphiresPull.push_back(
new Resplot(
"rphiPull_vs_pt", ptnbins, ptbinlims, 8*phiResBins, -5,5));
658 m_rzedresPull.push_back(
new Resplot(
"rzedPull_vs_pt", ptnbins, ptbinlims, 4*zfactor*zresBins, -5,5));
661 m_rptresPull.push_back(
new Resplot(
"rptPull_vs_pt", ptnbins, ptbinlims, 8*pTResBins, -5,5));
662 m_rd0resPull.push_back(
new Resplot(
"rd0Pull_vs_pt", ptnbins, ptbinlims, factor*8*a0resBins, -5,5));
687 m_rd0resPull.push_back(
new Resplot(
"rd0Pull_vs_nvtx", 12, 0, 36, factor*8*a0resBins, -5,5));
691 m_rzedresPull.push_back(
new Resplot(
"rzedPull_vs_ntracks", 60, 0, 600, 4*zfactor*zresBins, -5,5));
695 m_rd0resPull.push_back(
new Resplot(
"rd0Pull_vs_ntracks", 60, 0, 600, factor*8*a0resBins, -5,5));
697 m_rzedresPull.push_back(
new Resplot(
"rzedPull_vs_signed_pt", ptnbins2, ptbinlims2, 4*zfactor*zresBins, -5,5));
698 m_rzedresPull.push_back(
new Resplot(
"rzedPull_vs_ABS_pt", ptnbins, ptbinlims, 4*zfactor*zresBins, -5,5));
701 m_rd0resPull.push_back(
new Resplot(
"rd0Pull_vs_signed_pt", ptnbins2, ptbinlims2, factor*8*a0resBins, -5,5));
702 m_rd0resPull.push_back(
new Resplot(
"rd0Pull_vs_ABS_pt", ptnbins, ptbinlims, factor*8*a0resBins, -5,5));
704 m_rzedreslb =
new Resplot(
"rzed_vs_lb", 301, -0.5, 3009.5, 8*zfactor*zresBins, -2*zfactor*zresMax, 2*zfactor*zresMax );
706 m_rzedlb =
new Resplot(
"zed_vs_lb", 301, -0.5, 3009.5, 100, -300, 300 );
729 addHistogram(
new TH1F(
"nsct_rec",
"nsct_rec", NHits, -0.5,
float(NHits-0.5) ) );
732 addHistogram(
new TH1F(
"npix_rec",
"npix_rec", NHits, -0.5,
float(NHits-0.5) ) );
735 addHistogram(
new TH1F(
"nsi_rec",
"nsi_rec", NHits, -0.5,
float(NHits-0.5) ) );
736 addHistogram(
new TH1F(
"nsi_matched",
"nsi_matched", NHits, -0.5,
float(NHits-0.5) ) );
740 addHistogram(
new TH1F(
"ntrt_rec",
"ntrt_rec", NHits, -0.5,
float(NHits-0.5) ) );
742 addHistogram(
new TH1F(
"nstraw",
"nstraw", NHits*4, -0.5,
float(4*NHits-0.5) ) );
743 addHistogram(
new TH1F(
"nstraw_rec",
"nstraw_rec", NHits*4, -0.5,
float(4*NHits-0.5) ) );
745 addHistogram(
new TH1F(
"ntracks",
"ntracks", Ntracks+1, -0.5,
float(Ntracks+0.5) ) );
746 addHistogram(
new TH1F(
"ntracks_rec",
"ntracks_rec", Ntracks+1, -0.5,
float(Ntracks+0.5) ) );
750 m_h2 =
new Resplot(
"d0vphi", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
751 m_h2r =
new Resplot(
"d0vphi_rec", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
752 m_h2m =
new Resplot(
"d0vphi_match", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
754 m_h2a0 =
new Resplot(
"a0vphi", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
755 m_h2a0r =
new Resplot(
"a0vphi_rec", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
760 TH1F heffvlb(
"eff vs lb",
"eff vs lb", 301, -0.5, 3009.5 );
795 TH1F* eff_vs_mult =
new TH1F(
"eff_vs_mult",
"eff_vs_mult", 25, 0, 25 );
815 double nbins[23] = { 0, 29.5, 39.5, 48.5, 56.5, 63.5, 70.5, 77.5,
816 83.5, 91.5, 97.5, 106.5, 114.5, 121.5, 129.5, 140.5,
817 151.5, 164.5, 200.5, 250.5, 300.5, 400.5, 600 };
819 TH1F* n_vtx_tracks2 =
new TH1F(
"nvtxtracks2",
"nvtxtracks2", 22,
nbins);
821 delete n_vtx_tracks2;
823 m_n_vtx =
new TH1F(
"nvtx",
"nvtx", 81, -0.5, 80.5);
826 m_mu =
new TH1F(
"mu",
"mu", 90, 0, 90);
832 double etovpt_bins[39] = {
833 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,
834 1, 1.08571, 1.17877, 1.2798, 1.3895, 1.50859, 1.63789, 1.77828, 1.9307, 2.09618,
835 2.27585, 2.47091, 2.6827, 2.91263, 3.16228, 3.43332, 3.72759, 4.04709, 4.39397, 4.77058,
836 5.17947, 5.62341, 6.1054, 6.6287, 7.19686, 7.81371, 8.48343, 9.21055, 10
841 m_etovpt =
new TH1F(
"etovpt",
"ET / pT", 38, etovpt_bins );
845 m_et =
new TH1F(
"ET",
"ET; E_{T} [GeV]", ptnbins, ptbinlims );