CbmMuchSectorRadial.cxx

Go to the documentation of this file.

/* Copyright (C) 2012-2020 Petersburg Nuclear Physics Institute named by B.P.Konstantinov of National Research Centre "Kurchatov Institute", Gatchina
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Evgeny Kryshen [committer], Florian Uhlig , Apar Agarwal */
 
#include "CbmMuchSectorRadial.h"
 
#include "CbmMuchAddress.h"    // for CbmMuchAddress, kMuchChannel
#include "CbmMuchPadRadial.h"  // for CbmMuchPadRadial
 
#include <Logger.h>  // for Logger, LOG
 
#include <TMath.h>  // for Ceil, Floor, Pi, TwoPi
 
#include <vector>  // for vector
 
 
// -----   Default constructor   -------------------------------------------
CbmMuchSectorRadial::CbmMuchSectorRadial() : CbmMuchSector(), fR1(0.), fR2(0.), fPhi1(0.), fPhi2(0.), fPadDphi(0.) {}
// -------------------------------------------------------------------------
 
 
// -----  Standard constructor  --------------------------------------------
CbmMuchSectorRadial::CbmMuchSectorRadial(Int_t detId, Int_t secId, Double_t r1, Double_t r2, Double_t phi1,
                                         Double_t phi2)
  : CbmMuchSector(detId, secId, 0)
  , fR1(r1)
  , fR2(r2)
  , fPhi1(0.)
  , fPhi2(0.)
  , fPadDphi((r2 - r1) / r1)
{
  if (phi2 < phi1) LOG(error) << "Error in description of module paramaters (phi1 > phi2)";
  fNChannels = TMath::Ceil(((phi2 - phi1) / fPadDphi) - DBL_EPSILON);
  fPhi1      = phi1;
  fPhi2      = phi2;
}
// -----------------o--------------------------------------------------------
 
 
// -----   Public method GetChannel   --------------------------------------
Int_t CbmMuchSectorRadial::GetPadIndexByPhi(Double_t phi)
{
  if (fPhi2 > TMath::Pi() && phi < fPhi1) phi += TMath::TwoPi();
  if (phi < fPhi1 || phi > fPhi2) return -1;
  return TMath::Floor((fPhi2 - phi) / fPadDphi);
}
// -------------------------------------------------------------------------
 
 
// -----   Public method GetChannel   --------------------------------------
CbmMuchPadRadial* CbmMuchSectorRadial::GetPadByPhi(Double_t phi)
{
  Int_t i = GetPadIndexByPhi(phi);
  //  printf("i=%i\n",i);
  if (i == -1) return nullptr;
 
  return (CbmMuchPadRadial*) GetPadByChannelIndex(i);
}
// -------------------------------------------------------------------------
 
 
// -------------------------------------------------------------------------
void CbmMuchSectorRadial::AddPads()
{
  for (Int_t iChannel = 0; iChannel < fNChannels; iChannel++) {
    UInt_t address = CbmMuchAddress::SetElementId(fAddress, kMuchChannel, iChannel);
    Double_t r1, r2, phi1, phi2;
    r1              = fR1;
    r2              = fR2;
    phi1            = fPhi2 - iChannel * fPadDphi;
    phi2            = phi1 - fPadDphi;
    CbmMuchPad* pad = new CbmMuchPadRadial(address, r1, r2, phi1, phi2);
    fPads.push_back(pad);
    LOG(debug) << "iChannel=" << iChannel << " fPads.size()=" << fPads.size() << " fNChannels=" << fNChannels;
  }
}
// -------------------------------------------------------------------------
 
// -------------------------------------------------------------------------
void CbmMuchSectorRadial::DrawPads()
{
  for (Int_t iChannel = 0; iChannel < fNChannels; iChannel++) {
    CbmMuchPadRadial* pad = (CbmMuchPadRadial*) fPads[iChannel];
    pad->DrawPad();
    //     pad->Draw();
  }
}
// -------------------------------------------------------------------------
 
 
ClassImp(CbmMuchSectorRadial)