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

//Author: Ludmila Malinina, JINR (malinina@sunhe.jinr.ru)
#include "AliHBTLLWeightFctn.h"
#include "AliHBTLLWeights.h"
#include "AliHBTLLWeightsPID.h"

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

ClassImp( AliHBTLLWeightQInvFctn )  
/****************************************************************/
AliHBTLLWeightQInvFctn::AliHBTLLWeightQInvFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
           AliHBTTwoPairFctn1D(nbins,maxXval,minXval)
{
//ctor
}
/****************************************************************/
void  AliHBTLLWeightQInvFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//Processes Particles and tracks Same different even
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)     
  {
    Double_t weightPID=1.;
    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
    Double_t weight=weightHBT*weightPID;
    if(TMath::Abs(weight)<=10.)fNumerator->Fill(trackpair->GetQInv(),weight);
  }
} 
/****************************************************************/

void  AliHBTLLWeightQInvFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)  
  {
     fDenominator->Fill(trackpair->GetQInv());
  }
}
/**************************************************************/
TH1* AliHBTLLWeightQInvFctn::GetResult() 
                                                                               
{ 
//returns ratio of numerator and denominator                                    
 return GetRatio(Scale());                                                  
}                    
                                                              
/**************************************************************************************/
/**************************************************************************************/
/**************************************************************************************/
/**************************************************************************************/

ClassImp(AliHBTLLWeightQOutFctn)
/****************************************************************/
void AliHBTLLWeightQOutFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)     
  {
    Double_t weightPID=1.;
    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
    Double_t weight=weightHBT*weightPID;
    if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQOutCMSLC(),weight);
  }
} 
/****************************************************************/

void AliHBTLLWeightQOutFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)  
  {
     fDenominator->Fill(trackpair->GetQOutCMSLC());
  }
}
/**************************************************************/
TH1* AliHBTLLWeightQOutFctn::GetResult() 
                                                                               
{ 
//returns ratio of numerator and denominator                                    
 return GetRatio(Scale());                                                  
}                    
                                                              
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTLLWeightQLongFctn)
/****************************************************************/
void  AliHBTLLWeightQLongFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//Processes Particles and tracks Same different even
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)     
  {
    Double_t weightPID=1.;
    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
    Double_t weight=weightHBT*weightPID;
    if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQLongCMSLC(),weight);
  }
} 
/****************************************************************/

void  AliHBTLLWeightQLongFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)  
  {
     fDenominator->Fill(trackpair->GetQLongCMSLC());
  }
}
/**************************************************************/
TH1* AliHBTLLWeightQLongFctn::GetResult() 
                                                                               
{ 
//returns ratio of numerator and denominator                                    
 return GetRatio(Scale());                                                  
}                    
                                                              
/*************************************************************************************/ 
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTLLWeightQSideFctn)
/****************************************************************/
void  AliHBTLLWeightQSideFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//Processes Particles and tracks Same different even
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)     
  {
    Double_t weightPID=1.;
    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
    Double_t weight=weightHBT*weightPID;
    if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQSideCMSLC(),weight);
  }
} 
/****************************************************************/

void  AliHBTLLWeightQSideFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)  
  {
     fDenominator->Fill(trackpair->GetQSideCMSLC());
  }
}
/**************************************************************/
TH1* AliHBTLLWeightQSideFctn::GetResult() 
                                                                               
{ 
//returns ratio of numerator and denominator                                    
 return GetRatio(Scale());                                                  
}                    
                                                              
/*************************************************************************************/ 
/*************************************************************************************/ 

ClassImp(AliHBTLLWeightTwoKStarFctn)
/****************************************************************/
void  AliHBTLLWeightTwoKStarFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//Processes Particles and tracks Same different even
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)     
  {
    Double_t weightPID=1.;
    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
    Double_t weight=weightHBT*weightPID;
    if(TMath::Abs(weight)<=10.) fNumerator->Fill(2.0*(trackpair->GetKStar()),weight);
  }
} 
/****************************************************************/

void  AliHBTLLWeightTwoKStarFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
  trackpair = CheckPair(trackpair);
  partpair  = CheckPair(partpair);
  if ( trackpair && partpair)  
  {
     fDenominator->Fill(2.0*(trackpair->GetKStar()));
  }
}
/**************************************************************/
TH1* AliHBTLLWeightTwoKStarFctn::GetResult() 
                                                                               
{ 
//returns ratio of numerator and denominator                                    
 return GetRatio(Scale());                                                  
}                    
                                                              
