Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-09-28 23:48:26

0001 import FWCore.ParameterSet.Config as cms
0002 from DQMServices.Core.DQMEDHarvester import DQMEDHarvester
0003 from DQM.SiPixelPhase1Common.HistogramManager_cfi import *
0004 import DQM.SiPixelPhase1Common.TriggerEventFlag_cfi as trigger
0005 
0006 SiPixelPhase1TrackResidualsResidualsX = DefaultHistoTrack.clone(
0007   name = "residual_x",
0008   title = "Track Residuals X",
0009   range_min = -0.1, range_max = 0.1, range_nbins = 100,
0010   xlabel = "(x_rec - x_pred) [cm]",
0011   dimensions = 1,
0012   specs = VPSet(
0013     StandardSpecification2DProfile,
0014     Specification().groupBy("PXBarrel/PXLayer").saveAll(),
0015     Specification().groupBy("PXForward/PXDisk").saveAll(),
0016     
0017     Specification().groupBy("PXBarrel/PXLayer/LumiBlock")
0018                    .reduce("MEAN")
0019                    .groupBy("PXBarrel/PXLayer", "EXTEND_X")
0020                    .save(),
0021 
0022     Specification().groupBy("PXForward/PXDisk/LumiBlock")
0023                    .reduce("MEAN")
0024                    .groupBy("PXForward/PXDisk", "EXTEND_X")
0025                    .save(),
0026 
0027     Specification(PerLayer1D).groupBy("PXBarrel/Shell/PXLayer").save(),
0028     Specification(PerLayer1D).groupBy("PXForward/HalfCylinder/PXRing/PXDisk").save()
0029   )
0030 )
0031 
0032 SiPixelPhase1TrackResidualsResidualsY = SiPixelPhase1TrackResidualsResidualsX.clone(
0033   name = "residual_y",
0034   title = "Track Residuals Y",
0035   xlabel = "(y_rec - y_pred) [cm]",
0036 )
0037 
0038 SiPixelPhase1TrackResidualsResOnEdgeX = DefaultHistoTrack.clone(
0039   name = "residual_OnEdge_x",
0040   title = "Track Residuals X (OnEdge Clusters)",
0041   range_min = -0.1, range_max = 0.1, range_nbins = 100,
0042   xlabel = "(x_rec - x_pred) [cm]",
0043   dimensions = 1,
0044   specs = VPSet(
0045     Specification().groupBy("PXBarrel/PXLayer").saveAll(),
0046     Specification().groupBy("PXForward/PXDisk").saveAll(),
0047     Specification(PerLayer1D).groupBy("PXBarrel/Shell/PXLayer").save(),
0048     Specification(PerLayer1D).groupBy("PXForward/HalfCylinder/PXRing/PXDisk").save()
0049   )
0050 )
0051 
0052 SiPixelPhase1TrackResidualsResOnEdgeY = SiPixelPhase1TrackResidualsResOnEdgeX.clone(
0053   name = "residual_OnEdge_y",
0054   title = "Track Residuals Y (OnEdge Clusters)",
0055   xlabel = "(y_rec - y_pred) [cm]",
0056 )
0057 
0058 
0059 SiPixelPhase1TrackResidualsResOtherBadX = DefaultHistoTrack.clone(
0060   name = "residual_OtherBad_x",
0061   title = "Track Residuals X (OtherBad Clusters)",
0062   range_min = -0.1, range_max = 0.1, range_nbins = 100,
0063   xlabel = "(x_rec - x_pred) [cm]",
0064   dimensions = 1,
0065   specs = VPSet(
0066     Specification().groupBy("PXBarrel/PXLayer").saveAll(),
0067     Specification().groupBy("PXForward/PXDisk").saveAll(),
0068     Specification(PerLayer1D).groupBy("PXBarrel/Shell/PXLayer").save(),
0069     Specification(PerLayer1D).groupBy("PXForward/HalfCylinder/PXRing/PXDisk").save()
0070   )
0071 )
0072 
0073 SiPixelPhase1TrackResidualsResOtherBadY = SiPixelPhase1TrackResidualsResOtherBadX.clone(
0074   name = "residual_OtherBad_y",
0075   title = "Track Residuals Y (OtherBad Clusters)",
0076   xlabel = "(y_rec - y_pred) [cm]",
0077 )
0078 
0079 
0080 SiPixelPhase1TrackNormResX = DefaultHistoTrack.clone(
0081   topFolderName= cms.string("PixelPhase1/Tracks/ResidualsExtra"), 
0082   name = "NormRes_x",
0083   title = "Normalized Residuals X",
0084   range_min = -5, range_max = 5, range_nbins = 100,
0085   xlabel = "(x_rec - x_pred)/x_err",
0086   dimensions = 1,
0087   specs = VPSet(
0088     Specification().groupBy("PXBarrel/PXLayer").saveAll(),
0089     Specification().groupBy("PXForward/PXDisk").saveAll(),
0090     Specification(PerLayer1D).groupBy("PXBarrel/Shell/PXLayer").save(),
0091     Specification(PerLayer1D).groupBy("PXForward/HalfCylinder/PXRing/PXDisk").save()
0092   )
0093 )
0094 
0095 SiPixelPhase1TrackNormResY = SiPixelPhase1TrackNormResX.clone(
0096   name = "NormRes_y",
0097   title = "Normalized Residuals Y",
0098   range_min = -5, range_max = 5, range_nbins = 100,
0099   xlabel = "(y_rec - y_pred)/y_err",
0100 )
0101 
0102 SiPixelPhase1TrackDRnRX = DefaultHistoTrack.clone(
0103   topFolderName= cms.string("PixelPhase1/Tracks/ResidualsExtra"),
0104   name = "DRnR_x",
0105   title = "Distribution of RMS of Normalized Residuals X",
0106   range_min = -5, range_max = 5, range_nbins = 100,
0107   xlabel = "#sigma_{(x_rec - x_pred)/x_err}",
0108   dimensions = 1,
0109   specs = VPSet(
0110     StandardSpecification2DProfile
0111   )
0112 )
0113 
0114 SiPixelPhase1TrackDRnRY = SiPixelPhase1TrackDRnRX.clone(
0115   name = "DRnR_y",
0116   title = "Distribution of RMS of Normalized Residuals Y",
0117   range_min = -5, range_max = 5, range_nbins = 100,
0118   xlabel = "#sigma_{(y_rec - y_pred)/y_err}",
0119 )
0120 
0121 
0122 SiPixelPhase1TrackResidualsConf = cms.VPSet(
0123   SiPixelPhase1TrackResidualsResidualsX,
0124   SiPixelPhase1TrackResidualsResidualsY,
0125   SiPixelPhase1TrackResidualsResOnEdgeX,
0126   SiPixelPhase1TrackResidualsResOnEdgeY,
0127   SiPixelPhase1TrackResidualsResOtherBadX,
0128   SiPixelPhase1TrackResidualsResOtherBadY,
0129   SiPixelPhase1TrackNormResX,
0130   SiPixelPhase1TrackNormResY,
0131   SiPixelPhase1TrackDRnRX,
0132   SiPixelPhase1TrackDRnRY
0133 )
0134 
0135 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0136 SiPixelPhase1TrackResidualsAnalyzer = DQMEDAnalyzer('SiPixelPhase1TrackResiduals',
0137         trajectoryInput = cms.string("generalTracks"),
0138         Tracks        = cms.InputTag("generalTracks"),
0139         vertices = cms.InputTag("offlinePrimaryVertices"),
0140         histograms = SiPixelPhase1TrackResidualsConf,
0141         geometry = SiPixelPhase1Geometry,
0142         triggerflags = trigger.SiPixelPhase1Triggers,
0143         VertexCut = cms.untracked.bool(True)
0144 )
0145 
0146 SiPixelPhase1TrackResidualsHarvester = DQMEDHarvester("SiPixelPhase1Harvester",
0147         histograms = SiPixelPhase1TrackResidualsConf,
0148         geometry = SiPixelPhase1Geometry
0149 )