]>
Commit | Line | Data |
---|---|---|
76ce4b5b | 1 | //////////////////////////////////////////////////////////////////////////////// |
2 | // // | |
3 | // AliFemtoBasicEventCut - the basic cut for events. // | |
4 | // Only cuts on event multiplicity and z-vertex position // | |
5 | // // | |
6 | //////////////////////////////////////////////////////////////////////////////// | |
7 | ||
8 | #include "AliFemtoBasicEventCut.h" | |
9 | //#include <cstdio> | |
10 | ||
11 | #ifdef __ROOT__ | |
12 | ClassImp(AliFemtoBasicEventCut) | |
13 | #endif | |
14 | ||
15 | AliFemtoBasicEventCut::AliFemtoBasicEventCut() : | |
16 | AliFemtoEventCut(), | |
17 | fEventMult(), | |
18 | fVertZPos(), | |
19 | fAcceptBadVertex(false), | |
20 | fNEventsPassed(0), | |
21 | fNEventsFailed(0), | |
76ce4b5b | 22 | fSelectTrigger(0) |
23 | { | |
24 | // Default constructor | |
25 | fEventMult[0] = 0; | |
26 | fEventMult[1] = 100000; | |
27 | fVertZPos[0] = -100.0; | |
28 | fVertZPos[1] = 100.0; | |
d381642d | 29 | fPsiEP[0] = -1000.0; |
30 | fPsiEP[1] = 1000.0; | |
76ce4b5b | 31 | } |
32 | //------------------------------ | |
33 | AliFemtoBasicEventCut::~AliFemtoBasicEventCut(){ | |
34 | // Default destructor | |
35 | } | |
36 | //------------------------------ | |
973a91f8 | 37 | bool AliFemtoBasicEventCut::Pass(const AliFemtoEvent* event){ |
ce7b3d98 | 38 | |
76ce4b5b | 39 | // Pass events if they fall within the multiplicity and z-vertex |
40 | // position range. Fail otherwise | |
41 | // int mult = event->NumberOfTracks(); | |
42 | int mult = (int) event->UncorrectedNumberOfPrimaries(); | |
43 | double vertexZPos = event->PrimVertPos().z(); | |
5e2038a9 | 44 | |
5e2038a9 | 45 | // Double_t qxEPVZERO = 0, qyEPVZERO = 0; |
46 | // Double_t qVZERO = -999; | |
47 | double epvzero = event->ReactionPlaneAngle(); | |
48 | ||
49 | // cout << "AliFemtoBasicEventCut:: epvzero: " << fPsiEP[0] << " < " << epvzero << " < " << fPsiEP[1] << endl; | |
76ce4b5b | 50 | // cout << "AliFemtoBasicEventCut:: mult: " << fEventMult[0] << " < " << mult << " < " << fEventMult[1] << endl; |
51 | // cout << "AliFemtoBasicEventCut:: VertexZPos: " << fVertZPos[0] << " < " << vertexZPos << " < " << fVertZPos[1] << endl; | |
52 | // cout << "AliFemtoBasicEventCut:: VertexZErr: " << event->PrimVertCov()[4] << endl; | |
4eac0b05 | 53 | |
5e2038a9 | 54 | // cout << "AliFemtoBasicEventCut:: MagneticField: " << event->MagneticField() << endl; |
55 | // cout << "AliFemtoBasicEventCut:: IsCollisionCandidate: " << event->IsCollisionCandidate() << endl; | |
56 | // cout << "AliFemtoBasicEventCut:: TriggerCluster: " << event->TriggerCluster() << endl; | |
57 | // cout << "AliFemtoBasicEventCut:: fSelectTrigger: " << fSelectTrigger << endl; | |
58 | // cout << "AliFemtoBasicEventCut:: " << endl; | |
76ce4b5b | 59 | bool goodEvent = |
60 | ((mult >= fEventMult[0]) && | |
61 | (mult <= fEventMult[1]) && | |
62 | (vertexZPos > fVertZPos[0]) && | |
63 | (vertexZPos < fVertZPos[1]) && | |
5e2038a9 | 64 | (epvzero > fPsiEP[0]) && |
65 | (epvzero < fPsiEP[1]) && | |
76ce4b5b | 66 | ((!fAcceptBadVertex) || (event->ZDCParticipants() > 1.0)) && |
5e2038a9 | 67 | ((!fSelectTrigger) || (event->TriggerCluster() == fSelectTrigger)) |
68 | ); | |
69 | ||
70 | // cout << "AliFemtoBasicEventCut:: goodEvent" <<goodEvent << endl; | |
71 | ||
76ce4b5b | 72 | goodEvent ? fNEventsPassed++ : fNEventsFailed++ ; |
7dad907a | 73 | // cout << "AliFemtoBasicEventCut:: return : " << goodEvent << endl; |
76ce4b5b | 74 | // (fAcceptBadVertex || (event->PrimVertCov()[4] > -1000.0)) && |
4eac0b05 | 75 | |
76ce4b5b | 76 | return (goodEvent); |
77 | } | |
78 | //------------------------------ | |
79 | AliFemtoString AliFemtoBasicEventCut::Report(){ | |
80 | // Prepare report | |
81 | string stemp; | |
82 | char ctemp[100]; | |
83 | snprintf(ctemp , 100, "\nMultiplicity:\t %d-%d",fEventMult[0],fEventMult[1]); | |
84 | stemp = ctemp; | |
85 | snprintf(ctemp , 100, "\nVertex Z-position:\t %E-%E",fVertZPos[0],fVertZPos[1]); | |
86 | stemp += ctemp; | |
87 | snprintf(ctemp , 100, "\nNumber of events which passed:\t%ld Number which failed:\t%ld",fNEventsPassed,fNEventsFailed); | |
88 | stemp += ctemp; | |
89 | AliFemtoString returnThis = stemp; | |
90 | return returnThis; | |
91 | } | |
92 | void AliFemtoBasicEventCut::SetAcceptBadVertex(bool b) | |
93 | { | |
94 | fAcceptBadVertex = b; | |
95 | } | |
96 | bool AliFemtoBasicEventCut::GetAcceptBadVertex() | |
97 | { | |
98 | return fAcceptBadVertex; | |
99 | } | |
4eac0b05 | 100 |