From 6ac8f429c5ab6cff4455e3dd83f02ab4cf54b880 Mon Sep 17 00:00:00 2001
From: Mark Sutton <mark.sutton@cern.ch>
Date: Thu, 19 Nov 2020 11:04:45 +0000
Subject: [PATCH] Add 2D histograms for the showing correlation between the
 numbers of pixel and sct hits

In the TrigInDteAnalysis it is useful to have such a histogram
---
 .../Analysis/src/ConfAnalysis.cxx                      | 10 ++++++++++
 .../TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.h  |  5 ++++-
 .../share/TIDAdata_cuts-offline.dat                    |  3 ++-
 .../TrigInDetAnalysisUser/share/TIDAhisto-panel.dat    |  8 ++++----
 4 files changed, 20 insertions(+), 6 deletions(-)

diff --git a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.cxx b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.cxx
index e8947c581b8..aee47ba68d5 100644
--- a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.cxx
+++ b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.cxx
@@ -245,6 +245,9 @@ void ConfAnalysis::initialiseInternal() {
   mres.push_back(  rnscth_pt_rec = new Resplot( "nscth_pt_rec", ptnbins, ptbinlims,  22, -0.5, 21.5 ) );
 
 
+  mres.push_back(  m_rnsct_vs_npix     = new Resplot( "nsct_vs_npix",     12, -0.5,  11.5,   22, -0.5, 21.5 ) );
+  mres.push_back(  m_rnsct_vs_npix_rec = new Resplot( "nsct_vs_npix_rec", 12, -0.5,  11.5,   22, -0.5, 21.5 ) );
+
   mres.push_back(  rChi2prob = new Resplot( "Chi2prob", ptnbins, ptbinlims,  20, 0,   1 ) );
   mres.push_back(  rChi2     = new Resplot( "Chi2",     ptnbins, ptbinlims, 200, 0, 100 ) );
   mres.push_back(  rChi2dof  = new Resplot( "Chi2dof",  ptnbins, ptbinlims, 100, 0,  10 ) );
@@ -1378,6 +1381,9 @@ void ConfAnalysis::execute(const std::vector<TIDA::Track*>& reftracks,
 
     rnpix_lb->Fill( gevent->lumi_block(), npixt*1.0 );
 
+
+    m_rnsct_vs_npix->Fill( npixt, nsctt );
+
     double                 etovpt_val = 0;
     const TrackTrigObject* tobj       = 0;
     
@@ -1917,6 +1923,10 @@ void ConfAnalysis::execute(const std::vector<TIDA::Track*>& reftracks,
     double ntrtr   = testtracks[i]->trHits(); 
     double nstrawr = testtracks[i]->strawHits(); 
 
+
+    m_rnsct_vs_npix_rec->Fill( npixr, nsctr );
+
+
 #if 0
     double dpTr_b  = testtracks[i]->dpT()/1000;
     double detar_b = testtracks[i]->deta();
diff --git a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.h b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.h
index 6cb614022c2..798ec0a275e 100644
--- a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.h
+++ b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/ConfAnalysis.h
@@ -328,7 +328,10 @@ private:
   bool m_initialised;
   bool m_initialiseFirstEvent;
 
-  //  bool m_lfirst;
+
+  Resplot* m_rnsct_vs_npix;
+  Resplot* m_rnsct_vs_npix_rec;
+
 
 };
 
diff --git a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAdata_cuts-offline.dat b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAdata_cuts-offline.dat
index 4877a21535d..82f449ead7d 100644
--- a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAdata_cuts-offline.dat
+++ b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAdata_cuts-offline.dat
@@ -20,8 +20,9 @@ Rmatch = 0.05;
 
 dumpflag = 1;
 
+
 a0v = 10.5;
 z0v = 300.5;
 
 VertexSelection    = 0;
-VertexSelectionRec = 0;
\ No newline at end of file
+VertexSelectionRec = 0;
diff --git a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAhisto-panel.dat b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAhisto-panel.dat
index 7aa47db4793..df08b3920fa 100644
--- a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAhisto-panel.dat
+++ b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/share/TIDAhisto-panel.dat
@@ -31,11 +31,11 @@ diff_panel = {
      "ript_vs_eta/sigma",        "Residual 1/p_{T} vs #eta",     "xaxis:lin:auto",     "Offline #eta",         "yaxis:log:auto",  "1/p_{T} resolution [GeV^{-1}]",
      "ript_vs_pt/sigma",         "Residual 1/p_{T} vs p_{T}",    "xaxis:log:auto",     "Offline p_{T} [GeV]",   "yaxis:lin:auto", "1/p_{T} resolution [GeV^{-1}]",
 
-     "rzed_vs_eta/sigma",        "Residual z vs #eta",           "xaxis:lin",          "Offline #eta",          "yaxis:lin:auto", "z_{0} resolution [mm]",
-     "rzed_vs_pt/sigma",         "Residual z vs p_{T}",          "xaxis:log:auto",     "Offline p_{T} [GeV]",   "yaxis:lin:auto", "z_{0} resolution [mm]",
+     "rzed_vs_eta/sigma",        "Residual z vs #eta",           "xaxis:lin",          "Offline #eta",          "yaxis:lin:0:0.9", "z_{0} resolution [mm]",
+     "rzed_vs_pt/sigma",         "Residual z vs p_{T}",          "xaxis:log:auto",     "Offline p_{T} [GeV]",   "yaxis:lin:0:0.9", "z_{0} resolution [mm]",
 
-     "rd0_vs_eta/sigma",         "Residual d vs #eta",           "xaxis:lin",           "Offline #eta",          "yaxis:lin:auto",  "d_{0} resolution [mm]",                    
-     "rd0_vs_pt/sigma",          "Residual d vs p_{T}",          "xaxis:log:auto",      "Offline p_{T} [GeV]",   "yaxis:lin:auto",  "d_{0} resolution [mm]",
+     "rd0_vs_eta/sigma",         "Residual d vs #eta",           "xaxis:lin",           "Offline #eta",          "yaxis:lin:0:0.12",  "d_{0} resolution [mm]",              
+     "rd0_vs_pt/sigma",          "Residual d vs p_{T}",          "xaxis:log:auto",      "Offline p_{T} [GeV]",   "yaxis:lin:0:0.12",  "d_{0} resolution [mm]",
 
      "rd0_vs_signed_pt/sigma",   "Residual d vs signed p_{T}",   "xaxis:lin:-100:100",  "Offline p_{T} [GeV]",   "yaxis:lin:auto",  "d_{0} resolution [mm]",    
      "rzed_vs_signed_pt/sigma",  "Residual z vs signed p_{T}",   "xaxis:lin:-100:100",  "Offline p_{T} [GeV]",   "yaxis:lin:auto",  "z_{0} resolution [mm]"    
-- 
GitLab