cms.PSet(
        depth=cms.vint32(1, 2, 3, 4),
        threshold = _thresholdsHB,
        detectorEnum = cms.int32(1)
        ),
    cms.PSet(
        depth=cms.vint32(1, 2, 3, 4, 5, 6, 7),
        threshold = _thresholdsHE,
        detectorEnum = cms.int32(2)
        )
    )

# offline 2018 -- uncollapsed
from Configuration.Eras.Modifier_run2_HCAL_2018_cff import run2_HCAL_2018
run2_HCAL_2018.toModify(particleFlowRecHitHBHE.producers[0].qualityTests[0], cuts = cuts2018)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(particleFlowRecHitHBHE.producers[0].qualityTests[0], cuts = cuts2018)

# offline 2018 -- collapsed
run2_HECollapse_2018 =  cms.Modifier()
run2_HECollapse_2018.toModify(particleFlowRecHitHBHE.producers[0].qualityTests[0], cuts = cuts2017)

# offline 2019
from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(particleFlowRecHitHBHE.producers[0].qualityTests[0], cuts = cuts2019)

# offline phase2 restore what has been studied in the TDR
from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
phase2_hcal.toModify(particleFlowRecHitHBHE.producers[0].qualityTests[0], cuts = cutsPhase2)
Exemple #2
0
                        recHitEnergyNorms=recHitEnergyNorms2018)
run2_HCAL_2018.toModify(
    particleFlowClusterHBHE.seedFinder,
    thresholdsByDetector=seedFinderThresholdsByDetector2018)
run2_HCAL_2018.toModify(
    particleFlowClusterHBHE.initialClusteringStep,
    thresholdsByDetector=initialClusteringStepThresholdsByDetector2018)
run2_HCAL_2018.toModify(
    particleFlowClusterHBHE.pfClusterBuilder.positionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2018)
run2_HCAL_2018.toModify(
    particleFlowClusterHBHE.pfClusterBuilder.allCellsPositionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2018)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder,
                      recHitEnergyNorms=recHitEnergyNorms2018)
run2_HE_2018.toModify(particleFlowClusterHBHE.seedFinder,
                      thresholdsByDetector=seedFinderThresholdsByDetector2018)
run2_HE_2018.toModify(
    particleFlowClusterHBHE.initialClusteringStep,
    thresholdsByDetector=initialClusteringStepThresholdsByDetector2018)
run2_HE_2018.toModify(
    particleFlowClusterHBHE.pfClusterBuilder.positionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2018)
run2_HE_2018.toModify(
    particleFlowClusterHBHE.pfClusterBuilder.allCellsPositionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2018)

# offline 2018 -- collapsed
run2_HECollapse_2018 = cms.Modifier()
run2_HECollapse_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder,
               readoutFrameSize = cms.int32(3), 
               binOfMaximum     = cms.int32(2)
              )
)

from Configuration.Eras.Modifier_run2_HB_2018_cff import run2_HB_2018
run2_HB_2018.toModify( hcalSimParameters,
    hb = dict(
               readoutFrameSize = cms.int32(8), 
               binOfMaximum     = cms.int32(4)
              )
)
from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify( hcalSimParameters,
    he = dict(
               readoutFrameSize = cms.int32(8), 
               binOfMaximum     = cms.int32(4)
              )
)


from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify( hcalSimParameters,
    hb = dict(
        doSiPMSmearing = cms.bool(True),
        sipmTau = cms.double(10.),
    )
)

_newFactors = cms.vdouble(
    210.55, 197.93, 186.12, 189.64, 189.63,
    189.96, 190.03, 190.11, 190.18, 190.25,
              )
    )

#######################

# offline 2018 -- uncollapsed
from Configuration.Eras.Modifier_run2_HCAL_2018_cff import run2_HCAL_2018
run2_HCAL_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder, recHitEnergyNorms = recHitEnergyNorms2018)
run2_HCAL_2018.toModify(particleFlowClusterHBHE.seedFinder, thresholdsByDetector = seedFinderThresholdsByDetector2018)
run2_HCAL_2018.toModify(particleFlowClusterHBHE.initialClusteringStep, thresholdsByDetector = initialClusteringStepThresholdsByDetector2018)
run2_HCAL_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder.positionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2018)
run2_HCAL_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2018)


