| 8 |
|
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(500000) |
| 9 |
|
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(1000)) |
| 10 |
|
|
| 11 |
– |
|
| 12 |
– |
### For data |
| 13 |
– |
|
| 11 |
|
process.load("Configuration.StandardSequences.Services_cff") |
| 12 |
|
process.load('Configuration.StandardSequences.GeometryExtended_cff') |
| 13 |
|
process.load('Configuration.StandardSequences.MagneticField_38T_cff') |
| 15 |
|
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") |
| 16 |
|
|
| 17 |
|
### Conditions tags |
| 18 |
< |
|
| 19 |
< |
#process.GlobalTag.globaltag = 'START36_V10::ALL' |
| 23 |
< |
process.GlobalTag.globaltag = 'START38_V14::All' |
| 24 |
< |
#process.GlobalTag.globaltag = 'MC_38Y_V13::ALL' |
| 18 |
> |
process.GlobalTag.globaltag = 'START311_V2::All' |
| 19 |
> |
#process.GlobalTag.globaltag = 'START38_V14::All' |
| 20 |
|
|
| 21 |
|
### Input files |
| 22 |
|
|
| 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' |
| 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' |
| 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 |
< |
### Jet correction services |
| 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 |
< |
from JetMETCorrections.Configuration.DefaultJEC_cff import * |
| 47 |
> |
### Jet correction services |
| 48 |
|
process.load('JetMETCorrections.Configuration.DefaultJEC_cff') |
| 49 |
|
|
| 50 |
< |
### To include ak5 GenJets |
| 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 |
|
|
| 55 |
– |
### To get b-tags from ak5PFJets |
| 56 |
– |
from RecoJets.JetAssociationProducers.ak5JTA_cff import * |
| 57 |
– |
ak5JetTracksAssociatorAtVertex.jets = cms.InputTag("ak5PFJetsL2L3") |
| 58 |
– |
ak5JetTracksAssociatorAtCaloFace.jets = cms.InputTag("ak5PFJetsL2L3") |
| 59 |
– |
ak5JetExtender.jets = cms.InputTag("ak5PFJetsL2L3") |
| 60 |
– |
|
| 70 |
|
### DPS tree producer |
| 71 |
|
|
| 72 |
|
process.MPIntuple = cms.EDAnalyzer('MPIntuple', |
| 82 |
|
doPFJets = cms.untracked.bool(True), |
| 83 |
|
doGenJets = cms.untracked.bool(True), |
| 84 |
|
triggerHLT = cms.untracked.bool(True), |
| 85 |
< |
hltName = cms.untracked.string("REDIGI38X"), |
| 86 |
< |
rootfilename = cms.untracked.string("nuTuple_WW.root"), |
| 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 |
|
|
| 83 |
– |
process.p = cms.Path(process.genParticlesForJets * process.ak5GenJets * process.ak5PFJetsL2L3 * process.ak5JetTracksAssociatorAtVertex * process.btagging * process.MPIntuple) |
| 100 |
|
#process.p = cms.Path(process.genParticlesForJets * process.ak5GenJets * process.ak5PFJetsL2L3 * process.MPIntuple) |
