17 |
|
#include "Math/GSLIntegrator.h" |
18 |
|
#include "Math/WrappedTF1.h" |
19 |
|
|
20 |
< |
// TString fileName = "skimmed/QCDbin1_OldRelIso_rebin.root"; |
21 |
< |
TString fileName = "skimmed/QCDbin4_OldRelIso_FastSim.root"; |
20 |
> |
// TString fileName = "skimmed/QCDbin4_OldRelIso_rebin.root"; |
21 |
> |
TString fileName = "skimmed/QCDbin4_OldRelIso_FastSim.root"; |
22 |
|
|
23 |
|
// fit region weight |
24 |
|
double wxmin = 0.07; // 0.13 for FullSim 0.07 for FsatSim |
81 |
|
} |
82 |
|
h_all->Draw(); |
83 |
|
|
84 |
+ |
Int_t nbxMax = h_all->GetXaxis()->FindBin(1.0)+1; |
85 |
+ |
h_all->GetXaxis()->SetRange(1,h_all->GetXaxis()->FindBin(bxf)); |
86 |
|
Int_t nbMax=h_all->GetMaximumBin(); |
87 |
|
cout << "Bin number at peak= " << h_all->GetMaximumBin() << endl; |
88 |
+ |
h_all->GetXaxis()->SetRange(1,nbxMax); |
89 |
|
TAxis *axis0 = h_all->GetXaxis(); |
90 |
|
double bwidth = axis0->GetBinWidth(nbMax); |
91 |
|
double nbin = axis0->GetLast(); |