---

CalRecon.h

Go to the documentation of this file.

//   Reconstruct Events in CsI Data
//    by Jay Norris and Heather Arrighi Muise
// $Id: CalRecon.h,v 1.5 2001/01/22 15:53:55 burnett Exp $

#ifndef __CAL_RECON_H
#define __CAL_RECON_H 1

#include "reconstruction/Recon.h"

#include "geometry/Vector.h"
#include "geometry/Point.h"

#include <vector>

class CsIData;
class ReconVisitor;
class CalModules;
namespace gui { class DisplayRep; }

class CalRecon : public Recon
{

  public:
      CalRecon ();

      ~CalRecon ();


      void getParameters ();

      void accept (ReconVisitor& rv);

      void reconstruct (const CsIData* csiData);

      void clear ();

      void draw (gui::DisplayRep& v);

      void moments_anal (int counter);

      double RLinePnt (Point xPos);

      Point center () { return(m_cen); }

      Vector direction () { return(m_dir); }

          std::vector<double> Get_Elayer() { return E_layer; }


  private:
      unsigned int num_layers;
      int s_xNum;

      double mod_width, s_xtal_gap, s_rad_len, s_wall_thickness, s_wall_gap;

          double light_att;
      const char* calmat;

      float chisq;

      double  moment[3], moment1, moment2, moment3;

      double m_Rsq_mean;
      double eTotal;

      int numXtals, ntothits;

      LbldData::iterator i_CsI_Energy, i_CsI_Etop, i_CsI_Ebottom, i_CsI_Esides, i_CsI_S_or_B;      LbldData::iterator i_CsI_Xtals, i_CsI_mean_X, i_CsI_mean_Y, i_CsI_mean_Z;
      LbldData::iterator i_CsI_dir_X, i_CsI_dir_Y, i_CsI_dir_Z, i_CsI_chisq, i_CsI_mips, i_CsI_Xtals_trunc;
      LbldData::iterator i_CsI_Xtals_keep, i_CsI_Rsq_mean, i_CsI_D3Mean;
      LbldData::iterator i_CsI_E_moment, i_CsI_E_disp, i_CsI_moment1, i_CsI_moment2, i_CsI_moment3;
      LbldData::iterator i_CsI_MaxMod, i_CsI_ModXtals;
      std::vector<LbldData::iterator>i_CsI_eLayer;
      LbldData::iterator i_Diode_Energy;

    // Data Members for Associations

      CalModules *calModules;

      Point bCenter;

      Vector principal_axis[3];

      Point m_cen;

      Vector m_dir;

      //        new list to keep track of positions using
      //        left and right response of the ends of the logs
      std::vector<int>Use_R;
      std::vector<int>numlayhits;
      std::vector<double>E_hits;
      std::vector<double>E_layer;

      std::vector< std::vector< class Point >* > LRpos;
      std::vector <char*> layerLabel;

};


#endif

---