CbmStsHitProducerIdealAlgo.cxx

/* Copyright (C) 2019-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Florian Uhlig [committer] */

// -----------------------------------------------------------------------------
// -----                                                                   -----
// -----                     CbmStsHitProducerIdealAlgo                    -----
// -----                                                                   -----
// -----------------------------------------------------------------------------

#include "CbmStsHitProducerIdealAlgo.h"

#include "CbmTrdParSetGas.h"

#include <Logger.h>

#include <TList.h>

// -------------------------------------------------------------------------
CbmStsHitProducerIdealAlgo::CbmStsHitProducerIdealAlgo() : CbmAlgo() {}

CbmStsHitProducerIdealAlgo::~CbmStsHitProducerIdealAlgo()
{
  /// Clear buffers
}

// -------------------------------------------------------------------------
Bool_t CbmStsHitProducerIdealAlgo::Init()
{
  LOG(info) << "Initializing tutorial StsHitProducerIdeal algo";

  return kTRUE;
}

void CbmStsHitProducerIdealAlgo::Reset() {}

void CbmStsHitProducerIdealAlgo::Finish() {}

// -------------------------------------------------------------------------
Bool_t CbmStsHitProducerIdealAlgo::InitContainers()
{
  LOG(info) << "Init parameter containers for CbmStsHitProducerIdealAlgo";
  Bool_t initOK = ReInitContainers();

  return initOK;
}
Bool_t CbmStsHitProducerIdealAlgo::ReInitContainers()
{
  LOG(info) << "**********************************************";
  LOG(info) << "ReInit parameter containers for CbmStsHitProducerIdealAlgo";

  fTrdGasPar = static_cast<CbmTrdParSetGas*>(fParCList->FindObject("CbmTrdParSetGas"));
  if (nullptr == fTrdGasPar) return kFALSE;

  fTrdGasPar->Dump();
  Bool_t initOK = InitParameters();

  return initOK;
}

TList* CbmStsHitProducerIdealAlgo::GetParList()
{

  if (nullptr == fParCList) { fParCList = new TList(); }

  fTrdGasPar = new CbmTrdParSetGas("CbmTrdParSetGas");
  fParCList->Add(fTrdGasPar);

  return fParCList;
}

Bool_t CbmStsHitProducerIdealAlgo::InitParameters() { return kTRUE; }
// -------------------------------------------------------------------------

std::vector<CbmStsHit> CbmStsHitProducerIdealAlgo::ProcessInputData(const std::vector<CbmStsPoint>& pointVect)
{
  fTrdGasPar->Print();

  // Declare some variables
  //  CbmStsPoint* point{nullptr};
  Int_t detID {0};  // Detector ID
  Double_t x {0.};
  Double_t y {0.};
  Double_t z {0.1};      // Position
  Double_t dx {0.0001};  // Position error
  TVector3 pos {};
  TVector3 dpos {};  // Position and error vectors

  std::vector<CbmStsHit> hitVect {};


  //  for(auto point: pointVect) {
  for (unsigned long iPoint = 0; iPoint < pointVect.size(); ++iPoint) {

    // Detector ID
    detID = pointVect.at(iPoint).GetDetectorID();

    // Determine hit position (centre plane of station)
    x = 0.5 * (pointVect.at(iPoint).GetXOut() + pointVect.at(iPoint).GetXIn());
    y = 0.5 * (pointVect.at(iPoint).GetYOut() + pointVect.at(iPoint).GetYIn());
    z = 0.5 * (pointVect.at(iPoint).GetZOut() + pointVect.at(iPoint).GetZIn());

    // Create new hit
    pos.SetXYZ(x, y, z);
    dpos.SetXYZ(dx, dx, 0.);

    hitVect.emplace_back(detID, pos, dpos, 0., iPoint, iPoint, 0., 0.);
  }  // Loop over MCPoints

  return hitVect;
}

// -------------------------------------------------------------------------