23 |
|
cout << "makePlots(4) for glsbottom chi150" << endl; |
24 |
|
cout << "makePlots(5) for glsbottom chi300" << endl; |
25 |
|
cout << "makePlots(6) for glsbottom sbpair" << endl; |
26 |
+ |
cout << "The ntuples are to be found in $NTUPLE_DIR" << endl; |
27 |
|
return; |
28 |
|
} |
29 |
|
|
39 |
|
//return; |
40 |
|
|
41 |
|
// Open the ntuple |
42 |
< |
cout << "We are opening " << Form("Jun21/ntuple_%s.root",modelName) << endl; |
43 |
< |
TFile *file = TFile::Open(Form("Jun21/ntuple_%s.root",modelName)); |
42 |
> |
cout << "We are opening " << Form("$NTUPLE_DIR/ntuple_%s.root",modelName) << endl; |
43 |
> |
TFile *file = TFile::Open(Form("$NTUPLE_DIR/ntuple_%s.root",modelName)); |
44 |
|
file->cd(); |
45 |
|
TTree *tree = (TTree*)gDirectory->Get("tree"); |
46 |
|
|
154 |
|
TH2F* empty = new TH2F("empty","empty",nx,xmin,xmin+nx*xbinsize,ny,ymin,ymin+ny*ybinsize); |
155 |
|
|
156 |
|
|
157 |
< |
tree->Draw("lspmass:glmass>>ul","explimsrb/(3950.*effsrb)"); |
157 |
> |
tree->Draw("lspmass:glmass>>ul","obslimsrb/(3950.*effsrb)"); |
158 |
|
tree->Draw("lspmass:glmass>>ulbest","bestsr"); |
159 |
|
tree->Draw("lspmass:glmass>>acc","effsrb"); |
160 |
< |
tree->Draw("lspmass:glmass>>excl","explimsrb/(3950.*effsrb)<xsec"); |
161 |
< |
tree->Draw("lspmass:glmass>>exclup","explimsrb/(3950.*effsrb)<xsecup"); |
162 |
< |
tree->Draw("lspmass:glmass>>excldwn","explimsrb/(3950.*effsrb)<xsecdwn"); |
160 |
> |
tree->Draw("lspmass:glmass>>excl","obslimsrb/(3950.*effsrb)<xsec"); |
161 |
> |
tree->Draw("lspmass:glmass>>exclup","obslimsrb/(3950.*effsrb)<xsecup"); |
162 |
> |
tree->Draw("lspmass:glmass>>excldwn","obslimsrb/(3950.*effsrb)<xsecdwn"); |
163 |
|
tree->Draw("lspmass:glmass>>hxsec", "xsec"); |
164 |
|
tree->Draw("lspmass:glmass>>hxsecup", "xsecup"); |
165 |
|
tree->Draw("lspmass:glmass>>hxsecdwn","xsecdwn"); |