Merge pull request #35507 from CMSTrackerDPG/SiPixelLorentzAngleCalib…

…rationInPCL Prompt calibration loop module for SiPixelLorentzAngle
cms-sw · Oct 14, 2021 · 2a970a2 · 2a970a2
2 parents 0708a9b + a9c114b
commit 2a970a2
Show file tree

Hide file tree

Showing 19 changed files with 1,499 additions and 13 deletions.
diff --git a/CalibTracker/SiPixelLorentzAngle/BuildFile.xml b/CalibTracker/SiPixelLorentzAngle/BuildFile.xml
@@ -8,5 +8,9 @@
 <use name="Geometry/Records"/>
 <use name="SimTracker/TrackerHitAssociation"/>
 <use name="MagneticField/Records"/>
+<use name="DQMServices/Core"/>
+<use name="CondFormats/SiPixelObjects"/>
+<use name="CondFormats/SiPixelTransient"/>
+<use name="CondCore/DBOutputService"/>
 <use name="root"/>
 <use name="clhep"/>
diff --git a/CalibTracker/SiPixelLorentzAngle/interface/SiPixelLorentzAngleCalibrationStruct.h b/CalibTracker/SiPixelLorentzAngle/interface/SiPixelLorentzAngleCalibrationStruct.h
@@ -0,0 +1,35 @@
+#ifndef CalibTracker_SiPixelLorentzAngle_SiPixelLorentzAngleCalibrationStruct_h
+#define CalibTracker_SiPixelLorentzAngle_SiPixelLorentzAngleCalibrationStruct_h
+
+#include "DQMServices/Core/interface/DQMStore.h"
+#include <unordered_map>
+
+struct SiPixelLorentzAngleCalibrationHistograms {
+public:
+  SiPixelLorentzAngleCalibrationHistograms() = default;
+
+  using MonitorMap = std::unordered_map<uint32_t, dqm::reco::MonitorElement*>;
+
+  int nlay;
+  std::vector<int> nModules_;
+  std::vector<std::string> BPixnewmodulename_;
+  std::vector<int> BPixnewDetIds_;
+  std::vector<int> BPixnewModule_;
+  std::vector<int> BPixnewLayer_;
+
+  std::vector<std::string> FPixnewmodulename_;
+  std::vector<int> FPixnewDetIds_;
+  std::vector<int> FPixnewDisk_;
+  std::vector<int> FPixnewBlade_;
+  std::unordered_map<uint32_t, std::vector<uint32_t> > detIdsList;
+
+  MonitorMap h_drift_depth_adc_;
+  MonitorMap h_drift_depth_adc2_;
+  MonitorMap h_drift_depth_noadc_;
+  MonitorMap h_drift_depth_;
+  MonitorMap h_mean_;
+
+  dqm::reco::MonitorElement* h_tracks_;
+};
+
+#endif
diff --git a/CalibTracker/SiPixelLorentzAngle/python/SiPixelLorentzAnglePCLHarvester_cfi.py b/CalibTracker/SiPixelLorentzAngle/python/SiPixelLorentzAnglePCLHarvester_cfi.py
@@ -0,0 +1,11 @@
+import FWCore.ParameterSet.Config as cms
+from DQMServices.Core.DQMEDHarvester import DQMEDHarvester
+from CalibTracker.SiPixelLorentzAngle.SiPixelLorentzAnglePCLWorker_cfi import SiPixelLorentzAnglePCLWorker as worker
+
+SiPixelLorentzAnglePCLHarvester = DQMEDHarvester(
+    "SiPixelLorentzAnglePCLHarvester",
+    newmodulelist = cms.vstring(worker.newmodulelist.value()), # taken from worker configuration, need to stay in synch
+    dqmDir = cms.string(worker.folder.value()), # taken from worker configuration, need to stay in synch
+    record = cms.string("SiPixelLorentzAngleRcd"),
+    fitProbCut = cms.double(0.5)
+)
diff --git a/CalibTracker/SiPixelLorentzAngle/python/SiPixelLorentzAnglePCLWorker_cfi.py b/CalibTracker/SiPixelLorentzAngle/python/SiPixelLorentzAnglePCLWorker_cfi.py
@@ -0,0 +1,27 @@
+import FWCore.ParameterSet.Config as cms
+
+from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
+SiPixelLorentzAnglePCLWorker = DQMEDAnalyzer(
+    "SiPixelLorentzAnglePCLWorker",
+    folder = cms.string('AlCaReco/SiPixelLorentzAngle'),
+    notInPCL = cms.bool(False),
+    fileName = cms.string('testrun.root'),
+    newmodulelist = cms.vstring("BPix_BmI_SEC7_LYR2_LDR12F_MOD1",
+                                "BPix_BmI_SEC8_LYR2_LDR14F_MOD1",
+                                "BPix_BmO_SEC3_LYR2_LDR5F_MOD1",
+                                "BPix_BmO_SEC3_LYR2_LDR5F_MOD2",
+                                "BPix_BmO_SEC3_LYR2_LDR5F_MOD3",
+                                "BPix_BpO_SEC1_LYR2_LDR1F_MOD1",
+                                "BPix_BpO_SEC1_LYR2_LDR1F_MOD2",
+                                "BPix_BpO_SEC1_LYR2_LDR1F_MOD3"),
+    src = cms.InputTag("TrackRefitter"),
+    binsDepth    = cms.int32(50),
+    binsDrift =    cms.int32(200),
+    ptMin = cms.double(3),
+    normChi2Max = cms.double(2),
+    clustSizeYMin = cms.int32(4),
+    clustSizeYMinL4 = cms.int32(3),
+    clustSizeXMax = cms.int32(5),
+    residualMax = cms.double(0.005),
+    clustChargeMaxPerLength = cms.double(50000)
+)