---

acdReconAlg.h

Go to the documentation of this file.

// Interface Dependencies -------------------------------------------------

#ifndef __ACD_RECON_H
#define __ACD_RECON_H 1
#include "GaudiKernel/Algorithm.h"

#include "data/IVetoData.h"
#include "ntupleWriterSvc/INTupleWriterSvc.h"
#include "GlastSvc/GlastDetSvc/IGlastDetSvc.h"

#include "instrument/DetectorConverter.h"
#include "instrument/DetectorVolumes.h"
#include "instrument/Scintillator.h"
#include "geometry/Vector.h"
#include <map>

class Ray;
class IVetoData;

class acdReconAlg : public Algorithm
{
  public:
      acdReconAlg(const std::string& name, ISvcLocator* pSvcLocator); 

      StatusCode initialize();
      StatusCode execute();
      StatusCode finalize();
 
  protected:
 
  private:

      // reset all member variables for each iteration
      void clear ();

      void getParameters ();

      StatusCode reconstruct (const IVetoData* v);

      StatusCode acdTileDOCA();

      StatusCode writeNTuple();

      double DOCA (const Point &x0, const Vector &dir, std::vector<double> &doca_values);

      double hitTileDist(const Point &x0, const Vector &dir);

      static double threshold_energy;
      static int xNumTowers, yNumTowers, xNumTopTiles, yNumTopTiles, numSideRows;

    // Data Members for Associations

      const IVetoData *m_AcdData;

      // record of the tile with the minimum Distance of Closest Approach
      IVetoData::Tile m_hit_tile;

      std::string m_tupleName;
      INTupleWriterSvc *ntupleWriteSvc;

      IGlastDetSvc *m_glastDetSvc;

      // Items that will be output to the ntuple
      double m_totEnergy, m_tileCount, m_gammaDOCA, m_DOCA, m_act_dist;
      std::vector<double> m_rowDOCA_vec;


      // nested class that is also a visitor
      class TileParams;
      TileParams* m_tileParams; // the ACD parameters

};

//   Tile Rep Class - used to retrieve parameters of the ACD tiles
class   acdReconAlg::TileParams : public DetectorConverter  {
public:
    TileParams(){ m_tileDims.clear(); }
    void   forward ( const Scintillator& d )
    {
        if ( d.empty()) return;
        unsigned int id = d.type();
        const GlastDetector::_Volume * v = d.volume();
        const DetectorBox* db = dynamic_cast<const DetectorBox*> (d.volume());
        const Box& b = *db->box();
        double len = b.length();
        double wid = b.width();
        double height = b.height();
        m_tileDims[id] = Vector(len, wid, height);
    }
    void clear(){
        m_tileDims.clear();
    }
    void update(){};

    double length(int id) { return m_tileDims[id].x(); };
    double width(int id) {return m_tileDims[id].y(); };
    double height(int id) {return m_tileDims[id].z(); };

private:
    std::map<int, Vector> m_tileDims;
};

#endif

---