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

////////////////////////////////////////////////////////////////////////////////
//                                                                            //
// AliFemtoCorrFctnTPCNcls - A correlation function that saves the correlation//
// function as a function of number of TPC clusters of the track              //
//                                                                            //
// Authors: Adam Kisiel kisiel@mps.ohio-state.edu                             //
//                                                                            //
////////////////////////////////////////////////////////////////////////////////

#include "AliFemtoCorrFctnTPCNcls.h"
//#include "AliFemtoHisto.hh"
#include <cstdio>

#ifdef __ROOT__ 
ClassImp(AliFemtoCorrFctnTPCNcls)
#endif

//____________________________
AliFemtoCorrFctnTPCNcls::AliFemtoCorrFctnTPCNcls(char* title, const int& nbins, const float& QinvLo, const float& QinvHi):
  AliFemtoCorrFctn(),
  fNclsTPCMinNumerator(0),
  fNclsTPCMinDenominator(0)
{
  // set up numerator
  char tTitNum[101] = "NumNclsTPCMin";
  strncat(tTitNum,title, 100);
  fNclsTPCMinNumerator = new TH2D(tTitNum,title,nbins,QinvLo,QinvHi,159,0.5,159.5);
  // set up denominator
  char tTitDen[101] = "DenNclsTPCMin";
  strncat(tTitDen,title, 100);
  fNclsTPCMinDenominator = new TH2D(tTitDen,title,nbins,QinvLo,QinvHi,159,0.5,159.5);

  // to enable error bar calculation...
  fNclsTPCMinNumerator->Sumw2();
  fNclsTPCMinDenominator->Sumw2();
}

//____________________________
AliFemtoCorrFctnTPCNcls::AliFemtoCorrFctnTPCNcls(const AliFemtoCorrFctnTPCNcls& aCorrFctn) :
  AliFemtoCorrFctn(),
  fNclsTPCMinNumerator(0),
  fNclsTPCMinDenominator(0)
{
  // copy constructor
  if (aCorrFctn.fNclsTPCMinNumerator)
    fNclsTPCMinNumerator = new TH2D(*aCorrFctn.fNclsTPCMinNumerator);
  if (aCorrFctn.fNclsTPCMinDenominator)
    fNclsTPCMinDenominator = new TH2D(*aCorrFctn.fNclsTPCMinDenominator);
}
//____________________________
AliFemtoCorrFctnTPCNcls::~AliFemtoCorrFctnTPCNcls(){
  // destructor
  delete fNclsTPCMinNumerator;
  delete fNclsTPCMinDenominator;
}
//_________________________
AliFemtoCorrFctnTPCNcls& AliFemtoCorrFctnTPCNcls::operator=(const AliFemtoCorrFctnTPCNcls& aCorrFctn)
{
  // assignment operator
  if (this == &aCorrFctn)
    return *this;

  if (aCorrFctn.fNclsTPCMinNumerator)
    fNclsTPCMinNumerator = new TH2D(*aCorrFctn.fNclsTPCMinNumerator);
  else
    fNclsTPCMinNumerator = 0;
  if (aCorrFctn.fNclsTPCMinDenominator)
    fNclsTPCMinDenominator = new TH2D(*aCorrFctn.fNclsTPCMinDenominator);
  else
    fNclsTPCMinDenominator = 0;

  return *this;
}
//_________________________
void AliFemtoCorrFctnTPCNcls::Finish(){
  // here is where we should normalize, fit, etc...
  // we should NOT Draw() the histos (as I had done it below),
  // since we want to insulate ourselves from root at this level
  // of the code.  Do it instead at root command line with browser.
  //  mShareNumerator->Draw();
  //mShareDenominator->Draw();
  //mRatio->Draw();

}

