1 /////////////////////////////////////////////////////////////////////////////
3 // AliFemtoQPairCut - a simple cut which selects pairs based on the values //
4 // of their respective q components /
6 /////////////////////////////////////////////////////////////////////////////
7 /***************************************************************************
10 ***************************************************************************
13 ***************************************************************************
16 * Revision 1.2.6.1 2007/11/01 17:10:38 akisiel
17 * Fix code rule conformace
19 * Revision 1.2 2007/05/22 09:01:42 akisiel
20 * Add the possibiloity to save cut settings in the ROOT file
22 * Revision 1.1 2007/05/16 10:25:06 akisiel
23 * Making the directory structure of AliFemtoUser flat. All files go into one common directory
25 * Revision 1.4 2007/05/03 09:46:10 akisiel
26 * Fixing Effective C++ warnings
28 * Revision 1.3 2007/04/27 07:25:59 akisiel
29 * Make revisions needed for compilation from the main AliRoot tree
31 * Revision 1.1.1.1 2007/04/25 15:38:41 panos
32 * Importing the HBT code dir
34 * Revision 1.1.1.1 2007/03/07 10:14:49 mchojnacki
35 * First version on CVS
38 **************************************************************************/
40 #include "AliFemtoQPairCut.h"
45 ClassImp(AliFemtoQPairCut)
49 AliFemtoQPairCut::AliFemtoQPairCut():
53 // Default constructor
54 fNPairsPassed = fNPairsFailed = 0;
55 fQlong[0]=-1.0; fQlong[1]=100.0;
56 fQout[0]=-1.0; fQout[1]=100.0;
57 fQside[0]=-1.0; fQside[1]=100.0;
58 fQinv[0]=-1.0; fQinv[1]=100.0;
61 AliFemtoQPairCut::~AliFemtoQPairCut()
66 bool AliFemtoQPairCut::Pass(const AliFemtoPair* pair)
68 // Select pairs based on their q values
70 //temp ? fNPairsPassed++ : fNPairsFailed++;
71 if ((fabs(pair->QLongCMS())<fQlong[0])||(fabs(pair->QLongCMS())>fQlong[1]))
76 if ((fabs(pair->QOutCMS())<fQout[0])||(fabs(pair->QOutCMS())>fQout[1]))
81 if ((fabs(pair->QSideCMS())<fQside[0])||(fabs(pair->QSideCMS())>fQside[1]))
86 if ((fabs(pair->KStar())<fQinv[0])||(fabs(pair->KStar())>fQinv[1]))
95 AliFemtoString AliFemtoQPairCut::Report()
98 string stemp = "AliFemtoQ Pair Cut \n";
100 snprintf(ctemp , 100, "Number of pairs which passed:\t%ld Number which failed:\t%ld\n",fNPairsPassed,fNPairsFailed);
102 AliFemtoString returnThis = stemp;
106 TList *AliFemtoQPairCut::ListSettings()
108 // return a list of settings in a writable form
109 TList *tListSetttings = new TList();
111 snprintf(buf, 200, "AliFemtoQPairCut.qout.maximum=%f", fQout[0]);
112 tListSetttings->AddLast(new TObjString(buf));
114 snprintf(buf, 200, "AliFemtoQPairCut.qout.minimum=%f", fQout[1]);
115 tListSetttings->AddLast(new TObjString(buf));
117 snprintf(buf, 200, "AliFemtoQPairCut.qside.maximum=%f", fQside[0]);
118 tListSetttings->AddLast(new TObjString(buf));
120 snprintf(buf, 200, "AliFemtoQPairCut.qside.minimum=%f", fQside[1]);
121 tListSetttings->AddLast(new TObjString(buf));
123 snprintf(buf, 200, "AliFemtoQPairCut.qlong.maximum=%f", fQlong[0]);
124 tListSetttings->AddLast(new TObjString(buf));
126 snprintf(buf, 200, "AliFemtoQPairCut.qlong.minimum=%f", fQlong[1]);
127 tListSetttings->AddLast(new TObjString(buf));
129 snprintf(buf, 200, "AliFemtoQPairCut.qinv.maximum=%f", fQinv[0]);
130 tListSetttings->AddLast(new TObjString(buf));
132 snprintf(buf, 200, "AliFemtoQPairCut.qinv.minimum=%f", fQinv[1]);
133 tListSetttings->AddLast(new TObjString(buf));
135 return tListSetttings;