CMS 3D CMS Logo

/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/Validation/RecoHI/python/JetValidationHeavyIons_cff.py

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 ### genjet cleaning for improved matching in HI environment
00004 
00005 from RecoHI.HiJetAlgos.HiGenCleaner_cff import *
00006 
00007 iterativeCone5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets'))
00008 iterativeCone7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone7HiGenJets'))
00009 ak5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets'))
00010 ak7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak7HiGenJets'))
00011 
00012 ### jet analyzer for two radii (0.5, 0.7) and three algorithms:
00013 ### iterative cone with PU, anti-kt with PU, anti-kt with fastjet PU
00014 
00015 JetAnalyzerICPU5Calo = cms.EDAnalyzer("CaloJetTester",
00016     src = cms.InputTag("iterativeConePu5CaloJets"),
00017     srcGen = cms.InputTag("iterativeCone5HiCleanedGenJets"),
00018     genEnergyFractionThreshold = cms.double(0.05),
00019     genPtThreshold = cms.double(1.0),
00020     RThreshold = cms.double(0.3),
00021     reverseEnergyFractionThreshold = cms.double(0.5)                                    
00022 )
00023 
00024 JetAnalyzerICPU7Calo = cms.EDAnalyzer("CaloJetTester",
00025     src = cms.InputTag("iterativeConePu7CaloJets"),
00026     srcGen = cms.InputTag("iterativeCone7HiCleanedGenJets"), 
00027     genEnergyFractionThreshold = cms.double(0.05),
00028     genPtThreshold = cms.double(1.0),
00029     RThreshold = cms.double(0.3),
00030     reverseEnergyFractionThreshold = cms.double(0.5)                                    
00031 )
00032 
00033 JetAnalyzerAkPU5Calo = cms.EDAnalyzer("CaloJetTester",
00034     src = cms.InputTag("akPu5CaloJets"),
00035     srcGen = cms.InputTag("ak5HiCleanedGenJets"),  
00036     genEnergyFractionThreshold = cms.double(0.05),
00037     genPtThreshold = cms.double(1.0),
00038     RThreshold = cms.double(0.3),
00039     reverseEnergyFractionThreshold = cms.double(0.5)                                    
00040 )
00041 
00042 JetAnalyzerAkPU7Calo = cms.EDAnalyzer("CaloJetTester",
00043     src = cms.InputTag("akPu7CaloJets"),
00044     srcGen = cms.InputTag("ak7HiCleanedGenJets"),       
00045     genEnergyFractionThreshold = cms.double(0.05),
00046     genPtThreshold = cms.double(1.0),
00047     RThreshold = cms.double(0.3),
00048     reverseEnergyFractionThreshold = cms.double(0.5)
00049 )
00050 
00051 JetAnalyzerAkFastPU5Calo = cms.EDAnalyzer("CaloJetTester",
00052     src = cms.InputTag("akFastPu5CaloJets"),
00053     srcGen = cms.InputTag("ak5HiCleanedGenJets"),  
00054     genEnergyFractionThreshold = cms.double(0.05),
00055     genPtThreshold = cms.double(1.0),
00056     RThreshold = cms.double(0.3),
00057     reverseEnergyFractionThreshold = cms.double(0.5)                                    
00058 )
00059 
00060 JetAnalyzerAkFastPU7Calo = cms.EDAnalyzer("CaloJetTester",
00061     src = cms.InputTag("akFastPu7CaloJets"),
00062     srcGen = cms.InputTag("ak7HiCleanedGenJets"),       
00063     genEnergyFractionThreshold = cms.double(0.05),
00064     genPtThreshold = cms.double(1.0),
00065     RThreshold = cms.double(0.3),
00066     reverseEnergyFractionThreshold = cms.double(0.5)
00067 )
00068 
00069 hiJetValidation = cms.Sequence(
00070     (iterativeCone5HiCleanedGenJets * JetAnalyzerICPU5Calo) 
00071     #+ (iterativeCone7HiCleanedGenJets * JetAnalyzerICPU7Calo)
00072     #+ (ak5HiCleanedGenJets * JetAnalyzerAkPU5Calo * JetAnalyzerAkFastPU5Calo
00073     #+ (ak7HiCleanedGenJets*JetAnalyzerAkPU7Calo *JetAnalyzerAkFastPU7Calo)
00074     )