/// 
/// For questions or suggestions about this module please contact the 
/// current responsible and CC in the oaAnalysis package manager.
///
/// 21-Jul-2010: Current responsible for this module is, 
/// Ben Still (b.still [*a*t*] qmul.ac.uk) 
///
#ifndef TP0DECALReconModule_hxx_seen
#define TP0DECALReconModule_hxx_seen

#include <TNamed.h>
#include <TTree.h>
#include <string>
#include <ICOMETEvent.hxx>
#include <IReconCluster.hxx>
#include <IIntegerDatum.hxx>
#include "IReconShower.hxx"
#include "IReconTrack.hxx"
#include "IReconPID.hxx"


#include "TClonesArray.h"

#include "IAnalysisReconModuleBase.hxx"

namespace COMET
{
  class IP0DReconECalModule;
};

/// IP0DReconECalModule
///
/// Summary information about ECAL clusters.
/// The primary level store the RunNo, EventID, number of tracks and showers and TCloneArrays 
/// of TP0DECALTracks and TP0DECALShowers which are TObjects that are root browsable.
///
///
class COMET::IP0DReconECalModule : public IAnalysisReconModuleBase
{

  public:
    /// IP0DReconECalTrack
    /// The main object that contains ECAL track information
    /// Feel free to add to.
    class IP0DReconECalTrack : public TObject
    {
      public:
        virtual ~IP0DReconECalTrack();

        /// The object UniqueID which can be used to match this object to the Global Recon object
        UInt_t UniqueID;

        //     /// The name of the algorithm that created this reconstruction object.
        //     std::string AlgorithmName;

        double Curvature;        ///< The track curvature around the X axis.


        TVector3 Direction;        ///< The track direction.


        double EDeposit;        ///< the energy deposited in the track.


        ///< the full name of this datum. The full name places the datum in the
        ///< data tree and has the syntax //xxxx/yyyy/zzzz.
        //     TString FullName;

        int NDOF;        ///< the number of degrees of freedom in the reconstruction.


        TLorentzVector Position;        ///< the track starting position.


        double Quality;        ///< the goodness of fit for the reconstruction.


        TVector3 Width;        ///< the width of the track.


      double AMR;                 ///< The Pid Variables

      double CWTrackWidth;         ///< The Pid Variables

      double Max_Ratio;          ///< The Pid Variables

      double NormChargeSD;         ///< The Pid Variables

      double PathChargeRatio;      ///< The Pid Variables

      double TrShval;        ///< The Pid Variables

      int NHits;     

      ///<

      TVector3 BenCWPosition;      ///<Ben CW position.


      int StackNo;      ///< Stack of interaction

      
        ClassDef(IP0DReconECalModule::IP0DReconECalTrack, 1);
    };

  public:
    /// IP0DReconECalShower
    /// The main object that contains ECAL cluster information
    /// Feel free to add to.
    class IP0DReconECalShower : public TObject
    {
      public:
        virtual ~IP0DReconECalShower();

        UInt_t UniqueID;        ///< The object UniqueID which can be used to match this object to the Global Recon object


        //     std::string AlgorithmName;        ///< The name of the algorithm that created this reconstruction object.


        TVector3 ConeAngle;        ///< The shower opening angle.


        TVector3 Direction;        ///< The shower direction.


        double EDeposit;        ///< the energy deposited in the shower.


        ///< the full name of this datum. The full name places the datum in the
        ///< data tree and has the syntax //xxxx/yyyy/zzzz.
        //     TString FullName;

        int NDOF;        ///< the number of degrees of freedom in the reconstruction.


        TLorentzVector Position;        ///< the shower starting position.


        double Quality;        ///< the goodness of fit for the reconstruction.


      ///< The Pid Variables - AMR/F ChargeStandardDeviation/F MaxWidth/F MeanCharge/F NumHits/I Track_Length/I
      double AMR;      ///< The Pid Variables         
      double CWTrackWidth;     ///< The Pid Variables  
      double Max_Ratio;       ///< The Pid Variables 
      double NormChargeSD;    ///< The Pid Variables   
      double PathChargeRatio;   ///< The Pid Variables 
      double TrShval;     ///< The Pid Variables  
      int NHits;
      /// 

      TVector3 BenCWPosition;      ///<Ben CW position.


      int StackNo;      ///< Stack of interaction


        ClassDef(IP0DReconECalModule::IP0DReconECalShower, 1);
    };
      
    class IP0DECALLayerCharge : public TObject
    {
      public:
  virtual ~IP0DECALLayerCharge();

  double FirstTime;	///< first time of hits in hit


  double LastTime;	///< last time of hits in hit


  double Charge[6];	///< summed charge in each of 6 layers; index 0 is the innermost layer


        TString module_name;	///< string dictating the module of the detector

      
        ClassDef(IP0DReconECalModule::IP0DECALLayerCharge, 1);
    };


  public:

    IP0DReconECalModule(const char *name = "P0DECal",
      const char *title = "P0D ECal Recon Module");
    virtual ~IP0DReconECalModule();

    virtual Bool_t ProcessFirstEvent(COMET::ICOMETEvent& event);

  protected:

    virtual void InitializeModule();

    virtual void InitializeBranches();
    virtual bool FillTree(COMET::ICOMETEvent&);

  public:

    //Tree Entries
    Int_t        fNReconTracks;    ///< Number of objects in clones array.
    Int_t        fNReconShowers;     ///< Number of objects in clones array.
    Int_t        fNLayerCharge;      ///< Number of objects in clones array.

    TClonesArray *fReconTracks;        ///< TClonesArray of IP0DReconECalTrack Objects
    TClonesArray *fReconShowers;      ///< TClonesArray of IP0DReconECalShower Objects
    TClonesArray *fLayerCharge;      ///< TClonesArray of IP0DECALLayerCharge Objects

  private:

    ClassDef(IP0DReconECalModule,2);

};
#endif