TRestDetectorReadoutPlane.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_TRestDetectorReadoutPlane
#define RestCore_TRestDetectorReadoutPlane
 
#include <TGraph.h>
#include <TH2Poly.h>
 
#include <iostream>
 
#include "TRestDetectorReadoutChannel.h"
#include "TRestDetectorReadoutModule.h"
 
class TRestDetectorReadoutPlane {
   private:
    Int_t fId = -1;  //<
 
    TVector3 fPosition = {0, 0, 0};  //<
 
    TVector3 fNormal = {0, 0, 1};  //<
 
    TVector3 fAxisX;  //<
 
    TVector3 fAxisY;  //<
 
    Double_t fChargeCollection = 1;  //<
 
    Double_t fHeight = 0;  //<
 
    Double_t fRotation = 0;  //<
 
    std::vector<TRestDetectorReadoutModule> fReadoutModules;  //<
 
    void UpdateAxes();
 
   public:
    // Setters
    void SetID(int id) { fId = id; }
 
    void SetPosition(const TVector3& position) { fPosition = position; }
 
    void SetNormal(const TVector3& normal);
 
    void SetChargeCollection(Double_t charge) { fChargeCollection = charge; }
 
    void SetHeight(Double_t d);
 
    void SetRotation(Double_t radians);
 
    void SetAxisX(const TVector3& axis);
 
    // Getters
    inline Int_t GetID() const { return fId; }
 
    inline TVector3 GetPosition() const { return fPosition; }
 
    Double_t GetRotation() const { return fRotation; }
 
    inline TVector3 GetCathodePosition() const { return fPosition + fHeight * fNormal; }
 
    inline TVector3 GetNormal() const { return fNormal; }
 
    inline TVector3 GetAxisX() const { return fAxisX; }
 
    inline TVector3 GetAxisY() const { return fAxisY; }
 
    inline Double_t GetChargeCollection() const { return fChargeCollection; }
 
    inline Double_t GetDriftDistance() const { return fHeight; }
 
    inline Double_t GetHeight() const { return fHeight; }
 
    Double_t GetDistanceTo(const TVector3& pos) const;
 
    TVector2 GetDistanceToModule(Int_t mod, const TVector2& position) {
        return GetModuleByID(mod)->GetDistanceToModule(position);
    }
 
    TRestDetectorReadoutModule& operator[](int mod) { return fReadoutModules[mod]; }
 
    TRestDetectorReadoutModule* GetModule(size_t mod) {
        if (mod >= GetNumberOfModules()) {
            return nullptr;
        }
        return &fReadoutModules[mod];
    }
 
    size_t GetNumberOfModules() { return fReadoutModules.size(); }
 
    void AddModule(TRestDetectorReadoutModule& rModule) { fReadoutModules.push_back(rModule); }
 
    void PrintMetadata() { Print(); }
 
    Int_t GetNumberOfChannels();
 
    TRestDetectorReadoutModule* GetModuleByID(Int_t modID);
 
    Int_t isZInsideDriftVolume(Double_t z);
 
    Int_t isZInsideDriftVolume(const TVector3& position);
 
    Bool_t isDaqIDInside(Int_t daqId);
 
    Int_t GetModuleIDFromPosition(const TVector3& position);
 
    Int_t GetModuleIDFromPosition(Double_t x, Double_t y, Double_t z);
 
    TVector2 GetPositionInPlane(const TVector3& point) const;
    Double_t GetDistanceToPlane(const TVector3& point) const;
 
    TVector3 GetPositionInWorld(const TVector2& point, Double_t height = 0) const;
 
    void Draw();
 
    void Print(Int_t DetailLevel = 0);
 
    Int_t FindChannel(Int_t module, const TVector2& position);
 
    Double_t GetX(Int_t modID, Int_t chID);
    Double_t GetY(Int_t modID, Int_t chID);
 
    TH2Poly* GetReadoutHistogram();
    void GetBoundaries(double& xmin, double& xmax, double& ymin, double& ymax);
 
    // Constructor
    TRestDetectorReadoutPlane();
    // Destructor
    virtual ~TRestDetectorReadoutPlane();
 
    ClassDef(TRestDetectorReadoutPlane, 5);
};
#endif