---

AcdReconDigiAlg.h

Go to the documentation of this file.

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

#ifndef __ACD_DIGIRECON_H
#define __ACD_DIGIRECON_H 1
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/ObjectVector.h"

#include "GlastEvent/Digi/AcdDigi.h"
#include "ntupleWriterSvc/INTupleWriterSvc.h"
#include "GlastSvc/GlastDetSvc/IGlastDetSvc.h"
#include "geometry/Point.h"

#include <string>
#include <vector>
#include <map>

class Ray;
class Point;
class Vector;

class AcdReconDigiAlg : public Algorithm
{
  public:
      AcdReconDigiAlg(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 ();

      void readThresholds();

      StatusCode reconstruct (const AcdDigiVector* 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);

      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;

      AcdDigiVector *m_AcdData;

      std::map<int, Point> m_tileCenters;
      std::map<int, Vector> m_tileDims;
      std::map<int, double> m_vetoThresh;

      int m_numSideRows;

      std::string m_xmlFileStr, m_vetoThreshFileStr;
};

#endif

---