#include "AliAnalysisTaskSE.h"
#include "TObjString.h"
#include "TString.h"
-#include "TObjArray.h"
////////////////
class TProfile;
class TRandom3;
class TRefArray;
-
+class TObjArray;
class AliAnalysisTaskJetBackgroundSubtract : public AliAnalysisTaskSE
{
virtual Bool_t Notify();
// Task specific methods...
- virtual void AddJetBranch(const char* c){fJBArray->Add(new TObjString(c));}
+ virtual void AddJetBranch(const char* c);
virtual void SetSubtractionMethod(Int_t i){fSubtraction = i;}
virtual Int_t GetSubtractionMethod(){return fSubtraction;}
+ virtual void SetKeepJets(Bool_t b = kTRUE){fKeepJets = b;}
virtual void SetBackgroundBranch(char* c){fBackgroundBranch = c;}
- virtual void SetNonStdFile(char* c){fNonStdFile = c;}
+ virtual void SetNonStdOutputFile(char* c){fNonStdFile = c;}
virtual void SetToReplace(char* c){fReplaceString1 = c;}
const char* GetToReplace(){return fReplaceString1.Data();}
virtual void SetReplacementMask(char* c){fReplaceString2 = c;}
AliAnalysisTaskJetBackgroundSubtract& operator=(const AliAnalysisTaskJetBackgroundSubtract&);
Bool_t RescaleJetMomentum(AliAODJet *jet,Float_t pT);
Bool_t RescaleJet4vector(AliAODJet *jet,TLorentzVector backgroundv);
-
+ Int_t MultFromJetRefs(TClonesArray* jets);
Double_t RecalcRho(TClonesArray* fbkgclusters,Double_t meanarea);
Double_t RhoRC(TClonesArray* fbkgclustersRC);
void ResetOutJets();
TString fReplaceString1; // To construct the new output name
TString fReplaceString2; // To construct the new output name
Int_t fSubtraction; // Parameter for subtraction mode
+ Bool_t fKeepJets; // keeps the jets with negative p_t rescaled to 0.1 GeV
TList *fInJetArrayList; //! transient list to make ease the handling of input jets
TList *fOutJetArrayList; //! transient list to make ease the reset of output jets
TH2F* fh2CentvsRho; //! centrality vs background density
TH2F* fh2CentvsSigma; //! centrality vs background sigma
+ TH2F* fh2MultvsRho; //! centrality vs background density
+ TH2F* fh2MultvsSigma; //! centrality vs background sigma
TH2F* fh2ShiftEta; //! extended correction Eta
TH2F* fh2ShiftPhi; //! extended correction Phi
TH2F* fh2ShiftEtaLeading; //! extended correction Eta leading jet
TList *fHistList; //! the histograms output list
- ClassDef(AliAnalysisTaskJetBackgroundSubtract, 4)
+ ClassDef(AliAnalysisTaskJetBackgroundSubtract, 7)
};
#endif