Configuration/python/postLS1Customs.py

0001 from __future__ import print_function
0002 import FWCore.ParameterSet.Config as cms
0003
0004 def customisePostLS1_Common(process):
0005
0006     # deal with L1 Emulation separately
0007     from L1Trigger.L1TCommon.customsPostLS1 import customiseSimL1EmulatorForStage1
0008     process = customiseSimL1EmulatorForStage1(process)
0009
0010     # deal with CSC separately
0011     from SLHCUpgradeSimulations.Configuration.muonCustoms import customise_csc_PostLS1
0012     process = customise_csc_PostLS1(process)
0013
0014     # deal with FastSim separately
0015     from SLHCUpgradeSimulations.Configuration.fastSimCustoms import customise_fastSimPostLS1
0016     process = customise_fastSimPostLS1(process)
0017
0018     # all the rest
0019     if hasattr(process,'g4SimHits'):
0020         process = customise_Sim(process)
0021     if hasattr(process,'DigiToRaw'):
0022         process = customise_DigiToRaw(process)
0023     if hasattr(process,'RawToDigi'):
0024         process = customise_RawToDigi(process)
0025     if hasattr(process,'reconstruction'):
0026         process = customise_Reco(process)
0027     if hasattr(process,'digitisation_step') or ( hasattr(process,'mix') and hasattr(process.mix,'digitizers')):
0028         process = customise_Digi_Common(process)
0029     if hasattr(process,'L1simulation_step'):
0030         process = customise_L1Emulator(process)
0031     if hasattr(process,'dqmoffline_step'):
0032         process = customise_DQM(process)
0033     if hasattr(process,'dqmHarvesting'):
0034         process = customise_harvesting(process)
0035     if hasattr(process,'validation_step'):
0036         process = customise_Validation(process)
0037     if hasattr(process,'datamixing_step'):
0038         process = customise_DataMix(process)
0039
0040     return process
0041
0042
0043 def customisePostLS1(process,displayDeprecationWarning=True):
0044     if displayDeprecationWarning :
0045         print("""
0046         #
0047         # -- Warning! You are using a deprecated customisation function. --
0048         #
0049         # It will probably run fine, but the customisations you are getting may be out of date.
0050         # You should update your configuration file by
0051         #   If using cmsDriver:
0052         #       1) remove the "--customise SLHCUpgradeSimulations/Configuration/postLS1Customs.customisePostLS1" option
0053         #       2) add the option "--era Run2_25ns"
0054         #   If using a pre-made configuration file:
0055         #       1) remove or comment out the "process = customisePostLS1(process)" line.
0056         #       2) add "from Configuration.Eras.Era_Run2_25ns_cff import Run2_25ns" to the TOP of the config file (above
0057         #          the process declaration).
0058         #       3) add "Run2_25ns" as a parameter to the process object, e.g. "process = cms.Process('HLT',Run2_25ns)"
0059         #
0060         # There is more information at https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideCmsDriverEras
0061         #
0062         """)
0063     # common customisation
0064     process = customisePostLS1_Common(process)
0065
0066     # 25ns specific customisation
0067     if hasattr(process,'digitisation_step'):
0068         process = customise_Digi_25ns(process)
0069     if hasattr(process,'dqmoffline_step'):
0070         process = customise_DQM_25ns(process)
0071
0072     return process
0073
0074
0075 def customisePostLS1_lowPU(process):
0076
0077     # common customisations
0078     process = customisePostLS1_Common(process)
0079
0080     # 50ns specific customisation
0081     if hasattr(process,'digitisation_step'):
0082         process = customise_Digi_50ns(process)
0083
0084     return process
0085
0086
0087 def customisePostLS1_50ns(process,displayDeprecationWarning=True):
0088     if displayDeprecationWarning :
0089         print("""
0090         #
0091         # -- Warning! You are using a deprecated customisation function. --
0092         #
0093         # It will probably run fine, but the customisations you are getting may be out of date.
0094         # You should update your configuration file by
0095         #   If using cmsDriver:
0096         #       1) remove the "--customise SLHCUpgradeSimulations/Configuration/postLS1Customs.customisePostLS1_50ns" option
0097         #       2) add the option "--era Run2_50ns"
0098         #   If using a pre-made configuration file:
0099         #       1) remove or comment out the "process = customisePostLS1_50ns(process)" line.
0100         #       2) add "from Configuration.Eras.Era_Run2_50ns_cff import Run2_50ns" to the TOP of the config file (above
0101         #          the process declaration).
0102         #       3) add "Run2_50ns" as a parameter to the process object, e.g. "process = cms.Process('HLT',Run2_50ns)"
0103         #
0104         # There is more information at https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideCmsDriverEras
0105         #
0106         """)
0107
0108     # common customisations
0109     process = customisePostLS1_Common(process)
0110
0111     # 50ns specific customisation
0112     if hasattr(process,'digitisation_step'):
0113         process = customise_Digi_50ns(process)
0114
0115     return process
0116
0117
0118 def customisePostLS1_HI(process,displayDeprecationWarning=True):
0119     if displayDeprecationWarning :
0120         print("""
0121         #
0122         # -- Warning! You are using a deprecated customisation function. --
0123         #
0124         # It will probably run fine, but the customisations you are getting may be out of date.
0125         # You should update your configuration file by
0126         #   If using cmsDriver:
0127         #       1) remove the "--customise SLHCUpgradeSimulations/Configuration/postLS1Customs.customisePostLS1_HI" option
0128         #       2) add the option "--era Run2_HI"
0129         #   If using a pre-made configuration file:
0130         #       1) remove or comment out the "process = customisePostLS1_HI(process)" line.
0131         #       2) add "from Configuration.Eras.Era_Run2_HI_cff import Run2_HI" to the TOP of the config file (above
0132         #          the process declaration).
0133         #       3) add "Run2_HI" as a parameter to the process object, e.g. "process = cms.Process('HLT',Run2_HI)"
0134         #
0135         # There is more information at https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideCmsDriverEras
0136         #
0137         """)
0138
0139     # common customisation
0140     process = customisePostLS1_Common(process)
0141
0142     # HI Specific additional customizations:
0143     from L1Trigger.L1TCommon.customsPostLS1 import customiseSimL1EmulatorForPostLS1_Additional_HI
0144     process = customiseSimL1EmulatorForPostLS1_Additional_HI(process)
0145
0146     # HI L1Menu:
0147     #from L1Trigger.Configuration.customise_overwriteL1Menu import L1Menu_CollisionsHeavyIons2015_v0
0148     #process = L1Menu_CollisionsHeavyIons2015_v0(process)
0149
0150     return process
0151
0152 def digiEventContent(process):
0153     #extend the event content
0154
0155     alist=['RAWDEBUG','FEVTDEBUG','FEVTDEBUGHLT','GENRAW','RAWSIMHLT','FEVT']
0156     for a in alist:
0157         b = a + 'output'
0158         if hasattr(process,b):
0159             getattr(process,b).outputCommands.append('keep *_simMuonCSCDigis_*_*')
0160             getattr(process,b).outputCommands.append('keep *_simMuonRPCDigis_*_*')
0161             getattr(process,b).outputCommands.append('keep *_simHcalUnsuppressedDigis_*_*')
0162
0163     return process
0164
0165
0166 def customise_DQM(process):
0167     #process.dqmoffline_step.remove(process.jetMETAnalyzer)
0168     # Turn off flag of gangedME11a
0169     process.l1tCsctf.gangedME11a = cms.untracked.bool(False)
0170     # Turn off "low bias voltage" region in HCAL noise filters
0171     if hasattr(process,'HBHENoiseFilterResultProducer'):
0172         process.HBHENoiseFilterResultProducer.IgnoreTS4TS5ifJetInLowBVRegion = cms.bool(False)
0173     return process
0174
0175
0176 def customise_DQM_25ns(process):
0177     #Empty place-holder
0178     # Switch to Stage1 Digi
0179     process.load("L1Trigger.L1TCommon.l1tRawToDigi_cfi")
0180     process.load("L1Trigger.L1TCommon.caloStage1LegacyFormatDigis_cfi")
0181     if hasattr(process, 'RawToDigi'):
0182         process.L1RawToDigiSeq = cms.Sequence(process.gctDigis+process.caloStage1Digis+process.caloStage1LegacyFormatDigis)
0183         process.RawToDigi.replace(process.gctDigis, process.L1RawToDigiSeq)
0184
0185     process.load('DQMOffline.L1Trigger.L1TriggerDqmOffline_cff')
0186     if hasattr(process, 'l1tMonitorStage1Online'):
0187         process.l1tRct.rctSource = 'caloStage1Digis'
0188         process.l1tRctfromRCT.rctSource = 'rctDigis'
0189         process.l1tPUM.regionSource = cms.InputTag("rctDigis")
0190         process.l1tStage1Layer2.stage1_layer2_ = cms.bool(True)
0191         process.l1tStage1Layer2.gctCentralJetsSource = cms.InputTag("caloStage1LegacyFormatDigis","cenJets")
0192         process.l1tStage1Layer2.gctForwardJetsSource = cms.InputTag("caloStage1LegacyFormatDigis","forJets")
0193         process.l1tStage1Layer2.gctTauJetsSource = cms.InputTag("caloStage1LegacyFormatDigis","tauJets")
0194         process.l1tStage1Layer2.gctIsoTauJetsSource = cms.InputTag("caloStage1LegacyFormatDigis","isoTauJets")
0195         process.l1tStage1Layer2.gctEnergySumsSource = cms.InputTag("caloStage1LegacyFormatDigis")
0196         process.l1tStage1Layer2.gctIsoEmSource = cms.InputTag("caloStage1LegacyFormatDigis","isoEm")
0197         process.l1tStage1Layer2.gctNonIsoEmSource = cms.InputTag("caloStage1LegacyFormatDigis","nonIsoEm")
0198         process.dqmL1ExtraParticlesStage1.etTotalSource = 'caloStage1LegacyFormatDigis'
0199         process.dqmL1ExtraParticlesStage1.nonIsolatedEmSource = 'caloStage1LegacyFormatDigis:nonIsoEm'
0200         process.dqmL1ExtraParticlesStage1.etMissSource = 'caloStage1LegacyFormatDigis'
0201         process.dqmL1ExtraParticlesStage1.htMissSource = 'caloStage1LegacyFormatDigis'
0202         process.dqmL1ExtraParticlesStage1.forwardJetSource = 'caloStage1LegacyFormatDigis:forJets'
0203         process.dqmL1ExtraParticlesStage1.centralJetSource = 'caloStage1LegacyFormatDigis:cenJets'
0204         process.dqmL1ExtraParticlesStage1.tauJetSource = 'caloStage1LegacyFormatDigis:tauJets'
0205         process.dqmL1ExtraParticlesStage1.isolatedEmSource = 'caloStage1LegacyFormatDigis:isoEm'
0206         process.dqmL1ExtraParticlesStage1.etHadSource = 'caloStage1LegacyFormatDigis'
0207         process.dqmL1ExtraParticlesStage1.hfRingEtSumsSource = 'caloStage1LegacyFormatDigis'
0208         process.dqmL1ExtraParticlesStage1.hfRingBitCountsSource = 'caloStage1LegacyFormatDigis'
0209         process.l1ExtraDQMStage1.stage1_layer2_ = cms.bool(True)
0210         process.l1ExtraDQMStage1.L1ExtraIsoTauJetSource_ = cms.InputTag("dqmL1ExtraParticlesStage1", "IsoTau")
0211
0212     if hasattr(process, 'l1Stage1HwValEmulatorMonitor'):
0213         process.l1TdeRCT.rctSourceData = 'caloStage1Digis'
0214         process.l1TdeRCTfromRCT.rctSourceData = 'rctDigis'
0215         process.l1compareforstage1.GCTsourceData = cms.InputTag("caloStage1LegacyFormatDigis")
0216         process.l1compareforstage1.GCTsourceEmul = cms.InputTag("valCaloStage1LegacyFormatDigis")
0217         process.l1compareforstage1.stage1_layer2_ = cms.bool(True)
0218         process.valStage1GtDigis.GctInputTag = 'caloStage1LegacyFormatDigis'
0219
0220     return process
0221
0222
0223 def customise_Validation(process):
0224     #process.validation_step.remove(process.PixelTrackingRecHitsValid)
0225     # We don't run the HLT
0226     #process.validation_step.remove(process.HLTSusyExoVal)
0227     #process.validation_step.remove(process.hltHiggsValidator)
0228     return process
0229
0230
0231 def customise_Sim(process):
0232     # enable 2015 HF shower library
0233     process.g4SimHits.HFShowerLibrary.FileName = 'SimG4CMS/Calo/data/HFShowerLibrary_npmt_noatt_eta4_16en_v4.root'
0234     return process
0235
0236 def customise_Digi_Common(process):
0237     process = digiEventContent(process)
0238     if hasattr(process,'mix') and hasattr(process.mix,'digitizers'):
0239         if hasattr(process.mix.digitizers,'hcal') and hasattr(process.mix.digitizers.hcal,'ho'):
0240             process.mix.digitizers.hcal.ho.photoelectronsToAnalog = cms.vdouble([4.0]*16)
0241             process.mix.digitizers.hcal.ho.siPMCode = cms.int32(1)
0242             process.mix.digitizers.hcal.ho.pixels = cms.int32(2500)
0243             process.mix.digitizers.hcal.ho.doSiPMSmearing = cms.bool(False)
0244         if hasattr(process.mix.digitizers,'hcal') and hasattr(process.mix.digitizers.hcal,'hf1'):
0245             process.mix.digitizers.hcal.hf1.samplingFactor = cms.double(0.60)
0246         if hasattr(process.mix.digitizers,'hcal') and hasattr(process.mix.digitizers.hcal,'hf2'):
0247             process.mix.digitizers.hcal.hf2.samplingFactor = cms.double(0.60)
0248     return process
0249
0250 def customise_DataMix(process):
0251     if hasattr(process,'mixData'):
0252         if hasattr(process.mixData,'ho'):
0253             process.mixData.ho.photoelectronsToAnalog = cms.vdouble([4.0]*16)
0254             process.mixData.ho.siPMCode = cms.int32(1)
0255             process.mixData.ho.pixels = cms.int32(2500)
0256             process.mixData.ho.doSiPMSmearing = cms.bool(False)
0257         if hasattr(process.mixData,'hf1'):
0258             process.mixData.hf1.samplingFactor = cms.double(0.60)
0259         if hasattr(process.mixData,'hf2'):
0260             process.mixData.hf2.samplingFactor = cms.double(0.60)
0261     return process
0262
0263
0264 def customise_Digi_50ns(process):
0265     if hasattr(process,'mix') and hasattr(process.mix,'digitizers'):
0266         if hasattr(process.mix.digitizers,'pixel'):
0267             # pixel dynamic inefficency - 13TeV - 50ns case
0268             process.mix.digitizers.pixel.theInstLumiScaleFactor = cms.double(246.4)
0269             process.mix.digitizers.pixel.theLadderEfficiency_BPix1 = cms.vdouble(
0270                 0.979259,
0271                 0.976677,
0272                 0.979259,
0273                 0.976677,
0274                 0.979259,
0275                 0.976677,
0276                 0.979259,
0277                 0.976677,
0278                 0.979259,
0279                 0.976677,
0280                 0.979259,
0281                 0.976677,
0282                 0.979259,
0283                 0.976677,
0284                 0.979259,
0285                 0.976677,
0286                 0.979259,
0287                 0.976677,
0288                 0.979259,
0289                 0.976677,
0290             )
0291             process.mix.digitizers.pixel.theLadderEfficiency_BPix2 = cms.vdouble(
0292                 0.994321,
0293                 0.993944,
0294                 0.994321,
0295                 0.993944,
0296                 0.994321,
0297                 0.993944,
0298                 0.994321,
0299                 0.993944,
0300                 0.994321,
0301                 0.993944,
0302                 0.994321,
0303                 0.993944,
0304                 0.994321,
0305                 0.993944,
0306                 0.994321,
0307                 0.993944,
0308                 0.994321,
0309                 0.993944,
0310                 0.994321,
0311                 0.993944,
0312                 0.994321,
0313                 0.993944,
0314                 0.994321,
0315                 0.993944,
0316                 0.994321,
0317                 0.993944,
0318                 0.994321,
0319                 0.993944,
0320                 0.994321,
0321                 0.993944,
0322                 0.994321,
0323                 0.993944,
0324             )
0325             process.mix.digitizers.pixel.theLadderEfficiency_BPix3 = cms.vdouble(
0326                 0.996787,
0327                 0.996945,
0328                 0.996787,
0329                 0.996945,
0330                 0.996787,
0331                 0.996945,
0332                 0.996787,
0333                 0.996945,
0334                 0.996787,
0335                 0.996945,
0336                 0.996787,
0337                 0.996945,
0338                 0.996787,
0339                 0.996945,
0340                 0.996787,
0341                 0.996945,
0342                 0.996787,
0343                 0.996945,
0344                 0.996787,
0345                 0.996945,
0346                 0.996787,
0347                 0.996945,
0348                 0.996787,
0349                 0.996945,
0350                 0.996787,
0351                 0.996945,
0352                 0.996787,
0353                 0.996945,
0354                 0.996787,
0355                 0.996945,
0356                 0.996787,
0357                 0.996945,
0358                 0.996787,
0359                 0.996945,
0360                 0.996787,
0361                 0.996945,
0362                 0.996787,
0363                 0.996945,
0364                 0.996787,
0365                 0.996945,
0366                 0.996787,
0367                 0.996945,
0368                 0.996787,
0369                 0.996945,
0370             )
0371     if hasattr(process,'mixData'):
0372             # pixel dynamic inefficency - 13TeV - 50ns case
0373         process.mixData.theInstLumiScaleFactor = cms.double(246.4)
0374         process.mixData.theLadderEfficiency_BPix1 = cms.vdouble(
0375                 0.979259,
0376                 0.976677,
0377                 0.979259,
0378                 0.976677,
0379                 0.979259,
0380                 0.976677,
0381                 0.979259,
0382                 0.976677,
0383                 0.979259,
0384                 0.976677,
0385                 0.979259,
0386                 0.976677,
0387                 0.979259,
0388                 0.976677,
0389                 0.979259,
0390                 0.976677,
0391                 0.979259,
0392                 0.976677,
0393                 0.979259,
0394                 0.976677,
0395             )
0396         process.mixData.theLadderEfficiency_BPix2 = cms.vdouble(
0397                 0.994321,
0398                 0.993944,
0399                 0.994321,
0400                 0.993944,
0401                 0.994321,
0402                 0.993944,
0403                 0.994321,
0404                 0.993944,
0405                 0.994321,
0406                 0.993944,
0407                 0.994321,
0408                 0.993944,
0409                 0.994321,
0410                 0.993944,
0411                 0.994321,
0412                 0.993944,
0413                 0.994321,
0414                 0.993944,
0415                 0.994321,
0416                 0.993944,
0417                 0.994321,
0418                 0.993944,
0419                 0.994321,
0420                 0.993944,
0421                 0.994321,
0422                 0.993944,
0423                 0.994321,
0424                 0.993944,
0425                 0.994321,
0426                 0.993944,
0427                 0.994321,
0428                 0.993944,
0429             )
0430         process.mixData.theLadderEfficiency_BPix3 = cms.vdouble(
0431                 0.996787,
0432                 0.996945,
0433                 0.996787,
0434                 0.996945,
0435                 0.996787,
0436                 0.996945,
0437                 0.996787,
0438                 0.996945,
0439                 0.996787,
0440                 0.996945,
0441                 0.996787,
0442                 0.996945,
0443                 0.996787,
0444                 0.996945,
0445                 0.996787,
0446                 0.996945,
0447                 0.996787,
0448                 0.996945,
0449                 0.996787,
0450                 0.996945,
0451                 0.996787,
0452                 0.996945,
0453                 0.996787,
0454                 0.996945,
0455                 0.996787,
0456                 0.996945,
0457                 0.996787,
0458                 0.996945,
0459                 0.996787,
0460                 0.996945,
0461                 0.996787,
0462                 0.996945,
0463                 0.996787,
0464                 0.996945,
0465                 0.996787,
0466                 0.996945,
0467                 0.996787,
0468                 0.996945,
0469                 0.996787,
0470                 0.996945,
0471                 0.996787,
0472                 0.996945,
0473                 0.996787,
0474                 0.996945,
0475             )
0476     return process
0477
0478
0479 def customise_Digi_25ns(process):
0480     if hasattr(process,'mix') and hasattr(process.mix,'digitizers'):
0481         if hasattr(process.mix.digitizers,'pixel'):
0482             # pixel dynamic inefficency - 13TeV - 25ns case
0483             process.mix.digitizers.pixel.theInstLumiScaleFactor = cms.double(364)
0484             process.mix.digitizers.pixel.theLadderEfficiency_BPix1 = cms.vdouble(
0485                 1,
0486                 1,
0487                 1,
0488                 1,
0489                 1,
0490                 1,
0491                 1,
0492                 1,
0493                 1,
0494                 1,
0495                 1,
0496                 1,
0497                 1,
0498                 1,
0499                 1,
0500                 1,
0501                 1,
0502                 1,
0503                 1,
0504                 1,
0505             )
0506             process.mix.digitizers.pixel.theLadderEfficiency_BPix2 = cms.vdouble(
0507                 1,
0508                 1,
0509                 1,
0510                 1,
0511                 1,
0512                 1,
0513                 1,
0514                 1,
0515                 1,
0516                 1,
0517                 1,
0518                 1,
0519                 1,
0520                 1,
0521                 1,
0522                 1,
0523                 1,
0524                 1,
0525                 1,
0526                 1,
0527                 1,
0528                 1,
0529                 1,
0530                 1,
0531                 1,
0532                 1,
0533                 1,
0534                 1,
0535                 1,
0536                 1,
0537                 1,
0538                 1,
0539             )
0540             process.mix.digitizers.pixel.theLadderEfficiency_BPix3 = cms.vdouble(
0541                 1,
0542                 1,
0543                 1,
0544                 1,
0545                 1,
0546                 1,
0547                 1,
0548                 1,
0549                 1,
0550                 1,
0551                 1,
0552                 1,
0553                 1,
0554                 1,
0555                 1,
0556                 1,
0557                 1,
0558                 1,
0559                 1,
0560                 1,
0561                 1,
0562                 1,
0563                 1,
0564                 1,
0565                 1,
0566                 1,
0567                 1,
0568                 1,
0569                 1,
0570                 1,
0571                 1,
0572                 1,
0573                 1,
0574                 1,
0575                 1,
0576                 1,
0577                 1,
0578                 1,
0579                 1,
0580                 1,
0581                 1,
0582                 1,
0583                 1,
0584                 1,
0585             )
0586             process.mix.digitizers.pixel.theModuleEfficiency_BPix1 = cms.vdouble(
0587                 1,
0588                 1,
0589                 1,
0590                 1,
0591             )
0592             process.mix.digitizers.pixel.theModuleEfficiency_BPix2 = cms.vdouble(
0593                 1,
0594                 1,
0595                 1,
0596                 1,
0597                 )
0598             process.mix.digitizers.pixel.theModuleEfficiency_BPix3 = cms.vdouble(
0599                 1,
0600                 1,
0601                 1,
0602                 1,
0603             )
0604             process.mix.digitizers.pixel.thePUEfficiency_BPix1 = cms.vdouble(
0605                  1.00023,
0606                 -3.18350e-06,
0607                  5.08503e-10,
0608                 -6.79785e-14,
0609             )
0610             process.mix.digitizers.pixel.thePUEfficiency_BPix2 = cms.vdouble(
0611                  9.99974e-01,
0612                 -8.91313e-07,
0613                  5.29196e-12,
0614                 -2.28725e-15,
0615             )
0616             process.mix.digitizers.pixel.thePUEfficiency_BPix3 = cms.vdouble(
0617                  1.00005,
0618                 -6.59249e-07,
0619                  2.75277e-11,
0620                 -1.62683e-15,
0621             )
0622     if hasattr(process,'mixData'):
0623             # pixel dynamic inefficency - 13TeV - 25ns case
0624         process.mixData.theInstLumiScaleFactor = cms.double(364)
0625         process.mixData.theLadderEfficiency_BPix1 = cms.vdouble(
0626                 1,
0627                 1,
0628                 1,
0629                 1,
0630                 1,
0631                 1,
0632                 1,
0633                 1,
0634                 1,
0635                 1,
0636                 1,
0637                 1,
0638                 1,
0639                 1,
0640                 1,
0641                 1,
0642                 1,
0643                 1,
0644                 1,
0645                 1,
0646             )
0647         process.mixData.theLadderEfficiency_BPix2 = cms.vdouble(
0648                 1,
0649                 1,
0650                 1,
0651                 1,
0652                 1,
0653                 1,
0654                 1,
0655                 1,
0656                 1,
0657                 1,
0658                 1,
0659                 1,
0660                 1,
0661                 1,
0662                 1,
0663                 1,
0664                 1,
0665                 1,
0666                 1,
0667                 1,
0668                 1,
0669                 1,
0670                 1,
0671                 1,
0672                 1,
0673                 1,
0674                 1,
0675                 1,
0676                 1,
0677                 1,
0678                 1,
0679                 1,
0680             )
0681         process.mixData.theLadderEfficiency_BPix3 = cms.vdouble(
0682                 1,
0683                 1,
0684                 1,
0685                 1,
0686                 1,
0687                 1,
0688                 1,
0689                 1,
0690                 1,
0691                 1,
0692                 1,
0693                 1,
0694                 1,
0695                 1,
0696                 1,
0697                 1,
0698                 1,
0699                 1,
0700                 1,
0701                 1,
0702                 1,
0703                 1,
0704                 1,
0705                 1,
0706                 1,
0707                 1,
0708                 1,
0709                 1,
0710                 1,
0711                 1,
0712                 1,
0713                 1,
0714                 1,
0715                 1,
0716                 1,
0717                 1,
0718                 1,
0719                 1,
0720                 1,
0721                 1,
0722                 1,
0723                 1,
0724                 1,
0725                 1,
0726             )
0727         process.mixData.theModuleEfficiency_BPix1 = cms.vdouble(
0728                 1,
0729                 1,
0730                 1,
0731                 1,
0732             )
0733         process.mixData.theModuleEfficiency_BPix2 = cms.vdouble(
0734                 1,
0735                 1,
0736                 1,
0737                 1,
0738                 )
0739         process.mixData.theModuleEfficiency_BPix3 = cms.vdouble(
0740                 1,
0741                 1,
0742                 1,
0743                 1,
0744             )
0745         process.mixData.thePUEfficiency_BPix1 = cms.vdouble(
0746                  1.00023,
0747                 -3.18350e-06,
0748                  5.08503e-10,
0749                 -6.79785e-14,
0750             )
0751         process.mixData.thePUEfficiency_BPix2 = cms.vdouble(
0752                  9.99974e-01,
0753                 -8.91313e-07,
0754                  5.29196e-12,
0755                 -2.28725e-15,
0756             )
0757         process.mixData.thePUEfficiency_BPix3 = cms.vdouble(
0758                  1.00005,
0759                 -6.59249e-07,
0760                  2.75277e-11,
0761                 -1.62683e-15,
0762             )
0763     return process
0764
0765
0766 def customise_L1Emulator(process):
0767     return process
0768
0769
0770 def customise_RawToDigi(process):
0771     if hasattr(process,'gtEvmDigis'):
0772         process.RawToDigi.remove(process.gtEvmDigis)
0773     return process
0774
0775
0776 def customise_DigiToRaw(process):
0777     return process
0778
0779
0780 def customise_Reco(process):
0781     #lowering HO threshold with SiPM
0782     if hasattr(process,'particleFlowRecHitHO'):
0783         for prod in process.particleFlowRecHitHO.producers:
0784             prod.qualityTests = cms.VPSet(
0785                 cms.PSet(
0786                     name = cms.string("PFRecHitQTestThreshold"),
0787                     threshold = cms.double(0.05) # new threshold for SiPM HO
0788                     ),
0789                 cms.PSet(
0790                     name = cms.string("PFRecHitQTestHCALChannel"),
0791                     maxSeverities      = cms.vint32(11),
0792                     cleaningThresholds = cms.vdouble(0.0),
0793                     flags              = cms.vstring('Standard')
0794                     )
0795                 )
0796
0797     #Lower Thresholds also for Clusters!!!
0798
0799         for p in process.particleFlowClusterHO.seedFinder.thresholdsByDetector:
0800             p.seedingThreshold = cms.double(0.08)
0801
0802         for p in process.particleFlowClusterHO.initialClusteringStep.thresholdsByDetector:
0803             p.gatheringThreshold = cms.double(0.05)
0804
0805         for p in process.particleFlowClusterHO.pfClusterBuilder.recHitEnergyNorms:
0806             p.recHitEnergyNorm = cms.double(0.05)
0807
0808         process.particleFlowClusterHO.pfClusterBuilder.positionCalc.logWeightDenominator = cms.double(0.05)
0809         process.particleFlowClusterHO.pfClusterBuilder.allCellsPositionCalc.logWeightDenominator = cms.double(0.05)
0810
0811     # Muon reconstruction do not exclude bad chambers
0812     if hasattr(process, 'muonDetIdAssociator'):
0813         process.muonDetIdAssociator.includeBadChambers = cms.bool(True)
0814
0815     return process
0816
0817
0818 def customise_harvesting(process):
0819     #process.dqmHarvesting.remove(process.dataCertificationJetMET)
0820     #process.dqmHarvesting.remove(process.sipixelEDAClient)
0821     #process.dqmHarvesting.remove(process.sipixelCertification)
0822     return (process)
0823
0824
0825 def recoOutputCustoms(process):
0826
0827     alist=['AODSIM','RECOSIM','FEVTSIM','FEVTDEBUG','FEVTDEBUGHLT','RECODEBUG','RAWRECOSIMHLT','RAWRECODEBUGHLT']
0828     for a in alist:
0829         b = a + 'output'
0830         if hasattr(process,b):
0831             getattr(process,b).outputCommands.append('keep *_simMuonCSCDigis_*_*')
0832             getattr(process,b).outputCommands.append('keep *_simMuonRPCDigis_*_*')
0833             getattr(process,b).outputCommands.append('keep *_simHcalUnsuppressedDigis_*_*')
0834             getattr(process,b).outputCommands.append('keep *_rawDataCollector_*_*')
0835     return process