Base class for geometry parameters. More...

#include <ARICHAerogelHist.h>

Inheritance diagram for ARICHAerogelHist:

Public Member Functions
	ARICHAerogelHist ()
	Default constructor.

	~ARICHAerogelHist ()
	Default destructor.

	ARICHAerogelHist (const char name, const char title)
	Constructor with name, title.

Int_t	GetBinIDFromRingColumn (Int_t ring, Int_t column)
	Function which return histogram bin id from ring and column id's.

void	DrawHisto (TString opt, TString outDirName)
	Function to draw the histogram.

Protected Member Functions
void	SetInitialParametersByDefault ()
	Function which set initial values of input parameters.

void	SetUpVerticesMap ()
	Function for calculation vertices for one aerogel tile.

void	dumpVerticesMap ()
	Function to print vertices for one aerogel tile.

void	makeRotation (double xold, double yold, double &xnew, double &ynew, double phi)
	Function to rotate 2D point (x and y) around z axis by angle phi.

	ClassDef (ARICHAerogelHist, 2)
	ClassDef.

Protected Attributes
std::vector< Int_t >	m_nTiles
	Number of tiles per ring.

std::vector< double >	m_tileRmin
	Minimum radius of aerogel ring.

std::vector< double >	m_tileRmax
	Maximum radius of aerogel ring.

std::vector< double >	m_tileRcenter
	Center radius of aerogel ring.

std::vector< double >	m_tileDeltaPhiCenter
	Angle opening (phi) of the aerogel tile measured between two rays (0.0,0.0 : and centre of the ring from left/right).

std::vector< double >	m_aerogelAriGapDeltaPhiCenter
	Angle opening (phi) of the air gap between aerogel tiles.

Int_t	m_verboseLevel
	Verbose level.

Int_t	m_nCircularPoints
	Number of circular points.

double	m_aerogelTileGap
	Distance between aerogel tiles.

std::map< Int_t, std::vector< ROOT::Math::XYVector > >	m_verticesMap
	Aerogel vertices map.

TString	m_histName
	Histogram name.

TString	m_histTitle
	Histogram title.

Detailed Description

Base class for geometry parameters.

Definition at line 28 of file ARICHAerogelHist.h.

Constructor & Destructor Documentation

◆ ARICHAerogelHist() [1/2]

ARICHAerogelHist ( )

inline

Default constructor.

Definition at line 35 of file ARICHAerogelHist.h.

35{};

◆ ~ARICHAerogelHist()

~ARICHAerogelHist ( )

inline

Default destructor.

Definition at line 40 of file ARICHAerogelHist.h.

40{};

◆ ARICHAerogelHist() [2/2]

ARICHAerogelHist	(	const char *	name,
		const char *	title
	)

Constructor with name, title.

Parameters

name	name.
title	title.

Definition at line 29 of file ARICHAerogelHist.cc.

                                                                      : TH2Poly(), m_histName(name), m_histTitle(title)
{
  SetName(m_histName.Data());
  SetTitle(m_histTitle.Data());
  SetInitialParametersByDefault();
  SetUpVerticesMap();
 
  if (m_verboseLevel > 0) {
    std::cout << setw(10) << "ring" << setw(10) << "nTiles" << setw(10) << "tileRmin" << setw(10) << "tileRmax" << std::endl;
    for (unsigned int i = 0; i < m_nTiles.size(); i++)
      std::cout << setw(10) << i + 1 << setw(10) << m_nTiles[i] << setw(10) << m_tileRmin[i] << setw(10) << m_tileRmax[i] << std::endl;
    dumpVerticesMap();
  }
 
  unsigned int n;
  double* x;
  double* y;
 
  double xold;
  double yold;
  double xnew;
  double ynew;
  double phi;
 
  //Add "poly bins" to the histogram
  //Loop over rings with different radiuses
  unsigned int iR = 0;
  for (auto& m : m_verticesMap) {
    //Loop aerogel tiles (bins) within one ring.
    double phi0 = m_tileDeltaPhiCenter[iR] / 2.0;
    double deltaPhi = m_tileDeltaPhiCenter[iR] + m_aerogelAriGapDeltaPhiCenter[iR];
    for (Int_t j = 0; j < m_nTiles[iR]; j++) {
      phi = phi0 + deltaPhi * j;
      n = m.second.size();
      x = new double [n];
      y = new double [n];
      //Loop over polygonal points which defines bins.
      for (unsigned int i = 0; i < n; i++) {
        xold = m.second[i].X();
        yold = m.second[i].Y();
        makeRotation(xold, yold, xnew, ynew, phi);
        x[i] = xnew;
        y[i] = ynew;
      }
      AddBin(n, x, y);
      delete []x;
      delete []y;
    }
    iR++;
  }// for (auto& m: m_verticesMap) {
 
}

Member Function Documentation

◆ DrawHisto()

void DrawHisto	(	TString	opt = `"ZCOLOT text same"`,
		TString	outDirName = `"./"`
	)

Function to draw the histogram.

Parameters

opt	draw option string default value : "ZCOLOT text same".
outDirName	name of epe and pdf to save the plots.

Definition at line 104 of file ARICHAerogelHist.cc.

{
 
  TCanvas* c1 = new TCanvas("c1", "c1", 1000, 1000);
  c1->SetTitle(m_histTitle.Data());
  c1->SetName(m_histName.Data());
  c1->SetRightMargin(0.17);
  c1->SetLeftMargin(0.12);
  c1->SetTopMargin(0.15);
  c1->SetBottomMargin(0.15);
 
  TH2F* frame = new TH2F("h2", "h2", 40, -1200, 1200, 40, -1200, 1200);
  //TH2F *frame = new TH2F("h2","h2",40, 400, 650, 40, -50, 200);
  frame->GetXaxis()->SetTitle("x (mm)");
  frame->GetYaxis()->SetTitle("y (mm)");
  frame->GetXaxis()->CenterTitle();
  frame->GetYaxis()->CenterTitle();
  frame->GetYaxis()->SetTitleOffset(1.5);
  frame->SetStats(kFALSE);
  frame->SetTitle("");
  frame->Draw();
 
  Draw(opt.Data());
  c1->Modified();
  c1->Update();
  if (outDirName.Length() > 0) {
    std::cout << "outDirName.Length() " << outDirName.Length() << std::endl;
    TString outnamePDF = outDirName; outnamePDF += m_histName; outnamePDF += ".pdf";
    TString outnameEPS = outDirName; outnameEPS += m_histName; outnameEPS += ".eps";
    c1->SaveAs(outnamePDF.Data());
    c1->SaveAs(outnameEPS.Data());
  }
 
}

◆ dumpVerticesMap()

void dumpVerticesMap ( )

protected

Function to print vertices for one aerogel tile.

Definition at line 139 of file ARICHAerogelHist.cc.

{
 
  unsigned int i;
  for (auto& m : m_verticesMap) {
    std::cout << " --> Aerogel ring : " << m.first << '\n';
    for (i = 0; i < m.second.size(); i++) {
      std::cout << "                    " << setw(15) << m.second[i].X() << setw(15) << m.second[i].Y() << std::endl;
    }
  }
 
}

◆ GetBinIDFromRingColumn()

Int_t GetBinIDFromRingColumn	(	Int_t	ring,
		Int_t	column
	)

Function which return histogram bin id from ring and column id's.

Parameters

ring	ringID number.
column	columnID number.

Definition at line 82 of file ARICHAerogelHist.cc.

{
 
  Int_t binID = 0;
  unsigned int ringu = ring;
  //cout<<" ewewe "<<m_nTiles[ring-1]<<endl;
  if (ringu > m_nTiles.size() || ringu < 1)
    return -999;
  if (column > m_nTiles[ring - 1] || column < 1)
    return -999;
  for (Int_t i = 1; i < ring; i++)
    binID += m_nTiles[i - 1];
  return binID + column;
 
}

