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