from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder, recHitEnergyNorms = recHitEnergyNorms2018)
run2_HE_2018.toModify(particleFlowClusterHBHE.seedFinder, thresholdsByDetector = seedFinderThresholdsByDetector2018)
run2_HE_2018.toModify(particleFlowClusterHBHE.initialClusteringStep, thresholdsByDetector = initialClusteringStepThresholdsByDetector2018)
run2_HE_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder.positionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2018)
run2_HE_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2018)

# offline 2018 -- collapsed
run2_HECollapse_2018 =  cms.Modifier()
run2_HECollapse_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder, recHitEnergyNorms = recHitEnergyNorms2017)
run2_HECollapse_2018.toModify(particleFlowClusterHBHE.seedFinder, thresholdsByDetector = seedFinderThresholdsByDetector2017)
run2_HECollapse_2018.toModify(particleFlowClusterHBHE.initialClusteringStep, thresholdsByDetector = initialClusteringStepThresholdsByDetector2017)
run2_HECollapse_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder.positionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2017)
run2_HECollapse_2018.toModify(particleFlowClusterHBHE.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2017)


# offline 2019
Exemple #5
0
     UseRejectedHitsOnly = cms.bool(False),

#    Controls for hits to be included in the "bad" tower collection.
#    They have no effect unless UseRejectedHitsOnly=true                                
#    Hits passing the   HcalAcceptSeverityLevel
#    will be skipped as they are already in the default towers                                
     HcalAcceptSeverityLevelForRejectedHit = cms.uint32(9999),
#    List of ECAL problems that should be used in bad tower construction
#    Note that these can only be of type already rejected in default
#    reconstruction as specified in "EcalRecHitSeveritiesToBeExcluded"
     EcalSeveritiesToBeUsedInBadTowers = cms.vstring(),
                                

#    The code also checks the settings of the flags for the default
#    collection - if the recovered hits were used there, they
#    will be skipped for the "bad" tower collection regardless of these settings                                
     UseRejectedRecoveredHcalHits = cms.bool(True),
     UseRejectedRecoveredEcalHits = cms.bool(False),


# flag to allow/disallow missing inputs
    AllowMissingInputs = cms.bool(False),
	
# specify hcal upgrade phase - 0, 1, 2	
	HcalPhase = cms.int32(0)
    
)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(calotowermaker, HcalPhase = cms.int32(1))
               binOfMaximum     = cms.int32(2)
              )
)

from Configuration.Eras.Modifier_run2_HB_2018_cff import run2_HB_2018
run2_HB_2018.toModify( hcalSimParameters,
    hb = dict(
               readoutFrameSize = cms.int32(8), 
               binOfMaximum     = cms.int32(4)
              )
)
from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify( hcalSimParameters,
    he = dict(
               readoutFrameSize = cms.int32(8), 
               binOfMaximum     = cms.int32(4),
               threshold_currentTDC = cms.double(18.7)
              )
)


from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify( hcalSimParameters,
    hb = dict(
        doSiPMSmearing = cms.bool(True),
        threshold_currentTDC = cms.double(18.7),
        sipmTau = cms.double(10.),
    )
)

from Configuration.Eras.Modifier_run3_common_cff import run3_common
    ),
    cms.PSet(
        detector=cms.string("HCAL_ENDCAP"),
        depths=cms.vint32(1, 2, 3, 4, 5, 6, 7),
        logWeightDenominator=_thresholdsHE,
    ))

# offline 2018 -- uncollapsed
from Configuration.Eras.Modifier_run2_HCAL_2018_cff import run2_HCAL_2018
run2_HCAL_2018.toModify(
    particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2018)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(
    particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2018)

# offline 2018 -- collapsed
run2_HECollapse_2018 = cms.Modifier()
run2_HECollapse_2018.toModify(
    particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2017)

# offline 2019
from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(
    particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc,
    logWeightDenominatorByDetector=logWeightDenominatorByDetector2019)