//____________________________
AliFemtoString AliFemtoCorrFctnTPCNcls::Report(){
  // create report
  string stemp = "TPC Ncls Correlation Function Report:\n";
  char ctemp[100];
  snprintf(ctemp , 100, "Number of entries in numerator:\t%E\n",fNclsTPCMinNumerator->GetEntries());
  stemp += ctemp;
  snprintf(ctemp , 100, "Number of entries in denominator:\t%E\n",fNclsTPCMinDenominator->GetEntries());
  stemp += ctemp;
  //  stemp += mCoulombWeight->Report();
  AliFemtoString returnThis = stemp;
  return returnThis;
}
//____________________________
void AliFemtoCorrFctnTPCNcls::AddRealPair( AliFemtoPair* pair){
  // add real (effect) pair
  double tQinv = fabs(pair->QInv());   // note - qInv() will be negative for identical pairs...

  if (pair->Track1()->Track()->TPCncls()>pair->Track2()->Track()->TPCncls())
    fNclsTPCMinNumerator->Fill(tQinv, pair->Track2()->Track()->TPCncls());
  else
    fNclsTPCMinNumerator->Fill(tQinv, pair->Track1()->Track()->TPCncls());
}
//____________________________
void AliFemtoCorrFctnTPCNcls::AddMixedPair( AliFemtoPair* pair){
  // add mixed (background) pair
  double tQinv = fabs(pair->QInv());   // note - qInv() will be negative for identical pairs...

  if (pair->Track1()->Track()->TPCncls()>pair->Track2()->Track()->TPCncls())
    fNclsTPCMinDenominator->Fill(tQinv, pair->Track2()->Track()->TPCncls());
  else
    fNclsTPCMinDenominator->Fill(tQinv, pair->Track1()->Track()->TPCncls());
}


void AliFemtoCorrFctnTPCNcls::WriteHistos()
{
  // Write out result histograms
  fNclsTPCMinNumerator->Write();
  fNclsTPCMinDenominator->Write();
}

