The Class for ARICH reconstruction parameters. More...

#include <ARICHReconstructionPar.h>

Inheritance diagram for ARICHReconstructionPar:

Public Member Functions
	ARICHReconstructionPar ()
	Default constructor.

void	initializeDefault ()
	initializes "default" values of parameters

void	setBackgroundPDF (TF1 *bkgPDF)
	Set background PDF function for cherenkov angle distribution (vs.

void	setBackgroundPhiPDF (TF2 *bkgPhiPDF)
	Set background PDF function with flattening phi correction.

void	setThcResolution (TF1 *thcRes)
	Set function for cherenkov angle resolution vs.

void	setAerogelFOM (const std::vector< float > &aerogelFOM)
	Set aerogel figure of merit for cherenkov photon yield.

void	setFlatBkgPerPad (const float flatBkgPerPad)
	Set flat background per pad.

void	setParameters (const std::vector< float > &pars)
	Set additional cherenkov angle PDF parameters.

const std::vector< float > &	getParameters () const
	Get vector of additional cherenkov angle PDF parameters.

const TF1 *	getBackgroundPDF () const
	Get background PDF function.

const TF2 *	getBackgroundPhiPDF () const
	Get background PDF function (with phi correction)

double	getThcResolution (double momentum) const
	Get Cherenkov angle resolution (without smearing due to pad size!) at given track momentum.

double	getAerogelFOM (unsigned iLayer) const
	Get aerogel cherenkov photon yield for tile of selected layer.

float	getFlatBkgPerPad () const
	Get flat background per pad.

double	getBackgroundPerPad (double th_cer, const std::vector< double > &pars) const
	Get expected number of background hits for pad at given theta.

double	getPhiCorrectedBackgroundPerPad (double fi_cer_trk, double th_cer, const std::vector< double > &pars) const
	Get expected number of background hits for pad at given theta at given phi_Cer_trk (flat background in Cherenkov s)pace.

double	getExpectedBackgroundHits (const std::vector< double > &pars, double minThc=0.1, double maxThc=0.5) const
	Get number of expected background hits in ring (0.1<theta<0.5rad by default)

double	getNPadsInRing (double maxThc, double minThc=0.0, double trackTh=0.45) const
	Get average number of pads in ring.

void	print () const
	Print parameters values.

Private Member Functions
	ClassDef (ARICHReconstructionPar, 2)
	ClassDef, must be the last term before the closing {}.

Private Attributes
TF2 *	m_bkgPhiPDF = NULL
	background PDF function with phi correction

TF1 *	m_bkgPDF = NULL
	background PDF function (function of theta)

TF1 *	m_thcResolution = NULL
	cherenkov angle resolution (function of track momentum)

std::vector< float >	m_pars
	vector of other pdf parameters

std::vector< float >	m_aerogelFOM
	aerogel figure of merit (for photon yield)

float	m_flatBkgPerPad = 0.0
	expected background hits per pad (treated flat over detector surface)

Detailed Description

The Class for ARICH reconstruction parameters.

This class provides ARICH parameters for reconstruction (i.e. likelihood calculation)

Definition at line 25 of file ARICHReconstructionPar.h.

Constructor & Destructor Documentation

◆ ARICHReconstructionPar()

ARICHReconstructionPar ( )

inline

Default constructor.

Definition at line 30 of file ARICHReconstructionPar.h.

30{};

Member Function Documentation

◆ getAerogelFOM()

double getAerogelFOM ( unsigned iLayer ) const

inline

Get aerogel cherenkov photon yield for tile of selected layer.

Parameters

iLayer aerogel layer

Returns: aerogel cherenkov photon yield

Definition at line 132 of file ARICHReconstructionPar.h.

    {
      return m_aerogelFOM.at(iLayer);
    }

◆ getBackgroundPDF()

const TF1 * getBackgroundPDF ( ) const

inline

Get background PDF function.

Returns: background PDF function

Definition at line 103 of file ARICHReconstructionPar.h.

    {
      return m_bkgPDF;
    }

◆ getBackgroundPerPad()

double getBackgroundPerPad	(	double	th_cer,
		const std::vector< double > &	pars
	)		const

Get expected number of background hits for pad at given theta.

Parameters

th_cer	pad theta angle
pars	vector of parameters for PDF (beta, track hits HAPD window (1 or 0))

Returns: expected number of backgrond hits on pad

Definition at line 41 of file ARICHReconstructionPar.cc.

{
 
  int ipar = 0;
  for (auto par : pars) {
    m_bkgPDF->SetParameter(ipar, par);
    ipar++;
  }
  return m_bkgPDF->Eval(th_cer) + m_flatBkgPerPad;
}

◆ getBackgroundPhiPDF()

const TF2 * getBackgroundPhiPDF ( ) const

inline

Get background PDF function (with phi correction)

Returns: background PDF function - flat with phi dependance

Definition at line 112 of file ARICHReconstructionPar.h.

    {
      return m_bkgPhiPDF;
    }

◆ getExpectedBackgroundHits()

double getExpectedBackgroundHits	(	const std::vector< double > &	pars,
		double	minThc = `0.1`,
		double	maxThc = `0.5`
	)		const

Get number of expected background hits in ring (0.1<theta<0.5rad by default)

Parameters

pars	vector of parameters for PDF (beta, track hits HAPD window (1 or 0))
minThc	inner theta angle of ring
maxThc	outter theta angle of ring

Returns: expected number of backgrond hits in ring

Definition at line 64 of file ARICHReconstructionPar.cc.

{
 
  int ipar = 0;
  for (auto par : pars) {
    m_bkgPDF->SetParameter(ipar, par);
    ipar++;
  }
 
  // parameters of fit of  pol3 to the number of nPads per ring with given theta and width 5mrad
  const double surf[4] = { -2.19669e-02, 3.59010e+01, -2.77441e+01, 1.43564e+02};
  double step = 0.005;
  double thc = minThc + step / 2.;
  double bkg = 0;
  while (thc < maxThc) {
    bkg += (surf[0] + surf[1] * thc + surf[2] * pow(thc, 2) + surf[3] * pow(thc, 3)) * getBackgroundPerPad(thc, pars);
    thc += step;
  }
  return bkg;
 
}

◆ getFlatBkgPerPad()

float getFlatBkgPerPad ( ) const

inline

Get flat background per pad.

Returns: flat background per pad

Definition at line 141 of file ARICHReconstructionPar.h.

    {
      return m_flatBkgPerPad;
    }

◆ getNPadsInRing()

double getNPadsInRing	(	double	maxThc,
		double	minThc = `0.0`,
		double	trackTh = `0.45`
	)		const

Get average number of pads in ring.

Parameters

maxThc	outter theta angle of ring
minThc	inner theta angle of ring
trackTh	theta angle of track (global coordinates, z axis theta = 0)

Returns: average number of pads in ring

Definition at line 86 of file ARICHReconstructionPar.cc.

{
 
  double s1 = sqrt(tan(maxThc)) * pow((tan(maxThc + trackTh) + tan(maxThc - trackTh)) / 2, 3. / 2.);
  double s2 = 0;
  if (minThc > 0) s2 = sqrt(tan(minThc)) * pow((tan(minThc + trackTh) + tan(minThc - trackTh)) / 2, 3. / 2.);
  return 3.1416 * 0.18 * 0.18 * (s1 - s2) * 0.6 / 0.005 / 0.005; // pi*dist^2*s*avg_geo_acc/pad_size
}

◆ getParameters()

const std::vector< float > & getParameters ( ) const

inline

Get vector of additional cherenkov angle PDF parameters.

Returns: vector of parameters

Definition at line 94 of file ARICHReconstructionPar.h.

    {
      return m_pars;
    }

◆ getPhiCorrectedBackgroundPerPad()

double getPhiCorrectedBackgroundPerPad	(	double	fi_cer_trk,
		double	th_cer,
		const std::vector< double > &	pars
	)		const

Get expected number of background hits for pad at given theta at given phi_Cer_trk (flat background in Cherenkov s)pace.

Parameters

fi_cer_trk	pad phi angle
th_cer	pad theta angle
pars	vector of parameters for PDF (beta, track hits HAPD window (1 or 0))

Returns: expected number of backgrond hits on pad

Definition at line 52 of file ARICHReconstructionPar.cc.

{
 
  int ipar = 0;
  for (auto par : pars) {
    m_bkgPhiPDF->SetParameter(ipar, par);
    ipar++;
  }
  return m_flatBkgPerPad * m_bkgPhiPDF->Eval(fi_cer_trk, th_cer);
}

◆ getThcResolution()

double getThcResolution ( double momentum ) const

inline

Get Cherenkov angle resolution (without smearing due to pad size!) at given track momentum.

Parameters

momentum track momentum

Returns: cherenkov angle resolution

Definition at line 122 of file ARICHReconstructionPar.h.

    {
      return  m_thcResolution->Eval(momentum);
    }

◆ initializeDefault()

void initializeDefault ( )

initializes "default" values of parameters

Definition at line 18 of file ARICHReconstructionPar.cc.

{
 
  double parsBkg[10] = {0, 1, 8.764830e+00, -4.815370e+00, 3.268570e-06, 8.518560e+00, -9.609710e+00, 9.769060e-06, 1.095350e-02, -1.322410e-02};
  m_bkgPDF = new TF1("bkgPDF",
                     "0*[0] + (1-[1])*(exp([2]+[3]*x)*[4]) + [1]*(exp([5]+[6]*x)*[7]+x*[8]+x*x*[9])");
  m_bkgPDF->SetParameters(parsBkg);
 
  double parsRes[4] = {1.12147e-01, 3.81701e+00, 0.0093, -6.14984e-04};
  m_thcResolution = new TF1("thcResolution", "([p0]*exp(-[p1]*x)+[p2]+[p3]*x-((x>3)?[p3]*(x-3):0))*1.35");
  m_thcResolution->SetParameters(parsRes);
 
  double parsBkgPhi[8] = {0, 0, 6.68258e-02, 6.50048e+00, -1.04876e-01, -7.06652e-03, 1.34111e-01, 1.83823e-01};
  m_bkgPhiPDF = new TF2("m_bkgPhiPDF",
                        "[0]*0 + [1]*0 + 1/( ([2]*exp(y*[3]) + [4])*cos(x/2)*cos(x/2)*cos(x/2)*cos(x/2) + sqrt([5] + [6]*y + [7]*y*y ) )");
  m_bkgPhiPDF->SetParameters(parsBkgPhi);
 
  m_pars = {0.098304138, 5.973590e-02, 21, 0.88, 3.7, 0.357, 0.46, 1, 0.12, 0.12, 0.02, 1.5, 0.0288};
  m_flatBkgPerPad = 0.000892;
 
  m_aerogelFOM = {13.200797, 15.186757};
}

◆ print()

void print ( ) const

Print parameters values.

Definition at line 95 of file ARICHReconstructionPar.cc.

{
 
  cout << endl << "-----bkg PDF-----" << endl;
  m_bkgPDF->Print();
  int Npar = m_bkgPDF->GetNpar();
  std::vector<double> bkgPars(Npar, 0);
  m_bkgPDF->GetParameters(bkgPars.data());
  for (int i = 0; i < Npar; i++)cout << Form("bkg Pars %d = %e", i, bkgPars[i]) << endl;
 
  cout << endl << "-----bkg phi corr PDF-----" << endl;
  m_bkgPhiPDF->Print();
  int Npar2 = m_bkgPhiPDF->GetNpar();
  std::vector<double> bkgPhiCorPars(Npar2, 0);
  m_bkgPhiPDF->GetParameters(bkgPhiCorPars.data());
  for (int i = 0; i < Npar2; i++)cout << Form("bkg Pars %d = %e", i, bkgPhiCorPars[i]) << endl;
 
  cout << endl << "-----flat backgroud per pad-----"  << endl;
  cout << " flat background per pad is " << m_flatBkgPerPad << endl;
 
  cout << endl << "----additional parameters (for wide gaus, quartz etc.)-----"  << endl;
  for (int i = 0; i < (int)m_pars.size(); i++) {
    printf("m_pars[%d] = %e\n", i, m_pars[i]);
  }
 
  cout << endl << "----Aerogel FOM-----"  << endl;
  for (int i = 0; i < (int)m_aerogelFOM.size(); i++) {
    printf("m_aerogelFOM[%d] = %f\n", i, m_aerogelFOM[i]);
  }
 
  cout << endl << "----- Resolution PDF-----" << endl;
  m_thcResolution->Print();
  double resPars[4];
  m_thcResolution->GetParameters(resPars);
  for (int i = 0; i < 4; i++)cout << Form("resolution Pars %d = %e", i, resPars[i]) << endl;
 
}

◆ setAerogelFOM()

void setAerogelFOM ( const std::vector< float > & aerogelFOM )

inline

Set aerogel figure of merit for cherenkov photon yield.

Parameters

aerogelFOM vector of aerogel layers FOMs

Definition at line 67 of file ARICHReconstructionPar.h.

    {
      m_aerogelFOM =  aerogelFOM;
    }

◆ setBackgroundPDF()

void setBackgroundPDF ( TF1 * bkgPDF )

inline

Set background PDF function for cherenkov angle distribution (vs.

theta)

Parameters

bkgPDF background PDF

Definition at line 39 of file ARICHReconstructionPar.h.

    {
      m_bkgPDF = bkgPDF;
    }

◆ setBackgroundPhiPDF()

void setBackgroundPhiPDF ( TF2 * bkgPhiPDF )

inline

Set background PDF function with flattening phi correction.

Parameters

bkgPhiPDF background PDF f(phi, theta)

Definition at line 48 of file ARICHReconstructionPar.h.

    {
      m_bkgPhiPDF = bkgPhiPDF;
    }

◆ setFlatBkgPerPad()

void setFlatBkgPerPad ( const float flatBkgPerPad )

inline

Set flat background per pad.

Parameters

flatBkgPerPad value of flat background per pad

Definition at line 76 of file ARICHReconstructionPar.h.

    {
      m_flatBkgPerPad = flatBkgPerPad;
    }

◆ setParameters()

void setParameters ( const std::vector< float > & pars )

inline

Set additional cherenkov angle PDF parameters.

Parameters

pars	vector of parameters

Definition at line 85 of file ARICHReconstructionPar.h.

    {
      m_pars =  pars;
    }

◆ setThcResolution()

void setThcResolution ( TF1 * thcRes )

inline

Set function for cherenkov angle resolution vs.

momentum

Parameters

thcRes cherenkov angle resolution vs. momentum function

Definition at line 58 of file ARICHReconstructionPar.h.

    {
      m_thcResolution = thcRes;
    }

Member Data Documentation

◆ m_aerogelFOM

std::vector<float> m_aerogelFOM

private

aerogel figure of merit (for photon yield)

Definition at line 194 of file ARICHReconstructionPar.h.

◆ m_bkgPDF

TF1* m_bkgPDF = NULL

mutableprivate

background PDF function (function of theta)

Definition at line 191 of file ARICHReconstructionPar.h.

◆ m_bkgPhiPDF

TF2* m_bkgPhiPDF = NULL

mutableprivate

background PDF function with phi correction

Definition at line 190 of file ARICHReconstructionPar.h.

◆ m_flatBkgPerPad

float m_flatBkgPerPad = 0.0

private

expected background hits per pad (treated flat over detector surface)

Definition at line 195 of file ARICHReconstructionPar.h.

◆ m_pars

std::vector<float> m_pars

private

vector of other pdf parameters

Definition at line 193 of file ARICHReconstructionPar.h.

◆ m_thcResolution

TF1* m_thcResolution = NULL

private

cherenkov angle resolution (function of track momentum)

Definition at line 192 of file ARICHReconstructionPar.h.

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

arich/dbobjects/include/ARICHReconstructionPar.h
arich/dbobjects/src/ARICHReconstructionPar.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ARICHReconstructionPar()

Member Function Documentation

◆ getAerogelFOM()

◆ getBackgroundPDF()

◆ getBackgroundPerPad()

◆ getBackgroundPhiPDF()

◆ getExpectedBackgroundHits()

◆ getFlatBkgPerPad()

◆ getNPadsInRing()

◆ getParameters()

◆ getPhiCorrectedBackgroundPerPad()

◆ getThcResolution()

◆ initializeDefault()

◆ print()

◆ setAerogelFOM()

◆ setBackgroundPDF()

◆ setBackgroundPhiPDF()

◆ setFlatBkgPerPad()

◆ setParameters()

◆ setThcResolution()

Member Data Documentation

◆ m_aerogelFOM

◆ m_bkgPDF

◆ m_bkgPhiPDF

◆ m_flatBkgPerPad

◆ m_pars

◆ m_thcResolution