blinkseb/patch_produce_PAT_COMMON.py

## patch_produce_PAT_COMMON.py
diff --git a/crab/produce_PAT_COMMON.py b/crab/produce_PAT_COMMON.py
index 082490f..8e3f4e1 100644
--- a/crab/produce_PAT_COMMON.py
+++ b/crab/produce_PAT_COMMON.py
@@ -40,8 +40,18 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)
   # this function will modify the PAT sequences.
   from PhysicsTools.PatAlgos.tools.coreTools import removeSpecificPATObjects, removeMCMatching
   from PhysicsTools.PatAlgos.tools.pfTools import usePF2PAT, adaptPFIsoElectrons, adaptPVs, usePFIso
+  from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection
   #from PhysicsTools.PatAlgos.tools.metTools import *

+  process.load("RecoJets.Configuration.GenJetParticles_cff")
+  from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
+  process.ca8GenJetsNoNu = ca4GenJets.clone(
+      rParam = cms.double(0.8),
+      src = cms.InputTag("genParticlesForJetsNoNu")
+      )
+
+  from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets
+
   def usePF2PATForAnalysis(jetAlgo, postfix, useTypeIMET, usePFNoPU):

     # Jet corrections
@@ -137,6 +147,41 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)

     return getattr(process, "patPF2PATSequence" + p)

+  # This must be called after all the calls to usePF2PAT
+  def addC8Jets(p):
+    p2 = p.replace('AK5', 'CA8')
+
+    ###############################
+    ###### Bare CA 0.8 jets #######
+    ###############################
+    setattr(process, "ca8PFJets" + p2, ca4PFJets.clone(
+            rParam = cms.double(0.8),
+            src = cms.InputTag('pfNoElectron' + p),
+            doAreaFastjet = cms.bool(True),
+            doRhoFastjet = cms.bool(True),
+            Rho_EtaMax = cms.double(6.0),
+            Ghost_EtaMax = cms.double(7.0),
+            srcPVs = cms.InputTag("goodOfflinePrimaryVertices")
+            ))
+
+    # Add ca8PFJets to the main process path
+    getattr(process, "patPF2PATSequence" + p).replace(
+            getattr(process, "pfNoElectron" + p), getattr(process, "pfNoElectron" + p) * getattr(process, "ca8PFJets" + p2)
+            )
+
+    addJetCollection(process,
+            cms.InputTag('ca8PFJets' + p2),
+            'CA8', 'PFchs' if 'chs' in p else 'PF',
+            doJTA = True,
+            doBTagging = True,
+            jetCorrLabel = ('AK7PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute']),
+            doType1MET = True,
+            doL1Cleaning = False,
+            doL1Counters = False,
+            genJetCollection = cms.InputTag("ca8GenJetsNoNu"),
+            doJetID = True,
+            )
+
   if correctMETWithT1:
     process.load("PhysicsTools.PatUtils.patPFMETCorrections_cff")
     #from PhysicsTools.PatAlgos.tools.helpers import cloneProcessingSnippet
@@ -178,6 +223,11 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)

       process.QuarkGluonTagger.replace(getattr(process, 'QGTagger' + p), getattr(process, 'QGTagger' + p) + getattr(process, 'QGTagger' + chsP))

+  for p, algo in postfixes.items():
+    addC8Jets(p)
+    if runCHS:
+      addC8Jets(p + "chs")
+
   print "##########################"
   print "Calo jets" if processCaloJets else "No processing of calo jets"
   print "##########################"
@@ -233,6 +283,8 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)

   process.analysisSequence = cms.Sequence()

+  process.analysisSequence *= process.genParticlesForJetsNoNu * process.ca8GenJetsNoNu
+
   process.analysisSequence *= process.sequence_nochs
   if runCHS:
     process.analysisSequence *= process.sequence_chs

## produce_PAT_COMMON.py
import FWCore.ParameterSet.Config as cms

