---

LSQFit.h

Go to the documentation of this file.

//## begin module.cm preserve=no
//        %X% %Q% %Z% %W%
//## end module.cm

//## begin module.cp preserve=no
//## end module.cp

//## Module: LSQFit; Package specification
//## Subsystem: reconstruction
//      d:\code\glastsim\reconstruction
//## Source file: d:\code\glastsim\reconstruction\LSQFit.h

#ifndef __LSQFIT_H
#define __LSQFIT_H 1

#include <vector>

#include "geometry/Point.h"


//## begin module.declarations preserve=no
//## end module.declarations

//## begin module.additionalDeclarations preserve=yes
//## end module.additionalDeclarations


//## Class: FitData
//## Category: reconstruction
//## Subsystem: reconstruction
//## Persistence: Transient
//## Cardinality/Multiplicity: n


class FitData
{
  //## begin FitData.initialDeclarations preserve=yes
private:
        friend class LSQFit;
  //## end FitData.initialDeclarations

  public:
    //## Constructors (specified)
      //## Operation: FitData%907927192; C++
      FitData (unsigned int i, float dX, float dZ, float dwX)
        //## begin FitData::FitData%907927192.hasinit preserve=no
        //## end FitData::FitData%907927192.hasinit
        //## begin FitData::FitData%907927192.initialization preserve=yes
        : index( i ), x( dX), z(dZ), wt( dwX) 
        //## end FitData::FitData%907927192.initialization
      {
        //## begin FitData::FitData%907927192.body preserve=yes
        //## end FitData::FitData%907927192.body
      }

      //## Operation: FitData%1210066552; C++
      FitData ()
        //## begin FitData::FitData%1210066552.hasinit preserve=no
        //## end FitData::FitData%1210066552.hasinit
        //## begin FitData::FitData%1210066552.initialization preserve=yes
        //## end FitData::FitData%1210066552.initialization
      {
        //## begin FitData::FitData%1210066552.body preserve=yes
        //## end FitData::FitData%1210066552.body
      }

      //## Operation: FitData%666052249; C++
      FitData (const FitData& data)
        //## begin FitData::FitData%666052249.hasinit preserve=no
        //## end FitData::FitData%666052249.hasinit
        //## begin FitData::FitData%666052249.initialization preserve=yes
        //## end FitData::FitData%666052249.initialization
      {
        //## begin FitData::FitData%666052249.body preserve=yes
          index = data.index;
          x = data.x;
          z = data.z;
          wt = data.wt;
          res = data.res;
        //## end FitData::FitData%666052249.body
      }

    //## Destructor (specified)
      //## Operation: ~FitData%790131244; C++
      ~FitData ()
      {
        //## begin FitData::~FitData%790131244.body preserve=yes
        //## end FitData::~FitData%790131244.body
      }

    // Data Members for Class Attributes

      //## Attribute: index
      //        unique index tagging this data point
      //## begin FitData::index.attr preserve=no  public: unsigned int {VA} 
      unsigned int index;
      //## end FitData::index.attr

      //## Attribute: x
      //        2D space point for fit
      //## begin FitData::x.attr preserve=no  public: float {VA} 
      float x;
      //## end FitData::x.attr

      //## Attribute: z
      //        2D space point for fit
      //## begin FitData::z.attr preserve=no  public: float {VA} 
      float z;
      //## end FitData::z.attr

      //## Attribute: wt
      //        fit weights
      //## begin FitData::wt.attr preserve=no  public: float {VA} 
      float wt;
      //## end FitData::wt.attr

      //## Attribute: res
      //        fit Residuals
      //## begin FitData::res.attr preserve=no  public: float {VA} 
      float res;
      //## end FitData::res.attr

    // Additional Public Declarations
      //## begin FitData.public preserve=yes
      //## end FitData.public

  protected:
    // Additional Protected Declarations
      //## begin FitData.protected preserve=yes
      //## end FitData.protected

  private:
    // Additional Private Declarations
      //## begin FitData.private preserve=yes
      //## end FitData.private

