R3BNeulandTacquilaReader.cxx

Go to the documentation of this file.

/******************************************************************************
 *   Copyright (C) 2019 GSI Helmholtzzentrum für Schwerionenforschung GmbH    *
 *   Copyright (C) 2019-2025 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 "R3BNeulandTacquilaReader.h"
#include "FairLogger.h"
#include "FairRootManager.h"
#include "R3BNeulandTacquilaMappedData.h"
#include "TClonesArray.h"
#include "ext_data_struct_info.hh"
 
extern "C"
{
#include "ext_data_client.h"
#include "ext_h101_raw_nnp_tacquila.h"
}
 
R3BNeulandTacquilaReader::R3BNeulandTacquilaReader(EXT_STR_h101_raw_nnp_tacquila* data, UInt_t offset)
    : R3BReader("R3BNeulandTacquilaReader")
    , fNEvent(0)
    , fData(data)
    , fOffset(offset)
    , fLogger(FairLogger::GetLogger())
    , fArray(new TClonesArray("R3BNeulandTacquilaMappedData"))
{
}
 
R3BNeulandTacquilaReader::~R3BNeulandTacquilaReader() {}
 
auto R3BNeulandTacquilaReader::Init(ext_data_struct_info* a_struct_info) -> bool
{
    // Initialize input UCESB structure
    Int_t ok;
    EXT_STR_h101_raw_nnp_tacquila_ITEMS_INFO(ok, *a_struct_info, fOffset, EXT_STR_h101_raw_nnp_tacquila, 1);
    if (!ok)
    {
        // Throw error
        perror("ext_data_struct_info_item");
        LOG(error) << "Failed to setup structure information.";
        return kFALSE;
    }
 
    // Register output array in tree
    FairRootManager::Instance()->Register("NeulandTacquilaMappedData", "Neuland", fArray, kTRUE);
 
    return kTRUE;
}
 
Bool_t R3BNeulandTacquilaReader::R3BRead()
{
    // Convert plain raw data to multi-dimensional array
    auto* data = reinterpret_cast<EXT_STR_h101_raw_nnp_tacquila_onion*>(fData);
 
    // Loop over all planes, bars and PMT's
    Int_t bar;
    Int_t barId;
    for (Int_t p = 0; p < 8; p++)
    {
        for (Int_t pmt = 0; pmt < 2; pmt++)
        {
            for (Int_t k = 0; k < data->NNP[p]._[pmt]._; k++)
            {
                // Signal channel
                UInt_t tdc1 = data->NNP[p]._[pmt].TAC[k];
                UInt_t clock1 = data->NNP[p]._[pmt].CLK[k];
                UInt_t qdc = data->NNP[p]._[pmt].ADC[k];
 
                // Stop signal (17-th channel)
                UInt_t tdc2 = data->NNP[p]._[pmt].T[k];
                // UInt_t clock2 = data->NNP[p]._[pmt].S[k];
 
                bar = data->NNP[p]._[pmt].I[k];
 
                // Calculate global bar index
                barId = p * 50 + bar;
 
                // Store data in output
                if (tdc1 && tdc2)
                {
                    new ((*fArray)[fArray->GetEntriesFast()])
                        R3BNeulandTacquilaMappedData(63 - clock1, 4095 - tdc1, 4095 - tdc2, qdc, p + 1, bar, pmt + 1);
                }
            }
        }
    }
 
    fNEvent += 1;
 
    return kTRUE;
}
 
void R3BNeulandTacquilaReader::Reset()
{
    // Reset the output array
    fArray->Clear();
}
 
ClassImp(R3BNeulandTacquilaReader)