[PWGCF] calculate and apply V0s and cascades local density efficiency

calculate and apply V0s and cascades local density efficiency

Author: fuchuncui

PWGCF/Flow/Tasks/CMakeLists.txt

@@ -28,6 +28,11 @@ o2physics_add_dpl_workflow(flow-mc
                     SOURCES flowMc.cxx
                     PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::GFWCore
                     COMPONENT_NAME Analysis)
+                    
+o2physics_add_dpl_workflow(flow-qa
+                    SOURCES flowQa.cxx
+                    PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::AnalysisCCDB O2Physics::GFWCore
+                    COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(flow-gfw-task
                     SOURCES flowGfwTask.cxx

PWGCF/Flow/Tasks/flowEfficiencyCasc.cxx

@@ -80,15 +80,16 @@ struct flowEfficiencyCasc {
   using myCollisions = soa::Join<aod::StraCollisions, aod::StraEvSels>;
   using myMcCollisions = soa::Join<aod::StraMCCollisions, aod::StraMCCollMults>;
   using CascMCCandidates = soa::Join<aod::CascCollRefs, aod::CascCores, aod::CascExtras, aod::CascBBs, aod::CascCoreMCLabels>;
-  using V0MCCandidates =soa::Join<aod::V0CollRefs, aod::V0Cores, aod::V0Extras, aod::V0CoreMCLabels>;
+  using V0MCCandidates = soa::Join<aod::V0CollRefs, aod::V0Cores, aod::V0Extras, aod::V0CoreMCLabels>;
   using DaughterTracks = soa::Join<aod::DauTrackExtras, aod::DauTrackTPCPIDs>;
 
   // Define the output
   HistogramRegistry registry{"registry"};
 
   std::vector<float> cfgNSigma = cfgNSigmatpctof;
 
-  void init(InitContext const&) {
+  void init(InitContext const&)
+  {
     const AxisSpec axisCounter{1, 0, +1, ""};
     // create histograms
     registry.add("eventCounter", "eventCounter", kTH1F, {axisCounter});
@@ -177,7 +178,6 @@ struct flowEfficiencyCasc {
                  (!cfgcheckDauTPC || (std::fabs(bachelor.tpcNSigmaPi()) < cfgNSigma[0] && std::fabs(negdau.tpcNSigmaPr()) < cfgNSigma[1] && std::fabs(posdau.tpcNSigmaPi()) < cfgNSigma[0]))) {
         registry.fill(HIST("h2DRecXi"), casc.pt(), rectracknum);
       }
-
     }
 
     for (auto& v0 : V0s) {

PWGCF/Flow/Tasks/flowGFWOmegaXi.cxx

@@ -677,7 +677,7 @@ struct FlowGFWOmegaXi {
       return;
     if (eventSelected(collision, tracks.size(), cent))
       return;
-    TH1D *hLocalDensity = new TH1D("hphi", "hphi", 400, -constants::math::TwoPI, constants::math::TwoPI);
+    TH1D* hLocalDensity = new TH1D("hphi", "hphi", 400, -constants::math::TwoPI, constants::math::TwoPI);
     auto bc = collision.bc_as<aod::BCsWithTimestamps>();
     loadCorrections(bc.timestamp());
     float vtxz = collision.posZ();
@@ -814,7 +814,7 @@ struct FlowGFWOmegaXi {
           int density = hLocalDensity->Integral(phibin - cfgDeltaPhiLocDen, phibin + cfgDeltaPhiLocDen);
           setCurrentLocalDensityWeights(wloc, v0, density, 2);
         }
-        
+
         candNum[1] = candNum[1] + 1;
         registry.fill(HIST("InvMassLambda"), v0.pt(), v0.mLambda(), v0.eta(), cent);
         registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.phi(), v0.eta(), vtxz, wacc);
@@ -920,7 +920,7 @@ struct FlowGFWOmegaXi {
           phibin = hLocalDensity->FindBin(casc.phi() - constants::math::TwoPI);
         if (phibin > -900) {
           int density = hLocalDensity->Integral(phibin - cfgDeltaPhiLocDen, phibin + cfgDeltaPhiLocDen);
-          setCurrentLocalDensityWeights(wloc, casc, density, 3);          
+          setCurrentLocalDensityWeights(wloc, casc, density, 3);
         }
 
         candNum[2] = candNum[2] + 1;
@@ -939,7 +939,7 @@ struct FlowGFWOmegaXi {
         registry.fill(HIST("hEventCount"), 3.5, i + 0.5);
       }
     }
-    delete hLocalDensity; 
+    delete hLocalDensity;
     // Filling cumulant with ROOT TProfile and loop for all ptBins
     fillProfile(corrconfigs.at(15), HIST("c22"), cent);
     fillProfile(corrconfigs.at(16), HIST("c24"), cent);