Don't like ads? PRO users don't see any ads ;-)
Public Pastes
Guest

Untitled

By: a guest on May 5th, 2012  |  syntax: None  |  size: 4.03 KB  |  hits: 13  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. import FWCore.ParameterSet.Config as cms
  2.  
  3. from UserCode.L1TriggerDPG.l1Ntuple_cfg import *
  4.  
  5. # global tag
  6.  
  7. process.GlobalTag.globaltag = cms.string('GR_R_42_V14::All')
  8.  
  9. # load standard RECO for MET cleaning
  10. process.load('Configuration.StandardSequences.Reconstruction_cff')
  11.  
  12. ##################################good collisions############################################
  13. # This filter select ~73% of events in MinBias PD
  14. process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskAlgoTrigConfig_cff')
  15. process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
  16. process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
  17.  
  18. process.L1T1coll=process.hltLevel1GTSeed.clone()
  19. process.L1T1coll.L1TechTriggerSeeding = cms.bool(True)
  20. process.L1T1coll.L1SeedsLogicalExpression = cms.string('0 AND (40 OR 41) AND NOT (36 OR 37 OR 38 OR 39) AND NOT ((42 AND NOT 43) OR (43 AND NOT 42))')
  21.  
  22. process.primaryVertexFilter = cms.EDFilter("VertexSelector",
  23.    src = cms.InputTag("offlinePrimaryVertices"),
  24.    cut = cms.string("!isFake && ndof > 4 && abs(z) <= 15 && position.Rho <= 2"), # tracksSize() > 3 for the older cut
  25.    filter = cms.bool(True),   # otherwise it won't filter the events, just produce an empty vertex collection.
  26. )
  27.  
  28.  
  29. process.noscraping = cms.EDFilter("FilterOutScraping",
  30.   applyfilter = cms.untracked.bool(True),
  31.   debugOn = cms.untracked.bool(False),
  32.   numtrack = cms.untracked.uint32(10),
  33.   thresh = cms.untracked.double(0.25)
  34. )
  35.  
  36. ### HLT Filter ###
  37. from HLTrigger.HLTfilters.hltHighLevel_cfi import *
  38. process.HLTFilter  = hltHighLevel.clone(TriggerResultsTag = "TriggerResults::HLT", HLTPaths = ["HLT_L1Tech_BSC_minBias_threshold1_v1"])
  39.  
  40. # HCAL noise filter
  41. process.load('CommonTools/RecoAlgos/HBHENoiseFilter_cfi')
  42. process.hbhefilter = cms.Path(process.HBHENoiseFilter)
  43.  
  44.  
  45. # L1 ntuples
  46. process.l1NtupleProducer.generatorSource      = cms.InputTag("none")
  47. process.l1NtupleProducer.simulationSource     = cms.InputTag("none")
  48. process.l1NtupleProducer.gmtSource            = cms.InputTag("gtDigis")
  49. process.l1NtupleProducer.gtEvmSource          = cms.InputTag("none")
  50. process.l1NtupleProducer.gtSource             = cms.InputTag("gtDigis")
  51. process.l1NtupleProducer.gctIsoEmSource       = cms.InputTag("gctDigis","isoEm")
  52. process.l1NtupleProducer.gctNonIsoEmSource    = cms.InputTag("gctDigis","nonIsoEm")
  53. process.l1NtupleProducer.gctCentralJetsSource = cms.InputTag("gctDigis","cenJets")
  54. process.l1NtupleProducer.gctTauJetsSource     = cms.InputTag("gctDigis","tauJets")
  55. process.l1NtupleProducer.gctForwardJetsSource = cms.InputTag("gctDigis","forJets")
  56. process.l1NtupleProducer.gctEnergySumsSource  = cms.InputTag("gctDigis")
  57. process.l1NtupleProducer.rctSource            = cms.InputTag("gctDigis")
  58. process.l1NtupleProducer.dttfSource           = cms.InputTag("none")
  59. process.l1NtupleProducer.csctfTrkSource       = cms.InputTag("none")
  60. process.l1NtupleProducer.csctfLCTSource       = cms.InputTag("none")
  61. process.l1NtupleProducer.csctfStatusSource    = cms.InputTag("none")
  62.  
  63. process.l1RecoTreeProducer.jetptThreshold = cms.double(5)
  64.  
  65. process.gctDigis.numberOfGctSamplesToUnpack = cms.uint32(5)
  66. process.l1extraParticles.centralBxOnly = cms.bool(False)
  67.  
  68. process.p.remove(process.dttfDigis)
  69. process.p.remove(process.csctfDigis)
  70. process.p.remove(process.l1extraParticles)
  71. process.p.remove(process.l1ExtraTreeProducer)
  72. process.p.remove(process.l1MuonRecoTreeProducer)
  73.  
  74.  
  75. process.ntuple = cms.Path(
  76. #    process.L1T1coll
  77.     process.HLTFilter
  78.     +process.primaryVertexFilter
  79.     +process.noscraping
  80.     +process.HBHENoiseFilter
  81.     +process.l1NtupleProducer
  82.     +process.l1ExtraTreeProducer
  83.     +process.l1RecoTreeProducer
  84. )
  85. process.MessageLogger = cms.Service(
  86.     "MessageLogger",
  87.     destinations = cms.untracked.vstring('cout'),
  88.     cout         = cms.untracked.PSet(threshold = cms.untracked.string('INFO'))
  89.   )
  90.  
  91.  
  92.  
  93. # jet corrections
  94. process.load('JetMETCorrections.Configuration.DefaultJEC_cff')
  95. process.ak5CaloL1Offset.useCondDB = False
  96. process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1000) )
  97.  
  98.  
  99. readFiles.extend( [
  100.         'file:////vols/cms03/bm409/JetMay10ReReco_AOD.root'
  101. ] )
create a new version of this paste RAW Paste Data