[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / AliFemtoUser / AliFemtoPairCutPt.cxx

/////////////////////////////////////////////////////////////////////////////////////
//                                                                                 //
// AliFemtoPairCutPt - a pair cut which checks if the sum of the transverse        //
// momenta of two particles fit within given range 		                   //
// Authors: Malgorzata Janik, Warsaw University of Technology, majanik@cern.ch     //
//          Lukasz Graczykowski, Warsaw University of Technology, lgraczyk@cern.ch //
//  	       				                                           //
/////////////////////////////////////////////////////////////////////////////////////

#include "AliFemtoPairCutPt.h"
#include <string>
#include <cstdio>
#include <TMath.h>

#ifdef __ROOT__
ClassImp(AliFemtoPairCutPt)
#endif

//__________________
AliFemtoPairCutPt::AliFemtoPairCutPt():
  AliFemtoPairCut(),
  fSumPtMin(0),
  fSumPtMax(10000),
  fNPairsFailed(0),
  fNPairsPassed(0)
{

}
//__________________
AliFemtoPairCutPt::AliFemtoPairCutPt(double lo, double hi):
  AliFemtoPairCut(),
  fSumPtMin(lo),
  fSumPtMax(hi),
  fNPairsFailed(0),
  fNPairsPassed(0)
{
  fSumPtMin=lo;
  fSumPtMax=hi;
}
//__________________
AliFemtoPairCutPt::AliFemtoPairCutPt(const AliFemtoPairCutPt& c) : 
  AliFemtoPairCut(c),
  fSumPtMin(0),
  fSumPtMax(0),
  fNPairsFailed(0),
  fNPairsPassed(0)
{ 
  fSumPtMin = c.fSumPtMin;
  fSumPtMax = c.fSumPtMax;
}
AliFemtoPairCutPt& AliFemtoPairCutPt::operator=(const AliFemtoPairCutPt& c)
{
  if (this != &c) {
    fSumPtMin = c.fSumPtMin;
    fSumPtMax = c.fSumPtMax;
  }

  return *this;

}

//__________________
AliFemtoPairCutPt::~AliFemtoPairCutPt(){
  /* no-op */
}
//__________________
bool AliFemtoPairCutPt::Pass(const AliFemtoPair* pair){

  bool temp = true;

  double pt1 = pair->Track1()->Track()->Pt();
  double pt2 = pair->Track2()->Track()->Pt();

  double pt_sum = pt1 + pt2;

  if(pt_sum >= fSumPtMin && pt_sum <= fSumPtMax)
    temp = true;
  else
    temp = false;

  if(temp) 
    fNPairsPassed++;
  else fNPairsFailed++;


  return temp;
  
}
//__________________
AliFemtoString AliFemtoPairCutPt::Report(){
  // Prepare a report from the execution
  string stemp = "AliFemtoPairCutPt Pair Cut\n";  
  char ctemp[100];
  stemp += ctemp;
  snprintf(ctemp,100,"Number of pairs which passed:\t%ld  Number which failed:\t%ld\n",(long int) fNPairsPassed,(long int) fNPairsFailed);
  stemp += ctemp;
  AliFemtoString returnThis = stemp;
  return returnThis;}
//__________________

TList *AliFemtoPairCutPt::ListSettings()
{
  // return a list of settings in a writable form
  TList *tListSetttings = new TList();
  char buf[200];
  snprintf(buf, 200, "AliFemtoPairCutPt.sumptmin=%f", fSumPtMin);
  snprintf(buf, 200, "AliFemtoPairCutPr.sumptmax=%f", fSumPtMax);
  tListSetttings->AddLast(new TObjString(buf));

  return tListSetttings;
}

void AliFemtoPairCutPt::SetMinSumPt(Double_t sumptmin)
{
  fSumPtMin = sumptmin;
}

 
void AliFemtoPairCutPt::SetMaxSumPt(Double_t sumptmax)
{
  fSumPtMax = sumptmax;
}
Commit	Line	Data
76ce4b5b	1	/////////////////////////////////////////////////////////////////////////////////////
	2	// //
	3	// AliFemtoPairCutPt - a pair cut which checks if the sum of the transverse //
	4	// momenta of two particles fit within given range //
	5	// Authors: Malgorzata Janik, Warsaw University of Technology, majanik@cern.ch //
	6	// Lukasz Graczykowski, Warsaw University of Technology, lgraczyk@cern.ch //
	7	// //
	8	/////////////////////////////////////////////////////////////////////////////////////
	9
	10	#include "AliFemtoPairCutPt.h"
	11	#include <string>
	12	#include <cstdio>
	13	#include <TMath.h>
	14
	15	#ifdef __ROOT__
	16	ClassImp(AliFemtoPairCutPt)
	17	#endif
	18
	19	//__________________
	20	AliFemtoPairCutPt::AliFemtoPairCutPt():
	21	AliFemtoPairCut(),
	22	fSumPtMin(0),
	23	fSumPtMax(10000),
	24	fNPairsFailed(0),
	25	fNPairsPassed(0)
	26	{
	27
	28	}
	29	//__________________
	30	AliFemtoPairCutPt::AliFemtoPairCutPt(double lo, double hi):
	31	AliFemtoPairCut(),
	32	fSumPtMin(lo),
	33	fSumPtMax(hi),
	34	fNPairsFailed(0),
	35	fNPairsPassed(0)
	36	{
	37	fSumPtMin=lo;
	38	fSumPtMax=hi;
	39	}
	40	//__________________
	41	AliFemtoPairCutPt::AliFemtoPairCutPt(const AliFemtoPairCutPt& c) :
	42	AliFemtoPairCut(c),
	43	fSumPtMin(0),
	44	fSumPtMax(0),
	45	fNPairsFailed(0),
	46	fNPairsPassed(0)
	47	{
	48	fSumPtMin = c.fSumPtMin;
	49	fSumPtMax = c.fSumPtMax;
	50	}
	51	AliFemtoPairCutPt& AliFemtoPairCutPt::operator=(const AliFemtoPairCutPt& c)
	52	{
	53	if (this != &c) {
	54	fSumPtMin = c.fSumPtMin;
	55	fSumPtMax = c.fSumPtMax;
	56	}
	57
	58	return *this;
	59
	60	}
	61
	62	//__________________
	63	AliFemtoPairCutPt::~AliFemtoPairCutPt(){
	64	/* no-op */
65	}
66	//__________________
67	bool AliFemtoPairCutPt::Pass(const AliFemtoPair* pair){
68
69	bool temp = true;
70
71	double pt1 = pair->Track1()->Track()->Pt();
72	double pt2 = pair->Track2()->Track()->Pt();
73
74	double pt_sum = pt1 + pt2;
75
76	if(pt_sum >= fSumPtMin && pt_sum <= fSumPtMax)
77	temp = true;
78	else
79	temp = false;
80
81	if(temp)
82	fNPairsPassed++;
83	else fNPairsFailed++;
84
85
86	return temp;
87
88	}
89	//__________________
90	AliFemtoString AliFemtoPairCutPt::Report(){
91	// Prepare a report from the execution
92	string stemp = "AliFemtoPairCutPt Pair Cut\n";
93	char ctemp[100];
94	stemp += ctemp;
95	snprintf(ctemp,100,"Number of pairs which passed:\t%ld Number which failed:\t%ld\n",(long int) fNPairsPassed,(long int) fNPairsFailed);
96	stemp += ctemp;
97	AliFemtoString returnThis = stemp;
98	return returnThis;}
99	//__________________
100
101	TList *AliFemtoPairCutPt::ListSettings()
102	{
103	// return a list of settings in a writable form
104	TList *tListSetttings = new TList();
105	char buf[200];
106	snprintf(buf, 200, "AliFemtoPairCutPt.sumptmin=%f", fSumPtMin);
107	snprintf(buf, 200, "AliFemtoPairCutPr.sumptmax=%f", fSumPtMax);
108	tListSetttings->AddLast(new TObjString(buf));
109
110	return tListSetttings;
111	}
112
113	void AliFemtoPairCutPt::SetMinSumPt(Double_t sumptmin)
114	{
115	fSumPtMin = sumptmin;
116	}
117
118
119	void AliFemtoPairCutPt::SetMaxSumPt(Double_t sumptmax)
120	{
121	fSumPtMax = sumptmax;
122	}
123