R3BNeulandCalToHitParTask.cxx

Go to the documentation of this file.

/******************************************************************************
 *   Copyright (C) 2019 GSI Helmholtzzentrum für Schwerionenforschung GmbH    *
 *   Copyright (C) 2019-2023 Members of R3B Collaboration                     *
 *                                                                            *
 *             This software is distributed under the terms of the            *
 *                 GNU General Public Licence (GPL) version 3,                *
 *                    copied verbatim in the file "LICENSE".                  *
 *                                                                            *
 * In applying this license GSI does not waive the privileges and immunities  *
 * granted to it by virtue of its status as an Intergovernmental Organization *
 * or submit itself to any jurisdiction.                                      *
 ******************************************************************************/
 
#include "R3BNeulandCalToHitParTask.h"
#include "R3BDataMonitor.h"
#include "R3BException.h"
#include "R3BNeulandBasePar.h"
#include "R3BNeulandCalToHitPar.h"
#include "R3BNeulandCalibrationTask.h"
#include "R3BNeulandCommon.h"
#include <FairRootManager.h>
#include <FairRuntimeDb.h>
#include <R3BLogger.h>
#include <R3BNeulandLSQREngineAdaptor.h>
#include <R3BNeulandMillepede.h>
#include <R3BNeulandPredecessor.h>
#include <memory>
#include <string_view>
 
namespace R3B::Neuland
{
    // NOLINTBEGIN
    Cal2HitParTask::Cal2HitParTask(Cal2HitParMethod method,
                                   std::string_view cal_data_name,
                                   std::string_view base_par_name,
                                   std::string_view hit_par_name,
                                   std::string_view name,
                                   int iVerbose)
        : CalibrationTask(name, iVerbose)
        , cal_data_{ cal_data_name }
        , base_par_{ InputPar<CalibrationBasePar>(base_par_name) }
    // NOLINTEND
    {
        switch (method)
        {
            case Cal2HitParMethod::LSQT:
                R3BLOG(info, "Cal2HitPar method: LSQT.");
                engine_ = std::make_unique<Calibration::LSQREngineAdaptor>();
                hit_par_ = OutputPar<Cal2HitPar>(hit_par_name);
                break;
            case Cal2HitParMethod::millepede:
                R3BLOG(info, "Cal2HitPar method: Millepede.");
                engine_ = std::make_unique<Calibration::MillepedeEngine>();
                hit_par_ = InputPar<Cal2HitPar>(hit_par_name);
                break;
            case Cal2HitParMethod::predecessor:
                R3BLOG(info, "Cal2HitPar method: predecessor.");
                engine_ = std::make_unique<Calibration::Predecessor>();
                hit_par_ = OutputPar<Cal2HitPar>(hit_par_name);
                break;
        }
    }
 
    void Cal2HitParTask::HistogramInit(DataMonitor& histograms) { engine_->HistInit(histograms); }
 
    void Cal2HitParTask::ExtraInit(FairRootManager* /*rootMan*/)
    {
        cal_data_.init();
 
        const auto plane_num = base_par_->GetNumOfPlanes();
        if (plane_num == 0)
        {
            throw R3B::runtime_error("Plane number extracted from Map2CalPar is 0!");
        }
        engine_->SetTask(this);
        engine_->SetModuleSize(plane_num * BarsPerPlane);
        engine_->Init();
    }
 
    void Cal2HitParTask::SetExtraPar(FairRuntimeDb* rtdb) {}
 
    void Cal2HitParTask::TriggeredExec()
    {
        engine_->EventReset();
        engine_->AddSignals(cal_data_.get());
        // for (const auto& bar_signal : cal_data_)
        // {
        //     engine_->AddSignal(bar_signal);
        // }
        auto* eventHeader = GetEventHeader();
        if (eventHeader != nullptr)
        {
            engine_->EndOfEvent(eventHeader->GetEventno());
        }
        else
        {
            throw R3B::runtime_error("Event header is nullptr!");
        }
    }
 
    void Cal2HitParTask::EndOfTask()
    {
        if (hit_par_ == nullptr)
        {
            throw R3B::runtime_error("Hit parameter is nullptr during the calibration!");
        }
 
        R3BLOG(info, "Starting to do cosmic calibration calibration!");
        engine_->Calibrate(*hit_par_);
        engine_->EndOfTask();
    }
 
    auto Cal2HitParTask::CheckConditions() const -> bool
    {
        auto res = engine_->SignalFilter(cal_data_.get());
        return res;
    }
} // namespace R3B::Neuland