JetAnalysis/test/analyzeAlgo.py

import FWCore.ParameterSet.VarParsing as VarParsing

ivars = VarParsing.VarParsing('standard')

ivars.files = 'rfio:/castor/cern.ch/user/y/yilmaz/pat/CMSSW_3_6_0/UnquenchedDijet80to120_runs1to500.root'
ivars.output = 'algo_80to120_4.root'
ivars.maxEvents = 200

ivars.parseArguments()

import FWCore.ParameterSet.Config as cms

process = cms.Process('ALGO')

process.source = cms.Source("PoolSource",
                            fileNames = cms.untracked.vstring(ivars.files),
                            skipEvents = cms.untracked.uint32(200)
                            )

process.maxEvents = cms.untracked.PSet(
            input = cms.untracked.int32(ivars.maxEvents)
                    )

'''
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cout.placeholder = cms.untracked.bool(False)
process.MessageLogger.cout.threshold = cms.untracked.string('DEBUG')
process.MessageLogger.debugModules = cms.untracked.vstring('PileUpSubtractor','JetAlgorithmAnalyzer')
'''

process.load('Configuration.StandardSequences.Services_cff')
process.load('Configuration.StandardSequences.Geometry_cff')
process.load('RecoJets.Configuration.GenJetParticles_cff')
process.load('RecoHI.HiJetAlgos.HiGenJets_cff')
process.load('RecoHI.HiJetAlgos.HiRecoJets_cff')
process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
from PhysicsTools.PatAlgos.tools.heavyIonTools import *
configureHeavyIons(process)

process.RandomNumberGeneratorService.icPu5CaloJets = process.RandomNumberGeneratorService.generator.clone()
process.RandomNumberGeneratorService.icPu5CaloJets0 = process.RandomNumberGeneratorService.generator.clone()
process.RandomNumberGeneratorService.icPu5CaloJets1 = process.RandomNumberGeneratorService.generator.clone()
process.RandomNumberGeneratorService.icPu5CaloJets2 = process.RandomNumberGeneratorService.generator.clone()

process.TFileService = cms.Service('TFileService',
                                   fileName = cms.string('treefile4.root')
                                   )

from RecoJets.JetProducers.CaloJetParameters_cfi import *
from RecoJets.JetProducers.AnomalousCellParameters_cfi import *

process.icPu5CaloJets = cms.EDProducer('JetAlgorithmAnalyzer',
                                       CaloJetParameters,
                                       AnomalousCellParameters,
                                       jetAlgorithm = cms.string("IterativeCone"),
                                       rParam       = cms.double(0.5),
                                       avoidNegative = cms.bool(False)
                                       )


process.icPu5CaloJets.doPUOffsetCorr = True
process.icPu5CaloJets.doPVCorrection = False
process.icPu5CaloJets.inputEtMin = 0
process.icPu5CaloJets.jetPtMin = 10
process.icPu5CaloJets.centrality = cms.untracked.int32(0)

#process.icPu5CaloJets.reRunAlgo = cms.untracked.bool(True)

process.icPu5CaloJets0 = process.icPu5CaloJets.clone()
process.icPu5CaloJets1 = process.icPu5CaloJets.clone(centrality = cms.untracked.int32(15))
process.icPu5CaloJets2 = process.icPu5CaloJets.clone(centrality = cms.untracked.int32(30))

process.algoAna = cms.Sequence(process.icPu5CaloJets0) # + process.icPu5CaloJets1 + process.icPu5CaloJets2)

process.icPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("icPu5CaloJets"),
                                                    corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_IC5Calo"),
                                                                          L3Absolute = cms.string("L3Absolute_IC5Calo")))
process.icPu5clean = process.heavyIonCleanedGenJets.clone(src = cms.untracked.InputTag('iterativeCone5HiGenJets'))
process.icPu5match = process.patJetGenJetMatch.clone(src = cms.InputTag("icPu5CaloJets"),
                                                     matched = cms.InputTag("heavyIonCleanedGenJets"))
