// ALICE Off-Line framework:
//
// Piotr.Skowronski@cern.ch
-// more info: http://alisoft.cern.ch/people/skowron/analyzer/index.html
+// more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
//
////////////////////////////////////////////////////////////////////////////
//_________________________________________________________
#include <AliAnalysis.h>
-#include "AliAODPairCut.h"
-#include "AliAODParticleCut.h"
class AliHBTCut;
//class AliHBTPair;
class AliAOD;
class AliReader;
class AliEventBuffer;
+class AliEventCut;
class AliHBTOnePairFctn;
class AliHBTTwoPairFctn;
virtual ~AliHBTAnalysis();
Int_t Init();
- Int_t ProcessEvent(AliAOD* aodrec, AliAOD* aodsim = 0x0);
- Int_t Finish();
+ Int_t ProcessEvent(AliAOD* aodrec, AliAOD* aodsim = 0x0);
+ Int_t Finish();
- enum EProcessOption{kReconstructed,kSimulated,kSimulatedAndReconstructed};
- void SetProcessOption(EProcessOption option){fProcessOption = option;}//Init must be called to make effect
+ enum EProcessOption{kReconstructed,kSimulated,kSimulatedAndReconstructed};
+ void SetProcessOption(EProcessOption option){fProcessOption = option;}//Init must be called to make effect
- virtual void Process(Option_t* option = "TracksAndParticles");//Stand alone HBT analysis
+ void Process(Option_t* option = "TracksAndParticles");//Stand alone HBT analysis
- void SetGlobalPairCut(AliAODPairCut* cut);
+ void SetGlobalPairCut(AliAODPairCut* cut);
- void AddTrackFunction(AliHBTOnePairFctn* f);
- void AddParticleFunction(AliHBTOnePairFctn* f);
- void AddParticleAndTrackFunction(AliHBTTwoPairFctn* f);
+ void AddTrackFunction(AliHBTOnePairFctn* f);
+ void AddParticleFunction(AliHBTOnePairFctn* f);
+ void AddParticleAndTrackFunction(AliHBTTwoPairFctn* f);
- void AddParticleMonitorFunction(AliHBTMonOneParticleFctn* f); //z.ch.
- void AddTrackMonitorFunction(AliHBTMonOneParticleFctn* f); //z.ch.
- void AddParticleAndTrackMonitorFunction(AliHBTMonTwoParticleFctn* f);//z.ch.
+ void AddParticleMonitorFunction(AliHBTMonOneParticleFctn* f); //z.ch.
+ void AddTrackMonitorFunction(AliHBTMonOneParticleFctn* f); //z.ch.
+ void AddParticleAndTrackMonitorFunction(AliHBTMonTwoParticleFctn* f);//z.ch.
- void AddResolutionFunction(AliHBTTwoPairFctn* f){AddParticleAndTrackFunction(f);}
+ void AddResolutionFunction(AliHBTTwoPairFctn* f){AddParticleAndTrackFunction(f);}
- void SetReader(AliReader* r){fReader = r;}
+ void SetReader(AliReader* r){fReader = r;}
- void WriteFunctions();
-
- void SetBufferSize(Int_t buffsize){fBufferSize=buffsize;}
- void SetOwner(Bool_t owner=kTRUE){fIsOwner=owner;}
+ void WriteFunctions();
+ void SetOutputFileName(const char* fname);
+
+ void SetBufferSize(Int_t buffsize){fBufferSize=buffsize;}
+ void SetOwner(Bool_t owner=kTRUE){fIsOwner=owner;}
Bool_t IsOwner() const {return fIsOwner;}
Bool_t IsNonIdentAnalysis();
void ResetFunctions();
void SetDisplayInfo(Int_t howoften){fDisplayMixingInfo = howoften;}//defines every each line info about mixing is displayed
- void SetCutsOnParticles(); // -- aplies only to Process Tracks And Particles
- void SetCutsOnTracks();// -- aplies only to Process Tracks And Particles
- void SetCutsOnTracksAndParticles();// Default // -- aplies only to Process Tracks And Particles
+ void SetApparentVertex(Double_t x, Double_t y, Double_t z);//Sets apparent vertex
static void PressAnyKey();//small utility function that helps to make comfortable macros
- protected:
-
- Bool_t RunCoherencyCheck();
+
- void FilterOut(AliAOD* outpart1, AliAOD* outpart2, AliAOD* inpart,
- AliAOD* outtrack1, AliAOD* outtrack2, AliAOD* intrack)const;
- void FilterOut(AliAOD* out1, AliAOD* out2, AliAOD* in)const;
- void DeleteFunctions();
+ protected:
-
/**********************************************/
/* E V E N T P R O C E S S I N G */
/**********************************************/
// OLD legacy schema
- virtual void ProcessTracks();
- virtual void ProcessParticles();
- virtual void ProcessTracksAndParticles();
+ void ProcessTracks();
+ void ProcessParticles();
+ void ProcessTracksAndParticles();
- virtual void ProcessTracksAndParticlesNonIdentAnal();
- virtual void ProcessParticlesNonIdentAnal();
- virtual void ProcessTracksNonIdentAnal();
+ void ProcessTracksAndParticlesNonIdentAnal();
+ void ProcessParticlesNonIdentAnal();
+ void ProcessTracksNonIdentAnal();
+ Bool_t RunCoherencyCheck();
+
+ void FilterOut(AliAOD* outpart1, AliAOD* outpart2, AliAOD* inpart,
+ AliAOD* outtrack1, AliAOD* outtrack2, AliAOD* intrack)const;
+ void FilterOut(AliAOD* out1, AliAOD* out2, AliAOD* in)const;
+ void DeleteFunctions();
+
AliReader* fReader;//! Pointer to reader
AliHBTMonOneParticleFctn** fTrackMonitorFunctions; //! which are used for single particle analysis,
AliHBTMonTwoParticleFctn** fParticleAndTrackMonitorFunctions; //! cut monitoring, etc.
+ AliEventCut* fBkgEventCut;// We can narrow class of events used in
+
+ AliEventBuffer* fPartBuffer;//Sim Particles event buffer
+ AliEventBuffer* fTrackBuffer;//Rec Tracks event buffer
/**********************************************/
- /* Control parameters */
+ /* Control parameters */
/**********************************************/
-
- AliAODPairCut* fPairCut;//! Pair cut applied for all mixed particles
- Int_t fBufferSize; //!defines the size of buffer for mixed events; -1==MIX All
- Int_t fDisplayMixingInfo;//!defines every which particle mixing info is displayed
- Bool_t fIsOwner;//!defines of all functions are supposed to be deleted while by the way of analysis defaulr false
-
- AliEventBuffer* fPartBuffer;//Sim Particles event buffer
- AliEventBuffer* fTrackBuffer;//Rec Tracks event buffer
+ Int_t fBufferSize; //defines the size of buffer for mixed events; -1==MIX All
+ Int_t fDisplayMixingInfo;//!defines every which particle mixing info is displayed
+ Bool_t fIsOwner;//!defines of all functions are supposed to be deleted while by the way of analysis defaulr false
- EProcessOption fProcessOption;//Option that says waht analysis to do (Rec, Sim or SimAndRec)
- Bool_t fNoCorrfctns;//Internal flag indicating that no cfs are set by the user (only monitor ones)
+ EProcessOption fProcessOption;//Option that says waht analysis to do (Rec, Sim or SimAndRec)
+ Bool_t fNoCorrfctns;//Internal flag indicating that no cfs are set by the user (only monitor ones)
+ TString* fOutputFileName;//Fiele name where to dump results, if not specified reults are written to gDirectory
+
+ Double_t fVertexX;//X position of apparent vertex
+ Double_t fVertexY;//Y position of apparent vertex
+ Double_t fVertexZ;//Z position of apparent vertex
private:
-
- /**********************************************/
- /* C U T S */
- /**********************************************/
-
- Bool_t (AliHBTAnalysis::*fkPass)(AliAODPair* partpair, AliAODPair* trackpair) const;//Pointer to function that performes pair cut
- Bool_t (AliHBTAnalysis::*fkPass1)(AliVAODParticle* partpair, AliVAODParticle* trackpair) const;//Pointer to function that performes cut on first particle
- Bool_t (AliHBTAnalysis::*fkPass2)(AliVAODParticle* partpair, AliVAODParticle* trackpair) const;//Pointer to function that performes cut on second particle
- Bool_t (AliHBTAnalysis::*fkPassPairProp)(AliAODPair* partpair, AliAODPair* trackpair) const;//Pointer to function that performes pair cut
-
- Bool_t PassPartAndTrack (AliAODPair* partpair, AliAODPair* trackpair) const {return (fPairCut->Pass((AliAODPair*)partpair))?kTRUE:fPairCut->Pass((AliAODPair*)trackpair);}
- Bool_t PassPartAndTrack1(AliVAODParticle* part, AliVAODParticle* track) const;
- Bool_t PassPartAndTrack2(AliVAODParticle* part, AliVAODParticle* track) const;
- Bool_t PassPairPropPartAndTrack (AliAODPair* partpair, AliAODPair* trackpair) const {return (fPairCut->PassPairProp((AliAODPair*)partpair))?kTRUE:fPairCut->PassPairProp((AliAODPair*)trackpair);}
-
- Bool_t PassPart (AliAODPair* partpair, AliAODPair* /*trackpair*/) const {return fPairCut->Pass((AliAODPair*)partpair);}
- Bool_t PassPart1(AliVAODParticle* part, AliVAODParticle* /*track*/) const {return fPairCut->GetFirstPartCut()->Pass(part);}
- Bool_t PassPart2(AliVAODParticle* part, AliVAODParticle* /*track*/) const {return fPairCut->GetSecondPartCut()->Pass(part);}
- Bool_t PassPairPropPart (AliAODPair* partpair, AliAODPair* /*trackpair*/) const {return fPairCut->PassPairProp((AliAODPair*)partpair);}
-
- Bool_t PassTrack (AliAODPair* /*partpair*/, AliAODPair* trackpair) const {return fPairCut->Pass((AliAODPair*)trackpair);}
- Bool_t PassTrack1(AliVAODParticle* /*part*/, AliVAODParticle* track) const {return fPairCut->GetFirstPartCut()->Pass(track);}
- Bool_t PassTrack2(AliVAODParticle* /*part*/, AliVAODParticle* track) const {return fPairCut->GetSecondPartCut()->Pass(track);}
- Bool_t PassPairPropTrack (AliAODPair* /*partpair*/, AliAODPair* trackpair) const {return fPairCut->PassPairProp((AliAODPair*)trackpair);}
static const UInt_t fgkFctnArraySize;//!
static const UInt_t fgkDefaultMixingInfo;//!
ClassDef(AliHBTAnalysis,0)
};
-
-inline Bool_t AliHBTAnalysis::PassPartAndTrack1(AliVAODParticle* part,AliVAODParticle* track) const
-{
-//Checks first particle from both, particle and track pairs
- AliAODParticleCut* pc = fPairCut->GetFirstPartCut();
- return (pc->Pass(part))?kTRUE:pc->Pass(track);
-}
-/*************************************************************************************/
-inline Bool_t AliHBTAnalysis::PassPartAndTrack2(AliVAODParticle* part,AliVAODParticle* track) const
-{
-//Checks second particle from both, particle and track pairs
- AliAODParticleCut* pc = fPairCut->GetSecondPartCut();
- return (pc->Pass(part))?kTRUE:pc->Pass(track);
-}
-/*************************************************************************************/
-
#endif