| 1 |
yilmaz |
1.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 |
yilmaz |
1.2 |
|
| 29 |
yilmaz |
1.1 |
static bool subtract = 0;
|
| 30 |
|
|
|
| 31 |
|
|
static bool normLead = 0;
|
| 32 |
|
|
|
| 33 |
|
|
static int mixColor = 2;
|
| 34 |
|
|
static int dataColor = 1;
|
| 35 |
|
|
static int ppColor = 4;
|
| 36 |
|
|
|
| 37 |
|
|
static int centralBin = 8;
|
| 38 |
|
|
static int leadCut = 120;
|
| 39 |
|
|
static int subleadCut = 30;
|
| 40 |
|
|
|
| 41 |
|
|
static double sideMin = 0.1;
|
| 42 |
|
|
static double sideMax = TMath::Pi()/3 + 0.1;
|
| 43 |
|
|
|
| 44 |
|
|
static double sideCorrect = 1;
|
| 45 |
|
|
|
| 46 |
yilmaz |
1.3 |
static const char* LUM_PbPb = "L=150#mub^{-1}";
|
| 47 |
|
|
static const char* LUM_pPb = "L=17.3nb^{-1}";
|
| 48 |
yilmaz |
1.1 |
|
| 49 |
|
|
static bool plotSubtraction = 0;
|
| 50 |
|
|
|
| 51 |
|
|
static bool reweightCentrality = 1;
|
| 52 |
|
|
|
| 53 |
|
|
static const double pi = TMath::Pi();
|
| 54 |
|
|
|
| 55 |
|
|
static TCut lead = Form("pt1>%d",leadCut);
|
| 56 |
|
|
static TCut dijet = Form("pt1>%d && pt2>%d",leadCut,subleadCut);
|
| 57 |
|
|
|
| 58 |
|
|
static TCut deltaPhi = "abs(dphi)>2.0944";
|
| 59 |
|
|
static TCut side = Form("pt1>%d && pt2>%d && abs(dphi)>%f && abs(dphi)<%f",leadCut,subleadCut, sideMin, sideMax);
|
| 60 |
|
|
|
| 61 |
|
|
static TCut jetID = "trkMax1 > 4 || trkMax2 > 4";
|
| 62 |
|
|
static TCut noise = "noise < 0"; // Gplus Pile-Up cut also included
|
| 63 |
|
|
static TCut vtx = "abs(vz)<15";
|
| 64 |
|
|
|
| 65 |
yilmaz |
1.2 |
static TCut weight("weight*(pthat > 80)");
|
| 66 |
|
|
|
| 67 |
|
|
static TCut hfPlusBins[6] = {
|
| 68 |
|
|
"hfp>=30 && hfp<70",
|
| 69 |
|
|
"hfp>=20 && hfp<30",
|
| 70 |
|
|
"hfp>=15 && hfp<20",
|
| 71 |
|
|
"hfp>=10 && hfp<15",
|
| 72 |
|
|
"hfp>=5 && hfp<10",
|
| 73 |
|
|
"hfp<5"
|
| 74 |
|
|
};
|
| 75 |
|
|
|
| 76 |
|
|
// the same - no purpose
|
| 77 |
|
|
static TCut centBinsHFplus[6] = {
|
| 78 |
|
|
"cent>=30 && cent<70",
|
| 79 |
|
|
"cent>=20 && cent<30",
|
| 80 |
|
|
"cent>=15 && cent<20",
|
| 81 |
|
|
"cent>=10 && cent<15",
|
| 82 |
|
|
"cent>=5 && cent<10",
|
| 83 |
|
|
"cent<5"
|
| 84 |
|
|
};
|
| 85 |
yilmaz |
1.1 |
|
