Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95
import FWCore.ParameterSet.Config as cms

# ------------------------------------------------------------------------------
# configure a filter to run only on the events selected by TkAlMinBias AlcaReco
from HLTrigger.HLTfilters.hltHighLevel_cfi import *
ALCARECOCalMinBiasFilterForSiStripHitEff = hltHighLevel.clone(
    HLTPaths = ['pathALCARECOSiStripCalMinBias'],
    throw = True, ## throw on unknown path names
    TriggerResultsTag = ("TriggerResults","","RECO")
)

# ------------------------------------------------------------------------------
# This is the sequence for track refitting of the track saved by SiStripCalMinBias
# to have access to transient objects produced during RECO step and not saved
from Alignment.CommonAlignmentProducer.AlignmentTrackSelector_cfi import *
ALCARECOMonitoringTracks = AlignmentTrackSelector.clone(
    #    src = 'generalTracks',
    src = 'ALCARECOSiStripCalMinBias',
    filter = True,
    applyBasicCuts = True,
    ptMin = 0.8,
    nHitMin = 6,
    chi2nMax = 10.)

# FIXME: the beam-spot should be kept in the AlCaReco (if not already there) and dropped from here
from RecoVertex.BeamSpotProducer.BeamSpot_cff import *
from RecoTracker.IterativeTracking.InitialStep_cff import *
from RecoTracker.Configuration.RecoTrackerP5_cff import *
from RecoTracker.TrackProducer.TrackRefitter_cfi import *
from DQM.SiStripCommon.TkHistoMap_cff import *
from RecoTracker.MeasurementDet.MeasurementTrackerEventProducer_cfi import *

ALCARECOMonitoringTracksRefit = TrackRefitter.clone(
         src = "ALCARECOMonitoringTracks",
         NavigationSchool = cms.string("")
)

# ------------------------------------------------------------------------------
# refit and BS can be dropped if done together with RECO.
# track filter can be moved in acalreco if no otehr users
ALCARECOTrackFilterRefit = cms.Sequence(ALCARECOMonitoringTracks +
                                        MeasurementTrackerEvent +
                                        offlineBeamSpot +
                                        ALCARECOMonitoringTracksRefit)

# ------------------------------------------------------------------------------
# This is the module actually doing the calibration
from CalibTracker.SiStripHitEfficiency.siStripHitEfficiencyWorker_cfi import siStripHitEfficiencyWorker
ALCARECOSiStripHitEff =  siStripHitEfficiencyWorker.clone(
    dqmDir = "AlCaReco/SiStripHitEfficiency",
    lumiScalers= "scalersRawToDigi",
    addLumi = True,
    commonMode = "siStripDigis:CommonMode",
    addCommonMode= False,
    combinatorialTracks = "ALCARECOMonitoringTracksRefit",
    trajectories        = "ALCARECOMonitoringTracksRefit",
    siStripClusters     = "siStripClusters",
    siStripDigis        = "siStripDigis",
    trackerEvent        = "MeasurementTrackerEvent",
    # part 2
    Layer = 0, # =0 means do all layers
    Debug = True,
    # do not cut on the total number of tracks
    cutOnTracks = False,
    trackMultiplicity = 1000,
    # use or not first and last measurement of a trajectory (biases), default is false
    useFirstMeas = False,
    useLastMeas = False,
    useAllHitsFromTracksWithMissingHits = False,
    doMissingHitsRecovery = False,
    ## non-default settings
    ClusterMatchingMethod = 4,  # default 0  case0,1,2,3,4
    ClusterTrajDist       = 15, # default 64
)

from Configuration.Eras.Modifier_run3_common_cff import run3_common
run3_common.toModify(ALCARECOSiStripHitEff,
                     useAllHitsFromTracksWithMissingHits = True,
                     doMissingHitsRecovery = True)

# ----------------------------------------------------------------------------
MEtoEDMConvertSiStripHitEff = cms.EDProducer("MEtoEDMConverter",
                                             Name = cms.untracked.string('MEtoEDMConverter'),
                                             Verbosity = cms.untracked.int32(0), # 0 provides no output
                                             # 1 provides basic output
                                             # 2 provide more detailed output
                                             Frequency = cms.untracked.int32(50),
                                             MEPathToSave = cms.untracked.string('AlCaReco/SiStripHitEfficiency'))

# The actual sequence
seqALCARECOPromptCalibProdSiStripHitEfficiency = cms.Sequence(
    ALCARECOCalMinBiasFilterForSiStripHitEff *
    ALCARECOTrackFilterRefit *
    ALCARECOSiStripHitEff *
    MEtoEDMConvertSiStripHitEff)