calotowermaker_cfi.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

calotowermaker = cms.EDProducer("CaloTowersCreator",
    # Energy threshold for EB 5x5 crystal inclusion [GeV]
    EBSumThreshold = cms.double(0.2),
    # Weighting factor for HF short-fiber readouts
    HF2Weight = cms.double(1.0),
    # Weighting factor for EB   
    EBWeight = cms.double(1.0),
    # Label of HFRecHitCollection to use
    hfInput = cms.InputTag("hfreco"),
    # Energy threshold for EE crystals-in-tower inclusion [GeV]
    EESumThreshold = cms.double(0.45),
    # Energy threshold for HO cell inclusion [GeV]
    HOThreshold0 = cms.double(1.1),
    HOThresholdPlus1 = cms.double(3.5),
    HOThresholdMinus1 = cms.double(3.5),
    HOThresholdPlus2 = cms.double(3.5),
    HOThresholdMinus2 = cms.double(3.5),
    HBGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    # Energy threshold for HB cell inclusion [GeV]
    HBThreshold1 = cms.double(0.7), # depth 1
    HBThreshold2 = cms.double(0.7), # depth 2
    HBThreshold = cms.double(0.7), # depths 3-4
    EEWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    # Energy threshold for long-fiber HF readout inclusion [GeV]
    HF1Threshold = cms.double(0.5),
    HF2Weights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    HOWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    EEGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    # Weighting factor for HE 10-degree cells   
    HEDWeight = cms.double(1.0),
    # Weighting factor for EE   
    EEWeight = cms.double(1.0),
    # HO on/off flag for tower energy reconstruction
    UseHO = cms.bool(True),
    HBWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    # Weighting factor for HE 5-degree cells   
    HESWeight = cms.double(1.0),
    # Weighting factor for HF long-fiber readouts 
    HF1Weight = cms.double(1.0),
    HF2Grid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    HEDWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    HF1Grid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    EBWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    # Weighting factor for HO 
    HOWeight = cms.double(1.0),
                              
    # Energy threshold for EB crystal inclusion [GeV]
    EBThreshold = cms.double(0.07),
    # Energy threshold for EE crystal inclusion [GeV]
    EEThreshold = cms.double(0.3),
    # Flags specifying if the above thresholds
    # should be applied to Et (UseEtEXTreshold='True') or E ('False')
    # Flags for use of symmetric thresholds: |e|>threshold                          
    UseEtEBTreshold = cms.bool(False),
    UseSymEBTreshold = cms.bool(True),
    UseEtEETreshold = cms.bool(False),
    UseSymEETreshold = cms.bool(True),


    # Label of HBHERecHitCollection to use
    hbheInput = cms.InputTag("hbhereco"),
    # Global energy threshold on Hcal [GeV]
    HcalThreshold = cms.double(-1000.0),
    # Energy threshold for short-fiber HF readout inclusion [GeV]
    HF2Threshold = cms.double(0.85),

    # Energy threshold for 5-degree (phi) HE cell inclusion [GeV]
    HESThreshold1 = cms.double(0.8), # depth  1    
    HESThreshold  = cms.double(0.8), # depths 2-7
    HF1Weights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    # Label of HORecHitCollection to use
    hoInput = cms.InputTag("horeco"),
    HESGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    #
    HESWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
    # Energy threshold for 10-degree (phi) HE cel inclusion [GeV]
    HEDThreshold1 = cms.double(0.8), # depth  1  
    HEDThreshold  = cms.double(0.8), # depths 2-7
    # Global energy threshold on tower [GeV]
    EcutTower = cms.double(-1000.0),
    HEDGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
    # Label of EcalRecHitCollections to use
    ecalInputs = cms.VInputTag(cms.InputTag("ecalRecHit","EcalRecHitsEB"), cms.InputTag("ecalRecHit","EcalRecHitsEE")),
    # Weighting factor for HB   
    HBWeight = cms.double(1.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(-1.0, 1.0, 10.0, 100.0, 1000.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),

# parameters for handling of anomalous cells
# 
    # acceptable severity level
    HcalAcceptSeverityLevel = cms.uint32(9),
    #EcalAcceptSeverityLevel = cms.uint32(1),
    EcalRecHitSeveritiesToBeExcluded = cms.vstring('kTime','kWeird','kBad'),
                                
    # use of recovered hits
    UseHcalRecoveredHits = cms.bool(True),
      # The CaloTower code treats recovered cells as a separate category.
      # The flag to use (or not use) them should be explicitly set
      # regardless of the specified severity level in EcalAcceptSeverityLevel.                       
    UseEcalRecoveredHits = cms.bool(False),
                        

# 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),

#    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),

# Read HBHE thresholds from Global Tag
    usePFThresholdsFromDB = cms.bool(False),
# Read ECAL thresholds from Global Tag
    EcalRecHitThresh = cms.bool(False)
)

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
run3_HB.toModify(calotowermaker,
    HBThreshold1 = 0.1,
    HBThreshold2 = 0.2,
    HBThreshold = 0.3,
)

#--- Use DB conditions for HBHE thresholds for Run3 and phase2
from Configuration.Eras.Modifier_hcalPfCutsFromDB_cff import hcalPfCutsFromDB
hcalPfCutsFromDB.toModify( calotowermaker,
                           usePFThresholdsFromDB = True)