| Revision: | 1.3 |
| Committed: | Wed Jun 16 11:07:54 2010 UTC (14 years, 11 months ago) by yilmaz |
| 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.2: | +1 -1 lines |
| Log Message: | update |
| # | 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 | yilmaz | 1.3 | fileName = cms.string('treefile_' + ivars.output) |
| 47 | yilmaz | 1.1 | ) |
| 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.icPu5corr*process.icPu5match*process.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 |