PndGemHit.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                    PndGemHit header file                      -----
// -----                  Created by R.Kliemt          -----
// -------------------------------------------------------------------------


#ifndef PNDGEMHIT_H
#define PNDGEMHIT_H


#include "TVector3.h"
#include "TString.h"
#include "FairHit.h"

#include <stdio.h>
#include <iostream>

class PndGemHit : public FairHit
{
  friend std::ostream& operator<< (std::ostream& out, const PndGemHit& hit){
    out << "Gem hit" /*in detector " << hit.GetDetName()*/ << " at ("
            << hit.GetX() << ", " << hit.GetY() << ", " << hit.GetZ() << ") cm "
      << " with " << hit.GetCharge() << " e" << ", Cluster No. " << hit.GetRefIndex() << std::endl;

    return out;
  }

 public:

  PndGemHit();


  PndGemHit(Int_t detID, TVector3& pos, TVector3& dpos, Int_t index, Double_t charge, Int_t NDigiHits, TString fromStr = "GEMPoint");
  PndGemHit(Int_t detID, TVector3& pos, TVector3& dpos, Int_t digi1, Int_t digi2, Double_t dr, Double_t dp, Int_t index, TString fromStr = "GEMCluster");
  PndGemHit(Int_t detID, TVector3& pos, TVector3& dpos, Double_t charge, Double_t time, Int_t digi1, Int_t digi2, Double_t dr, Double_t dp, Int_t index, TString fromStr = "GEMCluster");

  //PndGemHit(PndGemHit& c);
  virtual ~PndGemHit();


  void SetCharge(Double_t charge){ fCharge    = charge;}
  void SetNDigiHits(Int_t pixel) { fNDigiHits = pixel;}
  void SetBotIndex(Int_t id)       { fBotIndex  = id;}
  void SetErrors(Double_t dr, Double_t dp) { fDr = dr; fDp = dp; };
  void SetDigiNr(Int_t digi1, Int_t digi2) { fDigiNr[0] = digi1; fDigiNr[1] = digi2; };

  Double_t GetCharge()    const { return fCharge;}
  Int_t    GetNDigiHits() const { return fNDigiHits;}
  TVector3 GetPosition()  const { return TVector3(fX,fY,fZ);}
  Int_t    GetBotIndex()  const { return fBotIndex;}
  Double_t GetEloss()     const { return (fCharge * 3.61e-9);}  // 3.6 eV/Electron in Silicon

  Double_t GetDr()        const { return fDr;  };
  Double_t GetDp()        const { return fDp;  };
  Int_t    GetDigiNr(Int_t iside) const { if ( iside*(iside-1) == 0 ) return fDigiNr[iside]; return -1;};
  
  Int_t    GetSystemId()   const { 
    return ( ( fDetectorID & (  31<<27) ) >> 27); }
  Int_t    GetStationNr()  const { 
    return ( ( fDetectorID & (8191<< 8) ) >>  8 ); }
  Int_t    GetSensorNr()   const {  // sensor number within station
    return ( ( fDetectorID & (   3<< 6) ) >>  6 ); }

  //   // CAUTION The errors in the GemHit are LOCAL, but the coordinates are in the LAB
  //
  //   Double_t GetDxLocal() const { return fDx;};
  //   Double_t GetDyLocal() const { return fDy;};
  //   Double_t GetDzLocal() const { return fDz;};
  //   void PositionErrorLocal(TVector3& dpos) const;
  
  //   /** overloaded accessors **/
  //   Double_t GetDx() {return GetD(0);};
  //   Double_t GetDy() {return GetD(1);};
  //   Double_t GetDz() {return GetD(2);};
  //   void PositionError(TVector3& dpos);
 
  virtual void Print(const Option_t* opt = 0) const;

 private:

  Int_t    fDigiNr[2];
  Double_t fDr;  // error in radius calculation
  Double_t fDp;  // error in angle calculation

  Double_t fCharge; 
  Int_t fNDigiHits; 
                    //if more then in the cluster cand, look for the bottom cluster.
  Int_t fBotIndex; 
  ClassDef(PndGemHit,1);

};

// inline void PndGemHit::PositionErrorLocal(TVector3& dpos) const {
//   dpos.SetXYZ(fDx, fDy, fDz);
// }
//
// inline void PndGemHit::PositionError(TVector3& dpos) {
//   dpos.SetXYZ(GetDx(), GetDy(), GetDz());
// }

#endif