1 /***************************************************************************
4 ***************************************************************************
7 ***************************************************************************
10 * Revision 1.1.1.1 2007/04/25 15:38:41 panos
11 * Importing the HBT code dir
13 * Revision 1.1.1.1 2007/03/07 10:14:49 mchojnacki
14 * First version on CVS
17 **************************************************************************/
19 #include "AliFemtoQPairCut.h"
24 ClassImp(AliFemtoQPairCut)
28 AliFemtoQPairCut::AliFemtoQPairCut()
30 fNPairsPassed = fNPairsFailed = 0;
31 fQlong[0]=-1.0; fQlong[1]=100.0;
32 fQout[0]=-1.0; fQout[1]=100.0;
33 fQside[0]=-1.0; fQside[1]=100.0;
34 fQinv[0]=-1.0; fQinv[1]=100.0;
37 AliFemtoQPairCut::~AliFemtoQPairCut()
42 bool AliFemtoQPairCut::Pass(const AliFemtoPair* pair)
45 //temp ? fNPairsPassed++ : fNPairsFailed++;
46 if ((fabs(pair->qLongCMS())<fQlong[0])||(fabs(pair->qLongCMS())>fQlong[1]))
51 if ((fabs(pair->qOutCMS())<fQout[0])||(fabs(pair->qOutCMS())>fQout[1]))
56 if ((fabs(pair->qSideCMS())<fQside[0])||(fabs(pair->qSideCMS())>fQside[1]))
61 if ((fabs(pair->KStar())<fQinv[0])||(fabs(pair->KStar())>fQinv[1]))
70 AliFemtoString AliFemtoQPairCut::Report()
72 string Stemp = "AliFemtoQ Pair Cut \n";
74 sprintf(Ctemp,"Number of pairs which passed:\t%ld Number which failed:\t%ld\n",fNPairsPassed,fNPairsFailed);
76 AliFemtoString returnThis = Stemp;