ARICH histogram with HAPD plane 3 options for bin segmentation are available type 0 - one bin per HAPD, type 1 - one bin per channel, type 2 - one bin per APD chip HAPD numbering is 1-420 (software numbering), channel is ASIC channel, APD is ASIC channel / 36. More...

#include <ARICHChannelHist.h>

Inheritance diagram for ARICHChannelHist:

[legend]

Collaboration diagram for ARICHChannelHist:

[legend]

Public Member Functions
	ARICHChannelHist ()
	Default constructor.

	ARICHChannelHist (const char name, const char title, int type=0, const std::vector< unsigned > &moduleIDs=std::vector< unsigned >())
	Constructor with name, title, and type (0 for channel bins (144/HAPD), 1 for HAPD bins (1/HAPD), 2 for APD bins (4/HAPD)) More...

void	fillBin (unsigned hapdID, unsigned chID, double weight=1.)
	Add entry to bin corresponding to hapd hapdID and channel chID. More...

void	setBinContent (unsigned hapdID, unsigned chID, double value)
	Set content of bin corresponding to hapd hapdID and channel chID. More...

void	fillBin (unsigned hapdID, double weight=1.)
	Add entry to bin corresponding to hapd hapdID. More...

void	setBinContent (unsigned hapdID, double value)
	Set content of bin corresponding to hapd hapdID. More...

void	Draw (Option_t *option="") override
	Draw the histogram. More...

void	fillFromTH1 (TH1 *hist)
	Fill the channelHist from the histogram Type 0 channelHist has to be filled with 420144bin TH1 (each channel) Type 1 channelHist can be filled with 420144 (channels) or 4204 (chips), or 420 (HAPDs) bin TH1 Type 2 channelHist can be filled with 420144 or 420*4 bin TH1 In case histogram has more bins than channelHist bin content is summed up Exampl: histogram with hits/channel can be used to make ChannelHist with hits/hapd. More...

void	setPoly (TH2Poly *poly)
	Fill pure TH2Poly from ARICHChannelHist, makes bins and fills content.

Protected Member Functions
	ClassDefOverride (ARICHChannelHist, 2)
	ClassDef.

Protected Attributes
TLine	lines [6]

TText	labels [6]

int	m_type = 0
	histogram type

std::vector< unsigned >	m_hapd2binMap
	map of bins

Detailed Description

ARICH histogram with HAPD plane 3 options for bin segmentation are available type 0 - one bin per HAPD, type 1 - one bin per channel, type 2 - one bin per APD chip HAPD numbering is 1-420 (software numbering), channel is ASIC channel, APD is ASIC channel / 36.

Definition at line 27 of file ARICHChannelHist.h.

Constructor & Destructor Documentation

◆ ARICHChannelHist()

ARICHChannelHist	(	const char *	name,
		const char *	title,
		int	type = `0`,
		const std::vector< unsigned > &	moduleIDs = `std::vector<unsigned>()`
	)

Constructor with name, title, and type (0 for channel bins (144/HAPD), 1 for HAPD bins (1/HAPD), 2 for APD bins (4/HAPD))

Parameters

name	name
title	title
type	type

Definition at line 21 of file ARICHChannelHist.cc.

                                                                          : TH2Poly()
 {
  
   m_type = type;
   SetName(name);
   SetTitle(title);
   m_hapd2binMap.assign(420, 0);
  
   // positions of HAPDs and channel mapping (avoid using DB classes...)
   const double rs[7] = {57.35, 65.81, 74.37, 82.868, 91.305, 99.794, 108.185};
   const unsigned nhapds[7] = {42, 48, 54, 60, 66, 72, 78};
   const unsigned chmap[144] = {88, 86, 96, 87, 75, 72, 97, 108, 73, 74, 98, 109, 84, 85, 120, 110, 76, 77, 132, 121, 136, 124, 99, 133, 125, 113, 122, 111, 101, 137, 134, 123, 89, 112, 100, 135, 52, 28, 3, 40, 41, 5, 15, 2, 17, 29, 27, 14, 4, 16, 1, 26, 53, 65, 0, 13, 48, 49, 39, 12, 61, 62, 25, 38, 63, 60, 24, 37, 64, 50, 51, 36, 91, 115, 140, 103, 102, 138, 128, 141, 126, 114, 116, 129, 139, 127, 142, 117, 90, 78, 143, 130, 95, 94, 104, 131, 82, 81, 118, 105, 80, 83, 119, 106, 79, 93, 92, 107, 55, 57, 47, 56, 68, 71, 46, 35, 70, 69, 45, 34, 59, 58, 23, 33, 67, 66, 11, 22, 7, 19, 44, 10, 18, 30, 21, 32, 42, 6, 9, 20, 54, 31, 43, 8};
   const double chns[12] = { -2.88, -2.37, -1.86, -1.35, -0.84, -0.33, 0.33, 0.84, 1.35, 1.86, 2.37, 2.88};
  
   float size = 0.5 / 2. - 0.01;
   if (m_type == 1) size = 7.0 / 2. - 0.5;
   if (m_type == 2) size = 3.3 / 2.;
  
   float X[5], Y[5], globX[5], globY[5];
   X[0] = -size;  Y[0] = -size;
   X[1] =  size;  Y[1] = -size;
   X[2] =  size;  Y[2] =  size;
   X[3] = -size;  Y[3] =  size;
   X[4] = -size;  Y[4] = -size;
  
   int nhapd = 1;
   unsigned iring = 0;
   unsigned ihapd = 0;
  
   std::vector<unsigned> ids;
   if (moduleIDs.size() > 0) ids = moduleIDs;
   else {
     for (int hapdID = 1; hapdID < 421; hapdID++) {
       ids.push_back(hapdID);
     }
   }
  
   // HAPD bins
   if (m_type == 1) {
     for (int hapdID = 1; hapdID < 421; hapdID++) {
       //for (unsigned hapdID : ids) {
       //m_hapd2binMap[hapdID - 1] = nhapd;
       //nhapd++;
       float r = rs[iring];
       float dphi = 2.*M_PI / nhapds[iring];
       float fi = dphi / 2. + ihapd * dphi;
       TVector2 centerPos(r * cos(fi), r * sin(fi));
       for (int i = 0; i < 5; i++) {
         float rotX = X[i] * cos(fi) - Y[i] * sin(fi);
         float rotY = X[i] * sin(fi) + Y[i] * cos(fi);
         globX[i] = rotX + centerPos.X();
         globY[i] = rotY + centerPos.Y();
       }
       if (std::find(ids.begin(), ids.end(), hapdID) != ids.end()) {
         m_hapd2binMap[hapdID - 1] = nhapd;
         nhapd++;
         TGraph* mybox = new TGraph(5, globX, globY);
         mybox->SetName((to_string(hapdID)).c_str());
         AddBin(mybox);
       }
       ihapd++;
       if (ihapd == nhapds[iring]) { iring++; ihapd = 0;}
     }
  
   } else if (m_type == 0) {
     for (int hapdID = 1; hapdID < 421; hapdID++) {
       //for (unsigned hapdID : ids) {
       //  m_hapd2binMap[hapdID - 1] = nhapd;
       //  nhapd++;
       float dphi = 2.*M_PI / nhapds[iring];
       float fi = dphi / 2. + ihapd * dphi;
       float r = rs[iring];
       for (int chID = 0; chID < 144; chID++) {
  
         unsigned chX = chmap[chID] % 12;
         unsigned chY = chmap[chID] / 12;
         TVector2 hapdPos(r * cos(fi), r * sin(fi));
         TVector2 locPos(chns[chX], chns[chY]);
         TVector2 centerPos = hapdPos + locPos.Rotate(fi);
  
         for (int i = 0; i < 5; i++) {
           float rotX = X[i] * cos(fi) - Y[i] * sin(fi);
           float rotY = X[i] * sin(fi) + Y[i] * cos(fi);
           globX[i] = rotX + centerPos.X();
           globY[i] = rotY + centerPos.Y();
         }
         if (std::find(ids.begin(), ids.end(), hapdID) != ids.end()) {
           m_hapd2binMap[hapdID - 1] = nhapd;
           if (chID == 143) nhapd++;
           TGraph* mybox = new TGraph(5, globX, globY);
           mybox->SetName((to_string(hapdID)).c_str());
           AddBin(mybox);
         }
       }
       ihapd++;
       if (ihapd == nhapds[iring]) { iring++; ihapd = 0;}
     }
   } else if (m_type == 2) {
     size += 0.2;
     for (int hapdID = 1; hapdID < 421; hapdID++) {
       float dphi = 2.*M_PI / nhapds[iring];
       float fi = dphi / 2. + ihapd * dphi;
       float r = rs[iring];
       TVector2 hapdPos(r * cos(fi), r * sin(fi));
       for (int chipID = 0; chipID < 4; chipID++) {
         TVector2 locPos(-size + (chipID / 2)*size * 2, size - (chipID % 2)*size * 2);
         TVector2 centerPos = hapdPos + locPos.Rotate(fi);
  
         for (int i = 0; i < 5; i++) {
           float rotX = X[i] * cos(fi) - Y[i] * sin(fi);
           float rotY = X[i] * sin(fi) + Y[i] * cos(fi);
           globX[i] = rotX + centerPos.X();
           globY[i] = rotY + centerPos.Y();
         }
  
         if (std::find(ids.begin(), ids.end(), hapdID) != ids.end()) {
           m_hapd2binMap[hapdID - 1] = nhapd;
           if (chipID == 3) nhapd++;
           TGraph* mybox = new TGraph(5, globX, globY);
           mybox->SetName((to_string(hapdID)).c_str());
           AddBin(mybox);
         }
       }
       ihapd++;
       if (ihapd == nhapds[iring]) { iring++; ihapd = 0;}
     }
   } else  std::cout << "Invalid histogram type! use 0 for channel bins or 1 for HAPD bins" << std::endl;
   SetOption("colz");
   SetStats(0);
   GetXaxis()->SetLimits(-115., 115.);
   GetYaxis()->SetLimits(-115., 115.);
 }

Member Function Documentation

◆ Draw()

void Draw ( Option_t * option = "" )

override

Draw the histogram.

Parameters

option root draw option

Definition at line 155 of file ARICHChannelHist.cc.

◆ fillBin() [1/2]

void fillBin	(	unsigned	hapdID,
		double	weight = `1.`
	)

Add entry to bin corresponding to hapd hapdID.

Parameters

hapdID ID number of HAPD module

Definition at line 198 of file ARICHChannelHist.cc.

◆ fillBin() [2/2]

void fillBin	(	unsigned	hapdID,
		unsigned	chID,
		double	weight = `1.`
	)

Add entry to bin corresponding to hapd hapdID and channel chID.

Parameters

hapdID	ID number of HAPD module
chID	channel ID (asic channel) or APD ID

Definition at line 177 of file ARICHChannelHist.cc.

◆ fillFromTH1()

void fillFromTH1 ( TH1 * hist )

Fill the channelHist from the histogram Type 0 channelHist has to be filled with 420*144bin TH1 (each channel) Type 1 channelHist can be filled with 420*144 (channels) or 420*4 (chips), or 420 (HAPDs) bin TH1 Type 2 channelHist can be filled with 420*144 or 420*4 bin TH1 In case histogram has more bins than channelHist bin content is summed up Exampl: histogram with hits/channel can be used to make ChannelHist with hits/hapd.

Parameters

hist historam

Definition at line 203 of file ARICHChannelHist.cc.

◆ setBinContent() [1/2]

void setBinContent	(	unsigned	hapdID,
		double	value
	)

Set content of bin corresponding to hapd hapdID.

Parameters

hapdID	ID number of HAPD module
value	bin content

Definition at line 193 of file ARICHChannelHist.cc.

◆ setBinContent() [2/2]

void setBinContent	(	unsigned	hapdID,
		unsigned	chID,
		double	value
	)

Set content of bin corresponding to hapd hapdID and channel chID.

Parameters

hapdID	ID number of HAPD module
chID	channel ID (asic channel) or APD ID
value	bin content

Definition at line 185 of file ARICHChannelHist.cc.

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ARICHChannelHist()

Member Function Documentation

◆ Draw()

◆ fillBin() [1/2]

◆ fillBin() [2/2]

◆ fillFromTH1()

◆ setBinContent() [1/2]

◆ setBinContent() [2/2]