| 1 |
import FWCore.ParameterSet.Config as cms
|
| 2 |
process = cms.Process("ReRec")
|
| 3 |
|
| 4 |
process.load("FWCore.MessageLogger.MessageLogger_cfi")
|
| 5 |
process.load("CondCore.DBCommon.CondDBSetup_cfi")
|
| 6 |
|
| 7 |
process.maxEvents = cms.untracked.PSet(
|
| 8 |
input = cms.untracked.int32(-1)
|
| 9 |
)
|
| 10 |
|
| 11 |
process.source = cms.Source("PoolSource",
|
| 12 |
fileNames = cms.untracked.vstring(
|
| 13 |
# pointing skim
|
| 14 |
#'/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/BSCNOBEAMHALO-Jan29Skim-v2/0019/0ACCA1F9-E00D-DF11-860C-0018F3D095F6.root'
|
| 15 |
#'/store/data/Commissioning10/MinimumBias/RAW-RECO/Apr1Skim_GOODCOLL-v1/0140/F4EA88D2-C83E-DF11-AF55-00261894392C.root',
|
| 16 |
#'/store/data/Commissioning10/MinimumBias/RAW-RECO/Apr1Skim_GOODCOLL-v1/0140/F4EA88D2-C83E-DF11-AF55-00261894392C.root',
|
| 17 |
#'/store/data/Commissioning10/MinimumBias/RAW-RECO/Apr1Skim_GOODCOLL-v1/0140/E27B88D1-8040-DF11-B3FC-00261894391B.root'
|
| 18 |
#'file:/data1/speer/Tk-DPG/data/Prompt_MB_Run135523_FCE0A77A-7360-DF11-8D7F-001D09F29321.root'
|
| 19 |
#'/store/relval/CMSSW_3_6_0_pre6/RelValMinBias/ALCARECO/START36_V4_SiStripCalMinBias-v1/0011/3AF6DD26-4F45-DF11-9FCF-0018F3D096BE.root'
|
| 20 |
'/store/data/Commissioning10/MinimumBias/RECO/GOODCOLL-Jun14thSkim_v1/0148/FED301EB-6284-DF11-9D47-002618943906.root',
|
| 21 |
|
| 22 |
)
|
| 23 |
)
|
| 24 |
|
| 25 |
|
| 26 |
process.options = cms.untracked.PSet(
|
| 27 |
wantSummary = cms.untracked.bool(True)
|
| 28 |
)
|
| 29 |
|
| 30 |
|
| 31 |
# output module#
|
| 32 |
#process.load("Configuration.EventContent.EventContentCosmics_cff")
|
| 33 |
# Carlotta: event content for collisions, let's try
|
| 34 |
process.load("Configuration.EventContent.EventContent_cff")
|
| 35 |
process.load("CalibTracker.SiStripESProducers.SiStripQualityESProducer_cfi")
|
| 36 |
#process.load("RecoLocalTracker.Configuration.RecoLocalTracker_Cosmics_cff")
|
| 37 |
process.load("RecoLocalTracker.Configuration.RecoLocalTracker_cff")
|
| 38 |
process.load("RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi")
|
| 39 |
|
| 40 |
# Conditions (Global Tag is used here):
|
| 41 |
process.load("Configuration.StandardSequences.Geometry_cff")
|
| 42 |
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
|
| 43 |
#process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")
|
| 44 |
process.GlobalTag.globaltag = 'GR_R_36X_V12A::All'
|
| 45 |
|
| 46 |
process.load("Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff")
|
| 47 |
|
| 48 |
|
| 49 |
#Geometry
|
| 50 |
process.load("Configuration.StandardSequences.Geometry_cff")
|
| 51 |
process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")
|
| 52 |
# Real data raw to digi
|
| 53 |
process.load("Configuration.StandardSequences.RawToDigi_Data_cff")
|
| 54 |
|
| 55 |
# reconstruction sequence
|
| 56 |
#process.load("Configuration.StandardSequences.ReconstructionCosmics_cff")
|
| 57 |
process.load("Configuration.StandardSequences.Reconstruction_cff")
|
| 58 |
|
| 59 |
# offline DQM
|
| 60 |
#process.load("DQMOffline.Configuration.DQMOfflineCosmics_cff")
|
| 61 |
process.load("DQMOffline.Configuration.DQMOffline_cff")
|
| 62 |
process.load("DQMServices.Components.MEtoEDMConverter_cff")
|
| 63 |
process.load("RecoTracker.TrackProducer.TrackRefitters_cff")
|
| 64 |
process.TrackRefitter.src = cms.InputTag("myTrackSelector")
|
| 65 |
#process.TrackRefitter.src = "myTrackSelector"
|
| 66 |
|
| 67 |
process.TFileService = cms.Service("TFileService",
|
| 68 |
fileName = cms.string('HitRes.root')
|
| 69 |
)
|
| 70 |
|
| 71 |
#process.source.lumisToProcess = cms.untracked.VLuminosityBlockRange('123592:2-123592:12','123596:1-123596:9999','123615:66-123615:9999','123732:57-123732:109','123815:8-123815:9999','123818:2-123818:42','123906:18-123906:28','123908:1-123908:12','124006:1-124006:9999','124008:1-124008:9999','124009:1-124009:68','124020:12-124020:9999','124022:69-124022:9999','124023:39-124023:9999','124024:1-124024:9999','124025:1-124025:13','124027:24-124027:9999','124030:1-124030:9999','1240120:1-124120:9999','124230:26-124230:9999','124275:1-124275:30')
|
| 72 |
|
| 73 |
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
|
| 74 |
|
| 75 |
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
|
| 76 |
#process.bitsel = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('(42 AND 43) OR (NOT 42 AND NOT 43)'))
|
| 77 |
|
| 78 |
process.bitsel = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('0 AND ((42 AND 43) OR (NOT 42 AND NOT 43))'))
|
| 79 |
|
| 80 |
|
| 81 |
process.myTrackSelector = cms.EDFilter("TrackSelector",
|
| 82 |
src = cms.InputTag("generalTracks"),
|
| 83 |
cut = cms.string('pt > 3'),
|
| 84 |
filter = cms.bool(True)
|
| 85 |
)
|
| 86 |
|
| 87 |
process.analysis = cms.EDFilter("HitRes",
|
| 88 |
usePXB = cms.bool(True),
|
| 89 |
usePXF = cms.bool(False),
|
| 90 |
useTIB = cms.bool(True),
|
| 91 |
useTOB = cms.bool(True),
|
| 92 |
useTID = cms.bool(False),
|
| 93 |
useTEC = cms.bool(False),
|
| 94 |
ROUList = cms.vstring('TrackerHitsTIBLowTof',
|
| 95 |
'TrackerHitsTIBHighTof',
|
| 96 |
'TrackerHitsTOBLowTof',
|
| 97 |
'TrackerHitsTOBHighTof'),
|
| 98 |
trajectories = cms.InputTag("TrackRefitter"),
|
| 99 |
associatePixel = cms.bool(False),
|
| 100 |
associateStrip = cms.bool(False),
|
| 101 |
associateRecoTracks = cms.bool(False),
|
| 102 |
tracks = cms.InputTag("TrackRefitter"),
|
| 103 |
barrelOnly = cms.bool(False)
|
| 104 |
)
|
| 105 |
|
| 106 |
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
|
| 107 |
process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
|
| 108 |
process.L1T1=process.hltLevel1GTSeed.clone()
|
| 109 |
process.L1T1.L1TechTriggerSeeding = cms.bool(True)
|
| 110 |
process.L1T1.L1SeedsLogicalExpression = cms.string('0 AND (40 OR 41) AND NOT (36 OR 37 OR 38 OR 39)')
|
| 111 |
process.bscnobeamhalo = cms.Path(process.L1T1)
|
| 112 |
|
| 113 |
process.hltHighLevel = cms.EDFilter("HLTHighLevel",
|
| 114 |
TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
|
| 115 |
HLTPaths = cms.vstring('HLT_L1_BscMinBiasOR_BptxPlusORMinus'), # provide list of HLT paths (or patterns) you want
|
| 116 |
eventSetupPathsKey = cms.string(''), # not empty => use read paths from AlCaRecoTriggerBitsRcd via this key
|
| 117 |
andOr = cms.bool(True), # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
|
| 118 |
throw = cms.bool(True) # throw exception on unknown path names
|
| 119 |
)
|
| 120 |
|
| 121 |
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
|
| 122 |
from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
|
| 123 |
process.bit40 = hltLevel1GTSeed.clone(L1TechTriggerSeeding = cms.bool(True), L1SeedsLogicalExpression = cms.string('0 AND (40 OR 41) AND NOT (36 OR 37 OR 38 OR 39)'))
|
| 124 |
|
| 125 |
|
| 126 |
|
| 127 |
process.oneGoodVertexFilter = cms.EDFilter("VertexSelector",
|
| 128 |
src = cms.InputTag("offlinePrimaryVertices"),
|
| 129 |
cut = cms.string("!isFake && ndof > 4 && abs(z) <= 15 && position.Rho <= 2"), # tracksSize() > 3 for the older cut
|
| 130 |
filter = cms.bool(True), # otherwise it won't filter the events, just produce an empty vertex collection.
|
| 131 |
)
|
| 132 |
|
| 133 |
|
| 134 |
process.noScraping= cms.EDFilter("FilterOutScraping",
|
| 135 |
applyfilter = cms.untracked.bool(True),
|
| 136 |
debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
|
| 137 |
numtrack = cms.untracked.uint32(10),
|
| 138 |
thresh = cms.untracked.double(0.25)
|
| 139 |
)
|
| 140 |
|
| 141 |
#For 7TeV collisions.
|
| 142 |
process.skimming = cms.EDFilter("PhysDecl",
|
| 143 |
HLTriggerResults = cms.InputTag("TriggerResults","","HLT"),
|
| 144 |
applyfilter = cms.untracked.bool(True)
|
| 145 |
)
|
| 146 |
|
| 147 |
process.load('HLTrigger.special.hltPhysicsDeclared_cfi')
|
| 148 |
process.hltPhysicsDeclared.L1GtReadoutRecordTag = 'gtDigis'
|
| 149 |
|
| 150 |
### Apply corrections to local reco - needed only if you run on TK-DECO data (i.e., MinBias) !!!! For cosmics in PEAK mode, comment them out
|
| 151 |
process.load("RecoLocalTracker.SiStripRecHitConverter.OutOfTime_cff")
|
| 152 |
process.OutOfTime.TOBlateBP=0.071
|
| 153 |
process.OutOfTime.TIBlateBP=0.036
|
| 154 |
|
| 155 |
from CondCore.DBCommon.CondDBSetup_cfi import *
|
| 156 |
|
| 157 |
#process.stripLorentzAngle = cms.ESSource("PoolDBESSource",CondDBSetup,
|
| 158 |
# connect = cms.string('sqlite_file:./SiStripLorentzAngle_Deco.db'),
|
| 159 |
# toGet = cms.VPSet(cms.PSet(record = cms.string('SiStripLorentzAngleRcd'),
|
| 160 |
# tag = cms.string('SiStripLorentzAngle_Deco') ))
|
| 161 |
# )
|
| 162 |
#process.es_prefer_stripLorentzAngle = cms.ESPrefer("PoolDBESSource", "stripLorentzAngle")
|
| 163 |
|
| 164 |
|
| 165 |
# Path and EndPath definitions
|
| 166 |
process.p = cms.Path(process.hltPhysicsDeclared*process.L1T1*process.oneGoodVertexFilter*process.noScraping*process.myTrackSelector*process.TrackRefitter*process.analysis)
|
| 167 |
|
| 168 |
|
| 169 |
#process.p = cms.Path(process.noScraping*process.myTrackSelector*process.TrackRefitter*process.analysis)
|
| 170 |
|
| 171 |
|
| 172 |
#FOR MC>>>>>>>>>
|
| 173 |
#process.p = cms.Path(process.L1T1*process.oneGoodVertexFilter*process.TrackRefitter*process.noScraping*process.analysis)
|
| 174 |
|
| 175 |
|
| 176 |
|
| 177 |
#process.p = cms.Path(process.physDecl*process.L1T1*process.oneGoodVertexFilter*process.TrackRefitter*process.noScraping)
|
| 178 |
|
| 179 |
#process.p = cms.Path(process.bitsel*process.TrackRefitter*process.analysis*process.bit40*process.bptxAnd*process.oneGoodVertexFilter)
|
| 180 |
#process.p = cms.Path(process.TrackRefitterP5)
|
| 181 |
#process.p = cms.Path(process.TrackRefitter)
|
| 182 |
#process.endjob_step = cms.Path(process.endOfProcess)
|
| 183 |
|
| 184 |
#process.outpath = cms.EndPath(process.output)
|
| 185 |
|
| 186 |
# Schedule definition
|
| 187 |
|
| 188 |
#process.schedule = cms.Schedule(process.endjob_step)
|
| 189 |
#process.schedule = cms.Schedule(process.p,process.endjob_step)
|
| 190 |
process.schedule = cms.Schedule(process.p)
|
| 191 |
#process.schedule = cms.Schedule(process.p,process.outpath)
|
| 192 |
|
| 193 |
|