# offline phase2 restore what has been studied in the TDR
Exemple #8
0
    HBPhase1=cms.bool(False),
    Plot_TP_ver=cms.bool(False))

from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(hcaldigisAnalyzer, simHits="fastSimProducer:HcalHits")

from Configuration.ProcessModifiers.premix_stage2_cff import premix_stage2
premix_stage2.toModify(
    hcaldigisAnalyzer,
    dataTPs="DMHcalTriggerPrimitiveDigis",
)

from Configuration.Eras.Modifier_run2_HCAL_2017_cff import run2_HCAL_2017
run2_HCAL_2017.toModify(hcaldigisAnalyzer, TestNumber=cms.bool(True))

from Configuration.Eras.Modifier_run2_HEPlan1_2017_cff import run2_HEPlan1_2017
run2_HEPlan1_2017.toModify(hcaldigisAnalyzer, hep17=cms.bool(True))

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(hcaldigisAnalyzer, HEPhase1=cms.bool(True))

from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
phase2_hcal.toModify(hcaldigisAnalyzer,
                     dataTPs=cms.InputTag(""),
                     digiTag=cms.InputTag("simHcalDigis"),
                     QIE10digiTag=cms.InputTag("simHcalDigis",
                                               "HFQIE10DigiCollection"),
                     QIE11digiTag=cms.InputTag("simHcalDigis",
                                               "HBHEQIE11DigiCollection"),
                     HBPhase1=cms.bool(True))
Exemple #9
0
    # NOTE: The following controls the creation of towers from
    #       rejected rechits.
    #       Always make sure that UseRejectedHitsOnly=false for
    #       normal reconstructions!!!
    UseRejectedHitsOnly=cms.bool(False),

    #    Controls for hits to be included in the "bad" tower collection.
    #    They have no effect unless UseRejectedHitsOnly=true
    #    Hits passing the   HcalAcceptSeverityLevel
    #    will be skipped as they are already in the default towers
    HcalAcceptSeverityLevelForRejectedHit=cms.uint32(9999),
    #    List of ECAL problems that should be used in bad tower construction
    #    Note that these can only be of type already rejected in default
    #    reconstruction as specified in "EcalRecHitSeveritiesToBeExcluded"
    EcalSeveritiesToBeUsedInBadTowers=cms.vstring(),

    #    The code also checks the settings of the flags for the default
    #    collection - if the recovered hits were used there, they
    #    will be skipped for the "bad" tower collection regardless of these settings
    UseRejectedRecoveredHcalHits=cms.bool(True),
    UseRejectedRecoveredEcalHits=cms.bool(False),

    # flag to allow/disallow missing inputs
    AllowMissingInputs=cms.bool(False),

    # specify hcal upgrade phase - 0, 1, 2
    HcalPhase=cms.int32(0))

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(calotowermaker, HcalPhase=cms.int32(1))
Exemple #10
0
from Configuration.Eras.Modifier_run2_HF_2017_cff import run2_HF_2017
run2_HF_2017.toModify(hcalSimParameters,
                      hf1=dict(readoutFrameSize=cms.int32(3),
                               binOfMaximum=cms.int32(2)),
                      hf2=dict(readoutFrameSize=cms.int32(3),
                               binOfMaximum=cms.int32(2)))

from Configuration.Eras.Modifier_run2_HB_2018_cff import run2_HB_2018
run2_HB_2018.toModify(hcalSimParameters,
                      hb=dict(readoutFrameSize=cms.int32(8),
                              binOfMaximum=cms.int32(4)))
from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(
    hcalSimParameters,
    he=dict(
        readoutFrameSize=cms.int32(8),
        binOfMaximum=cms.int32(4),
        threshold_currentTDC=cms.double(
            149.6)  #2.34)#9.35)#4.675)#9.35)#18.7)#37.4)
    ))

from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(
    hcalSimParameters,
    hb=dict(
        doSiPMSmearing=cms.bool(True),
        threshold_currentTDC=cms.double(149.6),  #18.7, 37.4, 56.1, 74.8
        sipmTau=cms.double(10.),
    ))

