TRestDetectorHitsSpecularProcess.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 "TRestDetectorHitsSpecularProcess.h"
 
using namespace std;
 
ClassImp(TRestDetectorHitsSpecularProcess);
 
TRestDetectorHitsSpecularProcess::TRestDetectorHitsSpecularProcess() { Initialize(); }
 
TRestDetectorHitsSpecularProcess::TRestDetectorHitsSpecularProcess(const char* configFilename) {
    Initialize();
    LoadConfigFromFile(configFilename);
}
 
TRestDetectorHitsSpecularProcess::~TRestDetectorHitsSpecularProcess() {}
 
void TRestDetectorHitsSpecularProcess::Initialize() {
    SetSectionName(this->ClassName());
    SetLibraryVersion(LIBRARY_VERSION);
 
    fInputEvent = nullptr;
    fOutputEvent = new TRestDetectorHitsEvent();
}
 
TRestEvent* TRestDetectorHitsSpecularProcess::ProcessEvent(TRestEvent* inputEvent) {
    fInputEvent = (TRestDetectorHitsEvent*)inputEvent;
    fOutputEvent->SetEventInfo(fInputEvent);
 
    for (unsigned int hit = 0; hit < fInputEvent->GetNumberOfHits(); hit++) {
        TVector3 position(fInputEvent->GetX(hit), fInputEvent->GetY(hit), fInputEvent->GetZ(hit));
        const auto type = fInputEvent->GetType(hit);
        const auto energy = fInputEvent->GetEnergy(hit);
        const auto time = fInputEvent->GetTime(hit);
 
        if (type != XYZ) {
            fOutputEvent->AddHit(position, energy, time, type);
            continue;
        }
 
        const TVector3 relativePosition = position - fPosition;
        const TVector3 reflectionVector = relativePosition - 2 * relativePosition.Dot(fNormal) * fNormal;
 
        position = fPosition + reflectionVector;
 
        fOutputEvent->AddHit(position, energy, time, type);
    }
    return fOutputEvent;
}
 
void TRestDetectorHitsSpecularProcess::InitFromConfigFile() {
    TRestEventProcess::InitFromConfigFile();
 
    fNormal = Get3DVectorParameterWithUnits("normal", {0, 0, 1});
    fNormal = fNormal.Unit();
    fPosition = Get3DVectorParameterWithUnits("position", {0, 0, 1});
}
 
void TRestDetectorHitsSpecularProcess::PrintMetadata() {
    BeginPrintProcess();
 
    RESTMetadata << " - Plane normal : ( " << fNormal.X() << ", " << fNormal.Y() << ", " << fNormal.Z()
                 << ") mm" << RESTendl;
    RESTMetadata << " - Plane position : ( " << fPosition.X() << ", " << fPosition.Y() << ", "
                 << fPosition.Z() << ")" << RESTendl;
 
    EndPrintProcess();
}