ARICHMergerHist Class Reference

ARICH merger histogram (72 blocks). More...

#include <ARICHMergerHist.h>

Inheritance diagram for ARICHMergerHist:

Public Member Functions
	ARICHMergerHist ()=default
	Default constructor.
	ARICHMergerHist (const char *name, const char *title, double globalRotationDeg=0.0, bool rightToLeftNumbering=true)
	Constructor.
void	fillFromTH1 (TH1 *hist)
	Fill from TH1.
void	SetBinContentFromSN (unsigned mergerID, double value)
	Set content using the real ARICH hardware MergerID.
void	Draw (Option_t *option="") override
	Draw with sector lines and labels.

Static Protected Member Functions
static void	pol2xy (double r, double angDeg, double &x, double &y)
	Convert polar coordinates to cartesian.

Protected Attributes
std::vector< int >	m_merger2bin
	Merger number (1..72) → TH2Poly bin index.
double	m_globalRotationDeg = 0.0
	Global rotation applied to the geometry (degrees).
bool	m_rightToLeft = true
	If true, numbering in each ring is right-to-left.
TLine	m_lines [6]
	Sector boundary lines drawn on top of the histogram.

Detailed Description

ARICH merger histogram (72 blocks).

Two numbering systems exist:

1) Internal position (1..72)

• This is the drawing / geometry order.
• Used by fillFromTH1().
• Bins correspond to 1_1, 2_1, 12_1,...,6_12 merger positions

2) Real ARICH hardware MergerID

• Used by SetBinContentFromSN(mergerID, value).

Definition at line 34 of file ARICHMergerHist.h.

Constructor & Destructor Documentation

ARICHMergerHist()

ARICHMergerHist	(	const char *	name,
		const char *	title,
		double	globalRotationDeg = 0.0,
		bool	rightToLeftNumbering = true )

Constructor.

Parameters

name	Histogram name
title	Histogram title
globalRotationDeg	Global rotation applied to the whole layout (degrees)
rightToLeftNumbering	If true, bins inside each ring are numbered right-to-left

Definition at line 69 of file ARICHMergerHist.cc.

  : TH2Poly()
{
  m_globalRotationDeg = globalRotationDeg;
  m_rightToLeft = rightToLeftNumbering;
  m_merger2bin.assign(72, 0);
 
  SetName(name);
  SetTitle(title);
 
  int polyBin = 1;
  int mergerNumber = 1;
 
  auto add_block = [&](double r, double angDeg) {
 
    double cx, cy;
    pol2xy(r, angDeg, cx, cy);
 
    const double rot = (angDeg + 90.0) * M_PI / 180.0;
    const double c = std::cos(rot);
    const double s = std::sin(rot);
 
    const double hx = BLOCK_W / 2.0;
    const double hy = BLOCK_H / 2.0;
 
    const double lx[5] = {-hx, +hx, +hx, -hx, -hx};
    const double ly[5] = {-hy, -hy, +hy, +hy, -hy};
 
    double gx[5], gy[5];
    for (int i = 0; i < 5; ++i) {
      const double rx = lx[i] * c - ly[i] * s;
      const double ry = lx[i] * s + ly[i] * c;
      gx[i] = rx + cx;
      gy[i] = ry + cy;
    }
 
    auto* g = new TGraph(5, gx, gy);
    AddBin(g);
 
    if (mergerNumber >= 1 && mergerNumber <= 72)
      m_merger2bin[mergerNumber - 1] = polyBin;
 
    ++polyBin;
    ++mergerNumber;
  };
 
  auto add_ring = [&](double r, const double * offsets, int nOffsets, double centerDeg) {
    if (m_rightToLeft) {
      for (int i = nOffsets - 1; i >= 0; --i)
        add_block(r, centerDeg + offsets[i]);
    } else {
      for (int i = 0; i < nOffsets; ++i)
        add_block(r, centerDeg + offsets[i]);
    }
  };
 
  for (int X = 1; X <= 6; ++X) {
    const double cdeg = SECTOR_CENTER_DEG[X] + m_globalRotationDeg;
 
    add_ring(R1, ROW1_OFFSETS, 3, cdeg);
    add_ring(R2, ROW2_OFFSETS, 4, cdeg);
    add_ring(R3, ROW3_OFFSETS, 4, cdeg);
    add_block(R4, cdeg + ROW4_OFFSETS[0]);
  }
 
  if (mergerNumber != 73)
    B2ERROR("ARICHMergerHist internal error: expected 72 bins, got "
            << mergerNumber - 1);
 
  SetOption("colz");
  SetStats(0);
  GetXaxis()->SetLimits(-140., 140.);
  GetYaxis()->SetLimits(-140., 140.);
}

