Update electronSelectionWithTpcEmcal.cxx

Author: rashigupt

PWGHF/HFL/TableProducer/electronSelectionWithTpcEmcal.cxx

@@ -118,10 +118,10 @@ struct HfElectronSelectionWithTpcEmcal {
   Configurable<float> m20EmcClusterElectronMin{"m20EmcClusterElectronMin", 0.0f, "min Electron  EMCal Cluster M20"};
   Configurable<float> tpcNsigmaElectronMin{"tpcNsigmaElectronMin", -0.5f, "min Electron TPCnsigma"};
   Configurable<float> tpcNsigmaElectronMax{"tpcNsigmaElectronMax", 3.0f, "max Electron TPCnsigma"};
-  Configurable<float> charmMinID{"charmMinID", 400.0, "Min Charm Hadron ID"};
-  Configurable<float> charmMaxID{"charmMaxID", 600.0, "Max Charm Hadron ID"};
-  Configurable<float> beautyMinID{"beautyMinID", 4000.0, "Min beauty Hadron ID"};
-  Configurable<float> beautyMaxID{"beautyMaxID", 6000.0, "Max beauty Hadron ID"};
+  Configurable<float> charmMinPdgCode{"charmMinPdgCode", 400.0, "Min Charm Hadron PdgCode"};
+  Configurable<float> charmMaxPdgCode{"charmMaxPdgCode", 600.0, "Max Charm Hadron PdgCode"};
+  Configurable<float> beautyMinPdgCode{"beautyMinPdgCode", 4000.0, "Min beauty Hadron PdgCode"};
+  Configurable<float> beautyMaxPdgCode{"beautyMaxPdgCode", 6000.0, "Max beauty Hadron PdgCode"};
 
   using TableCollisions = o2::soa::Filtered<o2::soa::Join<aod::Collisions, aod::Mults, aod::EvSels>>;
   using TableCollision = TableCollisions::iterator;
@@ -539,7 +539,7 @@ struct HfElectronSelectionWithTpcEmcal {
   {
     fillElectronTrack<false>(collision, tracks, emcClusters, matchedTracks, 0);
   }
-  PROCESS_SWITCH(HfElectronSelectionWithTpcEmcal, processData, "process Data info only", false);
+  PROCESS_SWITCH(HfElectronSelectionWithTpcEmcal, processData, "process Data info only", true);
   ///  Electron selection - for MC reco-level analysis
   void processMcRec(McTableCollision const& mcCollision,
                     McTableTracks const& mcTracks,
@@ -580,8 +580,8 @@ struct HfElectronSelectionWithTpcEmcal {
           //=================  eta->e ======================================
           if (std::abs(mother.pdgCode()) == kEta) {
             if (mother.isPhysicalPrimary()) {
-              if ((std::abs(gmother.pdgCode()) >= charmMinID && std::abs(gmother.pdgCode()) < charmMaxID) ||
-                  (std::abs(gmother.pdgCode()) >= beautyMinID && std::abs(gmother.pdgCode()) < beautyMaxID)) {
+              if ((std::abs(gmother.pdgCode()) >= charmMinPdgCode && std::abs(gmother.pdgCode()) < charmMaxPdgCode) ||
+                  (std::abs(gmother.pdgCode()) >= beautyMinPdgCode && std::abs(gmother.pdgCode()) < beautyMaxPdgCode)) {
                 continue;
               }
               isEmbEta = true;
@@ -591,17 +591,17 @@ struct HfElectronSelectionWithTpcEmcal {
 
           if (std::abs(mother.pdgCode()) == kPi0) {
             if (mother.isPhysicalPrimary()) {
-              if ((std::abs(gmother.pdgCode()) >= charmMinID && std::abs(gmother.pdgCode()) < charmMaxID) ||
-                  (std::abs(gmother.pdgCode()) >= beautyMinID && std::abs(gmother.pdgCode()) < beautyMaxID)) {
+              if ((std::abs(gmother.pdgCode()) >= charmMinPdgCode && std::abs(gmother.pdgCode()) < charmMaxPdgCode) ||
+                  (std::abs(gmother.pdgCode()) >= beautyMinPdgCode && std::abs(gmother.pdgCode()) < beautyMaxPdgCode)) {
                 continue;
               }
               isEmbPi0 = true; // pi0 -> e
             }
 
             if (std::abs(gmother.pdgCode()) == kEta) {
               if (gmother.isPhysicalPrimary()) {
-                if ((std::abs(ggmother.pdgCode()) >= charmMinID && std::abs(ggmother.pdgCode()) < charmMaxID) ||
-                    (std::abs(ggmother.pdgCode()) >= beautyMinID && std::abs(ggmother.pdgCode()) < beautyMaxID)) {
+                if ((std::abs(ggmother.pdgCode()) >= charmMinPdgCode && std::abs(ggmother.pdgCode()) < charmMaxPdgCode) ||
+                    (std::abs(ggmother.pdgCode()) >= beautyMinPdgCode && std::abs(ggmother.pdgCode()) < beautyMaxPdgCode)) {
                   continue;
                 }
                 isEmbEta = true; // eta->pi0-> e
@@ -613,8 +613,8 @@ struct HfElectronSelectionWithTpcEmcal {
           if (std::abs(mother.pdgCode()) == kGamma) {
             if (std::abs(gmother.pdgCode()) == kEta) {
               if (gmother.isPhysicalPrimary()) {
-                if ((std::abs(ggmother.pdgCode()) >= charmMinID && std::abs(ggmother.pdgCode()) < charmMaxID) ||
-                    (std::abs(ggmother.pdgCode()) >= beautyMinID && std::abs(ggmother.pdgCode()) < beautyMaxID)) {
+                if ((std::abs(ggmother.pdgCode()) >= charmMinPdgCode && std::abs(ggmother.pdgCode()) < charmMaxPdgCode) ||
+                    (std::abs(ggmother.pdgCode()) >= beautyMinPdgCode && std::abs(ggmother.pdgCode()) < beautyMaxPdgCode)) {
                   continue;
                 }
                 isEmbEta = true; // eta->gamma-> e
@@ -623,17 +623,17 @@ struct HfElectronSelectionWithTpcEmcal {
 
             if (std::abs(gmother.pdgCode()) == kPi0) {
               if (gmother.isPhysicalPrimary()) {
-                if ((std::abs(ggmother.pdgCode()) >= charmMinID && std::abs(ggmother.pdgCode()) < charmMaxID) ||
-                    (std::abs(ggmother.pdgCode()) >= beautyMinID && std::abs(ggmother.pdgCode()) < beautyMaxID)) {
+                if ((std::abs(ggmother.pdgCode()) >= charmMinPdgCode && std::abs(ggmother.pdgCode()) < charmMaxPdgCode) ||
+                    (std::abs(ggmother.pdgCode()) >= beautyMinPdgCode && std::abs(ggmother.pdgCode()) < beautyMaxPdgCode)) {
                   continue;
                 }
                 isEmbPi0 = true; // pi0-> gamma-> e
               }
 
               if (std::abs(ggmother.pdgCode()) == kEta) {
                 if (ggmother.isPhysicalPrimary()) {
-                  if ((std::abs(gggmother.pdgCode()) >= charmMinID && std::abs(gggmother.pdgCode()) < charmMaxID) ||
-                      (std::abs(gggmother.pdgCode()) >= beautyMinID && std::abs(gggmother.pdgCode()) < beautyMaxID)) {
+                  if ((std::abs(gggmother.pdgCode()) >= charmMinPdgCode && std::abs(gggmother.pdgCode()) < charmMaxPdgCode) ||
+                      (std::abs(gggmother.pdgCode()) >= beautyMinPdgCode && std::abs(gggmother.pdgCode()) < beautyMaxPdgCode)) {
                     continue;
                   }
                   isEmbEta = true; // eta->pi0->gamma-> e
@@ -659,7 +659,7 @@ struct HfElectronSelectionWithTpcEmcal {
     }
   }
 
-  PROCESS_SWITCH(HfElectronSelectionWithTpcEmcal, processMcGen, "Process MC Gen mode", true);
+  PROCESS_SWITCH(HfElectronSelectionWithTpcEmcal, processMcGen, "Process MC Gen mode", false);
 };
 WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
 {