| 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 |
| # | Content |
|---|---|
| 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 | ''' |
| 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 | 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('treefile_' + ivars.output) |
| 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 | process.icPu5CaloJets.inputEtMin = 0 |
| 64 | 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 | process.algoAna = cms.Sequence(process.icPu5CaloJets0) # + process.icPu5CaloJets1 + process.icPu5CaloJets2) |
| 74 | |
| 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 |