[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemto / AliFemtoKTPairCut.cxx

/////////////////////////////////////////////////////////////////////////////
//                                                                         //
// AliFemtoKTPairCut - a pair cut which selects pairs based on their       //
// transverse momentum kT                                                  //
//                                                                         //
/////////////////////////////////////////////////////////////////////////////
/***************************************************************************
 *
 * $Id: AliFemtoKTPairCut.cxx,v 1.1.2.2 2007/11/09 11:20:35 akisiel Exp $
 *
 * Author: Adam Kisiel, Ohio State, kisiel@mps.ohio-state.edu
 ***************************************************************************
 *
 * Description: part of STAR HBT Framework: AliFemtoMaker package
 *   a cut to remove "shared" and "split" pairs
 *
 ***************************************************************************
 *
 *
 **************************************************************************/

#include "AliFemtoKTPairCut.h"
#include <string>
#include <cstdio>

#ifdef __ROOT__
ClassImp(AliFemtoKTPairCut)
#endif

//__________________
AliFemtoKTPairCut::AliFemtoKTPairCut():
  AliFemtoPairCut(),
  fKTMin(0),
  fKTMax(1.0e6)
{
  fKTMin = 0;
   fKTMax = 1.0e6;
}
//__________________
AliFemtoKTPairCut::AliFemtoKTPairCut(double lo, double hi) :
  AliFemtoPairCut(),
  fKTMin(lo),
  fKTMax(hi)
{
}
//__________________
AliFemtoKTPairCut::AliFemtoKTPairCut(const AliFemtoKTPairCut& c) : 
  AliFemtoPairCut(c),
  fKTMin(0),
  fKTMax(1.0e6)
{ 
  fKTMin = c.fKTMin;
  fKTMax = c.fKTMax;
}

//__________________
AliFemtoKTPairCut::~AliFemtoKTPairCut(){
  /* no-op */
}
//__________________
bool AliFemtoKTPairCut::Pass(const AliFemtoPair* pair){
  bool temp = true;
  
  if (pair->KT() < fKTMin)
    temp = false;

  if (pair->KT() > fKTMax)
    temp = false;

  return temp;
}
//__________________
AliFemtoString AliFemtoKTPairCut::Report(){
  // Prepare a report from the execution
  string stemp = "AliFemtoKT Pair Cut \n";  char ctemp[100];
  sprintf(ctemp,"Accept pair with kT in range %lf , %lf",fKTMin,fKTMax);
  stemp += ctemp;
  AliFemtoString returnThis = stemp;
  return returnThis;}
//__________________

TList *AliFemtoKTPairCut::ListSettings()
{
  // return a list of settings in a writable form
  TList *tListSetttings =  new TList();
  char buf[200];
  snprintf(buf, 200, "AliFemtoKTPairCut.ktmax=%lf", fKTMax);
  tListSetttings->AddLast(new TObjString(buf));
  snprintf(buf, 200, "AliFemtoKTPairCut.ktmin=%lf", fKTMin);
  tListSetttings->AddLast(new TObjString(buf));

  return tListSetttings;
}

void AliFemtoKTPairCut::SetKTRange(double ktmin, double ktmax)
{
  fKTMin = ktmin;
  fKTMax = ktmax;
}
Commit	Line	Data
0b3bd1ac	1	/////////////////////////////////////////////////////////////////////////////
	2	// //
	3	// AliFemtoKTPairCut - a pair cut which selects pairs based on their //
	4	// transverse momentum kT //
	5	// //
	6	/////////////////////////////////////////////////////////////////////////////
	7	/***************************************************************************
	8	*
	9	* $Id: AliFemtoKTPairCut.cxx,v 1.1.2.2 2007/11/09 11:20:35 akisiel Exp $
	10	*
	11	* Author: Adam Kisiel, Ohio State, kisiel@mps.ohio-state.edu
	12	***************************************************************************
	13	*
	14	* Description: part of STAR HBT Framework: AliFemtoMaker package
	15	* a cut to remove "shared" and "split" pairs
	16	*
	17	***************************************************************************
	18	*
	19	*
	20	**************************************************************************/
	21
	22	#include "AliFemtoKTPairCut.h"
	23	#include <string>
	24	#include <cstdio>
	25
	26	#ifdef __ROOT__
	27	ClassImp(AliFemtoKTPairCut)
	28	#endif
	29
	30	//__________________
	31	AliFemtoKTPairCut::AliFemtoKTPairCut():
	32	AliFemtoPairCut(),
	33	fKTMin(0),
	34	fKTMax(1.0e6)
	35	{
	36	fKTMin = 0;
	37	fKTMax = 1.0e6;
	38	}
	39	//__________________
	40	AliFemtoKTPairCut::AliFemtoKTPairCut(double lo, double hi) :
	41	AliFemtoPairCut(),
	42	fKTMin(lo),
	43	fKTMax(hi)
	44	{
	45	}
	46	//__________________
	47	AliFemtoKTPairCut::AliFemtoKTPairCut(const AliFemtoKTPairCut& c) :
	48	AliFemtoPairCut(c),
	49	fKTMin(0),
	50	fKTMax(1.0e6)
	51	{
	52	fKTMin = c.fKTMin;
	53	fKTMax = c.fKTMax;
	54	}
	55
	56	//__________________
	57	AliFemtoKTPairCut::~AliFemtoKTPairCut(){
	58	/* no-op */
	59	}
	60	//__________________
	61	bool AliFemtoKTPairCut::Pass(const AliFemtoPair* pair){
	62	bool temp = true;
	63
	64	if (pair->KT() < fKTMin)
65	temp = false;
66
67	if (pair->KT() > fKTMax)
68	temp = false;
69
70	return temp;
71	}
72	//__________________
73	AliFemtoString AliFemtoKTPairCut::Report(){
74	// Prepare a report from the execution
75	string stemp = "AliFemtoKT Pair Cut \n"; char ctemp[100];
76	sprintf(ctemp,"Accept pair with kT in range %lf , %lf",fKTMin,fKTMax);
77	stemp += ctemp;
78	AliFemtoString returnThis = stemp;
79	return returnThis;}
80	//__________________
81
82	TList *AliFemtoKTPairCut::ListSettings()
83	{
84	// return a list of settings in a writable form
85	TList *tListSetttings = new TList();
86	char buf[200];
87	snprintf(buf, 200, "AliFemtoKTPairCut.ktmax=%lf", fKTMax);
88	tListSetttings->AddLast(new TObjString(buf));
89	snprintf(buf, 200, "AliFemtoKTPairCut.ktmin=%lf", fKTMin);
90	tListSetttings->AddLast(new TObjString(buf));
91
92	return tListSetttings;
93	}
94
95	void AliFemtoKTPairCut::SetKTRange(double ktmin, double ktmax)
96	{
97	fKTMin = ktmin;
98	fKTMax = ktmax;
99	}