Adding histograms for background fluctuation with random cones

[u/mrichter/AliRoot.git] / JETAN / AliAnalysisTaskJetCluster.h

diff --git a/JETAN/AliAnalysisTaskJetCluster.h b/JETAN/AliAnalysisTaskJetCluster.h
index 3fbf805498335ecaf9482c46e8f7c29e713c5129..adc639abcd5ee823ef1b632abfc5655ad2a25563 100644 (file)
--- a/JETAN/AliAnalysisTaskJetCluster.h
+++ b/JETAN/AliAnalysisTaskJetCluster.h
@@ -19,17 +19,19 @@
 class AliJetHeader;
 class AliESDEvent;
 class AliAODEvent;
 class AliJetHeader;
 class AliESDEvent;
 class AliAODEvent;

+class AliAODExtension;

 class AliAODJet;
 class AliGenPythiaEventHeader;
 class AliCFManager;
 class AliAODJet;
 class AliGenPythiaEventHeader;
 class AliCFManager;

-
+class AliAODJetEventBackground;
+class AliJetFinder;

 class TList;
 class TChain;
 class TH2F;
 class TH1F;
 class TH3F;
 class TProfile;
 class TList;
 class TChain;
 class TH2F;
 class TH1F;
 class TH3F;
 class TProfile;

-
+class TRefArray;

 
 
 class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
 
 
 class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE


@@ -37,7 +39,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
  public:
     AliAnalysisTaskJetCluster();
     AliAnalysisTaskJetCluster(const char* name);
  public:
     AliAnalysisTaskJetCluster();
     AliAnalysisTaskJetCluster(const char* name);

-    virtual ~AliAnalysisTaskJetCluster() {;}
+    virtual ~AliAnalysisTaskJetCluster();

     // Implementation of interface methods
     virtual void UserCreateOutputObjects();
     virtual void Init();
     // Implementation of interface methods
     virtual void UserCreateOutputObjects();
     virtual void Init();


@@ -46,15 +48,22 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     virtual void Terminate(Option_t *option);
     virtual Bool_t Notify();
 
     virtual void Terminate(Option_t *option);
     virtual Bool_t Notify();
 

+    
+

     virtual void SetUseGlobalSelection(Bool_t b){fUseGlobalSelection = b;}
     virtual void SetAODTrackInput(Bool_t b){fUseAODTrackInput = b;}
     virtual void SetAODMCInput(Bool_t b){fUseAODMCInput = b;}
     virtual void SetRecEtaWindow(Float_t f){fRecEtaWindow = f;}
     virtual void SetTrackTypeGen(Int_t i){fTrackTypeGen = i;}
     virtual void SetTrackTypeRec(Int_t i){fTrackTypeRec = i;}
     virtual void SetUseGlobalSelection(Bool_t b){fUseGlobalSelection = b;}
     virtual void SetAODTrackInput(Bool_t b){fUseAODTrackInput = b;}
     virtual void SetAODMCInput(Bool_t b){fUseAODMCInput = b;}
     virtual void SetRecEtaWindow(Float_t f){fRecEtaWindow = f;}
     virtual void SetTrackTypeGen(Int_t i){fTrackTypeGen = i;}
     virtual void SetTrackTypeRec(Int_t i){fTrackTypeRec = i;}

+    virtual void SetTrackPtCut(Float_t x){fTrackPtCut = x;}

     virtual void SetFilterMask(UInt_t i){fFilterMask = i;}
     virtual void SetFilterMask(UInt_t i){fFilterMask = i;}

-    // use for the CF

     
     

+    virtual void SetNSkipLeadingRan(Int_t x){fNSkipLeadingRan = x;}
+
+    virtual void SetJetOutputBranch(const char *c){fNonStdBranch = c;}
+    virtual void SetJetOutputFile(const char *c){fNonStdFile = c;}
+    virtual void SetJetOutputMinPt(Float_t x){fJetOutputMinPt = x;}

 
     // for Fast Jet
     fastjet::JetAlgorithm        GetAlgorithm()         const {return fAlgorithm;}
 
     // for Fast Jet
     fastjet::JetAlgorithm        GetAlgorithm()         const {return fAlgorithm;}


