| 7 |
|
# import of standard configurations |
| 8 |
|
process.load('Configuration/StandardSequences/Services_cff') |
| 9 |
|
process.load('FWCore/MessageService/MessageLogger_cfi') |
| 10 |
< |
process.load('Configuration/StandardSequences/GeometryIdeal_cff') |
| 11 |
< |
process.load('Configuration/StandardSequences/MagneticField_38T_cff') |
| 10 |
> |
#process.load('Configuration/StandardSequences/GeometryIdeal_cff') |
| 11 |
> |
#process.load('Configuration/StandardSequences/MagneticField_38T_cff') |
| 12 |
|
process.load('Configuration/StandardSequences/SimL1Emulator_cff') |
| 13 |
+ |
#process.load('Configuration/StandardSequences/GeometryDB_cff') |
| 14 |
+ |
#process.load('Configuration/StandardSequences/MagneticField_38T_cff') |
| 15 |
|
process.load("Configuration.StandardSequences.RawToDigi_Data_cff") |
| 16 |
+ |
process.load('Configuration.StandardSequences.L1Reco_cff') |
| 17 |
+ |
process.load('Configuration.StandardSequences.Reconstruction_cff') |
| 18 |
|
process.load('Configuration/StandardSequences/EndOfProcess_cff') |
| 19 |
|
process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff') |
| 20 |
|
process.load('Configuration/EventContent/EventContent_cff') |
| 21 |
< |
process.load('Configuration.StandardSequences.ReconstructionCosmics_cff') |
| 21 |
> |
process.load('Configuration/StandardSequences/Geometry_cff') |
| 22 |
> |
process.load('Configuration/StandardSequences/MagneticField_AutoFromDBCurrent_cff') |
| 23 |
> |
process.load("JetMETCorrections.Configuration.DefaultJEC_cff") |
| 24 |
|
|
| 25 |
|
# global tag |
| 26 |
< |
process.GlobalTag.globaltag = 'GR10_P_V2::All' |
| 26 |
> |
process.GlobalTag.globaltag = 'GR_P_V14::All' |
| 27 |
|
|
| 28 |
|
# output file |
| 29 |
|
process.TFileService = cms.Service("TFileService", |
| 30 |
|
fileName = cms.string('L1Tree.root') |
| 31 |
|
) |
| 32 |
|
|
| 27 |
– |
#is it MC or DATA |
| 28 |
– |
#WARNING: FOR MC turn isMC = True, otherwise the v4 of HF cleaning will be used, which includes timing cut. Timing is not modeled well in MC |
| 29 |
– |
isMC = False |
| 30 |
– |
useHBHEcleaning = True |
| 31 |
– |
useHBHEfilter = True |
| 32 |
– |
|
| 33 |
– |
HFPMTcleaningversion = 4 # version 1 = (loose), version 2 = (medium), version 3 = (tight) |
| 34 |
– |
# VERSION 4 is the currently recommended version, as of 28 May 2010. |
| 35 |
– |
|
| 36 |
– |
if useHBHEfilter == True: |
| 37 |
– |
process.load('CommonTools/RecoAlgos/HBHENoiseFilter_cfi') |
| 38 |
– |
process.hbhefilter = cms.Path(process.HBHENoiseFilter) |
| 39 |
– |
|
| 40 |
– |
# New SeverityLevelComputer that forces RecHits with UserDefinedBit0 set to be excluded from new rechit collection |
| 41 |
– |
import JetMETAnalysis.HcalReflagging.RemoveAddSevLevel as RemoveAddSevLevel |
| 42 |
– |
process.hcalRecAlgos=RemoveAddSevLevel.RemoveFlag(process.hcalRecAlgos,"HFLongShort") |
| 43 |
– |
|
| 44 |
– |
# UserDefinedBit0 is used by both the HF and HBHE reflaggers |
| 45 |
– |
process.hcalRecAlgos=RemoveAddSevLevel.AddFlag(process.hcalRecAlgos,"UserDefinedBit0",10) |
| 46 |
– |
|
| 47 |
– |
# HF RecHit reflagger |
| 48 |
– |
process.load("JetMETAnalysis/HcalReflagging/HFrechitreflaggerJETMET_cff") |
| 49 |
– |
if HFPMTcleaningversion==1: |
| 50 |
– |
process.hfrecoReflagged = process.HFrechitreflaggerJETMETv1.clone() |
| 51 |
– |
elif HFPMTcleaningversion==2: |
| 52 |
– |
process.hfrecoReflagged = process.HFrechitreflaggerJETMETv2.clone() |
| 53 |
– |
elif HFPMTcleaningversion==3: |
| 54 |
– |
process.hfrecoReflagged = process.HFrechitreflaggerJETMETv3.clone() |
| 55 |
– |
elif HFPMTcleaningversion==4: |
| 56 |
– |
if (isMC==False): |
| 57 |
– |
process.hfrecoReflagged = process.HFrechitreflaggerJETMETv4.clone() |
| 58 |
– |
else: |
| 59 |
– |
process.hfrecoReflagged = process.HFrechitreflaggerJETMETv2.clone() |
| 60 |
– |
elif HFPMTcleaningversion==5: |
| 61 |
– |
if (isMC==False): |
| 62 |
– |
process.hfrecoReflagged = process.HFrechitreflaggerJETMETv5.clone() |
| 63 |
– |
else: |
| 64 |
– |
process.hfrecoReflagged = process.HFrechitreflaggerJETMETv3.clone() |
| 65 |
– |
|
| 66 |
– |
|
| 67 |
– |
# HBHE RecHit reflagger |
| 68 |
– |
process.load("JetMETAnalysis/HcalReflagging/hbherechitreflaggerJETMET_cfi") |
| 69 |
– |
process.hbherecoReflagged = process.hbherechitreflaggerJETMET.clone() |
| 70 |
– |
process.hbherecoReflagged.debug=0 |
| 71 |
– |
|
| 72 |
– |
# Use the reflagged HF RecHits to make the CaloTowers |
| 73 |
– |
process.towerMaker.hfInput = "hfrecoReflagged" |
| 74 |
– |
process.towerMakerWithHO.hfInput = "hfrecoReflagged" |
| 75 |
– |
|
| 76 |
– |
# Path and EndPath definitions |
| 77 |
– |
|
| 78 |
– |
if (useHBHEcleaning==False): |
| 79 |
– |
process.reflagging_step = cms.Path(process.hfrecoReflagged) |
| 80 |
– |
else: |
| 81 |
– |
process.reflagging_step = cms.Path(process.hfrecoReflagged+process.hbherecoReflagged) |
| 82 |
– |
# Need to specify that new HBHE collection should be fed to calotower maker |
| 83 |
– |
process.towerMaker.hbheInput = "hbherecoReflagged" |
| 84 |
– |
process.towerMakerWithHO.hbheInput = "hbherecoReflagged" |
| 85 |
– |
|
| 86 |
– |
# Instead of rejecting the event, add a flag indicating the HBHE noise |
| 87 |
– |
process.load('CommonTools/RecoAlgos/HBHENoiseFilterResultProducer_cfi') |
| 88 |
– |
process.hbheflag = cms.Path(process.HBHENoiseFilterResultProducer) |
| 89 |
– |
|
| 90 |
– |
|
| 91 |
– |
process.rereco_step = cms.Path(process.filterSequence*process.caloTowersRec*(process.recoJets*process.recoJetIds+process.recoTrackJets)*process.recoJetAssociations*process.btagging*process.metreco) # re-reco jets and MET |
| 92 |
– |
|
| 93 |
– |
if useHBHEfilter == True: |
| 94 |
– |
process.schedule = cms.Schedule(process.hbhefilter, process.reflagging_step,process.rereco_step, USERCODE) |
| 95 |
– |
else: |
| 96 |
– |
process.schedule = cms.Schedule(process.reflagging_step, process.rereco_step, process.hbheflag, USERCODE) |
| 97 |
– |
|
| 98 |
– |
|
| 99 |
– |
|
| 33 |
|
# analysis |
| 34 |
|
process.load("L1Trigger.Configuration.L1Extra_cff") |
| 35 |
|
process.load("UserCode.L1TriggerDPG.l1NtupleProducer_cfi") |
| 36 |
|
process.load("UserCode.L1TriggerDPG.l1RecoTreeProducer_cfi") |
| 37 |
|
process.load("UserCode.L1TriggerDPG.l1ExtraTreeProducer_cfi") |
| 38 |
+ |
process.load("UserCode.L1TriggerDPG.l1MenuTreeProducer_cfi") |
| 39 |
|
process.load("UserCode.L1TriggerDPG.l1MuonRecoTreeProducer_cfi") |
| 40 |
+ |
process.load("EventFilter.L1GlobalTriggerRawToDigi.l1GtTriggerMenuLite_cfi") |
| 41 |
|
|
| 42 |
|
process.gctDigis.numberOfGctSamplesToUnpack = cms.uint32(5) |
| 43 |
|
process.l1extraParticles.centralBxOnly = cms.bool(False) |
| 44 |
|
|
| 45 |
|
process.p = cms.Path( |
| 46 |
< |
process.HBHENoiseFilter |
| 112 |
< |
+process.gtDigis |
| 46 |
> |
process.gtDigis |
| 47 |
|
+process.gtEvmDigis |
| 48 |
|
+process.gctDigis |
| 49 |
|
+process.dttfDigis |
| 51 |
|
+process.l1NtupleProducer |
| 52 |
|
+process.l1extraParticles |
| 53 |
|
+process.l1ExtraTreeProducer |
| 54 |
+ |
+process.l1GtTriggerMenuLite |
| 55 |
+ |
+process.l1MenuTreeProducer |
| 56 |
+ |
# +process.ak5CaloJetsCor |
| 57 |
|
+process.l1RecoTreeProducer |
| 58 |
|
+process.l1MuonRecoTreeProducer |
| 59 |
|
) |
| 67 |
|
secondaryFileNames = secFiles |
| 68 |
|
) |
| 69 |
|
|
| 70 |
+ |
|
| 71 |
+ |
|
| 72 |
|
readFiles.extend( [ |
| 73 |
< |
'/store/express/Commissioning10/ExpressPhysics/FEVT/v9/000/133/874/FEFC3201-644F-DF11-AED5-000423D98800.root' |
| 73 |
> |
# '/store/express/Commissioning10/ExpressPhysics/FEVT/v9/000/133/874/FEFC3201-644F-DF11-AED5-000423D98800.root' |
| 74 |
|
# '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/318/F23EF63B-7AD8-DE11-A6AC-0019B9F72F97.root', |
| 75 |
|
# '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/318/E0E8C7BF-7DD8-DE11-93F4-001617DC1F70.root', |
| 76 |
|
# '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/318/8E6024D2-7FD8-DE11-B2FD-001D09F295A1.root', |