_newFactors = cms.vdouble(
    210.55, 197.93, 186.12, 189.64, 189.63, 189.96, 190.03, 190.11, 190.18,
Exemple #11
0
HcalTimeSlewEP = cms.ESSource("HcalTimeSlewEP",
    appendToDataLabel = cms.string("HBHE"),                              
    # For method2/simulation
    timeSlewParametersM2 = cms.VPSet(
        cms.PSet(#Slow    
            tzero = cms.double(23.960177), slope = cms.double(-3.178648), tmax = cms.double(16.00)),
        cms.PSet(#Medium
            tzero = cms.double(13.307784), slope = cms.double(-1.556668), tmax = cms.double(10.00)),
        cms.PSet(#Fast
            tzero = cms.double(9.109694),  slope = cms.double(-1.075824), tmax = cms.double(6.25))
    ),
    # For method3                          
    timeSlewParametersM3 = cms.VPSet(        
        cms.PSet(#TestStand (Parameters not used)
            cap = cms.double(6.0), tspar0 = cms.double(12.2999), tspar1 = cms.double(-2.19142), tspar2 = cms.double(0.0),  tspar0_siPM = cms.double(0.0), tspar1_siPM = cms.double(0.0), tspar2_siPM = cms.double(0.0)),
        cms.PSet(#Data
            cap = cms.double(6.0), tspar0 = cms.double(15.5),    tspar1 = cms.double(-3.2),     tspar2 = cms.double(32.0), tspar0_siPM = cms.double(0.0), tspar1_siPM = cms.double(0.0), tspar2_siPM = cms.double(0.0)),
        cms.PSet(#MC
            cap = cms.double(6.0), tspar0 = cms.double(12.2999), tspar1 = cms.double(-2.19142), tspar2 = cms.double(0.0),  tspar0_siPM = cms.double(0.0), tspar1_siPM = cms.double(0.0), tspar2_siPM = cms.double(0.0)),
        cms.PSet(#HBHE
            cap = cms.double(6.0), tspar0 = cms.double(12.2999), tspar1 = cms.double(-2.19142), tspar2 = cms.double(0.0),  tspar0_siPM = cms.double(0.0), tspar1_siPM = cms.double(0.0), tspar2_siPM = cms.double(0.0))
    )
)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
#HBHE2018  -  Medium
run2_HE_2018.toModify(
    HcalTimeSlewEP, timeSlewParametersM2 = {
        1: dict(tzero = cms.double(11.977461), slope = cms.double(-1.5610227),tmax = cms.double(10.00))})
Exemple #12
0
    # Reconstruction algorithm configuration data to fetch from DB, if any
    algoConfigClass = cms.string(""),

    # Turn rechit status bit setters on/off
    setNegativeFlagsQIE8 = cms.bool(True),
    setNegativeFlagsQIE11 = cms.bool(False),
    setNoiseFlagsQIE8 = cms.bool(True),
    setNoiseFlagsQIE11 = cms.bool(False),
    setPulseShapeFlagsQIE8 = cms.bool(True),
    setPulseShapeFlagsQIE11 = cms.bool(False),
    setLegacyFlagsQIE8 = cms.bool(True),
    setLegacyFlagsQIE11 = cms.bool(False),

    # Parameter sets configuring rechit status bit setters
    flagParametersQIE8 = cms.PSet(
        hbheStatusFlag.qie8Config
    ),
    flagParametersQIE11 = cms.PSet(),

    pulseShapeParametersQIE8 = cms.PSet(
        pulseShapeFlag.qie8Parameters
    ),
    pulseShapeParametersQIE11 = cms.PSet()
)

# Disable the "triangle peak fit" and the corresponding HBHETriangleNoise flag
hbheprereco.pulseShapeParametersQIE8.TrianglePeakTS = cms.uint32(10000)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(hbheprereco, saveEffectivePedestal = cms.bool(True))
    Plot_TP_ver = cms.bool(False)
)

from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toModify(hcaldigisAnalyzer, simHits = "fastSimProducer:HcalHits")

