Merge pull request cms-sw#123 from kurtejung/fixHiGenEtaFilter2

fixing HiGen eta filter
mverwe · Jun 22, 2017 · a6071aa · a6071aa
2 parents 3669770 + 0979968
commit a6071aa
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/HeavyIonsAnalysis/JetAnalysis/python/HiGenAnalyzer_cfi.py b/HeavyIonsAnalysis/JetAnalysis/python/HiGenAnalyzer_cfi.py
@@ -2,7 +2,7 @@
 
 HiGenParticleAna = cms.EDAnalyzer('HiGenAnalyzer',
 				  doVertex = cms.untracked.bool(False),
-				  etaMax = cms.untracked.double(2),
+				  etaMax = cms.untracked.double(2.5),
 				  ptMin = cms.untracked.double(5),
 				  chargedOnly = cms.untracked.bool(False),
 				  stableOnly = cms.untracked.bool(True),

diff --git a/HeavyIonsAnalysis/TrackAnalysis/src/HiGenAnalyzer.cc b/HeavyIonsAnalysis/TrackAnalysis/src/HiGenAnalyzer.cc
@@ -322,6 +322,7 @@ HiGenAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
       nparticles++;
       if ((*it)->momentum().perp()<ptMin_) continue;
       //	if((*it)->status() == 1){
+      if (fabs((*it)->momentum().eta())>etaMax_) continue; 
       Int_t pdg_id = (*it)->pdg_id();
       Float_t eta = (*it)->momentum().eta();
       Float_t phi = (*it)->momentum().phi();
@@ -357,6 +358,7 @@ HiGenAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
       const reco::GenParticle& p = (*parts)[i];
       if (stableOnly_ && p.status()!=1) continue;
       if (p.pt()<ptMin_) continue;
+      if (fabs(p.eta())>etaMax_) continue;
       if (chargedOnly_&&p.charge()==0) continue;
       hev_.pt.push_back( p.pt());
       hev_.eta.push_back( p.eta());