HiEgammaAlgos/test/testSpike_cfg.py

0001
0002 import FWCore.ParameterSet.Config as cms
0003
0004 process = cms.Process("PAT")
0005
0006 ## MessageLogger
0007 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0008 process.Timing = cms.Service("Timing")
0009
0010 ## Options and Output Report
0011 process.options   = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0012
0013 ## Source
0014 process.source = cms.Source("PoolSource",
0015                             fileNames = cms.untracked.vstring(
0016     # 'rfio:/castor/cern.ch/cms/store/relval/CMSSW_3_6_3/RelValPhotonJets_Pt_10/GEN-SIM-RECO/START36_V10-v1/0005/780848E8-0978-DF11-9D44-00261894389C.root'  # pp
0017     'rfio:/castor/cern.ch/cms/store/relval/CMSSW_3_9_1/RelValHydjetQ_MinBias_2760GeV/GEN-SIM-RECO/MC_39Y_V3-v1/0062/E062CB6E-49E4-DF11-B10A-00261894397D.root'
0018     ),
0019                             duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
0020                             )
0021
0022 ## Maximal Number of Events
0023 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0024
0025 ## Geometry and Detector Conditions (needed for a few patTuple production steps)
0026 process.load("Configuration.StandardSequences.Geometry_cff")
0027 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0028 process.load('Configuration.StandardSequences.Reconstruction_cff')
0029 process.load("RecoHI.HiEgammaAlgos.HiEgamma_cff")
0030
0031 process.GlobalTag.globaltag = cms.string('MC_39Y_V3::All')
0032
0033 process.load("Configuration.StandardSequences.MagneticField_cff")
0034
0035 ## Standard PAT Configuration File
0036 #process.load("PhysicsTools.PatAlgos.patSequences_cff")
0037 process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
0038 from PhysicsTools.PatAlgos.tools.heavyIonTools import *
0039 configureHeavyIons(process)
0040
0041 from PhysicsTools.PatAlgos.tools.coreTools import *
0042
0043 process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(100)
0044
0045
0046
0047 # Centrality
0048 process.load("RecoHI.HiCentralityAlgos.HiCentrality_cfi")
0049
0050 # Heavyion
0051 process.GlobalTag.toGet = cms.VPSet(
0052     cms.PSet(record = cms.string("HeavyIonRcd"),
0053              tag = cms.string("CentralityTable_HFhits40_AMPT2760GeV_v3_mc"),
0054              connect = cms.untracked.string("frontier://FrontierPrep/CMS_COND_PHYSICSTOOLS")
0055              )
0056     )
0057
0058 #ecal filter
0059 process.goodPhotons = cms.EDFilter("PhotonSelector",
0060                            src = cms.InputTag("photons"),
0061                            cut = cms.string('et > 40 && hadronicOverEm < 0.1 && r9 > 0.8')
0062                            )
0063 # leading photon E_T filter
0064 process.photonFilter = cms.EDFilter("EtMinPhotonCountFilter",
0065                             src = cms.InputTag("goodPhotons"),
0066                             etMin = cms.double(40.0),
0067                             minNumber = cms.uint32(1)
0068                             )
0069
0070 process.TFileService = cms.Service("TFileService",
0071                                    fileName = cms.string('isoRecHitTree.root'),
0072                                    closeFileFast = cms.untracked.bool(True)
0073                                    )
0074
0075 #process.load("RecoHI.HiEgammaAlgos.hiEcalSpikeFilter_cfi")
0076
0077 process.load("CmsHi.PhotonAnalysis.isoConeInspector_cfi")
0078 process.load("CmsHi.PhotonAnalysis.ecalHistProducer_cfi")
0079
0080 process.isoConeMap.etCut = 5;  # For hydjet
0081
0082 process.load("RecoHI.HiEgammaAlgos.hiSpikeCleaner_cfi")
0083
0084 process.output = cms.OutputModule("PoolOutputModule",
0085                                   splitLevel = cms.untracked.int32(0),
0086                                   outputCommands = cms.untracked.vstring('keep *_*_*_*'),
0087                                   fileName = cms.untracked.string('edmOut.root')
0088                                   )
0089
0090
0091 process.p = cms.Path(
0092     #   process.goodPhotons *
0093     #    process.photonFilter*
0094     # process.hiEcalSpikeFilter
0095     # process.isoConeMap  *
0096     #  process.ecalHistProducer
0097     process.hiSpikeCleaner*
0098     process.output
0099     )
0100
0101
0102 #process.e = cms.EndPath(process.out)
0103
0104