from Configuration.ProcessModifiers.premix_stage2_cff import premix_stage2
premix_stage2.toModify(hcaldigisAnalyzer,
    dataTPs = "DMHcalTriggerPrimitiveDigis",
)

from Configuration.Eras.Modifier_run2_HCAL_2017_cff import run2_HCAL_2017
(run2_HCAL_2017 & ~fastSim).toModify(hcaldigisAnalyzer,
    TestNumber    = cms.bool(True)
)

from Configuration.Eras.Modifier_run2_HEPlan1_2017_cff import run2_HEPlan1_2017
run2_HEPlan1_2017.toModify(hcaldigisAnalyzer,
    hep17 = cms.bool(True)
)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(hcaldigisAnalyzer,
    HEPhase1 = cms.bool(True)
)
    
from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(hcaldigisAnalyzer,
    HBPhase1 = cms.bool(True)
)
logWeightDenominatorByDetectorPhase2 = cms.VPSet(
    cms.PSet( detector = cms.string("HCAL_BARREL1"),
              depths = cms.vint32(1, 2, 3, 4),
              logWeightDenominator = _thresholdsHB,
              ),
    cms.PSet( detector = cms.string("HCAL_ENDCAP"),
              depths = cms.vint32(1, 2, 3, 4, 5, 6, 7),
              logWeightDenominator = _thresholdsHE,
              )
    )

# offline 2018 -- uncollapsed
from Configuration.Eras.Modifier_run2_HCAL_2018_cff import run2_HCAL_2018
run2_HCAL_2018.toModify(particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2018)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2018)

# offline 2018 -- collapsed
run2_HECollapse_2018 =  cms.Modifier()
run2_HECollapse_2018.toModify(particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2017)

# offline 2019
from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetector2019)

# offline phase2 restore what has been studied in the TDR
from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
phase2_hcal.toModify(particleFlowClusterHCAL.pfClusterBuilder.allCellsPositionCalc, logWeightDenominatorByDetector= logWeightDenominatorByDetectorPhase2)
Exemple #15
0
     UseRejectedRecoveredEcalHits = cms.bool(False),


# flag to allow/disallow missing inputs
    AllowMissingInputs = cms.bool(False),
	
# specify hcal upgrade phase - 0, 1, 2	
	HcalPhase = cms.int32(0)
    
)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(calotowermaker, 
                      HcalPhase = cms.int32(1),
                      HESThreshold1 = cms.double(0.1),
                      HESThreshold  = cms.double(0.2),
                      HEDThreshold1 = cms.double(0.1),
                      HEDThreshold  = cms.double(0.2)
)

# needed to handle inner/outer assignment
from Configuration.ProcessModifiers.run2_HECollapse_2018_cff import run2_HECollapse_2018
run2_HECollapse_2018.toModify(calotowermaker,
    HcalPhase = cms.int32(0),
    HESThreshold1 = cms.double(0.8),
    HESThreshold  = cms.double(0.8),
    HEDThreshold1 = cms.double(0.8),
    HEDThreshold  = cms.double(0.8)
)

from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
Exemple #16
0
                      ))

from Configuration.Eras.Modifier_run2_HF_2017_cff import run2_HF_2017
run2_HF_2017.toModify(hcalSimParameters,
                      hf1=dict(readoutFrameSize=cms.int32(3),
                               binOfMaximum=cms.int32(2)),
                      hf2=dict(readoutFrameSize=cms.int32(3),
                               binOfMaximum=cms.int32(2)))

from Configuration.Eras.Modifier_run2_HB_2018_cff import run2_HB_2018
run2_HB_2018.toModify(hcalSimParameters,
                      hb=dict(readoutFrameSize=cms.int32(8),
                              binOfMaximum=cms.int32(4)))
from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(hcalSimParameters,
                      he=dict(readoutFrameSize=cms.int32(8),
                              binOfMaximum=cms.int32(4)))

from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(hcalSimParameters,
                 hb=dict(
                     doSiPMSmearing=cms.bool(True),
                     sipmTau=cms.double(10.),
                 ))

