TRestDetectorHitsSmearingProcess.cxx

/*************************************************************************
 * This file is part of the REST software framework.                     *
 *                                                                       *
 * Copyright (C) 2016 GIFNA/TREX (University of Zaragoza)                *
 * For more information see http://gifna.unizar.es/trex                  *
 *                                                                       *
 * REST is free software: you can redistribute it and/or modify          *
 * it under the terms of the GNU General Public License as published by  *
 * the Free Software Foundation, either version 3 of the License, or     *
 * (at your option) any later version.                                   *
 *                                                                       *
 * REST is distributed in the hope that it will be useful,               *
 * but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the          *
 * GNU General Public License for more details.                          *
 *                                                                       *
 * You should have a copy of the GNU General Public License along with   *
 * REST in $REST_PATH/LICENSE.                                           *
 * If not, see http://www.gnu.org/licenses/.                             *
 * For the list of contributors see $REST_PATH/CREDITS.                  *
 *************************************************************************/
 
 
#include "TRestDetectorHitsSmearingProcess.h"
 
using namespace std;
 
ClassImp(TRestDetectorHitsSmearingProcess);
 
TRestDetectorHitsSmearingProcess::TRestDetectorHitsSmearingProcess() { Initialize(); }
 
TRestDetectorHitsSmearingProcess::~TRestDetectorHitsSmearingProcess() { delete fOutputEvent; }
 
void TRestDetectorHitsSmearingProcess::Initialize() {
    SetSectionName(this->ClassName());
    SetLibraryVersion(LIBRARY_VERSION);
 
    fInputEvent = nullptr;
    fOutputEvent = new TRestDetectorHitsEvent();
}
 
void TRestDetectorHitsSmearingProcess::InitProcess() {
    delete fRandom;
    fRandom = new TRandom3(fSeed);
    fSeed = fRandom->TRandom::GetSeed();
 
    fReadout = dynamic_cast<TRestDetectorReadout*>(fRunInfo->GetMetadataClass("TRestDetectorReadout"));
 
    if (!fReadout && !fChannelType.empty()) {
        throw runtime_error(
            "TRestDetectorHitsSmearingProcess::InitProcess() : "
            "TRestDetectorReadout not found and channel type is requested");
    }
 
    transform(fChannelType.begin(), fChannelType.end(), fChannelType.begin(), ::tolower);
    const auto validChannelTypes = {"veto", "tpc"};
    if (find(validChannelTypes.begin(), validChannelTypes.end(), fChannelType) == validChannelTypes.end()) {
        throw runtime_error(
            "TRestDetectorHitsSmearingProcess::InitProcess() : Channel type not valid: valid types are: "
            "'tpc', 'veto'");
    }
}
 
TRestEvent* TRestDetectorHitsSmearingProcess::ProcessEvent(TRestEvent* inputEvent) {
    fInputEvent = (TRestDetectorHitsEvent*)inputEvent;
    fOutputEvent->SetEventInfo(fInputEvent);
 
    Double_t eDep = 0.0;
    for (int hit = 0; hit < static_cast<int>(fInputEvent->GetNumberOfHits()); hit++) {
        const auto hitType = fInputEvent->GetType(hit);
        if ((fChannelType == "veto" && hitType == REST_HitType::VETO) ||
            (fChannelType == "tpc" && hitType != REST_HitType::VETO)) {
            eDep += fInputEvent->GetEnergy(hit);
        }
    }
 
    Double_t gain = 1.0;
    if (eDep > 0) {
        Double_t eRes = fResolutionAtERef * TMath::Sqrt(fEnergyRef / eDep) / 2.35 / 100.0;
        gain = fRandom->Gaus(1.0, eRes);
    }
 
    for (int hit = 0; hit < static_cast<int>(fInputEvent->GetNumberOfHits()); hit++) {
        const auto hitType = fInputEvent->GetType(hit);
 
        // VETO hit types are processed according the "veto" channel type. All other hit types are processed
        // according to the "tpc" type (XYZ, XZ, etc.)
 
        if (fChannelType == "veto") {
            if (hitType == REST_HitType::VETO) {
                fOutputEvent->AddHit(fInputEvent->GetX(hit), fInputEvent->GetY(hit), fInputEvent->GetZ(hit),
                                     fInputEvent->GetEnergy(hit) * gain, fInputEvent->GetTime(hit),
                                     fInputEvent->GetType(hit));
            } else {
                // Do nothing
                fOutputEvent->AddHit(fInputEvent->GetX(hit), fInputEvent->GetY(hit), fInputEvent->GetZ(hit),
                                     fInputEvent->GetEnergy(hit), fInputEvent->GetTime(hit),
                                     fInputEvent->GetType(hit));
            }
        } else if (fChannelType == "tpc") {
            if (hitType == REST_HitType::VETO) {
                // Do nothing
                fOutputEvent->AddHit(fInputEvent->GetX(hit), fInputEvent->GetY(hit), fInputEvent->GetZ(hit),
                                     fInputEvent->GetEnergy(hit), fInputEvent->GetTime(hit),
                                     fInputEvent->GetType(hit));
            } else {
                fOutputEvent->AddHit(fInputEvent->GetX(hit), fInputEvent->GetY(hit), fInputEvent->GetZ(hit),
                                     fInputEvent->GetEnergy(hit) * gain, fInputEvent->GetTime(hit),
                                     fInputEvent->GetType(hit));
            }
        } else {
            // this should never happen
            throw runtime_error(
                "TRestDetectorHitsSmearingProcess::ProcessEvent() : "
                "Channel type not valid: valid types are: "
                "'tpc', 'veto'");
        }
    }
 
    return fOutputEvent;
}