- fh1FFIMK0ConeSmear = new TH1F("fh1FFIMK0ConeSmear","Smeared jet pt study for K0s-in-cone-jets; smeared jet #it{p}_{T}", 19,5.,100.);
- fh1FFIMLaConeSmear = new TH1F("fh1FFIMLaConeSmear","Smeared jet pt study for La-in-cone-jets; smeared jet #it{p}_{T}", 19,5.,100.);
- fh1FFIMALaConeSmear = new TH1F("fh1FFIMALaConeSmear","Smeared jet pt study for ALa-in-cone-jets; smeared jet #it{p}_{T}", 19,5.,100.);
-
- fh3MCrecK0Cone = new TH3F("fh3MCrecK0Cone", "MC rec {K^{0}}^{s} #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",19,5.,100., 400,0.3,0.7, 200,0.,20.);
- fh3MCrecLaCone = new TH3F("fh3MCrecLaCone", "MC rec {#Lambda #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2});#it{p}_{T}",19,5.,100., 200,1.05,1.25, 200,0.,20.);
- fh3MCrecALaCone = new TH3F("fh3MCrecALaCone", "MC rec {#Antilambda #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2});#it{p}_{T}",19,5.,100.,200,1.05,1.25, 200,0.,20.);
- fh3MCrecK0ConeSmear = new TH3F("fh3MCrecK0ConeSmear", "MC rec {K^{0}}^{s} #it{p}_{T} in cone around jet axis matching MC gen particle, with jet p_{T} smeared; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",19,5.,100., 400,0.3,0.7, 200,0.,20.);
- fh3MCrecLaConeSmear = new TH3F("fh3MCrecLaConeSmear", "MC rec {#Lambda #it{p}_{T} in cone around jet axis matching MC gen particle, with jet p_{T} smeared; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2});#it{p}_{T}",19,5.,100., 200,1.05,1.25, 200,0.,20.);
- fh3MCrecALaConeSmear = new TH3F("fh3MCrecALaConeSmear", "MC rec {#Antilambda #it{p}_{T} in cone around jet axis matching MC gen particle, with jet p_{T} smeared; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2});#it{p}_{T}",19,5.,100.,200,1.05,1.25, 200,0.,20.);
- fh3SecContinCone = new TH3F("fh3SecContinCone","secondary contamination of jet cones; jet #it{p}_{T}; track #it{p}_{T}, #eta",19,5.,100.,200,0.,20.,200,-1.,1.);
- fh3StrContinCone = new TH3F("fh3StrContinCone","strange particle contamination of jet cones; jet #it{p}_{T}; track #it{p}_{T}, #eta",19,5.,100.,200,0.,20.,200,-1.,1.);
+ fh1IMK0ConeSmear = new TH1F("fh1IMK0ConeSmear","Smeared jet pt study for K0s-in-cone-jets; smeared jet #it{p}_{T}", 19,5.,100.);
+ fh1IMLaConeSmear = new TH1F("fh1IMLaConeSmear","Smeared jet pt study for La-in-cone-jets; smeared jet #it{p}_{T}", 19,5.,100.);
+ fh1IMALaConeSmear = new TH1F("fh1IMALaConeSmear","Smeared jet pt study for ALa-in-cone-jets; smeared jet #it{p}_{T}", 19,5.,100.);
+
+ //8 new histograms: Cone, Incl, Lambda, Antilambda, Hijing, Injected:
+
+ Int_t binsrecMCHijingLaIncl[3] = {200, 200, 200};
+ Double_t xminrecMCHijingLaIncl[3] = {1.05, 0., -1.};
+ Double_t xmaxrecMCHijingLaIncl[3] = {1.25, 20., 1.};
+ fhnrecMCHijingLaIncl = new THnSparseF("fhnrecMCHijingLaIncl","La inv. mass; particle pT; particle #eta",3,binsrecMCHijingLaIncl,xminrecMCHijingLaIncl,xmaxrecMCHijingLaIncl);
+
+ Int_t binsrecMCHijingLaCone[4] = {19, 200, 200, 200};
+ Double_t xminrecMCHijingLaCone[4] = {5., 1.05, 0., -1.};
+ Double_t xmaxrecMCHijingLaCone[4] = {100., 1.25, 20., 1.};
+ fhnrecMCHijingLaCone = new THnSparseF("fhnrecMCHijingLaCone","La inv. mass; particle pT; particle #eta",4,binsrecMCHijingLaCone,xminrecMCHijingLaCone,xmaxrecMCHijingLaCone);
+
+ Int_t binsrecMCHijingALaIncl[3] = {200, 200, 200};
+ Double_t xminrecMCHijingALaIncl[3] = {1.05, 0., -1.};
+ Double_t xmaxrecMCHijingALaIncl[3] = {1.25, 20., 1.};
+ fhnrecMCHijingALaIncl = new THnSparseF("fhnrecMCHijingALaIncl","ALa inv. mass; particle pT; particle #eta",3,binsrecMCHijingALaIncl,xminrecMCHijingALaIncl,xmaxrecMCHijingALaIncl);
+
+ Int_t binsrecMCHijingALaCone[4] = {19, 200, 200, 200};
+ Double_t xminrecMCHijingALaCone[4] = {5., 1.05, 0., -1.};
+ Double_t xmaxrecMCHijingALaCone[4] = {100., 1.25, 20., 1.};
+ fhnrecMCHijingALaCone = new THnSparseF("fhnrecMCHijingALaCone","ALa inv. mass; particle pT; particle #eta",4,binsrecMCHijingALaCone,xminrecMCHijingALaCone,xmaxrecMCHijingALaCone);
+
+ Int_t binsrecMCInjectLaIncl[3] = {200, 200, 200};
+ Double_t xminrecMCInjectLaIncl[3] = {1.05, 0., -1.};
+ Double_t xmaxrecMCInjectLaIncl[3] = {1.25, 20., 1.};
+ fhnrecMCInjectLaIncl = new THnSparseF("fhnrecMCInjectLaIncl","La inv. mass; particle pT; particle #eta",3,binsrecMCInjectLaIncl,xminrecMCInjectLaIncl,xmaxrecMCInjectLaIncl);
+
+ Int_t binsrecMCInjectLaCone[4] = {19, 200, 200, 200};
+ Double_t xminrecMCInjectLaCone[4] = {5., 1.05, 0., -1.};
+ Double_t xmaxrecMCInjectLaCone[4] = {100., 1.25, 20., 1.};
+ fhnrecMCInjectLaCone = new THnSparseF("fhnrecMCInjectLaCone","La jet pT;inv. mass; particle pT; particle #eta",4,binsrecMCInjectLaCone,xminrecMCInjectLaCone,xmaxrecMCInjectLaCone);
+
+ Int_t binsrecMCInjectALaIncl[3] = {200, 200, 200};
+ Double_t xminrecMCInjectALaIncl[3] = {1.05, 0., -1.};
+ Double_t xmaxrecMCInjectALaIncl[3] = {1.25, 20., 1.};
+ fhnrecMCInjectALaIncl = new THnSparseF("fhnrecMCInjectALaIncl","ALa inv. mass; particle pT; particle #eta",3,binsrecMCInjectALaIncl,xminrecMCInjectALaIncl,xmaxrecMCInjectALaIncl);
+
+ Int_t binsrecMCInjectALaCone[4] = {19, 200, 200, 200};
+ Double_t xminrecMCInjectALaCone[4] = {5., 1.05, 0., -1.};
+ Double_t xmaxrecMCInjectALaCone[4] = {100., 1.25, 20., 1.};
+ fhnrecMCInjectALaCone = new THnSparseF("fhnrecMCInjectALaCone","ALa inv. mass; particle pT; particle #eta",4,binsrecMCInjectALaCone,xminrecMCInjectALaCone,xmaxrecMCInjectALaCone);
+
+
+ Int_t binsMCrecK0Cone[4] = {19, 200, 200, 200};
+ Double_t xminMCrecK0Cone[4] = {5.,0.3, 0., -1.};
+ Double_t xmaxMCrecK0Cone[4] = {100.,0.7, 20., 1.};
+ fhnMCrecK0Cone = new THnSparseF("fhnMCrecK0Cone", "MC rec {K^{0}}^{s} #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",4,binsMCrecK0Cone,xminMCrecK0Cone,xmaxMCrecK0Cone);
+
+ Int_t binsMCrecLaCone[4] = {19, 200, 200, 200};
+ Double_t xminMCrecLaCone[4] = {5.,0.3, 0., -1.};
+ Double_t xmaxMCrecLaCone[4] = {100.,0.7, 20., 1.};
+ fhnMCrecLaCone = new THnSparseF("fhnMCrecLaCone", "MC rec {#Lambda #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",4,binsMCrecLaCone,xminMCrecLaCone,xmaxMCrecLaCone);
+
+ Int_t binsMCrecALaCone[4] = {19, 200, 200, 200};
+ Double_t xminMCrecALaCone[4] = {5.,0.3, 0., -1.};
+ Double_t xmaxMCrecALaCone[4] = {100.,0.7, 20., 1.};
+ fhnMCrecALaCone = new THnSparseF("fhnMCrecALaCone", "MC rec {#bar{#Lambda} #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",4,binsMCrecALaCone,xminMCrecALaCone,xmaxMCrecALaCone);
+
+ Int_t binsMCrecK0ConeSmear[4] = {19, 200, 200, 200};
+ Double_t xminMCrecK0ConeSmear[4] = {5.,0.3, 0., -1.};
+ Double_t xmaxMCrecK0ConeSmear[4] = {100.,0.7, 20., 1.};
+ fhnMCrecK0ConeSmear = new THnSparseF("fhnMCrecK0ConeSmear", "MC rec {K^{0}}^{s} #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",4,binsMCrecK0ConeSmear,xminMCrecK0ConeSmear,xmaxMCrecK0ConeSmear);
+
+ Int_t binsMCrecLaConeSmear[4] = {19, 200, 200, 200};
+ Double_t xminMCrecLaConeSmear[4] = {5.,0.3, 0., -1.};
+ Double_t xmaxMCrecLaConeSmear[4] = {100.,0.7, 20., 1.};
+ fhnMCrecLaConeSmear = new THnSparseF("fhnMCrecLaConeSmear", "MC rec {#Lambda #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",4,binsMCrecLaConeSmear,xminMCrecLaConeSmear,xmaxMCrecLaConeSmear);
+
+ Int_t binsMCrecALaConeSmear[4] = {19, 200, 200, 200};
+ Double_t xminMCrecALaConeSmear[4] = {5.,0.3, 0., -1.};
+ Double_t xmaxMCrecALaConeSmear[4] = {100.,0.7, 20., 1.};
+ fhnMCrecALaConeSmear = new THnSparseF("fhnMCrecALaConeSmear", "MC rec {#bar{#Lambda} #it{p}_{T} in cone around jet axis matching MC gen particle; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",4,binsMCrecALaConeSmear,xminMCrecALaConeSmear,xmaxMCrecALaConeSmear);
+
+ Int_t binsK0sSecContinCone[3] = {19, 200, 200};
+ Double_t xminK0sSecContinCone[3] = {5.,0., -1.};
+ Double_t xmaxK0sSecContinCone[3] = {100.,20., 1.};
+ fhnK0sSecContinCone = new THnSparseF("fhnK0sSecContinCone", "Secondary contamination {K^{0}}^{s} #it{p}_{T} in cone around jet axis; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",3,binsK0sSecContinCone,xminK0sSecContinCone,xmaxK0sSecContinCone);
+
+ Int_t binsLaSecContinCone[3] = {19, 200, 200};
+ Double_t xminLaSecContinCone[3] = {5.,0., -1.};
+ Double_t xmaxLaSecContinCone[3] = {100.,20., 1.};
+ fhnLaSecContinCone = new THnSparseF("fhnLaSecContinCone", "Secondary contamination {#Lambda #it{p}_{T} in cone around jet axis; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",3,binsLaSecContinCone,xminLaSecContinCone,xmaxLaSecContinCone);
+
+ Int_t binsALaSecContinCone[3] = {19, 200, 200};
+ Double_t xminALaSecContinCone[3] = {5.,0., -1.};
+ Double_t xmaxALaSecContinCone[3] = {100.,20., 1.};
+ fhnALaSecContinCone = new THnSparseF("fhnALaSecContinCone", "Secondary contamination {#bar{#Lambda} #it{p}_{T} in cone around jet axis; jet #it{p}_{T}; inv mass (GeV/#it{c}^{2};#it{p}_{T}",3,binsALaSecContinCone,xminALaSecContinCone,xmaxALaSecContinCone);
+
+ Int_t binsK0sIncl[3] = {200, 200, 200};
+ Double_t xminK0sIncl[3] = {0.3, 0., -1.};
+ Double_t xmaxK0sIncl[3] = {0.7, 20., 1.};
+ fhnK0sIncl = new THnSparseF("fhnK0sIncl","K0s inv. mass; particle pT; particle #eta",3,binsK0sIncl,xminK0sIncl,xmaxK0sIncl);
+
+ Int_t binsK0sCone[4] = {19, 200, 200, 200};
+ Double_t xminK0sCone[4] = {5.,0.3, 0., -1.};
+ Double_t xmaxK0sCone[4] = {100.,0.7, 20., 1.};
+ fhnK0sCone = new THnSparseF("fhnK0sCone","jet pT; K0s inv. mass; particle pT; particle #eta",4,binsK0sCone,xminK0sCone,xmaxK0sCone);
+
+ Int_t binsLaIncl[3] = {200, 200, 200};
+ Double_t xminLaIncl[3] = {1.05, 0., -1.};
+ Double_t xmaxLaIncl[3] = {1.25, 20., 1.};
+ fhnLaIncl = new THnSparseF("fhnLaIncl","La inv. mass; particle pT; particle #eta",3,binsLaIncl,xminLaIncl,xmaxLaIncl);
+
+ Int_t binsLaCone[4] = {19, 200, 200, 200};
+ Double_t xminLaCone[4] = {5.,1.05, 0., -1.};
+ Double_t xmaxLaCone[4] = {100.,1.25, 20., 1.};
+ fhnLaCone = new THnSparseF("fhnLaCone","jet pT; La inv. mass; particle pT; particle #eta",4,binsLaCone,xminLaCone,xmaxLaCone);
+
+ Int_t binsALaIncl[3] = {200, 200, 200};
+ Double_t xminALaIncl[3] = {1.05, 0., -1.};
+ Double_t xmaxALaIncl[3] = {1.25, 20., 1.};
+ fhnALaIncl = new THnSparseF("fhnALaIncl","ALa inv. mass; particle pT; particle #eta",3,binsALaIncl,xminALaIncl,xmaxALaIncl);
+
+ Int_t binsALaCone[4] = {19, 200, 200, 200};
+ Double_t xminALaCone[4] = {5.,1.05, 0., -1.};
+ Double_t xmaxALaCone[4] = {100.,1.25, 20., 1.};
+ fhnALaCone = new THnSparseF("fhnALaCone","jet pT; ALa inv. mass; particle pT; particle #eta",4,binsALaCone,xminALaCone,xmaxALaCone);