ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/test/analyzeTowers.py
(Generate patch)

Comparing UserCode/CmsHi/JetAnalysis/test/analyzeTowers.py (file contents):
Revision 1.6 by yilmaz, Mon Jul 12 14:03:26 2010 UTC vs.
Revision 1.12 by nart, Thu Aug 12 21:13:13 2010 UTC

# Line 1 | Line 1
1 +
2   import FWCore.ParameterSet.VarParsing as VarParsing
3   ivars = VarParsing.VarParsing('standard')
4 + '''
5   ivars.files = [
6 <   # 'rfio:/castor/cern.ch/user/y/yilmaz/pat/CMSSW_3_7_0/SignalQuenchedDijet80to120_runs1to100.root'
7 <    'rfio:/castor/cern.ch/user/y/yilmaz/share/MinBias0708_runs1to50.root'
8 <              ]
6 >    'file:/net/hisrv0001/home/yetkin/hibat0007/analysis/jets/SignalUnquenchedDijet80to120_runs101to150.root',
7 >    'file:/net/hisrv0001/home/yetkin/hibat0007/analysis/jets/SignalUnquenchedDijet80to120_runs151to200.root',
8 >    'file:/net/hisrv0001/home/yetkin/hibat0007/analysis/jets/SignalUnquenchedDijet80to120_runs1to50.root',
9 >    'file:/net/hisrv0001/home/yetkin/hibat0007/analysis/jets/SignalUnquenchedDijet80to120_runs51to100.root'
10 >    ]
11 > '''
12  
13 < ivars.output = 'RandomCones_Hydjet_370.root'
13 > ivars.files = [
14 >    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100713/Hydjet_MinBias_2760GeV_runs11to20.root',
15 >    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100713/Hydjet_MinBias_2760GeV_runs21to40.root',
16 >    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100713/Hydjet_MinBias_2760GeV_runs61to80.root'
17 >    ]
18 >    
19 > '''
20 > ivars.files = [
21 > 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs11to20.root',
22 > 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs1to10.root',
23 > 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs1to5.root',
24 > 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs21to30.root',
25 > 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs31to40.root',
26 > 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs41to50.root',
27 > 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs51to100.root'
28 >    ]
29 > '''
30 >
31 > '''
32 > ivars.files = [
33 >    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs101to200.root',
34 >    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs1to100.root',
35 >    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs201to300.root',
36 >    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs301to400.root'
37 >    ]
38 > '''
39  
40 < ivars.maxEvents = 10000
40 > ivars.output = 'data0628_multiplePU_avoidT.root'
41 >
42 > ivars.maxEvents = -1
43  
44   ivars.register ('randomNumber',
45                  mult=ivars.multiplicity.singleton,
# Line 44 | Line 76 | process.output = cms.OutputModule("PoolO
76  
77   process.load('Configuration/StandardSequences/GeometryExtended_cff')
78   process.load('Configuration/StandardSequences/Services_cff')
79 + process.load('Configuration/StandardSequences/ReconstructionHeavyIons_cff')
80  
81   process.RandomNumberGeneratorService.bkg4Jets = cms.PSet(initialSeed = cms.untracked.uint32(ivars.randomNumber),
82                                                           engineName = cms.untracked.string('HepJamesRandom') )
83   process.RandomNumberGeneratorService.bkg5Jets = process.RandomNumberGeneratorService.bkg4Jets.clone()
84   process.RandomNumberGeneratorService.bkg6Jets = process.RandomNumberGeneratorService.bkg4Jets.clone()
85   process.RandomNumberGeneratorService.bkg7Jets = process.RandomNumberGeneratorService.bkg4Jets.clone()
86 + process.RandomNumberGeneratorService.bkgXJets = process.RandomNumberGeneratorService.bkg4Jets.clone()
87  
88   from RecoJets.JetProducers.CaloJetParameters_cfi import *
89   from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
90  
91 + process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
92 + from PhysicsTools.PatAlgos.tools.heavyIonTools import *
93 + configureHeavyIons(process)
94 +
95 + process.patJets.addJetCorrFactors = True
96 + process.patJets.addGenPartonMatch  = False
97 + process.patJets.addJetID = False
98 + process.patJets.addGenJetMatch = False
99 + process.patJets.embedGenJetMatch = False
100 +
101 + process.icPu5corr = process.patJetCorrFactors.clone()
102 + process.icPu5corr.jetSource = cms.InputTag("iterativeConePu5CaloJets","")
103 + process.icPu5corr.corrLevels.L2Relative = cms.string('L2Relative_IC5Calo')
104 + process.icPu5corr.corrLevels.L3Absolute = cms.string('L3Absolute_IC5Calo')
105 +
106 + process.icPu5patJets = process.patJets.clone()
107 + process.icPu5patJets.jetSource = cms.InputTag("iterativeConePu5CaloJets","")
108 + #process.icPu5patJets.genJetMatch = cms.InputTag("icPu5match")
109 + process.icPu5patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr") )
110 +
111 +
112   process.bkg4Jets = cms.EDProducer(
113 < #    "JetAlgorithmAnalyzer",
114 <  "BackgroundJetProducer",
115 <  CaloJetParameters,
116 <  AnomalousCellParameters,
117 <  jetAlgorithm = cms.string("IterativeCone"),
118 <  rParam       = cms.double(0.4),
119 <  avoidNegative = cms.bool(False),
120 <  subtractorName = cms.string("MultipleAlgoIterator")
121 <  )
113 >    "BackgroundJetProducer",
114 >    CaloJetParameters,
115 >    AnomalousCellParameters,
116 >    jetAlgorithm = cms.string("IterativeCone"),
117 >    rParam       = cms.double(0.4),
118 >    avoidNegative = cms.bool(True),
119 >    subtractorName = cms.string("MultipleAlgoIterator")
120 >    #    subtractorName = cms.string("ParametrizedSubtractor")
121 >    )
122 >
123 > process.bkgXJets = cms.EDProducer(
124 >    "JetAlgorithmAnalyzer",
125 >    CaloJetParameters,
126 >    AnomalousCellParameters,
127 >    jetAlgorithm = cms.string("IterativeCone"),
128 >    rParam       = cms.double(0.5),
129 >    avoidNegative = cms.bool(True),
130 >    #subtractorName = cms.string("ParametrizedSubtractor"),
131 >    subtractorName = cms.string("MultipleAlgoIterator"),
132 >    centralityTag =cms.InputTag("hiCentrality"),
133 >    evtPlaneTag = cms.InputTag("hiEvtPlane"),
134 >    doRecoEvtPlane = cms.untracked.bool(False),
135 >    interpolate = cms.untracked.bool(False),
136 >    centrality = cms.untracked.int32(-1),
137 >    sumRecHits = cms.bool(False)
138 >    )
139 >
140 > process.bkgXJets.doPUOffsetCorr = True
141 > process.bkgXJets.doPVCorrection = False
142 > process.bkgXJets.jetPtMin = 10
143 > process.bkgXJets.inputEtMin = -10.
144 > process.bkgXJets.radiusPU = 0.5
145 > process.bkgXJets.nSigmaPU = 1.
146  
147   process.bkg4Jets.doPUOffsetCorr = True
148   process.bkg4Jets.doPVCorrection = False
149   process.bkg4Jets.jetPtMin = 10
150 + process.bkg4Jets.inputEtMin = -10.
151   process.bkg4Jets.radiusPU = 0.4
152 + process.bkg4Jets.nSigmaPU = 1.
153  
154   process.bkg5Jets = process.bkg4Jets.clone()
155   process.bkg5Jets.rParam = 0.5
# Line 84 | Line 165 | process.bkg7Jets.radiusPU = 0.7
165  
166   process.bkgJets = cms.Sequence(process.bkg5Jets)
167  
168 +
169   process.ana = cms.EDAnalyzer('MinBiasTowerAnalyzer',
170                               jetTowersMean = cms.untracked.vdouble(21,21.4,21.4,22,22.5,21.3,17.4,16.1,11.5,0),
171 <                             jetTowersRMS = cms.untracked.vdouble(5.4,5.2,5.2,5.4,5.8,5.8,4.9,4.1,3.2,0),          
172 <                             fakeJetSrc = cms.untracked.InputTag('bkg5Jets'),
173 <                             patJetSrc = cms.untracked.InputTag('patJets')                            
174 <                             )
171 >                             jetTowersRMS = cms.untracked.vdouble(5.4,5.2,5.2,5.4,5.8,5.8,4.9,4.1,3.2,0),
172 >                             fakeJetSrc = cms.untracked.InputTag("bkg5Jets"),
173 >                             patJetSrc = cms.untracked.InputTag("icPu5patJets","","ANALYSIS"),  
174 >                             towersSrc =  cms.untracked.InputTag("towerMaker","","RECO"),
175 >                             nBins = cms.untracked.double(10),
176 >                             doRandomCone = cms.untracked.bool(True),
177 >                             doEvtPlane = cms.untracked.bool(True),
178 >                             centralitySrc = cms.untracked.InputTag("hiCentrality","","RECO"),
179 >                             excludeJets = cms.untracked.bool(False)
180  
181 +                             )
182  
183   process.load("CondCore.DBCommon.CondDBCommon_cfi")
184 < process.CondDBCommon.connect = "sqlite_file:/afs/cern.ch/user/n/nart/scratch0/CMSSW_3_7_0_patch4/src/CmsHi/JulyExercise/data/CentralityTables.db"
184 > process.CondDBCommon.connect = "sqlite_file:/net/hisrv0001/home/nart/scratch/july/CMSSW_3_7_0_patch4/src/RecoHI/HiCentralityAlgos/data/CentralityTables.db"
185   process.PoolDBESSource = cms.ESSource("PoolDBESSource",
186                                        process.CondDBCommon,
187                                        toGet = cms.VPSet(cms.PSet(record = cms.string('HeavyIonRcd'),
188 <                                                                 tag = cms.string('HFhits10_DataJulyExercise_AMPT2760GeV_MC_37Y_V5_v0')
188 >                                                                 tag = cms.string('HFhits10_MC_Hydjet2760GeV_MC_3XY_V24_NoZS_v0')
189                                                                   )
190                                                          )
191                                        )
192  
193 < process.p = cms.Path(process.bkgJets
194 <                     *process.ana
195 <                     )
196 < process.out_step = cms.EndPath(process.output)
197 <
198 <
193 > process.p = cms.Path(
194 >    process.bkgXJets +
195 >    process.hiEvtPlane +
196 >    process.ak5CaloJets +
197 >    process.kt4CaloJets +
198 >    process.bkgJets +
199 >    process.icPu5corr +
200 >    process.icPu5patJets +
201 >    process.ana
202 >    )
203 > #process.out_step = cms.EndPath(process.output)
204  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines