Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2025-01-18 10:18:38

0001 #include "RiceStyle.h"
0002 #include "common_bench/benchmark.h"
0003 
0004 using namespace std;
0005 
0006 int setbenchstatus(double eff){
0007     ///////////// Set benchmark status!
0008         // create our test definition
0009         // test_tag
0010         common_bench::Test rho_reco_eff_test{
0011           {
0012             {"name", "rho_reconstruction_efficiency"},
0013             {"title", "rho Reconstruction Efficiency for rho -> pi+pi- in the B0"},
0014             {"description", "u-channel rho->pi+pi- reconstruction efficiency "
0015                        "when both pions should be within B0 acceptance"},
0016             {"quantity", "efficiency"},
0017             {"target", "0.9"}
0018           }
0019         };      //these 2 need to be consistent 
0020         double eff_target = 0.9;    //going to find a way to use the same variable
0021 
0022         if(eff<0 || eff>1){
0023           rho_reco_eff_test.error(-1);
0024         }else if(eff > eff_target){
0025           rho_reco_eff_test.pass(eff);
0026         }else{
0027           rho_reco_eff_test.fail(eff);
0028         }
0029 
0030         // write out our test data
0031         common_bench::write_test(rho_reco_eff_test, "./benchmark_output/u_rho_eff.json");
0032     return 0;
0033 }
0034 
0035 void plot_rho_physics_benchmark(TString filename="./sim_output/plot_combined.root"){
0036   Ssiz_t dotPosition = filename.Last('.');
0037   TString figure_directory = filename(0, dotPosition);
0038   figure_directory += "_figures";   
0039   
0040   TFile* file = new TFile(filename);
0041   TString vm_label="#rho^{0}";
0042   TString daug_label="#pi^{#plus}#pi^{#minus}";
0043   //mass distribution
0044   TH1D* h_VM_mass_MC = (TH1D*) file->Get("h_VM_mass_MC");
0045   TH1D* h_VM_mass_REC = (TH1D*) file->Get("h_VM_mass_REC");
0046   TH1D* h_VM_mass_REC_justpions = (TH1D*) file->Get("h_VM_mass_REC_justpions");
0047   //mass distribution within B0
0048   TH1D* h_VM_mass_MC_etacut = (TH1D*) file->Get("h_VM_mass_MC_etacut");
0049   TH1D* h_VM_mass_REC_etacut = (TH1D*) file->Get("h_VM_mass_REC_etacut");
0050   //efficiencies
0051   TProfile2D* h_effEtaPtPi = (TProfile2D*) file->Get("h_effEtaPtPi");
0052   TProfile2D* h_effEtaPtPip = (TProfile2D*) file->Get("h_effEtaPtPip"); 
0053   TProfile2D* h_effEtaPtPim = (TProfile2D*) file->Get("h_effEtaPtPim"); 
0054   TProfile2D* h_effPhiEtaPi = (TProfile2D*) file->Get("h_effPhiEtaPi"); 
0055   TProfile2D* h_effPhiEtaPip = (TProfile2D*) file->Get("h_effPhiEtaPip"); 
0056   TProfile2D* h_effPhiEtaPim = (TProfile2D*) file->Get("h_effPhiEtaPim");   
0057 
0058 
0059   TLatex* r42 = new TLatex(0.18, 0.91, "ep 10#times100 GeV");
0060   r42->SetNDC();
0061   r42->SetTextSize(22);
0062   r42->SetTextFont(43);
0063   r42->SetTextColor(kBlack);
0064 
0065   TLatex* r43 = new TLatex(0.9,0.91, "#bf{EPIC}");
0066   r43->SetNDC();
0067   //r43->SetTextSize(0.04);
0068   r43->SetTextSize(22);
0069 
0070   TLatex* r44 = new TLatex(0.53, 0.78, "10^{-3}<Q^{2}<10 GeV^{2}, W>2 GeV");
0071   r44->SetNDC();
0072   r44->SetTextSize(20);
0073   r44->SetTextFont(43);
0074   r44->SetTextColor(kBlack);
0075 
0076   TLatex* r44_2 = new TLatex(0.5, 0.83, ""+vm_label+" #rightarrow "+daug_label+" eSTARlight");
0077   r44_2->SetNDC();
0078   r44_2->SetTextSize(30);
0079   r44_2->SetTextFont(43);
0080   r44_2->SetTextColor(kBlack);
0081 
0082   TCanvas* c2 = new TCanvas("c2","c2",1,1,600,600);
0083   gPad->SetTicks();
0084   gPad->SetLeftMargin(0.18);
0085   gPad->SetBottomMargin(0.18);
0086   gPad->SetTopMargin(0.10);
0087   gPad->SetRightMargin(0.01);
0088   TH1D* base2 = makeHist("base2", "", "#pi^{#plus}#pi^{#minus} inv. mass (GeV)", "counts", 100,0.05,2.05,kBlack);
0089   base2->GetYaxis()->SetRangeUser(0.5, 1.2*(h_VM_mass_MC->GetMaximum()));
0090   base2->GetXaxis()->SetTitleColor(kBlack);
0091   fixedFontHist1D(base2,1.,1.2);
0092   base2->GetYaxis()->SetTitleSize(base2->GetYaxis()->GetTitleSize()*1.5);
0093   base2->GetXaxis()->SetTitleSize(base2->GetXaxis()->GetTitleSize()*1.5);
0094   base2->GetYaxis()->SetLabelSize(base2->GetYaxis()->GetLabelSize()*1.5);
0095   base2->GetXaxis()->SetLabelSize(base2->GetXaxis()->GetLabelSize()*1.5);
0096   base2->GetXaxis()->SetNdivisions(4,4,0);
0097   base2->GetYaxis()->SetNdivisions(5,5,0);
0098   base2->GetYaxis()->SetTitleOffset(1.3);
0099   base2->Draw();
0100 
0101   TH1D* h_VM_mass_REC_justprotons = (TH1D*)h_VM_mass_REC->Clone("h_VM_mass_REC_justprotons");
0102   for(int ibin=1; ibin<h_VM_mass_REC_justprotons->GetNbinsX(); ibin++){
0103     h_VM_mass_REC_justprotons->SetBinContent(ibin,h_VM_mass_REC_justprotons->GetBinContent(ibin) - h_VM_mass_REC_justpions->GetBinContent(ibin));
0104   }
0105 
0106   h_VM_mass_MC->SetFillColorAlpha(kBlack,0.2);
0107   h_VM_mass_REC->SetFillColorAlpha(kMagenta,0.2);
0108   h_VM_mass_MC->SetLineColor(kBlack);
0109   h_VM_mass_REC->SetLineColor(kMagenta);
0110   h_VM_mass_REC_justpions->SetLineColor(kViolet+10);
0111   h_VM_mass_REC_justprotons->SetLineColor(kRed);
0112   h_VM_mass_MC->SetLineWidth(2);
0113   h_VM_mass_REC->SetLineWidth(2);
0114   h_VM_mass_REC_justpions->SetLineWidth(2);
0115   h_VM_mass_REC_justprotons->SetLineWidth(2);
0116 
0117   h_VM_mass_REC->Scale(3.0);
0118   h_VM_mass_REC_justpions->Scale(3.0);
0119   h_VM_mass_REC_justprotons->Scale(3.0);
0120 
0121   h_VM_mass_MC->Draw("HIST E same");
0122   h_VM_mass_REC->Draw("HIST E same");
0123   h_VM_mass_REC_justpions->Draw("HIST same");
0124   h_VM_mass_REC_justprotons->Draw("HIST same");
0125 
0126   r42->Draw("same");
0127   r43->Draw("same");
0128   r44->Draw("same");
0129   r44_2->Draw("same");
0130 
0131   TLegend *w8 = new TLegend(0.58,0.63,0.93,0.76);
0132   w8->SetLineColor(kWhite);
0133   w8->SetFillColor(0);
0134   w8->SetTextSize(17);
0135   w8->SetTextFont(45);
0136   w8->AddEntry(h_VM_mass_MC, ""+vm_label+" MC ", "L");
0137   w8->AddEntry(h_VM_mass_REC, vm_label+" reco.#times3", "L");
0138   w8->AddEntry(h_VM_mass_REC_justpions, vm_label+" reco.#times3 (#pi^{#minus}#pi^{#plus})", "L");
0139   w8->AddEntry(h_VM_mass_REC_justprotons, vm_label+" reco.#times3 (#pi^{#minus}p)", "L");
0140   w8->Draw("same");
0141 
0142   TString figure1name = figure_directory+"/benchmark_rho_mass.pdf";
0143   c2->Print(figure1name);
0144 
0145 
0146   ///////////////////// Figure 4
0147   TCanvas* c4 = new TCanvas("c4","c4",1,1,600,600);
0148   gPad->SetTicks();
0149   gPad->SetLeftMargin(0.18);
0150   gPad->SetBottomMargin(0.18);
0151   gPad->SetTopMargin(0.10);
0152   gPad->SetRightMargin(0.01);
0153   TH1D* base4 = makeHist("base4", "", "#pi^{#plus}#pi^{#minus} inv. mass (GeV)", "counts", 100,0.05,2.05,kBlack);
0154   base4->GetYaxis()->SetRangeUser(0.5, 1.2*(h_VM_mass_MC_etacut->GetMaximum()));
0155   base4->GetXaxis()->SetTitleColor(kBlack);
0156   fixedFontHist1D(base4,1.,1.2);
0157   base4->GetYaxis()->SetTitleSize(base4->GetYaxis()->GetTitleSize()*1.5);
0158   base4->GetXaxis()->SetTitleSize(base4->GetXaxis()->GetTitleSize()*1.5);
0159   base4->GetYaxis()->SetLabelSize(base4->GetYaxis()->GetLabelSize()*1.5);
0160   base4->GetXaxis()->SetLabelSize(base4->GetXaxis()->GetLabelSize()*1.5);
0161   base4->GetXaxis()->SetNdivisions(4,4,0);
0162   base4->GetYaxis()->SetNdivisions(5,5,0);
0163   base4->GetYaxis()->SetTitleOffset(1.3);
0164   base4->Draw();
0165 
0166   h_VM_mass_MC_etacut->SetFillColorAlpha(kBlack,0.2);
0167   h_VM_mass_REC_etacut->SetFillColorAlpha(kMagenta,0.2);
0168   h_VM_mass_MC_etacut->SetLineColor(kBlack);
0169   h_VM_mass_REC_etacut->SetLineColor(kMagenta);
0170   h_VM_mass_MC_etacut->SetLineWidth(2);
0171   h_VM_mass_REC_etacut->SetLineWidth(2);
0172 
0173   h_VM_mass_MC_etacut->Draw("HIST E same");
0174   h_VM_mass_REC_etacut->Draw("HIST E same");
0175 
0176   int minbineff = h_VM_mass_MC_etacut->FindBin(0.6);
0177   int maxbineff = h_VM_mass_MC_etacut->FindBin(1.0);
0178   double thiseff = 100.0*(1.0*h_VM_mass_REC_etacut->Integral(minbineff,maxbineff))/(1.0*h_VM_mass_MC_etacut->Integral(minbineff,maxbineff));
0179 
0180   r42->Draw("same");
0181   r43->Draw("same");
0182   r44->Draw("same");
0183   r44_2->Draw("same");
0184 
0185   TLegend *w10 = new TLegend(0.58,0.62,0.93,0.7);
0186   w10->SetLineColor(kWhite);
0187   w10->SetFillColor(0);
0188   w10->SetTextSize(17);
0189   w10->SetTextFont(45);
0190   w10->AddEntry(h_VM_mass_MC_etacut, vm_label+" MC ", "L");
0191   w10->AddEntry(h_VM_mass_REC_etacut, vm_label+" reco. (#pi^{#minus}#pi^{#plus})", "L");
0192   w10->Draw("same");
0193 
0194   TLatex* anglelabel = new TLatex(0.59, 0.73, "9<#theta_{#pi^{#pm},MC}<13 mrad");
0195   anglelabel->SetNDC();
0196   anglelabel->SetTextSize(20);
0197   anglelabel->SetTextFont(43);
0198   anglelabel->SetTextColor(kBlack);
0199   anglelabel->Draw("same");
0200 
0201   TLatex* efflabel = new TLatex(0.59, 0.55, "reco. eff (0.6<m^{2}<1 GeV^{2})");
0202   efflabel->SetNDC();
0203   efflabel->SetTextSize(20);
0204   efflabel->SetTextFont(43);
0205   efflabel->SetTextColor(kBlack);
0206   efflabel->Draw("same");
0207   TLatex* effnlabel = new TLatex(0.59, 0.51, Form("          = %.0f%%",thiseff));
0208   effnlabel->SetNDC();
0209   effnlabel->SetTextSize(20);
0210   effnlabel->SetTextFont(43);
0211   effnlabel->SetTextColor(kBlack);
0212   effnlabel->Draw("same");
0213 
0214   TString figure2name = figure_directory+"/benchmark_rho_mass_cuts.pdf";
0215   c4->Print(figure2name);
0216 
0217   ///////////////////// Figure 5
0218   TCanvas* c5 = new TCanvas("c5","c5",1,1,700,560);
0219   TPad* p5 = new TPad("p5","Pad5",0.,0.,1.,1.);
0220   p5->Divide(3,2,0,0);
0221   p5->Draw();
0222   gStyle->SetPalette(kBlueRedYellow);
0223   gStyle->SetOptStat(0);    
0224 
0225   h_effEtaPtPi  ->GetXaxis()->SetLabelSize(h_effEtaPtPi  ->GetXaxis()->GetLabelSize()*1.8);
0226   h_effEtaPtPip ->GetXaxis()->SetLabelSize(h_effEtaPtPip ->GetXaxis()->GetLabelSize()*1.8);
0227   h_effEtaPtPim ->GetXaxis()->SetLabelSize(h_effEtaPtPim ->GetXaxis()->GetLabelSize()*1.8);
0228   h_effEtaPtPi  ->GetYaxis()->SetLabelSize(h_effEtaPtPi  ->GetYaxis()->GetLabelSize()*1.8);
0229   h_effEtaPtPim ->GetZaxis()->SetLabelSize(h_effEtaPtPim ->GetZaxis()->GetLabelSize()*0.5);
0230   h_effEtaPtPim ->GetZaxis()->SetTitleSize(h_effEtaPtPim ->GetZaxis()->GetTitleSize()*0.5);
0231   h_effPhiEtaPi ->GetXaxis()->SetLabelSize(h_effPhiEtaPi ->GetXaxis()->GetLabelSize()*1.8);
0232   h_effPhiEtaPip->GetXaxis()->SetLabelSize(h_effPhiEtaPip->GetXaxis()->GetLabelSize()*1.8);
0233   h_effPhiEtaPim->GetXaxis()->SetLabelSize(h_effPhiEtaPim->GetXaxis()->GetLabelSize()*1.8);
0234   h_effPhiEtaPi ->GetYaxis()->SetLabelSize(h_effPhiEtaPi ->GetYaxis()->GetLabelSize()*1.8);
0235   h_effPhiEtaPim->GetZaxis()->SetLabelSize(h_effPhiEtaPim->GetZaxis()->GetLabelSize()*0.5);
0236   h_effPhiEtaPim->GetZaxis()->SetTitleSize(h_effPhiEtaPim->GetZaxis()->GetTitleSize()*0.5);
0237 
0238   fixedFontHist1D(h_effEtaPtPi,1.,1.2);
0239   fixedFontHist1D(h_effEtaPtPip,1.,1.2);
0240   fixedFontHist1D(h_effEtaPtPim,1.,1.2);
0241   fixedFontHist1D(h_effPhiEtaPi,1.,1.2);
0242   fixedFontHist1D(h_effPhiEtaPip,1.,1.2);
0243   fixedFontHist1D(h_effPhiEtaPim,1.,1.2);
0244 
0245   p5->cd(1);
0246   TVirtualPad* p51 = p5->cd(1);
0247   p51->SetTopMargin(0.08);
0248   p51->SetRightMargin(0);
0249   p51->SetLeftMargin(0.21);
0250   p51->SetBottomMargin(0.2);
0251   h_effEtaPtPi->GetXaxis()->SetRangeUser(3.9,6.05);
0252   h_effEtaPtPi->GetYaxis()->SetRangeUser(0,1.7);
0253   h_effEtaPtPi->GetZaxis()->SetRangeUser(0,1);
0254   h_effEtaPtPi->GetXaxis()->SetNdivisions(5);
0255   h_effEtaPtPi->GetYaxis()->SetNdivisions(5);
0256   h_effEtaPtPi->SetContour(99);
0257   h_effEtaPtPi->Draw("COLZ");
0258   TLatex* pilabel = new TLatex(0.81, 0.75, "#pi^{#pm}");
0259   pilabel->SetNDC();
0260   pilabel->SetTextSize(40);
0261   pilabel->SetTextFont(43);
0262   pilabel->SetTextColor(kBlack);
0263   pilabel->Draw("same");
0264   TLatex* r44fig5c = new TLatex(0.21, 0.93, "ep 10#times100 GeV         #rho^{0}#rightarrow#pi^{#plus}#pi^{#minus}");
0265   r44fig5c->SetNDC();
0266   r44fig5c->SetTextSize(15);
0267   r44fig5c->SetTextFont(43);
0268   r44fig5c->SetTextColor(kBlack);
0269   r44fig5c->Draw("same");
0270 
0271   p5->cd(2);
0272   TVirtualPad* p52 = p5->cd(2);
0273   p52->SetTopMargin(0.08);
0274   p52->SetRightMargin(0);
0275   p52->SetLeftMargin(0);
0276   p52->SetBottomMargin(0.2);
0277   h_effEtaPtPip->GetXaxis()->SetRangeUser(4.05,6.05);
0278   h_effEtaPtPip->GetYaxis()->SetRangeUser(0,1.7);
0279   h_effEtaPtPip->GetZaxis()->SetRangeUser(0,1);
0280   h_effEtaPtPip->GetXaxis()->SetNdivisions(5);
0281   h_effEtaPtPip->GetYaxis()->SetNdivisions(5);
0282   h_effEtaPtPip->SetContour(99);
0283   h_effEtaPtPip->Draw("COLZ");
0284   TLatex* piplabel = new TLatex(0.81, 0.75, "#pi^{#plus}");
0285   piplabel->SetNDC();
0286   piplabel->SetTextSize(40);
0287   piplabel->SetTextFont(43);
0288   piplabel->SetTextColor(kBlack);
0289   piplabel->Draw("same");
0290   TLatex* r44fig5a = new TLatex(0.01, 0.93, "eSTARlight        10^{-3}<Q^{2}<10 GeV^{2}");
0291   r44fig5a->SetNDC();
0292   r44fig5a->SetTextSize(15);
0293   r44fig5a->SetTextFont(43);
0294   r44fig5a->SetTextColor(kBlack);
0295   r44fig5a->Draw("same");
0296 
0297   p5->cd(3);
0298   TVirtualPad* p53 = p5->cd(3);
0299   p53->SetTopMargin(0.08);
0300   p53->SetRightMargin(0.2);
0301   p53->SetLeftMargin(0);
0302   p53->SetBottomMargin(0.2);
0303   h_effEtaPtPim->SetTitle(";#eta;;efficiency");
0304   h_effEtaPtPim->GetXaxis()->SetRangeUser(4.05,6.05);
0305   h_effEtaPtPim->GetYaxis()->SetRangeUser(0,1.7);
0306   h_effEtaPtPim->GetZaxis()->SetRangeUser(0,1);
0307   h_effEtaPtPim->GetXaxis()->SetNdivisions(5);
0308   h_effEtaPtPim->GetYaxis()->SetNdivisions(5);
0309   h_effEtaPtPim->SetContour(99);
0310   h_effEtaPtPim->Draw("COLZ");
0311   TLatex* pimlabel = new TLatex(0.62, 0.75, "#pi^{#minus}");
0312   pimlabel->SetNDC();
0313   pimlabel->SetTextSize(40);
0314   pimlabel->SetTextFont(43);
0315   pimlabel->SetTextColor(kBlack);
0316   pimlabel->Draw("same");
0317   TLatex* r43fig5 = new TLatex(0.66,0.93, "#bf{EPIC}");
0318   r43fig5->SetNDC();
0319   r43fig5->SetTextSize(15);
0320   r43fig5->SetTextFont(43);
0321   r43fig5->SetTextColor(kBlack);
0322   r43fig5->Draw("same");
0323   TLatex* r44fig5b = new TLatex(0.01, 0.93, "W>2 GeV");
0324   r44fig5b->SetNDC();
0325   r44fig5b->SetTextSize(15);
0326   r44fig5b->SetTextFont(43);
0327   r44fig5b->SetTextColor(kBlack);
0328   r44fig5b->Draw("same");
0329 
0330   p5->cd(4);
0331   TVirtualPad* p54 = p5->cd(4);
0332   p54->SetTopMargin(0.05);
0333   p54->SetRightMargin(0);
0334   p54->SetLeftMargin(0.2);
0335   p54->SetBottomMargin(0.21);
0336   h_effPhiEtaPi->GetXaxis()->SetRangeUser(0,6.2);
0337   h_effPhiEtaPi->GetYaxis()->SetRangeUser(4,6);
0338   h_effPhiEtaPi->GetZaxis()->SetRangeUser(0,1);
0339   h_effPhiEtaPi->GetXaxis()->SetNdivisions(8);
0340   h_effPhiEtaPi->GetYaxis()->SetNdivisions(5);
0341   h_effPhiEtaPi->SetContour(99);
0342   h_effPhiEtaPi->Draw("COLZ");
0343   TLatex* pilabela = new TLatex(0.3, 0.82, "#pi^{#pm}");
0344   TLatex* pilabelb = new TLatex(0.5, 0.84, "(p_{T}>0.2 GeV/c)");
0345   pilabela->SetNDC();
0346   pilabelb->SetNDC();
0347   pilabela->SetTextSize(40);
0348   pilabelb->SetTextSize(15);
0349   pilabela->SetTextFont(43);
0350   pilabelb->SetTextFont(43);
0351   pilabela->SetTextColor(kWhite);
0352   pilabelb->SetTextColor(kWhite);
0353   pilabela->Draw("same");
0354   pilabelb->Draw("same");
0355 
0356   p5->cd(5);
0357   TVirtualPad* p55 = p5->cd(5);
0358   p55->SetTopMargin(0.05);
0359   p55->SetRightMargin(0);
0360   p55->SetLeftMargin(0);
0361   p55->SetBottomMargin(0.2);
0362   h_effPhiEtaPip->GetXaxis()->SetRangeUser(0.15,6.2);
0363   h_effPhiEtaPip->GetYaxis()->SetRangeUser(4,6);
0364   h_effPhiEtaPip->GetZaxis()->SetRangeUser(0,1);
0365   h_effPhiEtaPip->GetXaxis()->SetNdivisions(8);
0366   h_effPhiEtaPip->GetYaxis()->SetNdivisions(5);
0367   h_effPhiEtaPip->SetContour(99);
0368   h_effPhiEtaPip->Draw("COLZ");
0369   TLatex* piplabela = new TLatex(0.2, 0.82, "#pi^{#plus}");
0370   TLatex* piplabelb = new TLatex(0.4, 0.84, "(p_{T}>0.2 GeV/c)");
0371   piplabela->SetNDC();
0372   piplabelb->SetNDC();
0373   piplabela->SetTextSize(40);
0374   piplabelb->SetTextSize(15);
0375   piplabela->SetTextFont(43);
0376   piplabelb->SetTextFont(43);
0377   piplabela->SetTextColor(kWhite);
0378   piplabelb->SetTextColor(kWhite);
0379   piplabela->Draw("same");
0380   piplabelb->Draw("same");
0381 
0382   p5->cd(6);
0383   TVirtualPad* p56 = p5->cd(6);
0384   p56->SetTopMargin(0.05);
0385   p56->SetRightMargin(0.2);
0386   p56->SetLeftMargin(0);
0387   p56->SetBottomMargin(0.2);
0388   h_effPhiEtaPim->SetTitle(";#phi (rad);;efficiency");
0389   h_effPhiEtaPim->GetXaxis()->SetRangeUser(0.15,6.2);
0390   h_effPhiEtaPim->GetYaxis()->SetRangeUser(4,6);
0391   h_effPhiEtaPim->GetZaxis()->SetRangeUser(0,1);
0392   h_effPhiEtaPim->GetXaxis()->SetNdivisions(8);
0393   h_effPhiEtaPim->GetYaxis()->SetNdivisions(5);
0394   h_effPhiEtaPim->SetContour(99);
0395   h_effPhiEtaPim->Draw("COLZ");
0396   TLatex* pimlabela = new TLatex(0.1, 0.82, "#pi^{#minus}");
0397   TLatex* pimlabelb = new TLatex(0.25, 0.84, "(p_{T}>0.2 GeV/c)");
0398   pimlabela->SetNDC();
0399   pimlabelb->SetNDC();
0400   pimlabela->SetTextSize(40);
0401   pimlabelb->SetTextSize(15);
0402   pimlabela->SetTextFont(43);
0403   pimlabelb->SetTextFont(43);
0404   pimlabela->SetTextColor(kWhite);
0405   pimlabelb->SetTextColor(kWhite);
0406   pimlabela->Draw("same");
0407   pimlabelb->Draw("same");
0408 
0409   TString figure3name = figure_directory+"/benchmark_rho_efficiencies.pdf";
0410   c5->Print(figure3name);
0411 
0412   double rhorecoeff = thiseff/100.0;
0413   setbenchstatus(rhorecoeff);
0414 }