PndSdsHit.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                    PndSdsHit header file                      -----
// -----                  Created by T.Stockmanns            -----
// -------------------------------------------------------------------------


#ifndef PNDSDSHIT_H
#define PNDSDSHIT_H


#include "TVector3.h"
#include "TString.h"
#include "FairHit.h"
#include "PndDetectorList.h"
#include "TMatrixD.h"
//#include "PndSingleLinkedData.h" 

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

class PndSdsHit : public FairHit
{
  friend std::ostream& operator<< (std::ostream& out, const PndSdsHit& hit){
    out << "PndSdsHit in " << hit.GetSensorID() << " at" << std::endl;
    out << "(" << hit.GetX() << ", " << hit.GetY() << ", " << hit.GetZ() << ") cm "
        << " with " << hit.GetCharge() << " e" << " at Time: " << hit.GetTimeStamp() << " +/- " << hit.GetTimeStampError() << " , Cluster No. " << hit.GetClusterIndex();
    if (hit.GetBotIndex() > -1)
        out << " and bottom " << hit.GetBotIndex();
    if (hit.GetRefIndex() > -1)
      out << ", mc point id = " << hit.GetRefIndex();
    else 
      out << ", noise hit without mc poit";
    out << std::endl;

    out <<"hit.GetClusterIndex() "<<hit.GetClusterIndex() <<std::endl;
    out <<"hit.GetBotIndex() "<< hit.GetBotIndex() <<std::endl;
    out <<"hit.GetRefIndex() "<< hit.GetRefIndex()<<std::endl;
    out <<"hit.GetNDigiHits() "<< hit.GetNDigiHits()<<std::endl;
    out <<"hit.GetCharge() "<< hit.GetCharge()<<"("<<hit.GetEloss()<<" GeV)"<<std::endl;
    out <<"hit.GetSensorID() "<< hit.GetSensorID()<<std::endl;
    out <<"Error values in FairHit part: (dx,dy,dz) = ("<< hit.GetDx()<<","<<hit.GetDy()<<","<<hit.GetDz()<<")"<<std::endl;
    return out;
  }

 public:

  PndSdsHit();


  PndSdsHit(Int_t detID, Int_t sensorID,
            TVector3& pos, TVector3& dpos, Int_t clindex, Double_t charge, Int_t NDigiHits, Int_t mcindex);

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


  void SetSensorID(Int_t sensorID)  { fSensorID   = sensorID;}
  void SetCharge(Double_t charge){ fCharge    = charge;}
  void SetNDigiHits(Int_t pixel) { fNDigiHits = pixel;}
  void SetClusterIndex(Int_t datasource, Int_t id, Int_t fileId = -1, Int_t evtId = -1) {
          fClusterIndex = id;
          SetLink(FairLink(fileId, evtId, datasource, id));
  }
  void SetBotIndex(Int_t id)     { fBotIndex  = id;}
  void SetCov(TMatrixD cov);
  
  //TString     GetDetName()      const { return fDetName;}
  Int_t     GetSensorID()                       const { return fSensorID;}
  Double_t      GetCharge()       const { return fCharge;}
  Int_t           GetNDigiHits()    const { return fNDigiHits;}
  TVector3  GetPosition()           const { return TVector3(fX,fY,fZ);}
  Int_t     GetClusterIndex() const { return fClusterIndex;}
  Int_t     GetTopIndex()     const { return GetClusterIndex();}
  Int_t     GetBotIndex()     const { return fBotIndex;}
  Double_t  GetEloss()        const { return (fCharge * 3.61e-9);}  // 3.6 eV/Electron in Silicon
  TMatrixD  GetCov()          const { return fCov;}


//  FIXME: CAUTION The errors in the SsdHit 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);

#ifndef __CINT__ // for BOOST serialization
  template<class Archive>
  void serialize(Archive & ar, const unsigned int version)
        {
                ar & boost::serialization::base_object<FairHit>(*this);
                ar & fSensorID;
                ar & fCharge;
                ar & fNDigiHits;
                ar & fClusterIndex;
                ar & fBotIndex;
                //ar & fCov;
        }
#endif // for BOOST serialization

  virtual void Print(const Option_t* opt = 0) const;

 private:
//   Double_t GetD(Int_t i);
#ifndef __CINT__ // for BOOST serialization
                friend class boost::serialization::access;

#endif // for BOOST serialization

  //TString fDetName;  // Detector name
  Int_t fSensorID; 
  Double_t fCharge; 
  Int_t fNDigiHits; 
  Int_t fClusterIndex; 
  Int_t fBotIndex; 
  TMatrixD fCov; 
  ClassDef(PndSdsHit,8);

};

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

#endif