/ElectroWeakAnalysis/WMuNu/test/WMuNuAODValidator.py

https://github.com/aivanov-cern/cmssw · Python · 84 lines · 47 code · 12 blank · 25 comment · 0 complexity · 668947cfce4ce660950107abfb485000 MD5 · raw file 

    

  1. import FWCore.ParameterSet.Config as cms
    2. 

  2. 

  3. # Process, how many events, inout files, ...
    4. 

  4. process = cms.Process("wmnsel")
    5. 

  5. process.maxEvents = cms.untracked.PSet(
    6. 

  6.       input = cms.untracked.int32(200)
    7. 

  7.       #input = cms.untracked.int32(100)
    8. 

  8. )
    9. 

  9. process.source = cms.Source("PoolSource",
    10. 

  10.       debugVerbosity = cms.untracked.uint32(0),
    11. 

  11.       debugFlag = cms.untracked.bool(False),
    12. 

  12.       #fileNames = cms.untracked.vstring("file:/data4/RelValWM_CMSSW_3_1_0-STARTUP31X_V1-v1_GEN-SIM-RECO/40BFAA1A-5466-DE11-B792-001D09F29533.root")
    13. 

  13.       #fileNames = cms.untracked.vstring("file:/data4/Wmunu-Summer09-MC_31X_V2_preproduction_311-v1/0011/F4C91F77-766D-DE11-981F-00163E1124E7.root")
    14. 

  14.        fileNames = cms.untracked.vstring("file:/data4/Wmunu_Summer09-MC_31X_V3_AODSIM-v1/0009/F82D4260-507F-DE11-B5D6-00093D128828.root")
    15. 

  15. 

  16. )
    17. 

  17. 

  18. # Debug/info printouts
    19. 

  19. process.MessageLogger = cms.Service("MessageLogger",
    20. 

  20.       debugModules = cms.untracked.vstring('wmnSelFilter'),
    21. 

  21.       cout = cms.untracked.PSet(
    22. 

  22.             default = cms.untracked.PSet( limit = cms.untracked.int32(10) ),
    23. 

  23.             threshold = cms.untracked.string('INFO')
    24. 

  24.             #threshold = cms.untracked.string('DEBUG')
    25. 

  25.       ),
    26. 

  26.       destinations = cms.untracked.vstring('cout')
    27. 

  27. )
    28. 

  28. 

  29. # Selector and parameters
    30. 

  30. process.wmnSelFilter = cms.EDFilter("WMuNuValidator",
    31. 

  31.       # Fast selection flag (no histograms or book-keeping) ->
    32. 

  32.       FastOption = cms.untracked.bool(False),
    33. 

  33. 

  34.       # Input collections ->
    35. 

  35.       TrigTag = cms.untracked.InputTag("TriggerResults::HLT"),
    36. 

  36.       MuonTag = cms.untracked.InputTag("muons"),
    37. 

  37.       METTag = cms.untracked.InputTag("met"),
    38. 

  38.       METIncludesMuons = cms.untracked.bool(False),
    39. 

  39.       JetTag = cms.untracked.InputTag("sisCone5CaloJets"),
    40. 

  40.       
    41. 

  41.       # Main cuts ->
    42. 

  42.       MuonTrig = cms.untracked.string("HLT_Mu9"),
    43. 

  43.       PtCut = cms.untracked.double(25.0),
    44. 

  44.       EtaCut = cms.untracked.double(2.1),
    45. 

  45.       IsRelativeIso = cms.untracked.bool(True),
    46. 

  46.       IsCombinedIso = cms.untracked.bool(False),
    47. 

  47.       IsoCut03 = cms.untracked.double(0.1),
    48. 

  48.       MtMin = cms.untracked.double(50.0),
    49. 

  49.       MtMax = cms.untracked.double(200.0),
    50. 

  50.       MetMin = cms.untracked.double(-999999.),
    51. 

  51.       MetMax = cms.untracked.double(999999.),
    52. 

  52.       AcopCut = cms.untracked.double(2.),
    53. 

  53. 

  54.       # Muon quality cuts ->
    55. 

  55.       DxyCut = cms.untracked.double(0.2),
    56. 

  56.       NormalizedChi2Cut = cms.untracked.double(10.),
    57. 

  57.       TrackerHitsCut = cms.untracked.int32(11),
    58. 

  58.       IsAlsoTrackerMuon = cms.untracked.bool(True),
    59. 

  59.       
    60. 

  60.       # To suppress Zmm ->
    61. 

  61.       PtThrForZ1 = cms.untracked.double(20.0),
    62. 

  62.       PtThrForZ2 = cms.untracked.double(10.0),
    63. 

  63.       
    64. 

  64.       # To further suppress ttbar ->
    65. 

  65.       EJetMin = cms.untracked.double(40.),
    66. 

  66.       NJetMax = cms.untracked.int32(999999)
    67. 

  67. )
    68. 

  68. 

  69. # Output
    70. 

  70. #process.load("Configuration.EventContent.EventContent_cff")
    71. 

  71. #process.wmnOutput = cms.OutputModule("PoolOutputModule",
    72. 

  72. #      process.AODSIMEventContent,
    73. 

  73. #      SelectEvents = cms.untracked.PSet(
    74. 

  74. #            SelectEvents = cms.vstring('wmnsel')
    75. 

  75. #      ),
    76. 

  76. #      fileName = cms.untracked.string('root_files/wmnsel.root')
    77. 

  77. #)
    78. 

  78. 

  79. # Output histograms
    80. 

  80. process.TFileService = cms.Service("TFileService", fileName = cms.string('WMuNu_histograms.root') )
    81. 

  81. 

  82. # Steering the process
    83. 

  83. process.wmnsel = cms.Path(process.wmnSelFilter)
    84. 

  84. #process.end = cms.EndPath(process.wmnOutput)
    85.