Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2025-04-17 02:41:47

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 before CMSSW_3_1_x
0019 ## process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
0020 # or Get DB parameters 
0021 # process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
0022 # process.GlobalTag.globaltag = "GR09_31X_V2H::All"
0023 process.load("CondCore.DBCommon.CondDBCommon_cfi")
0024 process.CondDBCommon.connect = 'oracle://cms_orcon_prod/CMS_COND_31X_ECAL' ##P5 stuff
0025 #process.CondDBCommon.connect = 'oracle://cms_orcoff_prod/CMS_COND_31X_ECAL' ##lxplus
0026 process.CondDBCommon.DBParameters.authenticationPath = '/nfshome0/popcondev/conddb' ###P5 stuff
0027 #process.CondDBCommon.DBParameters.authenticationPath = '/afs/cern.ch/cms/DB/conddb' ##lxplus
0028 
0029 process.PoolDBESSource = cms.ESSource("PoolDBESSource",
0030                              process.CondDBCommon,
0031                               toGet = cms.VPSet(
0032                    cms.PSet(
0033                        record = cms.string('EcalPedestalsRcd'),
0034                        tag = cms.string('EcalPedestals_mc')
0035                    ),
0036                    cms.PSet(
0037                        record = cms.string('EcalADCToGeVConstantRcd'),
0038                        tag = cms.string('EcalADCToGeVConstant_2010_V2_Bon_mc')
0039                    ),
0040                    cms.PSet(
0041                        record = cms.string('EcalIntercalibConstantsRcd'),
0042                        tag = cms.string('EcalIntercalibConstants_2010_V3_Bon_startup_mc')
0043                    ),
0044                    cms.PSet(
0045                        record = cms.string('EcalGainRatiosRcd'),
0046                        tag = cms.string('EcalGainRatios_mc')
0047                    ),
0048                    cms.PSet(
0049                        record = cms.string('EcalMappingElectronicsRcd'),
0050                        tag = cms.string('EcalMappingElectronics_EEMap')
0051                    )
0052     )
0053 )
0054 
0055 
0056 #########################
0057 process.source = cms.Source("EmptySource",
0058        firstRun = cms.untracked.uint32(100000000) ### need to use latest run to pick-up update values from DB 
0059 )
0060 
0061 
0062 process.maxEvents = cms.untracked.PSet(
0063     input = cms.untracked.int32(1)
0064 )
0065 
0066 ##process.TPGParamProducer = cms.EDFilter("EcalTPGParamBuilder",
0067 process.TPGParamProducer = cms.EDAnalyzer("EcalTPGParamBuilder",
0068 
0069     #### inputs/ouputs control ####
0070     writeToDB  = cms.bool(False),
0071     allowDBEE  = cms.bool(True),
0072 
0073     DBsid   = cms.string('cms_omds_lb'), ## real DB 
0074     ##DBsid   = cms.string('int2r'), ## test DB
0075     ## P5 online DB
0076     DBuser  = cms.string('cms_ecal_conf'),
0077     DBpass  = cms.string('0r4cms_3c4lc0nf'),
0078     ## test DB
0079     ##DBuser  = cms.string('cms_ecal_conf_test'),
0080     ##DBpass  = cms.string('0r4cms_3c4l'),
0081     DBport  = cms.uint32(10121),
0082 
0083     TPGWritePed = cms.uint32(1),
0084     TPGWriteLin = cms.uint32(1),
0085     TPGWriteSli = cms.uint32(1),
0086     TPGWriteWei = cms.uint32(1),
0087     TPGWriteLut = cms.uint32(1),
0088     TPGWriteFgr = cms.uint32(1),
0089     TPGWriteSpi = cms.uint32(1),
0090     TPGWriteDel = cms.uint32(1),
0091     TPGWriteBxt = cms.uint32(0),
0092     TPGWriteBtt = cms.uint32(0), #do not change
0093 
0094     writeToFiles = cms.bool(True),
0095     outFile = cms.string('TPG_beamv5_MC_startup.txt'), #modif
0096    #### TPG config tag and version (if not given it will be automatically given ) ####
0097     TPGtag = cms.string('BEAMV5MCSTARTUP'),
0098     TPGversion = cms.uint32(1),
0099                                         
0100    #### TPG calculation parameters ####
0101     useTransverseEnergy = cms.bool(True),    ## true when TPG computes transverse energy, false for energy
0102     Et_sat_EB = cms.double(128.0),            ## Saturation value (in GeV) of the TPG before the compressed-LUT (rem: with 35.84 the TPG_LSB = crystal_LSB)
0103     Et_sat_EE = cms.double(128.0),            ## Saturation value (in GeV) of the TPG before the compressed-LUT (rem: with 35.84 the TPG_LSB = crystal_LSB)
0104 
0105     sliding = cms.uint32(0),                 ## Parameter used for the FE data format, should'nt be changed
0106 
0107     weight_timeShift = cms.double(0.),       ## weights are computed shifting the timing of the shape by this amount in ns: val>0 => shape shifted to the right
0108     weight_sampleMax = cms.uint32(3),        ## position of the maximum among the 5 samples used by the TPG amplitude filter
0109     weight_unbias_recovery = cms.bool(True), ## true if weights after int conversion are forced to have sum=0. Pb, in that case it can't have sum f*w = 1
0110 
0111     forcedPedestalValue = cms.int32(-3),     ## use this value instead of getting it from DB or MC
0112                                              ## -1: means use value from DB or MC.
0113                                              ## -2: ped12 = 0 used to cope with FENIX bug
0114                                              ## -3: used with sFGVB: baseline subtracted is pedestal-offset*sin(theta)/G with G=mult*2^-(shift+2) 
0115     pedestal_offset =  cms.uint32(150),      ## pedestal offset used with option forcedPedestalValue = -3
0116 
0117     useInterCalibration = cms.bool(True),    ## use or not values from DB. If not, 1 is assumed
0118 
0119     timing_delays_EB = cms.string('Delays_EB.txt'), # timing delays for latency EB / TT 
0120     timing_delays_EE = cms.string('Delays_EE.txt'), # timing delays for latency EE / strip                                         
0121     timing_phases_EB = cms.string('Phases_EB.txt'), # TCC phase setting for EB / TT
0122     timing_phases_EE = cms.string('Phases_EE.txt'), # TCC phase setting for EE / strip                                                                                 
0123 
0124     SFGVB_Threshold = cms.uint32(9),              ## (adc) SFGVB threshold in FE
0125     SFGVB_lut = cms.uint32(0xfffefee8),           ## SFGVB LUT in FE                                
0126     SFGVB_SpikeKillingThreshold = cms.int32(8),   ## (GeV) Spike killing threshold applied in TPG ET in TCC
0127                                         
0128     forceEtaSlice = cms.bool(False),         ## when true, same linearization coeff for all crystals belonging to a given eta slice (tower)
0129 
0130     LUT_option = cms.string('Linear'),       ## compressed LUT option can be: "Identity", "Linear", "EcalResolution"
0131     LUT_threshold_EB = cms.double(0.250),    ## All Trigger Primitives <= threshold (in GeV) will be set to 0 
0132     LUT_threshold_EE = cms.double(0.250),    ## All Trigger Primitives <= threshold (in GeV) will be set to 0 
0133     LUT_stochastic_EB = cms.double(0.03),    ## Stochastic term of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
0134     LUT_noise_EB = cms.double(0.2),          ## noise term (GeV) of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
0135     LUT_constant_EB = cms.double(0.005),     ## constant term of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
0136     LUT_stochastic_EE = cms.double(0.03),    ## Stochastic term of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
0137     LUT_noise_EE = cms.double(0.2),          ## noise term (GeV) of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
0138     LUT_constant_EE = cms.double(0.005),     ## constant term of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
0139 
0140     TTF_lowThreshold_EB = cms.double(1.0),   ## EB Trigger Tower Flag low threshold in GeV
0141     TTF_highThreshold_EB = cms.double(2.0),  ## EB Trigger Tower Flag high threshold in GeV
0142     TTF_lowThreshold_EE = cms.double(1.0),   ## EE Trigger Tower Flag low threshold in GeV
0143     TTF_highThreshold_EE = cms.double(2.0),  ## EE Trigger Tower Flag high threshold in GeV
0144 
0145     FG_lowThreshold_EB = cms.double(3.9),      ## EB Fine Grain Et low threshold in GeV
0146     FG_highThreshold_EB = cms.double(3.9),     ## EB Fine Grain Et high threshold in GeV
0147     FG_lowRatio_EB = cms.double(0.9),          ## EB Fine Grain low-ratio
0148     FG_highRatio_EB = cms.double(0.9),         ## EB Fine Grain high-ratio
0149     FG_lut_EB = cms.uint32(0x08),              ## EB Fine Grain Look-up table. Put something != 0 if you really know what you do!
0150     FG_Threshold_EE = cms.double(18.75),       ## EE Fine threshold in GeV
0151     FG_lut_strip_EE = cms.uint32(0xfffefee8),  ## EE Fine Grain strip Look-up table
0152     FG_lut_tower_EE = cms.uint32(0)            ## EE Fine Grain tower Look-up table
0153 )
0154 
0155 process.p = cms.Path(process.TPGParamProducer)
0156