Adding some histograms for jet shape

[u/mrichter/AliRoot.git] / PWG4 / JetTasks / AliAnalysisTaskJetSpectrum2.h

diff --git a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.h b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.h
index 90db179c365da464545381771987197d3424b6f5..89d8c2ccc86f1e7a33dfc5fb16b6861712d76b98 100644 (file)
--- a/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.h
+++ b/PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.h
@@ -50,6 +50,8 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     virtual void SetAODMCInput(Bool_t b){fUseAODMCInput = b;}
     virtual void SetLimitGenJetEta(Bool_t b){fLimitGenJetEta = b;}
     virtual void SetRecEtaWindow(Float_t f){fRecEtaWindow = f;}
+    virtual void SetMinJetPt(Float_t f){fMinJetPt = f;}
+    virtual void SetDeltaPhiWindow(Float_t f){fDeltaPhiWindow = f;}
     virtual void SetAnalysisType(Int_t i){fAnalysisType = i;}
     virtual void SetBranchGen(const char* c){fBranchGen = c;}
     virtual void SetBranchRec(const char* c){fBranchRec = c;}
@@ -93,7 +95,9 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     AliAODEvent  *fAOD; // where we take the jets from can be input or output AOD
     THnSparseF   *fhnJetContainer[kMaxStep*2];   // like particle container in corrfw with different steps need AliCFContainer with Scale(), and clone() to do the same
     THnSparseF   *fhnCorrelation;           // response matrix for unfolding 
+    THnSparseF   *fhnCorrelationPhiZRec;       // response matrix for unfolding in max Z rec bins
 
+    TF1          *f1PtScale;                // correction function to correct to the average true jet energy depending on p_T,rec
 
     TString       fBranchRec;  // AOD branch name for reconstructed
     TString       fBranchGen;  // AOD brnach for genereated
@@ -111,6 +115,9 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     Float_t       fAvgTrials;             // Average nimber of trials
     Float_t       fExternalWeight;        // external weight
     Float_t       fRecEtaWindow;          // eta window used for corraltion plots between rec and gen 
+    Float_t       fMinJetPt;              // limits the jet p_T in addition to what already is done in the jet finder, this is important for jet matching for JF with lo threshold
+    Float_t       fDeltaPhiWindow;        // minium angle between dijets
+
 
     TProfile*     fh1Xsec;   // pythia cross section and trials
     TH1F*         fh1Trials; // trials are added
@@ -122,6 +129,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TH1F*         fh1PtTrackRec;    // track pt
     TH1F*         fh1SumPtTrackRec; // sum over all track pT    
     TH1F*         fh1SumPtTrackAreaRec; // sum over all track pT    
+    TH1F*         fh1TmpRho;            // just temporary histo for calculation
 
 
     TH1F*         fh1PtRecIn[kMaxJets];  // Jet pt for all this info is also in the THNsparse      
@@ -132,7 +140,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TH1F*         fh1PtTracksRecIn;  // track pt for all tracks
     TH1F*         fh1PtTracksLeadingRecIn;  // track pt for all tracks
     TH1F*         fh1PtTracksGenIn;  // track pt for all tracks
-
+    
 
     TH2F*         fh2NRecJetsPt;            // Number of found jets above threshold
     TH2F*         fh2NRecTracksPt;          // Number of found tracks above threshold
@@ -141,20 +149,34 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TH2F*         fh2TracksLeadingPhiEta;   // track correlation with leading track
     TH2F*         fh2TracksLeadingPhiPt;    // track correlation with leading track
     TH2F*         fh2TracksLeadingJetPhiPt; // track correlation with leading track
+    TH2F*         fh2JetPtJetPhi;           // Phi distribution of accepted jets 
+    TH2F*         fh2TrackPtTrackPhi;       // phi distribution of accepted tracks
     TH2F*         fh2PhiPt[kMaxJets];    // delta phi correlation of tracks with the jet      
     TH2F*         fh2PhiEta[kMaxJets];   // eta   phi correlation of tracks with the jet      
-
+    TH2F*         fh2RhoPtRec[kMaxJets];    // jet shape variable rho
+    TH2F*         fh2PsiPtRec[kMaxJets];    // jet shape variable psi
+    TH2F*         fh2RhoPtGen[kMaxJets];    //
+    TH2F*         fh2PsiPtGen[kMaxJets];    //
     TH2F*         fh2FragRec[kMaxJets];     // fragmentation function
     TH2F*         fh2FragLnRec[kMaxJets];   // fragmetation in xi
-
     TH2F*         fh2FragGen[kMaxJets];     // fragmentation function
     TH2F*         fh2FragLnGen[kMaxJets];   // fragmetation in xi
 
+    // Dijet histos
+    TH2F*   fh2DijetDeltaPhiPt;      // dijet delta phi vs pt
+    TH2F*   fh2DijetAsymPt;          // dijet asym vs pt
+    TH2F*   fh2DijetAsymPtCut;       // dijet asym vs pt after delta phi cut
+    TH2F*   fh2DijetDeltaPhiDeltaEta; // dijet delta phi delta eta
+    TH2F*   fh2DijetPt2vsPt1;          // dijet pt2 vs pt1 
+    TH2F*   fh2DijetDifvsSum;          // dijet dif vs sum
+    TH1F*   fh1DijetMinv;            // dijet inv mass
+    TH1F*   fh1DijetMinvCut;         // dijet inv after delta phi cut
+
 
     TList *fHistList; // Output list
    
 
-    ClassDef(AliAnalysisTaskJetSpectrum2, 1) // Analysis task for standard jet analysis
+    ClassDef(AliAnalysisTaskJetSpectrum2, 4) // Analysis task for standard jet analysis
 };
  
 #endif