/ElectroWeakAnalysis/WMuNu/test/WMuNuCandidateHistogrammer.py

https://github.com/aivanov-cern/cmssw · Python · 171 lines · 107 code · 29 blank · 35 comment · 0 complexity · 15243bd7de5441252b5405eab0f19cfc MD5 · raw file 

    

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

  2. 

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

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

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

  6.       input = cms.untracked.int32(-1)
    7. 

  7. )
    8. 

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

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

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

  11.        
    12. 

  12.      # fileNames = cms.untracked.vstring(
    13. 

  13.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_1.root',
    14. 

  14.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_2.root',
    15. 

  15.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_3.root',
    16. 

  16.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_4.root',
    17. 

  17.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_5.root',
    18. 

  18.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_6.root',
    19. 

  19.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_7.root',
    20. 

  20.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_8.root',
    21. 

  21.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_9.root',
    22. 

  22.      #         '/store/user/cepeda/mytestSkim_PTR_Wmunu_10pb/EWK_WMuNu_SubSkim_31Xv3_10.root'
    23. 

  23.      #)                  
    24. 

  24.       
    25. 

  25.      fileNames = cms.untracked.vstring(
    26. 

  26.          "file:EWK_WMuNu_SubSkim_31Xv3.root"
    27. 

  27.       #   "file:AOD_with_WCandidates.root"
    28. 

  28.      )
    29. 

  29. )
    30. 

  30. 

  31. # Debug/info printouts
    32. 

  32. process.MessageLogger = cms.Service("MessageLogger",
    33. 

  33.       debugModules = cms.untracked.vstring('corMetWMuNus','selcorMet'),
    34. 

  34.       cout = cms.untracked.PSet(
    35. 

  35.             default = cms.untracked.PSet( limit = cms.untracked.int32(-1) ),
    36. 

  36.             #threshold = cms.untracked.string('INFO')
    37. 

  37.             threshold = cms.untracked.string('DEBUG')
    38. 

  38.       ),
    39. 

  39.       destinations = cms.untracked.vstring('cout')
    40. 

  40. )
    41. 

  41. 

  42. process.selcorMet = cms.EDFilter("WMuNuSelector",
    43. 

  43.       # Fill Basc Histograms? ->
    44. 

  44.       plotHistograms = cms.untracked.bool(True),
    45. 

  45. 

  46.       # Input collections ->
    47. 

  47.       MuonTag = cms.untracked.InputTag("muons"),
    48. 

  48.       TrigTag = cms.untracked.InputTag("TriggerResults::HLT8E29"),
    49. 

  49.       JetTag = cms.untracked.InputTag("antikt5CaloJets"),
    50. 

  50.       WMuNuCollectionTag = cms.untracked.InputTag("corMetWMuNus"),
    51. 

  51. 

  52.       # Preselection! 
    53. 

  53.       MuonTrig = cms.untracked.string("HLT_Mu9"),
    54. 

  54.       PtThrForZ1 = cms.untracked.double(20.0),
    55. 

  55.       PtThrForZ2 = cms.untracked.double(10.0),
    56. 

  56.       EJetMin = cms.untracked.double(40.),
    57. 

  57.       NJetMax = cms.untracked.int32(999999),
    58. 

  58. 

  59.       # Main cuts ->
    60. 

  60.       PtCut = cms.untracked.double(25.0),
    61. 

  61.       EtaCut = cms.untracked.double(2.1),
    62. 

  62.       IsRelativeIso = cms.untracked.bool(True),
    63. 

  63.       IsCombinedIso = cms.untracked.bool(False),
    64. 

  64.       IsoCut03 = cms.untracked.double(0.1),
    65. 

  65.       MtMin = cms.untracked.double(50.0),
    66. 

  66.       MtMax = cms.untracked.double(200.0),
    67. 

  67.       MetMin = cms.untracked.double(-999999.),
    68. 

  68.       MetMax = cms.untracked.double(999999.),
    69. 

  69.       AcopCut = cms.untracked.double(2.),
    70. 

  70. 

  71.       # Muon quality cuts ->
    72. 

  72.       DxyCut = cms.untracked.double(0.2),
    73. 

  73.       NormalizedChi2Cut = cms.untracked.double(10.),
    74. 

  74.       TrackerHitsCut = cms.untracked.int32(11),
    75. 

  75.       IsAlsoTrackerMuon = cms.untracked.bool(True),
    76. 

  76. 

  77.       # Select only W-, W+ ( default is all Ws)  
    78. 

  78.       SelectByCharge=cms.untracked.int32(0)
    79. 

  79. 

  80. )
    81. 

  81. process.selpfMet = cms.EDFilter("WMuNuSelector",
    82. 

  82.       # Fill Basc Histograms? ->
    83. 

  83.       plotHistograms = cms.untracked.bool(True),
    84. 

  84. 

  85.       # Preselection! 
    86. 

  86.       MuonTrig = cms.untracked.string("HLT_Mu9"),
    87. 

  87.       PtThrForZ1 = cms.untracked.double(20.0),
    88. 

  88.       PtThrForZ2 = cms.untracked.double(10.0),
    89. 

  89.       EJetMin = cms.untracked.double(40.),
    90. 

  90.       NJetMax = cms.untracked.int32(999999),
    91. 

  91. 

  92.       # Input collections ->
    93. 

  93.       MuonTag = cms.untracked.InputTag("muons"),
    94. 

  94.       TrigTag = cms.untracked.InputTag("TriggerResults::HLT8E29"),
    95. 

  95.       JetTag = cms.untracked.InputTag("antikt5CaloJets"),
    96. 

  96.       WMuNuCollectionTag = cms.untracked.InputTag("pfMetWMuNus"),
    97. 

  97. 

  98.       # Main cuts ->
    99. 

  99.       UseTrackerPt = cms.untracked.bool(True),
    100. 

  100.       PtCut = cms.untracked.double(25.0),
    101. 

  101.       EtaCut = cms.untracked.double(2.1),
    102. 

  102.       IsRelativeIso = cms.untracked.bool(True),
    103. 

  103.       IsCombinedIso = cms.untracked.bool(False),
    104. 

  104.       IsoCut03 = cms.untracked.double(0.1),
    105. 

  105.       MtMin = cms.untracked.double(50.0),
    106. 

  106.       MtMax = cms.untracked.double(200.0),
    107. 

  107.       MetMin = cms.untracked.double(-999999.),
    108. 

  108.       MetMax = cms.untracked.double(999999.),
    109. 

  109.       AcopCut = cms.untracked.double(2.),
    110. 

  110. 

  111.       # Muon quality cuts ->
    112. 

  112.       DxyCut = cms.untracked.double(0.2),
    113. 

  113.       NormalizedChi2Cut = cms.untracked.double(10.),
    114. 

  114.       TrackerHitsCut = cms.untracked.int32(11),
    115. 

  115.       IsAlsoTrackerMuon = cms.untracked.bool(True),
    116. 

  116. 

  117.       # Select only W-, W+ ( default is all Ws)
    118. 

  118.       SelectByCharge=cms.untracked.int32(0)
    119. 

  119. 

  120. )
    121. 

  121. process.seltcMet = cms.EDFilter("WMuNuSelector",
    122. 

  122.       # Fill Basc Histograms? ->
    123. 

  123.       plotHistograms = cms.untracked.bool(True),
    124. 

  124. 

  125.       # Input collections ->
    126. 

  126.       MuonTag = cms.untracked.InputTag("muons"),
    127. 

  127.       TrigTag = cms.untracked.InputTag("TriggerResults::HLT8E29"),
    128. 

  128.       JetTag = cms.untracked.InputTag("antikt5CaloJets"),
    129. 

  129.       WMuNuCollectionTag = cms.untracked.InputTag("tcMetWMuNus"),
    130. 

  130. 

  131.       # Preselection! 
    132. 

  132.       MuonTrig = cms.untracked.string("HLT_Mu9"),
    133. 

  133.       PtThrForZ1 = cms.untracked.double(20.0),
    134. 

  134.       PtThrForZ2 = cms.untracked.double(10.0),
    135. 

  135.       EJetMin = cms.untracked.double(40.),
    136. 

  136.       NJetMax = cms.untracked.int32(999999),
    137. 

  137. 

  138.       # Main cuts ->
    139. 

  139.       UseTrackerPt = cms.untracked.bool(True),
    140. 

  140.       PtCut = cms.untracked.double(25.0),
    141. 

  141.       EtaCut = cms.untracked.double(2.1),
    142. 

  142.       IsRelativeIso = cms.untracked.bool(True),
    143. 

  143.       IsCombinedIso = cms.untracked.bool(False),
    144. 

  144.       IsoCut03 = cms.untracked.double(0.1),
    145. 

  145.       MtMin = cms.untracked.double(50.0),
    146. 

  146.       MtMax = cms.untracked.double(200.0),
    147. 

  147.       MetMin = cms.untracked.double(-999999.),
    148. 

  148.       MetMax = cms.untracked.double(999999.),
    149. 

  149.       AcopCut = cms.untracked.double(2.),
    150. 

  150. 

  151.       # Muon quality cuts ->
    152. 

  152.       DxyCut = cms.untracked.double(0.2),
    153. 

  153.       NormalizedChi2Cut = cms.untracked.double(10.),
    154. 

  154.       TrackerHitsCut = cms.untracked.int32(11),
    155. 

  155.       IsAlsoTrackerMuon = cms.untracked.bool(True),
    156. 

  156. 

  157.       # Select only W-, W+ ( default is all Ws)
    158. 

  158.       SelectByCharge=cms.untracked.int32(0)
    159. 

  159. 

  160. )
    161. 

  161. 

  162. process.TFileService = cms.Service("TFileService", fileName = cms.string('WMuNuBasicPlots.root') )
    163. 

  163. 

  164. 

  165. # Steering the process
    166. 

  166. process.path1 = cms.Path(process.selcorMet)
    167. 

  167. process.path2 = cms.Path(process.selpfMet)
    168. 

  168. process.path3 = cms.Path(process.seltcMet)
    169.