|
22 | 22 | #include <TDirectory.h> |
23 | 23 | #include <THn.h> |
24 | 24 | #include <TFile.h> |
| 25 | +#include <TVector2.h> |
25 | 26 |
|
26 | 27 | #include "Framework/runDataProcessing.h" |
27 | 28 | #include "Framework/AnalysisTask.h" |
@@ -165,6 +166,7 @@ struct CorrelationTask { |
165 | 166 | registry.add("invMass", "2-prong invariant mass (GeV/c^2)", {HistType::kTH3F, {axisInvMassHistogram, axisPtTrigger, axisMultiplicity}}); |
166 | 167 | if (doprocessSame2Prong2Prong || doprocessSame2Prong2ProngML) { |
167 | 168 | registry.add("invMassTwoPart", "2D 2-prong invariant mass (GeV/c^2)", {HistType::kTHnSparseF, {axisInvMassHistogram, axisInvMassHistogram, axisPtTrigger, axisPtAssoc, axisMultiplicity}}); |
| 169 | + registry.add("invMassTwoPartDPhi", "2D 2-prong invariant mass (GeV/c^2)", {HistType::kTHnSparseF, {axisInvMassHistogram, axisInvMassHistogram, axisPtTrigger, axisPtAssoc, axisDeltaPhi}}); |
168 | 170 | } |
169 | 171 | } |
170 | 172 | registry.add("multiplicity", "event multiplicity", {HistType::kTH1F, {{1000, 0, 100, "/multiplicity/centrality"}}}); |
@@ -328,6 +330,7 @@ struct CorrelationTask { |
328 | 330 | if (cfgCorrelationMethod == 2 && track1.decay() == track2.decay()) |
329 | 331 | continue; |
330 | 332 | registry.fill(HIST("invMassTwoPart"), track1.invMass(), track2.invMass(), track1.pt(), track2.pt(), multiplicity); |
| 333 | + registry.fill(HIST("invMassTwoPartDPhi"), track1.invMass(), track2.invMass(), track1.pt(), track2.pt(), TVector2::Phi_0_2pi(track1.phi() - track2.phi() + TMath::Pi() / 2.0) - TMath::Pi() / 2.0); |
331 | 334 | } |
332 | 335 | } |
333 | 336 | } |
|
0 commit comments