[u/mrichter/AliRoot.git] / PWGGA / EMCALJetTasks / AliJetResponseMaker.h

#ifndef ALIJETRESPONSEMAKER_H
#define ALIJETRESPONSEMAKER_H

// $Id$

class TClonesArray;
class TH1F;
class TH2F;

#include "AliAnalysisTaskEmcalJet.h"

class AliJetResponseMaker : public AliAnalysisTaskEmcalJet {
 public:
  AliJetResponseMaker();
  AliJetResponseMaker(const char *name);
  virtual ~AliJetResponseMaker();

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

  void                        SetMCJetsName(const char *n)       { fMCJetsName    = n; }
  void                        SetMCTracksName(const char *n)     { fMCTracksName  = n; }
  void                        SetMaxDistance(Double_t d)         { fMaxDistance   = d; }

 protected:
  void                        DoJetLoop(TClonesArray *jets1, TClonesArray *jets2, Bool_t mc);
  Bool_t                      FillHistograms();
  Bool_t                      RetrieveEventObjects();

  TString                     fMCTracksName;              // name of MC particle collection
  TString                     fMCJetsName;                // name of MC jet collection
  Double_t                    fMaxDistance;               // maximum distance between matched particle and detector level jets
  TClonesArray               *fMCTracks;                  //!MC particles
  TClonesArray               *fMCJets;                    //!MC jets
  // Particle level jets
  TH2F                       *fHistMCJetPhiEta;           //!phi-eta distribution of jets
  TH1F                       *fHistMCJetsPt;              //!inclusive jet pt spectrum
  TH2F                       *fHistMCJetsNEFvsPt;         //!jet neutral energy fraction vs. jet pt
  TH2F                       *fHistMCJetsZvsPt;           //!constituent Pt over Jet Pt ratio vs. jet pt
  // Detector level jets
  TH2F                       *fHistJetPhiEta;             //!phi-eta distribution of jets
  TH1F                       *fHistJetsPt;                //!inclusive jet pt spectrum
  TH2F                       *fHistJetsNEFvsPt;           //!jet neutral energy fraction vs. jet pt
  TH2F                       *fHistJetsZvsPt;             //!constituent Pt over Jet Pt ratio vs. jet pt
  // Detector-particle level matching
  TH1F                       *fHistClosestDistance;       //!distance between closest particle to detector level jet
  TH1F                       *fHistClosestDeltaPhi;       //!delta phi between closest particle to detector level jet
  TH1F                       *fHistClosestDeltaEta;       //!delta eta between closest particle to detector level jet
  TH1F                       *fHistClosestDeltaPt;        //!delta pt between closest particle to detector level jet
  TH1F                       *fHistNonMatchedMCJetPt;     //!non-matched mc jet pt distribution
  TH1F                       *fHistNonMatchedJetPt;       //!non-matched jet pt distribution
  TH2F                       *fHistPartvsDetecPt;         //!particle vs detector level jet pt
  TH1F                       *fHistMissedMCJets;          //!mc jets not measured

 private:
  AliJetResponseMaker(const AliJetResponseMaker&);            // not implemented
  AliJetResponseMaker &operator=(const AliJetResponseMaker&); // not implemented

  ClassDef(AliJetResponseMaker, 4) // Jet response matrix producing task
};
#endif
Commit	Line	Data
2949a2ec	1	#ifndef ALIJETRESPONSEMAKER_H
	2	#define ALIJETRESPONSEMAKER_H
	3
7549c451	4	// $Id$
2949a2ec	5
	6	class TClonesArray;
	7	class TH1F;
	8	class TH2F;
	9
6fd5039f	10	#include "AliAnalysisTaskEmcalJet.h"
2949a2ec	11
6fd5039f	12	class AliJetResponseMaker : public AliAnalysisTaskEmcalJet {
2949a2ec	13	public:
	14	AliJetResponseMaker();
	15	AliJetResponseMaker(const char *name);
	16	virtual ~AliJetResponseMaker();
	17
a7ab01d5	18	void UserCreateOutputObjects();
	19	void Terminate(Option_t *option);
	20
1b76c28f	21	void SetMCJetsName(const char *n) { fMCJetsName = n; }
	22	void SetMCTracksName(const char *n) { fMCTracksName = n; }
	23	void SetMaxDistance(Double_t d) { fMaxDistance = d; }
2949a2ec	24
2949a2ec	25	protected:
99c67c1b	26	void DoJetLoop(TClonesArray jets1, TClonesArray jets2, Bool_t mc);
6fd5039f	27	Bool_t FillHistograms();
6fd5039f	28	Bool_t RetrieveEventObjects();
2949a2ec	29
a7ab01d5	30	TString fMCTracksName; // name of MC particle collection
2949a2ec	31	TString fMCJetsName; // name of MC jet collection
1b76c28f	32	Double_t fMaxDistance; // maximum distance between matched particle and detector level jets
a7ab01d5	33	TClonesArray *fMCTracks; //!MC particles
2949a2ec	34	TClonesArray *fMCJets; //!MC jets
	35	// Particle level jets
	36	TH2F *fHistMCJetPhiEta; //!phi-eta distribution of jets
	37	TH1F *fHistMCJetsPt; //!inclusive jet pt spectrum
2949a2ec	38	TH2F *fHistMCJetsNEFvsPt; //!jet neutral energy fraction vs. jet pt
	39	TH2F *fHistMCJetsZvsPt; //!constituent Pt over Jet Pt ratio vs. jet pt
	40	// Detector level jets
	41	TH2F *fHistJetPhiEta; //!phi-eta distribution of jets
	42	TH1F *fHistJetsPt; //!inclusive jet pt spectrum
2949a2ec	43	TH2F *fHistJetsNEFvsPt; //!jet neutral energy fraction vs. jet pt
2949a2ec	44	TH2F *fHistJetsZvsPt; //!constituent Pt over Jet Pt ratio vs. jet pt
1b76c28f	45	// Detector-particle level matching
	46	TH1F *fHistClosestDistance; //!distance between closest particle to detector level jet
	47	TH1F *fHistClosestDeltaPhi; //!delta phi between closest particle to detector level jet
	48	TH1F *fHistClosestDeltaEta; //!delta eta between closest particle to detector level jet
	49	TH1F *fHistClosestDeltaPt; //!delta pt between closest particle to detector level jet
99c67c1b	50	TH1F *fHistNonMatchedMCJetPt; //!non-matched mc jet pt distribution
99c67c1b	51	TH1F *fHistNonMatchedJetPt; //!non-matched jet pt distribution
1b76c28f	52	TH2F *fHistPartvsDetecPt; //!particle vs detector level jet pt
2bddb6ae	53	TH1F *fHistMissedMCJets; //!mc jets not measured
2949a2ec	54
	55	private:
	56	AliJetResponseMaker(const AliJetResponseMaker&); // not implemented
	57	AliJetResponseMaker &operator=(const AliJetResponseMaker&); // not implemented
	58
2bddb6ae	59	ClassDef(AliJetResponseMaker, 4) // Jet response matrix producing task
2949a2ec	60	};
2949a2ec	61	#endif