1 #ifndef ALIHBTANALYSIS_H
2 #define ALIHBTANALYSIS_H
3 //_________________________________________________________
4 ////////////////////////////////////////////////////////////////////////////
6 // class AliHBTAnalysis
8 // Central Object Of HBTAnalyser:
9 // This class performs main looping within HBT Analysis
10 // User must plug a reader of Type AliReader
11 // User plugs in coorelation and monitor functions
12 // as well as monitor functions
14 // HBT Analysis Tool, which is integral part of AliRoot,
15 // ALICE Off-Line framework:
17 // Piotr.Skowronski@cern.ch
18 // more info: http://alisoft.cern.ch/people/skowron/analyzer/index.html
20 ////////////////////////////////////////////////////////////////////////////
21 //_________________________________________________________
23 #include <AliAnalysis.h>
24 #include "AliAODPairCut.h"
25 #include "AliAODParticleCut.h"
33 class AliHBTOnePairFctn;
34 class AliHBTTwoPairFctn;
36 class AliHBTMonOneParticleFctn;
37 class AliHBTMonTwoParticleFctn;
41 class AliHBTAnalysis: public AliAnalysis
45 AliHBTAnalysis(const AliHBTAnalysis& in);
46 AliHBTAnalysis& operator=(const AliHBTAnalysis& /*right*/);
47 virtual ~AliHBTAnalysis();
50 Int_t ProcessEvent(AliAOD* aodrec, AliAOD* aodsim = 0x0);
53 virtual void Process(Option_t* option = "TracksAndParticles");
55 void SetGlobalPairCut(AliAODPairCut* cut);
57 void AddTrackFunction(AliHBTOnePairFctn* f);
58 void AddParticleFunction(AliHBTOnePairFctn* f);
59 void AddParticleAndTrackFunction(AliHBTTwoPairFctn* f);
61 void AddParticleMonitorFunction(AliHBTMonOneParticleFctn* f); //z.ch.
62 void AddTrackMonitorFunction(AliHBTMonOneParticleFctn* f); //z.ch.
63 void AddParticleAndTrackMonitorFunction(AliHBTMonTwoParticleFctn* f);//z.ch.
65 void AddResolutionFunction(AliHBTTwoPairFctn* f){AddParticleAndTrackFunction(f);}
67 void SetReader(AliReader* r){fReader = r;}
69 void WriteFunctions();
71 void SetBufferSize(Int_t buffsize){fBufferSize=buffsize;}
72 void SetOwner(Bool_t owner=kTRUE){fIsOwner=owner;}
73 Bool_t IsOwner() const {return fIsOwner;}
74 Bool_t IsNonIdentAnalysis();
75 void ResetFunctions();
76 void SetDisplayInfo(Int_t howoften){fDisplayMixingInfo = howoften;}//defines every each line info about mixing is displayed
78 void SetCutsOnParticles(); // -- aplies only to Process Tracks And Particles
79 void SetCutsOnTracks();// -- aplies only to Process Tracks And Particles
80 void SetCutsOnTracksAndParticles();// Default // -- aplies only to Process Tracks And Particles
82 static void PressAnyKey();//small utility function that helps to make comfortable macros
85 Bool_t RunCoherencyCheck();
87 void FilterOut(AliAOD* outpart1, AliAOD* outpart2, AliAOD* inpart,
88 AliAOD* outtrack1, AliAOD* outtrack2, AliAOD* intrack)const;
89 void FilterOut(AliAOD* out1, AliAOD* out2, AliAOD* in)const;
90 void DeleteFunctions();
92 virtual void ProcessTracks();
93 virtual void ProcessParticles();
94 virtual void ProcessTracksAndParticles();
96 virtual void ProcessTracksAndParticlesNonIdentAnal();
97 virtual void ProcessParticlesNonIdentAnal();
98 virtual void ProcessTracksNonIdentAnal();
100 AliReader* fReader;//! Pointer to reader
102 UInt_t fNTrackFunctions; //! Number of Tracks functions
103 UInt_t fNParticleFunctions; //! Number of particles functions
104 UInt_t fNParticleAndTrackFunctions; //! Number of resolution functions
106 UInt_t fNTrackMonitorFunctions; //! Number of Track Monitor functions
107 UInt_t fNParticleMonitorFunctions; //! Number of Particles Monitor functions
108 UInt_t fNParticleAndTrackMonitorFunctions; //! Number of Resolution Monitor functions
110 AliHBTOnePairFctn** fTrackFunctions; //!array of pointers to functions that analyze rekonstructed tracks
111 AliHBTOnePairFctn** fParticleFunctions; //!array of pointers to functions that analyze generated particles
112 AliHBTTwoPairFctn** fParticleAndTrackFunctions; //!array of pointers to functions that analyze both
113 //reconstructed tracks and generated particles
114 //i.e. - resolution analyzers
115 AliHBTMonOneParticleFctn** fParticleMonitorFunctions; //! array of pointers to monitoring functions
116 AliHBTMonOneParticleFctn** fTrackMonitorFunctions; //! which are used for single particle analysis,
117 AliHBTMonTwoParticleFctn** fParticleAndTrackMonitorFunctions; //! cut monitoring, etc.
120 /**********************************************/
121 /* Control parameters */
122 /**********************************************/
124 AliAODPairCut* fPairCut;//! Pair cut applied for all mixed particles
126 Int_t fBufferSize; //!defines the size of buffer for mixed events; -1==MIX All
127 Int_t fDisplayMixingInfo;//!defines every which particle mixing info is displayed
128 Bool_t fIsOwner;//!defines of all functions are supposed to be deleted while by the way of analysis defaulr false
131 Bool_t (AliHBTAnalysis::*fkPass)(AliAODPair* partpair, AliAODPair* trackpair) const;//Pointer to function that performes pair cut
132 Bool_t (AliHBTAnalysis::*fkPass1)(AliVAODParticle* partpair, AliVAODParticle* trackpair) const;//Pointer to function that performes cut on first particle
133 Bool_t (AliHBTAnalysis::*fkPass2)(AliVAODParticle* partpair, AliVAODParticle* trackpair) const;//Pointer to function that performes cut on second particle
134 Bool_t (AliHBTAnalysis::*fkPassPairProp)(AliAODPair* partpair, AliAODPair* trackpair) const;//Pointer to function that performes pair cut
136 Bool_t PassPartAndTrack (AliAODPair* partpair, AliAODPair* trackpair) const {return (fPairCut->Pass((AliAODPair*)partpair))?kTRUE:fPairCut->Pass((AliAODPair*)trackpair);}
137 Bool_t PassPartAndTrack1(AliVAODParticle* part, AliVAODParticle* track) const;
138 Bool_t PassPartAndTrack2(AliVAODParticle* part, AliVAODParticle* track) const;
139 Bool_t PassPairPropPartAndTrack (AliAODPair* partpair, AliAODPair* trackpair) const {return (fPairCut->PassPairProp((AliAODPair*)partpair))?kTRUE:fPairCut->PassPairProp((AliAODPair*)trackpair);}
141 Bool_t PassPart (AliAODPair* partpair, AliAODPair* /*trackpair*/) const {return fPairCut->Pass((AliAODPair*)partpair);}
142 Bool_t PassPart1(AliVAODParticle* part, AliVAODParticle* /*track*/) const {return fPairCut->GetFirstPartCut()->Pass(part);}
143 Bool_t PassPart2(AliVAODParticle* part, AliVAODParticle* /*track*/) const {return fPairCut->GetSecondPartCut()->Pass(part);}
144 Bool_t PassPairPropPart (AliAODPair* partpair, AliAODPair* /*trackpair*/) const {return fPairCut->PassPairProp((AliAODPair*)partpair);}
146 Bool_t PassTrack (AliAODPair* /*partpair*/, AliAODPair* trackpair) const {return fPairCut->Pass((AliAODPair*)trackpair);}
147 Bool_t PassTrack1(AliVAODParticle* /*part*/, AliVAODParticle* track) const {return fPairCut->GetFirstPartCut()->Pass(track);}
148 Bool_t PassTrack2(AliVAODParticle* /*part*/, AliVAODParticle* track) const {return fPairCut->GetSecondPartCut()->Pass(track);}
149 Bool_t PassPairPropTrack (AliAODPair* /*partpair*/, AliAODPair* trackpair) const {return fPairCut->PassPairProp((AliAODPair*)trackpair);}
151 static const UInt_t fgkFctnArraySize;//!
152 static const UInt_t fgkDefaultMixingInfo;//!
153 static const Int_t fgkDefaultBufferSize;//!
155 ClassDef(AliHBTAnalysis,0)
158 inline Bool_t AliHBTAnalysis::PassPartAndTrack1(AliVAODParticle* part,AliVAODParticle* track) const
160 //Checks first particle from both, particle and track pairs
161 AliAODParticleCut* pc = fPairCut->GetFirstPartCut();
162 return (pc->Pass(part))?kTRUE:pc->Pass(track);
164 /*************************************************************************************/
165 inline Bool_t AliHBTAnalysis::PassPartAndTrack2(AliVAODParticle* part,AliVAODParticle* track) const
167 //Checks second particle from both, particle and track pairs
168 AliAODParticleCut* pc = fPairCut->GetSecondPartCut();
169 return (pc->Pass(part))?kTRUE:pc->Pass(track);
171 /*************************************************************************************/