00001
00002
00003 #include <iostream>
00004 #include <string>
00005 #include <typeinfo>
00006 #include <iostream>
00007 #include <vector>
00008
00009
00010 #include "Riostream.h"
00011 #include "TRint.h"
00012 #include "TROOT.h"
00013 #include "TStyle.h"
00014 #include "Riostream.h"
00015
00016
00017 #include "TH2.h"
00018 #include "TH1.h"
00019 #include "TF1.h"
00020
00021 #include "TGraphErrors.h"
00022 #include "TMultiGraph.h"
00023 #include "TCanvas.h"
00024 #include "TProfile.h"
00025 #include "TStopwatch.h"
00026
00027 #include "TFile.h"
00028 #include "TTree.h"
00029 #include "TBranch.h"
00030 #include "TMath.h"
00031 #include "TLegend.h"
00032 #include "TColor.h"
00033 #include "TText.h"
00034 #include "TLine.h"
00035 #include "TKey.h"
00036 #include "TLatex.h"
00037 #include "TCanvas.h"
00038
00039
00040
00041
00042
00043 #define NUM_SPADIC_CHA 8
00044 #define SPADIC_TRACE_SIZE 45
00045 #define SPADIC_ADC_MAX 255
00046 #define noBaselineTB 5
00047
00048 using namespace std;
00049
00050
00051 void plotResults() {
00052
00053 gROOT->SetStyle("Plain");
00054 gStyle->SetOptStat(kFALSE);
00055 gStyle->SetOptTitle(kFALSE);
00056 gStyle->SetPalette(1,0);
00057
00058 gStyle->SetPadTickX(1);
00059 gStyle->SetPadTickY(1);
00060 const TString picExt = ".eps";
00061 const Bool_t dummy = true;
00062 TString ofolder = "data2011/Analysis/Pics/Results/";
00063
00064 TString ifolder = "data2011/Analysis/";
00065 TString filename[24] = {
00066 "MS_HVScan_1.root",
00067 "MS_HVScan_2.root",
00068 "MS_HVScan_3.root",
00069 "MS_HVScan_4.root",
00070 "MS_HVScan_5.root",
00071 "MS_MomentumScan_2GeV.root",
00072 "MS_MomentumScan_4GeV.root",
00073 "MS_MomentumScan_6GeV.root",
00074 "MS_MomentumScan_8GeV.root",
00075 "MS_MomentumScan_10GeV.root",
00076 "MS_parallelBeam.root",
00077 "MS_RadiatorScan_1.root",
00078 "MS_RadiatorScan_2.root",
00079 "MS_RadiatorScan_3.root",
00080 "MS_RadiatorScan_4.root",
00081 "MS_RadiatorScan_5.root",
00082 "MS_RadiatorScan_6.root",
00083 "MS_RadiatorScan_7.root",
00084 "MS_RadiatorScan_8.root",
00085 "MS_Rotated.root",
00086 "MS_Tilted.root",
00087 "MS_BlackFoam.root",
00088 "MS_WhiteFoam.root",
00089 "MS_WhiteFoam_2GeV.root"
00090 };
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100 TString mergeName[2] = {"NotMerged","Merged"};
00101 TString filterName[4] = {"wOwU","woOwU","woOwoU","wOwoU"};
00102 Int_t usedSuId[8] = { 2, 12, 6, 18, 17, 16, 4, 5 };
00103 TString usedSuName[8] = {"MS336/10 ","MS336/11 ","MS444/10 ","MS444/11 "," "," "," "," "};
00104 TString recoSetName[9] = {
00105 "ClusterInteg",
00106 "3PadInteg",
00107 "1PadInteg",
00108 "ClusterIntegWindow",
00109 "3PadIntegWindow",
00110 "1PadIntegWindow",
00111 "ClusterAmpl",
00112 "3PadAmpl",
00113 "1PadAmpl"
00114 };
00115 TString macroPre[3] = {"EP_","PR_",""};
00116 TCanvas* c1_1 = new TCanvas("c1_1","c1_1", 800, 600);
00117 c1_1->Divide(1,1);
00118 TCanvas* c1_2 = new TCanvas("c1_2","c1_2", 800,1200);
00119 c1_2->Divide(1,2);
00120 TCanvas* c2_1 = new TCanvas("c2_1","c2_1",1600, 600);
00121 c2_1->Divide(2,1);
00122 TCanvas* c2_2 = new TCanvas("c2_2","c2_2",1600,1200);
00123 c2_2->Divide(2,2);
00124 TString title, newtitle, appending;
00125 TProfile* hEfficiency[40];
00126 TProfile* hPR[40];
00127
00128 TLegend* leg[4];
00129 leg[0] = new TLegend(0.12,0.15,0.88,0.4);
00130 leg[1] = new TLegend(0.12,0.15,0.88,0.4);
00131 leg[2] = new TLegend(0.12,0.15,0.88,0.4);
00132 leg[3] = new TLegend(0.12,0.15,0.88,0.4);
00133 for (Int_t i = 0; i < 4; i++) {
00134
00135 leg[i]->SetLineColor(0);
00136 leg[i]->SetLineStyle(0);
00137 leg[i]->SetLineWidth(0);
00138 leg[i]->SetFillColor(0);
00139 leg[i]->SetFillStyle(0);
00140 leg[i]->SetTextSize(0.03);
00141 }
00142
00143 Int_t suid(0), rs(1), imerge(0), ifilter(0), irebin(0), macro(0);
00144 for (imerge = 0; imerge < 2; imerge++) {
00145
00146 printf("::::::::::::: Anode Scan :::::::::::::\n");
00147 TString AHV[3][4] = {{"1800V","1800V","2000V","2000V"},
00148 {"1750V","1750V","1950V","1950V"},
00149 {"1700V","1700V","1900V","1900V"}
00150 };
00151 TString DV[4] = {"/600V","/600V","/400V","/400V"};
00152
00153 TString AV[4] = {"1750/","1750/","1950/","1950/"};
00154 TString DHV[3][4] = {{"600V","600V","400V","400V"},
00155 {"300V","300V","200V","200V"},
00156 { "60V", "60V", "40V", "40V"}
00157 };
00158
00159 if(dummy) {
00160 Double_t ax[4][3] = {{1800,1750,1700},
00161 {1800,1750,1700},
00162 {2000,1950,1900},
00163 {2000,1950,1900}
00164 };
00165 Double_t eax[4][3] = {{0.0}};
00166 Double_t ay[4][3] = {{0.0}};
00167 Double_t eay[4][3] = {{0.0}};
00168 Double_t dx[4][3] = {{600,300, 60},
00169 {600,300, 60},
00170 {400,200, 40},
00171 {400,200, 40}
00172 };
00173 Double_t edx[4][3] = {{0.0}};
00174 Double_t dy[4][3] = {{0.0}};
00175 Double_t edy[4][3] = {{0.0}};
00176 TProfile *h1dim[2][4];
00177 TGraphErrors* gAHV[4];
00178 TGraphErrors* gDHV[4];
00179 for (suid = 0; suid < 4; suid++) {
00180 title.Form("Anode_%i",usedSuId[suid]);
00181 h1dim[0][suid] = new TProfile(title,usedSuName[suid]+"Anode",7,1675,2025);
00182 h1dim[0][suid]->GetYaxis()->SetRangeUser(1e-3,100);
00183 h1dim[0][suid]->SetLineStyle(suid/2+1);
00184 h1dim[0][suid]->SetMarkerStyle(20);
00185 h1dim[0][suid]->SetMarkerColor(1+suid);
00186 h1dim[0][suid]->SetLineColor(1+suid);
00187 h1dim[0][suid]->SetXTitle("anode voltage [V]");
00188 h1dim[0][suid]->SetYTitle("pion efficiency [%]");
00189 title.Form("Drift_%i",usedSuId[suid]);
00190 h1dim[1][suid] = new TProfile(title,usedSuName[suid]+"Drift",65,0,650);
00191 h1dim[1][suid]->GetYaxis()->SetRangeUser(1e-3,100);
00192 h1dim[1][suid]->SetLineStyle(suid/2+1);
00193 h1dim[1][suid]->SetMarkerStyle(20);
00194 h1dim[1][suid]->SetMarkerColor(1+suid);
00195 h1dim[1][suid]->SetLineColor(1+suid);
00196 h1dim[1][suid]->SetXTitle("drift voltage [V]");
00197 h1dim[1][suid]->SetYTitle("pion efficiency [%]");
00198 }
00199
00200
00201 leg[0]->SetNColumns(4);
00202 leg[1]->SetNColumns(4);
00203 TString radi[4] = {" H"," E"," F"," I(p)"};
00204 TString header;
00205 for (Int_t i = 0; i < 4; i++)
00206 header += usedSuName[i] + radi[i];
00207
00208 leg[0]->SetHeader(header + ", 3GeV: ");
00209 leg[1]->SetHeader(header + ", 3GeV: ");
00210
00211 TFile inputFile1(ifolder + macroPre[macro] + filename[0],"READ");
00212 for (suid = 0; suid < 4; suid++) {
00213 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00214 newtitle.Form("Efficiency" + macroPre[macro] + filename[0] + "_%i",suid);
00215 if (hEfficiency[suid+0] = (TProfile*)inputFile1.Get(title)){
00216 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
00217 hEfficiency[suid+0]->SetMarkerStyle(20);
00218 hEfficiency[suid+0]->SetMarkerColor(1+suid);
00219 hEfficiency[suid+0]->SetLineColor(1+suid);
00220 ay[suid][0] = hEfficiency[suid+0]->GetBinContent(hEfficiency[suid+0]->FindBin(10));
00221
00222
00223
00224
00225
00226
00227 leg[0]->AddEntry(hEfficiency[suid+0],AHV[0][suid] + DV[suid],"lep");
00228 c2_1->cd(1)->SetLogy(1);
00229 if (suid == 0){
00230 hEfficiency[suid+0]->DrawCopy("L,P");
00231 }
00232 else{
00233 hEfficiency[suid+0]->DrawCopy("L,P,same");
00234 }
00235 }
00236 }
00237
00238 TFile inputFile2(ifolder + macroPre[macro] + filename[2],"READ");
00239 for (suid = 0; suid < 4; suid++) {
00240 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00241 newtitle.Form("Efficiency" + macroPre[macro] + filename[2] + "_%i",suid);
00242 if (hEfficiency[suid+4] = (TProfile*)inputFile2.Get(title)){
00243 hEfficiency[suid+4]->SetLineStyle(suid/2+1);
00244 hEfficiency[suid+4]->SetMarkerStyle(21);
00245 hEfficiency[suid+4]->SetMarkerColor(1+suid);
00246 hEfficiency[suid+4]->SetLineColor(1+suid);
00247 ay[suid][1] = hEfficiency[suid+4]->GetBinContent(hEfficiency[suid+4]->FindBin(10));
00248 dy[suid][0] = hEfficiency[suid+4]->GetBinContent(hEfficiency[suid+4]->FindBin(10));
00249
00250
00251
00252
00253
00254
00255 leg[0]->AddEntry(hEfficiency[suid+4],AHV[1][suid] + DV[suid],"lep");
00256 c2_1->cd(1)->SetLogy(1);
00257 hEfficiency[suid+4]->DrawCopy("L,P,same");
00258
00259
00260
00261
00262
00263
00264 leg[1]->AddEntry(hEfficiency[suid+0],AV[suid] + DHV[0][suid],"lep");
00265 c2_1->cd(2)->SetLogy(1);
00266 if (suid == 0)
00267 hEfficiency[suid+4]->DrawCopy("L,P");
00268 else
00269 hEfficiency[suid+4]->DrawCopy("L,P,same");
00270 }
00271 }
00272
00273 TFile inputFile3(ifolder + macroPre[macro] + filename[1],"READ");
00274 for (suid = 0; suid < 4; suid++) {
00275 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00276 newtitle.Form("Efficiency" + macroPre[macro] + filename[1] + "_%i",suid);
00277 if (hEfficiency[suid+8] = (TProfile*)inputFile3.Get(title)) {
00278 hEfficiency[suid+8]->SetLineStyle(suid/2+1);
00279 hEfficiency[suid+8]->SetMarkerStyle(22);
00280 hEfficiency[suid+8]->SetMarkerColor(1+suid);
00281 hEfficiency[suid+8]->SetLineColor(1+suid);
00282 ay[suid][2] = hEfficiency[suid+8]->GetBinContent(hEfficiency[suid+8]->FindBin(10));
00283
00284
00285
00286
00287
00288
00289 leg[0]->AddEntry(hEfficiency[suid+8],AHV[2][suid] + DV[suid],"lep");
00290 c2_1->cd(1)->SetLogy(1);
00291 hEfficiency[suid+8]->DrawCopy("L,P,same");
00292 }
00293 }
00294 c2_1->cd(1);
00295 leg[0]->Draw("same");
00296
00297
00298 printf("::::::::::::: Drift Scan :::::::::::::\n");
00299
00300
00301
00302
00303
00304
00305
00306
00307
00308
00309
00310
00311
00312
00313
00314
00315
00316
00317
00318
00319
00320
00321
00322
00323
00324 TFile inputFile5(ifolder + macroPre[macro] + filename[3],"READ");
00325 for (suid = 0; suid < 4; suid++) {
00326 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00327 newtitle.Form("Efficiency" + macroPre[macro] + filename[3] + "_%i",suid);
00328 if (hEfficiency[suid+4] = (TProfile*)inputFile5.Get(title)){
00329 hEfficiency[suid+4]->SetLineStyle(suid/2+1);
00330 hEfficiency[suid+4]->SetMarkerStyle(25);
00331 hEfficiency[suid+4]->SetMarkerColor(1+suid);
00332 hEfficiency[suid+4]->SetLineColor(1+suid);
00333 dy[suid][1] = hEfficiency[suid+4]->GetBinContent(hEfficiency[suid+4]->FindBin(10));
00334
00335
00336
00337
00338
00339
00340 leg[1]->AddEntry(hEfficiency[suid+4],AV[suid] + DHV[1][suid],"lep");
00341 c2_1->cd(2)->SetLogy(1);
00342 hEfficiency[suid+4]->DrawCopy("L,P,same");
00343 }
00344 }
00345
00346 TFile inputFile6(ifolder + macroPre[macro] + filename[4],"READ");
00347 for (suid = 0; suid < 4; suid++) {
00348 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00349 newtitle.Form("Efficiency" + macroPre[macro] + filename[4] + "_%i",suid);
00350 if (hEfficiency[suid+8] = (TProfile*)inputFile6.Get(title)){
00351 hEfficiency[suid+8]->SetLineStyle(suid/2+1);
00352 hEfficiency[suid+8]->SetMarkerStyle(26);
00353 hEfficiency[suid+8]->SetMarkerColor(1+suid);
00354 hEfficiency[suid+8]->SetLineColor(1+suid);
00355 dy[suid][2] = hEfficiency[suid+8]->GetBinContent(hEfficiency[suid+8]->FindBin(10));
00356
00357
00358
00359
00360
00361
00362 leg[1]->AddEntry(hEfficiency[suid+8],AV[suid] + DHV[2][suid],"lep");
00363 c2_1->cd(2)->SetLogy(1);
00364 hEfficiency[suid+8]->DrawCopy("L,P,same");
00365 }
00366 }
00367
00368 c2_1->cd(2);
00369 leg[1]->Draw("same");
00370 c2_1->SaveAs(ofolder + macroPre[macro] + "MS_HVScan_Comparison_" + mergeName[imerge] + picExt);
00371
00372 c1_1->Clear("D");
00373 c1_2->Clear("D");
00374 c2_1->Clear("D");
00375 c2_2->Clear("D");
00376
00377 for (Int_t i = 0; i < 4; i++)
00378 leg[i]->Clear();
00379 leg[0]->SetNColumns(2);
00380
00381 for (suid = 0; suid < 4; suid++) {
00382 leg[0]->AddEntry(h1dim[0][suid],usedSuName[suid] + radi[suid],"lep");
00383 c2_1->cd(1)->SetLogy(1);
00384 gAHV[suid] = new TGraphErrors(3,ax[suid],ay[suid],eax[suid],eay[suid]);
00385 gAHV[suid]->SetLineStyle(suid/2+1);
00386 gAHV[suid]->SetMarkerStyle(20);
00387 gAHV[suid]->SetMarkerColor(1+suid);
00388 gAHV[suid]->SetLineColor(1+suid);
00389 if (suid == 0){
00390 h1dim[0][suid]->DrawCopy();
00391 gAHV[suid]->Draw("same,P,L");
00392 }
00393
00394 else
00395 gAHV[suid]->Draw("same,P,L");
00396
00397 leg[0]->Draw("same");
00398 c2_1->cd(2)->SetLogy(2);
00399 gDHV[suid] = new TGraphErrors(3,dx[suid],dy[suid],edx[suid],edy[suid]);
00400 gDHV[suid]->SetLineStyle(suid/2+1);
00401 gDHV[suid]->SetMarkerStyle(20);
00402 gDHV[suid]->SetMarkerColor(1+suid);
00403 gDHV[suid]->SetLineColor(1+suid);
00404 if (suid == 0){
00405 h1dim[1][suid]->DrawCopy();
00406 gDHV[suid]->Draw("same,P,L");
00407
00408 }
00409 else
00410 gDHV[suid]->Draw("same,P,L");
00411
00412 leg[0]->Draw("same");
00413 }
00414 c2_1->SaveAs(ofolder + macroPre[macro] + "MS_HVScan_" + mergeName[imerge] + picExt);
00415
00416 c1_1->Clear("D");
00417 c1_2->Clear("D");
00418 c2_1->Clear("D");
00419 c2_2->Clear("D");
00420
00421 inputFile1.Close();
00422 inputFile2.Close();
00423 inputFile3.Close();
00424
00425 inputFile5.Close();
00426 inputFile6.Close();
00427
00428 for (Int_t i = 0; i < 4; i++)
00429 leg[i]->Clear();
00430 leg[0]->SetNColumns(1);
00431 leg[1]->SetNColumns(1);
00432 }
00433
00434
00435
00436 if (dummy) {
00437 printf("::::::::::::: Momentum Scan :::::::::::::\n");
00438
00439 TString momentum[5] = {"2GeV/c","4GeV/c","6GeV/c","8GeV/c","10GeV/c"};
00440 TString radi[4] = {" H:"," H(1/3):"," G(30):"," B:"};
00441 for (Int_t i = 0; i < 4; i++) {
00442 leg[i]->SetHeader(usedSuName[i] + radi[i]);
00443 }
00444 TProfile *h1dim[4];
00445 for (suid = 0; suid < 4; suid++) {
00446 title.Form("Momentum_%i",usedSuId[suid]);
00447 h1dim[suid] = new TProfile(title,usedSuName[suid],5,1.5,10.5);
00448 h1dim[suid]->GetYaxis()->SetRangeUser(1e-3,100);
00449 h1dim[suid]->SetLineStyle(suid/2+1);
00450 h1dim[suid]->SetMarkerStyle(20);
00451 h1dim[suid]->SetMarkerColor(1+suid);
00452 h1dim[suid]->SetLineColor(1+suid);
00453 h1dim[suid]->SetYTitle("pion efficiency [%]");
00454
00455
00456 h1dim[suid]->SetXTitle("momentum [GeV/c]");
00457 }
00458
00459 TFile inputFile1(ifolder + macroPre[macro] + filename[5],"READ");
00460 for (suid = 0; suid < 4; suid++) {
00461 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00462 newtitle.Form("Efficiency" + macroPre[macro] + filename[5] + "_%i",suid);
00463 if (hEfficiency[suid+0] = (TProfile*)inputFile1.Get(title)) {
00464 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
00465 hEfficiency[suid+0]->SetMarkerStyle(20);
00466 hEfficiency[suid+0]->SetMarkerColor(1+suid);
00467 hEfficiency[suid+0]->SetLineColor(1+suid);
00468 h1dim[suid]->Fill(2,hEfficiency[suid+0]->GetBinContent(hEfficiency[suid+0]->FindBin(10)));
00469 leg[suid]->AddEntry(hEfficiency[suid+0],momentum[0],"lep");
00470 c2_2->cd(suid+1)->SetLogy(1);
00471
00472 hEfficiency[suid+0]->DrawCopy("L,P");
00473
00474
00475 }
00476 }
00477 TFile inputFile2(ifolder + macroPre[macro] + filename[6],"READ");
00478 for (suid = 0; suid < 4; suid++) {
00479 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00480 newtitle.Form("Efficiency" + macroPre[macro] + filename[6] + "_%i",suid);
00481 if (hEfficiency[suid+4] = (TProfile*)inputFile2.Get(title)){
00482 hEfficiency[suid+4]->SetLineStyle(suid/2+1);
00483 hEfficiency[suid+4]->SetMarkerStyle(21);
00484 hEfficiency[suid+4]->SetMarkerColor(1+suid);
00485 hEfficiency[suid+4]->SetLineColor(1+suid);
00486 h1dim[suid]->Fill(4,hEfficiency[suid+4]->GetBinContent(hEfficiency[suid+4]->FindBin(10)));
00487 leg[suid]->AddEntry(hEfficiency[suid+4], momentum[1],"lep");
00488 c2_2->cd(suid+1)->SetLogy(1);
00489 hEfficiency[suid+4]->DrawCopy("L,P,same");
00490 }
00491 }
00492 TFile inputFile3(ifolder + macroPre[macro] + filename[7],"READ");
00493 for (suid = 0; suid < 4; suid++) {
00494 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00495 newtitle.Form("Efficiency" + macroPre[macro] + filename[7] + "_%i",suid);
00496 if (hEfficiency[suid+8] = (TProfile*)inputFile3.Get(title)){
00497 hEfficiency[suid+8]->SetLineStyle(suid/2+1);
00498 hEfficiency[suid+8]->SetMarkerStyle(22);
00499 hEfficiency[suid+8]->SetMarkerColor(1+suid);
00500 hEfficiency[suid+8]->SetLineColor(1+suid);
00501 h1dim[suid]->Fill(6,hEfficiency[suid+8]->GetBinContent(hEfficiency[suid+8]->FindBin(10)));
00502 leg[suid]->AddEntry(hEfficiency[suid+8], momentum[2],"lep");
00503 c2_2->cd(suid+1)->SetLogy(1);
00504 hEfficiency[suid+8]->DrawCopy("L,P,same");
00505 }
00506 }
00507 TFile inputFile4(ifolder + macroPre[macro] + filename[8],"READ");
00508 for (suid = 0; suid < 4; suid++) {
00509 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00510 newtitle.Form("Efficiency" + macroPre[macro] + filename[8] + "_%i",suid);
00511 if (hEfficiency[suid+12] = (TProfile*)inputFile4.Get(title)){
00512 hEfficiency[suid+12]->SetLineStyle(suid/2+1);
00513 hEfficiency[suid+12]->SetMarkerStyle(23);
00514 hEfficiency[suid+12]->SetMarkerColor(1+suid);
00515 hEfficiency[suid+12]->SetLineColor(1+suid);
00516 h1dim[suid]->Fill(8,hEfficiency[suid+12]->GetBinContent(hEfficiency[suid+12]->FindBin(10)));
00517 leg[suid]->AddEntry(hEfficiency[suid+12], momentum[3],"lep");
00518 c2_2->cd(suid+1)->SetLogy(1);
00519
00520 hEfficiency[suid+12]->DrawCopy("L,P,same");
00521 }
00522 }
00523 TFile inputFile5(ifolder + macroPre[macro] + filename[9],"READ");
00524 for (suid = 0; suid < 4; suid++) {
00525 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00526 newtitle.Form("Efficiency" + macroPre[macro] + filename[9] + "_%i",suid);
00527 if (hEfficiency[suid+16] = (TProfile*)inputFile5.Get(title)){
00528 hEfficiency[suid+16]->SetLineStyle(suid/2+1);
00529 hEfficiency[suid+16]->SetMarkerStyle(24);
00530 hEfficiency[suid+16]->SetMarkerColor(1+suid);
00531 hEfficiency[suid+16]->SetLineColor(1+suid);
00532 h1dim[suid]->Fill(10,hEfficiency[suid+16]->GetBinContent(hEfficiency[suid+16]->FindBin(10)));
00533 leg[suid]->AddEntry(hEfficiency[suid+16], momentum[4],"lep");
00534 c2_2->cd(suid+1)->SetLogy(1);
00535 hEfficiency[suid+16]->DrawCopy("L,P,same");
00536 }
00537 }
00538 for (Int_t i = 0; i < 4; i++){
00539 c2_2->cd(i+1);
00540 leg[i]->Draw("same");
00541 }
00542 c2_2->SaveAs(ofolder + macroPre[macro] + "MS_MomentumScan_Comparison_" + mergeName[imerge] + picExt);
00543
00544
00545
00546
00547 for (Int_t i = 0; i < 4; i++)
00548 leg[i]->Clear();
00549 leg[0]->SetNColumns(2);
00550 c1_1->cd(1)->SetLogy(1);
00551 for (suid = 0; suid < 4; suid++) {
00552 leg[0]->AddEntry(h1dim[suid],usedSuName[suid] + radi[suid],"lep");
00553 if (suid == 0)
00554 h1dim[suid]->DrawCopy("L,P");
00555 else
00556 h1dim[suid]->DrawCopy("L,P,same");
00557 }
00558 leg[0]->Draw("same");
00559 c1_1->SaveAs(ofolder + macroPre[macro] + "MS_MomentumScan_" + mergeName[imerge] + picExt);
00560
00561 for (Int_t i = 0; i < 4; i++)
00562 leg[i]->Clear();
00563 leg[0]->SetNColumns(1);
00564 c1_1->Clear("D");
00565 c1_2->Clear("D");
00566 c2_1->Clear("D");
00567 c2_2->Clear("D");
00568
00569 inputFile1.Close();
00570 inputFile2.Close();
00571 inputFile3.Close();
00572 inputFile4.Close();
00573 inputFile5.Close();
00574 }
00575
00576
00577 printf("::::::::::::: Radiator Scan :::::::::::::\n");
00578 if (dummy) {
00579
00580 TString rad[15] = {"no radiator","A","B","C","D","E","F","G(10)","G(20)","G(30)","H","H++","I(p)","I(u)","J"};
00581 TProfile *h1dim[4];
00582 for (suid = 0; suid < 4; suid++) {
00583 title.Form("Radiator_%i",usedSuId[suid]);
00584 h1dim[suid] = new TProfile(title,usedSuName[suid],14,-0.5,13.5);
00585 h1dim[suid]->GetYaxis()->SetRangeUser(1e-3,100);
00586 h1dim[suid]->SetLineStyle(suid/2+1);
00587
00588 h1dim[suid]->SetMarkerColor(1+suid);
00589 h1dim[suid]->SetMarkerStyle(24);
00590 h1dim[suid]->SetMarkerSize(1+0.25*suid);
00591 h1dim[suid]->SetLineColor(1+suid);
00592 h1dim[suid]->SetYTitle("pion efficiency [%]");
00593 for (Int_t j = 1; j <= 14; j++)
00594 h1dim[suid]->GetXaxis()->SetBinLabel(j, rad[j-1]);
00595 }
00596 TProfile *temp;
00597 Int_t radiator[8][4] = {{ 7, 7, 7, 7},
00598 { 0, 0, 0, 0},
00599 { 2, 3, 4, 5},
00600 { 3, 4, 5, 6},
00601 { 1, 8,14,12},
00602 {10,12, 1, 8},
00603 {10,13, 9, 2},
00604 {11,13, 9, 2}
00605 };
00606 for(Int_t r = 0; r < 8; r++) {
00607 TFile inputFile(ifolder + macroPre[macro] + filename[11+r],"READ");
00608 for (suid = 0; suid < 4; suid++) {
00609 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00610 if (temp = (TProfile*)inputFile.Get(title)){
00611 h1dim[suid]->Fill(radiator[r][suid],temp->GetBinContent(temp->FindBin(10)));
00612
00613 }
00614 }
00615 inputFile.Close();
00616 }
00617 leg[0]->SetNColumns(2);
00618 leg[0]->SetHeader("Regulare foil and mixed radiators, 3GeV:");
00619 leg[1]->SetNColumns(2);
00620 leg[1]->SetHeader("Irregulare and mixed radiators, 3GeV:");
00621 TFile inputFile0(ifolder + macroPre[macro] + filename[12],"READ");
00622
00623 for (suid = 0; suid < 4; suid++) {
00624 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00625 newtitle.Form("Efficiency" + macroPre[macro] + filename[12] + "_%i",suid);
00626 if (hEfficiency[suid+0] = (TProfile*)inputFile0.Get(title)){
00627 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
00628 hEfficiency[suid+0]->SetMarkerStyle(24);
00629 hEfficiency[suid+0]->SetMarkerColor(1+suid);
00630 hEfficiency[suid+0]->SetLineColor(1+suid);
00631
00632 leg[0]->AddEntry(hEfficiency[suid+0],usedSuName[suid] + rad[0],"lep");
00633 leg[1]->AddEntry(hEfficiency[suid+0],usedSuName[suid] + rad[0],"lep");
00634 c2_1->cd(1)->SetLogy(1);
00635 if (suid == 0)
00636 hEfficiency[suid+0]->DrawCopy("L,P");
00637 else
00638 hEfficiency[suid+0]->DrawCopy("L,P,same");
00639 c2_1->cd(2)->SetLogy(1);
00640 if (suid == 0)
00641 hEfficiency[suid+0]->DrawCopy("L,P");
00642 else
00643 hEfficiency[suid+0]->DrawCopy("L,P,same");
00644 }
00645 }
00646
00647 TFile inputFileG10(ifolder + macroPre[macro] + filename[11],"READ");
00648 for (suid = 0; suid < 4; suid++) {
00649 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00650 newtitle.Form("Efficiency" + macroPre[macro] + filename[11] + "_%i",suid);
00651 if (hEfficiency[suid+4] = (TProfile*)inputFileG10.Get(title)){
00652 hEfficiency[suid+4]->SetLineStyle(suid/2+1);
00653 hEfficiency[suid+4]->SetMarkerStyle(2);
00654 hEfficiency[suid+4]->SetMarkerColor(1+suid);
00655 hEfficiency[suid+4]->SetLineColor(1+suid);
00656
00657 leg[1]->AddEntry(hEfficiency[suid+4],usedSuName[suid] + rad[7],"lep");
00658 c2_1->cd(2)->SetLogy(1);
00659 hEfficiency[suid+4]->DrawCopy("L,P,same");
00660 }
00661 }
00662
00663 TFile inputFile3(ifolder + macroPre[macro] + filename[13],"READ");
00664 TFile inputFile4(ifolder + macroPre[macro] + filename[14],"READ");
00665 TFile inputFile5(ifolder + macroPre[macro] + filename[15],"READ");
00666 TFile inputFile6(ifolder + macroPre[macro] + filename[16],"READ");
00667 TFile inputFile7(ifolder + macroPre[macro] + filename[17],"READ");
00668 TFile inputFile8(ifolder + macroPre[macro] + filename[18],"READ");
00669
00670 suid = 1;
00671 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00672 newtitle.Form("Efficiency" + macroPre[macro] + filename[15] + "_%i",suid);
00673 if (hEfficiency[12] = (TProfile*)inputFile5.Get(title)){
00674 hEfficiency[12]->SetLineStyle(suid/2+1);
00675 hEfficiency[12]->SetMarkerStyle(3);
00676 hEfficiency[12]->SetMarkerColor(1+suid);
00677 hEfficiency[12]->SetLineColor(1+suid);
00678
00679 leg[1]->AddEntry(hEfficiency[12],usedSuName[suid] + rad[8],"lep");
00680 c2_1->cd(2)->SetLogy(1);
00681 hEfficiency[12]->DrawCopy("L,P,same");
00682 }
00683 suid = 3;
00684 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00685 newtitle.Form("Efficiency" + macroPre[macro] + filename[16] + "_%i",suid);
00686 if (hEfficiency[13] = (TProfile*)inputFile6.Get(title)){
00687 hEfficiency[13]->SetLineStyle(suid/2+1);
00688 hEfficiency[13]->SetMarkerStyle(3);
00689 hEfficiency[13]->SetMarkerColor(1+suid);
00690 hEfficiency[13]->SetLineColor(1+suid);
00691
00692 leg[1]->AddEntry(hEfficiency[13],usedSuName[suid] + rad[8],"lep");
00693 c2_1->cd(2)->SetLogy(1);
00694 hEfficiency[13]->DrawCopy("L,P,same");
00695 }
00696
00697 suid = 2;
00698 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00699 newtitle.Form("Efficiency" + macroPre[macro] + filename[17] + "_%i",suid);
00700 if (hEfficiency[20] = (TProfile*)inputFile7.Get(title)){
00701 hEfficiency[20]->SetLineStyle(suid/2+1);
00702 hEfficiency[20]->SetMarkerStyle(5);
00703 hEfficiency[20]->SetMarkerColor(1+suid);
00704 hEfficiency[20]->SetLineColor(1+suid);
00705
00706 leg[1]->AddEntry(hEfficiency[20],usedSuName[suid] + rad[9],"lep");
00707 c2_1->cd(2)->SetLogy(1);
00708 hEfficiency[20]->DrawCopy("L,P,same");
00709 }
00710
00711 suid = 0;
00712 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00713 newtitle.Form("Efficiency" + macroPre[macro] + filename[15] + "_%i",suid);
00714 if (hEfficiency[10] = (TProfile*)inputFile5.Get(title)){
00715 hEfficiency[10]->SetLineStyle(suid/2+1);
00716 hEfficiency[10]->SetMarkerStyle(29);
00717 hEfficiency[10]->SetMarkerColor(1+suid);
00718 hEfficiency[10]->SetLineColor(1+suid);
00719
00720 leg[1]->AddEntry(hEfficiency[10],usedSuName[suid] + rad[1],"lep");
00721 c2_1->cd(2)->SetLogy(1);
00722 hEfficiency[10]->DrawCopy("L,P,same");
00723 }
00724 suid = 2;
00725 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00726 newtitle.Form("Efficiency" + macroPre[macro] + filename[16] + "_%i",suid);
00727 if (hEfficiency[11] = (TProfile*)inputFile6.Get(title)){
00728 hEfficiency[11]->SetLineStyle(suid/2+1);
00729 hEfficiency[11]->SetMarkerStyle(29);
00730 hEfficiency[11]->SetMarkerColor(1+suid);
00731 hEfficiency[11]->SetLineColor(1+suid);
00732
00733 leg[1]->AddEntry(hEfficiency[11],usedSuName[suid] + rad[1],"lep");
00734 c2_1->cd(2)->SetLogy(1);
00735 hEfficiency[11]->DrawCopy("L,P,same");
00736 }
00737
00738 suid = 0;
00739 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00740 newtitle.Form("Efficiency" + macroPre[macro] + filename[13] + "_%i",suid);
00741 if (hEfficiency[8] = (TProfile*)inputFile3.Get(title)){
00742 hEfficiency[8]->SetLineStyle(suid/2+1);
00743 hEfficiency[8]->SetMarkerStyle(25);
00744 hEfficiency[8]->SetMarkerColor(1+suid);
00745 hEfficiency[8]->SetLineColor(1+suid);
00746
00747 leg[0]->AddEntry(hEfficiency[8],usedSuName[suid] + rad[2],"lep");
00748 c2_1->cd(1)->SetLogy(1);
00749 hEfficiency[8]->DrawCopy("L,P,same");
00750 }
00751 suid = 3;
00752 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00753 newtitle.Form("Efficiency" + macroPre[macro] + filename[17] + "_%i",suid);
00754 if (hEfficiency[9] = (TProfile*)inputFile7.Get(title)){
00755 hEfficiency[9]->SetLineStyle(suid/2+1);
00756 hEfficiency[9]->SetMarkerStyle(25);
00757 hEfficiency[9]->SetMarkerColor(1+suid);
00758 hEfficiency[9]->SetLineColor(1+suid);
00759
00760 leg[0]->AddEntry(hEfficiency[9],usedSuName[suid] + rad[2],"lep");
00761 c2_1->cd(1)->SetLogy(1);
00762 hEfficiency[9]->DrawCopy("L,P,same");
00763 }
00764
00765 suid = 1;
00766 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00767 newtitle.Form("Efficiency" + macroPre[macro] + filename[13] + "_%i",suid);
00768 if (hEfficiency[14] = (TProfile*)inputFile3.Get(title)){
00769 hEfficiency[14]->SetLineStyle(suid/2+1);
00770 hEfficiency[14]->SetMarkerStyle(26);
00771 hEfficiency[14]->SetMarkerColor(1+suid);
00772 hEfficiency[14]->SetLineColor(1+suid);
00773 h1dim[suid]->Fill(3,hEfficiency[14]->GetBinContent(hEfficiency[14]->FindBin(10)));
00774 leg[0]->AddEntry(hEfficiency[14],usedSuName[suid] + rad[3],"lep");
00775 c2_1->cd(1)->SetLogy(1);
00776 hEfficiency[14]->DrawCopy("L,P,same");
00777 }
00778 suid = 0;
00779 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00780 newtitle.Form("Efficiency" + macroPre[macro] + filename[14] + "_%i",suid);
00781 if (hEfficiency[15] = (TProfile*)inputFile4.Get(title)){
00782 hEfficiency[15]->SetLineStyle(suid/2+1);
00783 hEfficiency[15]->SetMarkerStyle(26);
00784 hEfficiency[15]->SetMarkerColor(1+suid);
00785 hEfficiency[15]->SetLineColor(1+suid);
00786
00787 leg[0]->AddEntry(hEfficiency[15],usedSuName[suid] + rad[3],"lep");
00788 c2_1->cd(1)->SetLogy(1);
00789 hEfficiency[15]->DrawCopy("L,P,same");
00790 }
00791
00792 suid = 2;
00793 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00794 newtitle.Form("Efficiency" + macroPre[macro] + filename[13] + "_%i",suid);
00795 if (hEfficiency[16] = (TProfile*)inputFile3.Get(title)){
00796 hEfficiency[16]->SetLineStyle(suid/2+1);
00797 hEfficiency[16]->SetMarkerStyle(27);
00798 hEfficiency[16]->SetMarkerColor(1+suid);
00799 hEfficiency[16]->SetLineColor(1+suid);
00800
00801 leg[0]->AddEntry(hEfficiency[16],usedSuName[suid] + rad[4],"lep");
00802 c2_1->cd(1)->SetLogy(1);
00803 hEfficiency[16]->DrawCopy("L,P,same");
00804 }
00805 suid = 1;
00806 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00807 newtitle.Form("Efficiency" + macroPre[macro] + filename[14] + "_%i",suid);
00808 if (hEfficiency[17] = (TProfile*)inputFile4.Get(title)){
00809 hEfficiency[17]->SetLineStyle(suid/2+1);
00810 hEfficiency[17]->SetMarkerStyle(27);
00811 hEfficiency[17]->SetMarkerColor(1+suid);
00812 hEfficiency[17]->SetLineColor(1+suid);
00813
00814 leg[0]->AddEntry(hEfficiency[17],usedSuName[suid] + rad[4],"lep");
00815 c2_1->cd(1)->SetLogy(1);
00816 hEfficiency[17]->DrawCopy("L,P,same");
00817 }
00818
00819 suid = 3;
00820 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00821 newtitle.Form("Efficiency" + macroPre[macro] + filename[13] + "_%i",suid);
00822 if (hEfficiency[18] = (TProfile*)inputFile3.Get(title)){
00823 hEfficiency[18]->SetLineStyle(suid/2+1);
00824 hEfficiency[18]->SetMarkerStyle(28);
00825 hEfficiency[18]->SetMarkerColor(1+suid);
00826 hEfficiency[18]->SetLineColor(1+suid);
00827
00828 leg[0]->AddEntry(hEfficiency[18],usedSuName[suid] + rad[5],"lep");
00829 c2_1->cd(1)->SetLogy(1);
00830 hEfficiency[18]->DrawCopy("L,P,same");
00831 }
00832 suid = 2;
00833 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00834 newtitle.Form("Efficiency" + macroPre[macro] + filename[14] + "_%i",suid);
00835 if (hEfficiency[19] = (TProfile*)inputFile4.Get(title)){
00836 hEfficiency[19]->SetLineStyle(suid/2+1);
00837 hEfficiency[19]->SetMarkerStyle(28);
00838 hEfficiency[19]->SetMarkerColor(1+suid);
00839 hEfficiency[19]->SetLineColor(1+suid);
00840
00841 leg[0]->AddEntry(hEfficiency[19],usedSuName[suid] + rad[5],"lep");
00842 c2_1->cd(1)->SetLogy(1);
00843 hEfficiency[19]->DrawCopy("L,P,same");
00844 }
00845
00846 suid = 3;
00847 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00848 newtitle.Form("Efficiency" + macroPre[macro] + filename[14] + "_%i",suid);
00849 if (hEfficiency[26] = (TProfile*)inputFile4.Get(title)){
00850 hEfficiency[26]->SetLineStyle(suid/2+1);
00851 hEfficiency[26]->SetMarkerStyle(20);
00852 hEfficiency[26]->SetMarkerColor(1+suid);
00853 hEfficiency[26]->SetLineColor(1+suid);
00854
00855 leg[0]->AddEntry(hEfficiency[26],usedSuName[suid] + rad[6],"lep");
00856 c2_1->cd(1)->SetLogy(1);
00857 hEfficiency[26]->DrawCopy("L,P,same");
00858 }
00859
00860 suid = 3;
00861 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00862 newtitle.Form("Efficiency" + macroPre[macro] + filename[15] + "_%i",suid);
00863 if (hEfficiency[23] = (TProfile*)inputFile5.Get(title)){
00864 hEfficiency[23]->SetLineStyle(suid/2+1);
00865 hEfficiency[23]->SetMarkerStyle(30);
00866 hEfficiency[23]->SetMarkerColor(1+suid);
00867 hEfficiency[23]->SetLineColor(1+suid);
00868
00869 leg[0]->AddEntry(hEfficiency[23],usedSuName[suid] + rad[12],"lep");
00870 c2_1->cd(1)->SetLogy(1);
00871 hEfficiency[23]->DrawCopy("L,P,same");
00872 leg[1]->AddEntry(hEfficiency[23],usedSuName[suid] + rad[12],"lep");
00873 c2_1->cd(2)->SetLogy(1);
00874 hEfficiency[23]->DrawCopy("L,P,same");
00875 }
00876 suid = 1;
00877 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00878 newtitle.Form("Efficiency" + macroPre[macro] + filename[16] + "_%i",suid);
00879 if (hEfficiency[24] = (TProfile*)inputFile6.Get(title)){
00880 hEfficiency[24]->SetLineStyle(suid/2+1);
00881 hEfficiency[24]->SetMarkerStyle(30);
00882 hEfficiency[24]->SetMarkerColor(1+suid);
00883 hEfficiency[24]->SetLineColor(1+suid);
00884
00885 leg[0]->AddEntry(hEfficiency[24],usedSuName[suid] + rad[12],"lep");
00886 c2_1->cd(1)->SetLogy(1);
00887 hEfficiency[24]->DrawCopy("L,P,same");
00888 leg[1]->AddEntry(hEfficiency[24],usedSuName[suid] + rad[12],"lep");
00889 c2_1->cd(2)->SetLogy(1);
00890 hEfficiency[24]->DrawCopy("L,P,same");
00891 }
00892
00893 suid = 1;
00894 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00895 newtitle.Form("Efficiency" + macroPre[macro] + filename[17] + "_%i",suid);
00896 if (hEfficiency[25] = (TProfile*)inputFile7.Get(title)){
00897 hEfficiency[25]->SetLineStyle(suid/2+1);
00898 hEfficiency[25]->SetMarkerStyle(29);
00899 hEfficiency[25]->SetMarkerColor(1+suid);
00900 hEfficiency[25]->SetLineColor(1+suid);
00901
00902 leg[0]->AddEntry(hEfficiency[25],usedSuName[suid] + rad[13],"lep");
00903 c2_1->cd(1)->SetLogy(1);
00904 hEfficiency[25]->DrawCopy("L,P,same");
00905 leg[1]->AddEntry(hEfficiency[25],usedSuName[suid] + rad[13],"lep");
00906 c2_1->cd(2)->SetLogy(1);
00907 hEfficiency[25]->DrawCopy("L,P,same");
00908 }
00909
00910 suid = 0;
00911 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00912 newtitle.Form("Efficiency" + macroPre[macro] + filename[16] + "_%i",suid);
00913 if (hEfficiency[21] = (TProfile*)inputFile6.Get(title)){
00914 hEfficiency[21]->SetLineStyle(suid/2+1);
00915 hEfficiency[21]->SetMarkerStyle(21);
00916 hEfficiency[21]->SetMarkerColor(1+suid);
00917 hEfficiency[21]->SetLineColor(1+suid);
00918
00919 leg[1]->AddEntry(hEfficiency[21],usedSuName[suid] + rad[10],"lep");
00920 c2_1->cd(2)->SetLogy(1);
00921 hEfficiency[21]->DrawCopy("L,P,same");
00922 }
00923
00924 suid = 0;
00925 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00926 newtitle.Form("Efficiency" + macroPre[macro] + filename[18] + "_%i",suid);
00927 if (hEfficiency[22] = (TProfile*)inputFile8.Get(title)){
00928 hEfficiency[22]->SetLineStyle(suid/2+1);
00929 hEfficiency[22]->SetMarkerStyle(22);
00930 hEfficiency[22]->SetMarkerColor(1+suid);
00931 hEfficiency[22]->SetLineColor(1+suid);
00932
00933 leg[1]->AddEntry(hEfficiency[22],usedSuName[suid] + rad[11],"lep");
00934 c2_1->cd(2)->SetLogy(1);
00935 hEfficiency[22]->DrawCopy("L,P,same");
00936 }
00937 c2_1->cd(1);
00938 leg[0]->Draw("same");
00939 c2_1->cd(2);
00940 leg[1]->Draw("same");
00941 c2_1->SaveAs(ofolder + macroPre[macro] + "MS_Radiator_ShortComparison_" + mergeName[imerge] + picExt);
00942
00943 for (Int_t i = 0; i < 4; i++)
00944 leg[i]->Clear();
00945 leg[0]->SetNColumns(2);
00946 c1_1->cd(1)->SetLogy(1);
00947 TLine designGoal(-0.5,1,13.5,1);
00948 designGoal.SetLineColor(1);
00949 designGoal.SetLineStyle(2);
00950 for (suid = 0; suid < 4; suid++) {
00951 leg[0]->AddEntry(h1dim[suid],usedSuName[suid],"lep");
00952 if (suid == 0){
00953 h1dim[suid]->DrawCopy("P");
00954 designGoal.Draw("same");
00955 }
00956 else
00957 h1dim[suid]->DrawCopy("P,same");
00958 }
00959 leg[0]->Draw("same");
00960 c1_1->SaveAs(ofolder + macroPre[macro] + "MS_Radiator_" + mergeName[imerge] + picExt);
00961
00962 for (Int_t i = 0; i < 4; i++)
00963 leg[i]->Clear();
00964 leg[0]->SetNColumns(1);
00965 leg[1]->SetNColumns(1);
00966 c1_1->Clear("D");
00967 c1_2->Clear("D");
00968 c2_1->Clear("D");
00969 c2_2->Clear("D");
00970
00971 inputFile0.Close();
00972 inputFileG10.Close();
00973 inputFile3.Close();
00974 inputFile4.Close();
00975 inputFile5.Close();
00976 inputFile6.Close();
00977 inputFile7.Close();
00978 inputFile8.Close();
00979 }
00980
00981 if (dummy) {
00982 TFile inputFile1(ifolder + macroPre[macro] + filename[11],"READ");
00983 TString rad[4] = {"G(10)","10 fiber mats","10 fiber mats","10 fiber mats"};
00984 leg[0]->SetHeader(rad[0] + ", 3GeV:");
00985 for (suid = 0; suid < 4; suid++) {
00986 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
00987 newtitle.Form("Efficiency" + macroPre[macro] + filename[11] + "_%i",suid);
00988 if (hEfficiency[suid+0] = (TProfile*)inputFile1.Get(title)){
00989 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
00990 hEfficiency[suid+0]->SetMarkerStyle(20);
00991 hEfficiency[suid+0]->SetMarkerColor(1+suid);
00992 hEfficiency[suid+0]->SetLineColor(1+suid);
00993 leg[0]->AddEntry(hEfficiency[suid+0],usedSuName[suid],"lep");
00994 c1_1->cd(1)->SetLogy(1);
00995 if (suid == 0)
00996 hEfficiency[suid+0]->DrawCopy("L,P");
00997 else
00998 hEfficiency[suid+0]->DrawCopy("L,P,same");
00999 }
01000 }
01001 c1_1->cd(1);
01002 leg[0]->Draw("same");
01003 c1_1->SaveAs(ofolder + macroPre[macro] + "MS_sameRadiator_Comparison_" + mergeName[imerge] + picExt);
01004
01005 for (Int_t i = 0; i < 4; i++)
01006 leg[i]->Clear();
01007
01008 c1_1->Clear("D");
01009 c1_2->Clear("D");
01010 c2_1->Clear("D");
01011 c2_2->Clear("D");
01012
01013 inputFile1.Close();
01014 }
01015
01016 if (dummy) {
01017 TFile inputFile1(ifolder + macroPre[macro] + filename[12],"READ");
01018 TString rad[4] = {"no radiator","no radiator","no radiator","no radiator"};
01019 leg[0]->SetHeader(rad[0] + ", 3GeV:");
01020 for (suid = 0; suid < 4; suid++) {
01021 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01022 newtitle.Form("Efficiency" + macroPre[macro] + filename[12] + "_%i",suid);
01023 if (hEfficiency[suid+0] = (TProfile*)inputFile1.Get(title)){
01024 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
01025 hEfficiency[suid+0]->SetMarkerStyle(24);
01026 hEfficiency[suid+0]->SetMarkerColor(1+suid);
01027 hEfficiency[suid+0]->SetLineColor(1+suid);
01028 leg[0]->AddEntry(hEfficiency[suid+0],usedSuName[suid],"lep");
01029 c1_1->cd(1)->SetLogy(1);
01030 if (suid == 0)
01031 hEfficiency[suid+0]->DrawCopy("L,P");
01032 else
01033 hEfficiency[suid+0]->DrawCopy("L,P,same");
01034 }
01035 }
01036 c1_1->cd(1);
01037 leg[0]->Draw("same");
01038 c1_1->SaveAs(ofolder + macroPre[macro] + "MS_noRadiator_Comparison_" + mergeName[imerge] + picExt);
01039
01040 for (Int_t i = 0; i < 4; i++)
01041 leg[i]->Clear();
01042
01043 c1_1->Clear("D");
01044 c1_2->Clear("D");
01045 c2_1->Clear("D");
01046 c2_2->Clear("D");
01047
01048 inputFile1.Close();
01049 }
01050
01051 if (dummy) {
01052 TFile inputFile1(ifolder + macroPre[macro] + filename[13],"READ");
01053 TString det[4][6] = {{"B reg. Pokalon (250/700/24)","C reg. PP (200/700/15)","A ALICE","H Foam","H Foam","H++ Foam"},
01054 {"C reg. PP (200/700/15)","D reg. PP (100/500/15)","G(20) Fiber","I(p) Mixed PP-foam Pokalon","I(u) Mixed PP-foam Pokalon","I(u) Mixed PP-foam Pokalon"},
01055 {"D reg. PP (100/500/15)","E reg. PP (120/500/20)","J(10) Fiber","A ALICE","G(30)Fiber","G(30) Fiber"},
01056 {"E reg. PP (120/500/20)","F reg. PP (220/250/20)","I(p) Mixed PP-foam Pokalon","G(20) Fiber","B reg. Pokalon (250/700/24)","B reg. Pokalon (250/700/24)"}
01057 };
01058 for (Int_t i = 0; i < 4; i++) {
01059 leg[i]->SetHeader(usedSuName[i] + "3GeV:");
01060 }
01061 for (suid = 0; suid < 4; suid++) {
01062 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01063 newtitle.Form("Efficiency" + macroPre[macro] + filename[13] + "_%i",suid);
01064 if (hEfficiency[suid+0] = (TProfile*)inputFile1.Get(title)){
01065 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
01066 hEfficiency[suid+0]->SetMarkerStyle(20);
01067 hEfficiency[suid+0]->SetMarkerColor(1+suid);
01068 hEfficiency[suid+0]->SetLineColor(1+suid);
01069 leg[suid]->AddEntry(hEfficiency[suid+0], det[suid][0],"lep");
01070 c2_2->cd(suid+1)->SetLogy(1);
01071 hEfficiency[suid+0]->DrawCopy("L,P");
01072 }
01073 }
01074 TFile inputFile2(ifolder + macroPre[macro] + filename[14],"READ");
01075 for (suid = 0; suid < 4; suid++) {
01076 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01077 newtitle.Form("Efficiency" + macroPre[macro] + filename[14] + "_%i",suid);
01078 if (hEfficiency[suid+4] = (TProfile*)inputFile2.Get(title)){
01079 hEfficiency[suid+4]->SetLineStyle(suid/2+1);
01080 hEfficiency[suid+4]->SetMarkerStyle(21);
01081 hEfficiency[suid+4]->SetMarkerColor(1+suid);
01082 hEfficiency[suid+4]->SetLineColor(1+suid);
01083 leg[suid]->AddEntry(hEfficiency[suid+4], det[suid][1],"lep");
01084 c2_2->cd(suid+1)->SetLogy(1);
01085 hEfficiency[suid+4]->DrawCopy("L,P,same");
01086 }
01087 }
01088 TFile inputFile3(ifolder + macroPre[macro] + filename[15],"READ");
01089 for (suid = 0; suid < 4; suid++) {
01090 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01091 newtitle.Form("Efficiency" + macroPre[macro] + filename[15] + "_%i",suid);
01092 if (hEfficiency[suid+8] = (TProfile*)inputFile3.Get(title)){
01093 hEfficiency[suid+8]->SetLineStyle(suid/2+1);
01094 hEfficiency[suid+8]->SetMarkerStyle(22);
01095 hEfficiency[suid+8]->SetMarkerColor(1+suid);
01096 hEfficiency[suid+8]->SetLineColor(1+suid);
01097 leg[suid]->AddEntry(hEfficiency[suid+8],det[suid][2],"lep");
01098 c2_2->cd(suid+1)->SetLogy(1);
01099 hEfficiency[suid+8]->DrawCopy("L,P,same");
01100 }
01101 }
01102 TFile inputFile4(ifolder + macroPre[macro] + filename[16],"READ");
01103 for (suid = 0; suid < 4; suid++) {
01104 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01105 newtitle.Form("Efficiency" + macroPre[macro] + filename[16] + "_%i",suid);
01106 if (hEfficiency[suid+12] = (TProfile*)inputFile4.Get(title)){
01107 hEfficiency[suid+12]->SetLineStyle(suid/2+1);
01108 hEfficiency[suid+12]->SetMarkerStyle(23);
01109 hEfficiency[suid+12]->SetMarkerColor(1+suid);
01110 hEfficiency[suid+12]->SetLineColor(1+suid);
01111 leg[suid]->AddEntry(hEfficiency[suid+12], det[suid][3],"lep");
01112 c2_2->cd(suid+1)->SetLogy(1);
01113 hEfficiency[suid+12]->DrawCopy("L,P,same");
01114 }
01115 }
01116 TFile inputFile5(ifolder + macroPre[macro] + filename[17],"READ");
01117 for (suid = 0; suid < 4; suid++) {
01118 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01119 newtitle.Form("Efficiency" + macroPre[macro] + filename[17] + "_%i",suid);
01120 if (hEfficiency[suid+16] = (TProfile*)inputFile5.Get(title)){
01121 hEfficiency[suid+16]->SetLineStyle(suid/2+1);
01122 hEfficiency[suid+16]->SetMarkerStyle(24);
01123 hEfficiency[suid+16]->SetMarkerColor(1+suid);
01124 hEfficiency[suid+16]->SetLineColor(1+suid);
01125 leg[suid]->AddEntry(hEfficiency[suid+16], det[suid][4],"lep");
01126 c2_2->cd(suid+1)->SetLogy(1);
01127 hEfficiency[suid+16]->DrawCopy("L,P,same");
01128 }
01129 }
01130 TFile inputFile6(ifolder + macroPre[macro] + filename[18],"READ");
01131 for (suid = 0; suid < 4; suid++) {
01132 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01133 newtitle.Form("Efficiency" + macroPre[macro] + filename[18] + "_%i",suid);
01134 if (hEfficiency[suid+20] = (TProfile*)inputFile6.Get(title)){
01135 hEfficiency[suid+20]->SetLineStyle(suid/2+1);
01136 hEfficiency[suid+20]->SetMarkerStyle(25);
01137 hEfficiency[suid+20]->SetMarkerColor(1+suid);
01138 hEfficiency[suid+20]->SetLineColor(1+suid);
01139 leg[suid]->AddEntry(hEfficiency[suid+20], det[suid][5],"lep");
01140 c2_2->cd(suid+1)->SetLogy(1);
01141 hEfficiency[suid+20]->DrawCopy("L,P,same");
01142 }
01143 }
01144 for (Int_t i = 0; i < 4; i++){
01145 c2_2->cd(i+1);
01146 leg[i]->Draw("same");
01147 }
01148 c2_2->SaveAs(ofolder + macroPre[macro] + "MS_RadiatorScan_Comparison_" + mergeName[imerge] + picExt);
01149
01150 for (Int_t i = 0; i < 4; i++)
01151 leg[i]->Clear();
01152 leg[0]->SetNColumns(4);
01153
01154 c1_1->Clear("D");
01155 c1_2->Clear("D");
01156 c2_1->Clear("D");
01157 c2_2->Clear("D");
01158
01159 inputFile1.Close();
01160 inputFile2.Close();
01161 inputFile3.Close();
01162 inputFile4.Close();
01163 inputFile5.Close();
01164 inputFile6.Close();
01165 }
01166
01167 printf("::::::::::::: Foam comparison ::::::::::\n");
01168 if (dummy) {
01169 TString det[4][4] = {{"H Foam 30cm","H Foam 6cm","G(30)Fiber","B reg. Pokalon (250/700/24)"},
01170 {"H Foam 30cm","H Foam 6cm","G(30)Fiber","B reg. Pokalon (250/700/24)"},
01171 {"Black Foam 26.5cm","H Foam 6cm","G(30)Fiber","B reg. Pokalon (250/700/24)"},
01172 {"White Foam 27cm","H Foam 6cm","G(30)Fiber","B reg. Pokalon (250/700/24)"}
01173 };
01174 leg[0]->SetNColumns(2);
01175 leg[0]->SetHeader("2GeV:");
01176 leg[1]->SetNColumns(2);
01177 leg[1]->SetHeader("4GeV:");
01178 for (suid = 0; suid < 2; suid++) {
01179 leg[0]->AddEntry((TObject*)0,usedSuName[suid],"");
01180 leg[1]->AddEntry((TObject*)0,usedSuName[suid],"");
01181 }
01182
01183 TFile inputFile2GeVRef(ifolder + macroPre[macro] + filename[5],"READ");
01184 if (!inputFile2GeVRef.IsOpen()) cout << "inputFile2GeVRef" << endl;
01185 for (suid = 0; suid < 2; suid++) {
01186 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01187 newtitle.Form("Efficiency" + macroPre[macro] + filename[5] + "_%i",suid);
01188 if (hEfficiency[suid+0] = (TProfile*)inputFile2GeVRef.Get(title)){
01189 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
01190 hEfficiency[suid+0]->SetMarkerStyle(21);
01191 hEfficiency[suid+0]->SetMarkerColor(1+suid);
01192 hEfficiency[suid+0]->SetLineColor(1+suid);
01193 leg[0]->AddEntry(hEfficiency[suid+0], det[0][suid],"lep");
01194 c2_1->cd(1)->SetLogy(1);
01195 if (suid == 0)
01196 hEfficiency[suid+0]->DrawCopy("L,P");
01197 else
01198 hEfficiency[suid+0]->DrawCopy("L,P,same");
01199 }
01200 }
01201 TFile inputFile4GeVRef(ifolder + macroPre[macro] + filename[6],"READ");
01202 if (!inputFile4GeVRef.IsOpen()) cout << "inputFile4GeVRef" << endl;
01203 for (suid = 0; suid < 2; suid++) {
01204 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01205 newtitle.Form("Efficiency" + macroPre[macro] + filename[6] + "_%i",suid);
01206 if (hEfficiency[suid+4] = (TProfile*)inputFile4GeVRef.Get(title)){
01207 hEfficiency[suid+4]->SetLineStyle(suid/2+1);
01208 hEfficiency[suid+4]->SetMarkerStyle(21);
01209 hEfficiency[suid+4]->SetMarkerColor(1+suid);
01210 hEfficiency[suid+4]->SetLineColor(1+suid);
01211 leg[1]->AddEntry(hEfficiency[suid+4], det[1][suid],"lep");
01212 c2_1->cd(2)->SetLogy(1);
01213 if (suid == 0)
01214 hEfficiency[suid+4]->DrawCopy("L,P");
01215 else
01216 hEfficiency[suid+4]->DrawCopy("L,P,same");
01217 }
01218 }
01219 TFile inputFileWhiteFoam2GeV(ifolder + macroPre[macro] + filename[23],"READ");
01220 if (!inputFileWhiteFoam2GeV.IsOpen()) cout << TString(ifolder + macroPre[macro] + filename[23]).Data() << endl;
01221 for (suid = 0; suid < 2; suid++) {
01222 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01223 newtitle.Form("Efficiency" + macroPre[macro] + filename[23] + "_%i",suid);
01224 if (hEfficiency[suid+8] = (TProfile*)inputFileWhiteFoam2GeV.Get(title)){
01225 hEfficiency[suid+8]->SetLineStyle(suid/2+1);
01226 hEfficiency[suid+8]->SetMarkerStyle(24);
01227 hEfficiency[suid+8]->SetMarkerColor(1+suid);
01228 hEfficiency[suid+8]->SetLineColor(1+suid);
01229 leg[0]->AddEntry(hEfficiency[suid+8], det[3][suid],"lep");
01230 c2_1->cd(1)->SetLogy(1);
01231 hEfficiency[suid+8]->DrawCopy("L,P,same");
01232 }
01233 }
01234 TFile inputFileWhiteFoam4GeV(ifolder + macroPre[macro] + filename[22],"READ");
01235 if (!inputFileWhiteFoam4GeV.IsOpen()) cout << TString(ifolder + macroPre[macro] + filename[22]).Data() << endl;
01236 for (suid = 0; suid < 2; suid++) {
01237 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01238 newtitle.Form("Efficiency" + macroPre[macro] + filename[22] + "_%i",suid);
01239 if (hEfficiency[suid+12] = (TProfile*)inputFileWhiteFoam4GeV.Get(title)){
01240 hEfficiency[suid+12]->SetLineStyle(suid/2+1);
01241 hEfficiency[suid+12]->SetMarkerStyle(24);
01242 hEfficiency[suid+12]->SetMarkerColor(1+suid);
01243 hEfficiency[suid+12]->SetLineColor(1+suid);
01244 leg[1]->AddEntry(hEfficiency[suid+12], det[3][suid],"lep");
01245 c2_1->cd(2)->SetLogy(1);
01246 hEfficiency[suid+12]->DrawCopy("L,P,same");
01247 }
01248 }
01249 TFile inputFileBlackFoam4GeV(ifolder + macroPre[macro] + filename[21],"READ");
01250 if (!inputFileBlackFoam4GeV.IsOpen()) cout << "inputFileBlackFoam4GeV" << endl;
01251 for (suid = 0; suid < 2; suid++) {
01252 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01253 newtitle.Form("Efficiency" + macroPre[macro] + filename[21] + "_%i",suid);
01254 if (hEfficiency[suid+16] = (TProfile*)inputFileBlackFoam4GeV.Get(title)){
01255 hEfficiency[suid+16]->SetLineStyle(suid/2+1);
01256 hEfficiency[suid+16]->SetMarkerStyle(20);
01257 hEfficiency[suid+16]->SetMarkerColor(1+suid);
01258 hEfficiency[suid+16]->SetLineColor(1+suid);
01259 leg[1]->AddEntry(hEfficiency[suid+16], det[2][suid],"lep");
01260 c2_1->cd(2)->SetLogy(1);
01261 hEfficiency[suid+16]->DrawCopy("L,P,same");
01262 }
01263 }
01264
01265 c2_1->cd(1);
01266 leg[0]->Draw("same");
01267 c2_1->cd(2);
01268 leg[1]->Draw("same");
01269 c2_1->SaveAs(ofolder + macroPre[macro] + "MS_Foam_Comparison_" + mergeName[imerge] + picExt);
01270 for (Int_t i = 0; i < 4; i++)
01271 leg[i]->Clear();
01272 leg[0]->SetNColumns(4);
01273
01274 c1_1->Clear("D");
01275 c1_2->Clear("D");
01276 c2_1->Clear("D");
01277 c2_2->Clear("D");
01278
01279 inputFile2GeVRef.Close();
01280 inputFile4GeVRef.Close();
01281 inputFileBlackFoam4GeV.Close();
01282 inputFileWhiteFoam4GeV.Close();
01283 inputFileWhiteFoam2GeV.Close();
01284 }
01285
01286 printf("::::::::::::: Modified Setup :::::::::::\n");
01287 if (dummy) {
01288 TString det[2][4] = {{"H Foam 30cm","H Foam 6cm","G(30)Fiber","B reg. Pokalon (250/700/24)"},
01289 {"H Foam 30cm","H Foam 6cm","A ALICE","B reg. Pokalon (250/700/24)"}
01290 };
01291 leg[0]->SetNColumns(4);
01292 leg[0]->SetHeader("Rotated, 4GeV: ");
01293 leg[1]->SetNColumns(4);
01294 leg[1]->SetHeader("Rotated & Tilted, 2GeV: ");
01295 for (suid = 0; suid < 4; suid++) {
01296 leg[0]->AddEntry((TObject*)0,usedSuName[suid],"");
01297 leg[1]->AddEntry((TObject*)0,usedSuName[suid],"");
01298 }
01299 TFile inputFileRot(ifolder + macroPre[macro] + filename[19],"READ");
01300 for (suid = 0; suid < 4; suid++) {
01301 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01302 newtitle.Form("Efficiency" + macroPre[macro] + filename[19] + "_%i",suid);
01303 if (hEfficiency[suid+0] = (TProfile*)inputFileRot.Get(title)){
01304 hEfficiency[suid+0]->SetLineStyle(suid/2+1);
01305 hEfficiency[suid+0]->SetMarkerStyle(20);
01306 hEfficiency[suid+0]->SetMarkerColor(1+suid);
01307 hEfficiency[suid+0]->SetLineColor(1+suid);
01308 leg[0]->AddEntry(hEfficiency[suid+0], det[0][suid],"lep");
01309 c2_1->cd(1)->SetLogy(1);
01310 if (suid == 0)
01311 hEfficiency[suid+0]->DrawCopy("L,P");
01312 else
01313 hEfficiency[suid+0]->DrawCopy("L,P,same");
01314 }
01315 }
01316 TFile inputFile4GeV(ifolder + macroPre[macro] + filename[14],"READ");
01317 for (suid = 0; suid < 4; suid++) {
01318 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01319 newtitle.Form("Efficiency" + macroPre[macro] + filename[14] + "_%i",suid);
01320 if (hEfficiency[suid+8] = (TProfile*)inputFile4GeV.Get(title)){
01321 hEfficiency[suid+8]->SetLineStyle(suid/2+1);
01322 hEfficiency[suid+8]->SetMarkerStyle(21);
01323 hEfficiency[suid+8]->SetMarkerColor(1+suid);
01324 hEfficiency[suid+8]->SetLineColor(1+suid);
01325 leg[0]->AddEntry(hEfficiency[suid+8], det[0][suid],"lep");
01326 c2_1->cd(1)->SetLogy(1);
01327 hEfficiency[suid+8]->DrawCopy("L,P,same");
01328 }
01329 }
01330 TFile inputFileRotTil(ifolder + macroPre[macro] + filename[20],"READ");
01331 for (suid = 0; suid < 4; suid++) {
01332 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01333 newtitle.Form("Efficiency" + macroPre[macro] + filename[20] + "_%i",suid);
01334 if (hEfficiency[suid+4] = (TProfile*)inputFileRotTil.Get(title)){
01335 hEfficiency[suid+4]->SetLineStyle(suid/2+1);
01336 hEfficiency[suid+4]->SetMarkerStyle(24);
01337 hEfficiency[suid+4]->SetMarkerColor(1+suid);
01338 hEfficiency[suid+4]->SetLineColor(1+suid);
01339 leg[1]->AddEntry(hEfficiency[suid+4], det[1][suid],"lep");
01340 c2_1->cd(2)->SetLogy(1);
01341 if (suid == 0)
01342 hEfficiency[suid+4]->DrawCopy("L,P");
01343 else
01344 hEfficiency[suid+4]->DrawCopy("L,P,same");
01345 }
01346 }
01347 TFile inputFile2GeV(ifolder + macroPre[macro] + filename[12],"READ");
01348 for (suid = 0; suid < 4; suid++) {
01349 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01350 newtitle.Form("Efficiency" + macroPre[macro] + filename[12] + "_%i",suid);
01351 if (hEfficiency[suid+12] = (TProfile*)inputFile2GeV.Get(title)){
01352 hEfficiency[suid+12]->SetLineStyle(suid/2+1);
01353 hEfficiency[suid+12]->SetMarkerStyle(25);
01354 hEfficiency[suid+12]->SetMarkerColor(1+suid);
01355 hEfficiency[suid+12]->SetLineColor(1+suid);
01356 leg[1]->AddEntry(hEfficiency[suid+12], det[0][suid],"lep");
01357 c2_1->cd(2)->SetLogy(1);
01358 hEfficiency[suid+12]->DrawCopy("L,P,same");
01359 }
01360 }
01361 c2_1->cd(1);
01362 leg[0]->Draw("same");
01363 c2_1->cd(2);
01364 leg[1]->Draw("same");
01365 c2_1->SaveAs(ofolder + macroPre[macro] + "MS_Modified_Setup_" + mergeName[imerge] + picExt);
01366 for (Int_t i = 0; i < 4; i++)
01367 leg[i]->Clear();
01368 leg[0]->SetNColumns(4);
01369
01370 c1_1->Clear("D");
01371 c1_2->Clear("D");
01372 c2_1->Clear("D");
01373 c2_2->Clear("D");
01374 }
01375
01376 printf("::::::::::::: Signal Scan :::::::::::::\n");
01377 if (dummy) {
01378 Int_t file = 11;
01379 for (Int_t i = 0; i < 4; i++) {
01380 leg[i]->SetNColumns(2);
01381 leg[i]->SetHeader(usedSuName[i] + ", G(10), 3GeV:");
01382 }
01383 TFile inputFile1(ifolder + macroPre[macro] + filename[file],"READ");
01384 TString rad[4] = {"10 fiber mats","10 fiber mats","10 fiber mats","10 fiber mats"};
01385 Int_t rese[6] = {0,1,3,4,6,7};
01386 TString ReSe[6] = {"full cluster full integral","3 pads full integral","full cluster time window integral","3 pads time window integral","full cluster amplitude"," 3 pads amplitude"};
01387 Int_t markerS[6] = {20,24,21,25,22,26};
01388 for (Int_t j = 0; j < 6; j++) {
01389 for (suid = 0; suid < 4; suid++) {
01390 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rese[j]].Data(),mergeName[imerge].Data(),filterName[ifilter].Data(),irebin);
01391 newtitle.Form("Efficiency" + macroPre[macro] + filename[file] + "_%i",suid);
01392 if (hEfficiency[j*4+suid] = (TProfile*)inputFile1.Get(title)){
01393 hEfficiency[j*4+suid]->SetLineStyle(suid/2+1);
01394 hEfficiency[j*4+suid]->SetMarkerStyle(markerS[j]);
01395 hEfficiency[j*4+suid]->SetMarkerColor(1+suid);
01396 hEfficiency[j*4+suid]->SetLineColor(1+suid);
01397 leg[suid]->AddEntry(hEfficiency[j*4+suid],ReSe[j],"lep");
01398 c2_2->cd(suid+1)->SetLogy(1);
01399 if (j == 0)
01400 hEfficiency[j*4+suid]->DrawCopy("L,P");
01401 else
01402 hEfficiency[j*4+suid]->DrawCopy("L,P,same");
01403 }
01404 }
01405 }
01406 for (suid = 0; suid < 4; suid++) {
01407 c2_2->cd(suid+1);
01408 leg[suid]->Draw("same");
01409 }
01410 c2_2->SaveAs(ofolder + macroPre[macro] + "MS_Signal_Comparison_" + mergeName[imerge] + picExt);
01411
01412 for (Int_t i = 0; i < 4; i++) {
01413 leg[i]->Clear();
01414 leg[i]->SetNColumns(1);
01415 }
01416 c1_1->Clear("D");
01417 c1_2->Clear("D");
01418 c2_1->Clear("D");
01419 c2_2->Clear("D");
01420
01421 inputFile1.Close();
01422 }
01423
01424 printf("::::::::::::: With or Without Scan :::::::::::::\n");
01425 if (dummy) {
01426 Int_t file = 11;
01427 suid = 1;
01428 TFile inputFile1(ifolder + macroPre[macro] + filename[file],"READ");
01429 TString FiSe[4] = {": with overflow; with underflow",": without overflow; with underflow",": without overflow; without underflow",": with overflow; without underflow"
01430
01431 };
01432 Int_t rese[6] = {0,1,3,4,6,7};
01433 TString ReSe[6] = {"full cluster full integral","3 pads full integral","full cluster time window integral","3 pads time window integral","full cluster amplitude"," 3 pads amplitude"};
01434 Int_t markerS[6] = {20,24,21,25,22,26};
01435 for (Int_t i = 0; i < 4; i++){
01436 leg[i]->SetHeader(usedSuName[suid] + ", G(10), 3GeV, " + FiSe[i]);
01437 leg[i]->SetNColumns(2);
01438 }
01439
01440 for (Int_t j = 0; j < 6; j++) {
01441 for (Int_t k = 0; k < 4; k++){
01442
01443 title.Form("EfficiencySuID%i_%s_%s_%s_Bin%i",usedSuId[suid],recoSetName[rese[j]].Data(),mergeName[imerge].Data(),filterName[k].Data(),irebin);
01444 newtitle.Form("Efficiency" + macroPre[macro] + filename[file] + "_%i",suid);
01445 if (hEfficiency[j*4+k] = (TProfile*)inputFile1.Get(title)){
01446 hEfficiency[j*4+k]->SetLineStyle(suid/2+1);
01447 hEfficiency[j*4+k]->SetMarkerStyle(markerS[j]);
01448 hEfficiency[j*4+k]->SetMarkerColor(1+k);
01449 hEfficiency[j*4+k]->SetLineColor(1+k);
01450 leg[k]->AddEntry(hEfficiency[j*4+k],ReSe[j],"lep");
01451 c2_2->cd(k+1)->SetLogy(1);
01452 if (j == 0)
01453 hEfficiency[j*4+k]->DrawCopy("L,P");
01454 else
01455 hEfficiency[j*4+k]->DrawCopy("L,P,same");
01456 }
01457 }
01458 }
01459 for (suid = 0; suid < 4; suid++) {
01460 c2_2->cd(suid+1);
01461 leg[suid]->Draw("same");
01462 }
01463 c2_2->SaveAs(ofolder + macroPre[macro] + "MS_with_without_Comparison_" + mergeName[imerge] + picExt);
01464
01465 for (Int_t i = 0; i < 4; i++) {
01466 leg[i]->Clear();
01467 leg[i]->SetNColumns(1);
01468 }
01469 c1_1->Clear("D");
01470 c1_2->Clear("D");
01471 c2_1->Clear("D");
01472 c2_2->Clear("D");
01473
01474 inputFile1.Close();
01475 }
01476
01477 printf("::::::::::::: Merged or notMerged Scan :::::::::::::\n");
01478
01479 }
01480
01481
01482 cout << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" << endl;
01483 suid = 0;
01484 rs = 1;
01485 imerge = 1;
01486 ifilter = 0;
01487 irebin = 0;
01488 macro = 1;
01489 for (imerge = 0; imerge < 2; imerge++) {
01490 if (dummy) {
01491
01492 printf("::::::::::::: Anode Scan :::::::::::::\n");
01493
01494 Double_t e[4] = {0.0};
01495 TString AHV[3][4] = {{"1800V","1800V","2000V","2000V"},
01496 {"1750V","1750V","1950V","1950V"},
01497 {"1700V","1700V","1900V","1900V"}
01498 };
01499 TString DV[4] = {" Drift: 600V"," Drift: 600V"," Drift: 400V"," Drift: 400V"};
01500 Double_t A[4][3] = {{1800,1750,1700},
01501 {1800,1750,1700},
01502 {2000,1950,1900},
01503 {2000,1950,1900}
01504 };
01505
01506 Float_t dx[2] = {1700,2000};
01507 Float_t dy[2] = {0,0};
01508 Float_t ddx[2] = {0,0};
01509 Float_t ddy[2] = {0,0};
01510 TGraphErrors dumm(2,dx,dy,ddx,ddy);
01511
01512 Double_t prAHV[4][3] = {{0.0}};
01513 Double_t dprAHV[4][3] = {{0.0}};
01514 TGraphErrors* gAHV[4];
01515
01516 TString AV[4] = {" Anode: 1750"," Anode: 1750"," Anode: 1950"," Anode: 1950"};
01517 TString DHV[3][4] = {{"600V","600V","400V","400V"},
01518 {"300V","300V","200V","200V"},
01519 { "60V", "60V", "40V", "40V"}
01520 };
01521 Double_t D[4][3] = {{600,300, 60},
01522 {600,300, 60},
01523 {400,200, 40},
01524 {400,200, 40}
01525 };
01526 Double_t prDHV[4][3] = {{0.0}};
01527 Double_t dprDHV[4][3] = {{0.0}};
01528 TGraphErrors* gDHV[4];
01529
01530
01531
01532
01533 TString radi[4] = {" H"," E"," F"," I(p)"};
01534
01535
01536 leg[0]->SetHeader(" 3GeV: ");
01537 leg[1]->SetHeader(" 3GeV: ");
01538
01539 appending.Form("_%s_%s_%s",recoSetName[rs].Data(),mergeName[imerge].Data(),filterName[ifilter].Data());
01540
01541
01542
01543 TFile inputFile1(ifolder + macroPre[macro] + filename[0],"READ");
01544
01545 title.Form("prResult" + appending);
01546 newtitle.Form("prResult" + macroPre[macro] + filename[0]);
01547 if (hPR[0] = (TProfile*)inputFile1.Get(title)){
01548
01549 for (suid = 0; suid < 4; suid++) {
01550 prAHV[suid][0] = hPR[0]->GetBinContent(hPR[0]->FindBin(suid));
01551 dprAHV[suid][0] = hPR[0]->GetBinError(hPR[0]->FindBin(suid));
01552
01553
01554
01555 }
01556 }
01557
01558 TFile inputFile2(ifolder + macroPre[macro] + filename[2],"READ");
01559
01560 title.Form("prResult" + appending);
01561 newtitle.Form("prResult" + macroPre[macro] + filename[2]);
01562 if (hPR[1] = (TProfile*)inputFile2.Get(title)){
01563
01564 for (suid = 0; suid < 4; suid++) {
01565 prAHV[suid][1] = hPR[1]->GetBinContent(hPR[1]->FindBin(suid));
01566 dprAHV[suid][1] = hPR[1]->GetBinError(hPR[1]->FindBin(suid));
01567
01568
01569
01570 }
01571 }
01572
01573 TFile inputFile3(ifolder + macroPre[macro] + filename[1],"READ");
01574
01575 title.Form("prResult" + appending);
01576 newtitle.Form("prResult" + macroPre[macro] + filename[1]);
01577 if(hPR[2] = (TProfile*)inputFile3.Get(title)){
01578
01579 for (suid = 0; suid < 4; suid++) {
01580 prAHV[suid][2] = hPR[2]->GetBinContent(hPR[2]->FindBin(suid));
01581 dprAHV[suid][2] = hPR[2]->GetBinError(hPR[2]->FindBin(suid));
01582
01583
01584
01585 }
01586 }
01587
01588 printf("::::::::::::: Drift Scan :::::::::::::\n");
01589
01590 TFile inputFile4(ifolder + macroPre[macro] + filename[2],"READ");
01591 title.Form("prResult" + appending);
01592 newtitle.Form("prResult" + macroPre[macro] + filename[2]);
01593 if(hPR[3] = (TProfile*)inputFile4.Get(title)){
01594
01595 for (suid = 0; suid < 4; suid++) {
01596 prDHV[suid][0] = hPR[3]->GetBinContent(hPR[3]->FindBin(suid));
01597 dprDHV[suid][0] = hPR[3]->GetBinError(hPR[3]->FindBin(suid));
01598
01599
01600
01601 }
01602 }
01603
01604 TFile inputFile5(ifolder + macroPre[macro] + filename[3],"READ");
01605 title.Form("prResult" + appending);
01606 newtitle.Form("prResult" + macroPre[macro] + filename[3]);
01607 if(hPR[4] = (TProfile*)inputFile5.Get(title)){
01608
01609 for (suid = 0; suid < 4; suid++) {
01610 prDHV[suid][1] = hPR[4]->GetBinContent(hPR[4]->FindBin(suid));
01611 dprDHV[suid][1] = hPR[4]->GetBinError(hPR[4]->FindBin(suid));
01612
01613
01614
01615 }
01616 }
01617
01618 TFile inputFile6(ifolder + macroPre[macro] + filename[4],"READ");
01619 title.Form("prResult" + appending);
01620 newtitle.Form("prResult" + macroPre[macro] + filename[4]);
01621 if(hPR[5] = (TProfile*)inputFile6.Get(title)){
01622
01623 for (suid = 0; suid < 4; suid++) {
01624 prDHV[suid][2] = hPR[5]->GetBinContent(hPR[5]->FindBin(suid));
01625 dprDHV[suid][2] = hPR[5]->GetBinError(hPR[5]->FindBin(suid));
01626
01627
01628
01629 }
01630 }
01631 c2_1->cd(1)->SetLogy(0);
01632 dumm.GetYaxis()->SetRangeUser(0,0.5);
01633 dumm.GetYaxis()->SetTitle("position resolution [mm]");
01634 dumm.GetXaxis()->SetTitle("Anode voltage [V]");
01635 dumm.Draw("ALP");
01636 for (Int_t i = 0; i < 4; i++){
01637
01638 c2_1->cd(1)->SetLogy(0);
01639 gAHV[i] = new TGraphErrors(3,A[i],prAHV[i],e,dprAHV[i]);
01640 leg[0]->AddEntry(gAHV[i],usedSuName[i] + DV[i] + radi[i],"lep");
01641 gAHV[i]->SetLineStyle(i/2+1);
01642 gAHV[i]->SetMarkerStyle(20+i);
01643 gAHV[i]->SetMarkerColor(1+i);
01644 gAHV[i]->SetLineColor(1+i);
01645 gAHV[i]->GetYaxis()->SetRangeUser(0,0.5);
01646 gAHV[i]->GetYaxis()->SetTitle("position resolution [mm]");
01647 gAHV[i]->GetXaxis()->SetRangeUser(1500,2200);
01648 gAHV[i]->GetXaxis()->SetTitle("Anode voltage [V]");
01649 if (i == 0)
01650 gAHV[i]->Draw("same,LP");
01651 else
01652 gAHV[i]->Draw("same,LP");
01653 leg[0]->Draw("same");
01654 c2_1->cd(2)->SetLogy(0);
01655 gDHV[i] = new TGraphErrors(3,D[i],prDHV[i],e,dprDHV[i]);
01656 leg[1]->AddEntry(gDHV[i],usedSuName[i] + AV[i] + radi[i],"lep");
01657 gDHV[i]->SetLineStyle(i/2+1);
01658 gDHV[i]->SetMarkerStyle(24+i);
01659 gDHV[i]->SetMarkerColor(1+i);
01660 gDHV[i]->SetLineColor(1+i);
01661 gDHV[i]->GetYaxis()->SetRangeUser(0,0.5);
01662 gDHV[i]->GetYaxis()->SetTitle("position resolution [mm]");
01663 gDHV[i]->GetXaxis()->SetRangeUser(0,650);
01664 gDHV[i]->GetXaxis()->SetTitle("Drift voltage [V]");
01665 if (i == 0)
01666 gDHV[i]->Draw("ALP");
01667 else
01668 gDHV[i]->Draw("same,LP");
01669 leg[1]->Draw("same");
01670 if (i == 3)
01671 c2_1->SaveAs(ofolder + macroPre[macro] + "MS_HVScan_Comparison_" + mergeName[imerge] + picExt);
01672 }
01673
01674 c1_1->Clear("D");
01675 c1_2->Clear("D");
01676 c2_1->Clear("D");
01677 c2_2->Clear("D");
01678
01679 inputFile1.Close();
01680 inputFile2.Close();
01681 inputFile3.Close();
01682 inputFile4.Close();
01683 inputFile5.Close();
01684 inputFile6.Close();
01685
01686 for (Int_t i = 0; i < 4; i++)
01687 leg[i]->Clear();
01688 leg[0]->SetNColumns(1);
01689 leg[1]->SetNColumns(1);
01690 }
01691
01692 if (dummy) {
01693 printf("::::::::::::: Momentum Scan :::::::::::::\n");
01694 TString radi[4] = {" H:"," H(1/3):"," G(30):"," B:"};
01695 Double_t M[5] = {2,4,6,8,10};
01696 Double_t dM[5] = {0.0};
01697 Double_t prM[4][5] = {{0.0}};
01698 Double_t dprM[4][5] = {{0.0}};
01699 TFile *inputFile[5];
01700 TGraphErrors *gM[5];
01701 for (Int_t file = 0; file < 5; file++) {
01702 inputFile[file] = new TFile(ifolder + macroPre[macro] + filename[file+5],"READ");
01703 title.Form("prResult" + appending);
01704 newtitle.Form("prResult" + macroPre[macro] + filename[file+5]);
01705 if(hPR[file] = (TProfile*)inputFile[file]->Get(title)){
01706 for (suid = 0; suid < 4; suid++) {
01707 prM[suid][file] = hPR[file]->GetBinContent(hPR[file]->FindBin(suid));
01708 dprM[suid][file] = hPR[file]->GetBinError(hPR[file]->FindBin(suid));
01709
01710
01711 }
01712 }
01713 }
01714 c1_1->cd(1)->SetLogy(0);
01715 for (suid = 0; suid < 4; suid++) {
01716 gM[suid] = new TGraphErrors(5,M,prM[suid],dM,dprM[suid]);
01717 leg[0]->AddEntry(gM[suid],usedSuName[suid] + radi[suid],"lep");
01718 gM[suid]->SetLineStyle(suid/2+1);
01719 gM[suid]->SetMarkerStyle(24+suid);
01720 gM[suid]->SetMarkerColor(1+suid);
01721 gM[suid]->SetLineColor(1+suid);
01722 gM[suid]->GetYaxis()->SetRangeUser(0,0.35);
01723 gM[suid]->GetYaxis()->SetTitle("position resolution [mm]");
01724 gM[suid]->GetXaxis()->SetRangeUser(0,12);
01725 gM[suid]->GetXaxis()->SetTitle("Momentum [GeV/c]");
01726 c1_1->cd(1)->SetLogy(0);
01727 if (suid == 0)
01728 gM[suid]->Draw("ALP");
01729 else
01730 gM[suid]->Draw("same,LP");
01731 if (suid == 3){
01732 leg[0]->Draw("same");
01733 c1_1->SaveAs(ofolder + macroPre[macro] + "MS_MomentumScan_Comparison_" + mergeName[imerge] + picExt);
01734 }
01735 }
01736
01737 for (Int_t i = 0; i < 4; i++) {
01738 leg[i]->Clear();
01739 leg[i]->SetNColumns(1);
01740 }
01741 for (Int_t file = 0; file < 5; file++)
01742 inputFile[file]->Close();
01743
01744
01745
01746
01747 }
01748 }
01749 if (dummy) {
01750
01751
01752 printf("::::::::::::: Signal Scan :::::::::::::\n");
01753 if (dummy) {
01754 Int_t file = 11;
01755
01756 leg[0]->SetNColumns(1);
01757 leg[0]->SetHeader(" G(10), 3GeV:");
01758 TProfile dumm("dummy","dummy",12,-0.5,11.5);
01759 dumm.GetYaxis()->SetRangeUser(0,0.35);
01760 Double_t S[12] = {0,1,2,3,4,5,6,7,8,9,10,11};
01761 Double_t dS[12] = {0};
01762 Double_t prS[4][2][12] = {{{0.0}}};
01763 Double_t dprS[4][2][12] = {{{0.0}}};
01764 TGraphErrors *gS[4][2];
01765 TFile inputFile1(ifolder + macroPre[macro] + filename[file],"READ");
01766 TString rad[4] = {"10 fiber mats","10 fiber mats","10 fiber mats","10 fiber mats"};
01767 Int_t rese[3] = {1,4,7};
01768 TString ReSe[3] = {"f. t. int. ","t. w. int. ","ampl. "};
01769 TString FiSe[4] = {"(wOwU)","(woOwU)","(woOwoU)","(wOwoU)"};
01770 TString MeSe[2] = {" not merged"," merged"};
01771 Int_t markerS[3] = {20,24,21};
01772 for (Int_t m = 0; m < 2; m++) {
01773 for (Int_t f = 0; f < 4; f++) {
01774 for (Int_t r = 0; r < 3; r++) {
01775 appending.Form("_%s_%s_%s",recoSetName[rese[r]].Data(),mergeName[m].Data(),filterName[f].Data());
01776 for (suid = 0; suid < 4; suid++) {
01777 title.Form("prwResult" + appending);
01778 newtitle.Form("prwResult" + macroPre[macro] + filename[file] + "_%i",suid);
01779 if(hPR[0] = (TProfile*)inputFile1.Get(title)){
01780 prS[suid][m][f*3+r] = hPR[0]->GetBinContent(hPR[0]->FindBin(suid));
01781 dprS[suid][m][f*3+r] = hPR[0]->GetBinError(hPR[0]->FindBin(suid));
01782 }
01783 }
01784 }
01785 }
01786 for (suid = 0; suid < 4; suid++) {
01787 gS[suid][m] = new TGraphErrors(12,S,prS[suid][m],dS,dprS[suid][m]);
01788 for (Int_t f = 0; f < 4; f++) {
01789 for (Int_t r = 0; r < 3; r++) {
01790
01791 dumm.GetXaxis()->SetBinLabel(f*3+r+1, ReSe[r] + FiSe[f]);
01792 }
01793 }
01794
01795 leg[0]->AddEntry(gS[suid][m],usedSuName[suid] + MeSe[m],"lep");
01796 gS[suid][m]->SetLineStyle(suid/2+1);
01797 if (m == 0)
01798 gS[suid][m]->SetMarkerStyle(24+suid);
01799 else
01800 gS[suid][m]->SetMarkerStyle(20+suid);
01801 gS[suid][m]->SetMarkerColor(1+suid);
01802 gS[suid][m]->SetLineColor(1+suid);
01803 gS[suid][m]->GetYaxis()->SetRangeUser(0,0.35);
01804 gS[suid][m]->GetYaxis()->SetTitle("position resolution [mm]");
01805 c1_1->cd(1)->SetLogy(0);
01806 if (suid == 0 && m == 0) {
01807 dumm.DrawCopy();
01808 gS[suid][m]->Draw("same,LP");
01809
01810 }
01811 else
01812 gS[suid][m]->Draw("same,LP");
01813 if (suid == 3){
01814 leg[0]->Draw("same");
01815
01816 }
01817 }
01818 }
01819 c1_1->SaveAs(ofolder + macroPre[macro] + "MS_Signal_Comparison_" + picExt);
01820
01821 for (Int_t i = 0; i < 4; i++) {
01822 leg[i]->Clear();
01823 leg[i]->SetNColumns(1);
01824 }
01825 c1_1->Clear("D");
01826 c1_2->Clear("D");
01827 c2_1->Clear("D");
01828 c2_2->Clear("D");
01829
01830 inputFile1.Close();
01831 }
01832
01833 }
01834
01835 c1_1->Close();
01836 c1_2->Close();
01837 c2_1->Close();
01838 c2_2->Close();
01839
01840
01841
01842 }