def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag):

  ## import skeleton process
  from PhysicsTools.PatAlgos.patTemplate_cfg import process

  # load the PAT config
  process.load("PhysicsTools.PatAlgos.patSequences_cff")
  process.load("JetMETCorrections.Configuration.JetCorrectionProducers_cff")

  # Do some CHS stuff
  process.ak5PFchsL1Fastjet  = process.ak5PFL1Fastjet.clone(algorithm = 'AK5PFchs')
  process.ak5PFchsL2Relative = process.ak5PFL2Relative.clone(algorithm = 'AK5PFchs')
  process.ak5PFchsL3Absolute = process.ak5PFL3Absolute.clone(algorithm = 'AK5PFchs')
  process.ak5PFchsResidual   = process.ak5PFResidual.clone(algorithm = 'AK5PFchs')
  process.ak5PFchsL1FastL2L3 = cms.ESProducer(
      'JetCorrectionESChain',
      correctors = cms.vstring('ak5PFchsL1Fastjet', 'ak5PFchsL2Relative','ak5PFchsL3Absolute')
      )


  from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector

  ## The good primary vertex filter ____________________________________________||
  ## This filter throw events with no primary vertex
  process.primaryVertexFilter = cms.EDFilter(
      "VertexSelector",
      src = cms.InputTag("offlinePrimaryVertices"),
      cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.Rho <= 2"),
      filter = cms.bool(True)
      )

  process.goodOfflinePrimaryVertices = cms.EDFilter("PrimaryVertexObjectFilter",
      filterParams = pvSelector.clone( minNdof = cms.double(4.0), maxZ = cms.double(24.0) ),
      src = cms.InputTag('offlinePrimaryVertices')
      )

  # Configure PAT to use PF2PAT instead of AOD sources
  # this function will modify the PAT sequences.
  from PhysicsTools.PatAlgos.tools.coreTools import removeSpecificPATObjects, removeMCMatching
  from PhysicsTools.PatAlgos.tools.pfTools import usePF2PAT, adaptPFIsoElectrons, adaptPVs, usePFIso
  from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection
  #from PhysicsTools.PatAlgos.tools.metTools import *

  process.load("RecoJets.Configuration.GenJetParticles_cff")
  from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
  process.ca8GenJetsNoNu = ca4GenJets.clone(
      rParam = cms.double(0.8),
      src = cms.InputTag("genParticlesForJetsNoNu")
      )

  from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets

  def usePF2PATForAnalysis(jetAlgo, postfix, useTypeIMET, usePFNoPU):

    # Jet corrections
    # No L2L3 Residual on purpose
    if usePFNoPU:
      jetCorrections = ("%sPFchs" % algo, ['L1FastJet', 'L2Relative', 'L3Absolute'])
      postfix += "chs"
    else:
      jetCorrections = ("%sPF" % algo, ['L1FastJet', 'L2Relative', 'L3Absolute'])

    #if not runOnMC:
    #  jetCorrections[1].append('L2L3Residual')

    p = postfix

    usePF2PAT(process, runPF2PAT=True, jetAlgo=jetAlgo, runOnMC=runOnMC, postfix=p, jetCorrections=jetCorrections, typeIMetCorrections = useTypeIMET)
    getattr(process, "pfPileUp" + p).Enable = True
    getattr(process, "pfPileUp" + p).Vertices = 'goodOfflinePrimaryVertices'
    getattr(process, "pfPileUp" + p).checkClosestZVertex = cms.bool(False)

    getattr(process, "pfJets" + p).doAreaFastjet = cms.bool(True)
    getattr(process, "pfJets" + p).doRhoFastjet = False
    getattr(process, 'patJetCorrFactors' + p).rho = cms.InputTag("kt6PFJets", "rho") # Do not use kt6PFJetsPFlowAK5, it's not ok for L1FastJet.

    # top projections in PF2PAT:
    getattr(process,"pfNoPileUp" + p).enable = cms.bool(usePFNoPU)
    getattr(process,"pfNoMuon" + p).enable = True
    getattr(process,"pfNoElectron" + p).enable = True
    getattr(process,"pfNoTau" + p).enable = False
    getattr(process,"pfNoJet" + p).enable = True

    getattr(process,"patElectrons" + p).embedTrack = True

    getattr(process,"patMuons" + p).embedTrack = True
    # enable delta beta correction for muon selection in PF2PAT?
    getattr(process,"pfIsolatedMuons" + p).doDeltaBetaCorrection = True

    getattr(process, "patJets" + p).embedPFCandidates = False
    # Keep only jets with pt > 2 Gev
    getattr(process, "selectedPatJets" + p).cut = "pt > 2";

    # Use a cone of 0.3 for photon isolation
    #adaptPFIsoPhotons(process, applyPostfix(process, "patPhotons", postfix), postfix, "03")

    # 2012 Photon ID

    # Electron conversion
    setattr(process, "patConversions" + p, cms.EDProducer("PATConversionProducer",
        # input collection
        electronSource = cms.InputTag("selectedPatElectrons" + p)
    ))

    # Switch electron isolation to dR = 0.3, for PF2PAT
    getattr(process, "pfIsolatedElectrons" + p).isolationValueMapsCharged = cms.VInputTag(cms.InputTag("elPFIsoValueCharged03PFId" + p))
    getattr(process, "pfIsolatedElectrons" + p).deltaBetaIsolationValueMap = cms.InputTag("elPFIsoValuePU03PFId" + p)
    getattr(process, "pfIsolatedElectrons" + p).isolationValueMapsNeutral = cms.VInputTag(cms.InputTag("elPFIsoValueNeutral03PFId" + p), cms.InputTag("elPFIsoValueGamma03PFId" + p))

    getattr(process, "pfElectrons" + p).isolationValueMapsCharged  = cms.VInputTag(cms.InputTag("elPFIsoValueCharged03PFId" + p))
    getattr(process, "pfElectrons" + p).deltaBetaIsolationValueMap = cms.InputTag("elPFIsoValuePU03PFId" + p)
    getattr(process, "pfElectrons" + p).isolationValueMapsNeutral  = cms.VInputTag(cms.InputTag( "elPFIsoValueNeutral03PFId" + p), cms.InputTag("elPFIsoValueGamma03PFId" + p))

    # ... And for PAT
    adaptPFIsoElectrons(process, getattr(process, "pfElectrons" + p), p, "03")


    # Setup quark gluon tagger
    #process.load('QuarkGluonTagger.EightTeV.QGTagger_RecoJets_cff')
    #cloneProcessingSnippet(process, process.QuarkGluonTagger, p)
    #getattr(process, "QGTagger" + p).srcJets = cms.InputTag("selectedPatJets" + p)
    #getattr(process, "QGTagger" + p).isPatJet = cms.untracked.bool(True)
    #getattr(process, "QGTagger" + p).useCHS = cms.untracked.bool(usePFNoPU)

    ## Remove the processing of primary vertices, as it's already what we do here
    #getattr(process, 'QGTagger' + p).srcPV = cms.InputTag('goodOfflinePrimaryVertices')
    #getattr(process, 'QuarkGluonTagger' + p).remove(getattr(process, 'goodOfflinePrimaryVerticesQG' + p))

    if not runOnMC:
      if 'L2L3Residual' in jetCorrections:
        getattr(process, 'patPFJetMETtype1p2Corr' + p).jetCorrLabel = 'L2L3Residual'
      getattr(process, 'patPFMet' + p).addGenMET = cms.bool(False)

    names = ["Taus"]
    if jetAlgo != "AK5":
      names += ["Electrons", "Muons"]
    if len(names) > 0:
      removeSpecificPATObjects(process, names = names, outputModules = ['out'], postfix = p)

    adaptPVs(process, pvCollection = cms.InputTag("goodOfflinePrimaryVertices"), postfix = p)

    getattr(process, "patDefaultSequence" + p).replace(getattr(process, "selectedPatElectrons" + p), getattr(process, "selectedPatElectrons" + p) + getattr(process, "patConversions" + p))

    #getattr(process, "patDefaultSequence" + p).replace(getattr(process, "selectedPatJets" + p), getattr(process, "selectedPatJets" + p) + getattr(process, "QuarkGluonTagger" + p))

    return getattr(process, "patPF2PATSequence" + p)

  # This must be called after all the calls to usePF2PAT
  def addC8Jets(p):
    p2 = p.replace('AK5', 'CA8')

    ###############################
    ###### Bare CA 0.8 jets #######
    ###############################
    setattr(process, "ca8PFJets" + p2, ca4PFJets.clone(
            rParam = cms.double(0.8),
            src = cms.InputTag('pfNoElectron' + p),
            doAreaFastjet = cms.bool(True),
            doRhoFastjet = cms.bool(True),
            Rho_EtaMax = cms.double(6.0),
            Ghost_EtaMax = cms.double(7.0),
            srcPVs = cms.InputTag("goodOfflinePrimaryVertices")
            ))

    # Add ca8PFJets to the main process path
    getattr(process, "patPF2PATSequence" + p).replace(
            getattr(process, "pfNoElectron" + p), getattr(process, "pfNoElectron" + p) * getattr(process, "ca8PFJets" + p2)
            )

    addJetCollection(process,
            cms.InputTag('ca8PFJets' + p2),
            'CA8', 'PFchs' if 'chs' in p else 'PF',
            doJTA = True,
            doBTagging = True,
            jetCorrLabel = ('AK7PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute']),
            doType1MET = True,
            doL1Cleaning = False,
            doL1Counters = False,
            genJetCollection = cms.InputTag("ca8GenJetsNoNu"),
            doJetID = True,
            )

  if correctMETWithT1:
    process.load("PhysicsTools.PatUtils.patPFMETCorrections_cff")
    #from PhysicsTools.PatAlgos.tools.helpers import cloneProcessingSnippet

  print "##########################"
  print "PF jets with PF2PAT"
  print "Using Type I met" if correctMETWithT1 else "NOT using Type I met"
  print "##########################"

  #runCHS = False

  #postfixes = {'PFlowAK5': 'AK5', 'PFlowAK7': 'AK7'}
  postfixes = {'PFlowAK5': 'AK5'}

  # Setup quark gluon tagger
  process.load('QuarkGluonTagger.EightTeV.QGTagger_RecoJets_cff')

  process.sequence_chs = cms.Sequence()
  process.sequence_nochs = cms.Sequence()
  for p, algo in postfixes.items():
    process.sequence_nochs += usePF2PATForAnalysis(jetAlgo=algo, postfix=p, usePFNoPU=False, useTypeIMET=correctMETWithT1)
    if runCHS:
      process.sequence_chs   += usePF2PATForAnalysis(jetAlgo=algo, postfix=p, usePFNoPU=True, useTypeIMET=correctMETWithT1)

    setattr(process, 'QGTagger' + p, process.QGTagger.clone())
    getattr(process, "QGTagger" + p).srcJets = cms.InputTag("selectedPatJets" + p)
    getattr(process, "QGTagger" + p).isPatJet = cms.untracked.bool(True)
    getattr(process, "QGTagger" + p).useCHS = cms.untracked.bool(False)

    process.QuarkGluonTagger.replace(process.QGTagger, getattr(process, 'QGTagger' + p))

    if runCHS:
      chsP = p + "chs"

      setattr(process, 'QGTagger' + chsP, process.QGTagger.clone())
      getattr(process, "QGTagger" + chsP).srcJets = cms.InputTag("selectedPatJets" + chsP)
      getattr(process, "QGTagger" + chsP).isPatJet = cms.untracked.bool(True)
      getattr(process, "QGTagger" + chsP).useCHS = cms.untracked.bool(True)

      process.QuarkGluonTagger.replace(getattr(process, 'QGTagger' + p), getattr(process, 'QGTagger' + p) + getattr(process, 'QGTagger' + chsP))

  for p, algo in postfixes.items():
    addC8Jets(p)
    if runCHS:
      addC8Jets(p + "chs")

  print "##########################"
  print "Calo jets" if processCaloJets else "No processing of calo jets"
  print "##########################"

  usePFIso(process, "")
  #adaptPFIsoPhotons(process,  process.patPhotons, "", "03")

  if processCaloJets:

    # No L2L3 Residual on purpose
    jetCorrections = ('AK5Calo', ['L1Offset', 'L2Relative', 'L3Absolute'])

    addJetCollection(process, cms.InputTag('ak7CaloJets'),
        'AK7',
        'Calo',
        doJTA            = True,
        doBTagging       = True,
        jetCorrLabel     = jetCorrections,
        doType1MET       = correctMETWithT1,
        doL1Cleaning     = False,
        doL1Counters     = False,
        genJetCollection = cms.InputTag("ak7GenJets"),
        doJetID          = True,
        jetIdLabel       = "ak7"
        )

    switchJetCollection(process, cms.InputTag('ak5CaloJets'),
        doJTA            = True,
        doBTagging       = True,
        jetCorrLabel     = jetCorrections,
        doType1MET       = correctMETWithT1,
        genJetCollection = cms.InputTag("ak5GenJets"),
        doJetID          = True,
        jetIdLabel       = "ak5"
        )

    process.selectedPatJets.cut = "pt > 2"
    process.selectedPatJetsAK7Calo.cut = "pt > 2"

  else:
    removeSpecificPATObjects(process, names = ['Jets', 'METs'])

  if not runOnMC:
    # Remove MC Matching
    removeMCMatching(process, names = ["All"])

  removeSpecificPATObjects(process, names = ['Electrons', 'Muons', 'Taus'])

  if runCHS:
    print "##########################"
    print "Running CHS sequence"
    print "##########################"

  process.analysisSequence = cms.Sequence()

  process.analysisSequence *= process.genParticlesForJetsNoNu * process.ca8GenJetsNoNu

  process.analysisSequence *= process.sequence_nochs
  if runCHS:
    process.analysisSequence *= process.sequence_chs

  # Quark Gluon tagging
  process.analysisSequence *= process.QuarkGluonTagger

  # Add default pat sequence to our path
  # This brings to life TcMET, Calo jets and Photons
  process.analysisSequence *= process.patDefaultSequence

  # Add our PhotonIsolationProducer to the analysisSequence. This producer compute pf isolations
  # for our photons
  process.photonPFIsolation = cms.EDProducer("PhotonIsolationProducer",
      src = cms.InputTag("selectedPatPhotons")
      )

  process.analysisSequence *= process.photonPFIsolation

  # Filtering

  # require physics declared
  process.load('HLTrigger.special.hltPhysicsDeclared_cfi')
  process.hltPhysicsDeclared.L1GtReadoutRecordTag = 'gtDigis'

  # require scraping filter
  process.scrapingVeto = cms.EDFilter("FilterOutScraping",
      applyfilter = cms.untracked.bool(True),
      debugOn = cms.untracked.bool(False),
      numtrack = cms.untracked.uint32(10),
      thresh = cms.untracked.double(0.25)
      )

  ## The iso-based HBHE noise filter ___________________________________________||
  process.load('CommonTools.RecoAlgos.HBHENoiseFilter_cfi')

  ## The CSC beam halo tight filter ____________________________________________||
  process.load('RecoMET.METAnalyzers.CSCHaloFilter_cfi')

  ## The HCAL laser filter _____________________________________________________||
  process.load("RecoMET.METFilters.hcalLaserEventFilter_cfi")

  ## The ECAL dead cell trigger primitive filter _______________________________||
  process.load('RecoMET.METFilters.EcalDeadCellTriggerPrimitiveFilter_cfi')

  ## The EE bad SuperCrystal filter ____________________________________________||
  process.load('RecoMET.METFilters.eeBadScFilter_cfi')

  ## The ECAL laser correction filter
  process.load('RecoMET.METFilters.ecalLaserCorrFilter_cfi')

  ## The Good vertices collection needed by the tracking failure filter ________||
  process.goodVertices = cms.EDFilter(
      "VertexSelector",
      filter = cms.bool(False),
      src = cms.InputTag("offlinePrimaryVertices"),
      cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.rho < 2")
      )

  ## The tracking failure filter _______________________________________________||
  process.load('RecoMET.METFilters.trackingFailureFilter_cfi')

  ## The tracking POG filters __________________________________________________||
  process.load('RecoMET.METFilters.trackingPOGFilters_cff')

  # Count events
  process.nEventsTotal    = cms.EDProducer("EventCountProducer")
  process.nEventsFiltered = cms.EDProducer("EventCountProducer")

  # Let it run
  process.p = cms.Path(
      process.nEventsTotal +

      # Filters
      process.primaryVertexFilter +
      process.scrapingVeto +
      process.HBHENoiseFilter +
      process.CSCTightHaloFilter +
      process.hcalLaserEventFilter +
      process.EcalDeadCellTriggerPrimitiveFilter +
      process.goodVertices + process.trackingFailureFilter +
      process.eeBadScFilter +
      process.ecalLaserCorrFilter +
      process.trkPOGFilters +

      process.goodOfflinePrimaryVertices +

      # Physics
      process.analysisSequence +

      process.nEventsFiltered
      )

  if runOnMC:
    process.p.remove(process.scrapingVeto)

  # Add PF2PAT output to the created file
  from PhysicsTools.PatAlgos.patEventContent_cff import patEventContentNoCleaning
  #process.load("CommonTools.ParticleFlow.PF2PAT_EventContent_cff")
  process.out.outputCommands = cms.untracked.vstring('drop *',
      'keep *_photonCore_*_*',
      'keep double_kt6*Jets*_rho_*',
      'keep *_goodOfflinePrimaryVertices_*_*',
      'keep recoPFCandidates_particleFlow_*_*',
      # Content of *patEventContentNoCleaning
      'keep *_selectedPatPhotons*_*_*', 'keep *_selectedPatElectrons*_*_*', 'keep *_selectedPatMuons*_*_*', 'keep *_selectedPatTaus*_*_*', 'keep *_selectedPatJets*_*_*', 'drop *_selectedPatJets_pfCandidates_*', 'drop *_*PF_caloTowers_*', 'drop *_*JPT_pfCandidates_*', 'drop *_*Calo_pfCandidates_*', 'keep *_patMETs*_*_*', 'keep *_selectedPatPFParticles*_*_*', 'keep *_selectedPatTrackCands*_*_*',
      'keep *_cleanPatPhotons*_*_*', 'keep *_cleanPatElectrons*_*_*', 'keep *_cleanPatMuons*_*_*', 'keep *_cleanPatTaus*_*_*',
      'keep *_cleanPatJets*_*_*', 'keep *_cleanPatHemispheres*_*_*', 'keep *_cleanPatPFParticles*_*_*',
      'keep *_cleanPatTrackCands*_*_*',
      'drop *_*PFlow_caloTowers_*',
      # Type I residual
      'drop *_selectedPatJetsForMET*_*_PAT',
      'keep *_patPFMet*_*_PAT', # Keep raw met
      # Trigger
      'keep *_TriggerResults_*_HLT',
      # Debug
      #'keep *_pf*_*_PAT'
      # Photon ID
      'keep *_patConversions*_*_*',
      'keep *_photonPFIsolation*_*_*',
      # Quark Gluon tagging
      'keep *_QGTagger*_*_*',
      'drop *_kt6PFJetsIsoQG_*_PAT',
      'drop *_kt6PFJetsQG_*_PAT',
      # MC truth
      'keep *_genParticles_*_*'
      )

  if runOnMC:
    process.out.outputCommands.extend(['keep *_addPileupInfo_*_*', 'keep *_generator_*_*'])

  # switch on PAT trigger
  # from PhysicsTools.PatAlgos.tools.trigTools import switchOnTrigger
  # switchOnTrigger( process )

  ## ------------------------------------------------------
  #  In addition you usually want to change the following
  #  parameters:
  ## ------------------------------------------------------
  #
  process.GlobalTag.globaltag = "%s::All" % (globalTag) ##  (according to https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideFrontierConditions)
  #                                         ##
  #process.source.fileNames =  cms.untracked.vstring('file:input_data.root')  ##  (e.g. 'file:AOD.root')
  #                                         ##
  process.maxEvents.input = 2500
  #                                         ##
  #   process.out.outputCommands = [ ... ]  ##  (e.g. taken from PhysicsTools/PatAlgos/python/patEventContent_cff.py)
  #                                         ##
  process.options.wantSummary = False   ##  (to suppress the long output at the end of the job)
  process.MessageLogger.cerr.FwkReport.reportEvery = 1000

  # Remove annoying ecalLaser messages
  process.MessageLogger.suppressError = cms.untracked.vstring ('ecalLaserCorrFilter')

  return process
	diff --git a/crab/produce_PAT_COMMON.py b/crab/produce_PAT_COMMON.py
	index 082490f..8e3f4e1 100644
	--- a/crab/produce_PAT_COMMON.py
	+++ b/crab/produce_PAT_COMMON.py
	@@ -40,8 +40,18 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)
	# this function will modify the PAT sequences.
	from PhysicsTools.PatAlgos.tools.coreTools import removeSpecificPATObjects, removeMCMatching
	from PhysicsTools.PatAlgos.tools.pfTools import usePF2PAT, adaptPFIsoElectrons, adaptPVs, usePFIso
	+ from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection
	#from PhysicsTools.PatAlgos.tools.metTools import *

	+ process.load("RecoJets.Configuration.GenJetParticles_cff")
	+ from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
	+ process.ca8GenJetsNoNu = ca4GenJets.clone(
	+ rParam = cms.double(0.8),
	+ src = cms.InputTag("genParticlesForJetsNoNu")
	+ )
	+
	+ from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets
	+
	def usePF2PATForAnalysis(jetAlgo, postfix, useTypeIMET, usePFNoPU):

	# Jet corrections
	@@ -137,6 +147,41 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)

	return getattr(process, "patPF2PATSequence" + p)

	+ # This must be called after all the calls to usePF2PAT
	+ def addC8Jets(p):
	+ p2 = p.replace('AK5', 'CA8')
	+
	+ ###############################
	+ ###### Bare CA 0.8 jets #######
	+ ###############################
	+ setattr(process, "ca8PFJets" + p2, ca4PFJets.clone(
	+ rParam = cms.double(0.8),
	+ src = cms.InputTag('pfNoElectron' + p),
	+ doAreaFastjet = cms.bool(True),
	+ doRhoFastjet = cms.bool(True),
	+ Rho_EtaMax = cms.double(6.0),
	+ Ghost_EtaMax = cms.double(7.0),
	+ srcPVs = cms.InputTag("goodOfflinePrimaryVertices")
	+ ))
	+
	+ # Add ca8PFJets to the main process path
	+ getattr(process, "patPF2PATSequence" + p).replace(
	+ getattr(process, "pfNoElectron" + p), getattr(process, "pfNoElectron" + p) * getattr(process, "ca8PFJets" + p2)
	+ )
	+
	+ addJetCollection(process,
	+ cms.InputTag('ca8PFJets' + p2),
	+ 'CA8', 'PFchs' if 'chs' in p else 'PF',
	+ doJTA = True,
	+ doBTagging = True,
	+ jetCorrLabel = ('AK7PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute']),
	+ doType1MET = True,
	+ doL1Cleaning = False,
	+ doL1Counters = False,
	+ genJetCollection = cms.InputTag("ca8GenJetsNoNu"),
	+ doJetID = True,
	+ )
	+
	if correctMETWithT1:
	process.load("PhysicsTools.PatUtils.patPFMETCorrections_cff")
	#from PhysicsTools.PatAlgos.tools.helpers import cloneProcessingSnippet
	@@ -178,6 +223,11 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)

	process.QuarkGluonTagger.replace(getattr(process, 'QGTagger' + p), getattr(process, 'QGTagger' + p) + getattr(process, 'QGTagger' + chsP))

	+ for p, algo in postfixes.items():
	+ addC8Jets(p)
	+ if runCHS:
	+ addC8Jets(p + "chs")
	+
	print "##########################"
	print "Calo jets" if processCaloJets else "No processing of calo jets"
	print "##########################"
	@@ -233,6 +283,8 @@ def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag)

	process.analysisSequence = cms.Sequence()

	+ process.analysisSequence = process.genParticlesForJetsNoNu process.ca8GenJetsNoNu
	+
	process.analysisSequence *= process.sequence_nochs
	if runCHS:
	process.analysisSequence *= process.sequence_chs
	import FWCore.ParameterSet.Config as cms

	def createProcess(runOnMC, runCHS, correctMETWithT1, processCaloJets, globalTag):

	## import skeleton process
	from PhysicsTools.PatAlgos.patTemplate_cfg import process

	# load the PAT config
	process.load("PhysicsTools.PatAlgos.patSequences_cff")
	process.load("JetMETCorrections.Configuration.JetCorrectionProducers_cff")

	# Do some CHS stuff
	process.ak5PFchsL1Fastjet = process.ak5PFL1Fastjet.clone(algorithm = 'AK5PFchs')
	process.ak5PFchsL2Relative = process.ak5PFL2Relative.clone(algorithm = 'AK5PFchs')
	process.ak5PFchsL3Absolute = process.ak5PFL3Absolute.clone(algorithm = 'AK5PFchs')
	process.ak5PFchsResidual = process.ak5PFResidual.clone(algorithm = 'AK5PFchs')
	process.ak5PFchsL1FastL2L3 = cms.ESProducer(
	'JetCorrectionESChain',
	correctors = cms.vstring('ak5PFchsL1Fastjet', 'ak5PFchsL2Relative','ak5PFchsL3Absolute')
	)


	from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector

	## The good primary vertex filter ____________________________________________\|\|
	## This filter throw events with no primary vertex
	process.primaryVertexFilter = cms.EDFilter(
	"VertexSelector",
	src = cms.InputTag("offlinePrimaryVertices"),
	cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.Rho <= 2"),
	filter = cms.bool(True)
	)

	process.goodOfflinePrimaryVertices = cms.EDFilter("PrimaryVertexObjectFilter",
	filterParams = pvSelector.clone( minNdof = cms.double(4.0), maxZ = cms.double(24.0) ),
	src = cms.InputTag('offlinePrimaryVertices')
	)

	# Configure PAT to use PF2PAT instead of AOD sources
	# this function will modify the PAT sequences.
	from PhysicsTools.PatAlgos.tools.coreTools import removeSpecificPATObjects, removeMCMatching
	from PhysicsTools.PatAlgos.tools.pfTools import usePF2PAT, adaptPFIsoElectrons, adaptPVs, usePFIso
	from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection
	#from PhysicsTools.PatAlgos.tools.metTools import *

	process.load("RecoJets.Configuration.GenJetParticles_cff")
	from RecoJets.JetProducers.ca4GenJets_cfi import ca4GenJets
	process.ca8GenJetsNoNu = ca4GenJets.clone(
	rParam = cms.double(0.8),
	src = cms.InputTag("genParticlesForJetsNoNu")
	)

	from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets

	def usePF2PATForAnalysis(jetAlgo, postfix, useTypeIMET, usePFNoPU):

	# Jet corrections
	# No L2L3 Residual on purpose
	if usePFNoPU:
	jetCorrections = ("%sPFchs" % algo, ['L1FastJet', 'L2Relative', 'L3Absolute'])
	postfix += "chs"
	else:
	jetCorrections = ("%sPF" % algo, ['L1FastJet', 'L2Relative', 'L3Absolute'])

	#if not runOnMC:
	# jetCorrections[1].append('L2L3Residual')

	p = postfix

	usePF2PAT(process, runPF2PAT=True, jetAlgo=jetAlgo, runOnMC=runOnMC, postfix=p, jetCorrections=jetCorrections, typeIMetCorrections = useTypeIMET)
	getattr(process, "pfPileUp" + p).Enable = True
	getattr(process, "pfPileUp" + p).Vertices = 'goodOfflinePrimaryVertices'
	getattr(process, "pfPileUp" + p).checkClosestZVertex = cms.bool(False)

	getattr(process, "pfJets" + p).doAreaFastjet = cms.bool(True)
	getattr(process, "pfJets" + p).doRhoFastjet = False
	getattr(process, 'patJetCorrFactors' + p).rho = cms.InputTag("kt6PFJets", "rho") # Do not use kt6PFJetsPFlowAK5, it's not ok for L1FastJet.

	# top projections in PF2PAT:
	getattr(process,"pfNoPileUp" + p).enable = cms.bool(usePFNoPU)
	getattr(process,"pfNoMuon" + p).enable = True
	getattr(process,"pfNoElectron" + p).enable = True
	getattr(process,"pfNoTau" + p).enable = False
	getattr(process,"pfNoJet" + p).enable = True

	getattr(process,"patElectrons" + p).embedTrack = True

	getattr(process,"patMuons" + p).embedTrack = True
	# enable delta beta correction for muon selection in PF2PAT?
	getattr(process,"pfIsolatedMuons" + p).doDeltaBetaCorrection = True

	getattr(process, "patJets" + p).embedPFCandidates = False
	# Keep only jets with pt > 2 Gev
	getattr(process, "selectedPatJets" + p).cut = "pt > 2";

	# Use a cone of 0.3 for photon isolation
	#adaptPFIsoPhotons(process, applyPostfix(process, "patPhotons", postfix), postfix, "03")

	# 2012 Photon ID

	# Electron conversion
	setattr(process, "patConversions" + p, cms.EDProducer("PATConversionProducer",
	# input collection
	electronSource = cms.InputTag("selectedPatElectrons" + p)
	))

	# Switch electron isolation to dR = 0.3, for PF2PAT
	getattr(process, "pfIsolatedElectrons" + p).isolationValueMapsCharged = cms.VInputTag(cms.InputTag("elPFIsoValueCharged03PFId" + p))
	getattr(process, "pfIsolatedElectrons" + p).deltaBetaIsolationValueMap = cms.InputTag("elPFIsoValuePU03PFId" + p)
	getattr(process, "pfIsolatedElectrons" + p).isolationValueMapsNeutral = cms.VInputTag(cms.InputTag("elPFIsoValueNeutral03PFId" + p), cms.InputTag("elPFIsoValueGamma03PFId" + p))

	getattr(process, "pfElectrons" + p).isolationValueMapsCharged = cms.VInputTag(cms.InputTag("elPFIsoValueCharged03PFId" + p))
	getattr(process, "pfElectrons" + p).deltaBetaIsolationValueMap = cms.InputTag("elPFIsoValuePU03PFId" + p)
	getattr(process, "pfElectrons" + p).isolationValueMapsNeutral = cms.VInputTag(cms.InputTag( "elPFIsoValueNeutral03PFId" + p), cms.InputTag("elPFIsoValueGamma03PFId" + p))

	# ... And for PAT
	adaptPFIsoElectrons(process, getattr(process, "pfElectrons" + p), p, "03")


	# Setup quark gluon tagger
	#process.load('QuarkGluonTagger.EightTeV.QGTagger_RecoJets_cff')
	#cloneProcessingSnippet(process, process.QuarkGluonTagger, p)
	#getattr(process, "QGTagger" + p).srcJets = cms.InputTag("selectedPatJets" + p)
	#getattr(process, "QGTagger" + p).isPatJet = cms.untracked.bool(True)
	#getattr(process, "QGTagger" + p).useCHS = cms.untracked.bool(usePFNoPU)

	## Remove the processing of primary vertices, as it's already what we do here
	#getattr(process, 'QGTagger' + p).srcPV = cms.InputTag('goodOfflinePrimaryVertices')
	#getattr(process, 'QuarkGluonTagger' + p).remove(getattr(process, 'goodOfflinePrimaryVerticesQG' + p))

	if not runOnMC:
	if 'L2L3Residual' in jetCorrections:
	getattr(process, 'patPFJetMETtype1p2Corr' + p).jetCorrLabel = 'L2L3Residual'
	getattr(process, 'patPFMet' + p).addGenMET = cms.bool(False)

	names = ["Taus"]
	if jetAlgo != "AK5":
	names += ["Electrons", "Muons"]
	if len(names) > 0:
	removeSpecificPATObjects(process, names = names, outputModules = ['out'], postfix = p)

	adaptPVs(process, pvCollection = cms.InputTag("goodOfflinePrimaryVertices"), postfix = p)

	getattr(process, "patDefaultSequence" + p).replace(getattr(process, "selectedPatElectrons" + p), getattr(process, "selectedPatElectrons" + p) + getattr(process, "patConversions" + p))

	#getattr(process, "patDefaultSequence" + p).replace(getattr(process, "selectedPatJets" + p), getattr(process, "selectedPatJets" + p) + getattr(process, "QuarkGluonTagger" + p))

	return getattr(process, "patPF2PATSequence" + p)

	# This must be called after all the calls to usePF2PAT
	def addC8Jets(p):
	p2 = p.replace('AK5', 'CA8')

	###############################
	###### Bare CA 0.8 jets #######
	###############################
	setattr(process, "ca8PFJets" + p2, ca4PFJets.clone(
	rParam = cms.double(0.8),
	src = cms.InputTag('pfNoElectron' + p),
	doAreaFastjet = cms.bool(True),
	doRhoFastjet = cms.bool(True),
	Rho_EtaMax = cms.double(6.0),
	Ghost_EtaMax = cms.double(7.0),
	srcPVs = cms.InputTag("goodOfflinePrimaryVertices")
	))

	# Add ca8PFJets to the main process path
	getattr(process, "patPF2PATSequence" + p).replace(
	getattr(process, "pfNoElectron" + p), getattr(process, "pfNoElectron" + p) * getattr(process, "ca8PFJets" + p2)
	)

	addJetCollection(process,
	cms.InputTag('ca8PFJets' + p2),
	'CA8', 'PFchs' if 'chs' in p else 'PF',
	doJTA = True,
	doBTagging = True,
	jetCorrLabel = ('AK7PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute']),
	doType1MET = True,
	doL1Cleaning = False,
	doL1Counters = False,
	genJetCollection = cms.InputTag("ca8GenJetsNoNu"),
	doJetID = True,
	)

	if correctMETWithT1:
	process.load("PhysicsTools.PatUtils.patPFMETCorrections_cff")
	#from PhysicsTools.PatAlgos.tools.helpers import cloneProcessingSnippet

	print "##########################"
	print "PF jets with PF2PAT"
	print "Using Type I met" if correctMETWithT1 else "NOT using Type I met"
	print "##########################"

	#runCHS = False

	#postfixes = {'PFlowAK5': 'AK5', 'PFlowAK7': 'AK7'}
	postfixes = {'PFlowAK5': 'AK5'}

	# Setup quark gluon tagger
	process.load('QuarkGluonTagger.EightTeV.QGTagger_RecoJets_cff')

	process.sequence_chs = cms.Sequence()
	process.sequence_nochs = cms.Sequence()
	for p, algo in postfixes.items():
	process.sequence_nochs += usePF2PATForAnalysis(jetAlgo=algo, postfix=p, usePFNoPU=False, useTypeIMET=correctMETWithT1)
	if runCHS:
	process.sequence_chs += usePF2PATForAnalysis(jetAlgo=algo, postfix=p, usePFNoPU=True, useTypeIMET=correctMETWithT1)

	setattr(process, 'QGTagger' + p, process.QGTagger.clone())
	getattr(process, "QGTagger" + p).srcJets = cms.InputTag("selectedPatJets" + p)
	getattr(process, "QGTagger" + p).isPatJet = cms.untracked.bool(True)
	getattr(process, "QGTagger" + p).useCHS = cms.untracked.bool(False)

	process.QuarkGluonTagger.replace(process.QGTagger, getattr(process, 'QGTagger' + p))

	if runCHS:
	chsP = p + "chs"

	setattr(process, 'QGTagger' + chsP, process.QGTagger.clone())
	getattr(process, "QGTagger" + chsP).srcJets = cms.InputTag("selectedPatJets" + chsP)
	getattr(process, "QGTagger" + chsP).isPatJet = cms.untracked.bool(True)
	getattr(process, "QGTagger" + chsP).useCHS = cms.untracked.bool(True)

	process.QuarkGluonTagger.replace(getattr(process, 'QGTagger' + p), getattr(process, 'QGTagger' + p) + getattr(process, 'QGTagger' + chsP))

	for p, algo in postfixes.items():
	addC8Jets(p)
	if runCHS:
	addC8Jets(p + "chs")

	print "##########################"
	print "Calo jets" if processCaloJets else "No processing of calo jets"
	print "##########################"

	usePFIso(process, "")
	#adaptPFIsoPhotons(process, process.patPhotons, "", "03")

	if processCaloJets:

	# No L2L3 Residual on purpose
	jetCorrections = ('AK5Calo', ['L1Offset', 'L2Relative', 'L3Absolute'])

	addJetCollection(process, cms.InputTag('ak7CaloJets'),
	'AK7',
	'Calo',
	doJTA = True,
	doBTagging = True,
	jetCorrLabel = jetCorrections,
	doType1MET = correctMETWithT1,
	doL1Cleaning = False,
	doL1Counters = False,
	genJetCollection = cms.InputTag("ak7GenJets"),
	doJetID = True,
	jetIdLabel = "ak7"
	)

	switchJetCollection(process, cms.InputTag('ak5CaloJets'),
	doJTA = True,
	doBTagging = True,
	jetCorrLabel = jetCorrections,
	doType1MET = correctMETWithT1,
	genJetCollection = cms.InputTag("ak5GenJets"),
	doJetID = True,
	jetIdLabel = "ak5"
	)

	process.selectedPatJets.cut = "pt > 2"
	process.selectedPatJetsAK7Calo.cut = "pt > 2"

	else:
	removeSpecificPATObjects(process, names = ['Jets', 'METs'])

	if not runOnMC:
	# Remove MC Matching
	removeMCMatching(process, names = ["All"])

	removeSpecificPATObjects(process, names = ['Electrons', 'Muons', 'Taus'])

	if runCHS:
	print "##########################"
	print "Running CHS sequence"
	print "##########################"

	process.analysisSequence = cms.Sequence()

	process.analysisSequence = process.genParticlesForJetsNoNu process.ca8GenJetsNoNu

	process.analysisSequence *= process.sequence_nochs
	if runCHS:
	process.analysisSequence *= process.sequence_chs

	# Quark Gluon tagging
	process.analysisSequence *= process.QuarkGluonTagger

	# Add default pat sequence to our path
	# This brings to life TcMET, Calo jets and Photons
	process.analysisSequence *= process.patDefaultSequence

	# Add our PhotonIsolationProducer to the analysisSequence. This producer compute pf isolations
	# for our photons
	process.photonPFIsolation = cms.EDProducer("PhotonIsolationProducer",
	src = cms.InputTag("selectedPatPhotons")
	)

	process.analysisSequence *= process.photonPFIsolation

	# Filtering

	# require physics declared
	process.load('HLTrigger.special.hltPhysicsDeclared_cfi')
	process.hltPhysicsDeclared.L1GtReadoutRecordTag = 'gtDigis'

	# require scraping filter
	process.scrapingVeto = cms.EDFilter("FilterOutScraping",
	applyfilter = cms.untracked.bool(True),
	debugOn = cms.untracked.bool(False),
	numtrack = cms.untracked.uint32(10),
	thresh = cms.untracked.double(0.25)
	)

	## The iso-based HBHE noise filter ___________________________________________\|\|
	process.load('CommonTools.RecoAlgos.HBHENoiseFilter_cfi')

	## The CSC beam halo tight filter ____________________________________________\|\|
	process.load('RecoMET.METAnalyzers.CSCHaloFilter_cfi')

	## The HCAL laser filter _____________________________________________________\|\|
	process.load("RecoMET.METFilters.hcalLaserEventFilter_cfi")

	## The ECAL dead cell trigger primitive filter _______________________________\|\|
	process.load('RecoMET.METFilters.EcalDeadCellTriggerPrimitiveFilter_cfi')

	## The EE bad SuperCrystal filter ____________________________________________\|\|
	process.load('RecoMET.METFilters.eeBadScFilter_cfi')

	## The ECAL laser correction filter
	process.load('RecoMET.METFilters.ecalLaserCorrFilter_cfi')

	## The Good vertices collection needed by the tracking failure filter ________\|\|
	process.goodVertices = cms.EDFilter(
	"VertexSelector",
	filter = cms.bool(False),
	src = cms.InputTag("offlinePrimaryVertices"),
	cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.rho < 2")
	)

	## The tracking failure filter _______________________________________________\|\|
	process.load('RecoMET.METFilters.trackingFailureFilter_cfi')

	## The tracking POG filters __________________________________________________\|\|
	process.load('RecoMET.METFilters.trackingPOGFilters_cff')

	# Count events
	process.nEventsTotal = cms.EDProducer("EventCountProducer")
	process.nEventsFiltered = cms.EDProducer("EventCountProducer")

	# Let it run
	process.p = cms.Path(
	process.nEventsTotal +

	# Filters
	process.primaryVertexFilter +
	process.scrapingVeto +
	process.HBHENoiseFilter +
	process.CSCTightHaloFilter +
	process.hcalLaserEventFilter +
	process.EcalDeadCellTriggerPrimitiveFilter +
	process.goodVertices + process.trackingFailureFilter +
	process.eeBadScFilter +
	process.ecalLaserCorrFilter +
	process.trkPOGFilters +

	process.goodOfflinePrimaryVertices +

	# Physics
	process.analysisSequence +

	process.nEventsFiltered
	)

	if runOnMC:
	process.p.remove(process.scrapingVeto)

	# Add PF2PAT output to the created file
	from PhysicsTools.PatAlgos.patEventContent_cff import patEventContentNoCleaning
	#process.load("CommonTools.ParticleFlow.PF2PAT_EventContent_cff")
	process.out.outputCommands = cms.untracked.vstring('drop *',
	'keep _photonCore__*',
	'keep double_kt6Jets_rho_*',
	'keep _goodOfflinePrimaryVertices__*',
	'keep recoPFCandidates_particleFlow__',
	# Content of *patEventContentNoCleaning
	'keep _selectedPatPhotons__', 'keep _selectedPatElectrons__', 'keep _selectedPatMuons__', 'keep _selectedPatTaus__', 'keep _selectedPatJets__', 'drop _selectedPatJets_pfCandidates_', 'drop _PF_caloTowers_', 'drop _JPT_pfCandidates_', 'drop _Calo_pfCandidates_', 'keep _patMETs__', 'keep _selectedPatPFParticles__', 'keep _selectedPatTrackCands__*',
	'keep _cleanPatPhotons__', 'keep _cleanPatElectrons__', 'keep _cleanPatMuons__', 'keep _cleanPatTaus__',
	'keep _cleanPatJets__', 'keep _cleanPatHemispheres__', 'keep _cleanPatPFParticles__',
	'keep _cleanPatTrackCands__',
	'drop _PFlow_caloTowers_*',
	# Type I residual
	'drop _selectedPatJetsForMET_*_PAT',
	'keep _patPFMet_*_PAT', # Keep raw met
	# Trigger
	'keep _TriggerResults__HLT',
	# Debug
	#'keep _pf_*_PAT'
	# Photon ID
	'keep _patConversions__',
	'keep _photonPFIsolation__',
	# Quark Gluon tagging
	'keep _QGTagger__',
	'drop _kt6PFJetsIsoQG__PAT',
	'drop _kt6PFJetsQG__PAT',
	# MC truth
	'keep _genParticles__*'
	)

	if runOnMC:
	process.out.outputCommands.extend(['keep _addPileupInfo__', 'keep _generator__'])

	# switch on PAT trigger
	# from PhysicsTools.PatAlgos.tools.trigTools import switchOnTrigger
	# switchOnTrigger( process )

	## ------------------------------------------------------
	# In addition you usually want to change the following
	# parameters:
	## ------------------------------------------------------
	#
	process.GlobalTag.globaltag = "%s::All" % (globalTag) ## (according to https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideFrontierConditions)
	# ##
	#process.source.fileNames = cms.untracked.vstring('file:input_data.root') ## (e.g. 'file:AOD.root')
	# ##
	process.maxEvents.input = 2500
	# ##
	# process.out.outputCommands = [ ... ] ## (e.g. taken from PhysicsTools/PatAlgos/python/patEventContent_cff.py)
	# ##
	process.options.wantSummary = False ## (to suppress the long output at the end of the job)
	process.MessageLogger.cerr.FwkReport.reportEvery = 1000

	# Remove annoying ecalLaser messages
	process.MessageLogger.suppressError = cms.untracked.vstring ('ecalLaserCorrFilter')

	return process