HLTEvF/python/HLTSiStripMonitoring_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from TrackingTools.RecoGeometry.RecoGeometries_cff import *
0004 hltESPDummyDetLayerGeometry = DummyDetLayerGeometry.clone(
0005     ComponentName =  "hltESPDummyDetLayerGeometry"
0006 )
0007
0008 # strip cluster monitor
0009 # in order to get the crossing angle through the sensor
0010 # the track trajectory is needed
0011 # => the track re-fit has to be performed
0012 # => some ESProducer have to be defined
0013
0014 ## NB: the following ESProducer should be the same used in the HLT menu
0015 ##     make sure they are not already defined somewhereelse in the final configuration
0016 from RecoLocalTracker.SiPixelRecHits.PixelCPETemplateReco_cfi import templates
0017 hltESPPixelCPETemplateReco = templates.clone(
0018   LoadTemplatesFromDB =  True ,
0019   ComponentName = "hltESPPixelCPETemplateReco" ,
0020   Alpha2Order =  True ,
0021   ClusterProbComputationFlag =  0 ,
0022   speed =  -2 ,
0023   UseClusterSplitter =  False
0024 )
0025
0026 from RecoLocalTracker.SiPixelRecHits.PixelCPEGeneric_cfi import PixelCPEGenericESProducer
0027 hltESPPixelCPEGeneric = PixelCPEGenericESProducer.clone(
0028   EdgeClusterErrorX =  50.0 ,
0029   DoCosmics =  False ,
0030   LoadTemplatesFromDB =  True ,
0031   UseErrorsFromTemplates =  True ,
0032   eff_charge_cut_highX =  1.0 ,
0033   TruncatePixelCharge =  True ,
0034   size_cutY =  3.0 ,
0035   size_cutX =  3.0 ,
0036   inflate_all_errors_no_trk_angle =  False ,
0037   IrradiationBiasCorrection =  False ,
0038   inflate_errors = False ,
0039   eff_charge_cut_lowX =  0.0 ,
0040   eff_charge_cut_highY =  1.0 ,
0041   ClusterProbComputationFlag = 0 ,
0042   EdgeClusterErrorY = 85.0 ,
0043   ComponentName = "hltESPPixelCPEGeneric" ,
0044   eff_charge_cut_lowY = 0.0 ,
0045   Alpha2Order =  True
0046 )
0047
0048 from RecoTracker.TransientTrackingRecHit.TTRHBuilderWithTemplate_cfi import TTRHBuilderAngleAndTemplate
0049 hltESPTTRHBuilderAngleAndTemplate = TTRHBuilderAngleAndTemplate.clone(
0050   StripCPE = "hltESPStripCPEfromTrackAngle" ,
0051   Matcher = "StandardMatcher" ,
0052   ComputeCoarseLocalPositionFromDisk = False ,
0053   PixelCPE = "hltESPPixelCPETemplateReco",
0054   ComponentName ="hltESPTTRHBuilderAngleAndTemplate"
0055 )
0056 hltESPTTRHBWithTrackAngle = TTRHBuilderAngleAndTemplate.clone(
0057   StripCPE = "hltESPStripCPEfromTrackAngle",
0058   Matcher = "StandardMatcher",
0059   ComputeCoarseLocalPositionFromDisk = False,
0060   PixelCPE = "hltESPPixelCPEGeneric",
0061   ComponentName = "hltESPTTRHBWithTrackAngle"
0062 )
0063
0064 from RecoLocalTracker.SiStripRecHitConverter.StripCPEESProducer_cfi import stripCPEESProducer
0065 hltESPStripCPEfromTrackAngle = stripCPEESProducer.clone(
0066   ComponentType = "StripCPEfromTrackAngle" ,
0067   ComponentName = "hltESPStripCPEfromTrackAngle",
0068   parameters = cms.PSet(
0069     mLC_P2 = cms.double(0.3),
0070     mLC_P1 = cms.double(0.618),
0071     mLC_P0 = cms.double(-0.326),
0072 #    useLegacyError = cms.bool( True ), # 50ns menu
0073 #    maxChgOneMIP = cms.double( -6000.0 ), # 50ns menu
0074     useLegacyError = cms.bool(False) , # 25ns menu
0075     maxChgOneMIP = cms.double(6000.0) , #25ns menu
0076     mTEC_P1 = cms.double( 0.471 ),
0077     mTEC_P0 = cms.double( -1.885 ),
0078     mTOB_P0 = cms.double( -1.026 ),
0079     mTOB_P1 = cms.double( 0.253 ),
0080     mTIB_P0 = cms.double( -0.742 ),
0081     mTIB_P1 = cms.double( 0.202 ),
0082     mTID_P0 = cms.double( -1.427 ),
0083     mTID_P1 = cms.double( 0.433 )
0084   )
0085 )
0086
0087 from RecoTracker.TkNavigation.NavigationSchoolESProducer_cfi import navigationSchoolESProducer
0088 navigationSchoolESProducer = navigationSchoolESProducer.clone(
0089   ComponentName = "SimpleNavigationSchool" ,
0090   SimpleMagneticField ="ParabolicMf"
0091 )
0092
0093 from RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi import MeasurementTracker
0094 hltESPMeasurementTracker = MeasurementTracker.clone(
0095   UseStripStripQualityDB =  True ,
0096   StripCPE = "hltESPStripCPEfromTrackAngle",
0097   UsePixelROCQualityDB = True ,
0098   DebugPixelROCQualityDB = False,
0099   UseStripAPVFiberQualityDB = True ,
0100   badStripCuts = cms.PSet(
0101     TOB = cms.PSet(
0102       maxConsecutiveBad = cms.uint32( 9999 ),
0103       maxBad = cms.uint32( 9999 )
0104     ),
0105     TID = cms.PSet(
0106       maxConsecutiveBad = cms.uint32( 9999 ),
0107       maxBad = cms.uint32( 9999 )
0108     ),
0109     TEC = cms.PSet(
0110       maxConsecutiveBad = cms.uint32( 9999 ),
0111       maxBad = cms.uint32( 9999 )
0112     ),
0113     TIB = cms.PSet(
0114       maxConsecutiveBad = cms.uint32( 9999 ),
0115       maxBad = cms.uint32( 9999 )
0116     )
0117   ),
0118   DebugStripModuleQualityDB = False ,
0119   ComponentName ="hltESPMeasurementTracker",
0120   DebugPixelModuleQualityDB = False,
0121   UsePixelModuleQualityDB = True,
0122   DebugStripAPVFiberQualityDB = False,
0123   HitMatcher = "StandardMatcher",
0124   DebugStripStripQualityDB = False,
0125   PixelCPE = "hltESPPixelCPEGeneric",
0126   SiStripQualityLabel = "" ,
0127   UseStripModuleQualityDB = True,
0128   MaskBadAPVFibers = True
0129 )
0130
0131 hltESPRungeKuttaTrackerPropagator = cms.ESProducer( "PropagatorWithMaterialESProducer",
0132   SimpleMagneticField = cms.string( "" ),
0133   PropagationDirection = cms.string( "alongMomentum" ),
0134   ComponentName = cms.string( "hltESPRungeKuttaTrackerPropagator" ),
0135   Mass = cms.double( 0.105 ),
0136   ptMin = cms.double( -1.0 ),
0137   MaxDPhi = cms.double( 1.6 ),
0138   useRungeKutta = cms.bool( True )
0139 )
0140
0141 from TrackingTools.KalmanUpdators.KFUpdatorESProducer_cfi import KFUpdatorESProducer
0142 hltESPKFUpdator = KFUpdatorESProducer.clone(
0143   ComponentName = "hltESPKFUpdator"
0144 )
0145
0146 hltESPChi2MeasurementEstimator30 = cms.ESProducer( "Chi2MeasurementEstimatorESProducer",
0147   appendToDataLabel = cms.string( "" ),
0148   MinimalTolerance = cms.double( 10.0 ),
0149   MaxDisplacement = cms.double( 100.0 ),
0150   ComponentName = cms.string( "hltESPChi2MeasurementEstimator30" ),
0151   nSigma = cms.double( 3.0 ),
0152   MaxSagitta = cms.double( -1.0 ),
0153   MaxChi2 = cms.double( 30.0 ),
0154   MinPtForHitRecoveryInGluedDet = cms.double( 1000000.0 )
0155 )
0156
0157 from TrackingTools.TrackFitters.KFTrajectoryFitter_cfi import KFTrajectoryFitter
0158 hltESPTrajectoryFitterRK = KFTrajectoryFitter.clone(
0159   minHits = cms.int32( 3 ),
0160   ComponentName = cms.string( "hltESPTrajectoryFitterRK" ),
0161   Estimator = cms.string( "hltESPChi2MeasurementEstimator30" ),
0162   Updator = cms.string( "hltESPKFUpdator" ),
0163   Propagator = cms.string( "hltESPRungeKuttaTrackerPropagator" ),
0164   RecoGeometry = cms.string( "hltESPDummyDetLayerGeometry" )
0165 )
0166
0167 hltSiStripExcludedFEDListProducer = cms.EDProducer( "SiStripExcludedFEDListProducer",
0168     ProductLabel = cms.InputTag( "rawDataCollector" )
0169 )
0170
0171 hltMeasurementTrackerEvent = cms.EDProducer( "MeasurementTrackerEventProducer",
0172     inactivePixelDetectorLabels = cms.VInputTag(  ),
0173     stripClusterProducer = cms.string( "hltSiStripRawToClustersFacility" ),
0174     pixelClusterProducer = cms.string( "hltSiPixelClusters" ),
0175     switchOffPixelsIfEmpty = cms.bool( True ),
0176     inactiveStripDetectorLabels = cms.VInputTag( ),
0177     skipClusters = cms.InputTag( "" ),
0178     measurementTracker = cms.string( "hltESPMeasurementTracker" )
0179 )
0180
0181 from Configuration.Eras.Modifier_pp_on_PbPb_run3_cff import pp_on_PbPb_run3
0182 pp_on_PbPb_run3.toModify(hltMeasurementTrackerEvent,
0183                          stripClusterProducer = cms.string( "hltHITrackingSiStripRawToClustersFacilityFullZeroSuppression" ),
0184                          pixelClusterProducer = cms.string( "hltSiPixelClustersAfterSplittingPPOnAA" ),
0185                          )
0186
0187 #####
0188 hltESPTrajectoryFitterRK = cms.ESProducer( "KFTrajectoryFitterESProducer",
0189   appendToDataLabel = cms.string( "" ),
0190   minHits = cms.int32( 3 ),
0191   ComponentName = cms.string( "hltESPTrajectoryFitterRK" ),
0192   Estimator = cms.string( "hltESPChi2MeasurementEstimator30" ),
0193   Updator = cms.string( "hltESPKFUpdator" ),
0194   Propagator = cms.string( "hltESPRungeKuttaTrackerPropagator" ),
0195   RecoGeometry = cms.string( "hltESPDummyDetLayerGeometry" )
0196 )
0197 hltESPTrajectorySmootherRK = cms.ESProducer( "KFTrajectorySmootherESProducer",
0198   errorRescaling = cms.double( 100.0 ),
0199   minHits = cms.int32( 3 ),
0200   ComponentName = cms.string( "hltESPTrajectorySmootherRK" ),
0201   appendToDataLabel = cms.string( "" ),
0202   Estimator = cms.string( "hltESPChi2MeasurementEstimator30" ),
0203   Updator = cms.string( "hltESPKFUpdator" ),
0204   Propagator = cms.string( "hltESPRungeKuttaTrackerPropagator" ),
0205   RecoGeometry = cms.string( "hltESPDummyDetLayerGeometry" )
0206 )
0207 hltESPFittingSmootherIT = cms.ESProducer( "KFFittingSmootherESProducer",
0208   EstimateCut = cms.double( -1.0 ),
0209   appendToDataLabel = cms.string( "" ),
0210   LogPixelProbabilityCut = cms.double( -16.0 ),
0211   MinDof = cms.int32( 2 ),
0212   NoOutliersBeginEnd = cms.bool( False ),
0213   Fitter = cms.string( "hltESPTrajectoryFitterRK" ),
0214   MinNumberOfHits = cms.int32( 3 ),
0215   Smoother = cms.string( "hltESPTrajectorySmootherRK" ),
0216   MaxNumberOfOutliers = cms.int32( 3 ),
0217   BreakTrajWith2ConsecutiveMissing = cms.bool( True ),
0218   MaxFractionOutliers = cms.double( 0.3 ),
0219   NoInvalidHitsBeginEnd = cms.bool( True ),
0220   ComponentName = cms.string( "hltESPFittingSmootherIT" ),
0221   RejectTracks = cms.bool( True )
0222 )
0223
0224 from DQMOffline.Trigger.SiStrip_OfflineMonitoring_cff import *
0225 hltTrackRefitterForSiStripMonitorTrack.TTRHBuilder             = 'hltESPTTRHBWithTrackAngle'
0226 hltTrackRefitterForSiStripMonitorTrack.Propagator              = 'hltESPRungeKuttaTrackerPropagator'
0227 hltTrackRefitterForSiStripMonitorTrack.Fitter                  = 'hltESPFittingSmootherIT'
0228 hltTrackRefitterForSiStripMonitorTrack.MeasurementTrackerEvent = 'hltMeasurementTrackerEvent'
0229 hltTrackRefitterForSiStripMonitorTrack.NavigationSchool        = 'navigationSchoolESProducer'
0230 hltTrackRefitterForSiStripMonitorTrack.src                     = 'hltMergedTracks' # hltIter2Merged
0231
0232 pp_on_PbPb_run3.toModify(hltTrackRefitterForSiStripMonitorTrack,
0233                          src = 'hltMergedTracksPPOnAA')
0234
0235 HLTSiStripMonitorTrack.TopFolderName = 'HLT/SiStrip'
0236 HLTSiStripMonitorTrack.TrackProducer = 'hltTrackRefitterForSiStripMonitorTrack'
0237 HLTSiStripMonitorTrack.TrackLabel    = ''
0238 HLTSiStripMonitorTrack.Cluster_src   = 'hltSiStripRawToClustersFacility'
0239 HLTSiStripMonitorTrack.AlgoName      = 'HLT'
0240 HLTSiStripMonitorTrack.Trend_On      = True
0241 HLTSiStripMonitorTrack.Mod_On        = False
0242 HLTSiStripMonitorTrack.OffHisto_On   = True
0243 HLTSiStripMonitorTrack.HistoFlag_On  = False
0244 HLTSiStripMonitorTrack.TkHistoMap_On = False
0245
0246 pp_on_PbPb_run3.toModify(HLTSiStripMonitorTrack,
0247                          Cluster_src = "hltHITrackingSiStripRawToClustersFacilityFullZeroSuppression")
0248
0249 pp_on_PbPb_run3.toModify(HLTSiStripMonitorCluster,
0250                          BPTXfilter = dict(l1Algorithms = ['L1_ZeroBias']))
0251
0252 HLTSiStripMonitorClusterAPVgainCalibration = HLTSiStripMonitorCluster.clone()
0253 from DQM.TrackingMonitorSource.pset4GenericTriggerEventFlag_cfi import *
0254 #HLTSiStripMonitorClusterAPVgainCalibration.BPTXfilter = genericTriggerEventFlag4fullTrackerAndHLTnoHIPnoOOTdb # HLT_ZeroBias_FirstCollisionAfterAbortGap_*
0255 HLTSiStripMonitorClusterAPVgainCalibration.BPTXfilter = cms.PSet(
0256    andOr         = cms.bool( False ),
0257 ### DCS selection
0258    dcsInputTag   = cms.InputTag( "scalersRawToDigi" ),
0259    dcsRecordInputTag   = cms.InputTag( "onlineMetaDataDigis" ),
0260    dcsPartitions = cms.vint32 ( 24, 25, 26, 27, 28, 29 ),
0261    andOrDcs      = cms.bool( False ),
0262    errorReplyDcs = cms.bool( True ),
0263 ### HLT selection
0264    andOrHlt      = cms.bool(True),# True:=OR; False:=AND
0265    hltInputTag   = cms.InputTag( "TriggerResults::HLT" ),
0266    hltPaths      = cms.vstring("HLT_ZeroBias_FirstCollisionAfterAbortGap_v*"),
0267    errorReplyHlt = cms.bool( False ),
0268 #   errorReplyHlt = cms.bool( True ),
0269 ### L1 selection
0270 #   andOrL1       = cms.bool( True ),
0271 #   l1Algorithms  = cms.vstring("L1_ZeroBias_FirstCollidingBunch"),
0272 ##   l1BeforeMask  = cms.bool( True ), # specifies, if the L1 algorithm decision should be read as before (true) or after (false) masking is applied.
0273 #   l1BeforeMask  = cms.bool( False ), # specifies, if the L1 algorithm decision should be read as before (true) or after (false) masking is applied.
0274 #   errorReplyL1  = cms.bool( False ),
0275    verbosityLevel = cms.uint32(1)
0276 )
0277 HLTSiStripMonitorClusterAPVgainCalibration.TopFolderName = cms.string('HLT/SiStrip/ZeroBias_FirstCollisionAfterAbortGap')
0278
0279 pp_on_PbPb_run3.toModify(HLTSiStripMonitorClusterAPVgainCalibration,
0280                          BPTXfilter = dict(hltPaths = ["HLT_HICentrality30100_FirstCollisionAfterAbortGap_v*"]),
0281                          TopFolderName = cms.string('HLT/SiStrip/HLT_HICentrality30100_FirstCollisionAfterAbortGap'))
0282
0283 sistripOnlineMonitorHLTsequence = cms.Sequence(
0284     hltMeasurementTrackerEvent
0285     * sistripMonitorHLTsequence # strip cluster monitoring
0286     * HLTSiStripMonitorClusterAPVgainCalibration
0287 )