PndSdsHit.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

// -------------------------------------------------------------------------
// -----                    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>

#ifndef __CINT__ // for BOOST serialization
#include <boost/serialization/serialization.hpp>
#include <boost/serialization/vector.hpp>
#endif

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 = nullptr) 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