customise_phase1_hcal Namespace Reference

Functions
def	customise
def	customise_3bx
def	customise_pu15_25ns
	Known alterations for Reco #####################################################. More...
def	customise_pu15_25ns_3bx
def	customise_pu15_25ns_wo_pairs
def	customise_pu50_25ns
def	customise_pu50_25ns_3bx
	if doing inefficiency at <PU>=50 More...
def	customise_pu50_25ns_wo_pairs
def	customise_wo_pairs
def	saveAllProds
	if doing inefficiency at <PU>=50 More...
def	turnOffLegacyPhiCorrection

Variables
tuple	averageNumber = cms.double(15.0)
	set the number of pileup More...

Function Documentation

def customise_phase1_hcal.customise

(

process

)

Definition at line 7 of file customise_phase1_hcal.py.

References reco.return().

    process.ctfWithMaterialTracks.TTRHBuilder = 'WithTrackAngle'
    process.PixelCPEGenericESProducer.UseErrorsFromTemplates = cms.bool(False)
    process.PixelCPEGenericESProducer.TruncatePixelCharge = cms.bool(False)
    process.PixelCPEGenericESProducer.LoadTemplatesFromDB = cms.bool(False)
    process.PixelCPEGenericESProducer.Upgrade = cms.bool(True)
    process.PixelCPEGenericESProducer.SmallPitch = False
    process.PixelCPEGenericESProducer.IrradiationBiasCorrection = False
    process.PixelCPEGenericESProducer.DoCosmics = False

    ## CPE for other steps
    process.siPixelRecHits.CPE = cms.string('PixelCPEGeneric')
    process.initialStepTracks.TTRHBuilder = cms.string('WithTrackAngle')
    process.lowPtTripletStepTracks.TTRHBuilder = cms.string('WithTrackAngle')
    process.pixelPairStepTracks.TTRHBuilder = cms.string('WithTrackAngle')
    process.detachedTripletStepTracks.TTRHBuilder = cms.string('WithTrackAngle')
    process.mixedTripletStepTracks.TTRHBuilder = cms.string('WithTrackAngle')
    process.pixelLessStepTracks.TTRHBuilder = cms.string('WithTrackAngle')
    process.tobTecStepTracks.TTRHBuilder = cms.string('WithTrackAngle')
    process.highPtTripletStepTracks.TTRHBuilder = cms.string('WithTrackAngle')

    # Need these lines to stop some errors about missing siStripDigis collections.
    # should add them to fakeConditions_Phase1_cff
    process.MeasurementTracker.inactiveStripDetectorLabels = cms.VInputTag()
    process.MeasurementTracker.UseStripModuleQualityDB     = cms.bool(False)
    process.MeasurementTracker.UseStripAPVFiberQualityDB   = cms.bool(False)
    process.MeasurementTracker.UseStripStripQualityDB      = cms.bool(False)
    process.MeasurementTracker.UsePixelModuleQualityDB     = cms.bool(False)
    process.MeasurementTracker.UsePixelROCQualityDB        = cms.bool(False)

    process.muons.TrackerKinkFinderParameters.TrackerRecHitBuilder = cms.string('WithTrackAngle')
    # The SeedMergerPSet should be added to the following file for Phase 1
    # RecoTracker/SpecialSeedGenerators/python/CombinatorialSeedGeneratorForCosmicsRegionalReconstruction_cfi.py
    # but pixel layers are not used here for cosmic TODO: need Maria and Jan to do appropriate thing here
    from RecoPixelVertexing.PixelTriplets.quadrupletseedmerging_cff import PixelSeedMergerQuadruplets
    process.regionalCosmicTrackerSeeds.SeedMergerPSet = cms.PSet(
                mergeTriplets = cms.bool(False),
                        ttrhBuilderLabel = cms.string( "PixelTTRHBuilderWithoutAngle" ),
                        addRemainingTriplets = cms.bool(False),
                        layerList = PixelSeedMergerQuadruplets
                        )
    process.regionalCosmicTracks.TTRHBuilder = cms.string('WithTrackAngle')


    ### back to standard job commands ##################################################
    process.DigiToRaw.remove(process.castorRawData)

    process.DigiToRaw.remove(process.siPixelRawData)
    process.RawToDigi.remove(process.siPixelDigis)


    if hasattr(process,'dqmoffline_step'):
        process.dqmoffline_step.remove(process.SiPixelTrackResidualSource)
        process.dqmoffline_step.remove(process.jetMETAnalyzer)
        process.dqmoffline_step.remove(process.hltMonMuBits)
        process.dqmoffline_step.remove(process.vbtfAnalyzer)
        process.dqmoffline_step.remove(process.hltResults)
        process.dqmoffline_step.remove(process.egHLTOffDQMSource)
        process.dqmoffline_step.remove(process.globalAnalyzer)
        process.dqmoffline_step.remove(process.jetMETHLTOfflineSource)
        ##
        process.dqmoffline_step.remove(process.TrackerCollisionTrackMon)
    if hasattr(process,'validation_step'):
        process.validation_step.remove(process.hltHITval)
        process.validation_step.remove(process.HLTSusyExoVal)
        process.validation_step.remove(process.relvalMuonBits)
        process.validation_step.remove(process.hltMuonValidator)
        #this takes forever - seems like an infinite loop that I didnt debug yet
        process.validation_step.remove(process.pixelDigisValid)
    else:
    ## removing large memory usage module if we don't need it
        process.pdigi.remove(process.mergedtruth)

