]>
Commit | Line | Data |
---|---|---|
a296aa84 | 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 | } |