ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/test/pphoton.py
Revision: 1.1
Committed: Fri Dec 2 22:04:54 2011 UTC (13 years, 5 months ago) by yilmaz
Content type: text/x-python
Branch: MAIN
CVS Tags: HiForest_V02_79, HiForest_V02_78, HiForest_V02_77, HiForest_V02_76, HiForest_V02_75, HiForest_V02_74, HiForest_V02_73, HiForest_V02_72, HiForest_V02_71, HiForest_V02_70, HiForest_V02_69, HiForest_V02_68, HiForest_V02_67, HiForest_V02_66, HiForest_V02_65, HiForest_V02_64, HiForest_V02_63, HiForest_V02_62, HiForest_V02_61, HiForest_V02_60, HiForest_V02_59, HiForest_V02_58, HiForest_V02_57, HiForest_V02_56, HiForest_V02_55, HiForest_V02_54, HiForest_V02_53, HiForest_V02_52, HiForest_V02_51, HiForest_V02_50, HiForest_V02_49, HiForest_V02_48, HiForest_V02_47, HiForest_V02_46, HiForest_V02_45, HiForest_V02_44, HiForest_V02_43, HiForest_V02_42, HiForest_V02_41, HiForest_V02_40, HiForest_V02_39, HiForest_V02_38, HiForest_V02_37, HiForest_V02_36, HiForest_V02_35, HiForest_V02_34, HiForest_V02_33, HiForest_V02_32, HiForest_V02_31, HiForest_V02_30, HiForest_V02_27, HiForest_V02_26, QM_2012, HiForest_V02_25, HiForest_V02_24, HiForest_V02_23, HiForest_V02_22, HiForest_V02_21, HiForest_V02_20, HiForest_V02_19, HiForest_V02_18, HiForest_V02_17, HiForest_V02_16, HiForest_V02_15, HiForest_V02_14, HiForest_V02_13, HiForest_V02_12, HiForest_V02_11, HiForest_V02_10, HiForest_V02_09, HiForest_V02_08, HiForest_V02_07, HiForest_V02_06, HiForest_V02_05, HiForest_V02_04, HiForest_V02_03, HiForest_V02_02, HiForest_V02_01, HiForest_V02_00
Branch point for: branch_44x
Log Message:
photon skim + forest

File Contents