◆ makeRotation()

void makeRotation	(	double	xold,
		double	yold,
		double &	xnew,
		double &	ynew,
		double	phi
	)

protected

Function to rotate 2D point (x and y) around z axis by angle phi.

Parameters

xold	old x coordinate.
yold	old y coordinate.
xnew	new x coordinate.
ynew	new y coordinate.
phi	roration angle.

Definition at line 98 of file ARICHAerogelHist.cc.

{
  xnew = xold * std::cos(phi) - yold * std::sin(phi);
  ynew = xold * std::sin(phi) + yold * std::cos(phi);
}

◆ SetInitialParametersByDefault()

void SetInitialParametersByDefault ( )

protected

Function which set initial values of input parameters.

Definition at line 206 of file ARICHAerogelHist.cc.

{
 
  //Verbose level
  m_verboseLevel = 0;
  //Distance between aerogel tiles
  m_aerogelTileGap = 1;
  //Number of circular points
  m_nCircularPoints = 10;
  //Number of tiles per ring.
  m_nTiles.push_back(22);
  m_nTiles.push_back(28);
  m_nTiles.push_back(34);
  m_nTiles.push_back(40);
  //Minimum radius of aerogel ring
  m_tileRmin.push_back(441.0);
  m_tileRmin.push_back(614.0);
  m_tileRmin.push_back(787.0);
  m_tileRmin.push_back(960.0);
  //Maximum radius of aerogel ring
  m_tileRmax.push_back(613.0);
  m_tileRmax.push_back(786.0);
  m_tileRmax.push_back(959.0);
  m_tileRmax.push_back(1117.0);
  for (unsigned int i = 0; i < m_nTiles.size(); i++) {
    m_tileRcenter.push_back((m_tileRmax[i] + m_tileRmin[i]) / 2.0);
    double r = m_tileRcenter[i];
    double l = 2 * TMath::Pi() * r / m_nTiles[i] - m_aerogelTileGap;
    double phi = l / r;
    m_tileDeltaPhiCenter.push_back(phi);
    m_aerogelAriGapDeltaPhiCenter.push_back(m_aerogelTileGap / r);
  }
 
}

◆ SetUpVerticesMap()

void SetUpVerticesMap ( )

protected

Function for calculation vertices for one aerogel tile.

Definition at line 152 of file ARICHAerogelHist.cc.

{
 
  for (unsigned int i = 0; i < m_nTiles.size(); i++) {
    std::vector<ROOT::Math::XYVector> vecTvec;
    int nTiles = m_nTiles[i];
    double rmin = m_tileRmin[i];
    double rmax = m_tileRmax[i];
    double lmin = 2 * TMath::Pi() * rmin / nTiles - m_aerogelTileGap;
    double lmax = 2 * TMath::Pi() * rmax / nTiles - m_aerogelTileGap;
    double phimin = lmin / rmin;
    double phimax = lmax / rmax;
    double x1 = rmin * TMath::Cos(phimin / 2.0);
    double y1 = rmin * TMath::Sin(phimin / 2.0);
    ROOT::Math::XYVector v1(x1, y1);
    vecTvec.push_back(v1);
    double x2 = rmax * TMath::Cos(phimax / 2.0);
    double y2 = rmax * TMath::Sin(phimax / 2.0);
    ROOT::Math::XYVector v2(x2, y2);
    vecTvec.push_back(v2);
    //Add circular points from outer radious (clockwise added)
    if (m_nCircularPoints > 0) {
      double dPhi = phimax / (m_nCircularPoints + 1);
      for (int j = 0; j < m_nCircularPoints; j++) {
        const double phi = phimax / 2.0 - dPhi * (j + 1);
        ROOT::Math::XYVector v(rmax * std::cos(phi), rmax * std::sin(phi));
        vecTvec.push_back(v);
      }
    }
 
    double x3 =  x2;
    double y3 = -y2;
    ROOT::Math::XYVector v3(x3, y3);
    vecTvec.push_back(v3);
    double x4 =  x1;
    double y4 = -y1;
    ROOT::Math::XYVector v4(x4, y4);
    vecTvec.push_back(v4);
    //Add circular points inner radious
    if (m_nCircularPoints > 0) {
      double dPhi = phimin / (m_nCircularPoints + 1);
      for (int j = 0; j < m_nCircularPoints; j++) {
        const double phi = -phimax / 2.0 + dPhi * (j + 1);
        ROOT::Math::XYVector v(rmax * std::cos(phi), rmax * std::sin(phi));
        vecTvec.push_back(v);
      }
    }
 
    vecTvec.push_back(v1);
    m_verticesMap[i] = vecTvec;
  }
 
}

