TRestGridMask.cxx

/*************************************************************************
 * This file is part of the REST software framework.                     *
 *                                                                       *
 * Copyright (C) 2016 GIFNA/TREX (University of Zaragoza)                *
 * For more information see https://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 https://www.gnu.org/licenses/.                            *
 * For the list of contributors see $REST_PATH/CREDITS.                  *
 *************************************************************************/
 
 
#include "TRestGridMask.h"
 
#include "TRandom3.h"
 
ClassImp(TRestGridMask);
 
TRestGridMask::TRestGridMask() : TRestPatternMask() { Initialize(); }
 
TRestGridMask::TRestGridMask(const char* cfgFileName, std::string name) : TRestPatternMask(cfgFileName) {
    Initialize();
 
    LoadConfigFromFile(fConfigFileName, name);
 
    if (GetVerboseLevel() >= TRestStringOutput::REST_Verbose_Level::REST_Info) PrintMetadata();
}
 
TRestGridMask::~TRestGridMask() {}
 
void TRestGridMask::Initialize() {
    SetSectionName(this->ClassName());
    SetType("Grid");
}
 
Int_t TRestGridMask::GetRegion(Double_t& x, Double_t& y) {
    if (TRestPatternMask::GetRegion(x, y)) return 0;
 
    Double_t xEval = fGridThickness / 2. + x;
 
    Int_t xcont = 0;
    if (xEval > 0) {
        while (xEval > fGridGap) {
            xEval -= fGridGap;
            xcont++;
        }
    } else {
        while (xEval < 0) {
            xEval += fGridGap;
            xcont--;
        }
    }
 
    if (xEval < fGridThickness) return 0;
 
    Double_t yEval = fGridThickness / 2. + y;
 
    Int_t ycont = 0;
    if (yEval > 0) {
        while (yEval > fGridGap) {
            yEval -= fGridGap;
            ycont++;
        }
    } else {
        while (yEval < 0) {
            yEval += fGridGap;
            ycont--;
        }
    }
 
    if (yEval < fGridThickness) return 0;
 
    xcont = xcont % fModulus;
    if (xcont < 0) xcont += fModulus;
 
    ycont = ycont % fModulus;
    if (ycont < 0) ycont += fModulus;
 
    return 1 + (fModulus * ycont + xcont) % fMaxRegions;
}
 
void TRestGridMask::PrintMetadata() {
    TRestPatternMask::PrintMetadata();
 
    PrintMaskMembers();
    RESTMetadata << "----" << RESTendl;
}
 
void TRestGridMask::PrintMask() {
    PrintCommonPatternMembers();
    RESTMetadata << "----" << RESTendl;
    PrintMaskMembers();
}
 
void TRestGridMask::PrintMaskMembers() {
    RESTMetadata << " - Grid gap : " << fGridGap << " mm" << RESTendl;
    RESTMetadata << " - Grid thickness : " << fGridThickness << " mm" << RESTendl;
}