Issue fix in Exotica HLT validation code

Author: borzari

HLTriggerOffline/Exotica/src/HLTExoticaPlotter.cc

@@ -126,9 +126,14 @@ void HLTExoticaPlotter::analyze(const bool &isPassTrigger,
     totalobjectssize++;
   totalobjectssize *= countobjects.size();
   // Fill the histos if pass the trigger (just the two with higher pt)
+  unsigned int jaux = 0;
+  // jaux is being used as a dedicated counter to avoid getting
+  // a non-existent element inside dxys
+  // more information in the issue https://github.com/cms-sw/cmssw/issues/32550
   for (size_t j = 0; j < matches.size(); ++j) {
     // Is this object owned by this trigger? If not we are not interested...
     if (_objectsType.find(matches[j].pdgId()) == _objectsType.end()) {
+      ++jaux;
       continue;
     }
 
@@ -149,8 +154,10 @@ void HLTExoticaPlotter::analyze(const bool &isPassTrigger,
     }
 
     if (!dxys.empty() &&
-        (objType == EVTColContainer::ELEC || objType == EVTColContainer::MUON || objType == EVTColContainer::MUTRK))
-      this->fillHist(isPassTrigger, source, objTypeStr, "Dxy", dxys[j]);
+        (objType == EVTColContainer::ELEC || objType == EVTColContainer::MUON || objType == EVTColContainer::MUTRK)){
+      this->fillHist(isPassTrigger, source, objTypeStr, "Dxy", dxys[jaux]);
+      ++jaux;
+    }
 
     if (countobjects[objType] == 0) {
       if (!(TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT")) || source != "gen") {

HLTriggerOffline/Exotica/src/HLTExoticaSubAnalysis.cc

@@ -659,9 +659,18 @@ void HLTExoticaSubAnalysis::analyze(const edm::Event &iEvent, const edm::EventSe
       }
     }
 
+    int jel = 0;
+    int jmu = 0;
+    int jmutrk = 0;
+
+    // jel, jmu and jmutrk are being used as a dedicated counters to avoid getting
+    // non-existent elements inside trkObjs[11], trkObjs[13] and trkObjs[130], respectively
+    // more information in the issue https://github.com/cms-sw/cmssw/issues/32550
+
     for (size_t j = 0; (j != matchesReco.size()) && isPassedLeadingCut; ++j) {
       const unsigned int objType = matchesReco[j].pdgId();
-      // std::cout << "(4) Gonna call with " << objType << std::endl;
+      //std::cout << "(4) Gonna call with " << objType << std::endl;
+
       const std::string objTypeStr = EVTColContainer::getTypeString(objType);
 
       float pt = matchesReco[j].pt();
@@ -703,10 +712,25 @@ void HLTExoticaSubAnalysis::analyze(const edm::Event &iEvent, const edm::EventSe
         this->fillHist("rec", objTypeStr, "SumEt", theSumEt[objType]);
       }
 
-      if (trkObjs[objType].size() >= j + 1) {
-        float dxyRec = trkObjs[objType].at(j)->dxy(cols->bs->position());
+      if(objType == 11){
+        float dxyRec = trkObjs[objType].at(jel)->dxy(cols->bs->position());
+        this->fillHist("rec", objTypeStr, "Dxy", dxyRec);
+        dxys.push_back(dxyRec);
+        ++jel;
+      }
+
+      if(objType == 13){
+        float dxyRec = trkObjs[objType].at(jmu)->dxy(cols->bs->position());
+        this->fillHist("rec", objTypeStr, "Dxy", dxyRec);
+        dxys.push_back(dxyRec);
+        ++jmu;
+      }
+
+      if(objType == 130){
+        float dxyRec = trkObjs[objType].at(jmutrk)->dxy(cols->bs->position());
         this->fillHist("rec", objTypeStr, "Dxy", dxyRec);
         dxys.push_back(dxyRec);
+        ++jmutrk;
       }
 
     }  // Closes loop in reco