Member Data Documentation

◆ m_aerogelAriGapDeltaPhiCenter

std::vector<double> m_aerogelAriGapDeltaPhiCenter

protected

Angle opening (phi) of the air gap between aerogel tiles.

Measured between ray (0.0,0.0 : and centre of the ring from left/right).

Definition at line 97 of file ARICHAerogelHist.h.

◆ m_aerogelTileGap

double m_aerogelTileGap

protected

Distance between aerogel tiles.

Definition at line 100 of file ARICHAerogelHist.h.

◆ m_histName

TString m_histName

protected

Histogram name.

Definition at line 105 of file ARICHAerogelHist.h.

◆ m_histTitle

TString m_histTitle

protected

Histogram title.

Definition at line 106 of file ARICHAerogelHist.h.

◆ m_nCircularPoints

Int_t m_nCircularPoints

protected

Number of circular points.

Definition at line 99 of file ARICHAerogelHist.h.

◆ m_nTiles

std::vector<Int_t> m_nTiles

protected

Number of tiles per ring.

Definition at line 90 of file ARICHAerogelHist.h.

◆ m_tileDeltaPhiCenter

std::vector<double> m_tileDeltaPhiCenter

protected

Angle opening (phi) of the aerogel tile measured between two rays (0.0,0.0 : and centre of the ring from left/right).

Definition at line 95 of file ARICHAerogelHist.h.

◆ m_tileRcenter

std::vector<double> m_tileRcenter

protected

Center radius of aerogel ring.

Definition at line 93 of file ARICHAerogelHist.h.

◆ m_tileRmax

std::vector<double> m_tileRmax

protected

Maximum radius of aerogel ring.

Definition at line 92 of file ARICHAerogelHist.h.

◆ m_tileRmin

std::vector<double> m_tileRmin

protected

Minimum radius of aerogel ring.

Definition at line 91 of file ARICHAerogelHist.h.

◆ m_verboseLevel

Int_t m_verboseLevel

protected

Verbose level.

Definition at line 98 of file ARICHAerogelHist.h.

◆ m_verticesMap

std::map<Int_t, std::vector<ROOT::Math::XYVector> > m_verticesMap

protected

Aerogel vertices map.

Definition at line 103 of file ARICHAerogelHist.h.

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

arich/utility/include/ARICHAerogelHist.h
arich/utility/src/ARICHAerogelHist.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ARICHAerogelHist() [1/2]

◆ ~ARICHAerogelHist()

◆ ARICHAerogelHist() [2/2]

Member Function Documentation

◆ DrawHisto()

◆ dumpVerticesMap()

◆ GetBinIDFromRingColumn()

◆ makeRotation()

◆ SetInitialParametersByDefault()

◆ SetUpVerticesMap()

Member Data Documentation

◆ m_aerogelAriGapDeltaPhiCenter

◆ m_aerogelTileGap

◆ m_histName

◆ m_histTitle

◆ m_nCircularPoints

◆ m_nTiles

◆ m_tileDeltaPhiCenter

◆ m_tileRcenter

◆ m_tileRmax

◆ m_tileRmin

◆ m_verboseLevel

◆ m_verticesMap