HBTAN/AliHBTMonSplitting.cxx

   1 #include "AliHBTMonSplitting.h"
   2
   3 ClassImp(AliHBTMonSplittingQosl)
   4
   5 AliHBTMonSplittingQosl::AliHBTMonSplittingQosl(Int_t nXbins, Double_t maxXval, Double_t minXval,
   6                                                Int_t nYbins, Double_t maxYval, Double_t minYval,
   7                                                 Int_t nZbins, Double_t maxZval, Double_t minZval):
   8  AliHBTTwoPairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
   9 {
  10   //ctor
  11   Rename("splitosl","Q_{out}-Q_{side}-Q_{long} Splitting Monitoring Function");
  12 }
  13
  14 void   AliHBTMonSplittingQosl::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
  15 {
  16   AliVAODParticle* p1 = partpair->Particle1();
  17   AliVAODParticle* p2 = partpair->Particle2();
  18
  19   if (p1->Px() != p2->Px()) return;
  20   if (p1->Py() != p2->Py()) return;
  21   if (p1->Pz() != p2->Pz()) return;
  22
  23   Double_t out = trackpair->GetQOutLCMS();
  24   Double_t side = trackpair->GetQSideLCMS();
  25   Double_t lon = trackpair->GetQLongLCMS();
  26
  27   fNumerator->Fill(out,side,lon);//here we fill in q's corresponding to track pair
  28                                           //weight calculated for the simulated one
  29 }
  30
  31 ClassImp(AliHBTMonSplittingDptDthetaDphi)
  32
  33 AliHBTMonSplittingDptDthetaDphi::AliHBTMonSplittingDptDthetaDphi(Int_t nXbins, Double_t maxXval, Double_t minXval,
  34                                                Int_t nYbins, Double_t maxYval, Double_t minYval,
  35                                                 Int_t nZbins, Double_t maxZval, Double_t minZval):
  36  AliHBTTwoPairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
  37 {
  38   //ctor
  39   Rename("splitdpdthedphi","\\Deltap_{t}-\\Delta\\theta-\\Delta\\phi Splitting Monitoring Function");
  40 }
  41
  42 void   AliHBTMonSplittingDptDthetaDphi::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
  43 {
  44   AliVAODParticle* p1 = partpair->Particle1();
  45   AliVAODParticle* p2 = partpair->Particle2();
  46
  47   if (p1->Px() != p2->Px()) return;
  48   if (p1->Py() != p2->Py()) return;
  49   if (p1->Pz() != p2->Pz()) return;
  50
  51   Double_t dpt = trackpair->GetDeltaPt();
  52   Double_t dphi = trackpair->GetDeltaPhi();
  53   Double_t dtheta = trackpair->GetDeltaTheta();
  54
  55   fNumerator->Fill(dpt, dphi, dtheta);//here we fill in q's corresponding to track pair
  56                                           //weight calculated for the simulated one
  57 }