#include <TList.h>
#include "TObjString.h"
#include "AliVEvent.h"
+#include "AliLog.h"
+#include "AliAnalysisManager.h"
-#define VERBOSE_STAT
+//#define VERBOSE_STAT
class AliESDEvent;
class TH2F;
class TCollection;
class AliTriggerAnalysis;
class AliAnalysisTaskSE;
-
+class AliOADBPhysicsSelection ;
+class AliOADBFillingScheme ;
+class AliOADBTriggerAnalysis;
+class TPRegexp;
class AliPhysicsSelection : public AliAnalysisCuts
{
public:
public:
- enum {kStatTriggerClass=1,kStatHWTrig,kStatV0ABG,kStatV0CBG,kStatMB1,kStatMB1Prime,kStatFMD,kStatFO1,kStatFO2,kStatV0A,kStatV0C,kStatSSD1,kStatFO1AndV0,kStatV0,kStatAny2Hits,kStatOffline,kStatBG,kStatAccepted};
+ enum {kStatTriggerClass=1,kStatHWTrig,kStatV0ABG,kStatV0CBG,kStatMB1,kStatMB1Prime,kStatLaserCut,kStatFO1,kStatFO2,kStatFO2L1,kStatV0A,kStatV0C,kStatZDCA,kStatZDCC,kStatZDCAC,kStatZDCTime,kStatV0,kStatOffline,kStatBG,kStatAccepted, kStatAcceptedPileUp};
#ifdef VERBOSE_STAT
- enum {kStatRowBG=0,kStatRowAcc,kStatRowBGFrac,kStatRowAccFrac,kStatRowErrGoodFrac,kStatRowGoodFrac,kStatRowErrGood,kStatRowGood}; // offset wrt fBGStatOffset
+ enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG,kStatRowAcc,kStatRowBGFrac,kStatRowAccFrac,kStatRowErrGoodFrac,kStatRowGoodFrac,kStatRowErrGood,kStatRowGood}; // offset wrt fBGStatOffset
#else
- enum {kStatRowBG=0,kStatRowAcc,kStatRowGood}; // offset wrt fBGStatOffset
+ enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG, kStatRowAcc,kStatRowGood}; // offset wrt fBGStatOffset
#endif
enum {kStatIdxAll=0,kStatIdxBin0=1};
-
+ enum ETriggerLogic { kCINT1 = 0, kCMBS2A, kCMBS2C, kCMBAC, kCMBACS2, kHighMultL1 };
typedef Bool_t (*Bin0Callback_t)(const AliESDEvent *);
Int_t GetCurrentRun() const {return fCurrentRun;}
UInt_t IsCollisionCandidate(const AliESDEvent* aEsd);
+ Bool_t Initialize(const AliESDEvent* aEsd);
Bool_t Initialize(Int_t runNumber);
void SetAnalyzeMC(Bool_t flag = kTRUE) { fMC = flag; }
void SetSkipTriggerClassSelection(Bool_t flag = kTRUE) { fSkipTriggerClassSelection = flag; }
- void SetSkipV0(Bool_t flag=kTRUE) { fSkipV0 = flag;}
+ void SetSkipV0(Bool_t flag=kTRUE) { AliError ("This method is deprecated! Please change the OADB configuration instead") ; fSkipV0 = flag;}
void AddBackgroundIdentification(AliAnalysisCuts* background) { fBackgroundIdentification = background; }
- virtual void Print(Option_t* option = "") const;
+ virtual void Print(const Option_t* option = "") const;
virtual Long64_t Merge(TCollection* list);
- void SaveHistograms(const char* folder = 0) const;
+ void SaveHistograms(const char* folder = 0);
const TList* GetCollisionTriggerClasses() const { return &fCollTrigClasses; }
const TList* GetBGTriggerClasses() const { return &fBGTrigClasses; }
- void AddCollisionTriggerClass(const char* className){ fCollTrigClasses.Add(new TObjString(className)); fUsingCustomClasses = kTRUE; }
- void AddBGTriggerClass(const char* className) { fBGTrigClasses.Add(new TObjString(className)); fUsingCustomClasses = kTRUE; }
-
+ void AddCollisionTriggerClass(const char* className);
+ void AddBGTriggerClass(const char* className) ;
+ void SetCustomOADBObjects(AliOADBPhysicsSelection * oadbPS, AliOADBFillingScheme * oadbFS, AliOADBTriggerAnalysis * oadbTA = 0) { fPSOADB = oadbPS; fFillOADB = oadbFS; fTriggerOADB = oadbTA; fUsingCustomClasses = kTRUE;}
+ const AliOADBPhysicsSelection * GetOADBPhysicsSelection() const {return fPSOADB; }
+ const AliOADBFillingScheme * GetOADBFillingScheme() const {return fFillOADB;}
+ const AliOADBTriggerAnalysis * GetOADBTriggerAnalysis() const {return fTriggerOADB;}
+
+ Int_t GetBGStatOffset() const { return fBGStatOffset; }
+
+
AliTriggerAnalysis* GetTriggerAnalysis(Int_t i = 0) { return (fTriggerAnalysis.GetEntries() > 0) ? (AliTriggerAnalysis*) fTriggerAnalysis.At(i) : 0; }
const TH2F* GetStatisticsHistogram(Int_t idx=kStatIdxAll) const { return fHistStatistics[idx]; }
const TH2F* GetBunchCrossingHistogram() const { return fHistBunchCrossing; }
+ virtual TObject *GetStatistics(const Option_t *option) const;
void SetBIFactors(const AliESDEvent * aESD);
void SetBin0Callback( const char * cb) {fBin0CallBack = cb;}
void SetBin0CallbackViaPointer( Bin0Callback_t cb) {fBin0CallBackPointer = cb;}// WARNING: THIS SHOULD NOT BE USED, WILL BE REMOVED SOON
+ static const char * GetOADBFileName() { static TString filename; filename.Form("%s/COMMON/PHYSICSSELECTION/data/physicsSelection.root", AliAnalysisManager::GetOADBPath()); return filename.Data();};
+
protected:
- UInt_t CheckTriggerClass(const AliESDEvent* aEsd, const char* trigger) const;
- Int_t GetTriggerScheme(UInt_t runNumber) const;
- const char * GetBXIDs(UInt_t runNumber, const char * trigger ) ;
- const char * GetFillingScheme(UInt_t runNumber) ;
+ UInt_t CheckTriggerClass(const AliESDEvent* aEsd, const char* trigger, Int_t& triggerLogic) const;
+ Bool_t EvaluateTriggerLogic(const AliESDEvent* aEsd, AliTriggerAnalysis* triggerAnalysis, const char* triggerLogic, Bool_t offline);
TH2F * BookHistStatistics(const char * tag) ;
Int_t GetStatRow(const char * triggerBXClass, UInt_t offlineTriggerType, UInt_t ** rowIDs) const;
+ const char * GetTriggerString(TObjString * obj);
+
Int_t fCurrentRun; // run number for which the object is initialized
Bool_t fMC; // flag if MC is analyzed
AliAnalysisCuts* fBackgroundIdentification; // class that performs additional background identification
TH2F* fHistStatistics[2]; // how many events are cut away why {all,bin 0}
+ // TH2F* fHistStatisticsTokens; // how many events are cut away why (new version, bins for all tokens. Only used tokens are filled)
TH2F* fHistBunchCrossing; // histograms of accepted bunch crossing numbers
TH1F* fHistTriggerPattern; // Pattern of the individual detectors in the MB1 trigger. Can reveal inconsistencies/inefficiencies in the trigger
TString fBin0CallBack; // callback used to determine if an event is in the bin0 (name of the task where the callback is implemented);
Bin0Callback_t fBin0CallBackPointer; //! don't stream this. TO BE REMOVED SOON
- ClassDef(AliPhysicsSelection, 10)
+ Bool_t fIsPP; // True if processing pp run, false if heavy ion
+
+ AliOADBPhysicsSelection * fPSOADB; // Physics selection OADB object
+ AliOADBFillingScheme * fFillOADB; // Filling scheme OADB object
+ AliOADBTriggerAnalysis * fTriggerOADB; // Trigger analysis OADB object
+
+ TPRegexp* fRegexp; //! regular expression for trigger tokens
+ TList* fCashedTokens; //! trigger token lookup list
+
+ ClassDef(AliPhysicsSelection, 15)
private:
AliPhysicsSelection(const AliPhysicsSelection&);