# User Rev Content
1 yilmaz 1.1 import FWCore.ParameterSet.Config as cms
2    
3     process = cms.Process('hiForestAna2011')
4    
5     process.options = cms.untracked.PSet(
6     wantSummary = cms.untracked.bool(True)
7     )
8    
9     #####################################################################################
10     # Input source
11     #####################################################################################
12    
13     process.source = cms.Source("PoolSource",
14     duplicateCheckMode = cms.untracked.string("noDuplicateCheck"),
15     fileNames = cms.untracked.vstring(
16     '/store/hidata/HIRun2011/HIHighPt/RECO/PromptReco-v1/000/181/611/7C1B26DC-0E10-E111-9CEF-003048F117F6.root'
17     # 'file:rereco/EED77885-BAF8-E011-B6B7-001D09F2532F.root'
18     ))
19    
20     # Number of events we want to process, -1 = all events
21     process.maxEvents = cms.untracked.PSet(
22     input = cms.untracked.int32(-1))
23    
24    
25     #####################################################################################
26     # Load some general stuff
27     #####################################################################################
28    
29     process.load("HLTrigger.HLTfilters.hltHighLevel_cfi")
30     process.photonFilter = process.hltHighLevel.clone()
31     process.photonFilter.HLTPaths = ["HLT_HISinglePhoton30_v*"]
32    
33    
34     process.singlePhotonPt35Filter = cms.EDFilter("PhotonSelector",
35     src = cms.InputTag("photons"),
36     cut = cms.string('pt > 35 && abs(eta) < 1.48 && sigmaIetaIeta > 0.002' ),
37     filter = cms.bool(True)
38     )
39    
40     process.superFilterSequence = cms.Sequence(process.photonFilter*process.singlePhotonPt35Filter)
41     process.superFilterPath = cms.Path(process.superFilterSequence)
42    
43     process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
44     process.load('Configuration.StandardSequences.Services_cff')
45     process.load('Configuration.StandardSequences.GeometryExtended_cff')
46     process.load('Configuration.StandardSequences.MagneticField_38T_cff')
47     process.load('Configuration.StandardSequences.ReconstructionHeavyIons_cff')
48     process.load('FWCore.MessageService.MessageLogger_cfi')
49     process.load('RecoLocalTracker.SiPixelRecHits.PixelCPEESProducers_cff')
50     # Data Global Tag 44x
51     process.GlobalTag.globaltag = 'GR_P_V27A::All'
52    
53     # MC Global Tag 44x
54     #process.GlobalTag.globaltag = 'STARTHI44_V4::All'
55    
56     # load centrality
57     from CmsHi.Analysis2010.CommonFunctions_cff import *
58     overrideCentrality(process)
59     process.HeavyIonGlobalParameters = cms.PSet(
60     centralityVariable = cms.string("HFtowers"),
61     nonDefaultGlauberModel = cms.string(""),
62     centralitySrc = cms.InputTag("hiCentrality")
63     )
64    
65     process.load('RecoHI.HiCentralityAlgos.HiCentrality_cfi')
66    
67     # EcalSeverityLevel ES Producer
68     process.load("RecoLocalCalo/EcalRecAlgos/EcalSeverityLevelESProducer_cfi")
69     process.load("RecoEcal.EgammaCoreTools.EcalNextToDeadChannelESProducer_cff")
70    
71    
72     #####################################################################################
73     # Define tree output
74     #####################################################################################
75    
76     process.TFileService = cms.Service("TFileService",
77     fileName=cms.string("HiForest.root"))
78    
79     #####################################################################################
80     # Jet energy correction
81     #####################################################################################
82    
83     #####################################################################################
84     # Additional Reconstruction and Analysis
85     #####################################################################################
86    
87     # MET: Calorimeter based MET
88     process.load("RecoMET.METProducers.CaloMET_cfi")
89    
90     # Define Analysis sequencues
91     process.load('CmsHi.JetAnalysis.EventSelection_cff')
92     #process.load('CmsHi.JetAnalysis.ExtraGenReco_cff')
93     #process.load('CmsHi.JetAnalysis.ExtraTrackReco_cff')
94     process.load('CmsHi.JetAnalysis.ExtraPfReco_cff')
95     process.load('CmsHi.JetAnalysis.ExtraJetReco_cff')
96     process.load('CmsHi.JetAnalysis.ExtraEGammaReco_cff')
97     process.load('CmsHi.JetAnalysis.PatAna_cff')
98     process.load('CmsHi.JetAnalysis.JetAnalyzers_cff')
99     process.load('CmsHi.JetAnalysis.EGammaAnalyzers_cff')
100     process.load("MitHig.PixelTrackletAnalyzer.trackAnalyzer_cff")
101     process.anaTrack.trackPtMin = 0
102     process.anaTrack.useQuality = True
103     process.anaTrack.doPFMatching = True
104     process.anaTrack.trackSrc = cms.InputTag("hiSelectedTracks")
105     process.load("MitHig.PixelTrackletAnalyzer.METAnalyzer_cff")
106     process.load("CmsHi.JetAnalysis.pfcandAnalyzer_cfi")
107     process.pfcandAnalyzer.skipCharged = False
108     process.pfcandAnalyzer.pfPtMin = 0
109     process.interestingTrackEcalDetIds.TrackCollection = cms.InputTag("hiSelectedTracks")
110     process.akPu3CaloJetAnalyzer.isMC = False
111     process.akPu5CaloJetAnalyzer.isMC = False
112     process.akPu3PFJetAnalyzer.isMC = False
113     process.akPu5PFJetAnalyzer.isMC = False
114     # Muons
115     process.load("MuTrig.HLTMuTree.hltMuTree_cfi")
116     process.muonTree = process.hltMuTree.clone()
117     process.muonTree.doGen = cms.untracked.bool(False)
118    
119     # Event tree
120     process.load("CmsHi/HiHLTAlgos.hievtanalyzer_cfi")
121     # Not working for the moment..
122     process.hiEvtAnalyzer.doMC = cms.bool(False)
123     process.hiEvtAnalyzer.doEvtPlane = cms.bool(True)
124    
125     process.ak5CaloJets = process.akPu5CaloJets.clone(doPUOffsetCorr = False)
126     process.ak5corr = process.icPu5corr.clone(
127     src = cms.InputTag("ak5CaloJets"),
128     payload = cms.string('AK5Calo')
129     )
130     process.ak5patJets = process.akPu5PFpatJets.clone(
131     jetSource = cms.InputTag("ak5CaloJets"),
132     jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak5corr"))
133     )
134    
135     #process.icPu5JetAnalyzer.hltTrgResults = cms.untracked.string('TriggerResults::RECO')
136     #process.akPu3PFJetAnalyzer.hltTrgResults = cms.untracked.string('TriggerResults::RECO')
137     process.icPu5JetAnalyzer.useCentrality = cms.untracked.bool(False) # doesn't fill cent info
138     process.akPu3PFJetAnalyzer.useCentrality = cms.untracked.bool(False) # doesn't fill cent info
139    
140     process.ak5CaloJetAnalyzer = process.icPu5JetAnalyzer.clone(
141     jetTag = 'ak5patJets',
142     genjetTag = 'ak5HiGenJets',
143     isMC = False
144     )
145    
146     process.ak3PFJetsX = process.akPu3PFJets.clone(doPUOffsetCorr = False)
147     process.ak3corrX = process.icPu5corr.clone(
148     src = cms.InputTag("ak3PFJetsX"),
149     payload = cms.string('AK3PF')
150     )
151     process.ak3patJetsX = process.akPu5PFpatJets.clone(
152     jetSource = cms.InputTag("ak3PFJetsX"),
153     jetCorrFactorsSource = cms.VInputTag(cms.InputTag("ak3corrX"))
154     )
155     process.ak3PFJetAnalyzer = process.icPu5JetAnalyzer.clone(
156     jetTag = 'ak3patJetsX',
157     genjetTag = 'ak3HiGenJets',
158     isMC = False
159     )
160    
161     process.ak5extra = cms.Sequence(process.ak5CaloJets*process.ak5corr*process.ak5patJets*process.ak5CaloJetAnalyzer)
162     process.ak3extra = cms.Sequence(process.ak3PFJetsX*process.ak3corrX*process.ak3patJetsX*process.ak3PFJetAnalyzer)
163    
164     #process.load("edwenger.Skims.EventFilter_cff")
165     #from edwenger.Skims.customise_cfi import *
166     #run2760GeVmode(process)
167    
168     # Filtering
169     process.hltJetHI.HLTPaths = ['HLT_HIJet35U','HLT_HIPhoton20']
170     #process.hltJetHI.TriggerResultsTag = cms.InputTag("TriggerResults::RECO")
171    
172     print "Add cleaning to analysis"
173     process.event_filter_seq = cms.Sequence(
174     process.hltJetHI *
175     process.siPixelRecHits *
176     process.collisionEventSelection *
177     process.HBHENoiseFilter *
178     process.hiEcalRecHitSpikeFilter
179    
180     # process.preTrgTest *
181     # process.minBiasBscFilter *
182     # process.postTrgTest *
183     # process.hfCoincFilter *
184     # process.purityFractionFilter
185    
186     )
187    
188     #Commented by Yen-Jie
189     #process.hiPixelAdaptiveVertex.useBeamConstraint = False
190    
191     process.load("RecoHI.HiMuonAlgos.HiRecoMuon_cff")
192     process.muons.JetExtractorPSet.JetCollectionLabel = cms.InputTag("iterativeConePu5CaloJets")
193    
194     process.load("UserCode.L1TriggerDPG.l1NtupleProducer_cfi")
195    
196     #process.hiGoodTightTracks.src = cms.InputTag("hiGlobalPrimTracks")
197     #process.hiGoodTightTracksDirect = process.hiGoodTightTracks.clone(keepAllTracks = True)
198     #process.hiGoodTracks = process.hiGoodTightTracks.clone()
199     process.akPu5PFJets.doAreaFastjet = False
200     process.akPu3PFJets.doAreaFastjet = False
201     process.iterativeConePu5CaloJets.doAreaFastjet = False
202     process.akPu5PFJets.doRhoFastjet = False
203     process.akPu3PFJets.doRhoFastjet = False
204     # process.iterativeCONEPu5CaloJets.doRhoFastjet = False
205    
206     #process.hiSelectedTrackHighPurity = cms.EDProducer("QualityFilter",
207     # TrackQuality = cms.string('highPurity'),
208     # recTracks = cms.InputTag("hiSelectedTracks")
209     #)
210    
211     process.hiSelectedTrackHighPurity = cms.EDFilter("TrackSelector",
212     src = cms.InputTag("hiSelectedTracks"),
213     cut = cms.string(
214     'quality("highPurity")')
215     )
216    
217     process.particleFlowClusterPS.thresh_Pt_Seed_Endcap = cms.double(99999.)
218     process.reco_extra = cms.Path( process.siPixelRecHits * process.siStripMatchedRecHits *
219     process.hiPixel3PrimTracks *
220     process.hiPixelTrackSeeds *
221     process.hiSelectedTrackHighPurity *
222     process.electronGsfTrackingHi *
223     process.hiElectronSequence *
224     process.HiParticleFlowReco *
225     process.iterativeConePu5CaloJets *
226     process.PFTowers
227     )
228    
229     process.reco_extra_jet = cms.Path( process.iterativeConePu5CaloJets *
230     process.akPu3PFJets * process.akPu5PFJets *
231     process.akPu3CaloJets * process.akPu5CaloJets *
232     process.photon_extra_reco)
233     #process.gen_step = cms.Path( process.hiGenParticles * process.hiGenParticlesForJets * process.genPartons * process.hiPartons * process.hiRecoGenJets)
234     process.pat_step = cms.Path(process.icPu5patSequence_data +
235     process.akPu3PFpatSequence_data + process.akPu5PFpatSequence_data +
236     process.akPu3patSequence_data + process.akPu5patSequence_data +
237     process.makeHeavyIonPhotons)
238     process.pat_step.remove(process.interestingTrackEcalDetIds)
239     process.pat_step.remove(process.photonMatch)
240     #+ process.patPhotons)
241    
242     process.patPhotons.addPhotonID = cms.bool(False)
243     #process.makeHeavyIonPhotons)
244     process.extrapatstep = cms.Path(process.selectedPatPhotons)
245     process.load("CmsHi.JetAnalysis.hcalNoise_cff")
246     process.load("CmsHi.JetAnalysis.RandomCones_cff")
247     process.akPu3PFConesAna.doMC = cms.untracked.bool(False)
248     process.akPu5PFConesAna.doMC = cms.untracked.bool(False)
249     process.akPu3CaloConesAna.doMC = cms.untracked.bool(False)
250     process.akPu5CaloConesAna.doMC = cms.untracked.bool(False)
251     process.icPu5CaloConesAna.doMC = cms.untracked.bool(False)
252    
253     process.multiPhotonAnalyzer.GammaEtaMax = cms.untracked.double(100)
254     process.multiPhotonAnalyzer.GammaPtMin = cms.untracked.double(0)
255     process.multiPhotonAnalyzer.gsfElectronCollection = cms.untracked.InputTag("ecalDrivenGsfElectrons")
256     process.ana_step = cms.Path( process.jetana_seq +
257     process.multiPhotonAnalyzer + process.anaTrack + process.pfcandAnalyzer +
258     process.met * process.anaMET +
259     process.muonTree +
260     process.hcalNoise +
261     process.hiEvtAnalyzer +
262     process.randomCones
263     )
264    
265     process.phltJetHI = cms.Path( process.hltJetHI )
266     process.pcollisionEventSelection = cms.Path(process.collisionEventSelection)
267     process.pHBHENoiseFilter = cms.Path( process.HBHENoiseFilter )
268     process.phiEcalRecHitSpikeFilter = cms.Path(process.hiEcalRecHitSpikeFilter )
269     #process.ppreTrgTest = cms.Path(process.preTrgTest )
270     #process.pminBiasBscFilter = cms.Path(process.minBiasBscFilter )
271     #process.ppostTrgTest = cms.Path(process.postTrgTest )
272     #process.phfCoincFilter = cms.Path(process.hfCoincFilter )
273     #process.ppurityFractionFilter = cms.Path(process.purityFractionFilter )
274    
275     # Customization
276     from CmsHi.JetAnalysis.customise_cfi import *
277     enableDataPat(process)
278     setPhotonObject(process,"cleanPhotons")
279     enableDataAnalyzers(process)
280    
281     # process.load('HLTrigger.HLTanalyzers.HLTBitAnalyser_cfi')
282     # process.hltbitanalysis.UseTFileService = cms.untracked.bool(True)
283     # process.hltanalysis = process.hltbitanalysis.clone(
284     # l1GtReadoutRecord = cms.InputTag("gtDigis"),
285     # l1GctHFBitCounts = cms.InputTag("gctDigis"),
286     # l1GctHFRingSums = cms.InputTag("gctDigis"),
287     # l1extramu = cms.string('l1extraParticles'),
288     # l1extramc = cms.string('l1extraParticles'),
289     # hltresults = cms.InputTag("TriggerResults","","HLT"),
290     # HLTProcessName = cms.string("HLT")
291     # )
292    
293     process.load('HLTrigger.HLTanalyzers.HLTBitAnalyser_cfi')
294     process.hltbitanalysis.UseTFileService = cms.untracked.bool(True)
295     process.hltanalysis = process.hltbitanalysis.clone(
296     dummyBranches = cms.untracked.vstring(),
297     l1GtReadoutRecord = cms.InputTag("gtDigis"),
298     l1GctHFBitCounts = cms.InputTag("gctDigis"),
299     l1GctHFRingSums = cms.InputTag("gctDigis"),
300     l1extramu = cms.string('l1extraParticles'),
301     l1extramc = cms.string('l1extraParticles'),
302     hltresults = cms.InputTag("TriggerResults","","HLT"),
303     HLTProcessName = cms.string("HLT")
304     )
305     process.hltanalysis.dummyBranches.extend( [
306     "HLT_HIL1DoubleMu0_HighQ_v1",
307     "HLT_HIL1DoubleMuOpen_v1",
308     "HLT_HIL2DoubleMu0_L1HighQL2NHitQ_v1",
309     "HLT_HIL2DoubleMu0_NHitQ_v1",
310     "HLT_HIL2DoubleMu0_v1",
311     "HLT_HIL2DoubleMu3_v1",
312     "HLT_HIL2Mu15_v1",
313     "HLT_HIL2Mu3_NHitQ_v1",
314     "HLT_HIL2Mu3_v1",
315     "HLT_HIL2Mu7_v1",
316     "HLT_HIL3DoubleMuOpen_Mgt2_OS_NoCowboy_v1",
317     "HLT_HIL3DoubleMuOpen_Mgt2_OS_v1",
318     "HLT_HIL3DoubleMuOpen_Mgt2_SS_v1",
319     "HLT_HIL3DoubleMuOpen_Mgt2_v1",
320     "HLT_HIL3DoubleMuOpen_v1",
321     "HLT_HIL3Mu3_v1",
322     "HLT_HIDiJet55_v1",
323     "HLT_HIDoublePhoton10_v1",
324     "HLT_HIDoublePhoton15_v1",
325     "HLT_HIDoublePhoton20_v1",
326     "HLT_HIFullTrack12_L1Central_v1",
327     "HLT_HIFullTrack12_L1Peripheral_v1",
328     "HLT_HIFullTrack14_L1Central_v1",
329     "HLT_HIFullTrack14_L1Peripheral_v1",
330     "HLT_HIFullTrack20_L1Central_v1",
331     "HLT_HIFullTrack20_L1Peripheral_v1",
332     "HLT_HIFullTrack25_L1Central_v1",
333     "HLT_HIFullTrack25_L1Peripheral_v1",
334     "HLT_HIJet55_v1",
335     "HLT_HIJet65_Jet55_v1",
336     "HLT_HIJet65_v1",
337     "HLT_HIJet80_v1",
338     "HLT_HIJet95_v1",
339     "HLT_HIJetE30_NoBPTX_v1",
340     "HLT_HIJetE50_NoBPTX3BX_NoHalo_v1",
341     "HLT_HIMET120_v1",
342     "HLT_HIMET200_v1",
343     "HLT_HIMET220_v1",
344     "HLT_HIPhoton10_Photon15_v1",
345     "HLT_HIPhoton15_Photon20_v1",
346     "HLT_HISinglePhoton15_v1",
347     "HLT_HISinglePhoton20_v1",
348     "HLT_HISinglePhoton30_v1",
349     "HLT_HISinglePhoton40_v1",
350     "HLT_HIBptxXOR_v1",
351     "HLT_HICentral10_v1",
352     "HLT_HICentralityVeto_v1",
353     "HLT_HIL1Algo_BptxXOR_BSC_OR_v1",
354     "HLT_HIMinBiasBSC_OR_v1",
355     "HLT_HIMinBiasBSC_v1",
356     "HLT_HIMinBiasHF_v1",
357     "HLT_HIMinBiasHfOrBSC_v1",
358     "HLT_HIMinBiasHf_OR_v1",
359     "HLT_HIMinBiasPixel_SingleTrack_v1",
360     "HLT_HIMinBiasZDCPixel_SingleTrack_v1",
361     "HLT_HIMinBiasZDC_Calo_PlusOrMinus_v1",
362     "HLT_HIMinBiasZDC_Calo_v1",
363     "HLT_HIMinBiasZDC_PlusOrMinusPixel_SingleTrack_v1",
364     "HLT_HIPhysics_v1",
365     "HLT_HIRandom_v1",
366     "HLT_HIUCC010_v1",
367     "HLT_HIUCC015_v1",
368     "HLT_HIUPCNeuEG2Pixel_SingleTrack_v1",
369     "HLT_HIUPCNeuEG5Pixel_SingleTrack_v1",
370     "HLT_HIUPCNeuHcalHfEG2Pixel_SingleTrack_v1",
371     "HLT_HIUPCNeuHcalHfEG5Pixel_SingleTrack_v1",
372     "HLT_HIUPCNeuHcalHfMuPixel_SingleTrack_v1",
373     "HLT_HIUPCNeuMuPixel_SingleTrack_v1",
374     "HLT_HIZeroBiasPixel_SingleTrack_v1",
375     "HLT_HIZeroBiasXOR_v1",
376     "HLT_HIZeroBias_v1",
377     ] )
378     process.hltanalysis.dummyBranches.extend( [
379     "HLT_HIL1DoubleMu0_HighQ_v2",
380     "HLT_HIL1DoubleMuOpen_v2",
381     "HLT_HIL2DoubleMu0_L1HighQL2NHitQ_v2",
382     "HLT_HIL2DoubleMu0_NHitQ_v2",
383     "HLT_HIL2DoubleMu0_v2",
384     "HLT_HIL2DoubleMu3_v2",
385     "HLT_HIL2Mu15_v2",
386     "HLT_HIL2Mu3_NHitQ_v2",
387     "HLT_HIL2Mu3_v2",
388     "HLT_HIL2Mu7_v2",
389     "HLT_HIL3DoubleMuOpen_Mgt2_OS_NoCowboy_v2",
390     "HLT_HIL3DoubleMuOpen_Mgt2_OS_v2",
391     "HLT_HIL3DoubleMuOpen_Mgt2_SS_v2",
392     "HLT_HIL3DoubleMuOpen_Mgt2_v2",
393     "HLT_HIL3DoubleMuOpen_v2",
394     "HLT_HIL3Mu3_v2",
395     "HLT_HIDiJet55_v2",
396     "HLT_HIDoublePhoton10_v2",
397     "HLT_HIDoublePhoton15_v2",
398     "HLT_HIDoublePhoton20_v2",
399     "HLT_HIFullTrack12_L1Central_v2",
400     "HLT_HIFullTrack12_L1Peripheral_v2",
401     "HLT_HIFullTrack14_L1Central_v2",
402     "HLT_HIFullTrack14_L1Peripheral_v2",
403     "HLT_HIFullTrack20_L1Central_v2",
404     "HLT_HIFullTrack20_L1Peripheral_v2",
405     "HLT_HIFullTrack25_L1Central_v2",
406     "HLT_HIFullTrack25_L1Peripheral_v2",
407     "HLT_HIJet55_v2",
408     "HLT_HIJet65_Jet55_v2",
409     "HLT_HIJet65_v2",
410     "HLT_HIJet80_v2",
411     "HLT_HIJet95_v2",
412     "HLT_HIJetE30_NoBPTX_v2",
413     "HLT_HIJetE50_NoBPTX3BX_NoHalo_v2",
414     "HLT_HIMET120_v2",
415     "HLT_HIMET200_v2",
416     "HLT_HIMET220_v2",
417     "HLT_HIPhoton10_Photon15_v2",
418     "HLT_HIPhoton15_Photon20_v2",
419     "HLT_HISinglePhoton15_v2",
420     "HLT_HISinglePhoton20_v2",
421     "HLT_HISinglePhoton30_v2",
422     "HLT_HISinglePhoton40_v2",
423     "HLT_HIBptxXOR_v2",
424     "HLT_HICentral10_v2",
425     "HLT_HICentralityVeto_v2",
426     "HLT_HIL1Algo_BptxXOR_BSC_OR_v2",
427     "HLT_HIMinBiasBSC_OR_v2",
428     "HLT_HIMinBiasBSC_v2",
429     "HLT_HIMinBiasHF_v2",
430     "HLT_HIMinBiasHfOrBSC_v2",
431     "HLT_HIMinBiasHf_OR_v2",
432     "HLT_HIMinBiasPixel_SingleTrack_v2",
433     "HLT_HIMinBiasZDCPixel_SingleTrack_v2",
434     "HLT_HIMinBiasZDC_Calo_PlusOrMinus_v2",
435     "HLT_HIMinBiasZDC_Calo_v2",
436     "HLT_HIMinBiasZDC_PlusOrMinusPixel_SingleTrack_v2",
437     "HLT_HIPhysics_v2",
438     "HLT_HIRandom_v2",
439     "HLT_HIUCC010_v2",
440     "HLT_HIUCC015_v2",
441     "HLT_HIUPCNeuEG2Pixel_SingleTrack_v2",
442     "HLT_HIUPCNeuEG5Pixel_SingleTrack_v2",
443     "HLT_HIUPCNeuHcalHfEG2Pixel_SingleTrack_v2",
444     "HLT_HIUPCNeuHcalHfEG5Pixel_SingleTrack_v2",
445     "HLT_HIUPCNeuHcalHfMuPixel_SingleTrack_v2",
446     "HLT_HIUPCNeuMuPixel_SingleTrack_v2",
447     "HLT_HIZeroBiasPixel_SingleTrack_v2",
448     "HLT_HIZeroBiasXOR_v2",
449     "HLT_HIZeroBias_v2",
450     ] )
451     process.hltanalysis.dummyBranches.extend( [
452     "HLT_HIL1DoubleMu0_HighQ_v3",
453     "HLT_HIL1DoubleMuOpen_v3",
454     "HLT_HIL2DoubleMu0_L1HighQL2NHitQ_v3",
455     "HLT_HIL2DoubleMu0_NHitQ_v3",
456     "HLT_HIL2DoubleMu0_v3",
457     "HLT_HIL2DoubleMu3_v3",
458     "HLT_HIL2Mu15_v3",
459     "HLT_HIL2Mu3_NHitQ_v3",
460     "HLT_HIL2Mu3_v3",
461     "HLT_HIL2Mu7_v3",
462     "HLT_HIL3DoubleMuOpen_Mgt2_OS_NoCowboy_v3",
463     "HLT_HIL3DoubleMuOpen_Mgt2_OS_v3",
464     "HLT_HIL3DoubleMuOpen_Mgt2_SS_v3",
465     "HLT_HIL3DoubleMuOpen_Mgt2_v3",
466     "HLT_HIL3DoubleMuOpen_v3",
467     "HLT_HIL3Mu3_v3",
468     "HLT_HIDiJet55_v3",
469     "HLT_HIDoublePhoton10_v3",
470     "HLT_HIDoublePhoton15_v3",
471     "HLT_HIDoublePhoton20_v3",
472     "HLT_HIFullTrack12_L1Central_v3",
473     "HLT_HIFullTrack12_L1Peripheral_v3",
474     "HLT_HIFullTrack14_L1Central_v3",
475     "HLT_HIFullTrack14_L1Peripheral_v3",
476     "HLT_HIFullTrack20_L1Central_v3",
477     "HLT_HIFullTrack20_L1Peripheral_v3",
478     "HLT_HIFullTrack25_L1Central_v3",
479     "HLT_HIFullTrack25_L1Peripheral_v3",
480     "HLT_HIJet55_v3",
481     "HLT_HIJet65_Jet55_v3",
482     "HLT_HIJet65_v3",
483     "HLT_HIJet80_v3",
484     "HLT_HIJet95_v3",
485     "HLT_HIJetE30_NoBPTX_v3",
486     "HLT_HIJetE50_NoBPTX3BX_NoHalo_v3",
487     "HLT_HIMET120_v3",
488     "HLT_HIMET200_v3",
489     "HLT_HIMET220_v3",
490     "HLT_HIPhoton10_Photon15_v3",
491     "HLT_HIPhoton15_Photon20_v3",
492     "HLT_HISinglePhoton15_v3",
493     "HLT_HISinglePhoton20_v3",
494     "HLT_HISinglePhoton30_v3",
495     "HLT_HISinglePhoton40_v3",
496     "HLT_HIBptxXOR_v3",
497     "HLT_HICentral10_v3",
498     "HLT_HICentralityVeto_v3",
499     "HLT_HIL1Algo_BptxXOR_BSC_OR_v3",
500     "HLT_HIMinBiasBSC_OR_v3",
501     "HLT_HIMinBiasBSC_v3",
502     "HLT_HIMinBiasHF_v3",
503     "HLT_HIMinBiasHfOrBSC_v3",
504     "HLT_HIMinBiasHf_OR_v3",
505     "HLT_HIMinBiasPixel_SingleTrack_v3",
506     "HLT_HIMinBiasZDCPixel_SingleTrack_v3",
507     "HLT_HIMinBiasZDC_Calo_PlusOrMinus_v3",
508     "HLT_HIMinBiasZDC_Calo_v3",
509     "HLT_HIMinBiasZDC_PlusOrMinusPixel_SingleTrack_v3",
510     "HLT_HIPhysics_v3",
511     "HLT_HIRandom_v3",
512     "HLT_HIUCC010_v3",
513     "HLT_HIUCC015_v3",
514     "HLT_HIUPCNeuEG2Pixel_SingleTrack_v3",
515     "HLT_HIUPCNeuEG5Pixel_SingleTrack_v3",
516     "HLT_HIUPCNeuHcalHfEG2Pixel_SingleTrack_v3",
517     "HLT_HIUPCNeuHcalHfEG5Pixel_SingleTrack_v3",
518     "HLT_HIUPCNeuHcalHfMuPixel_SingleTrack_v3",
519     "HLT_HIUPCNeuMuPixel_SingleTrack_v3",
520     "HLT_HIZeroBiasPixel_SingleTrack_v3",
521     "HLT_HIZeroBiasXOR_v3",
522     "HLT_HIZeroBias_v3",
523     ] )
524     process.hltanalysis.dummyBranches.extend( [
525     "HLT_HIL1DoubleMu0_HighQ_v4",
526     "HLT_HIL1DoubleMuOpen_v4",
527     "HLT_HIL2DoubleMu0_L1HighQL2NHitQ_v4",
528     "HLT_HIL2DoubleMu0_NHitQ_v4",
529     "HLT_HIL2DoubleMu0_v4",
530     "HLT_HIL2DoubleMu3_v4",
531     "HLT_HIL2Mu15_v4",
532     "HLT_HIL2Mu3_NHitQ_v4",
533     "HLT_HIL2Mu3_v4",
534     "HLT_HIL2Mu7_v4",
535     "HLT_HIL3DoubleMuOpen_Mgt2_OS_NoCowboy_v4",
536     "HLT_HIL3DoubleMuOpen_Mgt2_OS_v4",
537     "HLT_HIL3DoubleMuOpen_Mgt2_SS_v4",
538     "HLT_HIL3DoubleMuOpen_Mgt2_v4",
539     "HLT_HIL3DoubleMuOpen_v4",
540     "HLT_HIL3Mu3_v4",
541     "HLT_HIDiJet55_v4",
542     "HLT_HIDoublePhoton10_v4",
543     "HLT_HIDoublePhoton15_v4",
544     "HLT_HIDoublePhoton20_v4",
545     "HLT_HIFullTrack12_L1Central_v4",
546     "HLT_HIFullTrack12_L1Peripheral_v4",
547     "HLT_HIFullTrack14_L1Central_v4",
548     "HLT_HIFullTrack14_L1Peripheral_v4",
549     "HLT_HIFullTrack20_L1Central_v4",
550     "HLT_HIFullTrack20_L1Peripheral_v4",
551     "HLT_HIFullTrack25_L1Central_v4",
552     "HLT_HIFullTrack25_L1Peripheral_v4",
553     "HLT_HIJet55_v4",
554     "HLT_HIJet65_Jet55_v4",
555     "HLT_HIJet65_v4",
556     "HLT_HIJet80_v4",
557     "HLT_HIJet95_v4",
558     "HLT_HIJetE30_NoBPTX_v4",
559     "HLT_HIJetE50_NoBPTX3BX_NoHalo_v4",
560     "HLT_HIMET120_v4",
561     "HLT_HIMET200_v4",
562     "HLT_HIMET220_v4",
563     "HLT_HIPhoton10_Photon15_v4",
564     "HLT_HIPhoton15_Photon20_v4",
565     "HLT_HISinglePhoton15_v4",
566     "HLT_HISinglePhoton20_v4",
567     "HLT_HISinglePhoton30_v4",
568     "HLT_HISinglePhoton40_v4",
569     "HLT_HIBptxXOR_v4",
570     "HLT_HICentral10_v4",
571     "HLT_HICentralityVeto_v4",
572     "HLT_HIL1Algo_BptxXOR_BSC_OR_v4",
573     "HLT_HIMinBiasBSC_OR_v4",
574     "HLT_HIMinBiasBSC_v4",
575     "HLT_HIMinBiasHF_v4",
576     "HLT_HIMinBiasHfOrBSC_v4",
577     "HLT_HIMinBiasHf_OR_v4",
578     "HLT_HIMinBiasPixel_SingleTrack_v4",
579     "HLT_HIMinBiasZDCPixel_SingleTrack_v4",
580     "HLT_HIMinBiasZDC_Calo_PlusOrMinus_v4",
581     "HLT_HIMinBiasZDC_Calo_v4",
582     "HLT_HIMinBiasZDC_PlusOrMinusPixel_SingleTrack_v4",
583     "HLT_HIPhysics_v4",
584     "HLT_HIRandom_v4",
585     "HLT_HIUCC010_v4",
586     "HLT_HIUCC015_v4",
587     "HLT_HIUPCNeuEG2Pixel_SingleTrack_v4",
588     "HLT_HIUPCNeuEG5Pixel_SingleTrack_v4",
589     "HLT_HIUPCNeuHcalHfEG2Pixel_SingleTrack_v4",
590     "HLT_HIUPCNeuHcalHfEG5Pixel_SingleTrack_v4",
591     "HLT_HIUPCNeuHcalHfMuPixel_SingleTrack_v4",
592     "HLT_HIUPCNeuMuPixel_SingleTrack_v4",
593     "HLT_HIZeroBiasPixel_SingleTrack_v4",
594     "HLT_HIZeroBiasXOR_v4",
595     "HLT_HIZeroBias_v4",
596     ] )
597    
598     #process.hltanalysis.hltresults = cms.InputTag("TriggerResults","","RECO")
599     process.skimanalysis = cms.EDAnalyzer("FilterAnalyzer",
600     hltresults = cms.InputTag("TriggerResults","","hiForestAna2011")
601     )
602     process.hltAna = cms.Path(process.hltanalysis)
603     process.pAna = cms.EndPath(process.skimanalysis)
604    
605    
606    
607     process.load('CmsHi.JetAnalysis.rechitanalyzer_cfi')
608    
609     process.rechitanalyzer.HBHETreePtMin = cms.untracked.double(0.5)
610     process.rechitanalyzer.HFTreePtMin = cms.untracked.double(0.5)
611     process.rechitanalyzer.EBTreePtMin = cms.untracked.double(0.5)
612     process.rechitanalyzer.EETreePtMin = cms.untracked.double(0.5)
613     process.rechitanalyzer.TowerTreePtMin = cms.untracked.double(0.5)
614     process.rechitanalyzer.doHF = cms.untracked.bool(True)
615     process.rechitAna = cms.Path(process.rechitanalyzer)
616    
617     ########### random number seed
618     process.RandomNumberGeneratorService.akPu3PFConesAna = process.RandomNumberGeneratorService.generator.clone()
619     process.RandomNumberGeneratorService.icPu5CaloConesAna = process.RandomNumberGeneratorService.generator.clone()
620     process.RandomNumberGeneratorService.akPu5PFConesAna = process.RandomNumberGeneratorService.generator.clone()
621     process.RandomNumberGeneratorService.akPu3CaloConesAna = process.RandomNumberGeneratorService.generator.clone()
622     process.RandomNumberGeneratorService.akPu5CaloConesAna = process.RandomNumberGeneratorService.generator.clone()
623     process.RandomNumberGeneratorService.multiPhotonAnalyzer = process.RandomNumberGeneratorService.generator.clone()
624    
625     #####################################################################################
626     # Edm Output
627     #####################################################################################
628    
629     #process.out = cms.OutputModule("PoolOutputModule",
630     # fileName = cms.untracked.string("output.root")
631     # )
632     #process.save = cms.EndPath(process.out)
633    
634    
635     process.skimanalysis.superFilters = cms.vstring("superFilterPath")
636     for path in process.paths:
637     getattr(process,path)._seq = process.superFilterSequence*getattr(process,path)._seq
638    
639    
640    
641