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

// -------------------------------------------------------------------------
// -----                 PndSdsIdealRecoTask header file             -----
// -----                  Created 20/03/07  by R.Kliemt               -----
// -------------------------------------------------------------------------

#ifndef PNDSDSIDEALRECOTASK_H
#define PNDSDSIDEALRECOTASK_H

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

// framework includes
#include "PndSdsTask.h"

// PndSds includes
#include "PndSdsMCPoint.h"
#include "PndGeoHandling.h"

class TClonesArray;

class PndSdsIdealRecoTask : public PndSdsTask {
 public:
  PndSdsIdealRecoTask();

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

  PndSdsIdealRecoTask(PndSdsIdealRecoTask &other)
    : PndSdsTask(), fGeoH(other.fGeoH), fPointArray(other.fPointArray), fMctruthArray(other.fMctruthArray), fHitOutputArray(other.fHitOutputArray), fSigmaX(other.fSigmaX),
      fSigmaY(other.fSigmaY), fSigmaZ(other.fSigmaZ), fCurrentPndSdsMCPoint(other.fCurrentPndSdsMCPoint), fCurrentTransMat(other.fCurrentTransMat),
      fHitCovMatrix(other.fHitCovMatrix)
  {
  }

  virtual ~PndSdsIdealRecoTask();

  PndSdsIdealRecoTask &operator=(PndSdsIdealRecoTask &other)
  {
    fGeoH = other.fGeoH;
    fPointArray = other.fPointArray;
    fMctruthArray = other.fMctruthArray;
    fHitOutputArray = other.fHitOutputArray;
    fSigmaX = other.fSigmaX;
    fSigmaY = other.fSigmaY;
    fSigmaZ = other.fSigmaZ;
    fCurrentPndSdsMCPoint = other.fCurrentPndSdsMCPoint;
    fCurrentTransMat = other.fCurrentTransMat;
    fHitCovMatrix = other.fHitCovMatrix;
    return *this;
  };

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

  virtual void Exec(Option_t *opt);

 protected:
  PndGeoHandling *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
  PndSdsMCPoint *fCurrentPndSdsMCPoint;
  TGeoHMatrix *fCurrentTransMat;
  TMatrixT<Double_t> 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(PndSdsIdealRecoTask, 1);
};

#endif