---

SiClusters.h

Go to the documentation of this file.

// $Id: SiClusters.h,v 1.3 2000/12/08 23:01:08 igable Exp $

#ifndef SI_CLUSTER_H
#define SI_CLUSTER_H 1

#include "geometry/Point.h"
#include "data/SiData.h"

#include <iostream>

class SiClusters
{
 public:
  class Cluster;

 private:

  class Tray
    {
    private:
      friend class SiClusters;

      public:

          Tray (int m);
          ~Tray ();

          void flagHit (unsigned int idx);
          void flagHit (enum SiData::Axis , unsigned int );
          void unflagHit (unsigned int idx);
          void unflagHit (enum SiData::Axis , unsigned int );
          void clearHitFlags ();

        // Data Members for Class Attributes

          unsigned int module;
          std::vector<Cluster *> xhitList;
          std::vector<Cluster *> yhitList;
    };



  public:

  class Cluster
    {
    public:
      friend class SiClusters;
      Cluster (Point p, int n);
      Cluster (Point p, int n, int ns);

      int flag;
      int number;
      int num_noise;  // Number of noise strips in cluster
      Point pos;
    };

  public:

      explicit SiClusters (int  ntrays );
      // see note on size_type in .cxx file
      ~SiClusters ();


      //        setup call to transfer data from DetectorVisitor SiData
      //        and compute clusters
      void loadFrom (const SiData& data);

      //        -------------------- access functions for data
      //        analysis------------------------
      void clear ();

      //        number of  strip clusters in the given tray (tray=0..11)
      int nHits (enum SiData::Axis a, int trayNum) const;

      Point nextHit (int trayN, int* index);

      //        Function for cycling over allowed hit combinations
      Point nextHit (enum SiData::Axis a, int trayN, int* index);

      Point nearestHitInside( unsigned int trayNum, float dR,
                              const Point& x0, int* index);

      Point nearestHitInside( unsigned int trayNum, float dX, float dY,
                              const Point& x0, int* index);

      Point nearestHitInside( enum SiData::Axis a, unsigned int trayNum,
                              float dX, const Point& x0, int* index);

      Point nearestHitOutside( unsigned int trayNum, float dR,
                               const Point& x0, int* index);

      //        Functions for adding hits to trajectories
      Point nearestHitOutside( enum SiData::Axis a, unsigned int trayNum,
                               float dX, const Point& x0, int* index);

      Point meanHit (unsigned int trayNum, float* rmsX, float* rmsY);

      Point meanHit (enum SiData::Axis axis, unsigned int trayNum, float* rms);

      //        Functions for giving global pattern information
      Point meanHitInside (unsigned int trayNum, float dXY, const Point& x0);

      //        Functions for analyzing hit patterns
      int numberOfHitsNear( unsigned int trayNum, float dX, float dY,
                    const Point& x0);
      int numberOfHitsNear( unsigned int trayNum, float dR, const Point& x0);
          int numberOfHitsNear( enum SiData::Axis a, unsigned int trayNum, float dR,
                            const Point& x0);

      int clusterSize(enum SiData::Axis a, unsigned int index);

      Point clusterPosition(enum SiData::Axis a, unsigned int index);


      int clusterNoise(enum SiData::Axis a, unsigned int index);

      bool hitFlagged(enum SiData::Axis a, unsigned int index);

      void flagHit (unsigned int index);
      void flagHit (enum SiData::Axis a, unsigned int index);

      void unflagHit (unsigned int index);

      void unflagHit (enum SiData::Axis a, unsigned int index);

       //       Funcitons for managing hit bank for pattern recognition
      void clearHitFlags ();

      void printOn (std::ostream& cout = std::cout);

          static int tower(int indexhit);

   private:

      SiClusters::Tray* getTrayFor (int lyr, int mod);

       void transferData (enum SiData::Axis a, const SiData& data,
                          int lyr, int start, int stop);

    // Data Members for Associations

       std::vector< std::vector< Tray *> > trayList;

};


// Class SiClusters::Tray

inline SiClusters::Tray::Tray (int m) : module(m) {}

inline SiClusters::Tray::~Tray ()
{
    std::vector<Cluster *>::iterator it = xhitList.begin();
  for ( ; it != xhitList.end(); it++ ) {
    delete *it;
  }
  xhitList.clear();

  for ( it = yhitList.begin(); it != yhitList.end(); it++ ) {
    delete *it;
  }
  yhitList.clear();
}

// Class SiClusters::Cluster

inline SiClusters::Cluster::Cluster (Point p, int n, int ns)
    : flag(0), number(n), num_noise(ns), pos(p) {}

inline SiClusters::Cluster::Cluster (Point p, int n)
    : flag(0), number(n), num_noise(0), pos(p) {}

#endif

---