/*Compile with: g++ -o SMRDEffPlot $(root-config --libs --cflags) SMRDEffPlot.cxx */ //#include #include #include // interface to read and write #include // container to operate with sequences of characters #include //#include #include #include #include #include #include #include #include #include using namespace std; // group entities under a name int main( int argc, char** argv ){ //void Effpsplot() { if(argc!=3){ /* We print argv[0] assuming it is the program name */ printf( "Wrong number of inputs\n"); return 0; } char buff[100]; //Name of input file //sprintf(buff,"EffFinal_%d.root", first); sprintf(buff,"%s",argv[1]); //general style gStyle->SetOptStat(""); gStyle->SetHistMinimumZero(kTRUE); gStyle->SetOptFit(1); gStyle->SetOptTitle(0); gStyle->SetTitle(0); gStyle->SetPalette(1); // gStyle->SetFillColor(1); gStyle->SetPadColor(0); gStyle->SetCanvasColor(0); gStyle->SetStatColor(0); gStyle->SetTitleColor(1); gStyle->SetTitleSize(0.04); gStyle->SetPadBorderMode(0); gStyle->SetPadBottomMargin(0.15); gStyle->SetPadLeftMargin(0.12); gStyle->SetPadRightMargin(0.09); gStyle->SetFrameBorderMode(0); gStyle->SetCanvasBorderMode(0); //Data HISTOGRAMS TFile f1(buff); //open file // TFile f1("EffFinal_%i.root" % first); //open file f1.cd(); // f1.ls(); TH1F* HtrTin_D =(TH1F*)f1.Get("HtryT_in"); TH1F* HtrBin_D =(TH1F*)f1.Get("HtryB_in"); TH1F* HtrRin_D =(TH1F*)f1.Get("HtryR_in"); TH1F* HtrLin_D =(TH1F*)f1.Get("HtryL_in"); TH1F* HtrTou_D =(TH1F*)f1.Get("HtryT_ou"); TH1F* HtrBou_D =(TH1F*)f1.Get("HtryB_ou"); TH1F* HtrRou_D =(TH1F*)f1.Get("HtryR_ou"); TH1F* HtrLou_D =(TH1F*)f1.Get("HtryL_ou"); TH1F* HtrTc1_D =(TH1F*)f1.Get("HtryT_c1"); TH1F* HtrBc1_D =(TH1F*)f1.Get("HtryB_c1"); TH1F* HtrRc1_D =(TH1F*)f1.Get("HtryR_c1"); TH1F* HtrLc1_D =(TH1F*)f1.Get("HtryL_c1"); TH1F* HtrRc2_D =(TH1F*)f1.Get("HtryR_c2"); TH1F* HtrLc2_D =(TH1F*)f1.Get("HtryL_c2"); TH1F* HtrRc3_D =(TH1F*)f1.Get("HtryR_c3"); TH1F* HtrLc3_D =(TH1F*)f1.Get("HtryL_c3"); TH1F* HtrRc4_D =(TH1F*)f1.Get("HtryR_c4"); TH1F* HtrLc4_D =(TH1F*)f1.Get("HtryL_c4"); TH1F* HsuTin_D =(TH1F*)f1.Get("HsucT_in"); TH1F* HsuBin_D =(TH1F*)f1.Get("HsucB_in"); TH1F* HsuRin_D =(TH1F*)f1.Get("HsucR_in"); TH1F* HsuLin_D =(TH1F*)f1.Get("HsucL_in"); TH1F* HsuTou_D =(TH1F*)f1.Get("HsucT_ou"); TH1F* HsuBou_D =(TH1F*)f1.Get("HsucB_ou"); TH1F* HsuRou_D =(TH1F*)f1.Get("HsucR_ou"); TH1F* HsuLou_D =(TH1F*)f1.Get("HsucL_ou"); TH1F* HsuTc1_D =(TH1F*)f1.Get("HsucT_c1"); TH1F* HsuBc1_D =(TH1F*)f1.Get("HsucB_c1"); TH1F* HsuRc1_D =(TH1F*)f1.Get("HsucR_c1"); TH1F* HsuLc1_D =(TH1F*)f1.Get("HsucL_c1"); TH1F* HsuRc2_D =(TH1F*)f1.Get("HsucR_c2"); TH1F* HsuLc2_D =(TH1F*)f1.Get("HsucL_c2"); TH1F* HsuRc3_D =(TH1F*)f1.Get("HsucR_c3"); TH1F* HsuLc3_D =(TH1F*)f1.Get("HsucL_c3"); TH1F* HsuRc4_D =(TH1F*)f1.Get("HsucR_c4"); TH1F* HsuLc4_D =(TH1F*)f1.Get("HsucL_c4"); // TH1F* CQ_fail_D = (TH1F*)f1.Get("CQ_fail"); // TH1F* CQ_good_D = (TH1F*)f1.Get("CQ_good"); // TH1F* TQ_fail_D = (TH1F*)f1.Get("TQ_fail"); // TH1F* TQ_good_D = (TH1F*)f1.Get("TQ_good"); // TH1F* Cmb_Q_D = (TH1F*)f1.Get("Cmb_Q"); // TH1F* Cmb_Qwin_D = (TH1F*)f1.Get("Cmb_Qwin"); // TH1F* Cmb_T_D = (TH1F*)f1.Get("Cmb_T"); // // TH1F* LgrN_D = (TH1F*)f1.Get("LgrN"); // TH1F* Sng_Q_D = (TH1F*)f1.Get("Sng_Q"); // TH1F* Sng_T_D = (TH1F*)f1.Get("Sng_T"); // // TH1F* Sng_T_fine_D = (TH1F*)f1.Get("Sng_T_fine"); // // TH1F* TgrN_D = (TH1F*)f1.Get("TgrN"); // TH1F* cmbN_D = (TH1F*)f1.Get("cmbN"); // TH1F* sngN_D = (TH1F*)f1.Get("sngN"); // TH1F* HitR_Fail_D = (TH1F*)f1.Get("HitR_Fail"); // TH1F* HitT1_Fail_D = (TH1F*)f1.Get("HitT1_Fail"); // TH1F* HitT2_Fail_D = (TH1F*)f1.Get("HitT2_Fail"); // TH1F* HitT3_Fail_D = (TH1F*)f1.Get("HitT3_Fail"); // TH1F* TrgQt_D = (TH1F*)f1.Get("TrgQt"); // TH1F* TrgWin_D = (TH1F*)f1.Get("TrgWin"); // TH1F* tEnd_D = (TH1F*)f1.Get("tEnd"); // TH1F* qEnd_D = (TH1F*)f1.Get("qEnd"); // TH1F* TrgWinQ_D = (TH1F*)f1.Get("TrgWinQ"); // TH1F* HitWin_D = (TH1F*)f1.Get("HitWin"); // TH1F* Recon_TvsQ_D = (TH1F*)f1.Get("Recon_TvsQ"); // TH1F* Recon_TvsQmin_D = (TH1F*)f1.Get("Recon_TvsQmin"); TH1F* HfailTin = new TH1F("HfailT_in","",128,0.5,128.5); HfailTin->Add(HtrTin_D,HsuTin_D,1.,-1.); TH1F* HfailBin = new TH1F("HfailB_in","",128,0.5,128.5); HfailBin->Add(HtrBin_D,HsuBin_D,1.,-1.); TH1F* HfailRin = new TH1F("HfailR_in","",160,0.5,160.5); HfailRin->Add(HtrRin_D,HsuRin_D,1.,-1.); TH1F* HfailLin = new TH1F("HfailL_in","",160,0.5,160.5); HfailLin->Add(HtrLin_D,HsuLin_D,1.,-1.); TH1F* HfailTou = new TH1F("HfailT_ou","",128,0.5,128.5); HfailTou->Add(HtrTou_D,HsuTou_D,1.,-1.); TH1F* HfailBou = new TH1F("HfailB_ou","",128,0.5,128.5); HfailBou->Add(HtrBou_D,HsuBou_D,1.,-1.); TH1F* HfailRou = new TH1F("HfailR_ou","",160,0.5,160.5); HfailRou->Add(HtrRou_D,HsuRou_D,1.,-1.); TH1F* HfailLou = new TH1F("HfailL_ou","",160,0.5,160.5); HfailLou->Add(HtrLou_D,HsuLou_D,1.,-1.); TH1F* HfailTc1 = new TH1F("HfailT_c1","",128,0.5,128.5); HfailTc1->Add(HtrTc1_D,HsuTc1_D,1.,-1.); TH1F* HfailBc1 = new TH1F("HfailB_c1","",128,0.5,128.5); HfailBc1->Add(HtrBc1_D,HsuBc1_D,1.,-1.); TH1F* HfailRc1 = new TH1F("HfailR_c1","",160,0.5,160.5); HfailRc1->Add(HtrRc1_D,HsuRc1_D,1.,-1.); TH1F* HfailLc1 = new TH1F("HfailL_c1","",160,0.5,160.5); HfailLc1->Add(HtrLc1_D,HsuLc1_D,1.,-1.); TH1F* HfailRc2 = new TH1F("HfailR_c2","",60,100.5,160.5); HfailRc2->Add(HtrRc2_D,HsuRc2_D,1.,-1.); TH1F* HfailLc2 = new TH1F("HfailL_c2","",60,100.5,160.5); HfailLc2->Add(HtrLc2_D,HsuLc2_D,1.,-1.); TH1F* HfailRc3 = new TH1F("HfailR_c3","",40,120.5,160.5); HfailRc3->Add(HtrRc3_D,HsuRc3_D,1.,-1.); TH1F* HfailLc3 = new TH1F("HfailL_c3","",40,120.5,160.5); HfailLc3->Add(HtrLc3_D,HsuLc3_D,1.,-1.); TH1F* HfailRc4 = new TH1F("HfailR_c4","",40,120.5,160.5); HfailRc4->Add(HtrRc4_D,HsuRc4_D,1.,-1.); TH1F* HfailLc4 = new TH1F("HfailL_c4","",40,120.5,160.5); HfailLc4->Add(HtrLc4_D,HsuLc4_D,1.,-1.); TH1F* EffT_in = new TH1F("EffT_in","",128,0.5,128.5); EffT_in->Divide(HsuTin_D,HtrTin_D,100.,1.); TH1F* EffB_in = new TH1F("EffB_in","",128,0.5,128.5); EffB_in->Divide(HsuBin_D,HtrBin_D,100.,1.); TH1F* EffR_in = new TH1F("EffR_in","",160,0.5,160.5); EffR_in->Divide(HsuRin_D,HtrRin_D,100.,1.); TH1F* EffL_in = new TH1F("EffL_in","",160,0.5,160.5); EffL_in->Divide(HsuLin_D,HtrLin_D,100.,1.); TH1F* EffT_ou = new TH1F("EffT_ou","",128,0.5,128.5); EffT_ou->Divide(HsuTou_D,HtrTou_D,100.,1.); TH1F* EffB_ou = new TH1F("EffB_ou","",128,0.5,128.5); EffB_ou->Divide(HsuBou_D,HtrBou_D,100.,1.); TH1F* EffR_ou = new TH1F("EffR_ou","",160,0.5,160.5); EffR_ou->Divide(HsuRou_D,HtrRou_D,100.,1.); TH1F* EffL_ou = new TH1F("EffL_ou","",160,0.5,160.5); EffL_ou->Divide(HsuLou_D,HtrLou_D,100.,1.); TH1F* EffT_c1 = new TH1F("EffT_c1","",128,0.5,128.5); EffT_c1->Divide(HsuTc1_D,HtrTc1_D,100.,1.); Double_t Tc1t = HtrTc1_D->Integral(); Double_t Tc1s = HsuTc1_D->Integral(); Double_t Tc1Avg = Tc1s/Tc1t; Double_t Tc1Err = sqrt((Tc1Avg*(1.0-Tc1Avg))/Tc1t); char xTc1[50]; sprintf(xTc1,"Mean = %.3f Error = %.3f", Tc1Avg, Tc1Err); TH1F* EffB_c1 = new TH1F("EffB_c1","",128,0.5,128.5); EffB_c1->Divide(HsuBc1_D,HtrBc1_D,100.,1.); Double_t Bc1t = HtrBc1_D->Integral(); Double_t Bc1s = HsuBc1_D->Integral(); Double_t Bc1Avg = Bc1s/Bc1t; Double_t Bc1Err = sqrt(Bc1Avg*(1-Bc1Avg)/Bc1t); char xBc1[50]; sprintf(xBc1,"Mean = %.3f Error = %.3f", Bc1Avg, Bc1Err); TH1F* EffR_c1 = new TH1F("EffR_c1","",160,0.5,160.5); EffR_c1->Divide(HsuRc1_D,HtrRc1_D,100.,1.); Double_t Rc1t = HtrRc1_D->Integral(); Double_t Rc1s = HsuRc1_D->Integral(); Double_t Rc1Avg = Rc1s/Rc1t; Double_t Rc1Err = sqrt(Rc1Avg*(1-Rc1Avg)/Rc1t); char xRc1[50]; sprintf(xRc1,"Mean = %.3f Error = %.3f", Rc1Avg, Rc1Err); TH1F* EffL_c1 = new TH1F("EffL_c1","",160,0.5,160.5); EffL_c1->Divide(HsuLc1_D,HtrLc1_D,100.,1.); Double_t Lc1t = HtrLc1_D->Integral(); Double_t Lc1s = HsuLc1_D->Integral(); Double_t Lc1Avg = Lc1s/Lc1t; Double_t Lc1Err = sqrt(Lc1Avg*(1-Lc1Avg)/Lc1t); char xLc1[50]; sprintf(xLc1,"Mean = %.3f Error = %.3f", Lc1Avg, Lc1Err); TH1F* EffR_c2 = new TH1F("EffR_c2","",60,100.5,160.5); EffR_c2->Divide(HsuRc2_D,HtrRc2_D,100.,1.); Double_t Rc2t = HtrRc2_D->Integral(); Double_t Rc2s = HsuRc2_D->Integral(); Double_t Rc2Avg = Rc2s/Rc2t; Double_t Rc2Err = sqrt(Rc2Avg*(1-Rc2Avg)/Rc2t); char xRc2[50]; sprintf(xRc2,"Mean = %.3f Error = %.3f", Rc2Avg, Rc2Err); TH1F* EffL_c2 = new TH1F("EffL_c2","",60,100.5,160.5); EffL_c2->Divide(HsuLc2_D,HtrLc2_D,100.,1.); Double_t Lc2t = HtrLc2_D->Integral(); Double_t Lc2s = HsuLc2_D->Integral(); Double_t Lc2Avg = Lc2s/Lc2t; Double_t Lc2Err = sqrt(Lc2Avg*(1-Lc2Avg)/Lc2t); char xLc2[50]; sprintf(xLc2,"Mean = %.3f Error = %.3f", Lc2Avg, Lc2Err); TH1F* EffR_c3 = new TH1F("EffR_c3","",40,120.5,160.5); EffR_c3->Divide(HsuRc3_D,HtrRc3_D,100.,1.); Double_t Rc3t = HtrRc3_D->Integral(); Double_t Rc3s = HsuRc3_D->Integral(); Double_t Rc3Avg = Rc3s/Rc3t; Double_t Rc3Err = sqrt(Rc3Avg*(1-Rc3Avg)/Rc3t); char xRc3[50]; sprintf(xRc3,"Mean = %.3f Error = %.3f", Rc3Avg, Rc3Err); TH1F* EffL_c3 = new TH1F("EffL_c3","",40,120.5,160.5); EffL_c3->Divide(HsuLc3_D,HtrLc3_D,100.,1.); Double_t Lc3t = HtrLc3_D->Integral(); Double_t Lc3s = HsuLc3_D->Integral(); Double_t Lc3Avg = Lc3s/Lc3t; Double_t Lc3Err = sqrt(Lc3Avg*(1-Lc3Avg)/Lc3t); char xLc3[50]; sprintf(xLc3,"Mean = %.3f Error = %.3f", Lc3Avg, Lc3Err); TH1F* EffR_c4 = new TH1F("EffR_c4","",40,120.5,160.5); EffR_c4->Divide(HsuRc4_D,HtrRc4_D,100.,1.); Double_t Rc4t = HtrRc4_D->Integral(); Double_t Rc4s = HsuRc4_D->Integral(); Double_t Rc4Avg = Rc4s/Rc4t; Double_t Rc4Err = sqrt(Rc4Avg*(1-Rc4Avg)/Rc4t); char xRc4[50]; sprintf(xRc4,"Mean = %.3f Error = %.3f", Rc4Avg, Rc4Err); TH1F* EffL_c4 = new TH1F("EffL_c4","",40,120.5,160.5); EffL_c4->Divide(HsuLc4_D,HtrLc4_D,100.,1.); Double_t Lc4t = HtrLc4_D->Integral(); Double_t Lc4s = HsuLc4_D->Integral(); Double_t Lc4Avg = Lc4s/Lc4t; Double_t Lc4Err = sqrt(Lc4Avg*(1-Lc4Avg)/Lc4t); char xLc4[50]; sprintf(xLc4,"Mean = %.3f Error = %.3f", Lc4Avg, Lc4Err); TH1F* Effcen_TB = new TH1F("Effcen_TB","",51,0.5,1.01); Double_t bin = 0; for(Int_t i0=1;i0<129;i0++) { bin = (EffT_c1->GetBinContent(i0))/100.; if(bin>0.) Effcen_TB->Fill(bin); bin = (EffB_c1->GetBinContent(i0))/100.; if(bin>0.) Effcen_TB->Fill(bin); } TH1F* Effcen_RL = new TH1F("Effcen_RL","",51,0.5,1.01); bin = 0; for(Int_t i1=1;i1<161;i1++) { bin = (EffR_c1->GetBinContent(i1))/100.; if(bin>0.) Effcen_RL->Fill(bin); bin = (EffL_c1->GetBinContent(i1))/100.; if(bin>0.) Effcen_RL->Fill(bin); } for(Int_t i2=1;i2<61;i2++) { bin = (EffR_c2->GetBinContent(i2))/100.; if(bin>0.) Effcen_RL->Fill(bin); bin = (EffL_c2->GetBinContent(i2))/100.; if(bin>0.) Effcen_RL->Fill(bin); } for(Int_t i3=1;i3<41;i3++) { bin = (EffR_c3->GetBinContent(i3))/100.; if(bin>0.) Effcen_RL->Fill(bin); bin = (EffL_c3->GetBinContent(i3))/100.; if(bin>0.) Effcen_RL->Fill(bin); } for(Int_t i4=1;i4<41;i4++) { bin = (EffR_c4->GetBinContent(i4))/100.; if(bin>0.) Effcen_RL->Fill(bin); bin = (EffL_c4->GetBinContent(i4))/100.; if(bin>0.) Effcen_RL->Fill(bin); } // Efficiency Summary Effcen_TB->SetLineColor(1); sprintf(buff,"Central efficiency TB"); Effcen_TB->SetXTitle(buff); Effcen_RL->SetLineColor(1); sprintf(buff,"Central efficiency RL"); Effcen_RL->SetXTitle(buff); // Trial,Success,Failure,Efficiency Profiles inner HtrTin_D->SetLineColor(1); HtrTin_D->SetXTitle("Trial Profile Data Top-inner"); HtrBin_D->SetLineColor(1); HtrBin_D->SetXTitle("Trial Profile Data Bottom-inner"); HtrRin_D->SetLineColor(1); HtrRin_D->SetXTitle("Trial Profile Data Right-inner"); HtrLin_D->SetLineColor(1); HtrLin_D->SetXTitle("Trial Profile Data Left-inner"); HsuTin_D->SetLineColor(1); HsuTin_D->SetXTitle("Success Profile Data Top-inner"); HsuBin_D->SetLineColor(1); HsuBin_D->SetXTitle("Success Profile Data Bottom-inner"); HsuRin_D->SetLineColor(1); HsuRin_D->SetXTitle("Success Profile Data Right-inner"); HsuLin_D->SetLineColor(1); HsuLin_D->SetXTitle("Success Profile Data Left-inner"); HfailTin->SetLineColor(1); HfailTin->SetXTitle("Failure Profile Data Top-inner"); HfailBin->SetLineColor(1); HfailBin->SetXTitle("Failure Profile Data Bottom-inner"); HfailRin->SetLineColor(1); HfailRin->SetXTitle("Failure Profile Data Right-inner"); HfailLin->SetLineColor(1); HfailLin->SetXTitle("Failure Profile Data Left-inner"); EffT_in->SetLineColor(1); EffT_in->SetXTitle("Efficiency Profile Data Top-inner"); EffB_in->SetLineColor(1); EffB_in->SetXTitle("Efficiency Profile Data Bottom-inner"); EffR_in->SetLineColor(1); EffR_in->SetXTitle("Efficiency Profile Data Right-inner"); EffL_in->SetLineColor(1); EffL_in->SetXTitle("Efficiency Profile Data Left-inner"); // Trial,Success,Failure,Efficiency Profiles outer HtrTou_D->SetLineColor(1); HtrTou_D->SetXTitle("Trial Profile Data Top-outer"); HtrBou_D->SetLineColor(1); HtrBou_D->SetXTitle("Trial Profile Data Bottom-outer"); HtrRou_D->SetLineColor(1); HtrRou_D->SetXTitle("Trial Profile Data Right-outer"); HtrLou_D->SetLineColor(1); HtrLou_D->SetXTitle("Trial Profile Data Left-outer"); HsuTou_D->SetLineColor(1); HsuTou_D->SetXTitle("Success Profile Data Top-outer"); HsuBou_D->SetLineColor(1); HsuBou_D->SetXTitle("Success Profile Data Bottom-outer"); HsuRou_D->SetLineColor(1); HsuRou_D->SetXTitle("Success Profile Data Right-outer"); HsuLou_D->SetLineColor(1); HsuLou_D->SetXTitle("Success Profile Data Left-outer"); HfailTou->SetLineColor(1); HfailTou->SetXTitle("Failure Profile Data Top-outer"); HfailBou->SetLineColor(1); HfailBou->SetXTitle("Failure Profile Data Bottom-outer"); HfailRou->SetLineColor(1); HfailRou->SetXTitle("Failure Profile Data Right-outer"); HfailLou->SetLineColor(1); HfailLou->SetXTitle("Failure Profile Data Left-outer"); EffT_ou->SetLineColor(1); EffT_ou->SetXTitle("Efficiency Profile Data Top-outer"); EffB_ou->SetLineColor(1); EffB_ou->SetXTitle("Efficiency Profile Data Bottom-outer"); EffR_ou->SetLineColor(1); EffR_ou->SetXTitle("Efficiency Profile Data Right-outer"); EffL_ou->SetLineColor(1); EffL_ou->SetXTitle("Efficiency Profile Data Left-outer"); // Trial,Success,Efficiency Profiles center HtrTc1_D->SetLineColor(1); HtrTc1_D->SetXTitle("Trial Profile Data Top-center_1"); HtrBc1_D->SetLineColor(1); HtrBc1_D->SetXTitle("Trial Profile Data Bottom-center_1"); HtrRc1_D->SetLineColor(1); HtrRc1_D->SetXTitle("Trial Profile Data Right-center_1"); HtrLc1_D->SetLineColor(1); HtrLc1_D->SetXTitle("Trial Profile Data Left-center_1"); HsuTc1_D->SetLineColor(1); HsuTc1_D->SetXTitle("Success Profile Data Top-center_1"); HsuBc1_D->SetLineColor(1); HsuBc1_D->SetXTitle("Success Profile Data Bottom-center_1"); HsuRc1_D->SetLineColor(1); HsuRc1_D->SetXTitle("Success Profile Data Right-center_1"); HsuLc1_D->SetLineColor(1); HsuLc1_D->SetXTitle("Success Profile Data Left-center_1"); HfailTc1->SetLineColor(1); HfailTc1->SetXTitle("Failure Profile Data Top-center_1"); HfailBc1->SetLineColor(1); HfailBc1->SetXTitle("Failure Profile Data Bottom-center_1"); HfailRc1->SetLineColor(1); HfailRc1->SetXTitle("Failure Profile Data Right-center_1"); HfailLc1->SetLineColor(1); HfailLc1->SetXTitle("Failure Profile Data Left-center_1"); EffT_c1->SetLineColor(1); EffT_c1->SetXTitle("Efficiency Profile Data Top-center_1"); EffB_c1->SetLineColor(1); EffB_c1->SetXTitle("Efficiency Profile Data Bottom-center_1"); EffR_c1->SetLineColor(1); EffR_c1->SetXTitle("Efficiency Profile Data Right-center_1"); EffL_c1->SetLineColor(1); EffL_c1->SetXTitle("Efficiency Profile Data Left-center_1"); HtrRc2_D->SetLineColor(1); HtrRc2_D->SetXTitle("Trial Profile Data Right-center_2"); HtrLc2_D->SetLineColor(1); HtrLc2_D->SetXTitle("Trial Profile Data Left-center_2"); HsuRc2_D->SetLineColor(1); HsuRc2_D->SetXTitle("Success Profile Data Right-center_2"); HsuLc2_D->SetLineColor(1); HsuLc2_D->SetXTitle("Success Profile Data Left-center_2"); HfailRc2->SetLineColor(1); HfailRc2->SetXTitle("Failure Profile Data Right-center_2"); HfailLc2->SetLineColor(1); HfailLc2->SetXTitle("Failure Profile Data Left-center_2"); EffR_c2->SetLineColor(1); EffR_c2->SetXTitle("Efficiency Profile Data Right-center_2"); EffL_c2->SetLineColor(1); EffL_c2->SetXTitle("Efficiency Profile Data Left-center_2"); HtrRc3_D->SetLineColor(1); HtrRc3_D->SetXTitle("Trial Profile Data Right-center_3"); HtrLc3_D->SetLineColor(1); HtrLc3_D->SetXTitle("Trial Profile Data Left-center_3"); HsuRc3_D->SetLineColor(1); HsuRc3_D->SetXTitle("Success Profile Data Right-center_3"); HsuLc3_D->SetLineColor(1); HsuLc3_D->SetXTitle("Success Profile Data Left-center_3"); HfailRc3->SetLineColor(1); HfailRc3->SetXTitle("Failure Profile Data Right-center_3"); HfailLc3->SetLineColor(1); HfailLc3->SetXTitle("Failure Profile Data Left-center_3"); EffR_c3->SetLineColor(1); EffR_c3->SetXTitle("Efficiency Profile Data Right-center_3"); EffL_c3->SetLineColor(1); EffL_c3->SetXTitle("Efficiency Profile Data Left-center_3"); HtrRc4_D->SetLineColor(1); HtrRc4_D->SetXTitle("Trial Profile Data Right-center_4"); HtrLc4_D->SetLineColor(1); HtrLc4_D->SetXTitle("Trial Profile Data Left-center_4"); HsuRc4_D->SetLineColor(1); HsuRc4_D->SetXTitle("Success Profile Data Right-center_4"); HsuLc4_D->SetLineColor(1); HsuLc4_D->SetXTitle("Success Profile Data Left-center_4"); HfailRc4->SetLineColor(1); HfailRc4->SetXTitle("Failure Profile Data Right-center_4"); HfailLc4->SetLineColor(1); HfailLc4->SetXTitle("Failure Profile Data Left-center_4"); EffR_c4->SetLineColor(1); EffR_c4->SetXTitle("Efficiency Profile Data Right-center_4"); EffL_c4->SetLineColor(1); EffL_c4->SetXTitle("Efficiency Profile Data Left-center_4"); //Non-efficiency plots // CQ_fail_D->SetLineColor(1); // CQ_fail_D->SetXTitle("CQ_fail"); // CQ_good_D->SetLineColor(1); // CQ_good_D->SetXTitle("CQ_good"); // TQ_fail_D->SetLineColor(1); // TQ_fail_D->SetXTitle("TQ_fail"); // TQ_good_D->SetLineColor(1); // TQ_good_D->SetXTitle("TQ_good"); // Cmb_Q_D->SetLineColor(1); // Cmb_Q_D->SetXTitle("Combo Hit Charge"); // Cmb_Qwin_D->SetLineColor(1); // Cmb_Qwin_D->SetXTitle("Combo Hit Charge (window)"); // Cmb_T_D->SetLineColor(1); // Cmb_T_D->SetXTitle("Combo Hit Time"); // // LgrN_D->SetLineColor(1); // // LgrN_D->SetXTitle("LgrN"); // Sng_Q_D->SetLineColor(1); // Sng_Q_D->SetXTitle("Single Hit Charge"); // Sng_T_D->SetLineColor(1); // Sng_T_D->SetXTitle("Single Hit Time"); // // Sng_T_fine_D->SetLineColor(1); // // Sng_T_fine_D->SetXTitle("Single Hit Time-Fine"); // // TgrN_D->SetLineColor(1); // // TgrN_D->SetXTitle("TgrN"); // cmbN_D->SetLineColor(1); // cmbN_D->SetXTitle("Number of Combo Hits"); // sngN_D->SetLineColor(1); // sngN_D->SetXTitle("Number of Single Hits"); // HitR_Fail_D->SetLineColor(1); // HitR_Fail_D->SetXTitle("Number of Hits in a Slab for a Failed Trial"); // HitT1_Fail_D->SetLineColor(1); // HitT1_Fail_D->SetXTitle("Hit Time Difference for Failed Trials (qmin < 12)"); // HitT2_Fail_D->SetLineColor(1); // HitT2_Fail_D->SetXTitle("Hit Time Difference for Failed Trials (qmin 12-30)"); // HitT3_Fail_D->SetLineColor(1); // HitT3_Fail_D->SetXTitle("Hit Time Difference for Failed Trials (qmin > 30)"); // TrgQt_D->SetLineColor(1); // TrgQt_D->SetXTitle("Total Charge of Slabs Surrounding Trial Slab"); // TrgWin_D->SetLineColor(1); // TrgWin_D->SetXTitle("Time Difference of Hits between Slabs Surrounding Trial Slab (in ns)"); // tEnd_D->SetLineColor(1); // tEnd_D->SetXTitle("Time between MPPC firings in a slab (in ns)"); // qEnd_D->SetLineColor(1); // qEnd_D->SetXTitle("Hit Charges"); // TrgWinQ_D->SetLineColor(1); // TrgWinQ_D->SetXTitle("Time Difference of Hits between Slabs Surrounding Trial Slab (in ns) After Charge Cut"); // HitWin_D->SetLineColor(1); // HitWin_D->SetXTitle("Time Difference of success and average hit time of surrounding slabs (in ns)"); // // Recon_TvsQ_D->SetLineColor(1); // // Recon_TvsQ_D->SetXTitle("Hit Time Difference vs Sum of Hit Charges"); // // Recon_TvsQmin_D->SetLineColor(1); // // Recon_TvsQmin_D->SetXTitle("Hit Time Difference vs Minimum of Hit Charges"); sprintf(buff,"%s",argv[2]); TCanvas *c10 = new TCanvas("c10","",0,0,612,927); c10->Divide(1,4); c10->cd(1); HtrTin_D->Draw("e1"); c10->cd(2); HsuTin_D->Draw("e1"); c10->cd(3); HfailTin->Draw("e1"); c10->cd(4); EffT_in->Draw("e1"); c10->Print(buff,"Portrait"); TCanvas *c12 = new TCanvas("c12","",0,0,612,927); c12->Divide(1,4); c12->cd(1); HtrBin_D->Draw("e1"); c12->cd(2); HsuBin_D->Draw("e1"); c12->cd(3); HfailBin->Draw("e1"); c12->cd(4); EffB_in->Draw("e1"); c12->Print(buff,"Portrait"); TCanvas *c14 = new TCanvas("c14","",0,0,612,927); c14->Divide(1,4); c14->cd(1); HtrRin_D->Draw("e1"); c14->cd(2); HsuRin_D->Draw("e1"); c14->cd(3); HfailRin->Draw("e1"); c14->cd(4); EffR_in->Draw("e1"); c14->Print(buff,"Portrait"); TCanvas *c16 = new TCanvas("c16","",0,0,612,927); c16->Divide(1,4); c16->cd(1); HtrLin_D->Draw("e1"); c16->cd(2); HsuLin_D->Draw("e1"); c16->cd(3); HfailLin->Draw("e1"); c16->cd(4); EffL_in->Draw("e1"); c16->Print(buff,"Portrait"); TCanvas *c18 = new TCanvas("c18","",0,0,612,927); c18->Divide(1,4); c18->cd(1); HtrTou_D->Draw("e1"); c18->cd(2); HsuTou_D->Draw("e1"); c18->cd(3); HfailTou->Draw("e1"); c18->cd(4); EffT_ou->Draw("e1"); c18->Print(buff,"Portrait"); TCanvas *c20 = new TCanvas("c20","",0,0,612,927); c20->Divide(1,4); c20->cd(1); HtrBou_D->Draw("e1"); c20->cd(2); HsuBou_D->Draw("e1"); c20->cd(3); HfailBou->Draw("e1"); c20->cd(4); EffB_ou->Draw("e1"); c20->Print(buff,"Portrait"); TCanvas *c22 = new TCanvas("c22","",0,0,612,927); c22->Divide(1,4); c22->cd(1); HtrRou_D->Draw("e1"); c22->cd(2); HsuRou_D->Draw("e1"); c22->cd(3); HfailRou->Draw("e1"); c22->cd(4); EffR_ou->Draw("e1"); c22->Print(buff,"Portrait"); TCanvas *c24 = new TCanvas("c24","",0,0,612,927); c24->Divide(1,4); c24->cd(1); HtrLou_D->Draw("e1"); c24->cd(2); HsuLou_D->Draw("e1"); c24->cd(3); HfailLou->Draw("e1"); c24->cd(4); EffL_ou->Draw("e1"); c24->Print(buff,"Portrait"); TCanvas *c26 = new TCanvas("c26","",0,0,612,927); c26->Divide(1,4); c26->cd(1); HtrTc1_D->Draw("e1"); c26->cd(2); HsuTc1_D->Draw("e1"); c26->cd(3); HfailTc1->Draw("e1"); c26->cd(4); EffT_c1->Draw("e1"); TPaveText *Tc1 = new TPaveText(.2,.4,.55,.68,"NDC"); Tc1->SetTextSize(0.08); Tc1->SetFillColor(0); Tc1->SetTextAlign(12); /*TText* Tc1Entry =*/ Tc1->AddText(xTc1); Tc1->Draw("e1"); sprintf(buff,"%s(",argv[2]); c26->Print(buff,"Portrait"); TCanvas *c28 = new TCanvas("c28","",0,0,612,927); c28->Divide(1,4); c28->cd(1); HtrBc1_D->Draw("e1"); c28->cd(2); HsuBc1_D->Draw("e1"); c28->cd(3); HfailBc1->Draw("e1"); c28->cd(4); EffB_c1->Draw("e1"); TPaveText *Bc1 = new TPaveText(.2,.4,.55,.68,"NDC"); Bc1->SetTextSize(0.08); Bc1->SetFillColor(0); Bc1->SetTextAlign(12); /*TText* Bc1Entry =*/ Bc1->AddText(xBc1); Bc1->Draw("e1"); sprintf(buff,"%s(",argv[2]); c28->Print(buff,"Portrait"); TCanvas *c30 = new TCanvas("c30","",0,0,612,927); c30->Divide(1,4); c30->cd(1); HtrRc1_D->Draw("e1"); c30->cd(2); HsuRc1_D->Draw("e1"); c30->cd(3); HfailRc1->Draw("e1"); c30->cd(4); EffR_c1->Draw("e1"); TPaveText *Rc1 = new TPaveText(.2,.4,.55,.68,"NDC"); Rc1->SetTextSize(0.08); Rc1->SetFillColor(0); Rc1->SetTextAlign(12); /*TText* Rc1Entry =*/ Rc1->AddText(xRc1); Rc1->Draw("e1"); c30->Print(buff,"Portrait"); TCanvas *c32 = new TCanvas("c32","",0,0,612,927); c32->Divide(1,4); c32->cd(1); HtrLc1_D->Draw("e1"); c32->cd(2); HsuLc1_D->Draw("e1"); c32->cd(3); HfailLc1->Draw("e1"); c32->cd(4); EffL_c1->Draw("e1"); TPaveText *Lc1 = new TPaveText(.2,.4,.55,.68,"NDC"); Lc1->SetTextSize(0.08); Lc1->SetFillColor(0); Lc1->SetTextAlign(12); /*TText* Lc1Entry =*/ Lc1->AddText(xLc1); Lc1->Draw("e1"); c32->Print(buff,"Portrait"); TCanvas *c34 = new TCanvas("c34","",0,0,612,927); c34->Divide(1,4); c34->cd(1); HtrRc2_D->Draw("e1"); c34->cd(2); HsuRc2_D->Draw("e1"); c34->cd(3); HfailRc2->Draw("e1"); c34->cd(4); EffR_c2->Draw("e1"); TPaveText *Rc2 = new TPaveText(.2,.4,.55,.68,"NDC"); Rc2->SetTextSize(0.08); Rc2->SetFillColor(0); Rc2->SetTextAlign(12); /*TText* Rc2Entry =*/ Rc2->AddText(xRc2); Rc2->Draw("e1"); c34->Print(buff,"Portrait"); TCanvas *c36 = new TCanvas("c36","",0,0,612,927); c36->Divide(1,4); c36->cd(1); HtrLc2_D->Draw("e1"); c36->cd(2); HsuLc2_D->Draw("e1"); c36->cd(3); HfailLc2->Draw("e1"); c36->cd(4); EffL_c2->Draw("e1"); TPaveText *Lc2 = new TPaveText(.2,.4,.55,.68,"NDC"); Lc2->SetTextSize(0.08); Lc2->SetFillColor(0); Lc2->SetTextAlign(12); /*TText* Lc2Entry =*/ Lc2->AddText(xLc2); Lc2->Draw("e1"); c36->Print(buff,"Portrait"); TCanvas *c38 = new TCanvas("c38","",0,0,612,927); c38->Divide(1,4); c38->cd(1); HtrRc3_D->Draw("e1"); c38->cd(2); HsuRc3_D->Draw("e1"); c38->cd(3); HfailRc3->Draw("e1"); c38->cd(4); EffR_c3->Draw("e1"); TPaveText *Rc3 = new TPaveText(.2,.4,.55,.68,"NDC"); Rc3->SetTextSize(0.08); Rc3->SetFillColor(0); Rc3->SetTextAlign(12); /*TText* Rc3Entry =*/ Rc3->AddText(xRc3); Rc3->Draw("e1"); c38->Print(buff,"Portrait"); TCanvas *c40 = new TCanvas("c40","",0,0,612,927); c40->Divide(1,4); c40->cd(1); HtrLc3_D->Draw("e1"); c40->cd(2); HsuLc3_D->Draw("e1"); c40->cd(3); HfailLc3->Draw("e1"); c40->cd(4); EffL_c3->Draw("e1"); TPaveText *Lc3 = new TPaveText(.2,.4,.55,.68,"NDC"); Lc3->SetTextSize(0.08); Lc3->SetFillColor(0); Lc3->SetTextAlign(12); /*TText* Lc3Entry =*/ Lc3->AddText(xLc3); Lc3->Draw("e1"); c40->Print(buff,"Portrait"); TCanvas *c42 = new TCanvas("c42","",0,0,612,927); c42->Divide(1,4); c42->cd(1); HtrRc4_D->Draw("e1"); c42->cd(2); HsuRc4_D->Draw("e1"); c42->cd(3); HfailRc4->Draw("e1"); c42->cd(4); EffR_c4->Draw("e1"); TPaveText *Rc4 = new TPaveText(.2,.4,.55,.68,"NDC"); Rc4->SetTextSize(0.08); Rc4->SetFillColor(0); Rc4->SetTextAlign(12); /*TText* Rc4Entry =*/ Rc4->AddText(xRc4); Rc4->Draw("e1"); c42->Print(buff,"Portrait"); TCanvas *c44 = new TCanvas("c44","",0,0,612,927); c44->Divide(1,4); c44->cd(1); HtrLc4_D->Draw("e1"); c44->cd(2); HsuLc4_D->Draw("e1"); c44->cd(3); HfailLc4->Draw("e1"); c44->cd(4); EffL_c4->Draw("e1"); TPaveText *Lc4 = new TPaveText(.2,.4,.55,.68,"NDC"); Lc4->SetTextSize(0.08); Lc4->SetFillColor(0); Lc4->SetTextAlign(12); /*TText* Lc4Entry =*/ Lc4->AddText(xLc4); Lc4->Draw("e1"); c44->Print(buff,"Portrait"); // TCanvas *c46 = new TCanvas("c46","",0,0,612,927); // c46->Divide(1,4); // c46->cd(1); // TrgQt_D->Draw("e1"); // c46->cd(2); // TrgWin_D->Draw("e1"); // c46->cd(3); // TrgWinQ_D->Draw("e1"); // c46->cd(4); // HitWin_D->Draw("e1"); // c46->Print(buff,"Portrait"); // TCanvas *c48 = new TCanvas("c48","",0,0,612,927); // c48->Divide(1,3); // c48->cd(1); // tEnd_D->Draw("e1"); // c48->cd(2); // qEnd_D->Draw("e1"); // c48->cd(3); // HitR_Fail_D->Draw("e1"); // c48->Print(buff,"Portrait"); // TCanvas *c50 = new TCanvas("c50","",0,0,612,927); // c50->Divide(1,3); // c50->cd(1); // HitT1_Fail_D->Draw("e1"); // c50->cd(2); // HitT2_Fail_D->Draw("e1"); // c50->cd(3); // HitT3_Fail_D->Draw("e1"); // c50->Print(buff,"Portrait"); TCanvas *c99 = new TCanvas("c99","",0,0,612,927); c99->Divide(1,2); c99->cd(1); Effcen_TB->Draw("e1"); c99->cd(2); Effcen_RL->Draw("e1"); sprintf(buff,"%s(",argv[2]); c99->Print(buff,"Portrait"); f1.Close(); }