fix nummuons and numtaus partially addressing #44

cms-opendata-analyses · Jul 17, 2021 · 4510f5e · 4510f5e
1 parent ad87943
commit 4510f5e
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 6 deletions.
diff --git a/PhysObjectExtractor/src/MuonAnalyzer.cc b/PhysObjectExtractor/src/MuonAnalyzer.cc
@@ -187,8 +187,6 @@ MuonAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
 
   if(mymuons.isValid()){
 
-     // get the number of muons in the event
-     nummuon=(*mymuons).size();
      const float mu_min_pt = 3;
      math::XYZPoint pv(vertices->begin()->position());
 
@@ -231,7 +229,8 @@ MuonAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
         	    }
         	    muon_genpartidx.push_back(-1);
         	    muon_jetidx.push_back(-1);
-          } 
+		    nummuon++;
+	   } 
         }
   }
 

diff --git a/PhysObjectExtractor/src/TauAnalyzer.cc b/PhysObjectExtractor/src/TauAnalyzer.cc
@@ -194,9 +194,7 @@ TauAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
   tau_idantimutight.clear();
 
   if(mytaus.isValid()){
-    // get the number of taus in the event
-    numtau=mytaus->size();
-
+
     const float tau_min_pt = 15;
     for (reco::PFTauCollection::const_iterator itTau=mytaus->begin(); itTau!=mytaus->end(); ++itTau){
       if (itTau->pt() > tau_min_pt) {
@@ -222,6 +220,7 @@ TauAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
         tau_idantimutight.push_back(tausTightMuonRej->operator[](idx).second);
 
        	tau_reliso_all.push_back((itTau->isolationPFChargedHadrCandsPtSum() + itTau->isolationPFGammaCandsEtSum()) / itTau->pt());
+	numtau++;
       }   	
     }
   }