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
0008
0009
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 };
0020 double eff_target = 0.9;
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
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
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
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
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
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
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
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 }