Add specific tau PV producer for taus produced from MiniAOD input (ba…

…sed on work by Michal Bluj), and adapt MiniAOD tau sequence accordingly
cms-sw · Feb 28, 2018 · 5294f0e · mbluj · Mar 1, 2018 · 5294f0e
1 parent 99a3384
commit 5294f0e
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/RecoTauTag/Configuration/python/tools/adaptToRunAtMiniAOD.py b/RecoTauTag/Configuration/python/tools/adaptToRunAtMiniAOD.py
@@ -182,6 +182,7 @@ def adaptTauToMiniAODReReco(process, reclusterJets=True):
 	convertModuleToMiniAODInput(process, 'hpsPFTauPhotonPtSumOutsideSignalConedR03')
 
 	# Redefine tau PV producer
+	process.hpsPFTauPrimaryVertexProducer.__dict__['_TypedParameterizable__type'] = 'PFTauMiniAODPrimaryVertexProducer'
 	process.hpsPFTauPrimaryVertexProducer.PVTag = 'offlineSlimmedPrimaryVertices'
 	process.hpsPFTauPrimaryVertexProducer.qualityCuts.primaryVertexSrc = 'offlineSlimmedPrimaryVertices'
 	process.hpsPFTauPrimaryVertexProducer.packedCandidatesTag = cms.InputTag("packedPFCandidates")
@@ -298,6 +299,10 @@ def setOutputModule(mode=0):
 			)
 		)
 	output.outputCommands.append('keep *_selectedPatTaus_*_*')
+	output.outputCommands.append('keep *_combinatoricReco*_*_*')
+	output.outputCommands.append('keep *_ak4PFJetsRecoTauChargedHadrons_*_*')
+	output.outputCommands.append('keep *_ak4PFJetsLegacyHPSPiZeros_*_*')
+	output.outputCommands.append('keep *_patAK4PFJets_*_*')
 	if mode==1:
                 for prod in evtContent.RecoTauTagAOD.outputCommands:
 			if prod.find('ElectronRejection') > -1:

diff --git a/RecoTauTag/RecoTau/plugins/PFTauPrimaryVertexProducer.cc b/RecoTauTag/RecoTau/plugins/PFTauPrimaryVertexProducer.cc
@@ -195,12 +195,7 @@ void PFTauPrimaryVertexProducer::produce(edm::Event& iEvent,const edm::EventSetu
       }
       if (passed && cut_.get()){passed = (*cut_)(*tau);}
       if (passed){
-      	// Use two signal track collections: One with Refs for taus made from PFCandidates,
-      	// and one with pointers for taus made from PackedCandidates
 	std::vector<reco::TrackBaseRef> SignalTracks;
-	// std::vector<reco::Track*> SignalTrackPtrs;
-	// JAN - FIXME - THIS NEEDS TO BE CHANGED WITH MICHAL'S ADAPTIONS
-
 	for(reco::PFTauCollection::size_type jPFTau = 0; jPFTau < Tau->size(); jPFTau++) {
 	  if(useSelectedTaus_ || iPFTau==jPFTau){
 	    reco::PFTauRef RefPFTau(Tau, jPFTau);