TList* AliFemtoCorrFctnTPCNcls::GetOutputList()
{
  // Prepare the list of objects to be written to the output
  TList *tOutputList = new TList();

  tOutputList->Add(fNclsTPCMinNumerator);
  tOutputList->Add(fNclsTPCMinDenominator);
  
  return tOutputList;

}
Commit	Line	Data
76ce4b5b	1	////////////////////////////////////////////////////////////////////////////////
	2	// //
	3	// AliFemtoCorrFctnTPCNcls - A correlation function that saves the correlation//
	4	// function as a function of number of TPC clusters of the track //
	5	// //
	6	// Authors: Adam Kisiel kisiel@mps.ohio-state.edu //
	7	// //
	8	////////////////////////////////////////////////////////////////////////////////
	9
	10	#include "AliFemtoCorrFctnTPCNcls.h"
	11	//#include "AliFemtoHisto.hh"
	12	#include <cstdio>
	13
	14	#ifdef __ROOT__
	15	ClassImp(AliFemtoCorrFctnTPCNcls)
	16	#endif
	17
	18	//____________________________
	19	AliFemtoCorrFctnTPCNcls::AliFemtoCorrFctnTPCNcls(char* title, const int& nbins, const float& QinvLo, const float& QinvHi):
	20	AliFemtoCorrFctn(),
	21	fNclsTPCMinNumerator(0),
	22	fNclsTPCMinDenominator(0)
	23	{
	24	// set up numerator
	25	char tTitNum[101] = "NumNclsTPCMin";
	26	strncat(tTitNum,title, 100);
	27	fNclsTPCMinNumerator = new TH2D(tTitNum,title,nbins,QinvLo,QinvHi,159,0.5,159.5);
	28	// set up denominator
	29	char tTitDen[101] = "DenNclsTPCMin";
	30	strncat(tTitDen,title, 100);
	31	fNclsTPCMinDenominator = new TH2D(tTitDen,title,nbins,QinvLo,QinvHi,159,0.5,159.5);
	32
	33	// to enable error bar calculation...
	34	fNclsTPCMinNumerator->Sumw2();
	35	fNclsTPCMinDenominator->Sumw2();
	36	}
	37
	38	//____________________________
	39	AliFemtoCorrFctnTPCNcls::AliFemtoCorrFctnTPCNcls(const AliFemtoCorrFctnTPCNcls& aCorrFctn) :
	40	AliFemtoCorrFctn(),
	41	fNclsTPCMinNumerator(0),
	42	fNclsTPCMinDenominator(0)
	43	{
	44	// copy constructor
	45	if (aCorrFctn.fNclsTPCMinNumerator)
	46	fNclsTPCMinNumerator = new TH2D(*aCorrFctn.fNclsTPCMinNumerator);
	47	if (aCorrFctn.fNclsTPCMinDenominator)
	48	fNclsTPCMinDenominator = new TH2D(*aCorrFctn.fNclsTPCMinDenominator);
	49	}
	50	//____________________________
	51	AliFemtoCorrFctnTPCNcls::~AliFemtoCorrFctnTPCNcls(){
	52	// destructor
	53	delete fNclsTPCMinNumerator;
	54	delete fNclsTPCMinDenominator;
	55	}
	56	//_________________________
	57	AliFemtoCorrFctnTPCNcls& AliFemtoCorrFctnTPCNcls::operator=(const AliFemtoCorrFctnTPCNcls& aCorrFctn)
	58	{
	59	// assignment operator
	60	if (this == &aCorrFctn)
	61	return *this;
	62
	63	if (aCorrFctn.fNclsTPCMinNumerator)
	64	fNclsTPCMinNumerator = new TH2D(*aCorrFctn.fNclsTPCMinNumerator);
65	else
66	fNclsTPCMinNumerator = 0;
67	if (aCorrFctn.fNclsTPCMinDenominator)
68	fNclsTPCMinDenominator = new TH2D(*aCorrFctn.fNclsTPCMinDenominator);
69	else
70	fNclsTPCMinDenominator = 0;
71
72	return *this;
73	}
74	//_________________________
75	void AliFemtoCorrFctnTPCNcls::Finish(){
76	// here is where we should normalize, fit, etc...
77	// we should NOT Draw() the histos (as I had done it below),
78	// since we want to insulate ourselves from root at this level
79	// of the code. Do it instead at root command line with browser.
80	// mShareNumerator->Draw();
81	//mShareDenominator->Draw();
82	//mRatio->Draw();
83
84	}
85
86	//____________________________
87	AliFemtoString AliFemtoCorrFctnTPCNcls::Report(){
88	// create report
89	string stemp = "TPC Ncls Correlation Function Report:\n";
90	char ctemp[100];
91	snprintf(ctemp , 100, "Number of entries in numerator:\t%E\n",fNclsTPCMinNumerator->GetEntries());
92	stemp += ctemp;
93	snprintf(ctemp , 100, "Number of entries in denominator:\t%E\n",fNclsTPCMinDenominator->GetEntries());
94	stemp += ctemp;
95	// stemp += mCoulombWeight->Report();
96	AliFemtoString returnThis = stemp;
97	return returnThis;
98	}
99	//____________________________
100	void AliFemtoCorrFctnTPCNcls::AddRealPair( AliFemtoPair* pair){
101	// add real (effect) pair
102	double tQinv = fabs(pair->QInv()); // note - qInv() will be negative for identical pairs...
103
104	if (pair->Track1()->Track()->TPCncls()>pair->Track2()->Track()->TPCncls())
105	fNclsTPCMinNumerator->Fill(tQinv, pair->Track2()->Track()->TPCncls());
106	else
107	fNclsTPCMinNumerator->Fill(tQinv, pair->Track1()->Track()->TPCncls());
108	}
109	//____________________________
110	void AliFemtoCorrFctnTPCNcls::AddMixedPair( AliFemtoPair* pair){
111	// add mixed (background) pair
112	double tQinv = fabs(pair->QInv()); // note - qInv() will be negative for identical pairs...
113
114	if (pair->Track1()->Track()->TPCncls()>pair->Track2()->Track()->TPCncls())
115	fNclsTPCMinDenominator->Fill(tQinv, pair->Track2()->Track()->TPCncls());
116	else
117	fNclsTPCMinDenominator->Fill(tQinv, pair->Track1()->Track()->TPCncls());
118	}
119
120
121	void AliFemtoCorrFctnTPCNcls::WriteHistos()
122	{
123	// Write out result histograms
124	fNclsTPCMinNumerator->Write();
125	fNclsTPCMinDenominator->Write();
126	}
127
128	TList* AliFemtoCorrFctnTPCNcls::GetOutputList()
129	{
130	// Prepare the list of objects to be written to the output
131	TList *tOutputList = new TList();
132
133	tOutputList->Add(fNclsTPCMinNumerator);
134	tOutputList->Add(fNclsTPCMinDenominator);
135
136	return tOutputList;
137
138	}