dump_header can dump #counts for each detector now also for mud and nexus.

src/external/nexus/PNeXus.cpp

@@ -635,6 +635,21 @@ unsigned int PNeXusData1::GetHistoLength(unsigned int histoNo)
   return fHisto[histoNo].size();
 }
 
+//------------------------------------------------------------------------------------------
+// GetHistoLength (public)
+//------------------------------------------------------------------------------------------
+unsigned int PNeXusData1::GetHistoCounts(unsigned int histoNo)
+{
+  if (histoNo >= fHisto.size())
+    return 0;
+
+  unsigned int counts=0;
+  for (unsigned int i=0; i<fHisto[histoNo].size(); i++)
+    counts += fHisto[histoNo][i];
+
+  return counts;
+}
+
 //------------------------------------------------------------------------------------------
 // GetHisto (public)
 //------------------------------------------------------------------------------------------
@@ -1258,6 +1273,36 @@ int PNeXusDetector2::SetLastGoodBin(int *lgb)
   return result;
 }
 
+//------------------------------------------------------------------------------------------
+// GetHistoCounts (public)
+//------------------------------------------------------------------------------------------
+/**
+ * <p>Get number of counts in a given histogram histoNo
+ *
+ * \param idx_p period index
+ * \param idx_s spectrum index
+ *
+ * \return the number of entries in the selected histogram histoNo, or 0 if the histoNo is out of scope.
+ */
+unsigned int PNeXusDetector2::GetHistoCounts(int idx_p, int idx_s)
+{
+   unsigned counts = 0;
+
+   if (idx_p > 0)
+     if (idx_p > fNoOfPeriods)
+       return counts;
+
+   if (idx_s > 0)
+     if (idx_s > fNoOfSpectra)
+       return counts;
+
+   for (int i=0; i<fNoOfBins; i++) {
+     counts += GetHistoValue(idx_p,idx_s,i);
+   }
+
+   return counts;
+}
+
 //------------------------------------------------------------------------------------------
 // GetHistoValue (public)
 //------------------------------------------------------------------------------------------
@@ -2151,8 +2196,9 @@ int PNeXus::WriteFile(const char *fileName, const char *fileType, const unsigned
 /**
  * <p>Write the content of the NeXus file to stdout. Used for debugging purposes.
  *
+ * \param counts flag, if true, also dump the counts for all the histograms
  */
-void PNeXus::Dump()
+void PNeXus::Dump(const bool counts)
 {
   double dval;
   std::string str;
@@ -2259,6 +2305,8 @@ void PNeXus::Dump()
           std::cout << fNxEntry1->GetData()->GetHisto(i)->at(j) << ", ";
         }
         std::cout << "...";
+        if (counts)
+          std::cout << " total no of entries: " << fNxEntry1->GetData()->GetHistoCounts(i);
       }
     }
     if (fNxEntry1->GetData()->GetAlpha()->size() == 0) {
@@ -2408,6 +2456,8 @@ void PNeXus::Dump()
             std::cout << fNxEntry2->GetInstrument()->GetDetector()->GetHistoValue(i,j,k) << ", ";
           }
           std::cout << "...";
+          if (counts)
+            std::cout << " total number of entries: " << fNxEntry2->GetInstrument()->GetDetector()->GetHistoCounts(i,j);
         }
       }
     } else {
@@ -2421,6 +2471,8 @@ void PNeXus::Dump()
             std::cout << fNxEntry2->GetInstrument()->GetDetector()->GetHistoValue(0,j,k) << ", ";
           }
           std::cout << "...";
+          if (counts)
+            std::cout << " total number of entries: " << fNxEntry2->GetInstrument()->GetDetector()->GetHistoCounts(0,j);
         }
       } else {  // counts[ntc]
         std::cout << std::endl << "      (#bins=" << fNxEntry2->GetInstrument()->GetDetector()->GetNoOfBins() << ")";
@@ -2431,6 +2483,8 @@ void PNeXus::Dump()
           std::cout << fNxEntry2->GetInstrument()->GetDetector()->GetHistoValue(0,0,k) << ", ";
         }
         std::cout << "...";
+        if (counts)
+          std::cout << " total number of entries: " << fNxEntry2->GetInstrument()->GetDetector()->GetHistoCounts(0,0);
       }
     }
     std::cout << std::endl << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++";

src/external/nexus/PNeXus.h

@@ -254,6 +254,7 @@ class PNeXusData1 {
     virtual void GetHistoName(unsigned int idx, std::string &name, bool &ok);
     virtual unsigned int GetNoOfHistos() { return fHisto.size(); }
     virtual unsigned int GetHistoLength(unsigned int histoNo=0);
+    virtual unsigned int GetHistoCounts(unsigned int histoNo=0);
     virtual std::vector<unsigned int> *GetHisto(unsigned int histoNo);
     virtual std::vector<int> *GetGrouping() { return &fGrouping; }
     virtual std::vector<PNeXusAlpha1> *GetAlpha() { return &fAlpha; }
@@ -399,6 +400,7 @@ class PNeXusDetector2 {
     virtual int GetNoOfPeriods() { return fNoOfPeriods; }
     virtual int GetNoOfSpectra() { return fNoOfSpectra; }
     virtual int GetNoOfBins() { return fNoOfBins; }
+    virtual unsigned int GetHistoCounts(int idx_p, int idx_s);
     virtual int GetHistoValue(int idx_p, int idx_s, int idx_b);
     virtual int* GetHistos() { return fHisto; }
     virtual unsigned int GetSpectrumIndexSize() { return fSpectrumIndex.size(); }
@@ -586,7 +588,7 @@ class PNeXus {
 
     virtual void SetCreator(std::string str) { fCreator = str; }
 
-    virtual void Dump();
+    virtual void Dump(const bool counts);
 
   private:
     bool fValid;