1 /**************************************************************************
2 AliFemtoAnalysis - the most basic analysis there is.
3 Most others (e.g. AliFemtoVertexAnalysis) wrap this one.
4 **************************************************************************/
6 #ifndef ALIFEMTOANALYSIS_H
7 #define ALIFEMTOANALYSIS_H
12 #include "AliFemtoBaseAnalysis.h" // base analysis class
13 #include "AliFemtoPairCut.h"
14 #include "AliFemtoEventCut.h"
15 #include "AliFemtoParticleCut.h"
16 #include "AliFemtoCorrFctn.h"
17 #include "AliFemtoCorrFctnCollection.h"
18 #include "AliFemtoPicoEventCollection.h"
19 #include "AliFemtoParticleCollection.h"
20 #include "AliFemtoPicoEvent.h"
22 class AliFemtoPicoEventCollectionVectorHideAway;
24 class AliFemtoAnalysis : public AliFemtoBaseAnalysis {
26 // friend class AliFemtoLikeSignAnalysis;
30 AliFemtoAnalysis(const AliFemtoAnalysis& OriginalAnalysis); // copy constructor
31 virtual ~AliFemtoAnalysis();
33 AliFemtoAnalysis& operator=(const AliFemtoAnalysis& aAna);
37 virtual AliFemtoPairCut* PairCut();
38 virtual AliFemtoEventCut* EventCut();
39 virtual AliFemtoParticleCut* FirstParticleCut();
40 virtual AliFemtoParticleCut* SecondParticleCut();
42 AliFemtoCorrFctnCollection* CorrFctnCollection();
43 virtual AliFemtoCorrFctn* CorrFctn(int n); // Access to CFs within the collection
44 void AddCorrFctn(AliFemtoCorrFctn* AnotherCorrFctn);
46 void SetPairCut(AliFemtoPairCut* ThePairCut);
47 void SetEventCut(AliFemtoEventCut* TheEventCut);
48 void SetFirstParticleCut(AliFemtoParticleCut* TheFirstParticleCut);
49 void SetSecondParticleCut(AliFemtoParticleCut* TheSecondParticleCut);
51 void SetMinSizePartCollection(unsigned int minSize);
53 unsigned int NumEventsToMix() const;
54 void SetNumEventsToMix(const unsigned int& NumberOfEventsToMix);
55 AliFemtoPicoEvent* CurrentPicoEvent();
56 AliFemtoPicoEventCollection* MixingBuffer();
57 bool MixingBufferFull();
59 bool AnalyzeIdenticalParticles() const;
60 virtual AliFemtoString Report(); //! returns reports of all cuts applied and correlation functions being done
62 virtual void EventBegin(const AliFemtoEvent* TheEventToBegin); // startup for EbyE
63 virtual void ProcessEvent(const AliFemtoEvent* EventToProcess);
64 virtual void EventEnd(const AliFemtoEvent* TheEventToWrapUp); // cleanup for EbyE
65 int GetNeventsProcessed() const;
67 virtual void Finish();
71 void AddEventProcessed();
72 void MakePairs(const char* type,
73 AliFemtoParticleCollection* ParticlesPassingCut1,
74 AliFemtoParticleCollection* ParticlesPssingCut2=0);
76 AliFemtoPicoEventCollectionVectorHideAway* fPicoEventCollectionVectorHideAway; /* Mixing Buffer used for Analyses which wrap this one */
78 AliFemtoPairCut* fPairCut; /* cut applied to pairs */
79 AliFemtoCorrFctnCollection* fCorrFctnCollection; /* correlation functions of this analysis */
80 AliFemtoEventCut* fEventCut; /* cut to select events */
81 AliFemtoParticleCut* fFirstParticleCut; /* select particles of type #1 */
82 AliFemtoParticleCut* fSecondParticleCut; /* select particles of type #2 */
83 AliFemtoPicoEventCollection* fMixingBuffer; /* mixing buffer used in this simplest analysis */
84 AliFemtoPicoEvent* fPicoEvent; /* The current event, in the small (pico) form */
85 unsigned int fNumEventsToMix; /* How many "previous" events get mixed with this one, to make background */
86 unsigned int fNeventsProcessed; /* How many events processed so far */
88 unsigned int fMinSizePartCollection; /* Don't use event if it has fewer than this many particles passing ParticleCuts default 0*/
91 ClassDef(AliFemtoAnalysis, 0)
97 inline AliFemtoPairCut* AliFemtoAnalysis::PairCut() {return fPairCut;}
98 inline AliFemtoEventCut* AliFemtoAnalysis::EventCut() {return fEventCut;}
99 inline AliFemtoParticleCut* AliFemtoAnalysis::FirstParticleCut() {return fFirstParticleCut;}
100 inline AliFemtoParticleCut* AliFemtoAnalysis::SecondParticleCut() {return fSecondParticleCut;}
101 inline AliFemtoCorrFctnCollection* AliFemtoAnalysis::CorrFctnCollection() {return fCorrFctnCollection;}
102 inline unsigned int AliFemtoAnalysis::NumEventsToMix() const {return fNumEventsToMix;}
103 inline AliFemtoPicoEvent* AliFemtoAnalysis::CurrentPicoEvent() {return fPicoEvent;}
105 inline AliFemtoPicoEventCollection* AliFemtoAnalysis::MixingBuffer() {return fMixingBuffer;}
107 inline bool AliFemtoAnalysis::AnalyzeIdenticalParticles() const { return (fFirstParticleCut==fSecondParticleCut); }
110 inline void AliFemtoAnalysis::SetPairCut(AliFemtoPairCut* x) { fPairCut = x; x->SetAnalysis((AliFemtoBaseAnalysis*)this);}
111 inline void AliFemtoAnalysis::AddCorrFctn(AliFemtoCorrFctn* cf) {fCorrFctnCollection->push_back(cf); cf->SetAnalysis((AliFemtoBaseAnalysis*)this);}
112 inline void AliFemtoAnalysis::SetEventCut(AliFemtoEventCut* x) {fEventCut = x; x->SetAnalysis((AliFemtoBaseAnalysis*)this);}
113 inline void AliFemtoAnalysis::SetFirstParticleCut(AliFemtoParticleCut* x) {fFirstParticleCut = x; x->SetAnalysis((AliFemtoBaseAnalysis*)this);}
114 inline void AliFemtoAnalysis::SetSecondParticleCut(AliFemtoParticleCut* x) {fSecondParticleCut = x; x->SetAnalysis((AliFemtoBaseAnalysis*)this);}
116 inline void AliFemtoAnalysis::SetNumEventsToMix(const unsigned int& nmix){ fNumEventsToMix = nmix;}
117 inline bool AliFemtoAnalysis::MixingBufferFull(){return (fMixingBuffer->size() >= fNumEventsToMix);}
118 inline int AliFemtoAnalysis::GetNeventsProcessed() const {return fNeventsProcessed;}
120 inline void AliFemtoAnalysis::SetMinSizePartCollection(unsigned int minSize){fMinSizePartCollection = minSize;}