TRestDetectorReadoutModule.h

/*************************************************************************
 * 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.                  *
 *************************************************************************/
 
#ifndef RestCore_TRestDetectorReadoutModule
#define RestCore_TRestDetectorReadoutModule
 
#include <TMath.h>
#include <TVector2.h>
 
#include <iostream>
 
#include "TRestDetectorReadoutChannel.h"
#include "TRestDetectorReadoutMapping.h"
 
class TRestDetectorReadoutModule {
   private:
    Int_t fId = -1;  
 
    std::string fName;  
 
    TVector2 fOrigin = {0, 0};  
 
    TVector2 fSize = {0, 0};  
 
    Double_t fRotation = 0;  //<
 
    std::pair<Int_t, Int_t> fDaqIdRange = {
        -1, -1};  
 
    std::vector<TRestDetectorReadoutChannel>
        fReadoutChannel;  
 
    TRestDetectorReadoutMapping fMapping;  
 
    Double_t fTolerance;  
 
    Bool_t showWarnings;  
 
    void Initialize();
 
    inline TVector2 TransformToModuleCoordinates(const TVector2& xyPhysical) const {
        auto coords = xyPhysical - fOrigin;
        TVector2 rot = coords.Rotate(-fRotation);
 
        return rot;
    }
 
    inline TVector2 TransformToPlaneCoordinates(Double_t xMod, Double_t yMod) const {
        TVector2 coords(xMod, yMod);
 
        coords = coords.Rotate(fRotation);
        coords += fOrigin;
 
        return coords;
    }
 
   protected:
   public:
    // Setters
 
    inline void SetModuleID(Int_t modID) { fId = modID; }
 
    inline void SetSize(const TVector2& size) { fSize = size; }
 
    inline void SetOrigin(const TVector2& origin) { fOrigin = origin; }
 
    inline void SetRotation(Double_t rotation) { fRotation = rotation; }
 
    inline void SetName(const std::string& name) { fName = name; }
 
    inline void SetTolerance(Double_t tolerance) { fTolerance = tolerance; }
 
    inline Double_t GetTolerance() const { return fTolerance; }
 
    inline Int_t GetMinDaqID() const { return fDaqIdRange.first; }
 
    inline Int_t GetMaxDaqID() const { return fDaqIdRange.second; }
 
    inline Int_t DaqToReadoutChannel(Int_t daqChannel) {
        for (size_t n = 0; n < GetNumberOfChannels(); n++) {
            if (GetChannel(n)->GetDaqID() == daqChannel) {
                return n;
            }
        }
        return -1;
    }
 
    inline Int_t GetModuleID() const { return fId; }
 
    inline TVector2 GetOrigin() const { return fOrigin; }
 
    inline TVector2 GetSize() const { return fSize; }
 
    inline Double_t GetRotation() const { return fRotation; }
 
    TVector2 GetModuleCoordinates(const TVector2& p) const { return TransformToModuleCoordinates(p); }
 
    TVector2 GetPlaneCoordinates(const TVector2& p) { return TransformToPlaneCoordinates(p.X(), p.Y()); }
 
    inline const char* GetName() const { return fName.c_str(); }
 
    inline TRestDetectorReadoutMapping* GetMapping() { return &fMapping; }
 
    inline TRestDetectorReadoutChannel& operator[](int n) { return fReadoutChannel[n]; }
 
    inline TRestDetectorReadoutChannel* GetChannel(size_t n) {
        if (n >= GetNumberOfChannels()) {
            return nullptr;
        }
        return &fReadoutChannel[n];
    }
 
    inline size_t GetNumberOfChannels() const { return fReadoutChannel.size(); }
 
    inline void EnableWarnings() { showWarnings = true; }
 
    inline void DisableWarnings() { showWarnings = false; }
 
    void DoReadoutMapping(Int_t nodes = 0);
 
    Bool_t isInside(const TVector2& position);
 
    inline Bool_t isInside(Double_t x, Double_t y) { return isInside({x, y}); }
 
    Bool_t isInsideChannel(Int_t channel, Double_t x, Double_t y);
    Bool_t isInsideChannel(Int_t channel, const TVector2& position);
 
    Bool_t isInsidePixel(Int_t channel, Int_t pixel, const TVector2& position);
 
    Bool_t isDaqIDInside(Int_t daqID);
    Int_t FindChannel(const TVector2& position);
    TVector2 GetDistanceToModule(const TVector2& position);
 
    TVector2 GetPixelOrigin(Int_t channel, Int_t pixel);
    TVector2 GetPixelVertex(Int_t channel, Int_t pixel, Int_t vertex);
    TVector2 GetPixelCenter(Int_t channel, Int_t pixel);
    Bool_t GetPixelTriangle(Int_t channel, Int_t pixel);
 
    TVector2 GetPixelOrigin(TRestDetectorReadoutPixel* pix);
    TVector2 GetPixelVertex(TRestDetectorReadoutPixel* pix, Int_t vertex);
    TVector2 GetPixelCenter(TRestDetectorReadoutPixel* pix);
    Bool_t GetPixelTriangle(TRestDetectorReadoutPixel* pix);
 
    TVector2 GetVertex(int n) const;
 
    void AddChannel(TRestDetectorReadoutChannel& rChannel);
 
    void SetMinMaxDaqIDs();
 
    void Draw();
 
    void Print(Int_t DetailLevel = 0);
 
    // Constructor
    TRestDetectorReadoutModule();
    // Destructor
    virtual ~TRestDetectorReadoutModule();
 
    ClassDef(TRestDetectorReadoutModule, 3);
};
#endif