52 |
|
TH1D *hx0_Fill = new TH1D("hx0_Fill","x_{0} vs. LHC Fill",nentries,0.5,nentries+0.5); |
53 |
|
TH1D *hy0_Fill = new TH1D("hy0_Fill","y_{0} vs. LHC Fill",nentries,0.5,nentries+0.5); |
54 |
|
TH1D *hz0_Fill = new TH1D("hz0_Fill","z_{0} vs. LHC Fill",nentries,0.5,nentries+0.5); |
55 |
+ |
TH1D *hsigmaZ0_Fill = new TH1D("hsigmaZ0_Fill","#sigma_{z_{0}} vs. LHC Fill",nentries,0.5,nentries+0.5); |
56 |
|
|
57 |
|
for (Int_t i=0;i<nentries;i++) { |
58 |
|
chain.GetEntry(i); |
71 |
|
hz0_Fill->SetBinContent(i+1,z); |
72 |
|
hz0_Fill->SetBinError(i+1,zErr); |
73 |
|
hz0_Fill->GetXaxis()->SetBinLabel(i+1,xlabel); |
74 |
+ |
hsigmaZ0_Fill->SetBinContent(i+1,sigmaZ); |
75 |
+ |
hsigmaZ0_Fill->SetBinError(i+1,sigmaZErr); |
76 |
+ |
hsigmaZ0_Fill->GetXaxis()->SetBinLabel(i+1,xlabel); |
77 |
|
} |
78 |
|
|
79 |
|
gStyle->SetErrorX(0); |
129 |
|
|
130 |
|
c3->Print("z0_Fill.png"); |
131 |
|
|
132 |
+ |
gPad->Update(); |
133 |
+ |
|
134 |
+ |
TCanvas *c4 = new TCanvas("c4","c4",756,600); |
135 |
+ |
hsigmaZ0_Fill->GetYaxis()->SetRangeUser(hsigmaZ0_Fill->GetMinimum()-0.5,hsigmaZ0_Fill->GetMaximum()+0.9); |
136 |
+ |
hsigmaZ0_Fill->GetXaxis()->SetTitle("LHC Fill"); |
137 |
+ |
hsigmaZ0_Fill->GetYaxis()->SetTitle("#sigma_{z} of beam spots (cm)"); |
138 |
+ |
hsigmaZ0_Fill->Draw("E1"); |
139 |
+ |
ppt->Draw(); |
140 |
+ |
// l0->Draw(); |
141 |
+ |
|
142 |
+ |
c4->Print("sigmaZ0_Fill.png"); |
143 |
+ |
|
144 |
|
} |