TRestDetectorReadoutPixel.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 "TRestDetectorReadoutPixel.h"
 
#include <TRestMetadata.h>
 
using namespace std;
 
ClassImp(TRestDetectorReadoutPixel);
TRestDetectorReadoutPixel::TRestDetectorReadoutPixel() { Initialize(); }
 
TRestDetectorReadoutPixel::~TRestDetectorReadoutPixel() {}
 
void TRestDetectorReadoutPixel::Initialize() {}
 
TVector2 TRestDetectorReadoutPixel::GetCenter() const {
    TVector2 center(0, 0);
    TVector2 origin(fPixelOriginX, fPixelOriginY);
    TVector2 oppositeVertex = GetVertex(2);
 
    if (fTriangle) {
        center = (oppositeVertex - origin) / 4. + origin;
    } else {
        center = (origin + oppositeVertex) / 2.;
    }
    return center;
}
 
TVector2 TRestDetectorReadoutPixel::GetVertex(int n) const {
    TVector2 vertex(0, 0);
    TVector2 origin(fPixelOriginX, fPixelOriginY);
 
    if (n % 4 == 0)
        return origin;
    else if (n % 4 == 1) {
        vertex.Set(fPixelSizeX, 0);
        vertex = vertex.Rotate(fRotation * TMath::Pi() / 180.);
 
        vertex = vertex + origin;
    } else if (n % 4 == 2) {
        vertex.Set(fPixelSizeX, fPixelSizeY);
        vertex = vertex.Rotate(fRotation * TMath::Pi() / 180.);
 
        vertex = vertex + origin;
    } else if (n % 4 == 3) {
        vertex.Set(0, fPixelSizeY);
        vertex = vertex.Rotate(fRotation * TMath::Pi() / 180.);
 
        vertex = vertex + origin;
    }
    return vertex;
}
 
Bool_t TRestDetectorReadoutPixel::IsInside(const TVector2& inputPosition) {
    const auto pos = TransformToPixelCoordinates(inputPosition);
    Double_t const x = pos.X();
    if (pos.X() >= -fTolerance && pos.X() <= fPixelSizeX + fTolerance)  // Condition on X untouched
    {
        if (fTriangle && pos.Y() >= -fTolerance &&
            pos.Y() <= fPixelSizeY + fTolerance -
                           x * (fPixelSizeY / fPixelSizeX))  // if triangle, third condition depends on x
            return true;
        if (!fTriangle && pos.Y() >= -fTolerance &&
            pos.Y() <= fPixelSizeY + fTolerance)  // for a normal rectangular pixel, same
                                                  // simple conditions
            return true;
    }
    return false;
}
 
TVector2 TRestDetectorReadoutPixel::TransformToPixelCoordinates(const TVector2& pixel) const {
    TVector2 pos(pixel.X() - fPixelOriginX, pixel.Y() - fPixelOriginY);
    pos = pos.Rotate(-fRotation * TMath::Pi() / 180.);
    return pos;
}
 
void TRestDetectorReadoutPixel::Print() const {
    std::cout << "    ## Pixel  position : (" << GetOriginX() << "," << GetOriginY() << ") mm size : ("
              << GetSizeX() << "," << GetSizeY() << ") mm" << std::endl;
    std::cout << "       rotation : " << fRotation << " degrees"
              << " type : ";
    if (fTriangle)
        std::cout << "triangle" << std::endl;
    else
        std::cout << "rectangle" << std::endl;
}