@@ -67,6 +76,11 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     void SetStrategy(fastjet::Strategy f)                {fStrategy = f;}
     void SetRecombScheme(fastjet::RecombinationScheme f) {fRecombScheme = f;}
     void SetAreaType(fastjet::AreaType f)                {fAreaType = f;}
     void SetStrategy(fastjet::Strategy f)                {fStrategy = f;}
     void SetRecombScheme(fastjet::RecombinationScheme f) {fRecombScheme = f;}
     void SetAreaType(fastjet::AreaType f)                {fAreaType = f;}

+    void SetGhostArea(Double_t f) {fGhostArea = f;}
+    void SetActiveAreaRepeats(Int_t f) {fActiveAreaRepeats = f;}
+    void SetGhostEtamax(Double_t f) {fGhostEtamax = f;}
+
+

 
     // Helper
     //
 
     // Helper
     //


@@ -88,77 +102,97 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
 
     Int_t GetListOfTracks(TList *list,Int_t type);
 
 
     Int_t GetListOfTracks(TList *list,Int_t type);
 

-    AliAODEvent  *fAOD; // where we take the jets from can be input or output AOD
-
+    AliAODEvent     *fAOD;                // ! where we take the jets from can be input or output AOD
+    AliAODExtension *fAODExtension;       // ! AOD extension in case we write a non-sdt branch to a separate file and the aod is standard
+    TRefArray       *fRef;               // ! trefarray for track references within the jet

     Bool_t        fUseAODTrackInput;      // take track from input AOD not from ouptu AOD
     Bool_t        fUseAODMCInput;         // take MC from input AOD not from ouptu AOD
     Bool_t        fUseGlobalSelection;    // Limit the eta of the generated jets
     Bool_t        fUseAODTrackInput;      // take track from input AOD not from ouptu AOD
     Bool_t        fUseAODMCInput;         // take MC from input AOD not from ouptu AOD
     Bool_t        fUseGlobalSelection;    // Limit the eta of the generated jets

-    UInt_t        fFilterMask;             // filter bit for slecected tracks
+    UInt_t        fFilterMask;            // filter bit for slecected tracks

     Int_t         fTrackTypeRec;          // type of tracks used for FF 
     Int_t         fTrackTypeGen;          // type of tracks used for FF 
     Int_t         fTrackTypeRec;          // type of tracks used for FF 
     Int_t         fTrackTypeGen;          // type of tracks used for FF 

+    Int_t         fNSkipLeadingRan;        // number of leading tracks to be skipped in the randomized event

     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       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       fTrackPtCut;            // minimum track pt to be accepted
+    Float_t       fJetOutputMinPt;        // minimum p_t for jets to be written out
+
+    // output configurartion
+    TString       fNonStdBranch;      // the name of the non-std branch name, if empty no branch is filled
+    TString       fNonStdFile;        // The optional name of the output file the non-std brnach is written to
+    
+

     // Fast jet
     Double_t fRparam;
     fastjet::JetAlgorithm fAlgorithm; //fastjet::kt_algorithm
     fastjet::Strategy fStrategy;  //= fastjet::Best;
     fastjet::RecombinationScheme fRecombScheme; // = fastjet::BIpt_scheme;
     fastjet::AreaType fAreaType; 
     // Fast jet
     Double_t fRparam;
     fastjet::JetAlgorithm fAlgorithm; //fastjet::kt_algorithm
     fastjet::Strategy fStrategy;  //= fastjet::Best;
     fastjet::RecombinationScheme fRecombScheme; // = fastjet::BIpt_scheme;
     fastjet::AreaType fAreaType; 

