xml/XsilHistogram.hh

/*----------------------------------------------------------------------*/
/*                                                                      */
/* Module Name: xsilHistogram                                           */
/*                                                                      */
/* Module Description: Classes for LIGO-LW input/output to DMT objects. */
/*                                                                      */
/* Revision History:                                                    */
/* Rel   Date     Programmer    Comments                                */
/* 1.0   4/1/02   D. Sigg       Initial revision                        */
/*                                                                      */
/* Documentation References:                                            */
/*      Man Pages: none                                                 */
/*      References: none                                                */
/*                                                                      */
/* Author Information:                                                  */
/* Name          Telephone       Fax             e-mail                 */
/* Daniel Sigg   (509) 372-8132  (509) 372-8137  sigg_d@ligo.mit.edu    */
/*                                                                      */
/*                                                                      */
/*                      -------------------                             */
/*                                                                      */
/*                             LIGO                                     */
/*                                                                      */
/*        THE LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY.      */
/*                                                                      */
/*                     (C) The LIGO Project, 2001.                      */
/*                                                                      */
/*                                                                      */
/* Caltech                              MIT                             */
/* LIGO Project MS 51-33                LIGO Project NW-17 161          */
/* Pasadena CA 91125                    Cambridge MA 01239              */
/*                                                                      */
/* LIGO Hanford Observatory             LIGO Livingston Observatory     */
/* P.O. Box 1970 S9-02                  19100 LIGO Lane Rd.             */
/* Richland WA 99352                    Livingston, LA 70754            */
/*                                                                      */
/*----------------------------------------------------------------------*/

#ifndef _GDS_XSILHISTOGRAM_HH
#define _GDS_XSILHISTOGRAM_HH

#include <string>
#include <iosfwd>
#include "xml/Xsil.hh"
#include "Histogram1.hh"
#include "Histogram2.hh"

namespace xml {

   class xsilHistogram {
   protected:
      int            fLevel;
      const Histogram1* fH1;
      const Histogram2* fH2;
   public:
      xsilHistogram (int level = 1) 
      : fLevel(level), fH1 (0), fH2 (0) {
      }
      xsilHistogram (const Histogram1* h1, int level = 1) 
      : fLevel(level), fH1 (h1), fH2 (0) {              
      }
      xsilHistogram (const Histogram1& h1, int level = 1) 
      : fLevel(level), fH1 (&h1), fH2 (0) {              
      }
      xsilHistogram (const Histogram2* h2, int level = 1) 
      : fLevel(level), fH1 (0), fH2 (h2) {              
      }
      xsilHistogram (const Histogram2& h2, int level = 1) 
      : fLevel(level), fH1 (0), fH2 (&h2) {              
      }
      std::ostream& write (std::ostream &os) const;
   };


   inline std::ostream& operator<< (std::ostream &os, 
                     const xsilHistogram &h) {
      return h.write(os); }


   class xsilHandlerHistogram : public xsilHandler {
   protected:
      std::vector<Histogram1>*    fDat1;
      std::vector<Histogram2>*    fDat2;
   
      std::string                 fTitle;
      std::string                 fLabel[4];
      int                         fSubtype;
      unsigned long               fSec;
      unsigned long               fNsec;
      int                         fNBin[3];
      int                         fNData;
      double                      fSumWeight;
      double                      fSumWeightSqr;
      double                      fSumWeightDim[3];
      double                      fSumWeightDimSqr[3];
      double                      fSumWeightXY;
      double                      fLowEdge[3];
      double                      fSpacing[3];
      double*                     fBins[3];
      int                         fBinsLen[3];
      double*                     fErrors;
      int                         fErrorsDim[3];
      double*                     fData;
      int                         fDataDim[3];
   
   public:
      explicit xsilHandlerHistogram (std::vector<Histogram1>* objs1, 
                        std::vector<Histogram2>* objs2,
                        const attrlist *attr=0, bool ignore=false);
      virtual ~xsilHandlerHistogram();
      virtual bool HandleParameter (const std::string& name,
                        const attrlist& attr,
                        const bool& p, int N = 1) {
         return false; }
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr,
                        const char& p, int N = 1) {
         return false; }
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr,
                        const short& p, int N = 1) {
         return false; }
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr,
                        const int& p, int N=1);
   #ifndef __CINT__
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr,
                        const long long& p, int N=1) {
         return false; }
   #endif //__CINT__
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr, 
                        const float& p, int N=1) {
         return false; }
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr, 
                        const double& p, int N=1);
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr,
                        const std::complex<float>& p, int N = 1) {
         return false; }
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr,
                        const std::complex<double>& p, int N = 1) {
         return false; }
      virtual bool HandleParameter(const std::string& name, 
                        const attrlist& attr,
                        const std::string& p);
      virtual bool HandleTime (const std::string& name, 
                        const attrlist& attr,
                        unsigned long sec, unsigned long nsec);
      virtual bool HandleData (const std::string& name, float* x, 
                        int dim1, int dim2=0, int dim3=0, int dim4=0) {
         return false; }
      virtual bool HandleData (const std::string& name, 
                        std::complex<float>* x,
                        int dim1, int dim2=0, int dim3=0, int dim4=0) {
         return false; }
      virtual bool HandleData (const std::string& name, double* x, 
                        int dim1, int dim2=0, int dim3=0, int dim4=0);
      virtual bool HandleData (const std::string& name, 
                        std::complex<double>* x,
                        int dim1, int dim2=0, int dim3=0, int dim4=0) {
         return false; }
   };

   class xsilHandlerQueryHistogram : public xsilHandlerQuery {
   protected:
      std::vector<Histogram1>*    fDat1;
      std::vector<Histogram2>*    fDat2;
   
   public:
      explicit xsilHandlerQueryHistogram (std::vector<Histogram1>& h1)
      : fDat1 (&h1), fDat2 (0) {
      }
      explicit xsilHandlerQueryHistogram (std::vector<Histogram2>& h2)
      : fDat1 (0), fDat2 (&h2) {
      }
      explicit xsilHandlerQueryHistogram (std::vector<Histogram1>& h1,
                        std::vector<Histogram2>& h2)
      : fDat1 (&h1), fDat2 (&h2) {
      }
      virtual xsilHandler* GetHandler (const attrlist& attr);
   };

} // namespace xml

#endif // _GDS_XSILHISTOGRAM_HH

---