Back to home page

EIC code displayed by LXR
 
 

    

File indexing completed on 2026-10-08 08:36:29 

0001 #include "TFile.h"
0002 #include "TH2F.h"
0003 #include "TH1F.h"
0004 #include "TCanvas.h"
0005 #include "TLegend.h"
0006 #include "PlottingHeader.h"
0007 
0008 //********************************************************************************************************
0009 //********************************************************************************************************
0010 //********************************** MAIN FUNCTION *******************************************************
0011 //********************************************************************************************************
0012 //********************************************************************************************************
0013 void plotRadiationDose2023() {
0014     StyleSettingsBasics();
0015     // 1. Open the ROOT file
0016     //https://wiki.bnl.gov/EPIC/index.php?title=Radiation_Doses 
0017     TFile* fileDIS = TFile::Open("MEQ_ONLY_bryceCanyon_fluences_beamGas_275_central_7_8_2023.root");
0018     if (!fileDIS || fileDIS->IsZombie()) {
0019         printf("Error: Could not open file.\n");
0020         return;
0021     }
0022     TFile* filePBeamGas = TFile::Open("MEQ_ONLY_proton_bryceCanyon_fluences_beamGas_275_central_7_8_2023.root");
0023     if (!filePBeamGas || filePBeamGas->IsZombie()) {
0024         printf("Error: Could not open file.\n");
0025         return;
0026     }
0027 
0028     // Get 2D hists
0029     TH2F* h2DIS = (TH2F*)fileDIS->Get("neutron_MEQ_fluence");
0030     if (!h2DIS) {
0031         printf("Error: Could not find the 2D histogram.\n");
0032         fileDIS->Close();
0033         return;
0034     }
0035     TH2F* h2PBeamGas = (TH2F*)filePBeamGas->Get("neutron_MEQ_fluence");
0036     if (!h2PBeamGas) {
0037         printf("Error: Could not find the 2D histogram.\n");
0038         filePBeamGas->Close();
0039         return;
0040     }
0041 
0042     
0043     Double_t rangeZLFHCal[9] = {368, 378, 398, 418, 438, 458, 478, 488, 500 };
0044     TH1F* hLFHCalDIS[8];
0045     TH1F* hLFHCalPBeamGas[8];
0046     TH1F* hLFHCalAll[8];
0047     TH1F* hLFHCalAllFirst5years[8];
0048     TH1F* hLFHCalAllFirst15years[8];
0049 
0050     Double_t lumi5years   = 40; 
0051     Double_t lumi15years  = 150; 
0052     // Projections in various ranges
0053     Double_t max= 0;
0054     for (Int_t b = 0; b< 8; b++){
0055       hLFHCalDIS[b]        = (TH1F*)h2DIS->ProjectionY(Form("hLFHCalDIS_%i", b), h2DIS->GetXaxis()->FindBin(rangeZLFHCal[b]), h2DIS->GetXaxis()->FindBin(rangeZLFHCal[b+1]), "E");
0056       hLFHCalDIS[b]->Scale(1./(rangeZLFHCal[b+1]-rangeZLFHCal[b]));
0057       hLFHCalPBeamGas[b]   = (TH1F*)h2PBeamGas->ProjectionY(Form("hLFHCalPBeamGas%i", b), h2PBeamGas->GetXaxis()->FindBin(rangeZLFHCal[b]), h2PBeamGas->GetXaxis()->FindBin(rangeZLFHCal[b+1]), "E");
0058       hLFHCalPBeamGas[b]->Scale(1./(rangeZLFHCal[b+1]-rangeZLFHCal[b]));
0059       hLFHCalAll[b]        = (TH1F*)hLFHCalDIS[b]->Clone(Form("hLFHCalAll_%i", b));
0060       hLFHCalAll[b]->Sumw2();
0061       hLFHCalAll[b]->Add(hLFHCalPBeamGas[b]);
0062       if (max < hLFHCalAll[b]->GetMaximum()) max= hLFHCalAll[b]->GetMaximum();
0063       
0064       hLFHCalAllFirst5years[b]        = (TH1F*)hLFHCalAll[b]->Clone(Form("hLFHCalAllFirst5_%i", b));
0065       hLFHCalAllFirst5years[b]->Scale(lumi5years);
0066       hLFHCalAllFirst15years[b]       = (TH1F*)hLFHCalAll[b]->Clone(Form("hLFHCalAllFirst15_%i", b));
0067       hLFHCalAllFirst15years[b]->Scale(lumi15years);
0068     }
0069     
0070     Double_t textSizeRel = 0.035;
0071     Double_t min = hLFHCalPBeamGas[0]->GetMinimum()/3;
0072     
0073     hLFHCalDIS[0]->GetYaxis()->SetRangeUser(min, max*2);
0074     hLFHCalPBeamGas[0]->GetYaxis()->SetRangeUser(min, max*2);
0075     hLFHCalAll[0]->GetYaxis()->SetRangeUser(min, max*2);
0076     hLFHCalAllFirst5years[0]->GetYaxis()->SetRangeUser(min*lumi5years, max*2*lumi5years);
0077     hLFHCalAllFirst15years[0]->GetYaxis()->SetRangeUser(min*lumi15years, max*2*lumi15years);
0078     
0079     // 4. Draw the projected 1D histogram
0080     TCanvas* c1 = new TCanvas("c1", "Projection Canvas", 800, 600);
0081     DrawCanvasSettings(c1, 0.07, 0.02, 0.01, 0.07);
0082     
0083     TLegend* legend = GetAndSetLegend2( 0.72, 0.8, 0.95, 0.96,textSizeRel*0.85, 2, "", 42,0.2);
0084     
0085     Style_t markerStyle[10] = { 20, 21, 24, 25, 27, 28, 29, 33, 34, 43};
0086     Color_t color[10]       = { kBlack, kViolet+4, kBlue-3, kCyan+1, kGreen+1, 
0087                                 kYellow-4, kOrange, kRed-4, kPink-5, kMagenta+2 };
0088     TString labelLayer[8]   = {"seg. 1", "seg. 2", "seg. 3", "seg. 4", "seg. 5", "seg. 6", "seg. 7", "electronics" };
0089                                 
0090     for (Int_t b = 0; b< 8; b++){
0091       SetMarkerDefaults(hLFHCalDIS[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0092       hLFHCalDIS[b]->GetXaxis()->SetRangeUser(0,280);
0093       SetMarkerDefaults(hLFHCalPBeamGas[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0094       hLFHCalPBeamGas[b]->GetXaxis()->SetRangeUser(0,280);
0095       SetMarkerDefaults(hLFHCalAll[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0096       hLFHCalAll[b]->GetXaxis()->SetRangeUser(0,280);
0097 
0098       SetMarkerDefaults(hLFHCalAllFirst5years[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0099       hLFHCalAllFirst5years[b]->GetXaxis()->SetRangeUser(0,280);
0100 
0101       SetMarkerDefaults(hLFHCalAllFirst15years[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0102       hLFHCalAllFirst15years[b]->GetXaxis()->SetRangeUser(0,280);
0103       
0104       legend->AddEntry(hLFHCalDIS[b],labelLayer[b].Data(), "p");
0105     }
0106 
0107     c1->SetLogy();
0108     c1->cd();
0109     for (Int_t b = 0; b< 8; b++){
0110       if (b == 0) hLFHCalDIS[b]->Draw("pe");
0111       else hLFHCalDIS[b]->Draw("same,pe");
0112       
0113     }
0114     legend->SetHeader("DIS e-p");
0115     legend->Draw();
0116     
0117     DrawLines(20, 20, min,  min*4, 2, kGray+1, 7);
0118     DrawLines(40, 40, min,  min*4, 2, kGray+1, 7);
0119     DrawLines(100, 100, min,  min*4, 2, kGray+1, 7);
0120     DrawLines(250, 250, min,  min*4, 2, kGray+1, 7);
0121     
0122     DrawLatex(0.11, 0.10, "min R right", false,textSizeRel*0.85, 42,0.05, 1);
0123     DrawLatex(0.16, 0.14, "min R left", false,textSizeRel*0.85, 42,0.05, 1);
0124     DrawLatex(0.85, 0.10, "max R", false,textSizeRel*0.85, 42,0.05, 1);
0125     DrawLatex(0.13, 0.18, "insert", false,textSizeRel*0.85, 62,0.05, 1);
0126     DrawLatex(0.22, 0.18, "3x3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0127     DrawLatex(0.56, 0.18, "1.3x1.3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0128     
0129     DrawLatex(0.4, 0.92, "#it{LFHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0130     c1->SaveAs("LFHCal_DIS.pdf");
0131     
0132     c1->cd();
0133     for (Int_t b = 0; b< 8; b++){
0134       if (b == 0) hLFHCalPBeamGas[b]->Draw("pe");
0135       else hLFHCalPBeamGas[b]->Draw("same,pe");
0136     }
0137     legend->SetHeader("p beam gas");
0138     legend->Draw();
0139     DrawLines(20, 20, min,  min*4, 2, kGray+1, 7);
0140     DrawLines(40, 40, min,  min*4, 2, kGray+1, 7);
0141     DrawLines(100, 100, min,  min*4, 2, kGray+1, 7);
0142     DrawLines(250, 250, min,  min*4, 2, kGray+1, 7);
0143     
0144     DrawLatex(0.11, 0.10, "min R right", false,textSizeRel*0.85, 42,0.05, 1);
0145     DrawLatex(0.16, 0.14, "min R left", false,textSizeRel*0.85, 42,0.05, 1);
0146     DrawLatex(0.85, 0.10, "max R", false,textSizeRel*0.85, 42,0.05, 1);
0147     DrawLatex(0.13, 0.18, "insert", false,textSizeRel*0.85, 62,0.05, 1);
0148     DrawLatex(0.22, 0.18, "3x3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0149     DrawLatex(0.56, 0.18, "1.3x1.3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0150     
0151     DrawLatex(0.4, 0.92, "#it{LFHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0152     
0153     c1->SaveAs("LFHCal_PBeamGas.pdf");
0154     
0155     
0156     c1->cd();
0157     for (Int_t b = 0; b< 8; b++){
0158       if (b == 0) hLFHCalAll[b]->Draw("pe");
0159       else hLFHCalAll[b]->Draw("same,pe");
0160     }
0161     legend->SetHeader("DIS e-p + p beam gas");
0162     legend->Draw();
0163     DrawLines(20, 20, min,  min*4, 2, kGray+1, 7);
0164     DrawLines(40, 40, min,  min*4, 2, kGray+1, 7);
0165     DrawLines(100, 100, min,  min*4, 2, kGray+1, 7);
0166     DrawLines(250, 250, min,  min*4, 2, kGray+1, 7);
0167     
0168     DrawLatex(0.11, 0.10, "min R right", false,textSizeRel*0.85, 42,0.05, 1);
0169     DrawLatex(0.16, 0.14, "min R left", false,textSizeRel*0.85, 42,0.05, 1);
0170     DrawLatex(0.85, 0.10, "max R", false,textSizeRel*0.85, 42,0.05, 1);
0171     DrawLatex(0.13, 0.18, "insert", false,textSizeRel*0.85, 62,0.05, 1);
0172     DrawLatex(0.22, 0.18, "3x3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0173     DrawLatex(0.56, 0.18, "1.3x1.3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0174     
0175     DrawLatex(0.4, 0.92, "#it{LFHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0176     
0177     c1->SaveAs("LFHCal_All.pdf");
0178     
0179     
0180     c1->cd();
0181     for (Int_t b = 0; b< 8; b++){
0182       if (b == 0) hLFHCalAllFirst5years[b]->Draw("pe");
0183       else hLFHCalAllFirst5years[b]->Draw("same,pe");
0184     }
0185     legend->SetHeader("DIS e-p + p beam gas");
0186     legend->Draw();
0187     DrawLines(20, 20, min*lumi5years,  min*4*lumi5years, 2, kGray+1, 7);
0188     DrawLines(40, 40, min*lumi5years,  min*4*lumi5years, 2, kGray+1, 7);
0189     DrawLines(100, 100, min*lumi5years,  min*4*lumi5years, 2, kGray+1, 7);
0190     DrawLines(250, 250, min*lumi5years,  min*4*lumi5years, 2, kGray+1, 7);
0191     
0192     DrawLatex(0.11, 0.10, "min R right", false,textSizeRel*0.85, 42,0.05, 1);
0193     DrawLatex(0.16, 0.14, "min R left", false,textSizeRel*0.85, 42,0.05, 1);
0194     DrawLatex(0.85, 0.10, "max R", false,textSizeRel*0.85, 42,0.05, 1);
0195     DrawLatex(0.13, 0.18, "insert", false,textSizeRel*0.85, 62,0.05, 1);
0196     DrawLatex(0.22, 0.18, "3x3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0197     DrawLatex(0.56, 0.18, "1.3x1.3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0198     
0199     DrawLatex(0.4, 0.92, Form("#it{LFHCal} region, %.0f fb^{-1}", lumi5years), false,textSizeRel, 62,0.05, 1);
0200     
0201     c1->SaveAs("LFHCal_All_5Years.pdf");
0202     
0203     c1->cd();
0204     for (Int_t b = 0; b< 8; b++){
0205       if (b == 0) hLFHCalAllFirst15years[b]->Draw("pe");
0206       else hLFHCalAllFirst15years[b]->Draw("same,pe");
0207     }
0208     legend->SetHeader("DIS e-p + p beam gas");
0209     legend->Draw();
0210     DrawLines(20, 20, min*lumi15years,  min*4*lumi15years, 2, kGray+1, 7);
0211     DrawLines(40, 40, min*lumi15years,  min*4*lumi15years, 2, kGray+1, 7);
0212     DrawLines(100, 100, min*lumi15years,  min*4*lumi15years, 2, kGray+1, 7);
0213     DrawLines(250, 250, min*lumi15years,  min*4*lumi15years, 2, kGray+1, 7);
0214     
0215     DrawLatex(0.11, 0.10, "min R right", false,textSizeRel*0.85, 42,0.05, 1);
0216     DrawLatex(0.16, 0.14, "min R left", false,textSizeRel*0.85, 42,0.05, 1);
0217     DrawLatex(0.85, 0.10, "max R", false,textSizeRel*0.85, 42,0.05, 1);
0218     DrawLatex(0.13, 0.18, "insert", false,textSizeRel*0.85, 62,0.05, 1);
0219     DrawLatex(0.22, 0.18, "3x3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0220     DrawLatex(0.56, 0.18, "1.3x1.3 mm SiPMs", false,textSizeRel*0.85, 62,0.05, 1);
0221     
0222     DrawLatex(0.4, 0.92, Form("#it{LFHCal} region, %.0f fb^{-1}", lumi15years), false,textSizeRel, 62,0.05, 1);
0223     
0224     c1->SaveAs("LFHCal_All_15Years.pdf");
0225     
0226     // BHCAL
0227     
0228     Double_t rangeRBHCal[9] = {180, 220, 240, 245, 250  };
0229     TH1F* hBHCalDIS[8];
0230     TH1F* hBHCalPBeamGas[8];
0231     TH1F* hBHCalAll[8];
0232     TH1F* hBHCalAllFirst5years[8];
0233     TH1F* hBHCalAllFirst15years[8];
0234 
0235     // Projections in various ranges
0236     Double_t maxBHcal= 0;
0237     for (Int_t b = 0; b< 4; b++){
0238       hBHCalDIS[b]        = (TH1F*)h2DIS->ProjectionX(Form("hBHCalDIS_%i", b), h2DIS->GetYaxis()->FindBin(rangeRBHCal[b]), h2DIS->GetYaxis()->FindBin(rangeRBHCal[b+1]), "E");
0239       hBHCalDIS[b]->Scale(1./(rangeRBHCal[b+1]-rangeRBHCal[b]));
0240       hBHCalPBeamGas[b]   = (TH1F*)h2PBeamGas->ProjectionX(Form("hBHCalPBeamGas%i", b), h2PBeamGas->GetYaxis()->FindBin(rangeRBHCal[b]), h2PBeamGas->GetYaxis()->FindBin(rangeRBHCal[b+1]), "E");
0241       hBHCalPBeamGas[b]->Scale(1./(rangeRBHCal[b+1]-rangeRBHCal[b]));
0242       hBHCalAll[b]        = (TH1F*)hBHCalDIS[b]->Clone(Form("hBHCalAll_%i", b));
0243       hBHCalAll[b]->Sumw2();
0244       hBHCalAll[b]->Add(hBHCalPBeamGas[b]);
0245       if (maxBHcal < hBHCalAll[b]->GetMaximum()) maxBHcal= hBHCalAll[b]->GetMaximum();
0246       
0247       hBHCalAllFirst5years[b]        = (TH1F*)hBHCalAll[b]->Clone(Form("hBHCalAllFirst5_%i", b));
0248       hBHCalAllFirst5years[b]->Scale(lumi5years);
0249       hBHCalAllFirst15years[b]       = (TH1F*)hBHCalAll[b]->Clone(Form("hBHCalAllFirst15_%i", b));
0250       hBHCalAllFirst15years[b]->Scale(lumi15years);
0251     }
0252     
0253     Double_t minBHCal = hLFHCalPBeamGas[0]->GetMinimum();
0254     maxBHcal = maxBHcal*10;
0255     hBHCalDIS[0]->GetYaxis()->SetRangeUser(minBHCal, maxBHcal*2);
0256     hBHCalPBeamGas[0]->GetYaxis()->SetRangeUser(minBHCal, maxBHcal*2);
0257     hBHCalAll[0]->GetYaxis()->SetRangeUser(minBHCal, maxBHcal*2);
0258     hBHCalAllFirst5years[0]->GetYaxis()->SetRangeUser(minBHCal*lumi5years, maxBHcal*2*lumi5years);
0259     hBHCalAllFirst15years[0]->GetYaxis()->SetRangeUser(minBHCal*lumi15years, maxBHcal*2*lumi15years);
0260     
0261     TLegend* legendBHCal = GetAndSetLegend2( 0.72, 0.80, 0.95, 0.96,textSizeRel*0.85, 1, "", 42,0.2);
0262     
0263     for (Int_t b = 0; b< 4; b++){
0264       SetMarkerDefaults(hBHCalDIS[b], "Z (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0265       hBHCalDIS[b]->GetXaxis()->SetRangeUser(-350,350);
0266       SetMarkerDefaults(hBHCalPBeamGas[b], "Z (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0267       hBHCalPBeamGas[b]->GetXaxis()->SetRangeUser(-350,350);
0268       SetMarkerDefaults(hBHCalAll[b], "Z (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0269       hBHCalAll[b]->GetXaxis()->SetRangeUser(-350,350);
0270 
0271       SetMarkerDefaults(hBHCalAllFirst5years[b], "Z (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0272       hBHCalAllFirst5years[b]->GetXaxis()->SetRangeUser(-350,350);
0273 
0274       SetMarkerDefaults(hBHCalAllFirst15years[b], "Z (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0275       hBHCalAllFirst15years[b]->GetXaxis()->SetRangeUser(-350,350);
0276       
0277       legendBHCal->AddEntry(hBHCalDIS[b],Form("%.0f< R (cm) < %.0f", rangeRBHCal[b], rangeRBHCal[b+1]), "p");
0278     }
0279 
0280     c1->SetLogy();
0281     c1->cd();
0282     for (Int_t b = 0; b< 4; b++){
0283       if (b == 0) hBHCalDIS[b]->Draw("pe");
0284       else hBHCalDIS[b]->Draw("same,pe");
0285       
0286     }
0287     legendBHCal->SetHeader("DIS e-p");
0288     legendBHCal->Draw();
0289     DrawLines(-328, -328, minBHCal,  minBHCal*4, 2, kGray+1, 7);
0290     DrawLines(328, 328, minBHCal,  minBHCal*4, 2, kGray+1, 7);
0291     DrawLatex(0.4, 0.92, "#it{BHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0292     c1->SaveAs("BHCal_MEQ_DIS.pdf");
0293     
0294     c1->cd();
0295     for (Int_t b = 0; b< 4; b++){
0296       if (b == 0) hBHCalPBeamGas[b]->Draw("pe");
0297       else hBHCalPBeamGas[b]->Draw("same,pe");
0298     }
0299     legendBHCal->SetHeader("p beam gas");
0300     legendBHCal->Draw();
0301     DrawLines(-328, -328, minBHCal,  minBHCal*4, 2, kGray+1, 7);
0302     DrawLines(328, 328, minBHCal,  minBHCal*4, 2, kGray+1, 7);
0303     DrawLatex(0.4, 0.92, "#it{BHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0304     
0305     c1->SaveAs("BHCal_MEQ_PBeamGas.pdf");
0306     
0307     
0308     c1->cd();
0309     for (Int_t b = 0; b< 4; b++){
0310       if (b == 0) hBHCalAll[b]->Draw("pe");
0311       else hBHCalAll[b]->Draw("same,pe");
0312     }
0313     legendBHCal->SetHeader("DIS e-p + p beam gas");
0314     legendBHCal->Draw();
0315     DrawLatex(0.4, 0.92, "#it{BHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0316     DrawLines(-328, -328, minBHCal,  minBHCal*4, 2, kGray+1, 7);
0317     DrawLines(328, 328, minBHCal,  minBHCal*4, 2, kGray+1, 7);
0318     
0319     c1->SaveAs("BHCal_MEQ_All.pdf");
0320     
0321     
0322     c1->cd();
0323     for (Int_t b = 0; b< 4; b++){
0324       if (b == 0) hBHCalAllFirst5years[b]->Draw("pe");
0325       else hBHCalAllFirst5years[b]->Draw("same,pe");
0326     }
0327     legendBHCal->SetHeader("DIS e-p + p beam gas");
0328     legendBHCal->Draw();
0329     DrawLines(-328, -328, minBHCal*lumi5years,  minBHCal*4*lumi5years, 2, kGray+1, 7);
0330     DrawLines(328, 328, minBHCal*lumi5years,  minBHCal*4*lumi5years, 2, kGray+1, 7);
0331     DrawLatex(0.4, 0.92, Form("#it{BHCal} region, %.0f fb^{-1}", lumi5years), false,textSizeRel, 62,0.05, 1);
0332     
0333     c1->SaveAs("BHCal_MEQ_All_5Years.pdf");
0334     
0335     c1->cd();
0336     for (Int_t b = 0; b< 4; b++){
0337       if (b == 0) hBHCalAllFirst15years[b]->Draw("pe");
0338       else hBHCalAllFirst15years[b]->Draw("same,pe");
0339     }
0340     legendBHCal->SetHeader("DIS e-p + p beam gas");
0341     legendBHCal->Draw();    
0342     DrawLines(-328, -328, minBHCal*lumi15years,  minBHCal*4*lumi15years, 2, kGray+1, 7);
0343     DrawLines(328, 328, minBHCal*lumi15years,  minBHCal*4*lumi15years, 2, kGray+1, 7);
0344     DrawLatex(0.4, 0.92, Form("#it{BHCal} region, %.0f fb^{-1}", lumi15years), false,textSizeRel, 62,0.05, 1);
0345     
0346     c1->SaveAs("BHCal_MEQ_All_15Years.pdf");
0347     
0348     // NHCal 
0349     
0350     Double_t rangeZNHCal[7] = { 385, 400, 410, 420, 430, 440, 450 };
0351     TH1F* hNHCalDIS[8];
0352     TH1F* hNHCalPBeamGas[8];
0353     TH1F* hNHCalAll[8];
0354     TH1F* hNHCalAllFirst5years[8];
0355     TH1F* hNHCalAllFirst15years[8];
0356 
0357     // Projections in various ranges
0358     Double_t maxNHCal= 0;
0359     for (Int_t b = 0; b< 6; b++){
0360       hNHCalDIS[b]        = (TH1F*)h2DIS->ProjectionY(Form("hNHCalDIS_%i", b), h2DIS->GetXaxis()->FindBin(rangeZNHCal[b]), h2DIS->GetXaxis()->FindBin(rangeZNHCal[b+1]), "E");
0361       hNHCalDIS[b]->Scale(1./(TMath::Abs(rangeZNHCal[b+1]-rangeZNHCal[b])));
0362       hNHCalPBeamGas[b]   = (TH1F*)h2PBeamGas->ProjectionY(Form("hNHCalPBeamGas%i", b), h2PBeamGas->GetXaxis()->FindBin(rangeZNHCal[b]), h2PBeamGas->GetXaxis()->FindBin(rangeZNHCal[b+1]), "E");
0363       hNHCalPBeamGas[b]->Scale(1./(TMath::Abs(rangeZNHCal[b+1]-rangeZNHCal[b])));
0364       hNHCalAll[b]        = (TH1F*)hNHCalDIS[b]->Clone(Form("hNHCalAll_%i", b));
0365       hNHCalAll[b]->Sumw2();
0366       hNHCalAll[b]->Add(hNHCalPBeamGas[b]);
0367       if (maxNHCal < hNHCalAll[b]->GetMaximum()) maxNHCal= hNHCalAll[b]->GetMaximum();
0368       
0369       hNHCalAllFirst5years[b]        = (TH1F*)hNHCalAll[b]->Clone(Form("hNHCalAllFirst5_%i", b));
0370       hNHCalAllFirst5years[b]->Scale(lumi5years);
0371       hNHCalAllFirst15years[b]       = (TH1F*)hNHCalAll[b]->Clone(Form("hNHCalAllFirst15_%i", b));
0372       hNHCalAllFirst15years[b]->Scale(lumi15years);
0373     }
0374     
0375     Double_t minNHCal = hNHCalPBeamGas[0]->GetMinimum()/3;
0376     
0377     hNHCalDIS[0]->GetYaxis()->SetRangeUser(minNHCal, maxNHCal*2);
0378     hNHCalPBeamGas[0]->GetYaxis()->SetRangeUser(minNHCal, maxNHCal*2);
0379     hNHCalAll[0]->GetYaxis()->SetRangeUser(minNHCal, maxNHCal*2);
0380     hNHCalAllFirst5years[0]->GetYaxis()->SetRangeUser(minNHCal*lumi5years, maxNHCal*2*lumi5years);
0381     hNHCalAllFirst15years[0]->GetYaxis()->SetRangeUser(minNHCal*lumi15years, maxNHCal*2*lumi15years);
0382     
0383     TString labelLayerNHCal[8]   = {"electronics", "seg. 1", "seg. 2", "seg. 3", "seg. 4", "seg. 5", "seg. 6", "seg. 7"};
0384     TLegend* legendNHCal = GetAndSetLegend2( 0.72, 0.84, 0.95, 0.96,textSizeRel*0.85, 2, "", 42,0.2);
0385                                     
0386     for (Int_t b = 0; b< 6; b++){
0387       SetMarkerDefaults(hNHCalDIS[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0388       hNHCalDIS[b]->GetXaxis()->SetRangeUser(0,280);
0389       SetMarkerDefaults(hNHCalPBeamGas[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0390       hNHCalPBeamGas[b]->GetXaxis()->SetRangeUser(0,280);
0391       SetMarkerDefaults(hNHCalAll[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0392       hNHCalAll[b]->GetXaxis()->SetRangeUser(0,280);
0393 
0394       SetMarkerDefaults(hNHCalAllFirst5years[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0395       hNHCalAllFirst5years[b]->GetXaxis()->SetRangeUser(0,280);
0396 
0397       SetMarkerDefaults(hNHCalAllFirst15years[b], "R (cm)", "1 MEQ Neutron Fluence (cm^{-2})", markerStyle[b], 1, color[b], color[b], textSizeRel, 0.85*textSizeRel, 1., 1.);
0398       hNHCalAllFirst15years[b]->GetXaxis()->SetRangeUser(0,280);
0399       
0400       legendNHCal->AddEntry(hNHCalDIS[b],labelLayerNHCal[b].Data(), "p");
0401     }
0402 
0403     c1->SetLogy();
0404     c1->cd();
0405     for (Int_t b = 0; b< 6; b++){
0406       if (b == 0) hNHCalDIS[b]->Draw("pe");
0407       else hNHCalDIS[b]->Draw("same,pe");
0408       
0409     }
0410     legendNHCal->SetHeader("DIS e-p");
0411     legendNHCal->Draw();
0412     
0413     DrawLines(10, 10, minNHCal,  minNHCal*4, 2, kGray+1, 7);
0414     DrawLines(253, 253, minNHCal,  minNHCal*4, 2, kGray+1, 7);
0415         
0416     DrawLatex(0.4, 0.92, "#it{NHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0417     c1->SaveAs("NHCal_DIS.pdf");
0418     
0419     c1->cd();
0420     for (Int_t b = 0; b< 6; b++){
0421       if (b == 0) hNHCalPBeamGas[b]->Draw("pe");
0422       else hNHCalPBeamGas[b]->Draw("same,pe");
0423     }
0424     legendNHCal->SetHeader("p beam gas");
0425     legendNHCal->Draw();
0426     DrawLines(10, 10, minNHCal,  minNHCal*4, 2, kGray+1, 7);
0427     DrawLines(253, 253, minNHCal,  minNHCal*4, 2, kGray+1, 7);
0428     
0429     DrawLatex(0.4, 0.92, "#it{NHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0430     
0431     c1->SaveAs("NHCal_PBeamGas.pdf");
0432     
0433     
0434     c1->cd();
0435     for (Int_t b = 0; b< 6; b++){
0436       if (b == 0) hNHCalAll[b]->Draw("pe");
0437       else hNHCalAll[b]->Draw("same,pe");
0438     }
0439     legendNHCal->SetHeader("DIS e-p + p beam gas");
0440     legendNHCal->Draw();
0441     DrawLines(10, 10, minNHCal,  minNHCal*4, 2, kGray+1, 7);
0442     DrawLines(253, 253, minNHCal,  minNHCal*4, 2, kGray+1, 7);
0443 
0444     DrawLatex(0.4, 0.92, "#it{NHCal} region, 1 fb^{-1}", false,textSizeRel, 62,0.05, 1);
0445     
0446     c1->SaveAs("NHCal_All.pdf");
0447     
0448     
0449     c1->cd();
0450     for (Int_t b = 0; b< 6; b++){
0451       if (b == 0) hNHCalAllFirst5years[b]->Draw("pe");
0452       else hNHCalAllFirst5years[b]->Draw("same,pe");
0453     }
0454     legendNHCal->SetHeader("DIS e-p + p beam gas");
0455     legendNHCal->Draw();
0456     DrawLines(10, 10, minNHCal*lumi5years,  minNHCal*4*lumi5years, 2, kGray+1, 7);
0457     DrawLines(253, 253, minNHCal*lumi5years,  minNHCal*4*lumi5years, 2, kGray+1, 7);
0458     DrawLatex(0.4, 0.92, Form("#it{NHCal} region, %.0f fb^{-1}", lumi5years), false,textSizeRel, 62,0.05, 1);
0459     
0460     c1->SaveAs("NHCal_All_5Years.pdf");
0461     
0462     c1->cd();
0463     for (Int_t b = 0; b< 6; b++){
0464       if (b == 0) hNHCalAllFirst15years[b]->Draw("pe");
0465       else hNHCalAllFirst15years[b]->Draw("same,pe");
0466     }
0467     legendNHCal->SetHeader("DIS e-p + p beam gas");
0468     legendNHCal->Draw();
0469     DrawLines(10, 10, minNHCal*lumi15years,  minNHCal*4*lumi15years, 2, kGray+1, 7);
0470     DrawLines(253, 253, minNHCal*lumi15years,  minNHCal*4*lumi15years, 2, kGray+1, 7);
0471     
0472     DrawLatex(0.4, 0.92, Form("#it{NHCal} region, %.0f fb^{-1}", lumi15years), false,textSizeRel, 62,0.05, 1);
0473     
0474     c1->SaveAs("NHCal_All_15Years.pdf");
0475 
0476     delete c1;
0477 }
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!