EcalTPGTools/test/produceTPGParameters_Cosmics_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("ProdTPGParam")
0004
0005 # Calo geometry service model
0006 process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
0007 process.load("Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
0008 process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
0009
0010 # ecal mapping
0011 process.load("Geometry.EcalMapping.EcalMapping_cfi")
0012 process.eegeom = cms.ESSource("EmptyESSource",
0013     recordName = cms.string('EcalMappingRcd'),
0014     iovIsRunNotTime = cms.bool(True),
0015     firstValid = cms.vuint32(1)
0016 )
0017
0018 # Get hardcoded conditions the same used for standard digitization
0019 process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
0020
0021 #########################
0022 process.source = cms.Source("EmptySource")
0023
0024 process.maxEvents = cms.untracked.PSet(
0025     input = cms.untracked.int32(1)
0026 )
0027
0028 process.TPGParamProducer = cms.EDAnalyzer("EcalTPGParamBuilder",
0029
0030     #### inputs/ouputs control ####
0031     writeToDB  = cms.bool(False),
0032     allowDBEE  = cms.bool(False),
0033
0034     DBsid   = cms.string('ecalh4db'),
0035     DBuser  = cms.string('test09'),
0036     DBpass  = cms.string('oratest09'),
0037     DBport  = cms.uint32(1521),
0038
0039     writeToFiles = cms.bool(True),
0040     outFile = cms.string('TPG_cosmics.txt'),
0041    #### TPG config tag and version (if not given it will be automatically given ) ####
0042     TPGtag = cms.string('CRAFT'),
0043     TPGversion = cms.uint32(1),
0044
0045     #### TPG calculation parameters ####
0046
0047     #######################################################
0048     # Warning: for Cosmics, APD Gain = 200 (i.e. 50x4)    #
0049     #  ==> All numbers in GeV below must be divided by 4  #
0050     #      for the real energy correspondance             #
0051     #######################################################
0052
0053     useTransverseEnergy = cms.bool(False),   ## true when TPG computes transverse energy, false for energy
0054     Et_sat_EB = cms.double(35.84),          ## Saturation value (in GeV) of the TPG before the compressed-LUT (rem: with 35.84 the TPG_LSB = crystal_LSB)
0055     Et_sat_EE = cms.double(61.44),          ## Saturation value (in GeV) of the TPG before the compressed-LUT (rem: with 35.84 the TPG_LSB = crystal_LSB)
0056
0057     sliding = cms.uint32(0),                ## Parameter used for the FE data format, should'nt be changed
0058
0059     weight_sampleMax = cms.uint32(3),       ## position of the maximum among the 5 samples used by the TPG amplitude filter
0060
0061     forcedPedestalValue = cms.int32(160),   ## use this value instead of getting it from DB or MC (-1 means use DB or MC)
0062     forceEtaSlice = cms.bool(False),        ## when true, same linearization coeff for all crystals belonging to a given eta slice (tower). Implemented only for EB
0063
0064     LUT_option = cms.string('Identity'),    ## compressed LUT option can be: "Identity", "Linear", "EcalResolution"
0065     LUT_threshold_EB = cms.double(0.875),     ## All Trigger Primitives <= threshold (in GeV) will be set to 0
0066     LUT_threshold_EE = cms.double(0.875),     ## All Trigger Primitives <= threshold (in GeV) will be set to 0
0067     LUT_stochastic_EB = cms.double(0.03),   ## Stochastic term of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
0068     LUT_noise_EB = cms.double(0.2),         ## noise term (GeV) of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
0069     LUT_constant_EB = cms.double(0.005),    ## constant term of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
0070     LUT_stochastic_EE = cms.double(0.03),   ## Stochastic term of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
0071     LUT_noise_EE = cms.double(0.2),         ## noise term (GeV) of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
0072     LUT_constant_EE = cms.double(0.005),    ## constant term of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
0073
0074     TTF_lowThreshold_EB = cms.double(0.385),  ## EB Trigger Tower Flag low threshold in GeV
0075     TTF_highThreshold_EB = cms.double(0.385), ## EB Trigger Tower Flag high threshold in GeV
0076     TTF_lowThreshold_EE = cms.double(0.385),  ## EE Trigger Tower Flag low threshold in GeV
0077     TTF_highThreshold_EE = cms.double(0.385), ## EE Trigger Tower Flag high threshold in GeV
0078
0079     FG_lowThreshold_EB = cms.double(0.7),   ## EB Fine Grain Et low threshold in GeV
0080     FG_highThreshold_EB = cms.double(0.7),  ## EB Fine Grain Et high threshold in GeV
0081     FG_lowRatio_EB = cms.double(0.8),       ## EB Fine Grain low-ratio
0082     FG_highRatio_EB = cms.double(0.8),      ## EB Fine Grain high-ratio
0083     FG_lut_EB = cms.uint32(0),              ## EB Fine Grain Look-up table. Put something != 0 if you really know what you do!
0084     FG_Threshold_EE = cms.double(0.0),      ## EE Fine threshold in GeV
0085     FG_lut_strip_EE = cms.uint32(0),        ## EE Fine Grain strip Look-up table
0086     FG_lut_tower_EE = cms.uint32(0)         ## EE Fine Grain tower Look-up table
0087 )
0088
0089 process.p = cms.Path(process.TPGParamProducer)
0090