class to summarize SVD quantities per sensor and side More...

#include <SVDSummaryPlots.h>

Inheritance diagram for SVDSummaryPlots:

Public Types
enum	E_side { VIndex = 0 , UIndex = 1 }
	This enumeration assure the same semantic of the isU methods defined by Peter Kv. More...

Public Member Functions
	SVDSummaryPlots ()
	Default constructor.

	SVDSummaryPlots (TString name, TString title)
	this is the default constructor

	~SVDSummaryPlots ()
	clean everything in the destructor

TH2F *	getHistogram (int view)
	get a reference to the histogram for

float	getValue (const VxdID &vxdID, int view)
	get the value contained in the corresponding bin, given VxdID and view

float	getValue (int layer, int ladder, int sensor, int view)
	get the value contained in the corresponding bin, given

void	fill (int layer, int ladder, int sensor, int view, float value)
	fill the histogram for

void	fill (const VxdID &vxdID, int view, float value)
	fill the histogram for

void	setRunID (const TString &runID)
	set run ids in title

void	reset ()
	Reset histograms.

void	setStats (bool stats=true)
	set histograms stat

void	fill (const VxdID &vxdID, bool isU, float value)
	fill the histogram for

void	customizeString (std::string &base, bool isU)
	replaces layer ladder sensor view and apv with the current numbers

void	setMinimum (Int_t value=0)
	set histogram minimum

void	setMaximum (Int_t value=0)
	set histogram maximum

TH2Poly *	getPoly (int view, int min=-1111, int max=-1111)
	Create the TH2Poly version of the plot.

void	clean ()
	delete pointers

Private Member Functions
Int_t	findBinY (Int_t layer, Int_t sensor)
	find the Y bin given the layer and sensor number

void	customize (TH2F &histogram, int view)
	customize the histogram with the sensor, view

void	generateHistogramBins (TH2Poly *histogram)
	generate histogram bins

Int_t	findBinPoly (Int_t layer, Int_t ladder, Int_t sensor)
	find the poly-bin given the layer, ladder and sensor number

	ClassDef (SVDSummaryPlots, 3)
	needed by root

Private Attributes
TH2F *	m_histos [2]
	vector containing the U and V histograms

TH2Poly *	m_polyHistos [2] = {nullptr}
	vector containing the U and V poly-histograms

TH2F *	m_defaultHistogram = nullptr
	default histograms

TString	m_Title [2]
	Base title.

Detailed Description

class to summarize SVD quantities per sensor and side

Definition at line 26 of file SVDSummaryPlots.h.

Member Enumeration Documentation

◆ E_side

enum E_side

This enumeration assure the same semantic of the isU methods defined by Peter Kv.

Definition at line 61 of file SVDSummaryPlots.h.

61{ VIndex = 0, UIndex = 1 };

Constructor & Destructor Documentation

◆ SVDSummaryPlots() [1/2]

SVDSummaryPlots ( )

inline

Default constructor.

Definition at line 30 of file SVDSummaryPlots.h.

30 :

31 SVDSummaryPlots("", "") {};

Belle2::SVDSummaryPlots::SVDSummaryPlots

SVDSummaryPlots()

Default constructor.

Definition: SVDSummaryPlots.h:30

◆ SVDSummaryPlots() [2/2]

SVDSummaryPlots	(	TString	name,
		TString	title
	)

inline

this is the default constructor

Parameters

name	histogram name
title	histogram title

Definition at line 37 of file SVDSummaryPlots.h.

    {
      m_defaultHistogram = new TH2F(name.Data(), title.Data(),
                                    16, 0.5, 16.5,
                                    19, 0, 19);
 
      for (int view = VIndex ; view < UIndex + 1; view++) {
        std::string hname = name.Data();
        bool isU = (view == UIndex);
        customizeString(hname, isU);
        m_histos[view] = new TH2F(hname.c_str(), title.Data(),
                                  16, 0.5, 16.5,
                                  19, 0, 19);
        customize(*m_histos[view], view);
        m_Title[view] = m_histos[view]->GetTitle();
      }
    }

◆ ~SVDSummaryPlots()

~SVDSummaryPlots ( )

inline

clean everything in the destructor

Definition at line 57 of file SVDSummaryPlots.h.

57{ clean(); };

Belle2::SVDSummaryPlots::clean

void clean()

delete pointers

Definition: SVDSummaryPlots.h:246

Member Function Documentation

◆ clean()

void clean ( )

inline

delete pointers

Definition at line 246 of file SVDSummaryPlots.h.

    {
      delete m_histos[0];
      delete m_histos[1];
      delete m_defaultHistogram;
    }

◆ customize()

void customize	(	TH2F &	histogram,
		int	view
	)

