/ElectroWeakAnalysis/WENu/python/simpleCutBasedElectronIDSpring10_cfi.py

https://github.com/aivanov-cern/cmssw · Python · 164 lines · 83 code · 9 blank · 72 comment · 0 complexity · 3c8dcdf0043f2a3537197c33e62fc75f MD5 · raw file 

    

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

  2. 

  3. ## Electron ID Based on Simple Cuts: Spring10 MC tuned selections
    4. 

  4. #
    5. 

  5. #  Instructions on how to use this file
    6. 

  6. #  ====================================
    7. 

  7. #
    8. 

  8. #  The selections that are implemented in this python cfg are
    9. 

  9. #  explained in this twiki page:
    10. 

  10. #  https://twiki.cern.ch/twiki/bin/view/CMS/SimpleCutBasedEleID
    11. 

  11. #  In summary, they come in 6 different tightness levels. For
    12. 

  12. #  each tightness, the user can select whether they want
    13. 

  13. #  combined isolation or relative isolations.
    14. 

  14. #
    15. 

  15. #  In order to use this cfg file you have to include it from the
    16. 

  16. #  python directory that you have placed it, clone some selection
    17. 

  17. #  of your preference and run it in your sequence
    18. 

  18. #  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    19. 

  19. #  from ElectroWeakAnalysis.WENu.simpleCutBasedElectronID_cfi import *  
    20. 

  20. #
    21. 

  21. #  simpleEleId_95relIso = simpleCutBasedElectronID.clone()
    22. 

  22. #  simpleEleId_95relIso.electronQuality = '_95relIso_'
    23. 

  23. #  mySequence = cms.Sequence(...+...+..+simpleEleId95relIso+...)
    24. 

  24. #  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    25. 

  25. #  Warning: make sure that you use the correct tags for the
    26. 

  26. #  RecoEgamma/ElectronIdentification package
    27. 

  27. #  consult this twiki to obtain the latest information:
    28. 

  28. #
    29. 

  29. #  https://twiki.cern.ch/twiki/bin/view/CMS/SimpleCutBasedEleID
    30. 

  30. #
    31. 

  31. #  this version of the file needs 
    32. 

  32. #  V00-03-07-03   RecoEgamma/ElectronIdentification
    33. 

  33. 

  34. simpleCutBasedElectronID = cms.EDProducer("EleIdCutBasedExtProducer",
    35. 

  35. 

  36. #   import here your collections
    37. 

  37.     src = cms.InputTag("gsfElectrons"),
    38. 

  38.     #reducedBarrelRecHitCollection = cms.InputTag("ecalRecHit","EcalRecHitsEB"),
    39. 

  39.     #reducedEndcapRecHitCollection = cms.InputTag("ecalRecHit","EcalRecHitsEE"),
    40. 

  40.     # Spring10 uses these names:
    41. 

  41.     reducedBarrelRecHitCollection = cms.InputTag("reducedEcalRecHitsEB"),
    42. 

  42.     reducedEndcapRecHitCollection = cms.InputTag("reducedEcalRecHitsEE"),
    43. 

  43.     # if you want the vertices or the offline beam spot
    44. 

  44.     verticesCollection = cms.InputTag("offlineBeamSpot"),
    45. 

  45.     #dataMagneticFieldSetUp = cms.bool(False),
    46. 

  46.     #dcsTag = cms.InputTag("scalersRawToDigi"),                                          
    47. 

  47.     algorithm = cms.string('eIDCB'),
    48. 

  48. 

  49.     #electronIDType: robust  for the simple Cut-Based
    50. 

  50.     #electronQuality: see later
    51. 

  51.     #electronVersion: use V03 with the offline beam spot
    52. 

  52.     electronIDType  = cms.string('robust'),
    53. 

  53.     electronQuality = cms.string('test'),
    54. 

  54.     electronVersion = cms.string('V04'),
    55. 

  55. 

  56. ####
    57. 

  57. #### Selections with Relative Isolation                                          
    58. 

  58.     robust95relIsoEleIDCutsV04 = cms.PSet(
    59. 

  59.            barrel =  cms.vdouble(1.5e-01, 1.0e-02, 8.0e-01, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 1.5e-01, 
    60. 

  60.                                  2.0e+00, 1.2e-01, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.0, 0.0, ),
    61. 

  61. #           endcap =  cms.vdouble(7.0e-02, 3.0e-02, 7.0e-01, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 8.0e-02, 
    62. 

  62. #                                 6.0e-02, 5.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.0, 0.0, ),
    63. 

  63.            endcap =  cms.vdouble(7.0e-02, 3.0e-02, 7.0e-01, 1.0e-02, -1, -1, 9999., 9999., 9999., 9999., 9999., 8.0e-02, 
    64. 

  64.                                  6.0e-02, 5.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.0, 0.0, ),
    65. 

  65.     ),
    66. 

  66.     robust90relIsoEleIDCutsV04 = cms.PSet(
    67. 

  67.            barrel =  cms.vdouble(1.2e-01, 1.0e-02, 8.0e-01, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 1.2e-01, 
    68. 

  68.                                  9.0e-02, 1.0e-01, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    69. 

  69. #           endcap =  cms.vdouble(5.0e-02, 3.0e-02, 7.0e-01, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 5.0e-02, 
    70. 

  70. #                                 6.0e-02, 3.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),      
    71. 

  71.            endcap =  cms.vdouble(5.0e-02, 3.0e-02, 7.0e-01, 9.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 5.0e-02, 
    72. 

  72.                                  6.0e-02, 3.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    73. 

  73.     ),
    74. 

  74.     robust85relIsoEleIDCutsV04 = cms.PSet(
    75. 

  75.            barrel =  cms.vdouble(4.0e-02, 1.0e-02, 6.0e-02, 6.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9.0e-02, 
    76. 

  76.                                  8.0e-02, 1.0e-01, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    77. 

  77. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 4.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 5.0e-02, 
    78. 

  78. #                                 5.0e-02, 2.5e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    79. 

  79.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 4.0e-02, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 5.0e-02, 
    80. 

  80.                                  5.0e-02, 2.5e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    81. 

  81.     ),
    82. 

  82.     robust80relIsoEleIDCutsV04 = cms.PSet(
    83. 

  83.            barrel =  cms.vdouble(4.0e-02, 1.0e-02, 6.0e-02, 4.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9.0e-02, 
    84. 

  84.                                  7.0e-02, 1.0e-01, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    85. 

  85. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 3.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 4.0e-02, 
    86. 

  86. #                                 5.0e-02, 2.5e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),      
    87. 

  87.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 3.0e-02, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 4.0e-02, 
    88. 

  88.                                  5.0e-02, 2.5e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    89. 

  89.     ),
    90. 

  90.     # 70% point modified with restricting cuts to physical values                                                                   
    91. 

  91.     robust70relIsoEleIDCutsV04 = cms.PSet(
    92. 

  92.            barrel =  cms.vdouble(2.5e-02, 1.0e-02, 3.0e-02, 4.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 5.0e-02, 
    93. 

  93.                                  6.0e-02, 3.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    94. 

  94. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 2.5e-02, 
    95. 

  95. #                                 2.5e-02, 2.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    96. 

  96.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 5.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 2.5e-02, 
    97. 

  97.                                  2.5e-02, 2.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    98. 

  98.     ),
    99. 

  99.     # 60% point modified with restricting cuts to physical values                      
    100. 

  100.     robust60relIsoEleIDCutsV04 = cms.PSet(
    101. 

  101.            barrel =  cms.vdouble(2.5e-02, 1.0e-02, 2.5e-02, 4.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 4.0e-02, 
    102. 

  102.                                  4.0e-02, 3.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    103. 

  103. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 2.5e-02, 
    104. 

  104. #                                 2.0e-02, 2.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),      
    105. 

  105.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 5.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 2.5e-02, 
    106. 

  106.                                  2.0e-02, 2.0e-02, 9999., 9999., 9999., 9999., 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    107. 

  107.     ),
    108. 

  108. 

  109. ####
    110. 

  110. #### Selections with Combined Isolation
    111. 

  111. 

  112.     robust95cIsoEleIDCutsV04 = cms.PSet(
    113. 

  113.            barrel =  cms.vdouble(1.5e-01, 1.0e-02, 8.0e-01, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    114. 

  114.                                  9999., 9999., 9999., 9999., 9999., 1.5e-01, 0.0, -9999., 9999., 9999., 1, -1, 0.0, 0.0, ),
    115. 

  115. #           endcap =  cms.vdouble(7.0e-02, 3.0e-02, 7.0e-01, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    116. 

  116. #                                 9999., 9999., 9999., 9999., 9999., 1.0e-01, 0.0, -9999., 9999., 9999., 1, -1, 0.0, 0.0, ),           
    117. 

  117.            endcap =  cms.vdouble(7.0e-02, 3.0e-02, 7.0e-01, 1.0e-02, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    118. 

  118.                                  9999., 9999., 9999., 9999., 9999., 1.0e-01, 0.0, -9999., 9999., 9999., 1, -1, 0.0, 0.0, ),
    119. 

  119.     ),
    120. 

  120.     robust90cIsoEleIDCutsV04 = cms.PSet(
    121. 

  121.            barrel =  cms.vdouble(1.2e-01, 1.0e-02, 8.0e-01, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    122. 

  122.                                  9999., 9999., 9999., 9999., 9999., 1.0e-01, 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    123. 

  123. #           endcap =  cms.vdouble(5.0e-02, 3.0e-02, 7.0e-01, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    124. 

  124. #                                 9999., 9999., 9999., 9999., 9999., 7.0e-02, 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),        
    125. 

  125.            endcap =  cms.vdouble(5.0e-02, 3.0e-02, 7.0e-01, 9.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    126. 

  126.                                  9999., 9999., 9999., 9999., 9999., 7.0e-02, 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    127. 

  127.     ),
    128. 

  128.     robust85cIsoEleIDCutsV04 = cms.PSet(
    129. 

  129.            barrel =  cms.vdouble(4.0e-02, 1.0e-02, 6.0e-02, 6.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    130. 

  130.                                  9999., 9999., 9999., 9999., 9999., 9.0e-02, 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    131. 

  131. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 4.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    132. 

  132. #                                 9999., 9999., 9999., 9999., 9999., 6.0e-02, 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    133. 

  133.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 4.0e-02, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    134. 

  134.                                  9999., 9999., 9999., 9999., 9999., 6.0e-02, 0.0, -9999., 9999., 9999., 1, -1, 0.02, 0.02, ),
    135. 

  135.     ),
    136. 

  136.     robust80cIsoEleIDCutsV04 = cms.PSet(
    137. 

  137.            barrel =  cms.vdouble(4.0e-02, 1.0e-02, 6.0e-02, 4.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    138. 

  138.                                  9999., 9999., 9999., 9999., 9999., 7.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    139. 

  139. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 3.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    140. 

  140. #                                 9999., 9999., 9999., 9999., 9999., 6.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),        
    141. 

  141.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 3.0e-02, 7.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    142. 

  142.                                  9999., 9999., 9999., 9999., 9999., 6.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    143. 

  143.     ),
    144. 

  144.     # 70% point modified with restricting cuts to physical values                                          
    145. 

  145.     robust70cIsoEleIDCutsV04 = cms.PSet(
    146. 

  146.            barrel =  cms.vdouble(2.5e-02, 1.0e-02, 3.0e-02, 4.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    147. 

  147.                                  9999., 9999., 9999., 9999., 9999., 4.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    148. 

  148. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    149. 

  149. #                                 9999., 9999., 9999., 9999., 9999., 3.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    150. 

  150.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 5.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    151. 

  151.                                  9999., 9999., 9999., 9999., 9999., 3.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    152. 

  152.     ),
    153. 

  153.     # 60% point modified with restricting cuts to physical values
    154. 

  154.     robust60cIsoEleIDCutsV04 = cms.PSet(
    155. 

  155.            barrel =  cms.vdouble(2.5e-02, 1.0e-02, 2.5e-02, 4.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    156. 

  156.                                  9999., 9999., 9999., 9999., 9999., 3.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    157. 

  157. #           endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 9999., -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    158. 

  158. #                                 9999., 9999., 9999., 9999., 9999., 2.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    159. 

  159.            endcap =  cms.vdouble(2.5e-02, 3.0e-02, 2.0e-02, 5.0e-03, -1, -1, 9999., 9999., 9999., 9999., 9999., 9999., 
    160. 

  160.                                  9999., 9999., 9999., 9999., 9999., 2.0e-02, 0.0, -9999., 9999., 9999., 0, -1, 0.02, 0.02, ),
    161. 

  161.     ),
    162. 

  162. 

  163. )
    164.