GeaneTrackRep.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".              *
//****************************************************************************

//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      a GEANE (sd-system)  track representation
//      (q/p, v',w',v,w)
//      (v,w) refers to DetPlane system
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Sebastian Neubert    TUM            (original author)
//      Andrea Fontana       INFN
//
//
//-----------------------------------------------------------

#ifndef GeaneTRACKREP_HH
#define GeaneTRACKREP_HH

// Base Class Headers ----------------
#include "GFAbsTrackRep.h"
#include "FairTrackParP.h"

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

// Collaborating Class Declarations --
class FairGeanePro;

class GeaneTrackRep : public GFAbsTrackRep {
 public:
  // Constructors/Destructors ---------
  GeaneTrackRep();
  GeaneTrackRep(FairGeanePro *geane,
                const GFDetPlane &plane, // will be defined at origin of plane
                const TVector3 &mom, const TVector3 &poserr, const TVector3 &momerr, double q, int PDGCode);

  GeaneTrackRep(FairGeanePro *geane,
                const GFDetPlane &plane, // will be defined at origin of plane
                const TVector3 &mom, const TVector3 &poserr, const TVector3 &momerr, int q, int PDGCode);

  virtual ~GeaneTrackRep();

  virtual GFAbsTrackRep *clone() const { return new GeaneTrackRep(*this); }
  virtual GFAbsTrackRep *prototype() const { return new GeaneTrackRep(); }

  // Operators
  friend std::ostream &operator<<(std::ostream &s, const GeaneTrackRep &me);

  // Accessors -----------------------

  // Modifiers

  // Operations ----------------------

  virtual double extrapolate(const GFDetPlane &, TMatrixT<double> &statePred);
  // virtual void extrapolate(const GFDetPlane&,
  //               const TMatrixT<double>& stateFrom
  //               TMatrixT<double>& stateResult);

  virtual double extrapolate(const GFDetPlane &, TMatrixT<double> &statePred, TMatrixT<double> &covPred);

  // these two are overwritting ABsTrackRep methods
  void extrapolateToPoint(const TVector3 &pos, TVector3 &poca, TVector3 &dirInPoca);

  void extrapolateToLine(const TVector3 &point1, const TVector3 &point2, TVector3 &poca, TVector3 &dirInPoca, TVector3 &poca_onwire);

  TVector3 getPocaOnLine(const TVector3 &p1, const TVector3 &p2, bool back = false);

  virtual TVector3 getPos(const GFDetPlane &);
  virtual TVector3 getMom(const GFDetPlane &);
  virtual void getPosMom(const GFDetPlane &, TVector3 &pos, TVector3 &mom);
  virtual void getPosMomCov(const GFDetPlane &pl, TVector3 &pos, TVector3 &mom, TMatrixT<double> &cov);
  virtual double getCharge() const { return fState[0][0] > 0 ? 1. : -1.; }
  int getPropDir() { return _backw; }
  FairGeanePro *getPropagator() { return _geane; }
  int getPDG() { return _pdg; };
  double getSPU() { return _spu; }

  void setPropagator(FairGeanePro *g) { _geane = g; }
  void setPropDir(int d) { _backw = d; }
  void switchDirection() { _backw = -_backw; }

  // (-1,0,1) -> (backward prop,decide myself,forward)

 private:
  void checkState(); // checks if state vector is inside numerical limits

  // Private Data Members ------------
  FairGeanePro *_geane; 
  int _pdg;             // pdg code of the particle to be tracked
  int _backw;           // (-1,0,1) -> (backward prop,decide myself,forward)
  double _spu;          // sign of z-component of momentum

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

  // calculate jacobian of extrapolation
  // void Jacobian(const GFDetPlane& pl,
  //    const TMatrixT<double>& statePred,
  //    TMatrixT<double>& jacResult);

 public:
  ClassDef(GeaneTrackRep, 1)
};

#endif

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