#  HCAL Upgrade Geometry

    process.load("RecoLocalCalo.HcalRecProducers.HcalUpgradeReconstructor_cff")
    process.load("RecoJets.Configuration.CaloTowersRec_cff")
    process.load("RecoLocalCalo.HcalRecAlgos.hcalRecAlgoESProd_cfi")
    process.load("RecoLocalCalo.Configuration.RecoLocalCalo_cff")

    process.ecalGlobalUncalibRecHit.EBdigiCollection = cms.InputTag("simEcalDigis","ebDigis")
    process.ecalGlobalUncalibRecHit.EEdigiCollection = cms.InputTag("simEcalDigis","eeDigis")
    process.ecalRecHit.ebDetIdToBeRecovered = cms.InputTag("","")
    process.ecalRecHit.eeDetIdToBeRecovered = cms.InputTag("","")
    process.ecalRecHit.eeFEToBeRecovered = cms.InputTag("","")
    process.ecalRecHit.ebFEToBeRecovered = cms.InputTag("","")
    process.ecalRecHit.recoverEBFE = cms.bool(False)
    process.ecalRecHit.recoverEEFE = cms.bool(False)

    process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hbhe_cfi")
    process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_ho_cfi")
    process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_hf_cfi")

    process.hbheprereco.digiLabel = "simHcalUnsuppressedDigis"
    process.horeco.digiLabel = "simHcalUnsuppressedDigis"
    process.hfreco.digiLabel = "simHcalUnsuppressedDigis"
    process.hcalupgradereco.digiLabel = "simHcalUnsuppressedDigis"

### Known alterations for Reco #####################################################
    delattr(process,"hbhereco")
    process.hbhereco = process.hbheprereco.clone()

    process.metrecoPlusHCALNoise.remove(process.BeamHaloSummary)
    process.metrecoPlusHCALNoise.remove(process.GlobalHaloData)
    process.metrecoPlusHCALNoise.remove(process.HcalHaloData)

    process.metrecoPlusHCALNoise.remove(process.hcalnoise)
    process.jetGlobalReco = cms.Sequence(process.recoJets*process.recoTrackJets)
    process.jetHighLevelReco = cms.Sequence(process.recoJetAssociations*process.recoPFJets*process.recoJPTJets)

### Place to add in the reco steps one by one ######################################
    process.calolocalreco = cms.Sequence(process.ecalGlobalUncalibRecHit+
                                process.ecalDetIdToBeRecovered+
                                process.ecalRecHit+
                                process.ecalCompactTrigPrim+
                                process.ecalTPSkim+
                                process.ecalPreshowerRecHit+
                                #process.hbheprereco+
                                process.hbhereco+
                                process.horeco+process.hfreco+process.hcalupgradereco+process.towerMaker
                                #+process.zdcreco
                                )
    process.localreco  = cms.Sequence(process.trackerlocalreco+
                                process.muonlocalreco+
                                process.calolocalreco+
                                process.castorreco+
                                process.lumiProducer
                                )
    process.globalreco = cms.Sequence(process.offlineBeamSpot
                          *process.recopixelvertexing
                          *process.trackingGlobalReco
                          #*process.hcalGlobalRecoSequence 
                          *process.particleFlowCluster
                          *process.ecalClusters
                          *process.caloTowersRec
                          *process.vertexreco
                          *process.egammaGlobalReco
                          *process.pfTrackingGlobalReco
                          *process.jetGlobalReco
                          *process.muonrecoComplete
                          *process.muoncosmicreco
                          *process.CastorFullReco
                          )
    process.highlevelreco = cms.Sequence(process.egammaHighLevelRecoPrePF
                             *process.particleFlowReco
                             *process.egammaHighLevelRecoPostPF
                             *process.jetHighLevelReco
                             *process.tautagging
                             *process.metrecoPlusHCALNoise
                             *process.btagging
                             *process.recoPFMET
                             *process.PFTau
                             *process.regionalCosmicTracksSeq
                             *process.muoncosmichighlevelreco
                             *process.reducedRecHits
                             )
    process.reconstruction = cms.Sequence(  process.localreco       *
                                        process.globalreco      *
                                        process.highlevelreco   *
                                        process.logErrorHarvester
                                        )
    # Clean out some unused sequences
    process.hcalLocalRecoSequence=cms.Sequence()
    process.hcalGlobalRecoSequence=cms.Sequence()

    return(process)


def customise_pu15_25ns(process):

    process=customise(process)

    process.load("SLHCUpgradeSimulations.Geometry.mixLowLumPU_Phase1_R30F12_HCal_cff")

### set the number of pileup
    process.mix.input.nbPileupEvents = cms.PSet(
        averageNumber = cms.double(15.0)
        )
    return (process)

def customise_3bx(process):
    process.mix.maxBunch=1
    process.mix.minBunch=-1
    return process

def customise_pu15_25ns_3bx(process):
    process=customise_pu15_25ns(process)
    process=customise_3bx(process)
    return process

#pileup specific stuff here
def customise_pu50_25ns(process):

    process=customise(process)

    process.load("SLHCUpgradeSimulations.Geometry.mixLowLumPU_Phase1_R30F12_HCal_cff")

### set the number of pileup
    process.mix.input.nbPileupEvents = cms.PSet(
        averageNumber = cms.double(50.0)