inlineprivate

customize the histogram with the sensor, view

Definition at line 280 of file SVDSummaryPlots.h.

    {
 
      const int nY = 19;
      TString Ylabels[nY] = {"", "L3.x.1", "L3.x.2",
                             " ", "L4.x.1", "L4.x.2", "L4.x.3",
                             "  ", "L5.x.1", "L5.x.2", "L5.x.3", "L5.x.4",
                             "   ", "L6.x.1", "L6.x.2", "L6.x.3", "L6.x.4", "L6.x.5", "    "
                            };
 
      histogram.SetMarkerSize(1.1);
      histogram.GetXaxis()->SetTitle("ladder number");
      histogram.GetXaxis()->SetLabelSize(0.04);
      for (unsigned short i = 0; i < nY; i++)
        histogram.GetYaxis()->SetBinLabel(i + 1, Ylabels[i].Data());
 
      bool isU = view == UIndex;
      std::string title = histogram.GetTitle();
      customizeString(title, isU);
      histogram.SetTitle(title.c_str());
 
      std::string xAxis = histogram.GetXaxis()->GetTitle();
      customizeString(xAxis, isU);
      histogram.SetXTitle(xAxis.c_str());
 
      std::string yAxis = histogram.GetYaxis()->GetTitle();
      customizeString(yAxis, isU);
      histogram.SetYTitle(yAxis.c_str());
 
      std::string zAxis = histogram.GetZaxis()->GetTitle();
      customizeString(zAxis, isU);
      histogram.SetZTitle(zAxis.c_str());
 
    }

◆ customizeString()

void customizeString	(	std::string &	base,
		bool	isU
	)

inline

replaces layer ladder sensor view and apv with the current numbers

Definition at line 170 of file SVDSummaryPlots.h.

    {
      std::string view = isU ? "U" : "V" ;
      base = std::regex_replace(base, std::regex("[@]view"), view);
      std::string side = isU ? "P" : "N" ;
      base = std::regex_replace(base, std::regex("[@]side"), side);
    }

◆ fill() [1/3]

void fill	(	const VxdID &	vxdID,
		bool	isU,
		float	value
	)

inline

fill the histogram for

Parameters

vxdID	side
isU	with
value

Definition at line 158 of file SVDSummaryPlots.h.

    {
      int view = isU ? UIndex : VIndex;
 
      int layer = vxdID.getLayerNumber();
      int ladder = vxdID.getLadderNumber();
      int sensor = vxdID.getSensorNumber();
 
      fill(layer, ladder, sensor, view, value);
    }

◆ fill() [2/3]

void fill	(	const VxdID &	vxdID,
		int	view,
		float	value
	)

inline

fill the histogram for

Parameters

vxdID
view	with
value

Definition at line 117 of file SVDSummaryPlots.h.

    {
      int layer = vxdID.getLayerNumber();
      int ladder = vxdID.getLadderNumber();
      int sensor = vxdID.getSensorNumber();
 
      fill(layer, ladder, sensor, view, value);
    }

◆ fill() [3/3]

void fill	(	int	layer,
		int	ladder,
		int	sensor,
		int	view,
		float	value
	)

inline

fill the histogram for

Parameters

layer
ladder
sensor
view	with
value

Definition at line 109 of file SVDSummaryPlots.h.

    {
      getHistogram(view)->Fill(ladder, findBinY(layer, sensor), value);
    }

◆ findBinPoly()

Int_t findBinPoly	(	Int_t	layer,
		Int_t	ladder,
		Int_t	sensor
	)

inlineprivate

find the poly-bin given the layer, ladder and sensor number

Definition at line 359 of file SVDSummaryPlots.h.

    {
      const int initBinPerLayer[4] = {1, 15, 45, 93};
 
      const double nLadders[4] = {7, 10, 12, 16}; // per layer
      const double nSensors[4] = {2, 3, 4, 5}; // per ladder
 
      int bin = initBinPerLayer[layer - 3];
 
      for (int ld = 1; ld <= nLadders[layer - 3]; ld++) {
        bool binFound = false;
        for (int s = 1; s <= nSensors[layer - 3]; s++) {
          if (ld == ladder && s == sensor) binFound = true;
          if (binFound) break;
          bin += 1;
        }
        if (binFound) break;
      }
      return bin;
    }

◆ findBinY()

Int_t findBinY	(	Int_t	layer,
		Int_t	sensor
	)

inlineprivate

find the Y bin given the layer and sensor number

Definition at line 256 of file SVDSummaryPlots.h.

    {
 
      if (layer == 3)
        return sensor; //2
      if (layer == 4)
        return 2 + 1 + sensor; //6
      if (layer == 5)
        return 6 + 1 + sensor; // 11
      if (layer == 6)
        return 11 + 1 + sensor; // 17
      else
        return -1;
    }

