Adapt TrackToTrackComparisonHists to run in cosmics mode

Author: mmusich

DQM/TrackingMonitorSource/interface/TrackToTrackComparisonHists.h

@@ -103,6 +103,7 @@ class TrackToTrackComparisonHists : public DQMEDAnalyzer {
 
 private:
   //  edm::ParameterSet conf_;
+  const bool isCosmics_;
   std::string topDirName_;
   double dRmin_;
   double pTCutForPlateau_;

DQM/TrackingMonitorSource/src/TrackToTrackComparisonHists.cc

@@ -17,6 +17,7 @@
 TrackToTrackComparisonHists::TrackToTrackComparisonHists(const edm::ParameterSet& iConfig)
     : monitoredTrackInputTag_(iConfig.getParameter<edm::InputTag>("monitoredTrack")),
       referenceTrackInputTag_(iConfig.getParameter<edm::InputTag>("referenceTrack")),
+      isCosmics_(iConfig.getParameter<bool>("isCosmics")),
       topDirName_(iConfig.getParameter<std::string>("topDirName")),
       dRmin_(iConfig.getParameter<double>("dRmin")),
       pTCutForPlateau_(iConfig.getParameter<double>("pTCutForPlateau")),
@@ -119,15 +120,17 @@ void TrackToTrackComparisonHists::analyze(const edm::Event& iEvent, const edm::E
 
   edm::Handle<reco::VertexCollection> referencePVHandle;
   iEvent.getByToken(referencePVToken_, referencePVHandle);
-  if (!referencePVHandle.isValid()) {
+  if (!referencePVHandle.isValid() && !isCosmics_) {
     edm::LogError("TrackToTrackComparisonHists") << "referencePVHandle not found, skipping event";
     return;
   }
-  if (referencePVHandle->empty()) {
+  if (referencePVHandle->empty() && !isCosmics_) {
     edm::LogInfo("TrackToTrackComparisonHists") << "referencePVHandle->size is 0 ";
     return;
   }
-  reco::Vertex referencePV = referencePVHandle->at(0);
+  reco::Vertex referencePV = isCosmics_
+                                 ? reco::Vertex(referenceBS.position(), referenceBS.rotatedCovariance3D(), 0., 0., 0)
+                                 : referencePVHandle->at(0);
 
   //
   //  Get Monitored Track Info
@@ -151,15 +154,18 @@ void TrackToTrackComparisonHists::analyze(const edm::Event& iEvent, const edm::E
 
   edm::Handle<reco::VertexCollection> monitoredPVHandle;
   iEvent.getByToken(monitoredPVToken_, monitoredPVHandle);
-  if (!monitoredPVHandle.isValid()) {
-    edm::LogError("TrackToTrackComparisonHists") << "monitoredPVHandle not found, skipping event";
+  if (!monitoredPVHandle.isValid() && !isCosmics_) {
+    edm::LogError("TrackToTrackComparisonHists")
+        << "monitoredPVHandle not found, skipping event isCosmics value:" << isCosmics_;
     return;
   }
-  if (monitoredPVHandle->empty()) {
+  if (monitoredPVHandle->empty() && !isCosmics_) {
     edm::LogInfo("TrackToTrackComparisonHists") << "monitoredPVHandle->size is 0 ";
     return;
   }
-  reco::Vertex monitoredPV = monitoredPVHandle->at(0);
+  reco::Vertex monitoredPV = isCosmics_
+                                 ? reco::Vertex(monitoredBS.position(), monitoredBS.rotatedCovariance3D(), 0., 0., 0)
+                                 : monitoredPVHandle->at(0);
 
   edm::LogInfo("TrackToTrackComparisonHists")
       << "analyzing " << monitoredTrackInputTag_.process() << ":" << monitoredTrackInputTag_.label() << ":"
@@ -306,6 +312,7 @@ void TrackToTrackComparisonHists::fillDescriptions(edm::ConfigurationDescription
   edm::ParameterSetDescription desc;
 
   desc.add<bool>("requireValidHLTPaths", true);
+  desc.add<bool>("isCosmics", false);
 
   desc.add<edm::InputTag>("monitoredTrack", edm::InputTag("hltMergedTracks"));
   desc.add<edm::InputTag>("monitoredBeamSpot", edm::InputTag("hltOnlineBeamSpot"));

DQMOffline/Trigger/python/TrackToTrackMonitoringCosmics_cff.py

@@ -11,7 +11,17 @@
     referenceBeamSpot        = "offlineBeamSpot",
     topDirName               = "HLT/Tracking/ValidationWRTOffline/hltCtfWithMaterialTracksP5",
     referencePrimaryVertices = "offlinePrimaryVertices",
-    monitoredPrimaryVertices = "hltVerticesPFSelector"
+    monitoredPrimaryVertices = "hltPixelVertices",
+    isCosmics                = cms.bool(True),
+    dxyCutForPlateau         = 1e6,
+    histoPSet                = dict(
+        Dxy_rangeMin = -60,
+        Dxy_rangeMax = 60,
+        Dxy_nbin = 120,
+        Dz_rangeMin = -250,
+        Dz_rangeMax =  250,
+        Dz_nbin = 250,
+    ) 
 )
 
 hltToOfflineCosmicsTrackValidatorSequence = cms.Sequence(hltCtfWithMaterialTracksP5_2_ctfWithMaterialTracksP5)