  private:  //## implementation
    // Additional Implementation Declarations
      //## begin FitData.implementation preserve=yes
      //## end FitData.implementation

};

//## begin FitData.postscript preserve=yes
    bool operator<(const FitData&, const FitData&);
    bool operator==(const FitData&, const FitData&);


#include <stdlib.h>
#include "geometry/Vector.h"

#include <iostream>
#include <vector>

// End Implementation Dependencies -----------------------------------

class Ray;

//## end FitData.postscript



//## Class: LSQFit
//## Category: reconstruction
//## Subsystem: reconstruction
//## Persistence: Transient
//## Cardinality/Multiplicity: n

//## Uses: <unnamed>; Point { -> }
//## Uses: <unnamed>; ostream { -> }

class LSQFit 
{
  //## begin LSQFit.initialDeclarations preserve=yes
  //## end LSQFit.initialDeclarations

  public:
    //## Constructors (specified)
      //## Operation: LSQFit%-287297432; C++
      //        End Implementation Dependencies
      //        -----------------------------------
      LSQFit ();

      //## Operation: LSQFit%1816016125; C++
      LSQFit (const LSQFit& fit);

    //## Destructor (specified)
      //## Operation: ~LSQFit%1575098630; C++
      virtual ~LSQFit ();


    //## Other Operations (specified)
      //## Operation: clear%1564145804; C++
      //        ---------------------------------------------------------
      //        ---
      //        Member functions to access fit results
      //        ---------------------------------------------------------
      //        ---
      void clear ();

      //## Operation: position%-438469547; C++
      float position (float s)
      {
        //## begin LSQFit::position%-438469547.body preserve=yes
return (x0 + slopeX*s);
        //## end LSQFit::position%-438469547.body
      }

      //## Operation: slope%-1768319281; C++
      float slope ()
      {
        //## begin LSQFit::slope%-1768319281.body preserve=yes
return slopeX;
        //## end LSQFit::slope%-1768319281.body
      }

      //## Operation: chiSquare%-1141137480; C++
      float chiSquare () const
      {
        //## begin LSQFit::chiSquare%-1141137480.body preserve=yes
return chisq;
        //## end LSQFit::chiSquare%-1141137480.body
      }

      //## Operation: scatter%206744814; C++
      float scatter () const
      {
        //## begin LSQFit::scatter%206744814.body preserve=yes
return rmsResid;
        //## end LSQFit::scatter%206744814.body
      }

      //## Operation: maxResidual%-1294063902; C++
      float maxResidual (int* index);

      unsigned int numDataPoints () const;

      //## Operation: compareFits%-248894077; C++
      int compareFits (const LSQFit& fit);

      //## Operation: getHit%588524840; C++
      Point getHit (unsigned int i);

      //## Operation: getHitIndex%-1565627155; C++
      //        ---------------------------------------------------------
      //        ---
      //        Member functions to store data and do the fit
      //        ---------------------------------------------------------
      //        ---
      unsigned int getHitIndex (unsigned int i);

      //## Operation: doFit%-966633323; C++
      float doFit ();

      //## Operation: addData%860586691; C++
      int addData (float x, float z, float wx);

      //## Operation: addData%-2047337558; C++
      int addData (unsigned int idx, float x, float z, float wx);

      //## Operation: removeData%-2110859248; C++
      int removeData (unsigned int index);

      //## Operation: truncateFit%-10029783; C++
      //        ---------------------------------------------------------
      //        ---
      //        Displaying and Printing
      //        ---------------------------------------------------------
      //        ---
      float truncateFit ();

      //## Operation: printOn%353614648; C++
      //        ---------------------------------------------------------
      //        ----
      virtual void printOn (std::ostream& os = std::cout) const;

  protected:
    // Data Members for Class Attributes

      float x0;
      float slopeX;
      float chisq;
      float rmsResid;

    // Data Members for Associations

      std::vector< FitData > dataList;

};

#endif

---