PndRiemannHit.h

Go to the documentation of this file.

//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Conformal Mapping of a PndTpc Cluster (x,y)->(r,phi)-> riemann sphere
//
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Sebastian Neubert    TUM            (original author)
//
//
//-----------------------------------------------------------

#ifndef PNDRIEMANNHIT_HH
#define PNDRIEMANNHIT_HH

// Base Class Headers ----------------
#include "TObject.h"

// Collaborating Class Headers -------
#include <ostream> // remove if you do not need streaming op
#include "TVector3.h"
#include "TMatrixD.h"
#include "FairHit.h"
#include "TMath.h"

// Collaborating Class Declarations --
// class PndTpcCluster;
class PndRiemannTrack;

class PndRiemannHit : public TObject {
 public:
  // Constructors/Destructors ---------
  PndRiemannHit() : fX(0., 0., 0.), fSigmaX(0., 0., 0.), fCovX(3, 3), fHit(0), fHitID(-1), fS(-1.), fZ(-1.), fDeltaZ(-1.), fAlpha(-1.), fVerbose(0){};
  PndRiemannHit(double x, double y, double z, double dx, double dy, double dz);
  PndRiemannHit(FairHit *cl, int hitID = -1);
  ~PndRiemannHit();

  // Copy-/Assignment-Operator
  PndRiemannHit(const PndRiemannHit &myHit)
    : TObject(myHit), fX(myHit.fX), fSigmaX(myHit.fSigmaX), fCovX(myHit.fCovX), fHit(myHit.fHit), fHitID(myHit.fHitID), fS(myHit.fS), fZ(myHit.fZ), fDeltaZ(myHit.fDeltaZ),
      fAlpha(myHit.fAlpha), fVerbose(myHit.fVerbose)
  {
  }

  PndRiemannHit &operator=(const PndRiemannHit &myHit)
  {
    if (this != &myHit) {
      fX = myHit.fX;
      fSigmaX = myHit.fSigmaX;
      fCovX = myHit.fCovX;
      fHit = myHit.fHit;
      fHitID = myHit.fHitID;
      fS = myHit.fS;
      fZ = myHit.fZ;
      fDeltaZ = myHit.fDeltaZ;
      fAlpha = myHit.fAlpha;
      fVerbose = myHit.fVerbose;
    }
    return *this;
  }

  // Accessors -----------------------
  void setXYZ(double x, double y, double z);
  void setDXYZ(double dx, double dy, double dz);
  void setHit(FairHit *cl);
  const TVector3 &x() const { return fX; }
  const FairHit *hit() const { return fHit; }
  int hitID() const { return fHitID; }
  double s() const { return fS; }
  double z() const;
  double alpha() const { return fAlpha; }
  double sigmaXY() const;
  double sigmaX() const { return TMath::Sqrt(fCovX[0][0]); }
  double sigmaY() const { return TMath::Sqrt(fCovX[1][1]); }
  double sigmaW() const { return TMath::Sqrt(TMath::Power(2 * x().x() * sigmaX(), 2) + TMath::Power(2 * x().y() * sigmaY(), 2)); } // error of z-Coordinate (x2+y2) in RiemannSpace
  const TMatrixD &covX() const { return fCovX; }
  double covX(int row, int col) const { return fCovX[row][col]; }
  bool operator<(const PndRiemannHit &aHit) const
  { 
    if (s() >= 0 && aHit.s() >= 0) {
      return s() < aHit.s();
    } else if (s() <= 0 && aHit.s() <= 0) {
      return -s() < -aHit.s();
    } else
      return s() < aHit.s();
  }

  // Modifiers -----------------------
  // Operations ----------------------
  void calcPosOnTrk(PndRiemannTrack *trk);

 private:
  // Private Data Members ------------
  TVector3 fX;
  TVector3 fSigmaX; 
  TMatrixD fCovX;   
  FairHit *fHit;    //->  //no ownership over this pointer!
  int fHitID;
  double fS;      
  double fZ;      
  double fDeltaZ; 
  double fAlpha;  
  int fVerbose;

  // Private Methods -----------------

 public:
  ClassDef(PndRiemannHit, 2)
};

#endif

//--------------------------------------------------------------
// $Log$
//--------------------------------------------------------------