| 15 |
|
#include <TPaveStats.h> |
| 16 |
|
#include <TTree.h> |
| 17 |
|
|
| 18 |
< |
TString fileName = "skimmed/QCDbin4_OldRelIso_rebin.root"; |
| 19 |
< |
// TString fileName = "skimmed/QCDbin4_OldRelIso_FastSim.root"; |
| 18 |
> |
// Output file |
| 19 |
> |
TString fileName = "skimmed/QCDbinall_NewRelIso.root"; |
| 20 |
|
bool FastSim = false; |
| 21 |
< |
TString isoname = "Old"; |
| 21 |
> |
TString isoname = "New"; |
| 22 |
|
// Jet bin: |
| 23 |
< |
Int_t jbin = 4; // 1-3, 4 for >=4. Put in number less than 5 |
| 24 |
< |
double binwidth = 0.02; // 0.02,0.02,0.05,0.1 for New RelIso |
| 23 |
> |
Int_t jbin = 0; // 1-3, 4 for >=4. Put in number less than 5; 0 for all jets |
| 24 |
> |
double binwidth = 0.01; // 0.02,0.02,0.05,0.1 for New RelIso |
| 25 |
|
// Histo range |
| 26 |
|
double xMin = 0.0; |
| 27 |
|
double xMax = 0.0; |
| 147 |
|
bool goodD0sig = false; // d0sig condition |
| 148 |
|
|
| 149 |
|
if ( jbin > 4 ) {cout<<"Don't mess around!!! jbin should be less than 5!"<<endl;return;} |
| 150 |
< |
if (top_njets == jbin && jbin != 4) correctBin = true; |
| 151 |
< |
else if (top_njets >= jbin && jbin == 4) correctBin = true; |
| 150 |
> |
if ( jbin == 0 ){ |
| 151 |
> |
if ( top_njets > jbin ) { correctBin = true; } |
| 152 |
> |
} |
| 153 |
> |
else if ( jbin == 4 ){ |
| 154 |
> |
if ( top_njets >= jbin ) { correctBin = true; } |
| 155 |
> |
} |
| 156 |
> |
else { |
| 157 |
> |
if ( top_njets == jbin ) { correctBin = true; } |
| 158 |
> |
} |
| 159 |
> |
|
| 160 |
|
if ( d0sig < d0sigCut ) goodD0sig = true; |
| 161 |
|
|
| 162 |
|
// Counting |
