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