| 1 |
#include <iostream>
|
| 2 |
#include <TCut.h>
|
| 3 |
#ifndef SampleClassLoaded
|
| 4 |
#include "SampleClass.C"
|
| 5 |
#endif
|
| 6 |
|
| 7 |
|
| 8 |
using namespace std;
|
| 9 |
|
| 10 |
|
| 11 |
namespace PlottingSetup {
|
| 12 |
float luminosity=1078;//751.0;//486.0;//468.0//336.;//pb^{-1}
|
| 13 |
string jzbvariabledata="jzb[1]";
|
| 14 |
string jzbvariablemc="jzb[1]";
|
| 15 |
|
| 16 |
samplecollection allsamples("completesamplecollection");
|
| 17 |
int data=1;
|
| 18 |
int mc=0;
|
| 19 |
int mcwithsignal=2;
|
| 20 |
TCut passtrig("(passed_triggers||!is_data)");
|
| 21 |
TCut cutmass("abs(mll-91.2)<20");
|
| 22 |
//TCut cutmass("mll>2");
|
| 23 |
TCut basiccut("mll>2");//basically nothing.
|
| 24 |
TCut basicqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//don't use this for the "essential cut", because we want to plot nJets as well as mll in the inclusive case; we thus use it as an addition nJets cut.
|
| 25 |
//TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut
|
| 26 |
//TCut jetqualitycut("mll>0");
|
| 27 |
TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut);
|
| 28 |
TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
|
| 29 |
TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
|
| 30 |
//TCut basiccut("(passed_triggers||!is_data)");
|
| 31 |
|
| 32 |
TCut essential(passtrig);//add here any cuts you ALWAYS want
|
| 33 |
int dogaus=0;
|
| 34 |
int doKM=1;
|
| 35 |
int dogaus2sigma=2;
|
| 36 |
int dogaus3sigma=3;
|
| 37 |
int Kostasmethod=-99;
|
| 38 |
|
| 39 |
float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range!
|
| 40 |
float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range!
|
| 41 |
|
| 42 |
int method=-1;
|
| 43 |
|
| 44 |
}
|
