PndHypIdealRecoTask.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                 MvdIdealRecoTask header file             -----
// -----                  Created 20/03/07  by R.Kliemt           -----
// -----     Modified to HypIdealRecoTask 12.10.07 by A.Sanchez
// -------------------------------------------------------------------------

#ifndef PNDHYPIDEALRECOTASK_H
#define PNDHYPIDEALRECOTASK_H

// Root includes
#include "TVector3.h"
#include "TRandom.h"
#include "TString.h"
#include "TGeoManager.h"

// framework includes
#include "FairTask.h"

// Hyp includes
#include "PndHypPoint.h"
#include "PndHypGeoHandling.h"
class TClonesArray;

class PndHypIdealRecoTask : public FairTask {
 public:
  PndHypIdealRecoTask();

  PndHypIdealRecoTask(Double_t sx, Double_t sy, Double_t sz);

  ~PndHypIdealRecoTask();

  virtual void SetParContainers();
  virtual InitStatus Init();

  virtual void Exec(Option_t *opt);

 private:
  TString fBranchName;
  PndHypGeoHandling *fGeoH;
  TClonesArray *fPointArray;
  TClonesArray *fMctruthArray;

  TClonesArray *fHitOutputArray;
  Double_t fSigmaX; // Variance vector of sensor smearing
  Double_t fSigmaY; // Variance vector of sensor smearing
  Double_t fSigmaZ; // Variance vector of sensor smearing
  PndHypPoint *fCurrentHypPoint;
  TGeoHMatrix *fCurrentTransMat;
  TMatrixT<double> fHitCovMatrix;
  void InitTransMat();
  void smear(TVector3 &pos);
  void smearLocal(TVector3 &pos);
  void CalcDetPlane(TVector3 &o, TVector3 &u, TVector3 &v);

  void Register();

  void Reset();

  void ProduceHits();

  ClassDef(PndHypIdealRecoTask, 2);
};

#endif