| 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(500000)
|
| 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 |
### Conditions tags
|
| 18 |
process.GlobalTag.globaltag = 'START311_V2::All'
|
| 19 |
#process.GlobalTag.globaltag = 'START38_V14::All'
|
| 20 |
|
| 21 |
### Input files
|
| 22 |
|
| 23 |
#process.load('Signal')
|
| 24 |
|
| 25 |
process.source = cms.Source("PoolSource",
|
| 26 |
fileNames = cms.untracked.vstring(
|
| 27 |
|
| 28 |
# 'file:/uscms/home/naodell/nobackup/CA5A27DF-0AE8-DF11-9C93-E0CB4E553686.root'
|
| 29 |
|
| 30 |
'/store/mc/Fall10/QCD_Pt_120to170_TuneZ2_7TeV_pythia6/GEN-SIM-RECO/E7TeV_ProbDist_2010Data_BX156_START38_V12-v1/0000/FE4616C4-EDE3-DF11-A938-001A64789DA8.root'
|
| 31 |
# '/store/mc/Spring10/QCD_Pt-15_7TeV-pythia8/GEN-SIM-RECO/START3X_V26B-v2/0033/88571403-6E73-DF11-B10F-0015174ED57A.root',
|
| 32 |
# '/store/mc/Spring10/QCD_Pt-15_7TeV-pythia8/GEN-SIM-RECO/START3X_V26B-v2/0029/F45EEABC-EE71-DF11-8969-A4BADB1E65B1.root'
|
| 33 |
|
| 34 |
# '/store/mc/Fall10/QCD_Pt_30to50_Tune1_7TeV_pythia8/GEN-SIM-RECO/START38_V12-v1/0006/562C8C93-FCCA-DF11-96E6-E0CB4EA0A8FE.root'
|
| 35 |
# '/store/mc/Fall10/QCD_Pt_80to120_TuneZ2_7TeV_pythia6/GEN-SIM-RECO/START38_V12-v1/0000/FE6C9B24-41CB-DF11-9325-00E08178C0CD.root'
|
| 36 |
# '/store/mc/Fall10/WWtoAnything_TuneZ2_7TeV-pythia6-tauola/GEN-SIM-RECO/START38_V12-v1/0000/FEC1FFE9-06CD-DF11-A128-00E0817918C5.root'
|
| 37 |
)
|
| 38 |
)
|
| 39 |
|
| 40 |
# Select primary vertices
|
| 41 |
process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVerticesDA_cfi")
|
| 42 |
process.offlinePrimaryVerticesDAWithBS = process.offlinePrimaryVerticesDA.clone()
|
| 43 |
process.offlinePrimaryVerticesDAWithBS.useBeamConstraint = cms.bool(True)
|
| 44 |
process.offlinePrimaryVerticesDAWithBS.TkClusParameters.TkDAClusParameters.Tmin= cms.double(4.)
|
| 45 |
process.offlinePrimaryVerticesDAWithBS.TkClusParameters.TkDAClusParameters.vertexSize= cms.double(0.01)
|
| 46 |
|
| 47 |
### Jet correction services
|
| 48 |
process.load('JetMETCorrections.Configuration.DefaultJEC_cff')
|
| 49 |
|
| 50 |
### Extra jet collection for L1 corrections
|
| 51 |
process.load("RecoJets.Configuration.RecoPFJets_cff")
|
| 52 |
process.kt6PFJets.doRhoFastjet = True
|
| 53 |
process.kt6PFJets.Rho_EtaMax = cms.double(2.5)
|
| 54 |
process.kt6PFJets.rParam = cms.double(0.6)
|
| 55 |
|
| 56 |
process.ak5PFJetsL1 = process.ak5PFJets.clone()
|
| 57 |
process.ak5PFJetsL1.doAreaFastjet = True
|
| 58 |
process.ak5PFJetsL1.Rho_EtaMax = cms.double(2.5)
|
| 59 |
|
| 60 |
### To get b-tags from ak5PFJets
|
| 61 |
process.load('RecoJets.JetAssociationProducers.ak5JTA_cff')
|
| 62 |
process.ak5JetTracksAssociatorAtVertex.jets = cms.InputTag("ak5PFJetsL2L3")
|
| 63 |
process.ak5JetTracksAssociatorAtCaloFace.jets = cms.InputTag("ak5PFJetsL2L3")
|
| 64 |
process.ak5JetExtender.jets = cms.InputTag("ak5PFJetsL2L3")
|
| 65 |
|
| 66 |
### To include ak5 GenJets
|
| 67 |
process.load('RecoJets.Configuration.GenJetParticles_cff')
|
| 68 |
process.load('RecoJets.Configuration.RecoGenJets_cff')
|
| 69 |
|
| 70 |
### DPS tree producer
|
| 71 |
|
| 72 |
process.MPIntuple = cms.EDAnalyzer('MPIntuple',
|
| 73 |
|
| 74 |
RecoJetTag = cms.untracked.InputTag("ak5PFJets"),
|
| 75 |
RecoMETTag = cms.untracked.InputTag("pfMet"),
|
| 76 |
ElectronTag = cms.untracked.InputTag("gsfElectrons"),
|
| 77 |
MuonTag = cms.untracked.InputTag("muons"),
|
| 78 |
GenJetTag = cms.untracked.InputTag("ak5GenJets"),
|
| 79 |
PrimaryVtxTag = cms.untracked.InputTag("offlinePrimaryVertices"),
|
| 80 |
electronIDMap = cms.InputTag("eidTight"),
|
| 81 |
nJets = cms.untracked.int32(4),
|
| 82 |
doPFJets = cms.untracked.bool(True),
|
| 83 |
doGenJets = cms.untracked.bool(True),
|
| 84 |
triggerHLT = cms.untracked.bool(True),
|
| 85 |
hltName = cms.untracked.string("REDIGI38XPU"),
|
| 86 |
rootfilename = cms.untracked.string("nuTuple_QCDpy8_15to30_PU.root"),
|
| 87 |
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")
|
| 88 |
)
|
| 89 |
|
| 90 |
### Let it run
|
| 91 |
cmsSeq = cms.Sequence(
|
| 92 |
process.offlinePrimaryVerticesDAWithBS
|
| 93 |
* process.ak5PFJetsL2L3
|
| 94 |
* process.ak5JetTracksAssociatorAtVertex
|
| 95 |
* process.btagging
|
| 96 |
)
|
| 97 |
|
| 98 |
process.p = cms.Path(cmsSeq * process.MPIntuple)
|
| 99 |
|
| 100 |
#process.p = cms.Path(process.genParticlesForJets * process.ak5GenJets * process.ak5PFJetsL2L3 * process.MPIntuple)
|
