[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / AliEmcalJetTask.h

#ifndef ALIEMCALJETTASK_H
#define ALIEMCALJETTASK_H

// $Id$

class TClonesArray;
class AliVEvent;

namespace fastjet {
  class PseudoJet;
}

#include "AliLog.h"
#include "AliAnalysisTaskSE.h"

class AliEmcalJetTask : public AliAnalysisTaskSE {
 public:

  enum JetType {
    kNone=0,
    kKT=1<<0,
    kAKT=1<<1,
    kFullJet=1<<2,
    kChargedJet=1<<3,
    kNeutralJet=1<<4,
    kR020Jet=1<<5,
    kR030Jet=1<<6,
    kR040Jet=1<<7,
    kRX1Jet=1<<8,  // user defined radii, use SetRadius(Double_t)
    kRX2Jet=1<<9,
    kRX3Jet=1<<10
  };


  AliEmcalJetTask();
  AliEmcalJetTask(const char *name);
  virtual ~AliEmcalJetTask();

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

  void                   SelectConstituents(UInt_t constSel, UInt_t MCconstSel)  { fConstSel = constSel; fMCConstSel = MCconstSel; };
  void                   SetAlgo(Int_t a)                 { if (a==0) fJetType |= kKT; else fJetType |= kAKT; }  // for backward compatibility only
  void                   SetClusName(const char *n)       { fCaloName      = n     ; }
  void                   SetJetsName(const char *n)       { fJetsName      = n     ; }
  void                   SetJetType(UInt_t t)             { fJetType       = t     ; }
  void                   SetMarkConstituents(UInt_t m)    { fMarkConst     = m     ; }
  void                   SetMinJetArea(Double_t a)        { fMinJetArea    = a     ; }
  void                   SetMinJetClusPt(Double_t min)    { fMinJetClusPt  = min   ; }
  void                   SetMinJetPt(Double_t j)          { fMinJetPt      = j     ; }
  void                   SetMinJetTrackPt(Double_t min)   { fMinJetTrackPt = min   ; }
  void                   SetRadius(Double_t r)            { fRadius        = r     ; if ((fJetType & (kRX1Jet|kRX2Jet|kRX3Jet)) == 0) AliWarning("Radius value will be ignored if jet type is not set to a user defined radius (kRX1Jet,kRX2Jet,kRX3Jet)."); }
  void                   SetTrackEfficiency(Double_t t)   { fTrackEfficiency = t   ; }
  void                   SetTracksName(const char *n)     { fTracksName    = n     ; }
  void                   SetType(Int_t t)                 { if (t==0) fJetType |= kFullJet; 
                                                            else if (t==1) fJetType |= kChargedJet; 
                                                            else if (t==2) fJetType |= kNeutralJet; } // for backward compatibility only
  void                   SetEtaRange(Double_t emi, Double_t ema) {fEtaMin = emi; fEtaMax = ema; }
  void                   SetPhiRange(Double_t pmi, Double_t pma) {fPhiMin = pmi; fPhiMax = pma; }
  void                   SetJetEtaRange(Double_t emi, Double_t ema) {fJetEtaMin = emi; fJetEtaMax = ema; }
  void                   SetJetPhiRange(Double_t pmi, Double_t pma) {fJetPhiMin = pmi; fJetPhiMax = pma; }
  void                   SetGhostArea(Double_t gharea)    { fGhostArea      = gharea;  }
  void                   SetMinMCLabel(Int_t s)           { fMinMCLabel     = s     ;  }
  void                   SetRecombScheme(Int_t scheme)    { fRecombScheme   = scheme;  }
  void                   SelectCollisionCandidates(UInt_t offlineTriggerMask = AliVEvent::kMB)
  {
    if(!fIsPSelSet)
    {
      fIsPSelSet = kTRUE;
      fOfflineTriggerMask = offlineTriggerMask;
    }
    else
    {
      AliWarning("Manually setting the event selection for jet finders is not allowed! Using trigger=old_trigger | your_trigger");  
      fOfflineTriggerMask = fOfflineTriggerMask | offlineTriggerMask;
    }
  }
  void                   SetLegacyMode(Bool_t mode)       { fLegacyMode ^= mode; }
  void                   SetCodeDebug(Bool_t val)         { fCodeDebug = val; }

  UInt_t                 GetJetType()                     { return fJetType; }
  Bool_t                 GetLegacyMode()                  { return fLegacyMode; }

 protected:
  void                   FindJets();
  Bool_t                 DoInit();
  Bool_t                 GetSortedArray(Int_t indexes[], std::vector<fastjet::PseudoJet> array) const;

  TString                fTracksName;             // name of track collection
  TString                fCaloName;               // name of calo cluster collection
  TString                fJetsName;               // name of jet collection
  UInt_t                 fJetType;                // jet type (algorithm, radius, constituents)
  UInt_t                 fConstSel;               // select constituents from a previous jet finding
  UInt_t                 fMCConstSel;             // select MC constituents (label!=0) from a previous jet finding
  UInt_t                 fMarkConst;              // constituents are marked (via TObject::SetBit) as belonging to the # leading jets
  Double_t               fRadius;                 // jet radius
  Double_t               fMinJetTrackPt;          // min jet track momentum   (applied before clustering)
  Double_t               fMinJetClusPt;           // min jet cluster momentum (applied before clustering)
  Double_t               fPhiMin;                 // minimum phi for constituents (applied before clustering)
  Double_t               fPhiMax;                 // maximum phi for constituents (applied before clustering)
  Double_t               fEtaMin;                 // minimum eta for constituents (applied before clustering)
  Double_t               fEtaMax;                 // maximum eta for constituents (applied before clustering)
  Double_t               fMinJetArea;             // min area to keep jet in output
  Double_t               fMinJetPt;               // min jet pt to keep jet in output
  Double_t               fJetPhiMin;              // minimum phi to keep jet in output
  Double_t               fJetPhiMax;              // maximum phi to keep jet in output
  Double_t               fJetEtaMin;              // minimum eta to keep jet in output
  Double_t               fJetEtaMax;              // maximum eta to keep jet in output
  Double_t               fGhostArea;              // ghost area
  Int_t                  fMinMCLabel;             // minimum MC label value for the tracks/clusters being considered MC particles
  Int_t                  fRecombScheme;           // recombination scheme used by fastjet
  Double_t               fTrackEfficiency;        // artificial tracking inefficiency (0...1)
  Bool_t                 fIsInit;                 //!=true if already initialized
  Bool_t                 fIsPSelSet;              //!=true if physics selection was set
  Bool_t                 fIsMcPart;               //!=true if MC particles are given as input
  Bool_t                 fIsEmcPart;              //!=true if emcal particles are given as input (for clusters)
  Bool_t                 fLegacyMode;             //! if true, enable FJ 2.x behavior
  Bool_t                 fCodeDebug;              // use nontested code changes 
  TClonesArray          *fJets;                   //!jet collection
  AliVEvent             *fEvent;                  //!current event
  TClonesArray          *fTracks;                 //!tracks collection
  TClonesArray          *fClus;                   //!cluster collection

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

  ClassDef(AliEmcalJetTask, 10) // Jet producing task
};
#endif
Commit	Line	Data
914d486c	1	#ifndef ALIEMCALJETTASK_H
914d486c	2	#define ALIEMCALJETTASK_H
48d874ff	3
	4	// $Id$
	5
	6	class TClonesArray;
b65fac7a	7	class AliVEvent;
48d874ff	8
7030f36f	9	namespace fastjet {
	10	class PseudoJet;
	11	}
	12
d14d6e2c	13	#include "AliLog.h"
48d874ff	14	#include "AliAnalysisTaskSE.h"
48d874ff	15
914d486c	16	class AliEmcalJetTask : public AliAnalysisTaskSE {
48d874ff	17	public:
6a20534a	18
	19	enum JetType {
	20	kNone=0,
	21	kKT=1<<0,
	22	kAKT=1<<1,
	23	kFullJet=1<<2,
	24	kChargedJet=1<<3,
	25	kNeutralJet=1<<4,
	26	kR020Jet=1<<5,
	27	kR030Jet=1<<6,
	28	kR040Jet=1<<7,
	29	kRX1Jet=1<<8, // user defined radii, use SetRadius(Double_t)
	30	kRX2Jet=1<<9,
d1f2e0d9	31	kRX3Jet=1<<10
6a20534a	32	};
	33
	34
914d486c	35	AliEmcalJetTask();
	36	AliEmcalJetTask(const char *name);
	37	virtual ~AliEmcalJetTask();
48d874ff	38
b65fac7a	39	void UserCreateOutputObjects();
	40	void UserExec(Option_t *option);
	41	void Terminate(Option_t *option);
	42
6a20534a	43	void SelectConstituents(UInt_t constSel, UInt_t MCconstSel) { fConstSel = constSel; fMCConstSel = MCconstSel; };
6a20534a	44	void SetAlgo(Int_t a) { if (a==0) fJetType \|= kKT; else fJetType \|= kAKT; } // for backward compatibility only
b65fac7a	45	void SetClusName(const char *n) { fCaloName = n ; }
b65fac7a	46	void SetJetsName(const char *n) { fJetsName = n ; }
6a20534a	47	void SetJetType(UInt_t t) { fJetType = t ; }
7030f36f	48	void SetMarkConstituents(UInt_t m) { fMarkConst = m ; }
b65fac7a	49	void SetMinJetArea(Double_t a) { fMinJetArea = a ; }
111318e8	50	void SetMinJetClusPt(Double_t min) { fMinJetClusPt = min ; }
b65fac7a	51	void SetMinJetPt(Double_t j) { fMinJetPt = j ; }
111318e8	52	void SetMinJetTrackPt(Double_t min) { fMinJetTrackPt = min ; }
6a20534a	53	void SetRadius(Double_t r) { fRadius = r ; if ((fJetType & (kRX1Jet\|kRX2Jet\|kRX3Jet)) == 0) AliWarning("Radius value will be ignored if jet type is not set to a user defined radius (kRX1Jet,kRX2Jet,kRX3Jet)."); }
6421eeb0	54	void SetTrackEfficiency(Double_t t) { fTrackEfficiency = t ; }
111318e8	55	void SetTracksName(const char *n) { fTracksName = n ; }
6a20534a	56	void SetType(Int_t t) { if (t==0) fJetType \|= kFullJet;
	57	else if (t==1) fJetType \|= kChargedJet;
	58	else if (t==2) fJetType \|= kNeutralJet; } // for backward compatibility only
04905231	59	void SetEtaRange(Double_t emi, Double_t ema) {fEtaMin = emi; fEtaMax = ema; }
04905231	60	void SetPhiRange(Double_t pmi, Double_t pma) {fPhiMin = pmi; fPhiMax = pma; }
c9ad772c	61	void SetJetEtaRange(Double_t emi, Double_t ema) {fJetEtaMin = emi; fJetEtaMax = ema; }
c9ad772c	62	void SetJetPhiRange(Double_t pmi, Double_t pma) {fJetPhiMin = pmi; fJetPhiMax = pma; }
1f9c287f	63	void SetGhostArea(Double_t gharea) { fGhostArea = gharea; }
a7477843	64	void SetMinMCLabel(Int_t s) { fMinMCLabel = s ; }
8612dfc8	65	void SetRecombScheme(Int_t scheme) { fRecombScheme = scheme; }
6421eeb0	66	void SelectCollisionCandidates(UInt_t offlineTriggerMask = AliVEvent::kMB)
d14d6e2c	67	{
	68	if(!fIsPSelSet)
	69	{
	70	fIsPSelSet = kTRUE;
	71	fOfflineTriggerMask = offlineTriggerMask;
	72	}
	73	else
	74	{
	75	AliWarning("Manually setting the event selection for jet finders is not allowed! Using trigger=old_trigger \| your_trigger");
	76	fOfflineTriggerMask = fOfflineTriggerMask \| offlineTriggerMask;
	77	}
	78	}
5d0a5007	79	void SetLegacyMode(Bool_t mode) { fLegacyMode ^= mode; }
ffdca6aa	80	void SetCodeDebug(Bool_t val) { fCodeDebug = val; }
6a20534a	81
6a20534a	82	UInt_t GetJetType() { return fJetType; }
5d0a5007	83	Bool_t GetLegacyMode() { return fLegacyMode; }
6a20534a	84
48d874ff	85	protected:
04905231	86	void FindJets();
04905231	87	Bool_t DoInit();
ca5c29fa	88	Bool_t GetSortedArray(Int_t indexes[], std::vector<fastjet::PseudoJet> array) const;
48d874ff	89
900f5135	90	TString fTracksName; // name of track collection
900f5135	91	TString fCaloName; // name of calo cluster collection
7efbea04	92	TString fJetsName; // name of jet collection
6a20534a	93	UInt_t fJetType; // jet type (algorithm, radius, constituents)
	94	UInt_t fConstSel; // select constituents from a previous jet finding
	95	UInt_t fMCConstSel; // select MC constituents (label!=0) from a previous jet finding
7030f36f	96	UInt_t fMarkConst; // constituents are marked (via TObject::SetBit) as belonging to the # leading jets
48d874ff	97	Double_t fRadius; // jet radius
f5f3c8e9	98	Double_t fMinJetTrackPt; // min jet track momentum (applied before clustering)
f5f3c8e9	99	Double_t fMinJetClusPt; // min jet cluster momentum (applied before clustering)
04905231	100	Double_t fPhiMin; // minimum phi for constituents (applied before clustering)
	101	Double_t fPhiMax; // maximum phi for constituents (applied before clustering)
	102	Double_t fEtaMin; // minimum eta for constituents (applied before clustering)
	103	Double_t fEtaMax; // maximum eta for constituents (applied before clustering)
f5f3c8e9	104	Double_t fMinJetArea; // min area to keep jet in output
f5f3c8e9	105	Double_t fMinJetPt; // min jet pt to keep jet in output
c9ad772c	106	Double_t fJetPhiMin; // minimum phi to keep jet in output
	107	Double_t fJetPhiMax; // maximum phi to keep jet in output
	108	Double_t fJetEtaMin; // minimum eta to keep jet in output
	109	Double_t fJetEtaMax; // maximum eta to keep jet in output
1f9c287f	110	Double_t fGhostArea; // ghost area
a7477843	111	Int_t fMinMCLabel; // minimum MC label value for the tracks/clusters being considered MC particles
7071e730	112	Int_t fRecombScheme; // recombination scheme used by fastjet
6421eeb0	113	Double_t fTrackEfficiency; // artificial tracking inefficiency (0...1)
04905231	114	Bool_t fIsInit; //!=true if already initialized
d14d6e2c	115	Bool_t fIsPSelSet; //!=true if physics selection was set
04905231	116	Bool_t fIsMcPart; //!=true if MC particles are given as input
04905231	117	Bool_t fIsEmcPart; //!=true if emcal particles are given as input (for clusters)
5d0a5007	118	Bool_t fLegacyMode; //! if true, enable FJ 2.x behavior
ffdca6aa	119	Bool_t fCodeDebug; // use nontested code changes
48d874ff	120	TClonesArray *fJets; //!jet collection
b65fac7a	121	AliVEvent *fEvent; //!current event
04905231	122	TClonesArray *fTracks; //!tracks collection
04905231	123	TClonesArray *fClus; //!cluster collection
48d874ff	124
48d874ff	125	private:
914d486c	126	AliEmcalJetTask(const AliEmcalJetTask&); // not implemented
914d486c	127	AliEmcalJetTask &operator=(const AliEmcalJetTask&); // not implemented
48d874ff	128
6421eeb0	129	ClassDef(AliEmcalJetTask, 10) // Jet producing task
48d874ff	130	};
48d874ff	131	#endif