/ElectroWeakAnalysis/ZEE/test/pat_ZeeVBTF_ntuple_woCertJSON_forData.py

https://github.com/aivanov-cern/cmssw · Python · 469 lines · 253 code · 90 blank · 126 comment · 3 complexity · ab8aea1c8def169eb6aa4ed75293e93b MD5 · raw file 

    

  1. #   Initializations
    2. 

  2. #   ---------------
    3. 

  3. #
    4. 

  4. #   Get the base name of the Python script and use it to define input/output associated files
    5. 

  5. #
    6. 

  6. #   N.B:    Special variable __file__ is not available within cmsRun, and the following line fails:
    7. 

  7. #
    8. 

  8. #               pyFile=(__file__).replace('.py','')
    9. 

  9. #
    10. 

  10. #           As a result, the calling arguments must be scanned to locate the Python script.
    11. 

  11. #
    12. 

  12. import sys
    13. 

  13. 

  14. for inArg in sys.argv:
    15. 

  15.     if inArg.find('.py') == -1:
    16. 

  16.         continue
    17. 

  17.     else:
    18. 

  18.         pyBaseName=inArg.replace('.py','')
    19. 

  19. 

  20. #   Define input/output associated files.
    21. 

  21. #
    22. 

  22. pyFile = "%s.py" % (pyBaseName)             #   Name of the Python script itself
    23. 

  23. lstFile = "%s.list" % (pyBaseName)          #   Name of the list file       (contains a list of files to process)
    24. 

  24. resFile = "%s.results" % (pyBaseName)       #   Name of the results file    (stdout and stderr messages)
    25. 

  25. outFile = "%s.root" % (pyBaseName)          #   Name of the output file     (ROOT)
    26. 

  26. 

  27. 

  28. ##   Configuration for the production of the ICHEP VBTF ntuple
    29. 

  29. ##   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    30. 

  30. ##   MC, spring10
    31. 

  31. ##
    32. 

  32. ##   Stilianos Kesisoglou - Institute of Nuclear Physics
    33. 

  33. ##                                NCSR Demokritos
    34. 

  34. ##   25 June 2010
    35. 

  35. import FWCore.ParameterSet.Config as cms
    36. 

  36. 

  37. process = cms.Process("PAT")
    38. 

  38. 

  39. 

  40. 

  41. process.options = cms.untracked.PSet(
    42. 

  42.     Rethrow = cms.untracked.vstring('ProductNotFound')
    43. 

  43. )
    44. 

  44. 

  45. #process.MessageLogger = cms.Service(
    46. 

  46. #        "MessageLogger",
    47. 

  47. #            categories = cms.untracked.vstring('info', 'debug','cout')
    48. 

  48. #            )
    49. 

  49. 

  50. 

  51. process.load("FWCore.MessageService.MessageLogger_cfi")
    52. 

  52. 

  53. process.MessageLogger.cerr.threshold = cms.untracked.string("INFO")
    54. 

  54. 

  55. process.MessageLogger.cerr.FwkSummary = cms.untracked.PSet(
    56. 

  56.     reportEvery = cms.untracked.int32(1000000),
    57. 

  57.     limit = cms.untracked.int32(10000000)
    58. 

  58.        )
    59. 

  59.        
    60. 

  60. process.MessageLogger.cerr.FwkReport = cms.untracked.PSet(
    61. 

  61.       reportEvery = cms.untracked.int32(100000),
    62. 

  62.          limit = cms.untracked.int32(10000000)
    63. 

  63.       )
    64. 

  64.       
    65. 

  65. process.options = cms.untracked.PSet(
    66. 

  66.        wantSummary = cms.untracked.bool(True)
    67. 

  67.        )
    68. 

  68. 

  69. # 
    70. 

  70. # # source
    71. 

  71. # process.source = cms.Source("PoolSource", 
    72. 

  72. #      fileNames = cms.untracked.vstring(
    73. 

  73. # # SOME DATA FILE TO BE PUT HERE
    74. 

  74. # #
    75. 

  75. # #   DATA test (local) running on the OneElecPlusSC skim made from the EG DATA samples
    76. 

  76. # #   ----------------------------------------------------------------------------------
    77. 

  77. # #    'rfio:/tmp/ikesisog/TestFiles/skimDATA/Zee/00000000-0000-0000-0000-000000000000_woFakeJSON.root',
    78. 

  78. # #    'rfio:/tmp/ikesisog/TestFiles/skimDATA/Zee/00000000-0000-0000-0000-000000000000_wiFakeJSON.root',
    79. 

  79. # #
    80. 

  80. # #   DATA test (local) running directly on the EG DATA samples
    81. 

  81. # #   ----------------------------------------------------------
    82. 

  82. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/6AD8C6D3-2C92-DF11-AFDE-0030487C90C2.root',
    83. 

  83. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/34BD6E1C-3892-DF11-BD3F-001D09F27067.root',
    84. 

  84. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/04D0FF1D-2792-DF11-8B6C-003048F11942.root',
    85. 

  85. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/74CF44CB-1592-DF11-989E-001617C3B77C.root',
    86. 

  86. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/66072495-1A92-DF11-9461-0030487BC68E.root',
    87. 

  87. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/5E06C91D-1E92-DF11-9591-0016177CA778.root',
    88. 

  88. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/48B8741C-1E92-DF11-ADAD-001617C3B6DC.root',
    89. 

  89. #    'rfio:/tmp/ikesisog/TestFiles/dSetDATA/Zee/1EC74696-1A92-DF11-B9BB-0030487CD13A.root',
    90. 

  90. #     )
    91. 

  91. # )
    92. 

  92. # 
    93. 

  93. 

  94. 

  95. #   source - Use an input list to bypass the 255 file limit
    96. 

  96. import FWCore.Utilities.FileUtils as FileUtils
    97. 

  97. 

  98. filelist = FileUtils.loadListFromFile(lstFile) 
    99. 

  99. 

  100. process.source = cms.Source('PoolSource',
    101. 

  101.     fileNames = cms.untracked.vstring(*filelist)
    102. 

  102.     )
    103. 

  103. 

  104. 

  105. # 
    106. 

  106. # #   using locally provided JSON
    107. 

  107. # import PhysicsTools.PythonAnalysis.LumiList as LumiList
    108. 

  108. # import FWCore.ParameterSet.Types as CfgTypes
    109. 

  109. # myLumis = LumiList.LumiList(filename = 'Cert_132440-149442_7TeV_StreamExpress_Collisions10_JSON_v2.json').getCMSSWString().split(',')
    110. 

  110. # process.source.lumisToProcess = CfgTypes.untracked(CfgTypes.VLuminosityBlockRange())
    111. 

  111. # process.source.lumisToProcess.extend(myLumis)
    112. 

  112. # 
    113. 

  113. 

  114. process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
    115. 

  115. 

  116. ## Load additional processes
    117. 

  117. # 
    118. 

  118. process.load("Configuration.StandardSequences.Geometry_cff")
    119. 

  119. process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
    120. 

  120. 

  121. ## global tags:
    122. 

  122. #
    123. 

  123. process.GlobalTag.globaltag = cms.string('GR_R_39X_V5::All') # GLOBAL TAG FOR DATA
    124. 

  124. 

  125. process.load("Configuration.StandardSequences.MagneticField_cff")
    126. 

  126. 

  127. 

  128. ################################################################################################
    129. 

  129. ###    P r e p a r a t i o n      o f    t h e    P A T    O b j e c t s   f r o m    A O D  ###
    130. 

  130. ################################################################################################
    131. 

  131. 

  132. ## pat sequences to be loaded:
    133. 

  133. #process.load("PhysicsTools.PFCandProducer.PF2PAT_cff")
    134. 

  134. 

  135. process.load("PhysicsTools.PatAlgos.patSequences_cff")
    136. 

  136. 

  137. #process.load("PhysicsTools.PatAlgos.triggerLayer1.triggerProducer_cff")
    138. 

  138. 

  139. ##
    140. 

  140. #
    141. 

  141. 

  142. ## MET creation     <=== WARNING: YOU MAY WANT TO MODIFY THIS PART OF THE CODE       %%%%%%%%%%%%%
    143. 

  143. ##                                specify the names of the MET collections that you need here %%%%
    144. 

  144. ##                                                                                             #%%
    145. 

  145. ## if you don't specify anything the default MET is the raw Calo MET                           #%%
    146. 

  146. process.caloMET = process.patMETs.clone(                                                       #%%
    147. 

  147.     metSource = cms.InputTag("met","","RECO"),
    148. 

  148.     addTrigMatch = cms.bool(False),
    149. 

  149.     addMuonCorrections = cms.bool(False),
    150. 

  150.     addGenMET = cms.bool(False),
    151. 

  151. )
    152. 

  152. process.tcMET = process.patMETs.clone(                                                         #%%
    153. 

  153.     metSource = cms.InputTag("tcMet","","RECO"),
    154. 

  154.     addTrigMatch = cms.bool(False),
    155. 

  155.     addMuonCorrections = cms.bool(False),
    156. 

  156.     addGenMET = cms.bool(False),
    157. 

  157. )
    158. 

  158. process.pfMET = process.patMETs.clone(                                                         #%%
    159. 

  159.     metSource = cms.InputTag("pfMet","","RECO"),
    160. 

  160.     addTrigMatch = cms.bool(False),
    161. 

  161.     addMuonCorrections = cms.bool(False),
    162. 

  162.     addGenMET = cms.bool(False),
    163. 

  163. )
    164. 

  164. 

  165. ## specify here what you want to have on the plots! <===== MET THAT YOU WANT ON THE PLOTS  %%%%%%%
    166. 

  166. myMetCollection   = 'caloMET'
    167. 

  167. myPfMetCollection =   'pfMET'
    168. 

  168. myTcMetCollection =   'tcMET'
    169. 

  169. 

  170. ## modify the sequence of the MET creation:                                                    #%%
    171. 

  171. process.makePatMETs = cms.Sequence(process.caloMET*process.tcMET*process.pfMET)
    172. 

  172. 

  173. 

  174. ## modify the final pat sequence: keep only electrons + METS (muons are needed for met corrections)
    175. 

  175. 

  176. process.load("RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequence_cff")
    177. 

  177. #process.patElectronIsolation = cms.Sequence(process.egammaIsolationSequence)
    178. 

  178. 

  179. process.patElectrons.isoDeposits = cms.PSet()
    180. 

  180. process.patElectrons.userIsolation = cms.PSet()
    181. 

  181. process.patElectrons.addElectronID = cms.bool(True)
    182. 

  182. 

  183. process.patElectrons.electronIDSources = cms.PSet(
    184. 

  184.     simpleEleId95relIso = cms.InputTag("simpleEleId95relIso"),
    185. 

  185.     simpleEleId90relIso = cms.InputTag("simpleEleId90relIso"),
    186. 

  186.     simpleEleId85relIso = cms.InputTag("simpleEleId85relIso"),
    187. 

  187.     simpleEleId80relIso = cms.InputTag("simpleEleId80relIso"),
    188. 

  188.     simpleEleId70relIso = cms.InputTag("simpleEleId70relIso"),
    189. 

  189.     simpleEleId60relIso = cms.InputTag("simpleEleId60relIso"),
    190. 

  190.     simpleEleId95cIso   = cms.InputTag("simpleEleId95cIso"),
    191. 

  191.     simpleEleId90cIso   = cms.InputTag("simpleEleId90cIso"),
    192. 

  192.     simpleEleId85cIso   = cms.InputTag("simpleEleId85cIso"),
    193. 

  193.     simpleEleId80cIso   = cms.InputTag("simpleEleId80cIso"),
    194. 

  194.     simpleEleId70cIso   = cms.InputTag("simpleEleId70cIso"),
    195. 

  195.     simpleEleId60cIso   = cms.InputTag("simpleEleId60cIso")    
    196. 

  196. )
    197. 

  197. 

  198. ##
    199. 

  199. 

  200. process.patElectrons.addGenMatch = cms.bool(False)
    201. 

  201. process.patElectrons.embedGenMatch = cms.bool(False)
    202. 

  202. process.patElectrons.usePV = cms.bool(False)
    203. 

  203. 

  204. ##
    205. 

  205. 

  206. process.load("ElectroWeakAnalysis.ZEE.simpleEleIdSequence_cff")
    207. 

  207. 

  208. # you have to tell the ID that it is data. These are set to False for MC
    209. 

  209. process.simpleEleId95relIso.dataMagneticFieldSetUp = cms.bool(True)
    210. 

  210. process.simpleEleId90relIso.dataMagneticFieldSetUp = cms.bool(True)
    211. 

  211. process.simpleEleId85relIso.dataMagneticFieldSetUp = cms.bool(True)
    212. 

  212. process.simpleEleId80relIso.dataMagneticFieldSetUp = cms.bool(True)
    213. 

  213. process.simpleEleId70relIso.dataMagneticFieldSetUp = cms.bool(True)
    214. 

  214. process.simpleEleId60relIso.dataMagneticFieldSetUp = cms.bool(True)
    215. 

  215. process.simpleEleId95cIso.dataMagneticFieldSetUp = cms.bool(True)
    216. 

  216. process.simpleEleId90cIso.dataMagneticFieldSetUp = cms.bool(True)
    217. 

  217. process.simpleEleId85cIso.dataMagneticFieldSetUp = cms.bool(True)
    218. 

  218. process.simpleEleId80cIso.dataMagneticFieldSetUp = cms.bool(True)
    219. 

  219. process.simpleEleId70cIso.dataMagneticFieldSetUp = cms.bool(True)
    220. 

  220. process.simpleEleId60cIso.dataMagneticFieldSetUp = cms.bool(True)
    221. 

  221. #
    222. 

  222. process.patElectronIDs = cms.Sequence(process.simpleEleIdSequence)
    223. 

  223. 

  224. process.makePatElectrons = cms.Sequence(process.patElectronIDs*process.patElectrons)
    225. 

  225. 

  226. # process.makePatMuons may be needed depending on how you calculate the MET
    227. 

  227. 

  228. process.makePatCandidates = cms.Sequence(process.makePatElectrons+process.makePatMETs)
    229. 

  229. 

  230. process.patDefaultSequence = cms.Sequence(process.makePatCandidates)
    231. 

  231. 

  232. ## WARNING: you may want to modify this item:
    233. 

  233. HLT_process_name = "HLT"   # REDIGI for the production traditional MC / HLT for the powheg samples or data
    234. 

  234. 

  235. # Trigger Path(s)
    236. 

  236. HLT_path_name     = "HLT_Photon10_L1R"
    237. 

  237. 

  238. # Label of the last Trigger Filter in the Trigger Path
    239. 

  239. HLT_filter_name  =  "hltL1NonIsoHLTNonIsoSinglePhotonEt10HcalIsolFilter"
    240. 

  240. 

  241. HLT_path_name_extra0   = "HLT_Photon15_Cleaned_L1R"
    242. 

  242. HLT_filter_name_extra0 = cms.untracked.InputTag("hltL1NonIsoHLTNonIsoSinglePhotonEt15CleanedHcalIsolFilter","",HLT_process_name)
    243. 

  243. 

  244. HLT_path_name_extra1   = "HLT_Ele15_SW_CaloEleId_L1R"
    245. 

  245. HLT_filter_name_extra1 = cms.untracked.InputTag("hltL1NonIsoHLTNonIsoSingleElectronEt15CaloEleIdPixelMatchFilter","",HLT_process_name)
    246. 

  246. 

  247. HLT_path_name_extra2   = "HLT_Ele17_SW_CaloEleId_L1R"
    248. 

  248. HLT_filter_name_extra2 = cms.untracked.InputTag("hltL1NonIsoHLTNonIsoSingleElectronEt17CaloEleIdPixelMatchFilter","",HLT_process_name)
    249. 

  249. 

  250. HLT_path_name_extra3   = "HLT_Ele17_SW_TightEleId_L1R"
    251. 

  251. HLT_filter_name_extra3 = cms.untracked.InputTag("hltL1NonIsoHLTNonIsoSingleElectronEt17TightEleIdDphiFilter","",HLT_process_name)
    252. 

  252. 

  253. # HLT_path_name_extra4   = "HLT_Ele17_SW_TighterEleIdIsol_L1R_v2"
    254. 

  254. # HLT_filter_name_extra4 = cms.untracked.InputTag("hltL1NonIsoHLTNonIsoSingleElectronEt17TighterEleIdIsolTrackIsolFilter","",HLT_process_name)
    255. 

  255. 

  256. HLT_path_name_extra5   = "HLT_Ele22_SW_TighterCaloIdIsol_L1R_v1"
    257. 

  257. HLT_filter_name_extra5 = cms.untracked.InputTag("hltL1NonIsoHLTNonIsoSingleElectronEt22TighterCaloIdIsolTrackIsolFilter","",HLT_process_name)
    258. 

  258. 

  259. HLT_path_name_extra6   = "HLT_Ele22_SW_TighterCaloIdIsol_L1R_v2"
    260. 

  260. HLT_filter_name_extra6 = cms.untracked.InputTag("hltL1NonIsoHLTNonIsoSingleElectronEt22TighterCaloIdIsolTrackIsolFilter","",HLT_process_name)
    261. 

  261. 

  262. 

  263. rules_Filter = cms.PSet (
    264. 

  264.     ### the input collections needed:
    265. 

  265.     electronCollectionTag = cms.untracked.InputTag("patElectrons","","PAT"),
    266. 

  266.     metCollectionTag = cms.untracked.InputTag(myMetCollection,"","PAT"),
    267. 

  267.     pfMetCollectionTag = cms.untracked.InputTag(myPfMetCollection,"","PAT"),
    268. 

  268.     tcMetCollectionTag = cms.untracked.InputTag(myTcMetCollection,"","PAT"),
    269. 

  269.     triggerCollectionTag = cms.untracked.InputTag("TriggerResults","",HLT_process_name),
    270. 

  270.     triggerEventTag = cms.untracked.InputTag("hltTriggerSummaryAOD","",HLT_process_name),
    271. 

  271.     hltpath = cms.untracked.string(HLT_path_name), 
    272. 

  272.     hltpathFilter = cms.untracked.InputTag(HLT_filter_name,"",HLT_process_name),
    273. 

  273.     ebRecHits = cms.untracked.InputTag("reducedEcalRecHitsEB"),
    274. 

  274.     eeRecHits = cms.untracked.InputTag("reducedEcalRecHitsEE"),
    275. 

  275.     PrimaryVerticesCollection = cms.untracked.InputTag("offlinePrimaryVertices"),
    276. 

  276.     ### here the preselection is applied
    277. 

  277.     # fiducial cuts:
    278. 

  278.     BarrelMaxEta = cms.untracked.double(1.4442),
    279. 

  279.     EndCapMinEta = cms.untracked.double(1.5660),
    280. 

  280.     EndCapMaxEta = cms.untracked.double(2.5000),
    281. 

  281.     # demand ecal driven electron:
    282. 

  282.     useEcalDrivenElectrons = cms.untracked.bool(True),
    283. 

  283.     # demand offline spike cleaning with the Swiss Cross criterion:
    284. 

  284.     useSpikeRejection = cms.untracked.bool(False),
    285. 

  285.     spikeCleaningSwissCrossCut = cms.untracked.double(0.95),
    286. 

  286.     # demand geometrically matched to an HLT object with ET>15GeV
    287. 

  287.     useTriggerInfo = cms.untracked.bool(True),
    288. 

  288.     electronMatched2HLT = cms.untracked.bool(True),
    289. 

  289.     electronMatched2HLT_DR = cms.untracked.double(0.1),
    290. 

  290.     useHLTObjectETCut = cms.untracked.bool(True),
    291. 

  291.     hltObjectETCut = cms.untracked.double(15.0),
    292. 

  292.     useExtraTrigger = cms.untracked.bool(True),
    293. 

  293. #     vHltpathExtra = cms.untracked.vstring(HLT_path_name_extra0,HLT_path_name_extra1,HLT_path_name_extra2,HLT_path_name_extra3,HLT_path_name_extra4,HLT_path_name_extra5),
    294. 

  294. #     vHltpathFilterExtra = cms.untracked.VInputTag(HLT_filter_name_extra0,HLT_filter_name_extra1,HLT_filter_name_extra2,HLT_filter_name_extra3,HLT_filter_name_extra4,HLT_filter_name_extra5),
    295. 

  295.     vHltpathExtra = cms.untracked.vstring(HLT_path_name_extra0,HLT_path_name_extra1,HLT_path_name_extra2,HLT_path_name_extra3,HLT_path_name_extra5,HLT_path_name_extra6),
    296. 

  296.     vHltpathFilterExtra = cms.untracked.VInputTag(HLT_filter_name_extra0,HLT_filter_name_extra1,HLT_filter_name_extra2,HLT_filter_name_extra3,HLT_filter_name_extra5,HLT_filter_name_extra6),
    297. 

  297. 

  298.     # ET Cut in the SC
    299. 

  299.     ETCut = cms.untracked.double(25.0),                                  
    300. 

  300.     METCut = cms.untracked.double(0.0),
    301. 

  301.     # For DATA set it to True, for MC set it to False
    302. 

  302.     dataMagneticFieldSetUp = cms.untracked.bool(True),
    303. 

  303.     dcsTag = cms.untracked.InputTag("scalersRawToDigi")
    304. 

  304. )
    305. 

  305. 

  306. 

  307. rules_Filter_Elec0 = cms.PSet (
    308. 

  308.     # Other parameters of the code - leave them as they are
    309. 

  309.     useValidFirstPXBHit0            = cms.untracked.bool(False),
    310. 

  310.     useConversionRejection0         = cms.untracked.bool(False),
    311. 

  311.     useExpectedMissingHits0         = cms.untracked.bool(True),
    312. 

  312.     maxNumberOfExpectedMissingHits0 = cms.untracked.int32(0),
    313. 

  313.     # calculate some new cuts
    314. 

  314.     calculateValidFirstPXBHit0      = cms.untracked.bool(False),
    315. 

  315.     calculateConversionRejection0   = cms.untracked.bool(False),
    316. 

  316.     calculateExpectedMissingHits0   = cms.untracked.bool(False)
    317. 

  317. )
    318. 

  318. 

  319. rules_Filter_Elec1 = cms.PSet (
    320. 

  320.     # Other parameters of the code - leave them as they are
    321. 

  321.     useValidFirstPXBHit1            = cms.untracked.bool(False),
    322. 

  322.     useConversionRejection1         = cms.untracked.bool(False),
    323. 

  323.     useExpectedMissingHits1         = cms.untracked.bool(True),
    324. 

  324.     maxNumberOfExpectedMissingHits1 = cms.untracked.int32(0),
    325. 

  325.     # calculate some new cuts
    326. 

  326.     calculateValidFirstPXBHit1      = cms.untracked.bool(False),
    327. 

  327.     calculateConversionRejection1   = cms.untracked.bool(False),
    328. 

  328.     calculateExpectedMissingHits1   = cms.untracked.bool(False)
    329. 

  329. )
    330. 

  330. 

  331. rules_Filter_Elec2 = cms.PSet (
    332. 

  332.     # Other parameters of the code - leave them as they are
    333. 

  333.     useValidFirstPXBHit2            = cms.untracked.bool(False),
    334. 

  334.     useConversionRejection2         = cms.untracked.bool(False),
    335. 

  335.     useExpectedMissingHits2         = cms.untracked.bool(True),
    336. 

  336.     maxNumberOfExpectedMissingHits2 = cms.untracked.int32(0),
    337. 

  337.     # calculate some new cuts
    338. 

  338.     calculateValidFirstPXBHit2      = cms.untracked.bool(False),
    339. 

  339.     calculateConversionRejection2   = cms.untracked.bool(False),
    340. 

  340.     calculateExpectedMissingHits2   = cms.untracked.bool(False)
    341. 

  341. )
    342. 

  342. 

  343. 

  344. process.zeeFilter = cms.EDFilter('ZeeCandidateFilter',
    345. 

  345.     rules_Filter,
    346. 

  346.     rules_Filter_Elec0,    
    347. 

  347.     rules_Filter_Elec1,    
    348. 

  348.     rules_Filter_Elec2   
    349. 

  349. )
    350. 

  350. 

  351. 

  352. 

  353. ## the Z selection that you prefer
    354. 

  354. from ElectroWeakAnalysis.ZEE.simpleCutBasedSpring10SelectionBlocks_cfi import *
    355. 

  355. 

  356. selection_inverse = cms.PSet (
    357. 

  357.     deta_EB_inv = cms.untracked.bool(True),
    358. 

  358.     deta_EE_inv = cms.untracked.bool(True)
    359. 

  359.     )
    360. 

  360. 

  361. 

  362. rules_Plotter = cms.PSet (
    363. 

  363.     #   treat or not the elecrons with same criteria. if yes then rules_Plotter_Elec1/2 are neglected
    364. 

  364.     useSameSelectionOnBothElectrons = cms.untracked.bool(True),    
    365. 

  365.     # some extra information on the ntuple production:
    366. 

  366.     sihih0_EB = cms.untracked.double(1000.0),
    367. 

  367.     dphi0_EB  = cms.untracked.double(1000.0),
    368. 

  368.     deta0_EB  = cms.untracked.double(1000.0),
    369. 

  369.     hoe0_EB   = cms.untracked.double(1000.0),
    370. 

  370.     sihih0_EE = cms.untracked.double(1000.0),
    371. 

  371.     dphi0_EE  = cms.untracked.double(1000.0),
    372. 

  372.     deta0_EE  = cms.untracked.double(1000.0),
    373. 

  373.     hoe0_EE   = cms.untracked.double(1000.0),
    374. 

  374.     includeJetInformationInNtuples = cms.untracked.bool(True),
    375. 

  375.     caloJetCollectionTag = cms.untracked.InputTag('ak5CaloJetsL2L3'),
    376. 

  376.     pfJetCollectionTag = cms.untracked.InputTag('ak5PFJetsL2L3'),
    377. 

  377.     DRJetFromElectron = cms.untracked.double(0.3),
    378. 

  378.     #
    379. 

  379.     zeeCollectionTag = cms.untracked.InputTag("zeeFilter","selectedZeeCandidates","PAT"),
    380. 

  380.     ZEE_VBTFselectionFileName = cms.untracked.string("ZEE_VBTFselection.root"),
    381. 

  381.     ZEE_VBTFpreseleFileName = cms.untracked.string("ZEE_VBTFpreselection.root"),
    382. 

  382.     DatasetTag =  cms.untracked.int32(100),
    383. 

  383.     storeExtraInformation = cms.untracked.bool(True)
    384. 

  384. )
    385. 

  385. 

  386. rules_Plotter_Elec0 = cms.PSet (
    387. 

  387.     # The selection to be used here:
    388. 

  388.     usePrecalcID0                   = cms.untracked.bool(True),
    389. 

  389.     usePrecalcIDType0               = cms.untracked.string('simpleEleId80relIso'),
    390. 

  390.     usePrecalcIDSign0               = cms.untracked.string('='),
    391. 

  391.     usePrecalcIDValue0              = cms.untracked.double(7),    
    392. 

  392.     ## preselection criteria are independent of useSameSelectionOnBothElectrons
    393. 

  393.     #  set them to False if you don't want them
    394. 

  394.     useValidFirstPXBHit0            = cms.untracked.bool(False),
    395. 

  395.     useConversionRejection0         = cms.untracked.bool(False),
    396. 

  396.     useExpectedMissingHits0         = cms.untracked.bool(True),
    397. 

  397.     maxNumberOfExpectedMissingHits0 = cms.untracked.int32(0),    
    398. 

  398.     ##
    399. 

  399.     calculateValidFirstPXBHit0      = cms.untracked.bool(False),
    400. 

  400.     calculateConversionRejection0   = cms.untracked.bool(False),
    401. 

  401.     calculateExpectedMissingHits0   = cms.untracked.bool(False)
    402. 

  402. )
    403. 

  403. 

  404. rules_Plotter_Elec1 = cms.PSet (
    405. 

  405.     # The selection to be used here:
    406. 

  406.     usePrecalcID1                   = cms.untracked.bool(True),
    407. 

  407.     usePrecalcIDType1               = cms.untracked.string('simpleEleId80relIso'),
    408. 

  408.     usePrecalcIDSign1               = cms.untracked.string('='),
    409. 

  409.     usePrecalcIDValue1              = cms.untracked.double(7),    
    410. 

  410.     ## preselection criteria are independent of useSameSelectionOnBothElectrons
    411. 

  411.     #  set them to False if you don't want them
    412. 

  412.     useValidFirstPXBHit1            = cms.untracked.bool(False),
    413. 

  413.     useConversionRejection1         = cms.untracked.bool(False),
    414. 

  414.     useExpectedMissingHits1         = cms.untracked.bool(True),
    415. 

  415.     maxNumberOfExpectedMissingHits1 = cms.untracked.int32(0),    
    416. 

  416.     ##
    417. 

  417.     calculateValidFirstPXBHit1      = cms.untracked.bool(False),
    418. 

  418.     calculateConversionRejection1   = cms.untracked.bool(False),
    419. 

  419.     calculateExpectedMissingHits1   = cms.untracked.bool(False)
    420. 

  420. )
    421. 

  421. 

  422. rules_Plotter_Elec2 = cms.PSet (
    423. 

  423.     # The selection to be used here:
    424. 

  424.     usePrecalcID2                   = cms.untracked.bool(True),
    425. 

  425.     usePrecalcIDType2               = cms.untracked.string('simpleEleId80relIso'),
    426. 

  426.     usePrecalcIDSign2               = cms.untracked.string('='),
    427. 

  427.     usePrecalcIDValue2              = cms.untracked.double(7),    
    428. 

  428.     ## preselection criteria are independent of useSameSelectionOnBothElectrons
    429. 

  429.     #  set them to False if you don't want them
    430. 

  430.     useValidFirstPXBHit2            = cms.untracked.bool(False),
    431. 

  431.     useConversionRejection2         = cms.untracked.bool(False),
    432. 

  432.     useExpectedMissingHits2         = cms.untracked.bool(True),
    433. 

  433.     maxNumberOfExpectedMissingHits2 = cms.untracked.int32(0),    
    434. 

  434.     ##
    435. 

  435.     calculateValidFirstPXBHit2      = cms.untracked.bool(False),
    436. 

  436.     calculateConversionRejection2   = cms.untracked.bool(False),
    437. 

  437.     calculateExpectedMissingHits2   = cms.untracked.bool(False)
    438. 

  438. )
    439. 

  439. 

  440. 

  441. # we need to store jet information, hence we have to produce the jets:
    442. 

  442. process.load("JetMETCorrections.Configuration.DefaultJEC_cff")
    443. 

  443. 

  444. process.jetSequence = cms.Sequence( process.ak5CaloJetsL2L3  )
    445. 

  445. 

  446. process.pfjetAK5Sequence = cms.Sequence( process.ak5PFJetsL2L3 )
    447. 

  447. 

  448. process.ourJetSequence = cms.Sequence( process.jetSequence * process.pfjetAK5Sequence )
    449. 

  449. 

  450. 

  451. process.plotter = cms.EDAnalyzer('ZeePlots',
    452. 

  452.     rules_Plotter,
    453. 

  453.     rules_Plotter_Elec0,
    454. 

  454.     rules_Plotter_Elec1,
    455. 

  455.     rules_Plotter_Elec2   
    456. 

  456. )
    457. 

  457. 

  458. 

  459. # For DATA, 361p4 electrons need fixing first (misalignment corrections) *NOT* be used with MC
    460. 

  460. #process.load("RecoEgamma.EgammaTools.correctedElectronsProducer_cfi")
    461. 

  461. #process.p = cms.Path( process.gsfElectrons*process.ourJetSequence*process.patDefaultSequence*process.zeeFilter*process.plotter)
    462. 

  462. 

  463. # For DATA, 397 electrons are fixed (misalignment corrections included)
    464. 

  464. process.p = cms.Path( process.ourJetSequence*process.patDefaultSequence*process.zeeFilter*process.plotter)
    465. 

  465. 

  466. # For MonteCarlo,
    467. 

  467. #process.p = cms.Path( process.ourJetSequence * process.patDefaultSequence +process.zeeFilter + process.plotter)
    468.