| 1 |
import FWCore.ParameterSet.Config as cms
|
| 2 |
|
| 3 |
process = cms.Process("MPIntuple")
|
| 4 |
|
| 5 |
process.load("FWCore.MessageService.MessageLogger_cfi")
|
| 6 |
process.load("FWCore.MessageLogger.MessageLogger_cfi")
|
| 7 |
|
| 8 |
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(1000000)
|
| 9 |
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1000))
|
| 10 |
|
| 11 |
process.load("Configuration.StandardSequences.Services_cff")
|
| 12 |
process.load('Configuration.StandardSequences.GeometryExtended_cff')
|
| 13 |
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
|
| 14 |
process.load("Configuration.StandardSequences.Reconstruction_cff")
|
| 15 |
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
|
| 16 |
|
| 17 |
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
|
| 18 |
process.load('HLTrigger.HLTfilters.hltLevel1GTSeed_cfi')
|
| 19 |
process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
|
| 20 |
process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND NOT (36 OR 37 OR 38 OR 39)')
|
| 21 |
|
| 22 |
process.noscraping = cms.EDFilter("FilterOutScraping",
|
| 23 |
applyfilter = cms.untracked.bool(True),
|
| 24 |
debugOn = cms.untracked.bool(False),
|
| 25 |
numtrack = cms.untracked.uint32(10),
|
| 26 |
thresh = cms.untracked.double(0.25)
|
| 27 |
)
|
| 28 |
|
| 29 |
|
| 30 |
### Conditions tags
|
| 31 |
process.GlobalTag.globaltag = 'GR_R_311_V3::All'
|
| 32 |
#process.GlobalTag.globaltag = 'FT_R38X_V14A::All'
|
| 33 |
|
| 34 |
### Input files
|
| 35 |
process.load("Data")
|
| 36 |
|
| 37 |
#process.source = cms.Source("PoolSource",
|
| 38 |
# fileNames = cms.untracked.vstring(
|
| 39 |
# '/store/data/Run2010A/JetMET/RECO/Sep17ReReco_v2/0026/FE9C387C-24C8-DF11-8CAA-003048679274.root'
|
| 40 |
#)
|
| 41 |
#)
|
| 42 |
|
| 43 |
# Select primary vertices
|
| 44 |
process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVerticesDA_cfi")
|
| 45 |
process.offlinePrimaryVerticesDAWithBS = process.offlinePrimaryVerticesDA.clone()
|
| 46 |
process.offlinePrimaryVerticesDAWithBS.useBeamConstraint = cms.bool(True)
|
| 47 |
process.offlinePrimaryVerticesDAWithBS.TkClusParameters.TkDAClusParameters.Tmin= cms.double(4.)
|
| 48 |
process.offlinePrimaryVerticesDAWithBS.TkClusParameters.TkDAClusParameters.vertexSize= cms.double(0.01)
|
| 49 |
|
| 50 |
### Jet correction services
|
| 51 |
process.load('JetMETCorrections.Configuration.DefaultJEC_cff')
|
| 52 |
|
| 53 |
### Extra jet collection for L1 corrections
|
| 54 |
process.load("RecoJets.Configuration.RecoPFJets_cff")
|
| 55 |
process.kt6PFJets.doRhoFastjet = True
|
| 56 |
process.kt6PFJets.Rho_EtaMax = cms.double(2.5)
|
| 57 |
process.kt6PFJets.rParam = cms.double(0.6)
|
| 58 |
|
| 59 |
process.ak5PFJetsL1 = process.ak5PFJets.clone()
|
| 60 |
process.ak5PFJetsL1.doAreaFastjet = True
|
| 61 |
process.ak5PFJetsL1.Rho_EtaMax = cms.double(2.5)
|
| 62 |
|
| 63 |
### To get b-tags from ak5PFJets
|
| 64 |
process.load('RecoJets.JetAssociationProducers.ak5JTA_cff')
|
| 65 |
process.ak5JetTracksAssociatorAtVertex.jets = cms.InputTag("ak5PFJetsL2L3")
|
| 66 |
process.ak5JetTracksAssociatorAtCaloFace.jets = cms.InputTag("ak5PFJetsL2L3")
|
| 67 |
process.ak5JetExtender.jets = cms.InputTag("ak5PFJetsL2L3")
|
| 68 |
|
| 69 |
### DPS ntuple producer
|
| 70 |
process.MPIntuple = cms.EDAnalyzer('MPIntuple',
|
| 71 |
|
| 72 |
RecoJetTag = cms.untracked.InputTag("ak5PFJets"),
|
| 73 |
RecoMETTag = cms.untracked.InputTag("pfMet"),
|
| 74 |
ElectronTag = cms.untracked.InputTag("gsfElectrons"),
|
| 75 |
MuonTag = cms.untracked.InputTag("muons"),
|
| 76 |
GenJetTag = cms.untracked.InputTag("ak5GenJets"),
|
| 77 |
PrimaryVtxTag = cms.untracked.InputTag("offlinePrimaryVertices"),
|
| 78 |
electronIDMap = cms.InputTag("eidTight"),
|
| 79 |
nJets = cms.untracked.int32(2),
|
| 80 |
doPFJets = cms.untracked.bool(True),
|
| 81 |
doGenJets = cms.untracked.bool(False),
|
| 82 |
triggerHLT = cms.untracked.bool(True),
|
| 83 |
hltName = cms.untracked.string("HLT"),
|
| 84 |
rootfilename = cms.untracked.string("nuTuple_Mu_Run2010A.root"),
|
| 85 |
triggers = cms.untracked.vstring("HLT_L1_BscMinBiasOR_BptxPlusORMinus", "HLT_L1Jet6U", "HLT_L1Jet10U", "HLT_Jet15U", "HLT_Jet30U", "HLT_Jet50U", "HLT_Jet70U", "HLT_Jet100U", "HLT_FwdJet20U", "HLT_Jet15U_HcalNoiseFiltered", "HLT_BTagIP_Jet50U", "HLT_DiJetAve15U", "HLT_DiJetAve30U", "HLT_DiJetAve50U", "HLT_DoubleJet15U_ForwardBackward", "HLT_QuadJet15U", "HLT_HT100U", "HLT_Mu9", "HLT_Mu11", "HLT_Mu13", "HLT_Mu15", "HLT_Ele10_LW_L1R", "HLT_Ele15_SW_L1R", "")
|
| 86 |
)
|
| 87 |
|
| 88 |
### Let it run
|
| 89 |
cmsSeq = cms.Sequence(
|
| 90 |
process.hltLevel1GTSeed
|
| 91 |
* process.noscraping
|
| 92 |
* process.offlinePrimaryVerticesDAWithBS
|
| 93 |
* process.ak5PFJetsL2L3
|
| 94 |
* process.ak5JetTracksAssociatorAtVertex
|
| 95 |
* process.btagging
|
| 96 |
)
|
| 97 |
|
| 98 |
process.p = cms.Path(cmsSeq * process.MPIntuple)
|
| 99 |
#process.p = cms.Path(process.hltLevel1GTSeed * process.noscraping * process.ak5PFJetsL2L3 * process.MPIntuple)
|
| 100 |
|
