Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-11-18 00:40:21

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