_newFactors = cms.vdouble(
    210.55, 197.93, 186.12, 189.64, 189.63, 189.96, 190.03, 190.11, 190.18,
    190.25, 190.32, 190.40, 190.47, 190.54, 190.61, 190.69, 190.83, 190.94,
    190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94,
    190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94,
    190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94, 190.94,
Exemple #17
0
# set useConfigZSvalues = cms.int32(0) - default

import FWCore.ParameterSet.Config as cms

simHcalDigis = cms.EDProducer("HcalRealisticZS",
                              digiLabel=cms.string("simHcalUnsuppressedDigis"),
                              markAndPass=cms.bool(False),
                              HBlevel=cms.int32(8),
                              HElevel=cms.int32(9),
                              HOlevel=cms.int32(24),
                              HFlevel=cms.int32(-9999),
                              HBregion=cms.vint32(3, 6),
                              HEregion=cms.vint32(3, 6),
                              HOregion=cms.vint32(1, 8),
                              HFregion=cms.vint32(2, 3),
                              useConfigZSvalues=cms.int32(0))

from Configuration.Eras.Modifier_phase2_hcal_cff import phase2_hcal
phase2_hcal.toModify(simHcalDigis,
                     useConfigZSvalues=cms.int32(1),
                     HElevel=cms.int32(3))

from Configuration.Eras.Modifier_run2_HF_2017_cff import run2_HF_2017
run2_HF_2017.toModify(simHcalDigis, HFregion=cms.vint32(1, 2))

from Configuration.Eras.Modifier_run2_HB_2018_cff import run2_HB_2018
run2_HB_2018.toModify(simHcalDigis, HBregion=cms.vint32(2, 5))

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(simHcalDigis, HEregion=cms.vint32(2, 5))
    HOregion = cms.vint32(1,8),
    HFregion = cms.vint32(2,3),      
    useConfigZSvalues = cms.int32(0)
)

from Configuration.Eras.Modifier_run2_HF_2017_cff import run2_HF_2017
run2_HF_2017.toModify( simHcalDigis,
                             HFregion = cms.vint32(1,2)
)

from Configuration.Eras.Modifier_run2_HB_2018_cff import run2_HB_2018
run2_HB_2018.toModify( simHcalDigis,
                             HBregion = cms.vint32(2,5)
)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify( simHcalDigis,
                             HEregion = cms.vint32(2,5)
)

# Switch off HCAL ZS in digi for premixing stage1
from Configuration.ProcessModifiers.premix_stage1_cff import premix_stage1
premix_stage1.toModify(simHcalDigis,
    markAndPass = True,
    HBlevel = -999,
    HElevel = -999,
    HOlevel = -999,
    HFlevel = -999,
    useConfigZSvalues = 1
)
    UseRejectedRecoveredHcalHits=cms.bool(True),
    UseRejectedRecoveredEcalHits=cms.bool(False),

    #    If Hcal is masked, and Ecal is present, pretend Hcal = (this factor) * Ecal
    missingHcalRescaleFactorForEcal=cms.double(0),

    # flag to allow/disallow missing inputs
    AllowMissingInputs=cms.bool(False),

    # specify hcal upgrade phase - 0, 1, 2
    HcalPhase=cms.int32(0))

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(calotowermaker,
                      HcalPhase=cms.int32(1),
                      HESThreshold1=cms.double(0.1),
                      HESThreshold=cms.double(0.2),
                      HEDThreshold1=cms.double(0.1),
                      HEDThreshold=cms.double(0.2))

# needed to handle inner/outer assignment
from Configuration.ProcessModifiers.run2_HECollapse_2018_cff import run2_HECollapse_2018
run2_HECollapse_2018.toModify(calotowermaker,
                              HcalPhase=cms.int32(0),
                              HESThreshold1=cms.double(0.8),
                              HESThreshold=cms.double(0.8),
                              HEDThreshold1=cms.double(0.8),
                              HEDThreshold=cms.double(0.8))

