[u/mrichter/AliRoot.git] / HBTAN / AliHBTCorrelFctn.cxx

#include "AliHBTCorrelFctn.h"
//____________________________________________________________________________
//////////////////////////////////////////////////////////////////////////////
//
// class AliHBTQInvCorrelFctn
// class AliHBTQOutLCMSCorrelFctn
// class AliHBTQLongLCMSCorrelFctn
// class AliHBTQSideLCMSCorrelFctn
// class AliHBTInvMassCorrelFctn
// class AliHBTTwoKStarCorrelFctn
//
// Set of functions:
//   Q Invaraint Correlation Function
//   Invariant Mass Function
//
// more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
// Piotr.Skowronski@cern.ch
//
//////////////////////////////////////////////////////////////////////////////

ClassImp(AliHBTQInvCorrelFctn)

//Corroleation function is created from dividing two histograms of QInvariant:
//  of particles from the same evnt
//by 
//  of particles from different events

AliHBTQInvCorrelFctn::AliHBTQInvCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("qinvcf","Q_{inv} Correlation Function");
}
/*************************************************************************************/ 

TH1* AliHBTQInvCorrelFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTOutSideLongFctn)

AliHBTOutSideLongFctn::AliHBTOutSideLongFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
                                                   Int_t nYbins, Double_t maxYval, Double_t minYval,
                                                   Int_t nZbins, Double_t maxZval, Double_t minZval):
 AliHBTOnePairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval),
 fAbs(kTRUE)
{
//ctor
  fWriteNumAndDen = kTRUE;//change default behaviour
  Rename("qoslcf","Q_{out}-Q_{side}-Q_{long} Correlation Fctn");
}
/*************************************************************************************/ 

TH1* AliHBTOutSideLongFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}

void AliHBTOutSideLongFctn::GetValues(AliHBTPair* pair, Double_t& x, Double_t& y, Double_t& z) const
{ 
  //calculates values of that function
  //qout qside and qlong
  
  x=pair->GetQOutLCMS(); 
  y=pair->GetQSideLCMS(); 
  z=pair->GetQLongLCMS();
  if (fAbs)
   {
     x = TMath::Abs(x);
     y = TMath::Abs(y);
     z = TMath::Abs(z);
   }
} 

/*************************************************************************************/ 

ClassImp(AliHBTQOutLCMSCorrelFctn)
    
AliHBTQOutLCMSCorrelFctn::AliHBTQOutLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("qoutcf","Q_{out} Correlation Function");
}
/*************************************************************************************/ 
    
TH1* AliHBTQOutLCMSCorrelFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTQLongLCMSCorrelFctn)
    
AliHBTQLongLCMSCorrelFctn::AliHBTQLongLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("qlongcf","Q_{long} Correlation Function");
}
/*************************************************************************************/ 
    
TH1* AliHBTQLongLCMSCorrelFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTQSideLCMSCorrelFctn)
    
AliHBTQSideLCMSCorrelFctn::AliHBTQSideLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
 //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("qsidecf","Q_{side} Correlation Function");
}
/*************************************************************************************/ 
    
TH1* AliHBTQSideLCMSCorrelFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}


/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 
ClassImp(AliHBTQtLCMSCorrelFctn)
    
AliHBTQtLCMSCorrelFctn::AliHBTQtLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("Qtcf","Q_{t}(LCMS) Correlation Function");
}
/*************************************************************************************/ 
    
TH1* AliHBTQtLCMSCorrelFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 
ClassImp(AliHBTQtCorrelFctn)
    
AliHBTQtCorrelFctn::AliHBTQtCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("qtcf","Q_{t} Correlation Function");
}
/*************************************************************************************/ 
    
TH1* AliHBTQtCorrelFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTInvMassCorrelFctn)

AliHBTInvMassCorrelFctn::AliHBTInvMassCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
 //ctor 
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("InvMass CF","Invariant Mass Correlation Function");
}

TH1* AliHBTInvMassCorrelFctn::GetResult()
{
 //returns result
 return GetNumerator();
}
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTTwoKStarCorrelFctn)

AliHBTTwoKStarCorrelFctn::AliHBTTwoKStarCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
 //ctor 
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("twokstarcf","2K^{*} Correlation Function");
}

/*************************************************************************************/ 

TH1* AliHBTTwoKStarCorrelFctn::GetResult()
{  
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}

/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 
ClassImp(AliHBTAvSeparCorrelFctn)

AliHBTAvSeparCorrelFctn::AliHBTAvSeparCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
 AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
 //ctor 
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("avsepcf","Avarage separation Correlation Function");
}

/*************************************************************************************/ 

TH1* AliHBTAvSeparCorrelFctn::GetResult()
{  
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}

/*************************************************************************************/
Commit	Line	Data
1b446896	1	#include "AliHBTCorrelFctn.h"
db21c92b	2	//____________________________________________________________________________
db21c92b	3	//////////////////////////////////////////////////////////////////////////////
953577f3	4	//
db21c92b	5	// class AliHBTQInvCorrelFctn
78d7c6d3	6	// class AliHBTQOutLCMSCorrelFctn
	7	// class AliHBTQLongLCMSCorrelFctn
	8	// class AliHBTQSideLCMSCorrelFctn
