partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t weight = partpair->GetWeight();
- if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQOutCMSLC(),weight);
+ if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQOutLCMS(),weight);
}
/**************************************************************/
partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t weight = partpair->GetWeight();
- if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQSideCMSLC(),weight);
+ if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQSideLCMS(),weight);
}
/**************************************************************/
partpair = CheckPair(partpair);
if (partpair == 0x0) return;
Double_t weight = partpair->GetWeight();
- if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQLongCMSLC(),weight);
+ if(TMath::Abs(weight)<=10.) fNumerator->Fill(partpair->GetQLongLCMS(),weight);
}
/**************************************************************/
/*************************************************************/
/*************************************************************/
+ClassImp(AliHBTWeightTheorQtFctn)
+/*************************************************************/
+
+AliHBTWeightTheorQtFctn::AliHBTWeightTheorQtFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
+ AliHBTOnePairFctn1D(nbins,maxXval,minXval)
+{
+ //ctor
+ fWriteNumAndDen = kTRUE;//change default behaviour
+ Rename("wqttheorcf","Q_{t} Weight Theoretical Correlation Function");
+}
+/****************************************************************/
+
+void AliHBTWeightTheorQtFctn::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->GetQt(),weight);
+}
+/**************************************************************/
+
+TH1* AliHBTWeightTheorQtFctn::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,
/*************************************************************/
+void AliHBTWeightTheorOSLFctn::GetValues(AliHBTPair* pair, Double_t& x, Double_t& y, Double_t& z) const
+{
+//returns values of a functions for a given pair of particles
+ x=pair->GetQOutLCMS();
+ y=pair->GetQSideLCMS();
+ z=pair->GetQLongLCMS();
+
+ if (fAbs)
+ {
+ x = TMath::Abs(x);
+ y = TMath::Abs(y);
+ z = TMath::Abs(z);
+ }
+
+}
+
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->GetQOutCMSLC());
- Double_t side = TMath::Abs(partpair->GetQSideCMSLC());
- Double_t lon = TMath::Abs(partpair->GetQLongCMSLC());
+ Double_t out, side, lon;
+
+ GetValues(partpair, out, side, lon);
/*
if (out < 0.01)