CbmTaskMakeRecoEvents.h

Go to the documentation of this file.

/* Copyright (C) 2021-2025 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese [committer] */
 
 
#ifndef CBMTASKMAKERECOEVENTS_H
#define CBMTASKMAKERECOEVENTS_H 1
 
 
#include "CbmDefs.h"
#include "CbmDigiEvent.h"
#include "CbmEvent.h"
#include "CbmMvdDigi.h"
 
#include <FairTask.h>
 
#include <gsl/span>
#include <vector>
 
class TClonesArray;
 

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

  CbmTaskMakeRecoEvents(const CbmTaskMakeRecoEvents&) = delete;
 

  virtual ~CbmTaskMakeRecoEvents();
 

  virtual void Exec(Option_t* opt);
 

  virtual void Finish();
 

  CbmTaskMakeRecoEvents& operator=(const CbmTaskMakeRecoEvents&) = delete;
 
 
 private:  // methods
  virtual InitStatus Init();
 

  template<typename Digi>
  void FillTree(gsl::span<const Digi> inVec, std::vector<Digi>* outVec, CbmEvent* event, ECbmDataType digiType)
  {
    size_t startIndex = outVec->size();
    size_t stopIndex  = startIndex + inVec.size();
    outVec->insert(outVec->end(), inVec.begin(), inVec.end());
    for (size_t index = startIndex; index < stopIndex; index++)
      event->AddData(digiType, index);
  }

  void FillTreeMvd(gsl::span<const CbmMvdRawDigi> inVec, std::vector<CbmMvdDigi>* outVec, CbmEvent* event)
  {
    size_t index = outVec->size();
    for (auto& rawDigi : inVec) {
      outVec->push_back(CbmMvdDigi(rawDigi));
      event->AddData(ECbmDataType::kMvdDigi, index);
      index++;
    }
  }
 
 
 private:  // members
  const std::vector<CbmDigiEvent>* fDigiEvents = nullptr;
  TClonesArray* fRecoEvents                    = nullptr;
  std::vector<CbmBmonDigi>* fBmonDigis         = nullptr;
  std::vector<CbmMvdDigi>* fMvdDigis           = nullptr;
  std::vector<CbmStsDigi>* fStsDigis           = nullptr;
  std::vector<CbmRichDigi>* fRichDigis         = nullptr;
  std::vector<CbmMuchDigi>* fMuchDigis         = nullptr;
  std::vector<CbmTrdDigi>* fTrdDigis           = nullptr;
  std::vector<CbmTofDigi>* fTofDigis           = nullptr;
  std::vector<CbmFsdDigi>* fFsdDigis           = nullptr;
  double fTimeTot                              = 0.;  
  size_t fNumTs                                = 0;   
  size_t fNumEvents                            = 0;   
 
  ClassDef(CbmTaskMakeRecoEvents, 2);
};
 
#endif /* CBMTASKMAKERECOEVENTS_H */