1 /**************************************************************************
2 AliFemtoSimpleAnalysis - the most basic analysis there is.
3 Most others (e.g. AliFemtoVertexAnalysis) wrap this one.
4 **************************************************************************/
6 #ifndef ALIFEMTO_SIMPLE_ANALYSIS_H
7 #define ALIFEMTO_SIMPLE_ANALYSIS_H
12 #include "AliFemtoAnalysis.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 AliFemtoSimpleAnalysis : public AliFemtoAnalysis {
26 // friend class AliFemtoLikeSignAnalysis;
29 AliFemtoSimpleAnalysis();
30 AliFemtoSimpleAnalysis(const AliFemtoSimpleAnalysis& OriginalAnalysis); // copy constructor
31 virtual ~AliFemtoSimpleAnalysis();
33 AliFemtoSimpleAnalysis& operator=(const AliFemtoSimpleAnalysis& 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
61 virtual TList* ListSettings(); // return list of cut settings for the analysis
63 virtual void EventBegin(const AliFemtoEvent* TheEventToBegin); // startup for EbyE
64 virtual void ProcessEvent(const AliFemtoEvent* EventToProcess);
65 virtual void EventEnd(const AliFemtoEvent* TheEventToWrapUp); // cleanup for EbyE
66 int GetNeventsProcessed() const;
68 virtual void Finish();
71 ClassDef(AliFemtoSimpleAnalysis, 0)
77 void AddEventProcessed();
78 void MakePairs(const char* type,
79 AliFemtoParticleCollection* ParticlesPassingCut1,
80 AliFemtoParticleCollection* ParticlesPssingCut2=0);
82 AliFemtoPicoEventCollectionVectorHideAway* fPicoEventCollectionVectorHideAway; /* Mixing Buffer used for Analyses which wrap this one */
84 AliFemtoPairCut* fPairCut; /* cut applied to pairs */
85 AliFemtoCorrFctnCollection* fCorrFctnCollection; /* correlation functions of this analysis */
86 AliFemtoEventCut* fEventCut; /* cut to select events */
87 AliFemtoParticleCut* fFirstParticleCut; /* select particles of type #1 */
88 AliFemtoParticleCut* fSecondParticleCut; /* select particles of type #2 */
89 AliFemtoPicoEventCollection* fMixingBuffer; /* mixing buffer used in this simplest analysis */
90 AliFemtoPicoEvent* fPicoEvent; /* The current event, in the small (pico) form */
91 unsigned int fNumEventsToMix; /* How many "previous" events get mixed with this one, to make background */
92 unsigned int fNeventsProcessed; /* How many events processed so far */
94 unsigned int fMinSizePartCollection; /* Don't use event if it has fewer than this many particles passing ParticleCuts default 0*/
97 ClassDef(AliFemtoSimpleAnalysis, 0)
103 inline AliFemtoPairCut* AliFemtoSimpleAnalysis::PairCut() {return fPairCut;}
104 inline AliFemtoEventCut* AliFemtoSimpleAnalysis::EventCut() {return fEventCut;}
105 inline AliFemtoParticleCut* AliFemtoSimpleAnalysis::FirstParticleCut() {return fFirstParticleCut;}
106 inline AliFemtoParticleCut* AliFemtoSimpleAnalysis::SecondParticleCut() {return fSecondParticleCut;}
107 inline AliFemtoCorrFctnCollection* AliFemtoSimpleAnalysis::CorrFctnCollection() {return fCorrFctnCollection;}
108 inline unsigned int AliFemtoSimpleAnalysis::NumEventsToMix() const {return fNumEventsToMix;}
109 inline AliFemtoPicoEvent* AliFemtoSimpleAnalysis::CurrentPicoEvent() {return fPicoEvent;}
111 inline AliFemtoPicoEventCollection* AliFemtoSimpleAnalysis::MixingBuffer() {return fMixingBuffer;}
113 inline bool AliFemtoSimpleAnalysis::AnalyzeIdenticalParticles() const { return (fFirstParticleCut==fSecondParticleCut); }
116 inline void AliFemtoSimpleAnalysis::SetPairCut(AliFemtoPairCut* x) { fPairCut = x; x->SetAnalysis((AliFemtoAnalysis*)this);}
117 inline void AliFemtoSimpleAnalysis::AddCorrFctn(AliFemtoCorrFctn* cf) {fCorrFctnCollection->push_back(cf); cf->SetAnalysis((AliFemtoAnalysis*)this);}
118 inline void AliFemtoSimpleAnalysis::SetEventCut(AliFemtoEventCut* x) {fEventCut = x; x->SetAnalysis((AliFemtoAnalysis*)this);}
119 inline void AliFemtoSimpleAnalysis::SetFirstParticleCut(AliFemtoParticleCut* x) {fFirstParticleCut = x; x->SetAnalysis((AliFemtoAnalysis*)this);}
120 inline void AliFemtoSimpleAnalysis::SetSecondParticleCut(AliFemtoParticleCut* x) {fSecondParticleCut = x; x->SetAnalysis((AliFemtoAnalysis*)this);}
122 inline void AliFemtoSimpleAnalysis::SetNumEventsToMix(const unsigned int& nmix){ fNumEventsToMix = nmix;}
123 inline bool AliFemtoSimpleAnalysis::MixingBufferFull(){return (fMixingBuffer->size() >= fNumEventsToMix);}
124 inline int AliFemtoSimpleAnalysis::GetNeventsProcessed() const {return fNeventsProcessed;}
126 inline void AliFemtoSimpleAnalysis::SetMinSizePartCollection(unsigned int minSize){fMinSizePartCollection = minSize;}