db21c92b	9	// class AliHBTInvMassCorrelFctn
db21c92b	10	// class AliHBTTwoKStarCorrelFctn
953577f3	11	//
db21c92b	12	// Set of functions:
	13	// Q Invaraint Correlation Function
	14	// Invariant Mass Function
953577f3	15	//
db21c92b	16	// more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
	17	// Piotr.Skowronski@cern.ch
	18	//
	19	//////////////////////////////////////////////////////////////////////////////
1b446896	20
	21	ClassImp(AliHBTQInvCorrelFctn)
	22
88cb7938	23	//Corroleation function is created from dividing two histograms of QInvariant:
	24	// of particles from the same evnt
	25	//by
	26	// of particles from different events
1b446896	27
e6e21fa2	28	AliHBTQInvCorrelFctn::AliHBTQInvCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
e6e21fa2	29	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
953577f3	30	{
e6e21fa2	31	fWriteNumAndDen = kTRUE;//change default behaviour
953577f3	32	Rename("qinvcf","Q_{inv} Correlation Function");
953577f3	33	}
e6e21fa2	34	/*************************************************************************************/
953577f3	35
1b446896	36	TH1* AliHBTQInvCorrelFctn::GetResult()
fc13079c	37	{
	38	//returns the scaled ratio
	39	delete fRatio;
	40	fRatio = GetRatio(Scale());
	41	return fRatio;
953577f3	42	}
	43	/*************************************************************************************/
	44	/*************************************************************************************/
	45	/*************************************************************************************/
	46
3cc1fefe	47	ClassImp(AliHBTOutSideLongFctn)
	48
	49	AliHBTOutSideLongFctn::AliHBTOutSideLongFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
	50	Int_t nYbins, Double_t maxYval, Double_t minYval,
	51	Int_t nZbins, Double_t maxZval, Double_t minZval):
ed3b923a	52	AliHBTOnePairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval),
ed3b923a	53	fAbs(kTRUE)
3cc1fefe	54	{
f78d233e	55	//ctor
3cc1fefe	56	fWriteNumAndDen = kTRUE;//change default behaviour
	57	Rename("qoslcf","Q_{out}-Q_{side}-Q_{long} Correlation Fctn");
	58	}
f78d233e	59	/*************************************************************************************/
3cc1fefe	60
f78d233e	61	TH1* AliHBTOutSideLongFctn::GetResult()
	62	{
	63	//returns the scaled ratio
	64	delete fRatio;
	65	fRatio = GetRatio(Scale());
	66	return fRatio;
	67	}
ed3b923a	68
	69	void AliHBTOutSideLongFctn::GetValues(AliHBTPair* pair, Double_t& x, Double_t& y, Double_t& z) const
	70	{
	71	//calculates values of that function
	72	//qout qside and qlong
	73
78d7c6d3	74	x=pair->GetQOutLCMS();
	75	y=pair->GetQSideLCMS();
	76	z=pair->GetQLongLCMS();
ed3b923a	77	if (fAbs)
	78	{
	79	x = TMath::Abs(x);
	80	y = TMath::Abs(y);
	81	z = TMath::Abs(z);
	82	}
	83	}
	84
3cc1fefe	85	/*************************************************************************************/
3cc1fefe	86
78d7c6d3	87	ClassImp(AliHBTQOutLCMSCorrelFctn)
e6e21fa2	88
78d7c6d3	89	AliHBTQOutLCMSCorrelFctn::AliHBTQOutLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
e6e21fa2	90	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	91	{
	92	//ctor
	93	fWriteNumAndDen = kTRUE;//change default behaviour
	94	Rename("qoutcf","Q_{out} Correlation Function");
	95	}
	96	/*************************************************************************************/
	97
78d7c6d3	98	TH1* AliHBTQOutLCMSCorrelFctn::GetResult()
1b446896	99	{
fc13079c	100	//returns the scaled ratio
	101	delete fRatio;
	102	fRatio = GetRatio(Scale());
	103	return fRatio;
1b446896	104	}
953577f3	105	/*************************************************************************************/
	106	/*************************************************************************************/
	107	/*************************************************************************************/
	108
78d7c6d3	109	ClassImp(AliHBTQLongLCMSCorrelFctn)
e6e21fa2	110
78d7c6d3	111	AliHBTQLongLCMSCorrelFctn::AliHBTQLongLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
e6e21fa2	112	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	113	{
	114	//ctor
	115	fWriteNumAndDen = kTRUE;//change default behaviour
	116	Rename("qlongcf","Q_{long} Correlation Function");
	117	}
	118	/*************************************************************************************/
	119
78d7c6d3	120	TH1* AliHBTQLongLCMSCorrelFctn::GetResult()
953577f3	121	{
fc13079c	122	//returns the scaled ratio
	123	delete fRatio;
	124	fRatio = GetRatio(Scale());
	125	return fRatio;
953577f3	126	}
	127	/*************************************************************************************/
	128	/*************************************************************************************/
	129	/*************************************************************************************/
	130
