[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 AliJetContainer;

#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);

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