THcDCHit.h

#ifndef ROOT_THcDCHit
#define ROOT_THcDCHit

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
// THcDCHit                                                                 //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "TObject.h"
#include "THcDCWire.h"
#include "THcDriftChamberPlane.h"
#include "THcDriftChamber.h"
#include <cstdio>

class THcDCHit : public TObject {

public:
  THcDCHit( THcDCWire* wire=NULL, Int_t rawnorefcorrtime=0, Int_t rawtime=0, Double_t time=0.0,
    THcDriftChamberPlane* wp=0) :
    fWire(wire), fRawNoRefCorrTime(rawnorefcorrtime), fRawTime(rawtime), fTime(time), fWirePlane(wp),
      fDist(0.0), fLR(0), ftrDist(kBig) {
      if (wire) ConvertTimeToDist();
      fCorrected = 0;
    }
  virtual ~THcDCHit() {}

  virtual Double_t ConvertTimeToDist();
  Int_t  Compare ( const TObject* obj ) const;
  Bool_t IsSortable () const { return kTRUE; }
  virtual void Print( Option_t* opt="" ) const;

  // Get and Set Functions
  THcDCWire* GetWire() const { return fWire; }
  Double_t GetWireSigma() const { return fWire->GetSigma(); }
  Int_t    GetWireNum() const { return fWire->GetNum(); }
  Int_t    GetRawTime() const { return fRawTime; }
  Int_t    GetRawNoRefCorrTime() const { return fRawNoRefCorrTime; }
  Double_t GetTime()    const { return fTime; }
  Double_t GetDist()    const { return fDist; }
  Double_t GetPos()     const { return fWire->GetPos(); } //Position of hit wire
  Double_t GetCoord()   const { return fCoord; }
  Double_t GetdDist()   const { return fdDist; }
  Int_t    GetCorrectedStatus() const { return fCorrected; }

  THcDriftChamberPlane* GetWirePlane() const { return fWirePlane; }

  void     SetWire(THcDCWire * wire) { fWire = wire; ConvertTimeToDist(); }
  void     SetRawTime(Int_t time)     { fRawTime = time; }
  void     SetTime(Double_t time)     { fTime = time; }
  void     SetDist(Double_t dist)     { fDist = dist; }
  void     SetLeftRight(Int_t lr)   { fCoord = GetPos() + lr*fDist; fLR=lr;}
  Int_t    GetLR() { return fLR; }
  void     SetdDist(Double_t ddist)   { fdDist = ddist; }
  void     SetFitDist(Double_t dist)  { ftrDist = dist; }
  Int_t    GetPlaneNum() const { return fWirePlane->GetPlaneNum(); }
  Int_t    GetPlaneIndex() const { return fWirePlane->GetPlaneIndex(); }
  Int_t    GetChamberNum() const { return fWirePlane->GetChamberNum(); }
  void     SetCorrectedStatus(Int_t c) { fCorrected = c; }

  Int_t GetReadoutSide() {
    Int_t pln = fWirePlane->GetPlaneNum();
    Int_t tb = fWirePlane->GetReadoutTB();
    Int_t wn = fWire->GetNum();
    Int_t version = fWirePlane->GetVersion();

    //if new HMS
    if (version == 1) {
      if ((pln>=3 && pln<=4) || (pln>=9 && pln<=10)) {
	if (tb>0) {
	  if (wn < 60) {
	    fReadoutSide = 2;
	  } else {
	    fReadoutSide = 4;
	  }
	} else {
	  if (wn < 44) {
	    fReadoutSide = 4;
	  } else {
	    fReadoutSide = 2;
	  }
	}
      } else {
	if (tb>0) {
	  if (wn < 51) {
	    fReadoutSide = 2;
	  } else if (wn >= 51 && wn <= 64) {
	    fReadoutSide = 1;
	  } else {
	    fReadoutSide =4;
	  }
	} else {
	  if (wn < 33) {
	    fReadoutSide = 4;
	  } else if (wn >=33 && wn<=46) {
	    fReadoutSide = 1;
	  } else {
	    fReadoutSide = 2;
	  }
	}
      }
    } else{//appplies SHMS DC configuration
      //check if x board
      if ((pln>=3 && pln<=4) || (pln>=9 && pln<=10)) {
	if (tb>0) {
	  if (wn < 49) {
	    fReadoutSide = 4;
	  } else {
	    fReadoutSide = 2;
	  }
	} else {
	  if (wn < 33) {
	    fReadoutSide = 2;
	  } else {
	    fReadoutSide = 4;
	  }
	}
      } else { //else is u board
	if (tb>0) {
	  if (wn < 41) {
	    fReadoutSide = 4;
	  } else if (wn >= 41 && wn <= 63) {
	    fReadoutSide = 3;
	  } else if (wn >=64 && wn <=69) {
	    fReadoutSide = 1;
	  } else {
	    fReadoutSide = 2;
	  }
	} else {
	  if (wn < 39) {
	    fReadoutSide = 2;
	  } else if (wn >=39 && wn<=44) {
	    fReadoutSide = 1;
	  } else if (wn>=45 && wn<=67) {
	    fReadoutSide = 3;
	  } else {
	    fReadoutSide = 4;
	  }
	}
      }
    }
    return fReadoutSide;
  }


protected:
  static const Double_t kBig;  //!

  THcDCWire*  fWire;     // Wire on which the hit occurred
  Int_t       fRawNoRefCorrTime;  // raw TDC value (channels) 
  Int_t       fRawTime;  // TDC value (channels) reference time subtracted
  Double_t    fTime;     // Time corrected for time offset of wire (s)
  THcDriftChamberPlane* fWirePlane; //! Pointer to parent wire plane
  Double_t    fDist;     // (Perpendicular) Drift Distance
  Int_t       fLR;       // +1/-1 which side of wire
  Double_t    fCoord;    // Actual coordinate of hit
  Double_t    fdDist;    // uncertainty in fDist (for chi2 calc)
  Double_t    ftrDist;   // (Perpendicular) distance from the track
  Int_t       fCorrected; // Has this hit been corrected?
  Int_t       fReadoutSide; // Side where wire is read out. 1-4 is T/R/B/L from beam view for new chambers.

  THcDriftChamber* fChamber; //! Pointer to parent wire plane


private:
  THcDCHit( const THcDCHit& );
  THcDCHit& operator=( const THcDCHit& );

  ClassDef(THcDCHit,2)             // Drift Chamber Hit
};

////////////////////////////////////////////////////////////////////////////////

#endif