1 /* Copyright(c) 1998-2014, ALICE Experiment at CERN, All rights reserved. *
2 * See cxx source for full Copyright notice */
4 // Short comment describing what this class does needed!
6 //===========================================================
7 // Dummy comment, should be replaced by a real one
8 //===========================================================
10 #ifndef ALIJJETJTANALYSIS_H
11 #define ALIJJETJTANALYSIS_H
14 #include <TObjArray.h>
18 #include "AliJHistManager.h"
19 #include "AliJJetAnalysis.h"
24 class AliJJetJtAnalysis{
29 AliJJetJtAnalysis( AliJCard * card );
30 AliJJetJtAnalysis(const AliJJetJtAnalysis& ap);
31 AliJJetJtAnalysis& operator = (const AliJJetJtAnalysis& ap);
34 //type = charged jet or full jet, js = cone size
36 void FillHistosJets();
38 void AddJets(TObjArray * jets ){
42 fJetListOfList.Add( (TObject*)jets );
44 for( int i=0;i<jets->GetEntriesFast();i++ ){
45 //((AliJJet*)jets->At(i))->ReSum();
47 } // TODO clean before event
49 void SetJTracks(TClonesArray *tracks ){fTracks = tracks ;}
52 int GetNJets(){ return GetJetList()->GetEntriesFast(); }
53 TObjArray* GetJetList(){ return fJetList; }
54 //Double_t GetJetEtaRange(){ return fJetEtaRange; }
55 //void SetJetEtaRange(double eta){ fJetEtaRange=eta; }
56 void SetJetList(TObjArray* jetlist){ fJetList=jetlist; }
57 void SetInputList(TObjArray * ilist){ fInputList = ilist;}
58 //void SetTrackJetMap(std::vector<int> * v){ fTrackJetMap=v;}
59 //void SetJetPtBin( TVector * b){ fJetPtBins=b; }
60 void SetCard(AliJCard * card) {fCard = card;}
62 AliJJet & Jet( int i ){ return *(AliJJet*) fJetList->At(i); }
64 void UserCreateOutputObjects();
67 void ClearBeforeEvent();
68 void FillJtHistogram( TObjArray *Jets, int iContainer );
70 void FillBgJtWithSmallerR(const TClonesArray &Jets,int iContainer,
71 double nR, int iHist);
73 void FillBgJtWithDiffAxes (int iao, int ia, int iHist);
75 void WriteHistograms();
77 int GetBin(TVector *array, double val){
80 for(int i=1; i< array->GetNoElements(); i++){
81 if((*array)[i] <= val && val<(*array)[i+1]){
89 void SetJetFinderName(vector<TString> JetFinderName){ fJetFinderName = JetFinderName; }
91 // Need for event loop
92 void SetCentralityBin( int cbin) { cBin = cbin;}
93 void SetCentrality( float cent) { fcent = cent;}
94 void SetZVertex( float zvtx) { zVert = zvtx;}
95 void SetZVertexBin( int zbin) { zBin = zbin;}
96 void SetNumberOfJetFinders( int njfinder ) { nJetContainer = njfinder;}
97 AliJEfficiency* GetAliJEfficiency() { return fEfficiency;}
102 TObjArray * fInputList; // comment needed
103 TObjArray * fJetList; // comment needed
104 TObjArray fJetListOfList; // !comment needed
105 //TObjArray * fJetBgList; // comment needed
106 //TObjArray fJetBgListOfList; // !comment needed
107 vector<TClonesArray> fJetBgListOfList;
109 TVector *fJetTriggPtBorders;
110 TVector *fJetConstPtLowLimits;
111 TVector *fJetAssocPtBorders;
112 TVector *fDeltaRBorders;
115 AliJCard * fCard; // comment needed
116 AliJJetAnalysis *fJJetAnalysis;
117 vector<TString> fJetFinderName;
118 vector<double> fConeSizes;
120 AliJEfficiency *fEfficiency;
125 TClonesArray *fTracks;
128 AliJHistManager * fHMG;
130 AliJBin fJetFinderBin;
131 AliJBin fJetTriggerBin;
133 AliJBin fTrkLimPtBin;
144 AliJTH1D fhJtWeightBin;
145 AliJTH1D fhLogJtWeightBin;
146 AliJTH1D fhJtWithPtCutWeightBinBin;
147 AliJTH1D fhLogJtWithPtCutWeightBinBin;
148 AliJTH1D fhJtBinLimBin;
149 AliJTH1D fhJtWeightBinLimBin;
150 AliJTH1D fhLogJtWeightBinLimBin;
154 AliJTH1D fhJetBgPtBin;
159 AliJTH1D fhBgJtWeightBin;
160 AliJTH1D fhBgLogJtWeightBin;
161 AliJTH1D fhBgJtWithPtCutWeightBinBin;
162 AliJTH1D fhBgLogJtWithPtCutWeightBinBin;
163 AliJTH1D fhBgJtWithPtCutWeightBinBinSmallerR;
164 AliJTH1D fhBgLogJtWithPtCutWeightBinBinSmallerR;
165 AliJTH1D fhBgJtWithPtCutWeightBinBinDiffR;
166 AliJTH1D fhBgLogJtWithPtCutWeightBinBinDiffR;
167 AliJTH1D fhBgJtBinLimBin;
168 AliJTH1D fhBgJtWeightBinLimBin;
169 AliJTH1D fhBgLogJtWeightBinLimBin;
173 AliJTH1D fhFullJetEChJetBin;
174 AliJTH1D fhFullChdRChJetBin;
175 AliJTH2D fh2DFullEvsChEdN0;
176 AliJTH2D fh2DFullEvsChEdNnot0;
179 //double fJetPtMinCut;