◆ generateHistogramBins()

void generateHistogramBins ( TH2Poly * histogram )

inlineprivate

generate histogram bins

Definition at line 316 of file SVDSummaryPlots.h.

    {
      // parameters
      const double rLayer[4] = {40, 70, 110, 160}; // arbitrary units. True radii are {39, 80, 104, 135}; // mm
      const double nLadders[4] = {7, 10, 12, 16}; // per layer
      const double nSensors[4] = {2, 3, 4, 5}; // per ladder
      const double inclination[4] = {-17, -5, -13, -12}; // de.g. Layer inclination with respect to the x-axis
      const double delta[4] = {9, 8, 8, 8}; // arbitrary units. Width of the bins for each sensor
 
      double pi = TMath::Pi();
 
      for (int layer = 0; layer < 4; layer ++) {
        for (int ladder = 1; ladder <= nLadders[layer]; ladder++) {
          for (int sensor = 1; sensor <= nSensors[layer]; sensor++) {
            double r = rLayer[layer] + (delta[layer]) * (sensor - 1); //position of the layer + position of the sensor
            double phi = 2 * pi / nLadders[layer]; // angle for each sensor (step)
            double dphiThisPoint = (ladder - 1) * phi - phi / 2 + inclination[layer] * pi /
                                   180.; // position of the edge of the sensor (ladder-1 to start from zero)
            double dphiNextPoint = dphiThisPoint + phi; // position of the next edge of the sensor
            // bin coordinates
            double xr0 = r * TMath::Cos(dphiThisPoint);
            double xr1 = (r + delta[layer]) * TMath::Cos(dphiThisPoint);
            double xr2 = (r + delta[layer]) * TMath::Cos(dphiNextPoint);
            double xr3 = r * TMath::Cos(dphiNextPoint);
            double yr0 = r * TMath::Sin(dphiThisPoint);
            double yr1 = (r + delta[layer]) * TMath::Sin(dphiThisPoint);
            double yr2 = (r + delta[layer]) * TMath::Sin(dphiNextPoint);
            double yr3 = r * TMath::Sin(dphiNextPoint);
            // add bin to the histogram
            double xbin[4] = {xr0, xr1, xr2, xr3};
            double ybin[4] = {yr0, yr1, yr2, yr3};
            histogram->AddBin(4, xbin, ybin);
            // resetting bin coordinates
            for (int k = 0; k < 4; k++) {
              xbin[k] = 0.;
              ybin[k] = 0.;
            }
          }
        }
      }
    }

◆ getHistogram()

TH2F * getHistogram ( int view )

inline

get a reference to the histogram for

Parameters

view	please, use the enumeration SVDSummaryPlots::Vindex and SVDSummaryPlots::UIndex

Definition at line 66 of file SVDSummaryPlots.h.

    {
      TH2F* returnValue = m_defaultHistogram;
      try {
        returnValue = m_histos[view];
      } catch (...) {
        B2WARNING("Unexpected view: " << view);
 
        returnValue = m_defaultHistogram;
      }
 
      return returnValue;
    }

◆ getPoly()

TH2Poly * getPoly	(	int	view,
		int	min = `-1111`,
		int	max = `-1111`
	)

inline

Create the TH2Poly version of the plot.

Definition at line 198 of file SVDSummaryPlots.h.

    {
      TH2F* histogram = getHistogram(view);
      TString name = histogram->GetName();
      TString title = histogram->GetTitle();
 
      double xLow = -210.;
      double xHigh = 210.;
      double yLow = -210.;
      double yHigh = 210.;
 
      if (m_polyHistos[view] == NULL) {
        m_polyHistos[view] = new TH2Poly(name + "_poly", title, xLow, xHigh, yLow, yHigh);
        generateHistogramBins(m_polyHistos[view]);
        m_polyHistos[view]->GetXaxis()->SetTitle("x");
        m_polyHistos[view]->GetXaxis()->SetLabelSize(0);
        m_polyHistos[view]->GetXaxis()->SetTickLength(0);
        m_polyHistos[view]->GetXaxis()->SetTitleOffset(0.5);
 
        m_polyHistos[view]->GetYaxis()->SetTitle("y");
        m_polyHistos[view]->GetYaxis()->SetLabelSize(0);
        m_polyHistos[view]->GetYaxis()->SetTickLength(0);
        m_polyHistos[view]->GetYaxis()->SetTitleOffset(0.5);
      }
 
      m_polyHistos[view]->SetName(name + "poly");
      m_polyHistos[view]->SetTitle(title);
 
      m_polyHistos[view]->SetMinimum(min);
      m_polyHistos[view]->SetMaximum(max);
 
      const double nLadders[4] = {7, 10, 12, 16}; // per layer
      const double nSensors[4] = {2, 3, 4, 5}; // per ladder
 
      for (int layer = 0; layer < 4; layer++) {
        for (int ladder = 1; ladder < nLadders[layer] + 1; ladder++) {
          for (int sensor = 1; sensor < nSensors[layer] + 1; sensor++) {
            double value = getValue(layer + 3, ladder, sensor, view);
            int bin = findBinPoly(layer + 3, ladder, sensor);
            m_polyHistos[view]->SetBinContent(bin, value);
          }
        }
      }
 
      return m_polyHistos[view];
    }