from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
run3_HB.toModify(
    calotowermaker,
Exemple #20
0
    UseRejectedRecoveredEcalHits=cms.bool(False),

    #    If Hcal is masked, and Ecal is present, pretend Hcal = (this factor) * Ecal
    missingHcalRescaleFactorForEcal=cms.double(0),

    # flag to allow/disallow missing inputs
    AllowMissingInputs=cms.bool(False),

    # specify hcal upgrade phase - 0, 1, 2
    HcalPhase=cms.int32(0))

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018

run2_HE_2018.toModify(calotowermaker,
                      HcalPhase=1,
                      HESThreshold1=0.1,
                      HESThreshold=0.2,
                      HEDThreshold1=0.1,
                      HEDThreshold=0.2)

# needed to handle inner/outer assignment
from Configuration.ProcessModifiers.run2_HECollapse_2018_cff import run2_HECollapse_2018

run2_HECollapse_2018.toModify(calotowermaker,
                              HcalPhase=0,
                              HESThreshold1=0.8,
                              HESThreshold=0.8,
                              HEDThreshold1=0.8,
                              HEDThreshold=0.8)

from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
Exemple #21
0
    caloLabel=cms.InputTag("calotowermaker"),
    HF1Weights=cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    HESEScale=cms.double(50.0),
    HESGrid=cms.vdouble(0.0, 2.0, 4.0, 5.0, 9.0, 20.0, 30.0, 50.0, 100.0,
                        1000.0),
    HEDEScale=cms.double(50.0),
    HESWeights=cms.vdouble(1.7, 1.57, 1.54, 1.49, 1.41, 1.26, 1.19, 1.15, 1.12,
                           1.0),
    HEDGrid=cms.vdouble(0.0, 2.0, 4.0, 5.0, 9.0, 20.0, 30.0, 50.0, 100.0,
                        1000.0),
    EBEScale=cms.double(50.0),
    # Weighting factor for HB
    HBWeight=cms.double(1.0),
    EEEScale=cms.double(50.0),
    HOGrid=cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    # Energy dependent weights and energy scale to be used
    EBGrid=cms.vdouble(2.0, 4.0, 5.0, 9.0, 20.0, 30.0, 50.0, 100.0, 300.0),
    # momentum assignment
    # Method for momentum reconstruction
    MomConstrMethod=cms.int32(1),
    # Depth, fraction of the respective calorimeter [0,1]
    MomHBDepth=cms.double(0.2),
    MomHEDepth=cms.double(0.4),
    MomEBDepth=cms.double(0.3),
    MomEEDepth=cms.double(0.0),
    HcalPhase=cms.int32(0))

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018

run2_HE_2018.toModify(ct2ct, HcalPhase=cms.int32(1))
Exemple #22
0
    HESEScale = cms.double(50.0),
    HESGrid = cms.vdouble(0.0, 2.0, 4.0, 5.0, 9.0, 
        20.0, 30.0, 50.0, 100.0, 1000.0),
    HEDEScale = cms.double(50.0),
    HESWeights = cms.vdouble(1.7, 1.57, 1.54, 1.49, 1.41, 
        1.26, 1.19, 1.15, 1.12, 1.0),
    HEDGrid = cms.vdouble(0.0, 2.0, 4.0, 5.0, 9.0, 
        20.0, 30.0, 50.0, 100.0, 1000.0),
    EBEScale = cms.double(50.0),
    # Weighting factor for HB   
    HBWeight = cms.double(1.0),
    EEEScale = cms.double(50.0),
    HOGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    # Energy dependent weights and energy scale to be used
    EBGrid = cms.vdouble(2.0, 4.0, 5.0, 9.0, 20.0, 
        30.0, 50.0, 100.0, 300.0),
    # momentum assignment
    # Method for momentum reconstruction
    MomConstrMethod = cms.int32(1),                           
    # Depth, fraction of the respective calorimeter [0,1]
    MomHBDepth = cms.double(0.2),
    MomHEDepth = cms.double(0.4),   
    MomEBDepth = cms.double(0.3),
    MomEEDepth = cms.double(0.0),
	HcalPhase = cms.int32(0)
)

from Configuration.Eras.Modifier_run2_HE_2018_cff import run2_HE_2018
run2_HE_2018.toModify(ct2ct, HcalPhase = 1)