Fixing reserve issue and adding extra cut

Author: lucamicheletti

PWGDQ/Core/CutsLibrary.cxx

@@ -3067,6 +3067,11 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
     return cut;
   }
 
+  if (nameStr == "muonEtaCut") {
+    cut->AddCut(GetAnalysisCut("muonEtaCut"));
+    return cut;
+  }
+
   if (nameStr == "muonMinimalCuts") {
     cut->AddCut(GetAnalysisCut("muonMinimalCuts"));
     return cut;
@@ -6558,6 +6563,11 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName)
     return cut;
   }
 
+  if (nameStr == "muonEtaCut") {
+    cut->AddCut(VarManager::kEta, -4.0, -2.5);
+    return cut;
+  }
+
   if (nameStr == "muonMinimalCuts") {
     cut->AddCut(VarManager::kEta, -4.0, -2.5);
     cut->AddCut(VarManager::kMuonRAtAbsorberEnd, 17.6, 89.5);

PWGDQ/TableProducer/tableMakerMuonMchTrkEfficiency.cxx

@@ -420,7 +420,7 @@ struct tableMakerMuonMchTrkEfficiency {
     for (int i = 0; i < ncuts; i++)
       nselmuons.push_back(0);
 
-    rowCandidateBase.reserve(tracksMuon.size());
+    //rowCandidateBase.reserve(tracksMuon.size());
     for (auto& muon : tracksMuon) {
 
       VarManager::FillTrack<TMuonFillMap>(muon);
@@ -486,8 +486,8 @@ struct tableMakerMuonMchTrkEfficiency {
     for (int i = 0; i < ncuts; i++)
       nselmuons.push_back(0);
 
-    rowCandidateBase.reserve(tracksMuon.size());
-    rowCandidateGen.reserve(tracksMuon.size());
+    //rowCandidateBase.reserve(tracksMuon.size());
+    //rowCandidateGen.reserve(tracksMuon.size());
     for (auto& muon : tracksMuon) {
       ///
       /// First compute MC matched quantities using either the DQ skimmed or the Framework data models

PWGDQ/Tasks/taskFwdTrackPid.cxx

@@ -131,7 +131,7 @@ struct taskFwdTrackPid {
   template <bool TMatchedOnly, uint32_t TEventFillMap, uint32_t TTrackFillMap, typename TEvent, typename Muons, typename MftTracks>
   void runFwdTrackPid(TEvent const& event, Muons const& muons, MftTracks const& mftTracks)
   {
-    fwdPidAllList.reserve(1);
+    //fwdPidAllList.reserve(1);
     for (const auto& muon : muons) {
       if (muon.has_matchMFTTrack() && muon.trackType() == 0 && TMath::Abs(muon.fwdDcaX()) < fConfigMaxDCA && TMath::Abs(muon.fwdDcaY()) < fConfigMaxDCA) {
         auto mftTrack = muon.template matchMFTTrack_as<MyMftTracks>();
@@ -157,7 +157,7 @@ struct taskFwdTrackPid {
   template <bool TMatchedOnly, uint32_t TEventFillMap, uint32_t TEventMCFillMap, uint32_t TTrackFillMap, typename TEvent, typename Muons, typename MftTracks, typename TEventsMC, typename TTracksMC>
   void runFwdTrackPidMC(TEvent const& event, Muons const& muons, MftTracks const& mftTracks, TEventsMC const& /*eventsMC*/, TTracksMC const& /*tracksMC*/)
   {
-    fwdPidAllList.reserve(1);
+    //fwdPidAllList.reserve(1);
     for (const auto& muon : muons) {
       if (muon.has_matchMFTTrack() && muon.trackType() == 0 && TMath::Abs(muon.fwdDcaX()) < fConfigMaxDCA && TMath::Abs(muon.fwdDcaY()) < fConfigMaxDCA) {
         auto mftTrack = muon.template matchMFTTrack_as<MyMftTracksMC>();