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

/* $Id$ */

//_____________________________________________________________________________
///////////////////////////////////////////////////////////////////////////////
//
// class AliHBTWeightTheorQInvFctn
//
// This function allows to obtain Q_inv correlation function with weights
// calculated by Lednicky's alghorithm.
// Numerator is filled with weighted events. Weights are attributed to simulated particles.
// Weights are calculated with corresponding simulated particles momenta.
// Denominator is filled with mixing unweighted simulated particles.
// One needs only simulated pairs, so 
// this function is of class AliHBTOnePairFctn1D.
//-----------------------------------------------------------
// This class introduces the weights calculated according 
// with functions of efficiency of identification (TPC+TOF) 
// (calculated by B.V. Batyunia).
// Author: Ludmila Malinina, JINR (malinina@sunhe.jinr.ru)
//-----------------------------------------------------------

#include "AliHBTWeightTheorFctn.h"

//--for test--AliHBTWeightQInvFctn* yyy= new AliHBTWeightQInvFctn();
/*************************************************************/
/*************************************************************/
/*************************************************************/

ClassImp(AliHBTWeightTheorQInvFctn)  
/*************************************************************/

AliHBTWeightTheorQInvFctn::AliHBTWeightTheorQInvFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
  AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("wqinvtheorcf","Q_{inv} Weight Theoretical Correlation Function");
}
/**************************************************************/

void  AliHBTWeightTheorQInvFctn::ProcessSameEventParticles(AliHBTPair* partpair)
{
  //Processes Particles and tracks Same different event
  partpair  = CheckPair(partpair);
  if (partpair == 0x0) return;
  Double_t weight = partpair->GetWeight();
  if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQInv(),weight);
} 
/**************************************************************/

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

ClassImp(AliHBTWeightTheorQOutFctn)  
/*************************************************************/

AliHBTWeightTheorQOutFctn::AliHBTWeightTheorQOutFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
  AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("wqouttheorcf","Q_{out} Weight Theoretical Correlation Function");
}
/****************************************************************/

void  AliHBTWeightTheorQOutFctn::ProcessSameEventParticles(AliHBTPair* partpair)
{
  //Processes Particles and tracks Same different even
  partpair  = CheckPair(partpair);
  if (partpair == 0x0) return;
  Double_t weight = partpair->GetWeight();
  if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQOutLCMS(),weight);
} 

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

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

ClassImp(AliHBTWeightTheorQSideFctn)  
/*************************************************************/

AliHBTWeightTheorQSideFctn::AliHBTWeightTheorQSideFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
  AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("wqsidetheorcf","Q_{side} Weight Theoretical Correlation Function");
}
/****************************************************************/

void  AliHBTWeightTheorQSideFctn::ProcessSameEventParticles(AliHBTPair* partpair)
{
  //Processes Particles and tracks Same different even
  partpair  = CheckPair(partpair);
  if (partpair == 0x0) return;
  Double_t weight = partpair->GetWeight();
  if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQSideLCMS(),weight);
} 
/**************************************************************/

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

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

ClassImp(AliHBTWeightTheorQLongFctn)  
/*************************************************************/

AliHBTWeightTheorQLongFctn::AliHBTWeightTheorQLongFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
  AliHBTOnePairFctn1D(nbins,maxXval,minXval)
{
  //ctor
 fWriteNumAndDen = kTRUE;//change default behaviour
 Rename("wqlongtheorcf","Q_{long} Weight Theoretical Correlation Function");
}
/****************************************************************/

void  AliHBTWeightTheorQLongFctn::ProcessSameEventParticles(AliHBTPair* partpair)
{
  //Processes Particles and tracks Same different even
  partpair  = CheckPair(partpair);
  if (partpair == 0x0) return;
  Double_t weight = partpair->GetWeight();
  if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQLongLCMS(),weight);
} 
/**************************************************************/

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

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

ClassImp(AliHBTWeightTheorOSLFctn)

AliHBTWeightTheorOSLFctn::AliHBTWeightTheorOSLFctn(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)
{
  //ctor
  fWriteNumAndDen = kTRUE;//change default behaviour
  Rename("wqosltheorcf","Q_{out}-Q_{side}-Q_{long} Weight Theoretical Correlation Fctn");
}

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

void AliHBTWeightTheorOSLFctn::ProcessSameEventParticles(AliHBTPair* partpair)
{
//Fills numerator
  partpair  = CheckPair(partpair);
  if (partpair == 0x0) return;
  Double_t weight = partpair->GetWeight();
  Double_t out = TMath::Abs(partpair->GetQOutLCMS());
  Double_t side = TMath::Abs(partpair->GetQSideLCMS());
  Double_t lon = TMath::Abs(partpair->GetQLongLCMS());

/*  
  if (out < 0.01)
    if (side < 0.01)
      if (lon < 0.01)
       {
         Info("TheorOSL Num","================================================================");
         Info("TheorOSL Num","o:%f, s:%f, l:%f, w%f",out,side,lon,weight);
         Info("TheorOSL Num","First");
         partpair->Particle1()->Print();
         Info("TheorOSL Num","Second");
         partpair->Particle2()->Print();
         fflush(0);
       }
*/ 
  fNumerator->Fill(out,side,lon,weight);
}
/*************************************************************/