78d7c6d3	131	ClassImp(AliHBTQSideLCMSCorrelFctn)
e6e21fa2	132
78d7c6d3	133	AliHBTQSideLCMSCorrelFctn::AliHBTQSideLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
e6e21fa2	134	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	135	{
	136	//ctor
	137	fWriteNumAndDen = kTRUE;//change default behaviour
	138	Rename("qsidecf","Q_{side} Correlation Function");
	139	}
	140	/*************************************************************************************/
	141
78d7c6d3	142	TH1* AliHBTQSideLCMSCorrelFctn::GetResult()
953577f3	143	{
fc13079c	144	//returns the scaled ratio
	145	delete fRatio;
	146	fRatio = GetRatio(Scale());
	147	return fRatio;
953577f3	148	}
953577f3	149
1b446896	150
98bc79ce	151	/*************************************************************************************/
	152	/*************************************************************************************/
	153	/*************************************************************************************/
	154	ClassImp(AliHBTQtLCMSCorrelFctn)
	155
	156	AliHBTQtLCMSCorrelFctn::AliHBTQtLCMSCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
	157	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	158	{
	159	//ctor
	160	fWriteNumAndDen = kTRUE;//change default behaviour
00aa96f1	161	Rename("Qtcf","Q_{t}(LCMS) Correlation Function");
98bc79ce	162	}
	163	/*************************************************************************************/
	164
	165	TH1* AliHBTQtLCMSCorrelFctn::GetResult()
	166	{
	167	//returns the scaled ratio
	168	delete fRatio;
	169	fRatio = GetRatio(Scale());
	170	return fRatio;
	171	}
953577f3	172	/*************************************************************************************/
	173	/*************************************************************************************/
	174	/*************************************************************************************/
00aa96f1	175	ClassImp(AliHBTQtCorrelFctn)
	176
	177	AliHBTQtCorrelFctn::AliHBTQtCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
	178	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	179	{
	180	//ctor
	181	fWriteNumAndDen = kTRUE;//change default behaviour
	182	Rename("qtcf","Q_{t} Correlation Function");
	183	}
	184	/*************************************************************************************/
	185
	186	TH1* AliHBTQtCorrelFctn::GetResult()
	187	{
	188	//returns the scaled ratio
	189	delete fRatio;
	190	fRatio = GetRatio(Scale());
	191	return fRatio;
	192	}
	193	/*************************************************************************************/
	194	/*************************************************************************************/
	195	/*************************************************************************************/
1b446896	196
	197	ClassImp(AliHBTInvMassCorrelFctn)
	198
e6e21fa2	199	AliHBTInvMassCorrelFctn::AliHBTInvMassCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
e6e21fa2	200	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
1b446896	201	{
e6e21fa2	202	//ctor
	203	fWriteNumAndDen = kTRUE;//change default behaviour
	204	Rename("InvMass CF","Invariant Mass Correlation Function");
1b446896	205	}
	206
	207	TH1* AliHBTInvMassCorrelFctn::GetResult()
	208	{
e6e21fa2	209	//returns result
fc13079c	210	return GetNumerator();
1b446896	211	}
4866a439	212	/*************************************************************************************/
	213	/*************************************************************************************/
	214	/*************************************************************************************/
	215
	216	ClassImp(AliHBTTwoKStarCorrelFctn)
	217
e6e21fa2	218	AliHBTTwoKStarCorrelFctn::AliHBTTwoKStarCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
e6e21fa2	219	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
4866a439	220	{
e6e21fa2	221	//ctor
e6e21fa2	222	fWriteNumAndDen = kTRUE;//change default behaviour
4866a439	223	Rename("twokstarcf","2K^{*} Correlation Function");
	224	}
	225
	226	/*************************************************************************************/
	227
	228	TH1* AliHBTTwoKStarCorrelFctn::GetResult()
	229	{
fc13079c	230	//returns the scaled ratio
	231	delete fRatio;
	232	fRatio = GetRatio(Scale());
	233	return fRatio;
4866a439	234	}
	235
	236	/*************************************************************************************/
087f87e7	237	/*************************************************************************************/
	238	/*************************************************************************************/
	239	ClassImp(AliHBTAvSeparCorrelFctn)
	240
	241	AliHBTAvSeparCorrelFctn::AliHBTAvSeparCorrelFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
	242	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	243	{
	244	//ctor
	245	fWriteNumAndDen = kTRUE;//change default behaviour
	246	Rename("avsepcf","Avarage separation Correlation Function");
	247	}
	248
	249	/*************************************************************************************/
	250
	251	TH1* AliHBTAvSeparCorrelFctn::GetResult()
	252	{
	253	//returns the scaled ratio
	254	delete fRatio;
	255	fRatio = GetRatio(Scale());
	256	return fRatio;
	257	}
	258
	259	/*************************************************************************************/