process.icPu5parton = process.patJetPartonMatch.clone(src = cms.InputTag("icPu5CaloJets"))
process.icPu5patJets = process.patJets.clone(jetSource  = cms.InputTag("icPu5CaloJets"),
                                             genJetMatch = cms.InputTag("icPu5match"),
                                             genPartonMatch= cms.InputTag("icPu5parton"),
                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")))
process.icPu5patSequence = cms.Sequence(process.icPu5corr*process.icPu5match*process.icPu5parton*process.icPu5patJets)

process.load("HeavyIonsAnalysis.Configuration.analysisEventContent_cff")

process.output = cms.OutputModule("PoolOutputModule",
                                  process.jetTrkSkimContent,
                                  fileName = cms.untracked.string(ivars.output)
                                  )

process.output.outputCommands.append("keep *_*_*_ALGO")


process.load("CondCore.DBCommon.CondDBCommon_cfi")
process.CondDBCommon.connect = "sqlite_file:/home/yilmaz/cms361v2/src/RecoHI/HiCentralityAlgos/data/CentralityTables.db"

process.PoolDBESSource = cms.ESSource("PoolDBESSource",
                                      process.CondDBCommon,
                                      toGet = cms.VPSet(cms.PSet(record = cms.string('HeavyIonRcd'),
                                                                 tag = cms.string('HFhits40_MC_Hydjet2760GeV_MC_3XY_V24_v0')
                                                                 )
                                                        )
                                      )


process.runAllJets = cms.Path(process.hiGenParticlesForJets *
                              process.hiRecoGenJets +
                              process.iterativeConePu5CaloJets *
                              process.makeHeavyIonJets +
                              process.algoAna # *
#                              process.icPu5patSequence
                              )

process.out_step = cms.EndPath(process.output)

Revision:	1.2
Committed:	Mon Jun 7 14:27:28 2010 UTC (14 years, 11 months ago) by yilmaz
Content type:	text/x-python
Branch:	MAIN
Changes since 1.1:	+9 -1 lines
Log Message:	debug stuff
#	User	Rev	Content
1	yilmaz	1.1	import FWCore.ParameterSet.VarParsing as VarParsing
2
3			ivars = VarParsing.VarParsing('standard')
4
5			ivars.files = 'rfio:/castor/cern.ch/user/y/yilmaz/pat/CMSSW_3_6_0/UnquenchedDijet80to120_runs1to500.root'
6			ivars.output = 'algo_80to120_4.root'
7			ivars.maxEvents = 200
8
9			ivars.parseArguments()
10
11			import FWCore.ParameterSet.Config as cms
12
13			process = cms.Process('ALGO')
14
15			process.source = cms.Source("PoolSource",
16			fileNames = cms.untracked.vstring(ivars.files),
17			skipEvents = cms.untracked.uint32(200)
18			)
19
20			process.maxEvents = cms.untracked.PSet(
21			input = cms.untracked.int32(ivars.maxEvents)
22			)
23
24	yilmaz	1.2	'''
25			process.load("FWCore.MessageLogger.MessageLogger_cfi")
26			process.MessageLogger.cout.placeholder = cms.untracked.bool(False)
27			process.MessageLogger.cout.threshold = cms.untracked.string('DEBUG')
28			process.MessageLogger.debugModules = cms.untracked.vstring('PileUpSubtractor','JetAlgorithmAnalyzer')
29			'''
30
31	yilmaz	1.1	process.load('Configuration.StandardSequences.Services_cff')
32			process.load('Configuration.StandardSequences.Geometry_cff')
33			process.load('RecoJets.Configuration.GenJetParticles_cff')
34			process.load('RecoHI.HiJetAlgos.HiGenJets_cff')
35			process.load('RecoHI.HiJetAlgos.HiRecoJets_cff')
36			process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
37			from PhysicsTools.PatAlgos.tools.heavyIonTools import *
38			configureHeavyIons(process)
39
40			process.RandomNumberGeneratorService.icPu5CaloJets = process.RandomNumberGeneratorService.generator.clone()
41			process.RandomNumberGeneratorService.icPu5CaloJets0 = process.RandomNumberGeneratorService.generator.clone()
42			process.RandomNumberGeneratorService.icPu5CaloJets1 = process.RandomNumberGeneratorService.generator.clone()
43			process.RandomNumberGeneratorService.icPu5CaloJets2 = process.RandomNumberGeneratorService.generator.clone()
44
45			process.TFileService = cms.Service('TFileService',
46			fileName = cms.string('treefile4.root')
47			)
48
49			from RecoJets.JetProducers.CaloJetParameters_cfi import *
50			from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
51
52			process.icPu5CaloJets = cms.EDProducer('JetAlgorithmAnalyzer',
53			CaloJetParameters,
54			AnomalousCellParameters,
55			jetAlgorithm = cms.string("IterativeCone"),
56			rParam = cms.double(0.5),
57			avoidNegative = cms.bool(False)
58			)
59
60
61			process.icPu5CaloJets.doPUOffsetCorr = True
62			process.icPu5CaloJets.doPVCorrection = False
63	yilmaz	1.2	process.icPu5CaloJets.inputEtMin = 0
64	yilmaz	1.1	process.icPu5CaloJets.jetPtMin = 10
65			process.icPu5CaloJets.centrality = cms.untracked.int32(0)
66
67			#process.icPu5CaloJets.reRunAlgo = cms.untracked.bool(True)
68
69			process.icPu5CaloJets0 = process.icPu5CaloJets.clone()
70			process.icPu5CaloJets1 = process.icPu5CaloJets.clone(centrality = cms.untracked.int32(15))
71			process.icPu5CaloJets2 = process.icPu5CaloJets.clone(centrality = cms.untracked.int32(30))
72
73	yilmaz	1.2	process.algoAna = cms.Sequence(process.icPu5CaloJets0) # + process.icPu5CaloJets1 + process.icPu5CaloJets2)
74	yilmaz	1.1
75			process.icPu5corr = process.patJetCorrFactors.clone(jetSource = cms.InputTag("icPu5CaloJets"),
76			corrLevels = cms.PSet(L2Relative = cms.string("L2Relative_IC5Calo"),
77			L3Absolute = cms.string("L3Absolute_IC5Calo")))
78			process.icPu5clean = process.heavyIonCleanedGenJets.clone(src = cms.untracked.InputTag('iterativeCone5HiGenJets'))
79			process.icPu5match = process.patJetGenJetMatch.clone(src = cms.InputTag("icPu5CaloJets"),
80			matched = cms.InputTag("heavyIonCleanedGenJets"))
81			process.icPu5parton = process.patJetPartonMatch.clone(src = cms.InputTag("icPu5CaloJets"))
82			process.icPu5patJets = process.patJets.clone(jetSource = cms.InputTag("icPu5CaloJets"),
83			genJetMatch = cms.InputTag("icPu5match"),
84			genPartonMatch= cms.InputTag("icPu5parton"),
85			jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")))
86			process.icPu5patSequence = cms.Sequence(process.icPu5corrprocess.icPu5matchprocess.icPu5parton*process.icPu5patJets)
87
88			process.load("HeavyIonsAnalysis.Configuration.analysisEventContent_cff")
89
90			process.output = cms.OutputModule("PoolOutputModule",
91			process.jetTrkSkimContent,
92			fileName = cms.untracked.string(ivars.output)
93			)
94
95			process.output.outputCommands.append("keep __*_ALGO")
96
97
98			process.load("CondCore.DBCommon.CondDBCommon_cfi")
99			process.CondDBCommon.connect = "sqlite_file:/home/yilmaz/cms361v2/src/RecoHI/HiCentralityAlgos/data/CentralityTables.db"
100
101			process.PoolDBESSource = cms.ESSource("PoolDBESSource",
102			process.CondDBCommon,
103			toGet = cms.VPSet(cms.PSet(record = cms.string('HeavyIonRcd'),
104			tag = cms.string('HFhits40_MC_Hydjet2760GeV_MC_3XY_V24_v0')
105			)
106			)
107			)
108
109
110			process.runAllJets = cms.Path(process.hiGenParticlesForJets *
111			process.hiRecoGenJets +
112			process.iterativeConePu5CaloJets *
113			process.makeHeavyIonJets +
114			process.algoAna # *
115			# process.icPu5patSequence
116			)
117
118			process.out_step = cms.EndPath(process.output)
119