TRestAxionEventProcess.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 "TRestAxionEventProcess.h"
 
using namespace std;
 
ClassImp(TRestAxionEventProcess);
 
TRestAxionEventProcess::TRestAxionEventProcess() { fSingleThreadOnly = false; }
 
TRestAxionEventProcess::~TRestAxionEventProcess() {}
 
void TRestAxionEventProcess::BeginOfEventProcess(TRestEvent* inEv) {
    TRestEventProcess::BeginOfEventProcess(inEv);
 
    fAxionEvent = (TRestAxionEvent*)inEv;
    RESTDebug << "BoEP: Initial Position. X: " << fAxionEvent->GetPosition().X()
              << " Y: " << fAxionEvent->GetPosition().Y() << " Z: " << fAxionEvent->GetPosition().Z()
              << RESTendl;
    RESTDebug << "BoEP: Initial Direction. X: " << fAxionEvent->GetDirection().X()
              << " Y: " << fAxionEvent->GetDirection().Y() << " Z: " << fAxionEvent->GetDirection().Z()
              << RESTendl;
    fAxionEvent->RotateYX(fPosition, fInternalYaw, fInternalPitch);
    fAxionEvent->RotateYX(fExternalRotationCenter, fExternalYaw, fExternalPitch);
    fAxionEvent->Translate(TVector3(-fPosition.X(), -fPosition.Y(), -fPosition.Z()));
    RESTDebug << " ---- " << RESTendl;
    RESTDebug << "BoEP: Final Position. X: " << fAxionEvent->GetPosition().X()
              << " Y: " << fAxionEvent->GetPosition().Y() << " Z: " << fAxionEvent->GetPosition().Z()
              << RESTendl;
    RESTDebug << "BoEP: Final Direction. X: " << fAxionEvent->GetDirection().X()
              << " Y: " << fAxionEvent->GetDirection().Y() << " Z: " << fAxionEvent->GetDirection().Z()
              << RESTendl;
    RESTDebug << " ++++ " << RESTendl;
}
 
void TRestAxionEventProcess::EndOfEventProcess(TRestEvent* evInput) {
    TRestEventProcess::EndOfEventProcess(evInput);
 
    RESTDebug << "EoEP: Initial Position. X: " << fAxionEvent->GetPosition().X()
              << " Y: " << fAxionEvent->GetPosition().Y() << " Z: " << fAxionEvent->GetPosition().Z()
              << RESTendl;
    RESTDebug << "EoEP: Initial Direction. X: " << fAxionEvent->GetDirection().X()
              << " Y: " << fAxionEvent->GetDirection().Y() << " Z: " << fAxionEvent->GetDirection().Z()
              << RESTendl;
    RESTDebug << " ---- " << RESTendl;
    fAxionEvent->Translate(TVector3(fPosition.X(), fPosition.Y(), fPosition.Z()));
    if (fLocalAxis) {
        fAxionEvent->RotateXY(fExternalRotationCenter, -fExternalPitch, -fExternalYaw);
        fAxionEvent->RotateXY(fPosition, -fInternalPitch, -fInternalYaw);
    }
    RESTDebug << "EoEP: Final Position. X: " << fAxionEvent->GetPosition().X()
              << " Y: " << fAxionEvent->GetPosition().Y() << " Z: " << fAxionEvent->GetPosition().Z()
              << RESTendl;
    RESTDebug << "EoEP: Final Direction. X: " << fAxionEvent->GetDirection().X()
              << " Y: " << fAxionEvent->GetDirection().Y() << " Z: " << fAxionEvent->GetDirection().Z()
              << RESTendl;
    RESTDebug << " ++++ " << RESTendl;
}
 
void TRestAxionEventProcess::BeginPrintProcess() {
    TRestEventProcess::BeginPrintProcess();
 
    RESTMetadata << "Position: (" << fPosition.X() << ", " << fPosition.Y() << ", " << fPosition.Z() << ") mm"
                 << RESTendl;
    RESTMetadata << " " << RESTendl;
    RESTMetadata << "Internal Yaw angle (Y-axis): " << fInternalYaw * units("degrees") << " degrees"
                 << RESTendl;
    RESTMetadata << "Internal Pitch angle (X-axis): " << fInternalPitch * units("degrees") << " degrees"
                 << RESTendl;
    RESTMetadata << " " << RESTendl;
    RESTMetadata << "External rotation center: (" << fExternalRotationCenter.X() << ", "
                 << fExternalRotationCenter.Y() << ", " << fExternalRotationCenter.Z() << ") mm" << RESTendl;
    RESTMetadata << "External Yaw angle (Y-axis): " << fExternalYaw * units("degrees") << " degrees"
                 << RESTendl;
    RESTMetadata << "External Pitch angle (X-axis): " << fExternalPitch * units("degrees") << " degrees"
                 << RESTendl;
    RESTMetadata << " " << RESTendl;
    if (fLocalAxis)
        RESTMetadata << "Local axis is true " << RESTendl;
    else
        RESTMetadata << "Universal axis is true " << RESTendl;
    RESTMetadata << " --------------------------- " << RESTendl;
    RESTMetadata << " " << RESTendl;
}
 
void TRestAxionEventProcess::EndPrintProcess() { TRestEventProcess::EndPrintProcess(); }