ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/test/analyzeData.py
Revision: 1.2
Committed: Sat Oct 23 21:09:53 2010 UTC (14 years, 6 months ago) by nart
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, hi44X_02, hi413_03, hi441_1, hi441_0, hi413_11, hi413_10, hi413_09, hi413_08, hi413_07, hi413_06, hi413_05, hi413_04, hi413_02, hi39X_01, tag_d20110915, cmssw39x_base, cmssw39X_base
Branch point for: branch_44x, cmssw39x_branch
Changes since 1.1: +27 -4 lines
Log Message:
jet stuff

File Contents

# User Rev Content
1 nart 1.1
2     import FWCore.ParameterSet.VarParsing as VarParsing
3     ivars = VarParsing.VarParsing('standard')
4    
5     ivars.files = [
6     # ZS Data:
7 nart 1.2 # 'rfio:/castor/cern.ch/cms/store/data/Run2010B/MinimumBias/RECO/PromptReco-v2/000/146/510/82BB6441-6EC7-DF11-BED4-001617C3B66C.root',
8 nart 1.1 #NZS Data:
9     'rfio:/castor/cern.ch/cms/store/data/Run2010B/HeavyIonTest/RECO/PromptReco-v2/000/146/421/FE8B65B2-ABC6-DF11-8534-0030487CD7E0.root',
10     ]
11    
12    
13     ivars.output = 'Data.root'
14    
15     ivars.maxEvents = 10
16    
17     ivars.register ('randomNumber',
18     mult=ivars.multiplicity.singleton,
19     info="for testing")
20     ivars.randomNumber=5
21     ivars.parseArguments()
22    
23     import FWCore.ParameterSet.Config as cms
24    
25     process = cms.Process('ANALYSIS')
26    
27     process.maxEvents = cms.untracked.PSet(
28     input = cms.untracked.int32(ivars.maxEvents)
29     )
30    
31     # Input source
32     process.source = cms.Source("PoolSource",
33     fileNames = cms.untracked.vstring(ivars.files)
34 nart 1.2 )
35 nart 1.1
36     process.TFileService = cms.Service('TFileService',
37     fileName = cms.string('plots_' + ivars.output)
38     )
39    
40    
41     process.output = cms.OutputModule("PoolOutputModule",
42 nart 1.2 outputCommands = cms.untracked.vstring("keep *_*_*_*"),
43 nart 1.1 fileName = cms.untracked.string(ivars.output)
44     )
45    
46    
47     process.load('Configuration/StandardSequences/GeometryExtended_cff')
48     process.load('Configuration/StandardSequences/Services_cff')
49 nart 1.2 process.load('Configuration/StandardSequences/ReconstructionHeavyIons_cff')
50    
51     process.load('PhysicsTools.PatAlgos.patHeavyIonSequences_cff')
52     from PhysicsTools.PatAlgos.tools.heavyIonTools import *
53     configureHeavyIons(process)
54    
55     process.patJets.addJetCorrFactors = True
56     process.patJets.addGenPartonMatch = False
57     process.patJets.addJetID = False
58     process.patJets.addGenJetMatch = False
59     process.patJets.embedGenJetMatch = False
60    
61     process.icPu5corr = process.patJetCorrFactors.clone()
62     process.icPu5corr.jetSource = cms.InputTag("iterativeCone5CaloJets","")
63     process.icPu5corr.corrLevels.L2Relative = cms.string('L2Relative_IC5Calo')
64     process.icPu5corr.corrLevels.L3Absolute = cms.string('L3Absolute_IC5Calo')
65    
66     process.icPu5patJets = process.patJets.clone()
67     process.icPu5patJets.jetSource = cms.InputTag("iterativeCone5CaloJets","")
68     process.icPu5patJets.jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr") )
69 nart 1.1
70     process.ana = cms.EDAnalyzer('MinBiasTowerAnalyzer',
71     nBins = cms.untracked.int32(1),
72 nart 1.2 patJetSrc = cms.untracked.InputTag("icPu5patJets","","ANALYSIS"),
73 nart 1.1 towersSrc = cms.untracked.InputTag("towerMaker","","RECO"),
74     doRandomCone = cms.untracked.bool(False),
75     doEvtPlane = cms.untracked.bool(False),
76 nart 1.2 excludeJets = cms.untracked.bool(False),
77 nart 1.1 isSignal = cms.untracked.bool(True),
78     doMC = cms.untracked.bool(False),
79     doTowers = cms.untracked.bool(True),
80     doRecHits = cms.untracked.bool(True),
81     ecalEBRecHitSrc = cms.untracked.InputTag("ecalRecHit","EcalRecHitsEB","RECO"),
82     ecalEERecHitSrc = cms.untracked.InputTag("ecalRecHit","EcalRecHitsEE","RECO")
83     )
84    
85     process.pana = cms.Path(
86 nart 1.2 process.icPu5corr +
87     process.icPu5patJets +
88 nart 1.1 process.ana
89     )