Configuration/test/centralSkimsHI_SKIM.py

0001 # Auto generated configuration file
0002 # using:
0003 # Revision: 1.232.2.6
0004 # Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v
0005 # with command line options: centralSkimsHI -s SKIM:DiJet+Photon+ZMM+ZEE --conditions GR10_P_V12::All --scenario HeavyIons --filein=/store/hidata/HIRun2010/HIAllPhysics/RECO/PromptReco-v1/000/150/063/B497BEDB-8BE8-DF11-B09D-0030487A18F2.root --data --no_exec
0006 import FWCore.ParameterSet.Config as cms
0007
0008 process = cms.Process('SKIM')
0009
0010 # import of standard configurations
0011 process.load('Configuration.StandardSequences.Services_cff')
0012 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0013 process.load('FWCore.MessageService.MessageLogger_cfi')
0014 process.load('Configuration.StandardSequences.GeometryDB_cff')
0015 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0016 process.load('Configuration.StandardSequences.SkimsHeavyIons_cff')
0017 process.load('Configuration.StandardSequences.ReconstructionHeavyIons_cff')
0018 process.load('RecoHI.HiEgammaAlgos.HiElectronSequence_cff')
0019 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0020 process.load('Configuration.EventContent.EventContentHeavyIons_cff')
0021
0022 process.configurationMetadata = cms.untracked.PSet(
0023     version = cms.untracked.string('$Revision: 1.6 $'),
0024     annotation = cms.untracked.string('centralSkimsHI nevts:1'),
0025     name = cms.untracked.string('PyReleaseValidation')
0026 )
0027 process.maxEvents = cms.untracked.PSet(
0028     input = cms.untracked.int32(-1)
0029 )
0030
0031 # Input source
0032 process.source = cms.Source("PoolSource",
0033     fileNames = cms.untracked.vstring(
0034     '/store/hidata/HIRun2010/HICorePhysics/RECO/PromptReco-v3/000/151/353/DC04F2F3-32F2-DF11-BD72-0030487C5CE2.root'),
0035     secondaryFileNames = cms.untracked.vstring(
0036     '/store/hidata/HIRun2010/HICorePhysics/RAW/v1/000/151/353/00B95950-00F2-DF11-9485-001D09F252DA.root')
0037 )
0038
0039 process.options = cms.untracked.PSet(
0040     #wantSummary = cms.untracked.bool(True)
0041 )
0042
0043 #process.Timing = cms.Service("Timing")
0044
0045 # Output definition
0046
0047 # Additional output definition
0048 process.SKIMStreamDiJet = cms.OutputModule("PoolOutputModule",
0049     process.FEVTEventContent,
0050     SelectEvents = cms.untracked.PSet(
0051         SelectEvents = cms.vstring('diJetSkimPath')
0052     ),
0053     fileName = cms.untracked.string('DiJet.root'),
0054     dataset = cms.untracked.PSet(
0055         filterName = cms.untracked.string('DiJet'),
0056         dataTier = cms.untracked.string('RAW-RECO')
0057     )
0058 )
0059 process.SKIMStreamPhoton = cms.OutputModule("PoolOutputModule",
0060     process.FEVTEventContent,
0061     SelectEvents = cms.untracked.PSet(
0062         SelectEvents = cms.vstring('photonSkimPath')
0063     ),
0064     fileName = cms.untracked.string('Photon.root'),
0065     dataset = cms.untracked.PSet(
0066         filterName = cms.untracked.string('Photon'),
0067         dataTier = cms.untracked.string('RAW-RECO')
0068     )
0069 )
0070 process.SKIMStreamZEE = cms.OutputModule("PoolOutputModule",
0071     process.FEVTEventContent,
0072     SelectEvents = cms.untracked.PSet(
0073         SelectEvents = cms.vstring('zEESkimPath')
0074     ),
0075     fileName = cms.untracked.string('ZEE.root'),
0076     dataset = cms.untracked.PSet(
0077         filterName = cms.untracked.string('ZEE'),
0078         dataTier = cms.untracked.string('RAW-RECO')
0079     )
0080 )
0081 process.SKIMStreamZMM = cms.OutputModule("PoolOutputModule",
0082     process.FEVTEventContent,
0083     SelectEvents = cms.untracked.PSet(
0084         SelectEvents = cms.vstring('zMMSkimPath')
0085     ),
0086     fileName = cms.untracked.string('ZMM.root'),
0087     dataset = cms.untracked.PSet(
0088         filterName = cms.untracked.string('ZMM'),
0089         dataTier = cms.untracked.string('RAW-RECO')
0090     )
0091 )
0092
0093 # Other statements
0094 process.GlobalTag.globaltag = 'GR10_P_V12::All'
0095
0096 # Valid vertex filter
0097 process.primaryVertexFilter = cms.EDFilter("VertexSelector",
0098     src = cms.InputTag("hiSelectedVertex"),
0099     cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"),
0100     filter = cms.bool(True),   # otherwise it won't filter the events
0101     )
0102
0103 # Refine photon cuts
0104 process.goodPhotons.cut = cms.string('et > 20 && hadronicOverEm < 0.1 && r9 > 0.8 && sigmaIetaIeta > 0.002')
0105
0106 # Looser photon cuts for ZEE
0107 process.goodPhotonsForZEE = process.goodPhotons.clone(
0108     cut=cms.string('et > 20 && hadronicOverEm < 0.2 && r9 > 0.5 && sigmaIetaIeta > 0.002')
0109     )
0110 process.goodCleanPhotonsForZEE = process.goodPhotonsForZEE.clone(src=cms.InputTag("cleanPhotons"))
0111 process.twoPhotonFilter.src = cms.InputTag("goodPhotonsForZEE")
0112 process.photonCombiner.decay = cms.string('goodCleanPhotonsForZEE goodCleanPhotonsForZEE')
0113 process.fullZEESkimSequence = cms.Sequence(process.hltPhotonHI
0114                                            * process.primaryVertexFilter
0115                                            * process.goodPhotonsForZEE
0116                                            * process.twoPhotonFilter
0117                                            * process.hiPhotonCleaningSequence
0118                                            * process.goodCleanPhotonsForZEE
0119                                            * process.photonCombiner * process.photonPairCounter
0120                                            * process.siPixelRecHits * process.siStripMatchedRecHits
0121                                            * process.hiPrimSeeds * process.hiElectronSequence)
0122
0123
0124 # Higher trigger thresholds
0125 process.hltJetHI.HLTPaths = ["HLT_HIJet50U_Core"]
0126 process.hltPhotonHI.HLTPaths = ["HLT_HIPhoton20_Core"]
0127 process.hltZMMHI.HLTPaths = ["HLT_HIL2DoubleMu3_Core"]
0128
0129
0130 # Dijet requirement
0131 process.leadingCaloJet = cms.EDFilter( "LargestEtCaloJetSelector",
0132     src = cms.InputTag( "icPu5CaloJetsL2L3" ),
0133     filter = cms.bool( False ),
0134     maxNumber = cms.uint32( 1 )
0135     )
0136
0137 process.goodLeadingJet = cms.EDFilter("CaloJetSelector",
0138     src = cms.InputTag("leadingCaloJet"),
0139     cut = cms.string("et > 130")
0140     )
0141
0142 process.goodSecondJet = cms.EDFilter("CaloJetSelector",
0143     src = cms.InputTag("icPu5CaloJetsL2L3"),
0144     cut = cms.string("et > 50")
0145     )
0146
0147 process.backToBackDijets = cms.EDProducer("CandViewShallowCloneCombiner",
0148     checkCharge = cms.bool(False),
0149     cut = cms.string('abs(deltaPhi(daughter(0).phi,daughter(1).phi)) > 2.5'),
0150     decay = cms.string("goodLeadingJet goodSecondJet")
0151     )
0152
0153 process.dijetFilter = cms.EDFilter("CandViewCountFilter",
0154     src = cms.InputTag("backToBackDijets"),
0155     minNumber = cms.uint32(1)
0156     )
0157
0158 process.backToBackSequence = cms.Sequence(process.leadingCaloJet*process.goodLeadingJet*
0159                                           process.goodSecondJet*process.backToBackDijets*
0160                                           process.dijetFilter)
0161
0162 process.particleFlowSequence = cms.Sequence(process.siPixelRecHits * process.siStripMatchedRecHits *
0163                                             process.heavyIonTracking * process.HiParticleFlowReco)
0164
0165
0166 # Path and EndPath definitions
0167 process.zEESkimPath = cms.Path(process.fullZEESkimSequence)
0168
0169 process.photonSkimPath = cms.Path(process.photonSkimSequence*process.primaryVertexFilter)
0170
0171 process.diJetSkimPath = cms.Path(process.diJetSkimSequence*process.primaryVertexFilter*
0172                                  process.backToBackSequence*process.particleFlowSequence)
0173
0174 process.zMMSkimPath = cms.Path(process.zMMSkimSequence*process.primaryVertexFilter)
0175
0176 process.SKIMStreamDiJetOutPath = cms.EndPath(process.SKIMStreamDiJet)
0177
0178 process.SKIMStreamPhotonOutPath = cms.EndPath(process.SKIMStreamPhoton)
0179
0180 process.SKIMStreamZEEOutPath = cms.EndPath(process.SKIMStreamZEE)
0181
0182 process.SKIMStreamZMMOutPath = cms.EndPath(process.SKIMStreamZMM)
0183
0184
0185 # Schedule definition
0186 process.schedule = cms.Schedule(process.photonSkimPath,process.zMMSkimPath,process.zEESkimPath,process.diJetSkimPath,process.SKIMStreamDiJetOutPath,process.SKIMStreamPhotonOutPath,process.SKIMStreamZEEOutPath,process.SKIMStreamZMMOutPath)