/*************************************************************************************/
Commit	Line	Data
3506b383	1	//Author: Ludmila Malinina, JINR (malinina@sunhe.jinr.ru)
7f92929e	2	#include "AliHBTLLWeightFctn.h"
7f92929e	3	#include "AliHBTLLWeights.h"
b3fb9814	4	#include "AliHBTLLWeightsPID.h"
7f92929e	5
	6	//--for test--AliHBTLLWeightQInvFctn* yyy= new AliHBTLLWeightQInvFctn();
	7
	8	ClassImp( AliHBTLLWeightQInvFctn )
	9	/****************************************************************/
	10	AliHBTLLWeightQInvFctn::AliHBTLLWeightQInvFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
	11	AliHBTTwoPairFctn1D(nbins,maxXval,minXval)
	12	{
2f8eea63	13	//ctor
7f92929e	14	}
	15	/****************************************************************/
	16	void AliHBTLLWeightQInvFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
	17	{
2f8eea63	18	//Processes Particles and tracks Same different even
7f92929e	19	trackpair = CheckPair(trackpair);
	20	partpair = CheckPair(partpair);
	21	if ( trackpair && partpair)
	22	{
2f8eea63	23	Double_t weightPID=1.;
	24	Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
	25	Double_t weight=weightHBT*weightPID;
b3fb9814	26	if(TMath::Abs(weight)<=10.)fNumerator->Fill(trackpair->GetQInv(),weight);
7f92929e	27	}
7f92929e	28	}
7f92929e	29	/****************************************************************/
2f8eea63	30
2f8eea63	31	void AliHBTLLWeightQInvFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
7f92929e	32	{
	33	trackpair = CheckPair(trackpair);
	34	partpair = CheckPair(partpair);
	35	if ( trackpair && partpair)
	36	{
	37	fDenominator->Fill(trackpair->GetQInv());
	38	}
	39	}
	40	/**************************************************************/
	41	TH1* AliHBTLLWeightQInvFctn::GetResult()
	42
	43	{
	44	//returns ratio of numerator and denominator
af8cf344	45	return GetRatio(Scale());
7f92929e	46	}
7f92929e	47
b3fb9814	48	/**************************************************************************************/
	49	/**************************************************************************************/
	50	/**************************************************************************************/
	51	/**************************************************************************************/
	52
	53	ClassImp(AliHBTLLWeightQOutFctn)
	54	/****************************************************************/
	55	void AliHBTLLWeightQOutFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
	56	{
	57	trackpair = CheckPair(trackpair);
	58	partpair = CheckPair(partpair);
	59	if ( trackpair && partpair)
	60	{
	61	Double_t weightPID=1.;
	62	Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
	63	Double_t weight=weightHBT*weightPID;
	64	if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQOutCMSLC(),weight);
	65	}
	66	}
	67	/****************************************************************/
	68
	69	void AliHBTLLWeightQOutFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
	70	{
	71	trackpair = CheckPair(trackpair);
	72	partpair = CheckPair(partpair);
	73	if ( trackpair && partpair)
	74	{
	75	fDenominator->Fill(trackpair->GetQOutCMSLC());
	76	}
	77	}
	78	/**************************************************************/
	79	TH1* AliHBTLLWeightQOutFctn::GetResult()
	80
	81	{
	82	//returns ratio of numerator and denominator
	83	return GetRatio(Scale());
	84	}
	85
	86	/*************************************************************************************/
	87	/*************************************************************************************/
	88	/*************************************************************************************/
	89
	90	ClassImp(AliHBTLLWeightQLongFctn)
	91	/****************************************************************/
	92	void AliHBTLLWeightQLongFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
	93	{
	94	//Processes Particles and tracks Same different even
	95	trackpair = CheckPair(trackpair);
	96	partpair = CheckPair(partpair);
	97	if ( trackpair && partpair)
	98	{
	99	Double_t weightPID=1.;
	100	Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
	101	Double_t weight=weightHBT*weightPID;
	102	if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQLongCMSLC(),weight);
	103	}
	104	}
	105	/****************************************************************/
	106
	107	void AliHBTLLWeightQLongFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
	108	{
	109	trackpair = CheckPair(trackpair);
	110	partpair = CheckPair(partpair);
	111	if ( trackpair && partpair)
112	{
113	fDenominator->Fill(trackpair->GetQLongCMSLC());
114	}
115	}
116	/**************************************************************/
117	TH1* AliHBTLLWeightQLongFctn::GetResult()
118
119	{
120	//returns ratio of numerator and denominator
121	return GetRatio(Scale());
122	}
123
124	/*************************************************************************************/
125	/*************************************************************************************/
126	/*************************************************************************************/
127
128	ClassImp(AliHBTLLWeightQSideFctn)
129	/****************************************************************/
130	void AliHBTLLWeightQSideFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
131	{
132	//Processes Particles and tracks Same different even
133	trackpair = CheckPair(trackpair);
134	partpair = CheckPair(partpair);
135	if ( trackpair && partpair)
136	{
137	Double_t weightPID=1.;
138	Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
139	Double_t weight=weightHBT*weightPID;
140	if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQSideCMSLC(),weight);
141	}
142	}
143	/****************************************************************/
144
145	void AliHBTLLWeightQSideFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
146	{
147	trackpair = CheckPair(trackpair);
148	partpair = CheckPair(partpair);
149	if ( trackpair && partpair)
150	{
151	fDenominator->Fill(trackpair->GetQSideCMSLC());
152	}
153	}
154	/**************************************************************/
155	TH1* AliHBTLLWeightQSideFctn::GetResult()
156
157	{
158	//returns ratio of numerator and denominator
159	return GetRatio(Scale());
160	}
161
162	/*************************************************************************************/
163	/*************************************************************************************/
164
165	ClassImp(AliHBTLLWeightTwoKStarFctn)
166	/****************************************************************/
167	void AliHBTLLWeightTwoKStarFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
168	{
169	//Processes Particles and tracks Same different even
170	trackpair = CheckPair(trackpair);
171	partpair = CheckPair(partpair);
172	if ( trackpair && partpair)
173	{
174	Double_t weightPID=1.;
175	Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
176	Double_t weight=weightHBT*weightPID;
177	if(TMath::Abs(weight)<=10.) fNumerator->Fill(2.0*(trackpair->GetKStar()),weight);
178	}
179	}
180	/****************************************************************/
181
182	void AliHBTLLWeightTwoKStarFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
183	{
184	trackpair = CheckPair(trackpair);
185	partpair = CheckPair(partpair);
186	if ( trackpair && partpair)
187	{
188	fDenominator->Fill(2.0*(trackpair->GetKStar()));
189	}
190	}
191	/**************************************************************/
192	TH1* AliHBTLLWeightTwoKStarFctn::GetResult()
193
194	{
195	//returns ratio of numerator and denominator
196	return GetRatio(Scale());
197	}
198
199	/*************************************************************************************/
200