HBTAN/AliHBTLLWeightFctn.cxx

   1 //Author: Ludmila Malinina, JINR (malinina@sunhe.jinr.ru)
   2
   3 #include "AliHBTLLWeightFctn.h"
   4 #include "AliHBTLLWeights.h"
   5 #include "AliHBTLLWeightsPID.h"
   6
   7 //--for test--AliHBTLLWeightQInvFctn* yyy= new AliHBTLLWeightQInvFctn();
   8
   9 ClassImp( AliHBTLLWeightQInvFctn )
  10 /****************************************************************/
  11 AliHBTLLWeightQInvFctn::AliHBTLLWeightQInvFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
  12            AliHBTTwoPairFctn1D(nbins,maxXval,minXval)
  13 {
  14 //ctor
  15 }
  16 /****************************************************************/
  17 void  AliHBTLLWeightQInvFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
  18 {
  19 //Processes Particles and tracks Same different even
  20   trackpair = CheckPair(trackpair);
  21   partpair  = CheckPair(partpair);
  22   if ( trackpair && partpair)
  23   {
  24     Double_t weightPID=1.;
  25     Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
  26     Double_t weight=weightHBT*weightPID;
  27     if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQInv(),weight);
  28   }
  29 }
  30 /****************************************************************/
  31
  32 void  AliHBTLLWeightQInvFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
  33 {
  34   trackpair = CheckPair(trackpair);
  35   partpair  = CheckPair(partpair);
  36   if ( trackpair && partpair)
  37   {
  38      fDenominator->Fill(trackpair->GetQInv());
  39   }
  40 }
  41 /**************************************************************/
  42 TH1* AliHBTLLWeightQInvFctn::GetResult()
  43
  44 {
  45 //returns ratio of numerator and denominator
  46  return GetRatio(Scale());
  47 }
  48