[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / UserTasks / AliAnalysisTaskEmcalDiJetResponse.h

#ifndef ALIANALYSISTASKEMCALDIJETRESPONSE_H
#define ALIANALYSISTASKEMCALDIJETRESPONSE_H

class TH1;
class TH2;
class TH3;
class TH3F;
class THnSparse;
class TClonesArray;
class TArrayI;
class AliAnalysisUtils;
class AliAnalysisManager;
class AliGenPythiaEventHeader;

#include "AliJetContainer.h"

#include "AliAnalysisTaskEmcalDiJetBase.h"

class AliAnalysisTaskEmcalDiJetResponse : public AliAnalysisTaskEmcalDiJetBase {
 public:
  AliAnalysisTaskEmcalDiJetResponse();
  AliAnalysisTaskEmcalDiJetResponse(const char *name);
  virtual ~AliAnalysisTaskEmcalDiJetResponse();

  void                        UserCreateOutputObjects();
  void                        Terminate(Option_t *option);

  //Setters
  void SetMatchFullCharged(Bool_t b)        { fDoMatchFullCharged = b; }
  void SetResponseVar(Int_t v)              { fnUsedResponseVar   = v; }    

  //Getters

 protected:
  Bool_t                      Run()              ;
  void                        CorrelateJets(const Int_t type);
  void                        CorrelateAllJets(const Int_t type);
  void                        CorrelateTwoJets(const Int_t type);

  Bool_t                      FillHistograms()   ;
  void                        FillDiJetHistos(const AliEmcalJet *jet1 = 0, const AliEmcalJet *jet2 = 0, const Int_t mode = 0);
  void                        FillMatchHistos();
  Bool_t                      RetrieveEventObjects();

  void                        FillDiJetResponse(const AliEmcalJet *jetTrigMC = 0, const AliEmcalJet *jetAssocMC = 0, const AliEmcalJet *jetTrigDet = 0, const AliEmcalJet *jetAssocDet = 0, Int_t type = 0);

 private:
  Bool_t            fDoMatchFullCharged;          //  do full-charged matching histos
  THnSparse        *fhnDiJetResponseCharged;      //! sparse with di-jet properties (full-full)
  THnSparse        *fhnDiJetResponseFullCharged;  //! sparse with di-jet properties (full-full)
  TH1F             *fh1TriggersCharged[2];        //! charged jet triggers
  TH1F             *fh1TriggersFull[2];           //! full jet triggers
  TH1F             *fh1TriggersLostCharged;       //! lost charged jet triggers
  TH1F             *fh1TriggersLostFull;          //! lost full jet triggers
  TH3F             *fh3AssocLostPtDeltaPhiCharged;//! lost charged associated jet
  TH3F             *fh3AssocLostPtDeltaPhiFull;   //! lost full associated jet
  THnSparse        *fhnMatchingCharged;           //! sparse comparing matched particle and detector level charged jets
  THnSparse        *fhnMatchingFull;              //! sparse comparing matched particle and detector level charged jets
  Int_t             fnUsedResponseVar;            //  build response for kt (0) or dijet eta (1)


  AliAnalysisTaskEmcalDiJetResponse(const AliAnalysisTaskEmcalDiJetResponse&);            // not implemented
  AliAnalysisTaskEmcalDiJetResponse &operator=(const AliAnalysisTaskEmcalDiJetResponse&); // not implemented

  ClassDef(AliAnalysisTaskEmcalDiJetResponse, 2) // jet sample analysis task
};
#endif
Commit	Line	Data
6e8b6371	1	#ifndef ALIANALYSISTASKEMCALDIJETRESPONSE_H
	2	#define ALIANALYSISTASKEMCALDIJETRESPONSE_H
	3
	4	class TH1;
	5	class TH2;
	6	class TH3;
	7	class TH3F;
	8	class THnSparse;
	9	class TClonesArray;
	10	class TArrayI;
	11	class AliAnalysisUtils;
	12	class AliAnalysisManager;
	13	class AliGenPythiaEventHeader;
	14
	15	#include "AliJetContainer.h"
	16
	17	#include "AliAnalysisTaskEmcalDiJetBase.h"
	18
	19	class AliAnalysisTaskEmcalDiJetResponse : public AliAnalysisTaskEmcalDiJetBase {
	20	public:
	21	AliAnalysisTaskEmcalDiJetResponse();
	22	AliAnalysisTaskEmcalDiJetResponse(const char *name);
	23	virtual ~AliAnalysisTaskEmcalDiJetResponse();
	24
	25	void UserCreateOutputObjects();
	26	void Terminate(Option_t *option);
	27
	28	//Setters
6ab30d5f	29	void SetMatchFullCharged(Bool_t b) { fDoMatchFullCharged = b; }
6ab30d5f	30	void SetResponseVar(Int_t v) { fnUsedResponseVar = v; }
6e8b6371	31
	32	//Getters
	33
	34	protected:
	35	Bool_t Run() ;
	36	void CorrelateJets(const Int_t type);
6ab30d5f	37	void CorrelateAllJets(const Int_t type);
	38	void CorrelateTwoJets(const Int_t type);
	39
6e8b6371	40	Bool_t FillHistograms() ;
	41	void FillDiJetHistos(const AliEmcalJet jet1 = 0, const AliEmcalJet jet2 = 0, const Int_t mode = 0);
	42	void FillMatchHistos();
	43	Bool_t RetrieveEventObjects();
	44
6ab30d5f	45	void FillDiJetResponse(const AliEmcalJet jetTrigMC = 0, const AliEmcalJet jetAssocMC = 0, const AliEmcalJet jetTrigDet = 0, const AliEmcalJet jetAssocDet = 0, Int_t type = 0);
6ab30d5f	46
6e8b6371	47	private:
	48	Bool_t fDoMatchFullCharged; // do full-charged matching histos
	49	THnSparse *fhnDiJetResponseCharged; //! sparse with di-jet properties (full-full)
	50	THnSparse *fhnDiJetResponseFullCharged; //! sparse with di-jet properties (full-full)
	51	TH1F *fh1TriggersCharged[2]; //! charged jet triggers
	52	TH1F *fh1TriggersFull[2]; //! full jet triggers
	53	TH1F *fh1TriggersLostCharged; //! lost charged jet triggers
	54	TH1F *fh1TriggersLostFull; //! lost full jet triggers
67c36e97	55	TH3F *fh3AssocLostPtDeltaPhiCharged;//! lost charged associated jet
67c36e97	56	TH3F *fh3AssocLostPtDeltaPhiFull; //! lost full associated jet
6e8b6371	57	THnSparse *fhnMatchingCharged; //! sparse comparing matched particle and detector level charged jets
6e8b6371	58	THnSparse *fhnMatchingFull; //! sparse comparing matched particle and detector level charged jets
6ab30d5f	59	Int_t fnUsedResponseVar; // build response for kt (0) or dijet eta (1)
6e8b6371	60
	61
	62	AliAnalysisTaskEmcalDiJetResponse(const AliAnalysisTaskEmcalDiJetResponse&); // not implemented
	63	AliAnalysisTaskEmcalDiJetResponse &operator=(const AliAnalysisTaskEmcalDiJetResponse&); // not implemented
	64
6ab30d5f	65	ClassDef(AliAnalysisTaskEmcalDiJetResponse, 2) // jet sample analysis task
6e8b6371	66	};
6e8b6371	67	#endif