TH1* AliHBTWeightTheorOSLFctn::GetResult()
{
 //returns the scaled ratio
 delete fRatio;
 fRatio = GetRatio(Scale());
 return fRatio;
}
Commit	Line	Data
604034d8	1	/* $Id$ */
604034d8	2
eb828ce0	3	//_____________________________________________________________________________
	4	///////////////////////////////////////////////////////////////////////////////
	5	//
	6	// class AliHBTWeightTheorQInvFctn
	7	//
	8	// This function allows to obtain Q_inv correlation function with weights
	9	// calculated by Lednicky's alghorithm.
	10	// Numerator is filled with weighted events. Weights are attributed to simulated particles.
	11	// Weights are calculated with corresponding simulated particles momenta.
	12	// Denominator is filled with mixing unweighted simulated particles.
	13	// One needs only simulated pairs, so
	14	// this function is of class AliHBTOnePairFctn1D.
604034d8	15	//-----------------------------------------------------------
	16	// This class introduces the weights calculated according
	17	// with functions of efficiency of identification (TPC+TOF)
	18	// (calculated by B.V. Batyunia).
	19	// Author: Ludmila Malinina, JINR (malinina@sunhe.jinr.ru)
	20	//-----------------------------------------------------------
a6e49985	21
dd82cadc	22	#include "AliHBTWeightTheorFctn.h"
a6e49985	23
dd82cadc	24	//--for test--AliHBTWeightQInvFctn* yyy= new AliHBTWeightQInvFctn();
9714a029	25	/*************************************************************/
	26	/*************************************************************/
	27	/*************************************************************/
a6e49985	28
dd82cadc	29	ClassImp(AliHBTWeightTheorQInvFctn)
a6e49985	30	/*************************************************************/
a6e49985	31
dd82cadc	32	AliHBTWeightTheorQInvFctn::AliHBTWeightTheorQInvFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
604034d8	33	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
a6e49985	34	{
604034d8	35	//ctor
9714a029	36	fWriteNumAndDen = kTRUE;//change default behaviour
9714a029	37	Rename("wqinvtheorcf","Q_{inv} Weight Theoretical Correlation Function");
a6e49985	38	}
fc13079c	39	/**************************************************************/
fc13079c	40
dd82cadc	41	void AliHBTWeightTheorQInvFctn::ProcessSameEventParticles(AliHBTPair* partpair)
a6e49985	42	{
dd82cadc	43	//Processes Particles and tracks Same different event
a6e49985	44	partpair = CheckPair(partpair);
ce5d39b1	45	if (partpair == 0x0) return;
dd82cadc	46	Double_t weight = partpair->GetWeight();
a6e49985	47	if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQInv(),weight);
a6e49985	48	}
a6e49985	49	/**************************************************************/
fc13079c	50
dd82cadc	51	TH1* AliHBTWeightTheorQInvFctn::GetResult()
a6e49985	52	{
fc13079c	53	//returns the scaled ratio
	54	delete fRatio;
	55	fRatio = GetRatio(Scale());
	56	return fRatio;
a6e49985	57	}
a6e49985	58
9714a029	59	/*************************************************************/
	60	/*************************************************************/
	61	/*************************************************************/
	62
dd82cadc	63	ClassImp(AliHBTWeightTheorQOutFctn)
9714a029	64	/*************************************************************/
9714a029	65
dd82cadc	66	AliHBTWeightTheorQOutFctn::AliHBTWeightTheorQOutFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
9714a029	67	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	68	{
	69	//ctor
	70	fWriteNumAndDen = kTRUE;//change default behaviour
	71	Rename("wqouttheorcf","Q_{out} Weight Theoretical Correlation Function");
	72	}
	73	/****************************************************************/
fc13079c	74
dd82cadc	75	void AliHBTWeightTheorQOutFctn::ProcessSameEventParticles(AliHBTPair* partpair)
9714a029	76	{
	77	//Processes Particles and tracks Same different even
	78	partpair = CheckPair(partpair);
ce5d39b1	79	if (partpair == 0x0) return;
dd82cadc	80	Double_t weight = partpair->GetWeight();
78d7c6d3	81	if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQOutLCMS(),weight);
9714a029	82	}
	83
	84	/**************************************************************/
dd82cadc	85	TH1* AliHBTWeightTheorQOutFctn::GetResult()
9714a029	86	{
fc13079c	87	//returns the scaled ratio
	88	delete fRatio;
	89	fRatio = GetRatio(Scale());
	90	return fRatio;
9714a029	91	}
	92
	93	/*************************************************************/
	94	/*************************************************************/
	95	/*************************************************************/
	96
