@@ -54,10 +54,14 @@ struct flowTest {
5454
5555 void init (InitContext&)
5656 {
57- // pT histograms
57+ // QA and detailed studies
5858 histos.add <TH1>(" hImpactParameter" , " hImpactParameter" , HistType::kTH1D , {axisB});
5959 histos.add <TH2>(" hNchVsImpactParameter" , " hNchVsImpactParameter" , HistType::kTH2D , {axisB, axisNch});
6060 histos.add <TH1>(" hEventPlaneAngle" , " hEventPlaneAngle" , HistType::kTH1D , {axisPhi});
61+ histos.add <TH2>(" hTrackPhiVsEventPlaneAngle" , " hTrackPhiVsEventPlaneAngle" , HistType::kTH2D , {axisPhi, axisPhi});
62+ histos.add <TH2>(" hTrackDeltaPhiVsEventPlaneAngle" , " hTrackDeltaPhiVsEventPlaneAngle" , HistType::kTH2D , {axisPhi, axisPhi});
63+
64+ // analysis
6165 histos.add <TH2>(" hPtVsPhiGenerated" , " hPtVsPhiGenerated" , HistType::kTH2D , {axisPhi, axisPt});
6266 histos.add <TH2>(" hPtVsPhiGlobal" , " hPtVsPhiGlobal" , HistType::kTH2D , {axisPhi, axisPt});
6367 histos.add <TH3>(" hBVsPtVsPhiGenerated" , " hBVsPtVsPhiGenerated" , HistType::kTH3D , {axisB, axisPhi, axisPt});
@@ -83,9 +87,10 @@ struct flowTest {
8387 histos.add <TH3>(" hBVsPtVsPhiGlobalOmega" , " hBVsPtVsPhiGlobalOmega" , HistType::kTH3D , {axisB, axisPhi, axisPt});
8488 }
8589
86- using recoTracks = soa::Join<aod::TracksIU, aod::TracksExtra, aod::McTrackLabels>;
90+ using recoTracks = soa::Join<aod::TracksIU, aod::TracksExtra>;
91+ using recoTracksWithLabels = soa::Join<aod::TracksIU, aod::TracksExtra, aod::McTrackLabels>;
8792
88- void process (aod::McCollision const & mcCollision, soa::Join<aod::McParticles, aod::ParticlesToTracks> const & mcParticles, recoTracks const &)
93+ void process (aod::McCollision const & mcCollision, soa::Join<aod::McParticles, aod::ParticlesToTracks> const & mcParticles, recoTracksWithLabels const &)
8994 {
9095
9196 float imp = mcCollision.impactParameter ();
@@ -118,6 +123,9 @@ struct flowTest {
118123 deltaPhi += 2 . * TMath::Pi ();
119124 if (deltaPhi > 2 . * TMath::Pi ())
120125 deltaPhi -= 2 . * TMath::Pi ();
126+
127+ histos.fill (HIST (" hTrackDeltaPhiVsEventPlaneAngle" ), evPhi, deltaPhi);
128+ histos.fill (HIST (" hTrackPhiVsEventPlaneAngle" ), evPhi, mcParticle.phi ());
121129 histos.fill (HIST (" hPtVsPhiGenerated" ), deltaPhi, mcParticle.pt ());
122130 histos.fill (HIST (" hBVsPtVsPhiGenerated" ), imp, deltaPhi, mcParticle.pt ());
123131
@@ -134,7 +142,7 @@ struct flowTest {
134142 bool validAnalysisTrack = false ;
135143 bool validAnalysisTrackFake = false ;
136144 if (mcParticle.has_tracks ()) {
137- auto const & tracks = mcParticle.tracks_as <recoTracks >();
145+ auto const & tracks = mcParticle.tracks_as <recoTracksWithLabels >();
138146 for (auto const & track : tracks) {
139147 bool isITSFake = false ;
140148
0 commit comments