◆ getValue() [1/2]

float getValue	(	const VxdID &	vxdID,
		int	view
	)

inline

get the value contained in the corresponding bin, given VxdID and view

Definition at line 81 of file SVDSummaryPlots.h.

    {
      int layer = vxdID.getLayerNumber();
      int ladder = vxdID.getLadderNumber();
      int sensor = vxdID.getSensorNumber();
 
      return getValue(layer, ladder, sensor, view);
    }

◆ getValue() [2/2]

float getValue	(	int	layer,
		int	ladder,
		int	sensor,
		int	view
	)

inline

get the value contained in the corresponding bin, given

Parameters

layer
ladder
sensor
view

Definition at line 96 of file SVDSummaryPlots.h.

    {
      int bin = m_histos[view]->FindBin(ladder, findBinY(layer, sensor));
      return getHistogram(view)->GetBinContent(bin);
    }

◆ reset()

void reset ( )

inline

Reset histograms.

Definition at line 142 of file SVDSummaryPlots.h.

    {
      for (int view = VIndex ; view < UIndex + 1; view++)
        getHistogram(view)->Reset();
    }

◆ setMaximum()

void setMaximum ( Int_t value = 0 )

inline

set histogram maximum

Parameters

value

Definition at line 190 of file SVDSummaryPlots.h.

    {
      for (int view = VIndex ; view < UIndex + 1; view++)
        getHistogram(view)->SetMaximum(value);
    }

◆ setMinimum()

void setMinimum ( Int_t value = 0 )

inline

set histogram minimum

Parameters

value

Definition at line 181 of file SVDSummaryPlots.h.

    {
      for (int view = VIndex ; view < UIndex + 1; view++)
        getHistogram(view)->SetMinimum(value);
    }

◆ setRunID()

void setRunID ( const TString & runID )

inline

set run ids in title

Parameters

runID

Definition at line 129 of file SVDSummaryPlots.h.

    {
      for (int view = VIndex ; view < UIndex + 1; view++) {
        // add blank if needed before adding runID
        if (!m_Title[view].EndsWith(" "))
          m_Title[view].Append(" ");
 
        getHistogram(view)->SetTitle(m_Title[view] + runID);
      }
    }

◆ setStats()

void setStats ( bool stats = true )

inline

set histograms stat

Parameters

stats

Definition at line 151 of file SVDSummaryPlots.h.

    {
      for (int view = VIndex ; view < UIndex + 1; view++)
        getHistogram(view)->SetStats(stats);
    }

Member Data Documentation

◆ m_defaultHistogram

TH2F* m_defaultHistogram = nullptr

private

default histograms

Definition at line 275 of file SVDSummaryPlots.h.

◆ m_histos

TH2F* m_histos[2]

private

vector containing the U and V histograms

Definition at line 271 of file SVDSummaryPlots.h.

◆ m_polyHistos

TH2Poly* m_polyHistos[2] = {nullptr}

private

vector containing the U and V poly-histograms

Definition at line 273 of file SVDSummaryPlots.h.

◆ m_Title

TString m_Title[2]

private

Base title.

Definition at line 277 of file SVDSummaryPlots.h.

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

svd/dataobjects/include/SVDSummaryPlots.h

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ E_side

Constructor & Destructor Documentation

◆ SVDSummaryPlots() [1/2]

◆ SVDSummaryPlots() [2/2]

◆ ~SVDSummaryPlots()

Member Function Documentation

◆ clean()

◆ customize()

◆ customizeString()

◆ fill() [1/3]

◆ fill() [2/3]

◆ fill() [3/3]

◆ findBinPoly()

◆ findBinY()

◆ generateHistogramBins()

◆ getHistogram()

◆ getPoly()

◆ getValue() [1/2]

◆ getValue() [2/2]

◆ reset()

◆ setMaximum()

◆ setMinimum()

◆ setRunID()

◆ setStats()

Member Data Documentation

◆ m_defaultHistogram

◆ m_histos

◆ m_polyHistos

◆ m_Title