dd82cadc	97	ClassImp(AliHBTWeightTheorQSideFctn)
9714a029	98	/*************************************************************/
9714a029	99
dd82cadc	100	AliHBTWeightTheorQSideFctn::AliHBTWeightTheorQSideFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
9714a029	101	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	102	{
	103	//ctor
	104	fWriteNumAndDen = kTRUE;//change default behaviour
	105	Rename("wqsidetheorcf","Q_{side} Weight Theoretical Correlation Function");
	106	}
	107	/****************************************************************/
fc13079c	108
dd82cadc	109	void AliHBTWeightTheorQSideFctn::ProcessSameEventParticles(AliHBTPair* partpair)
9714a029	110	{
	111	//Processes Particles and tracks Same different even
	112	partpair = CheckPair(partpair);
ce5d39b1	113	if (partpair == 0x0) return;
dd82cadc	114	Double_t weight = partpair->GetWeight();
78d7c6d3	115	if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQSideLCMS(),weight);
9714a029	116	}
9714a029	117	/**************************************************************/
fc13079c	118
dd82cadc	119	TH1* AliHBTWeightTheorQSideFctn::GetResult()
9714a029	120	{
fc13079c	121	//returns the scaled ratio
	122	delete fRatio;
	123	fRatio = GetRatio(Scale());
	124	return fRatio;
9714a029	125	}
	126
	127	/*************************************************************/
	128	/*************************************************************/
	129	/*************************************************************/
	130
dd82cadc	131	ClassImp(AliHBTWeightTheorQLongFctn)
9714a029	132	/*************************************************************/
9714a029	133
dd82cadc	134	AliHBTWeightTheorQLongFctn::AliHBTWeightTheorQLongFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
9714a029	135	AliHBTOnePairFctn1D(nbins,maxXval,minXval)
	136	{
	137	//ctor
	138	fWriteNumAndDen = kTRUE;//change default behaviour
	139	Rename("wqlongtheorcf","Q_{long} Weight Theoretical Correlation Function");
	140	}
	141	/****************************************************************/
fc13079c	142
dd82cadc	143	void AliHBTWeightTheorQLongFctn::ProcessSameEventParticles(AliHBTPair* partpair)
9714a029	144	{
	145	//Processes Particles and tracks Same different even
	146	partpair = CheckPair(partpair);
ce5d39b1	147	if (partpair == 0x0) return;
dd82cadc	148	Double_t weight = partpair->GetWeight();
78d7c6d3	149	if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQLongLCMS(),weight);
9714a029	150	}
9714a029	151	/**************************************************************/
fc13079c	152
dd82cadc	153	TH1* AliHBTWeightTheorQLongFctn::GetResult()
dd82cadc	154	{
fc13079c	155	//returns the scaled ratio
	156	delete fRatio;
	157	fRatio = GetRatio(Scale());
	158	return fRatio;
dd82cadc	159	}
	160
	161	/*************************************************************/
	162	/*************************************************************/
	163	/*************************************************************/
	164
	165	ClassImp(AliHBTWeightTheorOSLFctn)
	166
	167	AliHBTWeightTheorOSLFctn::AliHBTWeightTheorOSLFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
	168	Int_t nYbins, Double_t maxYval, Double_t minYval,
	169	Int_t nZbins, Double_t maxZval, Double_t minZval):
	170	AliHBTOnePairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
	171	{
f78d233e	172	//ctor
dd82cadc	173	fWriteNumAndDen = kTRUE;//change default behaviour
	174	Rename("wqosltheorcf","Q_{out}-Q_{side}-Q_{long} Weight Theoretical Correlation Fctn");
	175	}
f78d233e	176
dd82cadc	177	/*************************************************************/
	178
	179	void AliHBTWeightTheorOSLFctn::ProcessSameEventParticles(AliHBTPair* partpair)
	180	{
	181	//Fills numerator
	182	partpair = CheckPair(partpair);
	183	if (partpair == 0x0) return;
	184	Double_t weight = partpair->GetWeight();
78d7c6d3	185	Double_t out = TMath::Abs(partpair->GetQOutLCMS());
	186	Double_t side = TMath::Abs(partpair->GetQSideLCMS());
	187	Double_t lon = TMath::Abs(partpair->GetQLongLCMS());
92d472ba	188
	189	/*
	190	if (out < 0.01)
	191	if (side < 0.01)
	192	if (lon < 0.01)
	193	{
	194	Info("TheorOSL Num","================================================================");
	195	Info("TheorOSL Num","o:%f, s:%f, l:%f, w%f",out,side,lon,weight);
	196	Info("TheorOSL Num","First");
	197	partpair->Particle1()->Print();
	198	Info("TheorOSL Num","Second");
	199	partpair->Particle2()->Print();
	200	fflush(0);
	201	}
	202	*/
dd82cadc	203	fNumerator->Fill(out,side,lon,weight);
	204	}
	205	/*************************************************************/
	206
	207	TH1* AliHBTWeightTheorOSLFctn::GetResult()
9714a029	208	{
fc13079c	209	//returns the scaled ratio
	210	delete fRatio;
	211	fRatio = GetRatio(Scale());
	212	return fRatio;
9714a029	213	}