| 1 |
#include <iostream>
|
| 2 |
#include "TCanvas.h"
|
| 3 |
#include "TError.h"
|
| 4 |
#include "TPad.h"
|
| 5 |
#include "TString.h"
|
| 6 |
#include "TRandom.h"
|
| 7 |
#include "TH1F.h"
|
| 8 |
#include "TMath.h"
|
| 9 |
#include "TF1.h"
|
| 10 |
|
| 11 |
#include "TFile.h"
|
| 12 |
#include "TTree.h"
|
| 13 |
#include "TH1D.h"
|
| 14 |
#include "TH2D.h"
|
| 15 |
#include "TCanvas.h"
|
| 16 |
#include "TLegend.h"
|
| 17 |
#include "TLatex.h"
|
| 18 |
#include "TString.h"
|
| 19 |
#include "TChain.h"
|
| 20 |
|
| 21 |
#include "TCut.h"
|
| 22 |
#include "TNtuple.h"
|
| 23 |
|
| 24 |
#include "THStack.h"
|
| 25 |
|
| 26 |
using namespace std;
|
| 27 |
|
| 28 |
static bool subtract = 0;
|
| 29 |
static int QID = 3;
|
| 30 |
|
| 31 |
static TString weightString;
|
| 32 |
|
| 33 |
static bool normLead = 0;
|
| 34 |
|
| 35 |
static int mixColor = 2;
|
| 36 |
static int dataColor = 1;
|
| 37 |
static int ppColor = 4;
|
| 38 |
|
| 39 |
static int centralBin = 8;
|
| 40 |
static int leadCut = 120;
|
| 41 |
static int subleadCut = 30;
|
| 42 |
|
| 43 |
static double sideMin = 0.1;
|
| 44 |
static double sideMax = TMath::Pi()/3 + 0.1;
|
| 45 |
|
| 46 |
static double sideCorrect = 1;
|
| 47 |
|
| 48 |
static const char* LUM = "#int L dt=150";
|
| 49 |
|
| 50 |
static bool plotSubtraction = 0;
|
| 51 |
|
| 52 |
static bool reweightCentrality = 1;
|
| 53 |
|
| 54 |
static const double pi = TMath::Pi();
|
| 55 |
|
| 56 |
static TCut lead = Form("pt1>%d",leadCut);
|
| 57 |
static TCut dijet = Form("pt1>%d && pt2>%d",leadCut,subleadCut);
|
| 58 |
|
| 59 |
static TCut deltaPhi = "abs(dphi)>2.0944";
|
| 60 |
static TCut side = Form("pt1>%d && pt2>%d && abs(dphi)>%f && abs(dphi)<%f",leadCut,subleadCut, sideMin, sideMax);
|
| 61 |
|
| 62 |
static TCut jetID = "trkMax1 > 4 || trkMax2 > 4";
|
| 63 |
static TCut noise = "noise < 0"; // Gplus Pile-Up cut also included
|
| 64 |
static TCut vtx = "abs(vz)<15";
|
| 65 |
|
| 66 |
TCut weight("weight*(pthat > 80)");
|
| 67 |
|
