Update flowTest.cxx

Author: ddobrigk

Common/Tasks/flowTest.cxx

@@ -59,8 +59,9 @@ struct flowTest {
     histos.add<TH2>("hNchVsImpactParameter", "hNchVsImpactParameter", HistType::kTH2D, {axisB, axisNch});
     histos.add<TH1>("hEventPlaneAngle", "hEventPlaneAngle", HistType::kTH1D, {axisPhi});
     histos.add<TH2>("hTrackPhiVsEventPlaneAngle", "hTrackPhiVsEventPlaneAngle", HistType::kTH2D, {axisPhi, axisPhi});
+    histos.add<TH2>("hTrackDeltaPhiVsEventPlaneAngle", "hTrackDeltaPhiVsEventPlaneAngle", HistType::kTH2D, {axisPhi, axisPhi});
 
-    // analysis
+    // analysis 
     histos.add<TH2>("hPtVsPhiGenerated", "hPtVsPhiGenerated", HistType::kTH2D, {axisPhi, axisPt});
     histos.add<TH2>("hPtVsPhiGlobal", "hPtVsPhiGlobal", HistType::kTH2D, {axisPhi, axisPt});
     histos.add<TH3>("hBVsPtVsPhiGenerated", "hBVsPtVsPhiGenerated", HistType::kTH3D, {axisB, axisPhi, axisPt});
@@ -123,6 +124,7 @@ struct flowTest {
         if (deltaPhi > 2. * TMath::Pi())
           deltaPhi -= 2. * TMath::Pi();
 
+        histos.fill(HIST("hTrackDeltaPhiVsEventPlaneAngle"), evPhi, deltaPhi);
         histos.fill(HIST("hTrackPhiVsEventPlaneAngle"), evPhi, mcParticle.phi());
         histos.fill(HIST("hPtVsPhiGenerated"), deltaPhi, mcParticle.pt());
         histos.fill(HIST("hBVsPtVsPhiGenerated"), imp, deltaPhi, mcParticle.pt());
@@ -137,28 +139,28 @@ struct flowTest {
         bool validITSTrackFake = false;
         bool validITSABTrack = false;
         bool validITSABTrackFake = false;
-        bool validAnalysisTrack = false;
+        bool validAnalysisTrack = false; 
         bool validAnalysisTrackFake = false;
         if (mcParticle.has_tracks()) {
           auto const& tracks = mcParticle.tracks_as<recoTracksWithLabels>();
           for (auto const& track : tracks) {
             bool isITSFake = false;
 
-            for (int bit = 0; bit < 7; bit++) {
-              if (bitcheck(track.mcMask(), bit)) {
+            for(int bit = 0; bit < 7; bit++){ 
+              if(bitcheck(track.mcMask(), bit)){ 
                 isITSFake = true;
               }
             }
 
-            if (track.tpcNClsFound() > analysisMinimumTPCClusters && track.itsNCls() > analysisMinimumITSClusters) {
+            if (track.tpcNClsFound() > analysisMinimumTPCClusters && track.itsNCls() > analysisMinimumITSClusters){
               validAnalysisTrack = true;
-              if (isITSFake) {
+              if(isITSFake){ 
                 validAnalysisTrackFake = true;
               }
             }
             if (track.hasTPC() && track.hasITS()) {
               validGlobal = true;
-              if (isITSFake) {
+              if(isITSFake){ 
                 validGlobalFake = true;
               }
             }
@@ -170,13 +172,13 @@ struct flowTest {
             }
             if (track.hasITS() && track.itsChi2NCl() > -1e-6) {
               validITSTrack = true;
-              if (isITSFake) {
+              if(isITSFake){ 
                 validITSTrackFake = true;
               }
             }
             if (track.hasITS() && track.itsChi2NCl() < -1e-6) {
               validITSABTrack = true;
-              if (isITSFake) {
+              if(isITSFake){ 
                 validITSABTrackFake = true;
               }
             }
@@ -188,13 +190,13 @@ struct flowTest {
           histos.fill(HIST("hPtVsPhiGlobal"), deltaPhi, mcParticle.pt());
           histos.fill(HIST("hBVsPtVsPhiGlobal"), imp, deltaPhi, mcParticle.pt());
         }
-        if (validGlobalFake) {
+        if (validGlobalFake){
           histos.fill(HIST("hBVsPtVsPhiGlobalFake"), imp, deltaPhi, mcParticle.pt());
         }
         if (validAnalysisTrack) {
           histos.fill(HIST("hBVsPtVsPhiAnalysis"), imp, deltaPhi, mcParticle.pt());
         }
-        if (validAnalysisTrackFake) {
+        if (validAnalysisTrackFake){
           histos.fill(HIST("hBVsPtVsPhiAnalysisFake"), imp, deltaPhi, mcParticle.pt());
         }
         // if any track present, fill