File indexing completed on 2024-04-06 12:33:30
0001 {
0002
0003 gROOT->Reset();
0004 gROOT->Clear();
0005
0006 gStyle->SetNdivisions(5);
0007 gStyle->SetCanvasBorderMode(0);
0008 gStyle->SetPadBorderMode(1);
0009 gStyle->SetOptTitle(1);
0010 gStyle->SetStatFont(42);
0011 gStyle->SetCanvasColor(10);
0012 gStyle->SetPadColor(0);
0013 gStyle->SetTitleFont(62,"xy");
0014 gStyle->SetLabelFont(62,"xy");
0015 gStyle->SetTitleFontSize(0.07);
0016 gStyle->SetTitleSize(0.045,"xy");
0017 gStyle->SetLabelSize(0.05,"xy");
0018 gStyle->SetHistFillStyle(1001);
0019 gStyle->SetHistFillColor(0);
0020 gStyle->SetHistLineStyle(1);
0021 gStyle->SetHistLineWidth(1);
0022 gStyle->SetHistLineColor(1);
0023 gStyle->SetTitleXOffset(1.1);
0024 gStyle->SetTitleYOffset(1.15);
0025
0026 gStyle->SetOptStat(kFALSE);
0027 gStyle->SetOptFit(0111);
0028 gStyle->SetStatH(0.1);
0029
0030
0031
0032
0033 const char * samples="SCENARIOS";
0034 TString samplestring(samples);
0035 std::cout<<"Going to analyze the following scenarios:"<<samplestring<<endl;
0036 TString delim(" ");
0037 TObjArray *samplelist=samplestring.Tokenize(delim);
0038 int nscenarios=samplelist->GetEntries();
0039 TFile *files[nscenarios];
0040 TDirectory *dir[nscenarios];
0041 TH1F *fakeratevseta[nscenarios], *fakeratevsPt[nscenarios], *fakeratevshits[nscenarios];
0042 TH1F *efficvseta[nscenarios], *efficvsPt[nscenarios], *efficvshits[nscenarios];
0043 TH1F *sigmadxyvseta[nscenarios], *sigmadxyvspt[nscenarios], *pullDxy[nscenarios];
0044 TH1F *sigmadzvseta[nscenarios], *sigmadzvspt[nscenarios], *pullDz[nscenarios];
0045 TH1F *sigmaptvseta[nscenarios], *sigmaptvspt[nscenarios], *pullPt[nscenarios];
0046 for(int i=0;i<nscenarios;i++){
0047 TObjString * string =(TObjString*)samplelist->At(i);
0048 files[i]=new TFile("val.SAMPLE."+string->GetString()+".root");
0049 files[i]->cd("DQMData/RecoTrackV/Track");
0050 TDirectory * dir[i]=gDirectory;
0051 TList *l= dir[i]->GetListOfKeys();
0052 string rdirName=l->At(0)->GetName();
0053
0054 dir[i]->GetObject((rdirName+"/fakerate").c_str(),fakeratevseta[i]);
0055 dir[i]->GetObject((rdirName+"/fakeratePt").c_str(),fakeratevsPt[i]);
0056 dir[i]->GetObject((rdirName+"/fakerate_vs_hit").c_str(),fakeratevshits[i]);
0057
0058 dir[i]->GetObject((rdirName+"/effic").c_str(),efficvseta[i]);
0059 dir[i]->GetObject((rdirName+"/efficPt").c_str(),efficvsPt[i]);
0060 dir[i]->GetObject((rdirName+"/effic_vs_hit").c_str(),efficvshits[i]);
0061
0062 dir[i]->GetObject((rdirName+"/sigmadxy").c_str(),sigmadxyvseta[i]);
0063 dir[i]->GetObject((rdirName+"/sigmadxyPt").c_str(),sigmadxyvspt[i]);
0064 dir[i]->GetObject((rdirName+"/pullDxy").c_str(),pullDxy[i]);
0065
0066 dir[i]->GetObject((rdirName+"/sigmadz").c_str(),sigmadzvseta[i]);
0067 dir[i]->GetObject((rdirName+"/sigmadzPt").c_str(),sigmadzvspt[i]);
0068 dir[i]->GetObject((rdirName+"/pullDz").c_str(),pullDz[i]);
0069
0070 dir[i]->GetObject((rdirName+"/sigmapt").c_str(),sigmaptvseta[i]);
0071 dir[i]->GetObject((rdirName+"/sigmaptPt").c_str(),sigmaptvspt[i]);
0072 dir[i]->GetObject((rdirName+"/pullPt").c_str(),pullPt[i]);
0073 }
0074
0075
0076
0077
0078
0079 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0080
0081
0082
0083 c1->SetGrid(1,1);
0084 for(int i=0;i<nscenarios;i++){
0085
0086 fakeratevseta[i]->SetTitle("Fake rate vs #eta for SAMPLE events");
0087 fakeratevseta[i]->SetMarkerStyle(20+i);
0088 fakeratevseta[i]->SetMarkerColor(2+i);
0089 fakeratevseta[i]->SetMarkerSize(0.9);
0090 fakeratevseta[i]->SetLineColor(1);
0091 fakeratevseta[i]->SetLineWidth(1);
0092 }
0093
0094 TLegend *leg1 = new TLegend(0.1,0.75,0.35,0.9);
0095 leg1->SetTextAlign(32);
0096 leg1->SetTextColor(1);
0097 leg1->SetTextSize(0.04);
0098
0099 for(int i=0;i<nscenarios;i++){
0100 fakeratevseta[i]->SetTitle("Fake rate vs #eta for SAMPLE events");
0101 fakeratevseta[i]->SetXTitle("#eta");
0102 fakeratevseta[i]->SetYTitle("Fake Rate");
0103 fakeratevseta[i]->SetAxisRange(-0.005,0.4,"Y");
0104 if(i==0)fakeratevseta[i]->Draw();
0105 else fakeratevseta[i]->Draw("same");
0106 leg1->AddEntry(fakeratevseta[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0107 }
0108
0109 leg1->Draw();
0110
0111 c1->Update();
0112 c1->SaveAs("Fake_eta_SAMPLE.eps");
0113 c1->SaveAs("Fake_eta_SAMPLE.gif");
0114 c1->WaitPrimitive();
0115
0116 delete c1;
0117 delete leg1;
0118
0119
0120
0121
0122 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0123
0124
0125
0126 c1->SetGrid(1,1);
0127 for(int i=0;i<nscenarios;i++){
0128
0129 fakeratevsPt[i]->SetTitle("Fake rate vs Pt for SAMPLE events");
0130 fakeratevsPt[i]->SetMarkerStyle(20+i);
0131 fakeratevsPt[i]->SetMarkerColor(2+i);
0132 fakeratevsPt[i]->SetMarkerSize(0.9);
0133 fakeratevsPt[i]->SetLineColor(1);
0134 fakeratevsPt[i]->SetLineWidth(1);
0135 }
0136
0137 TLegend *leg1 = new TLegend(0.1,0.75,0.35,0.9);
0138 leg1->SetTextAlign(32);
0139 leg1->SetTextColor(1);
0140 leg1->SetTextSize(0.04);
0141
0142 for(int i=0;i<nscenarios;i++){
0143 fakeratevsPt[i]->SetTitle("Fake rate vs Pt for SAMPLE events");
0144 fakeratevsPt[i]->SetXTitle("Pt");
0145 fakeratevsPt[i]->SetYTitle("Fake Rate");
0146 fakeratevsPt[i]->SetAxisRange(-0.005,0.5,"Y");
0147 if(i==0)fakeratevsPt[i]->Draw();
0148 else fakeratevsPt[i]->Draw("same");
0149 leg1->AddEntry(fakeratevsPt[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0150 }
0151
0152 leg1->Draw();
0153
0154 c1->Update();
0155 c1->SaveAs("Fake_Pt_SAMPLE.eps");
0156 c1->SaveAs("Fake_Pt_SAMPLE.gif");
0157 c1->WaitPrimitive();
0158
0159 delete c1;
0160 delete leg1;
0161
0162
0163
0164
0165 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0166
0167
0168
0169 c1->SetGrid(1,1);
0170 for(int i=0;i<nscenarios;i++){
0171
0172 fakeratevshits[i]->SetTitle("Fake rate vs hits for SAMPLE events");
0173 fakeratevshits[i]->SetMarkerStyle(20+i);
0174 fakeratevshits[i]->SetMarkerColor(2+i);
0175 fakeratevshits[i]->SetMarkerSize(0.9);
0176 fakeratevshits[i]->SetLineColor(1);
0177 fakeratevshits[i]->SetLineWidth(1);
0178 }
0179
0180 TLegend *leg1 = new TLegend(0.65,0.75,0.9,0.9);
0181 leg1->SetTextAlign(32);
0182 leg1->SetTextColor(1);
0183 leg1->SetTextSize(0.04);
0184
0185 for(int i=0;i<nscenarios;i++){
0186 fakeratevshits[i]->SetTitle("Fake rate vs hits for SAMPLE events");
0187 fakeratevshits[i]->SetXTitle("hits");
0188 fakeratevshits[i]->SetYTitle("Fake Rate");
0189 fakeratevshits[i]->SetAxisRange(-0.005,0.7,"Y");
0190 if(i==0)fakeratevshits[i]->Draw();
0191 else fakeratevshits[i]->Draw("same");
0192 leg1->AddEntry(fakeratevshits[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0193 }
0194
0195 leg1->Draw();
0196
0197 c1->Update();
0198 c1->SaveAs("Fake_hits_SAMPLE.eps");
0199 c1->SaveAs("Fake_hits_SAMPLE.gif");
0200 c1->WaitPrimitive();
0201
0202 delete c1;
0203 delete leg1;
0204
0205
0206
0207
0208 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0209
0210
0211
0212 c1->SetGrid(1,1);
0213 for(int i=0;i<nscenarios;i++){
0214
0215 efficvseta[i]->SetTitle("Efficiency vs #eta for SAMPLE events");
0216 efficvseta[i]->SetMarkerStyle(20+i);
0217 efficvseta[i]->SetMarkerColor(2+i);
0218 efficvseta[i]->SetMarkerSize(0.9);
0219 efficvseta[i]->SetLineColor(1);
0220 efficvseta[i]->SetLineWidth(1);
0221 }
0222
0223 TLegend *leg1 = new TLegend(0.1,0.1,0.35,0.35);
0224 leg1->SetTextAlign(32);
0225 leg1->SetTextColor(1);
0226 leg1->SetTextSize(0.04);
0227
0228 for(int i=0;i<nscenarios;i++){
0229 efficvseta[i]->SetXTitle("#eta");
0230 efficvseta[i]->SetYTitle("Efficiency");
0231 efficvseta[i]->SetAxisRange(-0.005,1.,"Y");
0232 if(i==0)efficvseta[i]->Draw();
0233 else efficvseta[i]->Draw("same");
0234 leg1->AddEntry(efficvseta[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0235 }
0236
0237 leg1->Draw();
0238
0239 c1->Update();
0240 c1->SaveAs("Efficiency_eta_SAMPLE.eps");
0241 c1->SaveAs("Efficiency_eta_SAMPLE.gif");
0242 c1->WaitPrimitive();
0243
0244 delete c1;
0245 delete leg1;
0246
0247
0248
0249
0250 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0251
0252
0253
0254 c1->SetGrid(1,1);
0255 for(int i=0;i<nscenarios;i++){
0256
0257 efficvsPt[i]->SetTitle("Efficiency vs Pt for SAMPLE events");
0258 efficvsPt[i]->SetMarkerStyle(20+i);
0259 efficvsPt[i]->SetMarkerColor(2+i);
0260 efficvsPt[i]->SetMarkerSize(0.9);
0261 efficvsPt[i]->SetLineColor(1);
0262 efficvsPt[i]->SetLineWidth(1);
0263 }
0264
0265 TLegend *leg1 = new TLegend(0.2,0.1,0.45,0.35);
0266 leg1->SetTextAlign(32);
0267 leg1->SetTextColor(1);
0268 leg1->SetTextSize(0.04);
0269
0270 for(int i=0;i<nscenarios;i++){
0271 efficvsPt[i]->SetXTitle("Pt");
0272 efficvsPt[i]->SetYTitle("Efficiency");
0273 efficvsPt[i]->SetAxisRange(-0.005,1.,"Y");
0274 if(i==0)efficvsPt[i]->Draw();
0275 else efficvsPt[i]->Draw("same");
0276 leg1->AddEntry(efficvsPt[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0277 }
0278
0279 leg1->Draw();
0280
0281 c1->Update();
0282 c1->SaveAs("Efficiency_Pt_SAMPLE.eps");
0283 c1->SaveAs("Efficiency_Pt_SAMPLE.gif");
0284 c1->WaitPrimitive();
0285
0286 delete c1;
0287 delete leg1;
0288
0289
0290
0291
0292 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0293
0294
0295
0296 c1->SetGrid(1,1);
0297 for(int i=0;i<nscenarios;i++){
0298
0299 efficvshits[i]->SetTitle("Efficiency vs hits for SAMPLE events");
0300 efficvshits[i]->SetMarkerStyle(20+i);
0301 efficvshits[i]->SetMarkerColor(2+i);
0302 efficvshits[i]->SetMarkerSize(0.9);
0303 efficvshits[i]->SetLineColor(1);
0304 efficvshits[i]->SetLineWidth(1);
0305 }
0306
0307
0308
0309
0310
0311 TLegend *leg1 = new TLegend(0.2,0.1,0.45,0.35);
0312 leg1->SetTextAlign(32);
0313 leg1->SetTextColor(1);
0314 leg1->SetTextSize(0.04);
0315
0316 for(int i=0;i<nscenarios;i++){
0317 efficvshits[i]->SetXTitle("hits");
0318 efficvshits[i]->SetYTitle("Efficiency");
0319 efficvshits[i]->SetAxisRange(0.,1.,"Y");
0320 if(i==0)efficvshits[i]->Draw();
0321 else efficvshits[i]->Draw("same");
0322 leg1->AddEntry(efficvshits[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0323 }
0324
0325 leg1->Draw();
0326
0327 c1->Update();
0328 c1->SaveAs("Efficiency_hits_SAMPLE.eps");
0329 c1->SaveAs("Efficiency_hits_SAMPLE.gif");
0330 c1->WaitPrimitive();
0331
0332 delete c1;
0333 delete leg1;
0334
0335
0336
0337
0338 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0339
0340
0341
0342 c1->SetGrid(1,1);
0343 for(int i=0;i<nscenarios;i++){
0344
0345 sigmadxyvseta[i]->SetTitle("Dxy resolution vs #eta for SAMPLE events");
0346 sigmadxyvseta[i]->SetMarkerStyle(20+i);
0347 sigmadxyvseta[i]->SetMarkerColor(2+i);
0348 sigmadxyvseta[i]->SetMarkerSize(0.9);
0349 sigmadxyvseta[i]->SetLineColor(1);
0350 sigmadxyvseta[i]->SetLineWidth(1);
0351 }
0352
0353
0354 TLegend *leg1 = new TLegend(0.2,0.1,0.45,0.35);
0355 leg1->SetTextAlign(32);
0356 leg1->SetTextColor(1);
0357 leg1->SetTextSize(0.04);
0358
0359 for(int i=0;i<nscenarios;i++){
0360 sigmadxyvseta[i]->SetXTitle("#eta");
0361 sigmadxyvseta[i]->SetYTitle("dxy resolution");
0362 sigmadxyvseta[i]->SetAxisRange(0.0005,0.1,"Y");
0363 if(i==0)sigmadxyvseta[i]->Draw();
0364 else sigmadxyvseta[i]->Draw("same");
0365 leg1->AddEntry(sigmadxyvseta[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0366 }
0367 c1->SetLogy();
0368 leg1->Draw();
0369
0370 c1->Update();
0371 c1->SaveAs("Dxy_res_eta_SAMPLE.eps");
0372 c1->SaveAs("Dxy_res_eta_SAMPLE.gif");
0373 c1->WaitPrimitive();
0374
0375 delete c1;
0376 delete leg1;
0377
0378
0379
0380
0381 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0382
0383
0384
0385 c1->SetGrid(1,1);
0386 for(int i=0;i<nscenarios;i++){
0387
0388 sigmadxyvspt[i]->SetTitle("Dxy resolution vs Pt for SAMPLE events");
0389 sigmadxyvspt[i]->SetMarkerStyle(20+i);
0390 sigmadxyvspt[i]->SetMarkerColor(2+i);
0391 sigmadxyvspt[i]->SetMarkerSize(0.9);
0392 sigmadxyvspt[i]->SetLineColor(1);
0393 sigmadxyvspt[i]->SetLineWidth(1);
0394 }
0395
0396
0397 TLegend *leg1 = new TLegend(0.1,0.1,0.35,0.35);
0398 leg1->SetTextAlign(32);
0399 leg1->SetTextColor(1);
0400 leg1->SetTextSize(0.04);
0401
0402 for(int i=0;i<nscenarios;i++){
0403 sigmadxyvspt[i]->SetXTitle("Pt");
0404 sigmadxyvspt[i]->SetYTitle("dxy resolution");
0405 sigmadxyvspt[i]->SetAxisRange(0.,30,"X");
0406 sigmadxyvspt[i]->SetAxisRange(0.0005,0.1,"Y");
0407 if(i==0)sigmadxyvspt[i]->Draw();
0408 else sigmadxyvspt[i]->Draw("same");
0409 leg1->AddEntry(sigmadxyvspt[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0410 }
0411
0412 leg1->Draw();
0413 c1->SetLogy();
0414 c1->Update();
0415 c1->SaveAs("Dxy_res_Pt_SAMPLE.eps");
0416 c1->SaveAs("Dxy_res_Pt_SAMPLE.gif");
0417 c1->WaitPrimitive();
0418
0419 delete c1;
0420 delete leg1;
0421
0422
0423
0424
0425 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0426
0427
0428
0429 c1->SetGrid(1,1);
0430 for(int i=0;i<nscenarios;i++){
0431
0432 pullDxy[i]->SetTitle("Pull dxy for SAMPLE events");
0433 pullDxy[i]->SetMarkerStyle(20+i);
0434 pullDxy[i]->SetMarkerColor(2+i);
0435 pullDxy[i]->SetMarkerSize(0.9);
0436 pullDxy[i]->SetLineColor(2+i);
0437 pullDxy[i]->SetLineWidth(1);
0438 }
0439
0440 TLegend *leg1 = new TLegend(0.1,0.6,0.35,0.85);
0441 leg1->SetTextAlign(32);
0442 leg1->SetTextColor(1);
0443 leg1->SetTextSize(0.04);
0444
0445 for(int i=0;i<nscenarios;i++){
0446 pullDxy[i]->SetXTitle("Pulldxy");
0447
0448
0449 if(i==0)pullDxy[i]->Draw();
0450 else pullDxy[i]->Draw("same");
0451 leg1->AddEntry(pullDxy[i],((TObjString *)(samplelist->At(i)))->GetString(),"l");
0452 }
0453
0454 leg1->Draw();
0455 c1->SetLogy();
0456 c1->Update();
0457 c1->SaveAs("Pulldxy_SAMPLE.eps");
0458 c1->SaveAs("Pulldxy_SAMPLE.gif");
0459 c1->WaitPrimitive();
0460
0461 delete c1;
0462 delete leg1;
0463
0464
0465
0466
0467 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0468
0469
0470
0471 c1->SetGrid(1,1);
0472 for(int i=0;i<nscenarios;i++){
0473
0474 sigmadzvseta[i]->SetTitle("Dz resolution vs #eta for SAMPLE events");
0475 sigmadzvseta[i]->SetMarkerStyle(20+i);
0476 sigmadzvseta[i]->SetMarkerColor(2+i);
0477 sigmadzvseta[i]->SetMarkerSize(0.9);
0478 sigmadzvseta[i]->SetLineColor(1);
0479 sigmadzvseta[i]->SetLineWidth(1);
0480 }
0481
0482
0483 TLegend *leg1 = new TLegend(0.2,0.1,0.45,0.35);
0484 leg1->SetTextAlign(32);
0485 leg1->SetTextColor(1);
0486 leg1->SetTextSize(0.04);
0487
0488 for(int i=0;i<nscenarios;i++){
0489 sigmadzvseta[i]->SetXTitle("#eta");
0490 sigmadzvseta[i]->SetYTitle("dz resolution");
0491 sigmadzvseta[i]->SetAxisRange(0.0005,0.5,"Y");
0492 if(i==0)sigmadzvseta[i]->Draw();
0493 else sigmadzvseta[i]->Draw("same");
0494 leg1->AddEntry(sigmadzvseta[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0495 }
0496
0497 leg1->Draw();
0498 c1->SetLogy();
0499
0500 c1->Update();
0501 c1->SaveAs("Dz_res_eta_SAMPLE.eps");
0502 c1->SaveAs("Dz_res_eta_SAMPLE.gif");
0503 c1->WaitPrimitive();
0504
0505 delete c1;
0506 delete leg1;
0507
0508
0509
0510
0511 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0512
0513
0514
0515 c1->SetGrid(1,1);
0516 for(int i=0;i<nscenarios;i++){
0517
0518 sigmadzvspt[i]->SetTitle("Dz resolution vs Pt for SAMPLE events");
0519 sigmadzvspt[i]->SetMarkerStyle(20+i);
0520 sigmadzvspt[i]->SetMarkerColor(2+i);
0521 sigmadzvspt[i]->SetMarkerSize(0.9);
0522 sigmadzvspt[i]->SetLineColor(1);
0523 sigmadzvspt[i]->SetLineWidth(1);
0524 }
0525
0526
0527 TLegend *leg1 = new TLegend(0.2,0.1,0.45,0.35);
0528 leg1->SetTextAlign(32);
0529 leg1->SetTextColor(1);
0530 leg1->SetTextSize(0.04);
0531
0532 for(int i=0;i<nscenarios;i++){
0533 sigmadzvspt[i]->SetXTitle("Pt");
0534 sigmadzvspt[i]->SetYTitle("dz resolution");
0535 sigmadzvspt[i]->SetAxisRange(0,30.,"X");
0536 sigmadzvspt[i]->SetAxisRange(0.0005,0.5,"Y");
0537 if(i==0)sigmadzvspt[i]->Draw();
0538 else sigmadzvspt[i]->Draw("same");
0539 leg1->AddEntry(sigmadzvspt[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0540 }
0541
0542 leg1->Draw();
0543 c1->SetLogy();
0544
0545 c1->Update();
0546 c1->SaveAs("Dz_res_Pt_SAMPLE.eps");
0547 c1->SaveAs("Dz_res_Pt_SAMPLE.gif");
0548 c1->WaitPrimitive();
0549
0550 delete c1;
0551 delete leg1;
0552
0553
0554
0555
0556 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0557
0558
0559
0560 c1->SetGrid(1,1);
0561 for(int i=0;i<nscenarios;i++){
0562
0563 pullDz[i]->SetTitle("Pull dz for SAMPLE events");
0564 pullDz[i]->SetMarkerStyle(20+i);
0565 pullDz[i]->SetMarkerColor(2+i);
0566 pullDz[i]->SetMarkerSize(0.9);
0567 pullDz[i]->SetLineColor(2+i);
0568 pullDz[i]->SetLineWidth(1);
0569 }
0570
0571 TLegend *leg1 = new TLegend(0.1,0.6,0.35,0.85);
0572 leg1->SetTextAlign(32);
0573 leg1->SetTextColor(1);
0574 leg1->SetTextSize(0.04);
0575
0576 for(int i=0;i<nscenarios;i++){
0577 pullDz[i]->SetXTitle("Pulldz");
0578
0579 if(i==0)pullDz[i]->Draw();
0580 else pullDz[i]->Draw("same");
0581 leg1->AddEntry(pullDz[i],((TObjString *)(samplelist->At(i)))->GetString(),"l");
0582 }
0583
0584 leg1->Draw();
0585 c1->SetLogy();
0586
0587 c1->Update();
0588 c1->SaveAs("Pulldz_SAMPLE.eps");
0589 c1->SaveAs("Pulldz_SAMPLE.gif");
0590 c1->WaitPrimitive();
0591
0592 delete c1;
0593 delete leg1;
0594
0595
0596
0597
0598 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0599
0600
0601
0602 c1->SetGrid(1,1);
0603 for(int i=0;i<nscenarios;i++){
0604
0605 sigmaptvseta[i]->SetTitle("Pt resolution vs #eta for SAMPLE events");
0606 sigmaptvseta[i]->SetMarkerStyle(20+i);
0607 sigmaptvseta[i]->SetMarkerColor(2+i);
0608 sigmaptvseta[i]->SetMarkerSize(0.9);
0609 sigmaptvseta[i]->SetLineColor(1);
0610 sigmaptvseta[i]->SetLineWidth(1);
0611 }
0612
0613
0614 TLegend *leg1 = new TLegend(0.2,0.1,0.45,0.35);
0615 leg1->SetTextAlign(32);
0616 leg1->SetTextColor(1);
0617 leg1->SetTextSize(0.04);
0618
0619 for(int i=0;i<nscenarios;i++){
0620 sigmaptvseta[i]->SetXTitle("#eta");
0621 sigmaptvseta[i]->SetYTitle("pt resolution");
0622 sigmaptvseta[i]->SetAxisRange(0.0005,0.1,"Y");
0623 if(i==0)sigmaptvseta[i]->Draw();
0624 else sigmaptvseta[i]->Draw("same");
0625 leg1->AddEntry(sigmaptvseta[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0626 }
0627
0628 leg1->Draw();
0629 c1->SetLogy();
0630
0631 c1->Update();
0632 c1->SaveAs("Pt_res_eta_SAMPLE.eps");
0633 c1->SaveAs("Pt_res_eta_SAMPLE.gif");
0634 c1->WaitPrimitive();
0635
0636 delete c1;
0637 delete leg1;
0638
0639
0640
0641
0642 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0643
0644
0645
0646 c1->SetGrid(1,1);
0647 for(int i=0;i<nscenarios;i++){
0648
0649 sigmaptvspt[i]->SetTitle("Pt resolution vs Pt for SAMPLE events");
0650 sigmaptvspt[i]->SetMarkerStyle(20+i);
0651 sigmaptvspt[i]->SetMarkerColor(2+i);
0652 sigmaptvspt[i]->SetMarkerSize(0.9);
0653 sigmaptvspt[i]->SetLineColor(1);
0654 sigmaptvspt[i]->SetLineWidth(1);
0655 }
0656
0657
0658 TLegend *leg1 = new TLegend(0.2,0.1,0.45,0.35);
0659 leg1->SetTextAlign(32);
0660 leg1->SetTextColor(1);
0661 leg1->SetTextSize(0.04);
0662
0663 for(int i=0;i<nscenarios;i++){
0664 sigmaptvspt[i]->SetXTitle("Pt");
0665 sigmaptvspt[i]->SetYTitle("pt resolution");
0666 sigmaptvspt[i]->SetAxisRange(0.,30.,"X");
0667 sigmaptvspt[i]->SetAxisRange(0.0005,0.1,"Y");
0668 if(i==0)sigmaptvspt[i]->Draw();
0669 else sigmaptvspt[i]->Draw("same");
0670 leg1->AddEntry(sigmaptvspt[i],((TObjString *)(samplelist->At(i)))->GetString(),"P");
0671 }
0672
0673 leg1->Draw();
0674 c1->SetLogy();
0675
0676 c1->Update();
0677 c1->SaveAs("Pt_res_Pt_SAMPLE.eps");
0678 c1->SaveAs("Pt_res_Pt_SAMPLE.gif");
0679 c1->WaitPrimitive();
0680
0681 delete c1;
0682 delete leg1;
0683
0684
0685
0686
0687 TCanvas *c1 = new TCanvas("c1", "c1",129,17,926,703);
0688
0689
0690
0691 c1->SetGrid(1,1);
0692 for(int i=0;i<nscenarios;i++){
0693
0694 pullPt[i]->SetTitle("Pull pt for SAMPLE events");
0695 pullPt[i]->SetMarkerStyle(20+i);
0696 pullPt[i]->SetMarkerColor(2+i);
0697 pullPt[i]->SetMarkerSize(0.9);
0698 pullPt[i]->SetLineColor(2+i);
0699 pullPt[i]->SetLineWidth(1);
0700 }
0701
0702 TLegend *leg1 = new TLegend(0.1,0.6,0.35,0.85);
0703 leg1->SetTextAlign(32);
0704 leg1->SetTextColor(1);
0705 leg1->SetTextSize(0.04);
0706
0707 for(int i=0;i<nscenarios;i++){
0708 pullPt[i]->SetXTitle("Pullpt");
0709
0710 if(i==0)pullPt[i]->Draw();
0711 else pullPt[i]->Draw("same");
0712 leg1->AddEntry(pullPt[i],((TObjString *)(samplelist->At(i)))->GetString(),"l");
0713 }
0714
0715 leg1->Draw();
0716 c1->SetLogy();
0717
0718 c1->Update();
0719 c1->SaveAs("Pullpt_SAMPLE.eps");
0720 c1->SaveAs("Pullpt_SAMPLE.gif");
0721 c1->WaitPrimitive();
0722
0723 delete c1;
0724 delete leg1;
0725 delete samplelist;
0726 gROOT->Reset();
0727 gROOT->Clear();
0728
0729 }
0730