Back to home page

sPhenix code displayed by LXR
 
 

    

File indexing completed on 2025-08-06 08:11:56 

0001 void GammaNeutronEnergy2DGraphMaker()
0002 {
0003   TFile *fin = TFile::Open("Gamma_Neutron_Energy_Graphs.root");
0004   gROOT->cd();
0005   TCanvas *c1 = new TCanvas();
0006   gStyle->SetOptStat(0);
0007   TH2F *n2h2 = new TH2F("Percentage of Energy Deposited by 2 GeV Particles vs Cone Size","Percentage of Energy Deposited by 2 GeV Particles vs Cone Size",1,0,1,1,0,1.1);
0008   TGraph *n2gr = (TGraph *) fin->Get("neutron2GeV")->Clone();
0009   n2h2->SetXTitle("Cone Size");
0010   n2h2->SetYTitle("Percentage of Energy Deposited");
0011   n2h2->Draw();
0012   n2gr->Draw("*");
0013   TGraph *g2gr = (TGraph *) fin->Get("gamma2GeV")->Clone();
0014   //TH2F *g2h2 = new TH2F("g2h2","g2h2",1,0,1,1,0,1.1);
0015   //  g2h2->Draw("same");
0016   g2gr->Draw("*");
0017   TLine *tl = new TLine();
0018   tl->SetLineStyle(2);
0019   tl->DrawLine(0.1,0,0.1,1.1);
0020  TLine *t2 = new TLine();
0021   t2->SetLineStyle(6);
0022   t2->DrawLine(0.024,0,0.024,1.1);
0023   TLegend *legrda = new TLegend(0.63,0.24,0.84,0.59,NULL,"brNDC");
0024   legrda->SetLineColor(1);
0025   legrda->SetLineStyle(1);
0026   legrda->SetLineWidth(1);
0027   legrda->SetFillColor(10);
0028   legrda->SetFillStyle(1001);
0029   legrda->SetBorderSize(0);
0030   //  legrda->SetTextSize(labelsize);
0031   //legrda->AddEntry(hjbkg,"HIJING bkg"); 
0032   legrda->AddEntry(gamma2GeV,          "2 GeV Gamma","p"); 
0033   legrda->AddEntry(neutron2GeV,          "2 GeV Neutron", "p"); 
0034   legrda->AddEntry(tl,          "HCAL tower size","l"); 
0035   legrda->AddEntry(t2,"EMCAL tower size","l");
0036   legrda->Draw();
0037   c1->Print("Gamma_vs_Neutron_2GeV.png");
0038   
0039   TCanvas *c2 = new TCanvas();
0040   gStyle->SetOptStat(0);
0041   TH2F *n50h2 = new TH2F("Percentage of Energy Deposited by 50 GeV Particles vs Cone Size","Percentage of Energy Deposited by 50 GeV Particles vs Cone Size",1,0,1,1,0,1.1);
0042   TGraph *n50gr = (TGraph *) fin->Get("neutron50GeV")->Clone();
0043   n50h2->SetXTitle("Cone Size");
0044   n50h2->SetYTitle("Percentage of Energy Deposited");
0045   n50h2->Draw();
0046   n50gr->Draw("*");
0047   TGraph *g50gr = (TGraph *) fin->Get("gamma50GeV")->Clone();
0048   //TH2F *g50h2 = new TH2F("g50h2","g50h2",1,0,1,1,0,1.1);
0049   //g50h2->Draw("same");
0050   g50gr->Draw("*");
0051   TLine *tl = new TLine();
0052   tl->SetLineStyle(2);
0053   tl->DrawLine(0.1,0,0.1,1.1);
0054  TLine *t2 = new TLine();
0055   t2->SetLineStyle(6);
0056   t2->DrawLine(0.024,0,0.024,1.1);
0057   TLegend *legrda2 = new TLegend(0.63,0.24,0.84,0.59,NULL,"brNDC");
0058   legrda2->SetLineColor(1);
0059   legrda2->SetLineStyle(1);
0060   legrda2->SetLineWidth(1);
0061   legrda2->SetFillColor(10);
0062   legrda2->SetFillStyle(1001);
0063   legrda2->SetBorderSize(0);
0064   //  legrda2->SetTextSize(labelsize);
0065   //legrda2->AddEntry(hjbkg,"HIJING bkg"); 
0066   legrda2->AddEntry(gamma50GeV,          "50 GeV Gamma","p"); 
0067   legrda2->AddEntry(neutron50GeV,          "50 GeV Neutron", "p"); 
0068   legrda2->AddEntry(tl,          "HCAL tower size","l");
0069   legrda2->AddEntry(t2,"EMCAL tower size","l"); 
0070   legrda2->Draw();
0071   c2->Print("Gamma_vs_Neutron_50GeV.png");
0072   
0073 
0074   TCanvas *c3 = new TCanvas();
0075   gStyle->SetOptStat(0);
0076   TH2F *ga2 = new TH2F("Gamma Energy Deposition vs Cone Size","Gamma Energy Deposition vs Cone Size",1,0,0.5,1,0,1.1);
0077   TGraph *gamma2GeV = (TGraph *) fin->Get("gamma2GeV")->Clone();
0078   ga2->SetXTitle("Cone Size");
0079   ga2->SetYTitle("Percentage of Energy Deposited");
0080   ga2->Draw();
0081   gamma2GeV->SetMarkerSize(1);
0082   gamma2GeV->SetMarkerStyle(4);
0083   gamma2GeV->Draw("P");
0084   int pz[7] = {4,8,12,16,32,40,50};
0085   int col[7] = {2,9,4,5,6,7,8};
0086   int mark[7] = {24,25,26,27,28,35,36};
0087   for (int i=0; i<7;i++)
0088     {
0089       char hname[100];
0090       sprintf(hname,"gamma%dGeV",pz[i]);
0091       //TH2F *gall = new TH2F(hname,hname,1,0,0.5,1,0,1.1);
0092       TGraph *gamall = (TGraph *) fin->Get(hname);
0093       //gall->Draw("same");
0094       gamall->SetMarkerColor(col[i]);
0095       gamall->SetMarkerStyle(4);
0096       gamall->SetMarkerSize(1);
0097       //gamall->SetLineWidth(3);
0098       gamall->Draw("P");
0099     }
0100   TLine *tl = new TLine();
0101   tl->SetLineStyle(2);
0102   tl->DrawLine(0.1,0,0.1,1.1);
0103   TLine *t2 = new TLine();
0104   t2->SetLineStyle(6);
0105   t2->DrawLine(0.024,0,0.024,1.1);
0106   TLegend *legrda3 = new TLegend(0.63,0.24,0.84,0.59,NULL,"brNDC");
0107   legrda3->SetLineColor(1);
0108   legrda3->SetLineStyle(1);
0109   legrda3->SetLineWidth(1);
0110   legrda3->SetFillColor(10);
0111   legrda3->SetFillStyle(1001);
0112   legrda3->SetBorderSize(0);
0113   //  legrda3->SetTextSize(labelsize);
0114   //legrda3->AddEntry(hjbkg,"HIJING bkg");
0115   legrda3->AddEntry(gamma2GeV,          "2 GeV Gamma","p");
0116     legrda3->AddEntry(gamma4GeV,          "4 GeV Gamma","p");
0117     legrda3->AddEntry(gamma8GeV,          "8 GeV Gamma","p");
0118     legrda3->AddEntry(gamma12GeV,          "12 GeV Gamma","p");
0119     legrda3->AddEntry(gamma16GeV,          "16 GeV Gamma","p");
0120     legrda3->AddEntry(gamma32GeV,          "32 GeV Gamma","p");
0121     legrda3->AddEntry(gamma40GeV,          "40 GeV Gamma","p");
0122     legrda3->AddEntry(gamma50GeV,          "50 GeV Gamma","p");
0123    legrda3->AddEntry(tl,          "HCAL tower size","l" );
0124    legrda3->AddEntry(t2,"EMCAL tower size", "l");
0125   legrda3->Draw();
0126   c3->Print("Gamma_All.png");
0127   
0128   TCanvas *c4 = new TCanvas();
0129   gStyle->SetOptStat(0);
0130   TH2F *an2 = new TH2F("Antineutron Energy Deposition vs Cone Size","Antineutron Energy Deposition vs Cone Size",1,0,0.5,1,0,1.1);
0131   TGraph *anti_neutron2GeV = (TGraph *) fin->Get("anti_neutron2GeV")->Clone();
0132   an2->SetXTitle("Cone Size");
0133   an2->SetYTitle("Percentage of Energy Deposited");
0134   an2->Draw();
0135   anti_neutron2GeV->SetMarkerSize(1);
0136   anti_neutron2GeV->SetMarkerStyle(4);
0137   anti_neutron2GeV->Draw("P");
0138   int pz[7] = {4,8,12,16,32,40,50};
0139   int col[7] = {2,9,4,5,6,7,8};
0140   int mark[7] = {24,25,26,27,28,35,36};
0141   for (int i=0; i<7;i++)
0142     {
0143       char hname[100];
0144       sprintf(hname,"anti_neutron%dGeV",pz[i]);
0145       //TH2F *gall = new TH2F(hname,hname,1,0,0.5,1,0,1.1);
0146       TGraph *anall = (TGraph *) fin->Get(hname);
0147       //gall->Draw("same");
0148       anall->SetMarkerColor(col[i]);
0149       anall->SetMarkerStyle(4);
0150       anall->SetMarkerSize(1);
0151       //gamall->SetLineWidth(3);
0152       anall->Draw("P");
0153     }
0154   TLine *tl = new TLine();
0155   tl->SetLineStyle(2);
0156   tl->DrawLine(0.1,0,0.1,1.1);
0157   TLine *t2 = new TLine();
0158   t2->SetLineStyle(6);
0159   t2->DrawLine(0.024,0,0.024,1.1);
0160   TLegend *legrda4 = new TLegend(0.63,0.24,0.84,0.59,NULL,"brNDC");
0161   legrda4->SetLineColor(1);
0162   legrda4->SetLineStyle(1);
0163   legrda4->SetLineWidth(1);
0164   legrda4->SetFillColor(10);
0165   legrda4->SetFillStyle(1001);
0166   legrda4->SetBorderSize(0);
0167   //  legrda4->SetTextSize(labelsize);
0168   //legrda4->AddEntry(hjbkg,"HIJING bkg");
0169 legrda4->AddEntry(anti_neutron2GeV,          "2 GeV Anti Neutron","p");
0170 legrda4->AddEntry(anti_neutron4GeV,          "4 GeV Anti Neutron","p");
0171 legrda4->AddEntry(anti_neutron8GeV,          "8 GeV Anti Neutron","p");
0172 legrda4->AddEntry(anti_neutron12GeV,          "12 GeV Anti Neutron","p");
0173 legrda4->AddEntry(anti_neutron16GeV,          "16 GeV Anti Neutron","p");
0174 legrda4->AddEntry(anti_neutron32GeV,          "32 GeV Anti Neutron","p");
0175 legrda4->AddEntry(anti_neutron40GeV,          "40 GeV Anti Neutron","p");
0176 legrda4->AddEntry(anti_neutron50GeV,          "50 GeV Anti Neutron","p");
0177    legrda4->AddEntry(tl,          "HCAL tower size","l"); 
0178    legrda4->AddEntry(t2,"EMCAL Tower size","l");
0179   legrda4->Draw();
0180   c4->Print("Antineutron_All.png");
0181 
0182   TCanvas *c5 = new TCanvas();
0183   gStyle->SetOptStat(0);
0184   TH2F *n2 = new TH2F("Neutron Energy Deposition vs Cone Size","Neutron Energy Deposition vs Cone Size",1,0,0.5,1,0,1.1);
0185   TGraph *neutron2GeV = (TGraph *) fin->Get("neutron2GeV")->Clone();
0186   n2->SetXTitle("Cone Size");
0187   n2->SetYTitle("Percentage of Energy Deposited");
0188   n2->Draw();
0189   neutron2GeV->SetMarkerSize(1);
0190   neutron2GeV->SetMarkerStyle(4);
0191   neutron2GeV->Draw("P");
0192   int pz[7] = {4,8,12,16,32,40,50};
0193   int col[7] = {2,9,4,5,6,7,8};
0194   int mark[7] = {24,25,26,27,28,35,36};
0195   for (int i=0; i<7;i++)
0196     {
0197       char hname[100];
0198       sprintf(hname,"neutron%dGeV",pz[i]);
0199       //TH2F *gall = new TH2F(hname,hname,1,0,0.5,1,0,1.1);
0200       TGraph *nall = (TGraph *) fin->Get(hname);
0201       //gall->Draw("same");
0202       nall->SetMarkerColor(col[i]);
0203       nall->SetMarkerStyle(4);
0204       nall->SetMarkerSize(1);
0205       //gamall->SetLineWidth(3);
0206       nall->Draw("P");
0207     }
0208  TLine *tl = new TLine();
0209   tl->SetLineStyle(2);
0210   tl->DrawLine(0.1,0,0.1,1.1);
0211  TLine *t2 = new TLine();
0212   t2->SetLineStyle(6);
0213   t2->DrawLine(0.024,0,0.024,1.1);
0214   TLegend *legrda5 = new TLegend(0.63,0.24,0.84,0.59,NULL,"brNDC");
0215   legrda5->SetLineColor(1);
0216   legrda5->SetLineStyle(1);
0217   legrda5->SetLineWidth(1);
0218   legrda5->SetFillColor(10);
0219   legrda5->SetFillStyle(1001);
0220   legrda5->SetBorderSize(0);
0221   //  legrda->SetTextSize(labelsize);
0222   //legrda->AddEntry(hjbkg,"HIJING bkg");
0223 legrda5->AddEntry(neutron2GeV,          "2 GeV Neutron","p");
0224 legrda5->AddEntry(neutron4GeV,          "4 GeV Neutron","p");
0225 legrda5->AddEntry(neutron8GeV,          "8 GeV Neutron","p");
0226 legrda5->AddEntry(neutron12GeV,          "12 GeV Neutron","p");
0227 legrda5->AddEntry(neutron16GeV,          "16 GeV Neutron","p");
0228 legrda5->AddEntry(neutron32GeV,          "32 GeV Neutron","p");
0229 legrda5->AddEntry(neutron40GeV,          "40 GeV Neutron","p");
0230 legrda5->AddEntry(neutron50GeV,          "50 GeV Neutron","p");
0231    legrda5->AddEntry(tl,          "HCAL tower size","l"); 
0232    legrda5->AddEntry(t2,"EMCAL tower size","l");
0233   legrda5->Draw();
0234   c5->Print("Neutron_All.png");
0235 
0236   fin->Close();
0237 }
0238
[ 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!