HcalCalibAlgos/python/isoTrackCalib.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("hcalCalib")
0004
0005 process.load("Configuration.StandardSequences.GeometryHCAL_cff")
0006
0007 # no CMS input files
0008 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(0) )
0009 process.source = cms.Source("EmptySource")
0010
0011
0012 process.hcalCalib = cms.EDAnalyzer("HcalCalibrator",
0013
0014
0015
0016 #-------- File containing the list of the input root files
0017   inputFileList = cms.untracked.string("test/inputFiles_IsoTrkCalib.list"),
0018
0019 #-------- calibration type: ISO_TRACK or DI_JET
0020   calibType = cms.untracked.string("ISO_TRACK"),
0021
0022 #-------- calibration method: L3, MATRIX_INV_OF_ETA_AVE, L3_AND_MTRX_INV
0023   calibMethod = cms.untracked.string("L3_AND_MTRX_INV"),
0024
0025 # --------min target energy: ref jet (had) energy, or isotrack momentum
0026   minTargetE = cms.untracked.double(40.0),
0027
0028 # ------- max target energy: ref jet (had) energy, or isotrack momentum
0029   maxTargetE = cms.untracked.double(60.0),
0030
0031 # ------- min energy in the cell
0032   minCellE = cms.untracked.double(0.0),
0033
0034 # ------- minimum e/p to accept an IsoTrack for calibration (not used for DiJEts)
0035   minEOverP = cms.untracked.double(0.2),
0036
0037 # ------- maximum e/p to accept an IsoTrack for calibration (not used for DiJEts)
0038   maxEOverP = cms.untracked.double(999.0),
0039
0040 # ------- MIP cut: maximum EM energy associated with the track (not used for DiJEts)
0041   maxTrkEmE = cms.untracked.double(1.0),
0042
0043 #--------Maximum Et allowed for third jet in the event (affects only DiJets)
0044   maxEtThirdJet = cms.untracked.double(5.0),
0045
0046 # ------- Minimum deltaPhi for the dijets in degrees: 0-180  (affects only DiJets)
0047   minDPhiDiJets = cms.untracked.double(150.0),
0048
0049 # ------- logical flag to sum depths in HCAL
0050   sumDepths = cms.untracked.bool(True),
0051
0052 # ------- logical flag to sum depths 1,2 in HB for towers 15,16.
0053 # --------If sumDepths="true" this flag has no effect
0054   sumSmallDepths = cms.untracked.bool(True),
0055
0056 # ------- logical flag to combine phi in HCAL
0057   combinePhi = cms.untracked.bool(True),
0058
0059 # ------- cluster size in HB for isotracks: 3 or 5 (means 3x3, 5x5)
0060   hbClusterSize = cms.untracked.int32(3),
0061
0062 # ------- cluster size in HE for isotracks: 3 or 5 (means 3x3, 5x5)
0063   heClusterSize = cms.untracked.int32(3),
0064
0065 # -------- flag to use cone clustering -> overrides the above cluster sizes
0066 # -------- and uses cone size as specified below
0067  useConeClustering = cms.untracked.bool(True),
0068
0069 # -------- size of the cone (when useConeClustering=True)
0070  maxConeDist = cms.untracked.double(26.2),
0071
0072 # ------- max ABS(iEta) used in the calibration: for matrix inversion sets the range for performing inversion
0073 # ------- For all methods: controls the range of correction factors are saved in the output file
0074   calibAbsIEtaMax  = cms.untracked.int32(24),
0075
0076 # ------- max ABS(iEta) used in the calibration: for matrix inversion sets the range for performing inversion
0077 # ------- For all methods: controls the range of correction factors are saved in the output file
0078   calibAbsIEtaMin  = cms.untracked.int32(1),
0079
0080 # ------- max EmFraction of probe jet allowed (does not affect isotracks)
0081   maxProbeJetEmFrac  = cms.untracked.double(0.1),
0082
0083 # ------- max EmFraction of tag jet allowed (does not affect isotracks)
0084   maxTagJetEmFrac = cms.untracked.double(1.0),
0085
0086 # ------- max abs(eta) of tag jet allowed (does not affect isotracks)
0087   maxTagJetAbsEta = cms.untracked.double(1.0),
0088
0089 # ------- min Et of tag jet allowed (does not affect isotracks)
0090   minTagJetEt = cms.untracked.double(30.0),
0091
0092 # ------- min abs(eta) of probe jet allowed (does not affect isotracks)
0093   minProbeJetAbsEta = cms.untracked.double(1.4),
0094
0095 # ------- file containing the phi symmetry corrections (to apply on the fly in isoTrack calibration)
0096 # ------- USER PROVIDED!
0097   phiSymCorFileName = cms.untracked.string("phiSymCor.txt"),
0098
0099 # ------- Flag to read phi symmetry corrections and apply on the fly
0100 # ------- Set to false if already applied in Analyzer
0101   applyPhiSymCorFlag = cms.untracked.bool(False),
0102
0103 # ------- output file with correction factors
0104   outputCorCoefFileName = cms.untracked.string("calibConst_IsoTrk.txt"),
0105
0106 # ------- output file with some hisograms
0107   histoFileName = cms.untracked.string("histoFile.root")
0108
0109 )
0110
0111
0112 process.p = cms.Path(process.hcalCalib)