-
-    TProfile*     fh1Xsec;   // pythia cross section and trials
-    TH1F*         fh1Trials; // trials are added
-    TH1F*         fh1PtHard;  // Pt har of the event...       
-    TH1F*         fh1PtHardNoW;  // Pt har of the event without weigt       
-    TH1F*         fh1PtHardTrials;  // Number of trials 
-
-    TH1F*         fh1NJetsRec; // number of reconstructed jets
-    TH1F*         fh1NConstRec;// number of constiutens in leading jet
-    TH1F*         fh1NConstLeadingRec;// number of constiutens in leading jet
-    TH1F*         fh1PtJetsRecIn;  // Jet pt for all jets
-    TH1F*         fh1PtJetsLeadingRecIn;  // Jet pt for all jets
-    TH1F*         fh1PtJetConstRec;// pt of constituents
+    Double_t fGhostArea;
+    Int_t fActiveAreaRepeats;
+    Double_t fGhostEtamax;
+    TProfile*     fh1Xsec;   //! pythia cross section and trials
+    TH1F*         fh1Trials; //! trials are added
+    TH1F*         fh1PtHard;  //! Pt har of the event...       
+    TH1F*         fh1PtHardNoW;  //! Pt har of the event without weigt       
+    TH1F*         fh1PtHardTrials;  //! Number of trials 
+
+    TH1F*         fh1NJetsRec; //! number of reconstructed jets
+    TH1F*         fh1NConstRec;//! number of constiutens in leading jet
+    TH1F*         fh1NConstLeadingRec;//! number of constiutens in leading jet
+    TH1F*         fh1PtJetsRecIn;  //! Jet pt for all jets
+    TH1F*         fh1PtJetsLeadingRecIn;  //! Jet pt for all jets
+    TH1F*         fh1PtJetConstRec;//! pt of constituents

     TH1F*         fh1PtJetConstLeadingRec;// pt of constituents
     TH1F*         fh1PtJetConstLeadingRec;// pt of constituents

-    TH1F*         fh1PtTracksRecIn;  // track pt for all tracks
-    TH1F*         fh1PtTracksLeadingRecIn;  // track pt for all tracks
+    TH1F*         fh1PtTracksRecIn;  //! track pt for all tracks
+    TH1F*         fh1PtTracksLeadingRecIn;  //! track pt for all tracks

 
     // Randomized track histos
 
     // Randomized track histos