Member Function Documentation

Draw()

void Draw ( Option_t * option = "" )

override

Draw with sector lines and labels.

Parameters

option

ROOT draw option string

Definition at line 179 of file ARICHMergerHist.cc.

{
  TH2Poly::Draw(option);
 
  // -------------------------------------------------
  // Sector lines
  // -------------------------------------------------
  const double rInner = 50.0;
  const double rOuter = 130.0;
 
  for (int k = 0; k < 6; ++k) {
    const double aDeg = k * 60.0 + m_globalRotationDeg;
 
    double x1, y1, x2, y2;
    pol2xy(rInner, aDeg, x1, y1);
    pol2xy(rOuter, aDeg, x2, y2);
 
    m_lines[k] = TLine(x1, y1, x2, y2);
    m_lines[k].SetLineWidth(2);
    m_lines[k].SetLineColor(kGray + 2);
    m_lines[k].Draw();
  }
 
  // -------------------------------------------------
  // Sector labels S-1 ... S-6
  // -------------------------------------------------
  const double rLabel = 45.0;
 
  for (int s = 0; s < 6; ++s) {
    const double centerDeg = SECTOR_CENTER_DEG[s + 1] + m_globalRotationDeg;
 
    double x, y;
    pol2xy(rLabel, centerDeg, x, y);
 
    auto* txt = new TLatex(x, y, Form("S-%d", s + 1));
    txt->SetBit(kCanDelete);
    txt->SetTextAlign(22);
    txt->SetTextFont(62);
    txt->SetTextSize(0.030);
    txt->Draw();
  }
}

fillFromTH1()

void fillFromTH1

(

TH1 *

hist

)

Fill from TH1.

Histogram must have exactly 72 bins.

Parameters

hist	Input TH1 with 72 bins

Definition at line 150 of file ARICHMergerHist.cc.

{
  if (!hist) return;
 
  if (hist->GetNbinsX() != 72) {
    B2ERROR("ARICHMergerHist::fillFromTH1: histogram must have 72 bins (got "
            << hist->GetNbinsX() << ")");
    return;
  }
 
  for (int i = 1; i <= 72; ++i) {
    const int b = m_merger2bin[i - 1];
    SetBinContent(b, hist->GetBinContent(i));
  }
}

pol2xy()

void pol2xy ( double r, double angDeg, double & x, double & y )

staticprotected

Convert polar coordinates to cartesian.

Parameters

r	Radius
angDeg	Angle in degrees
x	Output x coordinate
y	Output y coordinate

Definition at line 62 of file ARICHMergerHist.cc.

{
  const double th = angDeg * M_PI / 180.0;
  x = r * std::cos(th);
  y = r * std::sin(th);
}

SetBinContentFromSN()

void SetBinContentFromSN	(	unsigned	mergerID,
		double	value )

Set content using the real ARICH hardware MergerID.

Definition at line 169 of file ARICHMergerHist.cc.

{
  if (mergerID < 1 || mergerID > 100) return;
 
  int internal = MERGER_TO_INTERNAL.at(mergerID);
  int polybin  = m_merger2bin[internal - 1];
 
  SetBinContent(polybin, value);
}

Member Data Documentation

m_globalRotationDeg

double m_globalRotationDeg = 0.0

protected

Global rotation applied to the geometry (degrees).

Definition at line 63 of file ARICHMergerHist.h.

m_lines

TLine m_lines[6]

protected

Sector boundary lines drawn on top of the histogram.

Definition at line 65 of file ARICHMergerHist.h.

m_merger2bin

std::vector<int> m_merger2bin

protected

Merger number (1..72) → TH2Poly bin index.

Definition at line 62 of file ARICHMergerHist.h.

m_rightToLeft

bool m_rightToLeft = true

protected

If true, numbering in each ring is right-to-left.

Definition at line 64 of file ARICHMergerHist.h.

---

The documentation for this class was generated from the following files:

• arich/utility/include/ARICHMergerHist.h
• arich/utility/src/ARICHMergerHist.cc