void AliHBTCorrFitFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//Fills the numerator using pair from the same event
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if(partpair == 0x0) return;
trackpair = CheckPair(trackpair);
if(trackpair == 0x0) return;
{
// Fills the denominator using mixed pairs
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQInv());
void AliHBTTwoPairFctn1D::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
// Fills the numerator using pairs from the same event
- partpair = CheckPair(partpair);
- if( partpair == 0x0) return;
+ trackpair = CheckPair(trackpair);
+ if( trackpair == 0x0) return;
Double_t x = GetValue(trackpair,partpair);
fNumerator->Fill(x);
void AliHBTTwoPairFctn1D::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
// Fills the denumerator usin mixed pairs
- partpair = CheckPair(partpair);
- if( partpair == 0x0) return;
+ trackpair = CheckPair(trackpair);
+ if( trackpair == 0x0) return;
Double_t x = GetValue(trackpair,partpair);
fDenominator->Fill(x);
void AliHBTTwoPairFctn2D::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//processes pair of particles coming from a same events (real pair)
- partpair = CheckPair(partpair); //check cuts
- if (partpair == 0x0) return;
+ trackpair = CheckPair(trackpair);
+ if( trackpair == 0x0) return;
Double_t x,y;
GetValues(trackpair,partpair,x,y);
void AliHBTTwoPairFctn2D::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//processes pair of particles coming from a different events (mixed pair)
- partpair = CheckPair(partpair);
- if (partpair == 0x0) return;
+ trackpair = CheckPair(trackpair);
+ if( trackpair == 0x0) return;
Double_t x,y;
GetValues(trackpair,partpair,x,y);
void AliHBTTwoPairFctn3D::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
// Fills th numerator using pairs from the same event
- partpair = CheckPair(partpair);
- if( partpair == 0x0) return;
+ trackpair = CheckPair(trackpair);
+ if( trackpair == 0x0) return;
Double_t x,y,z;
GetValues(trackpair,partpair,x,y,z);
void AliHBTTwoPairFctn3D::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
// Fills the denumerator using mixed pairs
- partpair = CheckPair(partpair);
- if( partpair == 0x0) return;
+ trackpair = CheckPair(trackpair);
+ if( trackpair == 0x0) return;
Double_t x,y,z;
GetValues(trackpair,partpair,x,y,z);
//check if pair and both particles meets the cut criteria
if(fPairCut->Rejected(pair)) return 0x0; //if the pair is BAD
-// It is notallowed to change the order here beacause analysis enforce the order
+// It is not allowed to change the order here beacause analysis enforce the order
// {//it is BAD
// pair = pair->GetSwappedPair();
// //
///////////////////////////////////////////////////////
+AliHBTNonId3DCorrFctn x;
+
ClassImp(AliHBTNonId3DCorrFctn)
AliHBTNonId3DCorrFctn::AliHBTNonId3DCorrFctn(const char* name, const char* title):
void WriteFunction();
- TH1* GetResult();
+ TH1* GetResult();
protected:
- Double_t GetValue(AliHBTPair* pair) {return pair->GetQInv();}
+ Double_t GetValue(AliHBTPair* pair) const {return pair->GetQInv();}
void BuildHistos(Int_t nbins, Float_t max, Float_t min);
TH1D* fNumOutP;
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
fNumerator->Fill(trackpair->GetQInv());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
fDenominator->Fill(trackpair->GetQInv());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
fNumerator->Fill(trackpair->GetQInv(),partpair->GetWeight());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
fDenominator->Fill(trackpair->GetQInv());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t weight = partpair->GetWeight();
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t out = TMath::Abs(trackpair->GetQOutLCMS());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t out = TMath::Abs(trackpair->GetQOutLCMS());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t out = TMath::Abs(trackpair->GetQOutLCMS());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t x = TMath::Abs(trackpair->GetDeltaPt());
trackpair = CheckPair(trackpair);
if (trackpair == 0x0) return;
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t x = TMath::Abs(trackpair->GetDeltaPt());
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
// Fills the denominator using mixed pairs
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQInv());
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQOutLCMS());
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQLongLCMS());
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQSideLCMS());
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(2.0*(trackpair->GetKStar()));
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQOutLCMS(),trackpair->GetQSideLCMS());
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQOutLCMS(),trackpair->GetQLongLCMS());
{
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
fDenominator->Fill(trackpair->GetQSideLCMS(),trackpair->GetQLongLCMS());
//Fills numerator
//process particles from same events (fills numerator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
// Double_t weightPID=1.;
Double_t side = TMath::Abs(trackpair->GetQSideLCMS());
Double_t lon = TMath::Abs(trackpair->GetQLongLCMS());
- if (out < 0.005)
- if (side < 0.005)
- if (lon < 0.005)
- {
- trackpair->Particle1()->Print();
- trackpair->Particle2()->Print();
- Info("","Delta Theta %f, Delta Phi %f",
- trackpair->GetDeltaTheta(),trackpair->GetDeltaPhi());
- }
-
fNumerator->Fill(out,side,lon,weight);//here we fill in q's corresponding to track pair
//weight calculated for the simulated one
}
{
//process particles from diff events (fills denominator)
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
if ( trackpair && partpair)
{
Double_t out = TMath::Abs(trackpair->GetQOutLCMS());
void AliHBTWeightNonId3DCorrFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
{
//Fills the numerator using pair from the same event
- partpair = CheckPair(partpair);
+
+// partpair = CheckPair(partpair);//why checking cuts on simulated?
+
trackpair = CheckPair(trackpair);
if(partpair == 0x0) return;
if(trackpair == 0x0) return;
{
trackpair = CheckPair(trackpair);
- partpair = CheckPair(partpair);
+// partpair = CheckPair(partpair);
double tKStar = TMath::Abs(partpair->GetKStar());
double tKStarOut = trackpair->GetKStarOut();