-    TH1F*         fh1NJetsRecRan; // number of reconstructed jets from randomized
-    TH1F*         fh1NConstRecRan;// number of constiutens in leading jet
-    TH1F*         fh1PtJetsLeadingRecInRan;  // Jet pt for all jets
-    TH1F*         fh1NConstLeadingRecRan;// number of constiutens in leading jet
-    TH1F*         fh1PtJetsRecInRan;  // Jet pt for all jets
-
-    TH1F*         fh1PtTracksGenIn;  // track pt for all tracks
-
-
-    TH2F*         fh2NRecJetsPt;            // Number of found jets above threshold
-    TH2F*         fh2NRecTracksPt;          // Number of found tracks above threshold
-    TH2F*         fh2NConstPt;           // number of constituents vs. pt
-    TH2F*         fh2NConstLeadingPt;           // number of constituents vs. pt
-    TH2F*         fh2JetPhiEta;             // jet phi eta
-    TH2F*         fh2LeadingJetPhiEta;      // leading jet phi eta
-    TH2F*         fh2JetEtaPt;              // leading jet eta
-    TH2F*         fh2LeadingJetEtaPt;              // leading jet eta
-    TH2F*         fh2TrackEtaPt;              // track eta
-    TH2F*         fh2LeadingTrackEtaPt;       // leading track eta
-    TH2F*         fh2JetsLeadingPhiEta;     // jet phi eta
-    TH2F*         fh2JetsLeadingPhiPt;      // jet correlation with leading jet
-    TH2F*         fh2TracksLeadingPhiEta;   // track correlation with leading track
-    TH2F*         fh2TracksLeadingPhiPt;    // track correlation with leading track
-    TH2F*         fh2TracksLeadingJetPhiPt; // track correlation with leading Jet
-    TH2F*         fh2JetsLeadingPhiPtW;      // jet correlation with leading jet
-    TH2F*         fh2TracksLeadingPhiPtW;   // track correlation with leading track
-    TH2F*         fh2TracksLeadingJetPhiPtW; // track correlation with leading Jet
-    TH2F*         fh2NRecJetsPtRan;            // Number of found jets above threshold
-    TH2F*         fh2NConstPtRan;           // number of constituents vs. pt
-    TH2F*         fh2NConstLeadingPtRan;           // number of constituents vs. pt
-    TH2F*         fh2TracksLeadingJetPhiPtRan; // track correlation with leading Jet
-    TH2F*         fh2TracksLeadingJetPhiPtWRan; // track correlation with leading Jet
+    TH1F*         fh1NJetsRecRan; //! number of reconstructed jets from randomized
+    TH1F*         fh1NConstRecRan;//! number of constiutens in leading jet
+    TH1F*         fh1PtJetsLeadingRecInRan;  //! Jet pt for all jets
+    TH1F*         fh1NConstLeadingRecRan;//! number of constiutens in leading jet
+    TH1F*         fh1PtJetsRecInRan;  //! Jet pt for all jets
+
+    TH1F*         fh1PtTracksGenIn;  //! track pt for all tracks
+    TH1F*         fh1Nch;            //! charged particle mult
+    TH1F*         fh1BiARandomCones[3]; //! Residual distribtion from reandom cones on real event
+    TH1F*         fh1BiARandomConesRan[3]; //! Residual distribtion from reandom cones on random event
+
+
+    TH2F*         fh2NRecJetsPt;            //! Number of found jets above threshold
+    TH2F*         fh2NRecTracksPt;          //! Number of found tracks above threshold
+    TH2F*         fh2NConstPt;           //! number of constituents vs. pt
+    TH2F*         fh2NConstLeadingPt;           //! number of constituents vs. pt
+    TH2F*         fh2JetPhiEta;             //! jet phi eta
+    TH2F*         fh2LeadingJetPhiEta;      //! leading jet phi eta
+    TH2F*         fh2JetEtaPt;              //! leading jet eta
+    TH2F*         fh2LeadingJetEtaPt;              //! leading jet eta
+    TH2F*         fh2TrackEtaPt;              //! track eta
+    TH2F*         fh2LeadingTrackEtaPt;       //! leading track eta
+    TH2F*         fh2JetsLeadingPhiEta;     //! jet phi eta
+    TH2F*         fh2JetsLeadingPhiPt;      //! jet correlation with leading jet
+    TH2F*         fh2TracksLeadingPhiEta;   //! track correlation with leading track
+    TH2F*         fh2TracksLeadingPhiPt;    //! track correlation with leading track
+    TH2F*         fh2TracksLeadingJetPhiPt; //! track correlation with leading Jet
+    TH2F*         fh2JetsLeadingPhiPtW;      //! jet correlation with leading jet
+    TH2F*         fh2TracksLeadingPhiPtW;   //! track correlation with leading track
+    TH2F*         fh2TracksLeadingJetPhiPtW; //! track correlation with leading Jet
+    TH2F*         fh2NRecJetsPtRan;            //! Number of found jets above threshold
+    TH2F*         fh2NConstPtRan;           //! number of constituents vs. pt
+    TH2F*         fh2NConstLeadingPtRan;           //! number of constituents vs. pt
+    TH2F*         fh2PtNch;               //! p_T of cluster vs. multiplicity,
+    TH2F*         fh2PtNchRan;            //! p_T of cluster vs. multiplicity,random
+    TH2F*         fh2PtNchN;               //! p_T of cluster vs. multiplicity, weigthed with constituents
+    TH2F*         fh2PtNchNRan;            //! p_T of cluster vs. multiplicity, weigthed with constituents random
+    TH2F*         fh2TracksLeadingJetPhiPtRan; //! track correlation with leading Jet
+    TH2F*         fh2TracksLeadingJetPhiPtWRan; //! track correlation with leading Jet
+

     TList *fHistList; // Output list
    
 
     TList *fHistList; // Output list
    
 

-    ClassDef(AliAnalysisTaskJetCluster, 2) // Analysis task for standard jet analysis
+    ClassDef(AliAnalysisTaskJetCluster, 6) 

 };
  
 #endif
 };
  
 #endif