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.1 by yilmaz, Mon Jul 5 22:18:11 2010 UTC vs.
Revision 1.13 by nart, Mon Aug 16 20:51:21 2010 UTC

# Line 1 | Line 1
1 +
2 + import FWCore.ParameterSet.VarParsing as VarParsing
3 + ivars = VarParsing.VarParsing('standard')
4 + '''
5 + ivars.files = [
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.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 + 'dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/yetkin/reco/CMSSW_3_8_1_patch1/MC_38Y_V8/Hydjet_MinBias_2760GeV_d20100812/Hydjet_Quenched_MinBias_2760GeV_cfi_py_GEN_SIM_DIGI_L1_DIGI2RAW_HLT_1.root',
22 + 'dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/yetkin/reco/CMSSW_3_8_1_patch1/MC_38Y_V8/Hydjet_MinBias_2760GeV_d20100812/Hydjet_Quenched_MinBias_2760GeV_cfi_py_GEN_SIM_DIGI_L1_DIGI2RAW_HLT_10.root',
23 + 'dcache:/pnfs/cmsaf.mit.edu/t2bat/cms/store/user/yetkin/reco/CMSSW_3_8_1_patch1/MC_38Y_V8/Hydjet_MinBias_2760GeV_d20100812/Hydjet_Quenched_MinBias_2760GeV_cfi_py_GEN_SIM_DIGI_L1_DIGI2RAW_HLT_100.root'
24 + ]
25 +
26 + '''
27 + ivars.files = [
28 + 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs11to20.root',
29 + 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs1to10.root',
30 + 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs1to5.root',
31 + 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs21to30.root',
32 + 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs31to40.root',
33 + 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs41to50.root',
34 + 'file:/net/hisrv0001/home/yetkin/hibat0007/aod/JulyExercise/MinBias0709/MinBias0709_runs51to100.root'
35 +    ]
36 + '''
37 +
38 + '''
39 + ivars.files = [
40 +    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs101to200.root',
41 +    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs1to100.root',
42 +    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs201to300.root',
43 +    'file:/net/hisrv0001/home/yetkin/pstore02/reco/NoZSP/local2/Hydjet_MinBias_2760GeV_d20100628/Hydjet_MinBias_2760GeV_runs301to400.root'
44 +    ]
45 + '''
46 +
47 + ivars.output = 'Hydjet_Quenched_2760GeV.root'
48 +
49 + ivars.maxEvents = -1
50 +
51 + ivars.register ('randomNumber',
52 +                mult=ivars.multiplicity.singleton,
53 +                info="for testing")
54 + ivars.randomNumber=5
55 + ivars.parseArguments()
56 +
57 + import FWCore.ParameterSet.Config as cms
58 +
59 +
60   process = cms.Process('ANALYSIS')
61  
62   process.maxEvents = cms.untracked.PSet(
63 <    input = cms.untracked.int32(2)
64 < )
63 >    input = cms.untracked.int32(ivars.maxEvents)  
64 >    )
65 >
66  
67   # Input source
68   process.source = cms.Source("PoolSource",
69 <                            fileNames = cms.untracked.vstring('rfio:/castor/cern.ch/user/y/yilmaz/pat/CMSSW_3_7_0/Hydjet_MinBias_2760GeV_000.root')
10 <                            )
69 >                            fileNames = cms.untracked.vstring(ivars.files)
70  
71 < process.ana = cms.EDAnalyzer('MinBiasTowerAnalyzer')
71 >                          )
72  
73   process.TFileService = cms.Service('TFileService',
74 <                                   fileName = cms.string('plots.root')
74 >                                   fileName = cms.string('plots_' + ivars.output)
75                                     )
76  
77 < process.p = cms.Path(process.ana)
77 >
78 > process.output = cms.OutputModule("PoolOutputModule",
79 >                                   outputCommands = cms.untracked.vstring("keep *_*_*_*"),
80 >                                  fileName = cms.untracked.string(ivars.output)
81 >                                  )
82 >
83 >
84 > process.load('Configuration/StandardSequences/GeometryExtended_cff')
85 > process.load('Configuration/StandardSequences/Services_cff')
86 > process.load('Configuration/StandardSequences/ReconstructionHeavyIons_cff')
87 >
88 > process.RandomNumberGeneratorService.bkg4Jets = cms.PSet(initialSeed = cms.untracked.uint32(ivars.randomNumber),
89 >                                                         engineName = cms.untracked.string('HepJamesRandom') )
90 > process.RandomNumberGeneratorService.bkg5Jets = process.RandomNumberGeneratorService.bkg4Jets.clone()
91 > process.RandomNumberGeneratorService.bkg6Jets = process.RandomNumberGeneratorService.bkg4Jets.clone()
92 > process.RandomNumberGeneratorService.bkg7Jets = process.RandomNumberGeneratorService.bkg4Jets.clone()
93 > process.RandomNumberGeneratorService.bkgXJets = process.RandomNumberGeneratorService.bkg4Jets.clone()
94 >
95 > from RecoJets.JetProducers.CaloJetParameters_cfi import *
96 > from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
97 >
98 > process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
99 > from PhysicsTools.PatAlgos.tools.heavyIonTools import *
100 > configureHeavyIons(process)
101 >
102 > process.patJets.addJetCorrFactors = True
103 > process.patJets.addGenPartonMatch  = False
104 > process.patJets.addJetID = False
105 > process.patJets.addGenJetMatch = False
106 > process.patJets.embedGenJetMatch = False
107 >
108 > process.icPu5corr = process.patJetCorrFactors.clone()
109 > process.icPu5corr.jetSource = cms.InputTag("iterativeConePu5CaloJets","")
110 > process.icPu5corr.corrLevels.L2Relative = cms.string('L2Relative_IC5Calo')
111 > process.icPu5corr.corrLevels.L3Absolute = cms.string('L3Absolute_IC5Calo')
112 >
113 > process.icPu5patJets = process.patJets.clone()
114 > process.icPu5patJets.jetSource = cms.InputTag("iterativeConePu5CaloJets","")
115 > #process.icPu5patJets.genJetMatch = cms.InputTag("icPu5match")
116 > process.icPu5patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr") )
117 >
118 >
119 > process.bkg4Jets = cms.EDProducer(
120 >    "BackgroundJetProducer",
121 >    CaloJetParameters,
122 >    AnomalousCellParameters,
123 >    jetAlgorithm = cms.string("IterativeCone"),
124 >    rParam       = cms.double(0.4),
125 >    avoidNegative = cms.bool(True),
126 >    subtractorName = cms.string("MultipleAlgoIterator")
127 >    #    subtractorName = cms.string("ParametrizedSubtractor")
128 >    )
129 >
130 > process.bkgXJets = cms.EDProducer(
131 >    "JetAlgorithmAnalyzer",
132 >    CaloJetParameters,
133 >    AnomalousCellParameters,
134 >    jetAlgorithm = cms.string("IterativeCone"),
135 >    rParam       = cms.double(0.5),
136 >    avoidNegative = cms.bool(True),
137 >    #subtractorName = cms.string("ParametrizedSubtractor"),
138 >    subtractorName = cms.string("MultipleAlgoIterator"),
139 >    centralityTag =cms.InputTag("hiCentrality"),
140 >    evtPlaneTag = cms.InputTag("hiEvtPlane"),
141 >    doRecoEvtPlane = cms.untracked.bool(False),
142 >    interpolate = cms.untracked.bool(False),
143 >    centrality = cms.untracked.int32(-1),
144 >    sumRecHits = cms.bool(False)
145 >    )
146 >
147 > process.bkgXJets.doPUOffsetCorr = True
148 > process.bkgXJets.doPVCorrection = False
149 > process.bkgXJets.jetPtMin = 10
150 > process.bkgXJets.inputEtMin = -10.
151 > process.bkgXJets.radiusPU = 0.5
152 > process.bkgXJets.nSigmaPU = 1.
153 >
154 > process.bkg4Jets.doPUOffsetCorr = True
155 > process.bkg4Jets.doPVCorrection = False
156 > process.bkg4Jets.jetPtMin = 10
157 > process.bkg4Jets.inputEtMin = -10.
158 > process.bkg4Jets.radiusPU = 0.4
159 > process.bkg4Jets.nSigmaPU = 1.
160 >
161 > process.bkg5Jets = process.bkg4Jets.clone()
162 > process.bkg5Jets.rParam = 0.5
163 > process.bkg5Jets.radiusPU = 0.5
164 >
165 > process.bkg6Jets = process.bkg4Jets.clone()
166 > process.bkg6Jets.rParam = 0.6
167 > process.bkg6Jets.radiusPU = 0.6
168 >
169 > process.bkg7Jets = process.bkg4Jets.clone()
170 > process.bkg7Jets.rParam = 0.7
171 > process.bkg7Jets.radiusPU = 0.7
172 >
173 > process.bkgJets = cms.Sequence(process.bkg5Jets)
174 >
175 >
176 > process.ana = cms.EDAnalyzer('MinBiasTowerAnalyzer',
177 >                             jetTowersMean = cms.untracked.vdouble(21,21.4,21.4,22,22.5,21.3,17.4,16.1,11.5,0),
178 >                             jetTowersRMS = cms.untracked.vdouble(5.4,5.2,5.2,5.4,5.8,5.8,4.9,4.1,3.2,0),
179 >                             fakeJetSrc = cms.untracked.InputTag("bkg5Jets"),
180 >                             patJetSrc = cms.untracked.InputTag("icPu5patJets","","ANALYSIS"),  
181 >                             towersSrc =  cms.untracked.InputTag("towerMaker","","RECO"),
182 >                             nBins = cms.untracked.double(10),
183 >                             doRandomCone = cms.untracked.bool(True),
184 >                             doEvtPlane = cms.untracked.bool(True),
185 >                             centralitySrc = cms.untracked.InputTag("hiCentrality","","RECO"),
186 >                             excludeJets = cms.untracked.bool(False),
187 >                             doRecHits = cms.untracked.bool(True),
188 >                             ecalEBRecHitSrc = cms.untracked.InputTag("ecalRecHit","EcalRecHitsEB","RECO"),
189 >                             ecalEERecHitSrc = cms.untracked.InputTag("ecalRecHit","EcalRecHitsEE","RECO")
190 >                             )
191 >
192 > process.load("CondCore.DBCommon.CondDBCommon_cfi")
193 > process.CondDBCommon.connect = "sqlite_file:/net/hisrv0001/home/nart/scratch/july/CMSSW_3_7_0_patch4/src/RecoHI/HiCentralityAlgos/data/CentralityTables.db"
194 > process.PoolDBESSource = cms.ESSource("PoolDBESSource",
195 >                                      process.CondDBCommon,
196 >                                      toGet = cms.VPSet(cms.PSet(record = cms.string('HeavyIonRcd'),
197 >                                                                 tag = cms.string('HFhits10_MC_Hydjet2760GeV_MC_3XY_V24_NoZS_v0')
198 >                                                                 )
199 >                                                        )
200 >                                      )
201 >
202 > process.p = cms.Path(
203 >    process.bkgXJets +
204 >    process.hiEvtPlane +
205 >    process.ak5CaloJets +
206 >    process.kt4CaloJets +
207 >    process.bkgJets +
208 >    process.icPu5corr +
209 >    process.icPu5patJets +
210 >    process.ana
211 >    )
212 > #process.out_step = cms.EndPath(process.output)
213  

Diff Legend

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