Class computes probability density for Muid calculation. More...

#include <B2BIIFixMdstModule.h>

MuidProb methods:
static const int	kRange = 16
	Array size of range.

static const int	kRchisq = 50
	Array size of reduced chi-squared.

static const int	kEEcl = 100
	Array size of ECL energy.

static const int	kPTrk = 50
	Array size of CDC momentum.

static const double	kRchisqMax = 10.0
	Overflow value of reduced chi-squared.

static const double	kEEclMax = 1.0
	Overflow value of ECL energy.

static const double	kPTrkMax = 5.0
	Overflow value of CDC momentum.

double	fRange [2][4][15][kRange]
	Range pdf.

double	fRchisq [2][4][kRchisq+1]
	Reduced chi-squared pdf.

double	fRchisqN [2][4][15]
	Non-overflow normalization.

double	fRchisqD1 [2][4][kRchisq+1] = {{{0.0}}}
	First derivatives of Reduced chi-squared pdf.

double	fRchisqD2 [2][4][kRchisq+1] = {{{0.0}}}
	Second derivatives of Reduced chi-squared pdf.

double	fRchisqD3 [2][4][kRchisq+1] = {{{0.0}}}
	Third derivatives of Reduced chi-squared pdf.

double	fEEcl [kPTrk][kEEcl+1]
	ECL energy pdf.

double	fEEclD1 [kPTrk][kEEcl+1] = {{0.0}}
	First derivatives of ECL energy pdf.

double	fEEclD2 [kPTrk][kEEcl+1] = {{0.0}}
	Second derivatives of ECL energy pdf.

double	fEEclD3 [kPTrk][kEEcl+1] = {{0.0}}
	Third derivatives of ECL energy pdf.

	MuidProb (const char *, int &)
	Constructor.

	~MuidProb ()
	Destructor.

double	prob (int, int, int, int, double) const
	Compute probability density.

double	probRange (int, int, int, int) const
	Compute probability density for range.

double	probRchisq (int, int, int, double) const
	Compute probability density for reduced chi-squared.

double	probECL (double, double) const
	Compute probability density for ECL energy deposit.

void	readDB (const char *, int &)
	Read in probability density functions from database.

Detailed Description

Class computes probability density for Muid calculation.

Definition at line 285 of file B2BIIFixMdstModule.h.

Constructor & Destructor Documentation

◆ MuidProb()

MuidProb	(	const char *	dbtemplate,
		int &	expNo
	)

Constructor.

Definition at line 287 of file B2BIIFixMdstModule_muid.cc.

{
 
  readDB(dbtemplate, expNo);
 
}

◆ ~MuidProb()

~MuidProb ( )

inline

Destructor.

Definition at line 294 of file B2BIIFixMdstModule.h.

294{};

Member Function Documentation

◆ prob()

double prob	(	int	ECMaxLyr,
		int	outcome,
		int	lyr_ext,
		int	lyr_dif,
		double	chi2_red
	)		const

Compute probability density.

Definition at line 376 of file B2BIIFixMdstModule_muid.cc.

{
 
  // ECMaxLyr: 0 if Endcap_MX_layer == 13; 1 if Endcap_MX_layer == 11
  // outcome:  0=??, 1=Barrel Stop, 2=Endcap Stop, 3=Barrel Exit, 4=Endcap Exit
  // lyr_ext:  last layer that Ext track touched
  // lyr_dif:  difference between last Ext layer and last hit layer
  // chi2_red: reduced chi**2 of the transverse deviations of all associated
  //           hits from the corresponding Ext track crossings
 
  return probRange(ECMaxLyr, outcome, lyr_ext, lyr_dif) *
         probRchisq(ECMaxLyr, outcome, lyr_ext, chi2_red);
 
}

◆ probECL()

double probECL	(	double	eEcl,
		double	p
	)		const

Compute probability density for ECL energy deposit.

Definition at line 456 of file B2BIIFixMdstModule_muid.cc.

{
 
  // eEcl: ECL energy (GeV)
  // p:    CDC momentum (GeV/c)
 
  double prob;
  int jTrk = (p < kPTrkMax) ? (int)(p / (kPTrkMax / kPTrk)) : kPTrk - 1;
  if (eEcl >= kEEclMax) {
    prob = fEEcl[jTrk][kEEcl];       // overflow bin
  } else {
    int    i  = (int)(eEcl / (kEEclMax / kEEcl));
    double dx = eEcl - i * (kEEclMax / kEEcl);
    prob = fEEcl[jTrk][i] +
           dx * (fEEclD1[jTrk][i] +
                 dx * (fEEclD2[jTrk][i] +
                       dx * fEEclD3[jTrk][i]));
    if (prob < 0.0) { prob = 0.0; }     // don't let spline fit go negative
  }
  return prob;
 
}

◆ probRange()

double probRange	(	int	ECMaxLyr,
		int	outcome,
		int	lyr_ext,
		int	lyr_dif
	)		const

Compute probability density for range.

Definition at line 394 of file B2BIIFixMdstModule_muid.cc.

{
 
 
  // ECMaxLyr: 0 if max endcap layer is 13; 1 if it's 11.
  // outcome:  0=??, 1=Barrel Stop, 2=Endcap Stop, 3=Barrel Exit, 4=Endcap Exit
  // lyr_ext:  last layer that Ext track touched
  // lyr_dif:  difference between last Ext layer and last hit layer
 
  if (outcome <=  0) { return 0.0; }
  if (outcome >   4) { return 0.0; }
  if (lyr_ext <   0) { return 0.0; }
  if (lyr_ext >  14) { return 0.0; }
  if (lyr_dif <   0) { return 0.0; }
  if (lyr_dif >= kRange) { lyr_dif = kRange - 1; }
 
  return fRange[ECMaxLyr][outcome - 1][lyr_ext][lyr_dif];
 
}

◆ probRchisq()

double probRchisq	(	int	ECMaxLyr,
		int	outcome,
		int	lyr_ext,
		double	chi2_red
	)		const

Compute probability density for reduced chi-squared.

Definition at line 417 of file B2BIIFixMdstModule_muid.cc.

{
 
  // ECMaxLyr: 0 if max endcap layer is 13; 1 if it's 11.
  // outcome:  0=??, 1=Barrel Stop, 2=Endcap Stop, 3=Barrel Exit, 4=Endcap Exit
  // lyr_ext:  last layer that Ext track touched
  // chi2_red: reduced chi**2 of the transverse deviations of all associated
  //           hits from the corresponding Ext track crossings
 
  // Extract the probability density for a particular value of reduced
  // chi-squared by spline interpolation of the binned histogram values.
  // This avoids binning artifacts that were seen when the histogram
  // was sampled directly.
 
  if (outcome  <= 0) { return 0.0; }
  if (outcome  >  4) { return 0.0; }
  if (lyr_ext  <  0) { return 0.0; }
  if (lyr_ext  > 14) { return 0.0; }
  if (chi2_red <  0.0) { return 0.0; }
 
  double prob;
  double area = fRchisqN[ECMaxLyr][outcome - 1][lyr_ext];
  if (chi2_red >= kRchisqMax) {
    prob = area * fRchisq[ECMaxLyr][outcome - 1][kRchisq - 1] + (1.0 - area);
  } else {
    int    i  = (int)(chi2_red / (kRchisqMax / kRchisq));
    double dx = chi2_red - i * (kRchisqMax / kRchisq);
    prob = fRchisq[ECMaxLyr][outcome - 1][i] +
           dx * (fRchisqD1[ECMaxLyr][outcome - 1][i] +
                 dx * (fRchisqD2[ECMaxLyr][outcome - 1][i] +
                       dx * fRchisqD3[ECMaxLyr][outcome - 1][i]));
    prob = (prob < 0.0) ? 0.0 : area * prob;
  }
  return prob;
}

◆ readDB()

void readDB	(	const char *	dbtemplate,
		int &	expNo
	)

Read in probability density functions from database.

Definition at line 294 of file B2BIIFixMdstModule_muid.cc.

{
 
  double dx;    // bin size
 
  if (expNo == 0) return;   // do nothing if pdf file was not found previously
 
  if ((expNo < 5) || ((expNo & 1) == 0)) {
    expNo = 5;    // Runs 1-4 and all even-numbered runs map to Exp #5
  }
 
  char dbname[] = "muid_xxxx_e000000.dat";
  std::string pathname = "";
  bool tmp = Belle::set_belfnm_verbose(false);
  while ((pathname == "") && (expNo >= 5)) {
    std::sprintf(dbname, "%s%s%s%06d%s", "muid_", dbtemplate,
                 "_e", expNo, ".dat");
    pathname = Belle::belfnm(dbname, 0, "share/belle_legacy/data-files/muid");
    expNo -= 2;
  }
  (void)Belle::set_belfnm_verbose(tmp);
  expNo += 2;
  if (pathname == "") {
    B2ERROR("%MuidProb: Failed to open database file."
            << "Module B2BIIFixMdst muid will do nothing.");
    expNo = 0;
    return;
  }
 
  std::ifstream probDB;
  char ident[90];
  probDB.open(pathname.c_str());
  probDB.get(ident, 80, '\n');
  B2INFO("%MuidProb: Reading file " << pathname
         << std::endl << "%MuidProb: " << ident
         << std::endl << "%MuidProb: ");
 
  for (int l = 0; l < 2; l++) {
    B2INFO("range ");
    for (int k = 0; k < 4; k++) {
      for (int j = 0; j < 15; j++) {
        for (int i = 0; i < kRange; i++) {
          probDB >> fRange[l][k][j][i];
        }
      }
    }
    B2INFO("chi**2 ");
    dx = kRchisqMax / kRchisq;
    for (int k = 0; k < 4; k++) {
      for (int i = 0; i < kRchisq; i++) {
        probDB >> fRchisq[l][k][i];
      }
      Spline::muidSpline(kRchisq, dx, &fRchisq[l][k][0], &fRchisqD1[l][k][0],
                         &fRchisqD2[l][k][0], &fRchisqD3[l][k][0]);
    }
    for (int k = 0; k < 4; k++) {
      for (int j = 0; j < 15; j++) {
        probDB >> fRchisqN[l][k][j];
      }
    }
  }
  B2INFO("eEcl ");
  dx = kEEclMax / kEEcl;
  for (int k = 0; k < kPTrk; k++) {
    double probSum = 0.0; // normalization denominator
    for (int i = 0; i < kEEcl; i++) {
      probDB >> fEEcl[k][i];
      probSum += fEEcl[k][i];
    }
    fEEcl[k][kEEcl] = 1.0 - probSum;            // overflow bin
    if (fEEcl[k][kEEcl] < 0.0) { fEEcl[k][kEEcl] = 0.0; }
    if (fEEcl[k][kEEcl] > 1.0) { fEEcl[k][kEEcl] = 1.0; }
    Spline::muidSpline(kEEcl, dx, &fEEcl[k][0], &fEEclD1[k][0],
                       &fEEclD2[k][0], &fEEclD3[k][0]);
  }
 
  probDB.close();
 
}

Member Data Documentation

◆ fEEcl

double fEEcl[kPTrk][kEEcl+1]

private

ECL energy pdf.

Definition at line 334 of file B2BIIFixMdstModule.h.

◆ fEEclD1

double fEEclD1[kPTrk][kEEcl+1] = {{0.0}}

private

First derivatives of ECL energy pdf.

Definition at line 335 of file B2BIIFixMdstModule.h.

◆ fEEclD2

double fEEclD2[kPTrk][kEEcl+1] = {{0.0}}

private

Second derivatives of ECL energy pdf.

Definition at line 336 of file B2BIIFixMdstModule.h.

◆ fEEclD3

double fEEclD3[kPTrk][kEEcl+1] = {{0.0}}

private

Third derivatives of ECL energy pdf.

Definition at line 337 of file B2BIIFixMdstModule.h.

◆ fRange

double fRange[2][4][15][kRange]

private

Range pdf.

Definition at line 328 of file B2BIIFixMdstModule.h.

◆ fRchisq

double fRchisq[2][4][kRchisq+1]

private

Reduced chi-squared pdf.

Definition at line 329 of file B2BIIFixMdstModule.h.

◆ fRchisqD1

double fRchisqD1[2][4][kRchisq+1] = {{{0.0}}}

private

First derivatives of Reduced chi-squared pdf.

Definition at line 331 of file B2BIIFixMdstModule.h.

◆ fRchisqD2

double fRchisqD2[2][4][kRchisq+1] = {{{0.0}}}

private

Second derivatives of Reduced chi-squared pdf.

Definition at line 332 of file B2BIIFixMdstModule.h.

◆ fRchisqD3

double fRchisqD3[2][4][kRchisq+1] = {{{0.0}}}

private

Third derivatives of Reduced chi-squared pdf.

Definition at line 333 of file B2BIIFixMdstModule.h.

◆ fRchisqN

double fRchisqN[2][4][15]

private

Non-overflow normalization.

Definition at line 330 of file B2BIIFixMdstModule.h.

◆ kEEcl

const int kEEcl = 100

staticprivate

Array size of ECL energy.

Definition at line 313 of file B2BIIFixMdstModule.h.

◆ kEEclMax

const double kEEclMax = 1.0

staticprivate

Overflow value of ECL energy.

Definition at line 319 of file B2BIIFixMdstModule.h.

◆ kPTrk

const int kPTrk = 50

staticprivate

Array size of CDC momentum.

Definition at line 314 of file B2BIIFixMdstModule.h.

◆ kPTrkMax

const double kPTrkMax = 5.0

staticprivate

Overflow value of CDC momentum.

Definition at line 320 of file B2BIIFixMdstModule.h.

◆ kRange

const int kRange = 16

staticprivate

Array size of range.

Definition at line 311 of file B2BIIFixMdstModule.h.

◆ kRchisq

const int kRchisq = 50

staticprivate

Array size of reduced chi-squared.

Definition at line 312 of file B2BIIFixMdstModule.h.

◆ kRchisqMax

const double kRchisqMax = 10.0

staticprivate

Overflow value of reduced chi-squared.

Definition at line 318 of file B2BIIFixMdstModule.h.

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

b2bii/modules/B2BIIMdstInput/include/B2BIIFixMdstModule.h
b2bii/modules/B2BIIMdstInput/src/B2BIIFixMdstModule_muid.cc

MuidProb methods:

Detailed Description

Constructor & Destructor Documentation

◆ MuidProb()

◆ ~MuidProb()

Member Function Documentation

◆ prob()

◆ probECL()

◆ probRange()

◆ probRchisq()

◆ readDB()

Member Data Documentation

◆ fEEcl

◆ fEEclD1

◆ fEEclD2

◆ fEEclD3

◆ fRange

◆ fRchisq

◆ fRchisqD1

◆ fRchisqD2

◆ fRchisqD3

◆ fRchisqN

◆ kEEcl

◆ kEEclMax

◆ kPTrk

◆ kPTrkMax

◆ kRange

◆ kRchisq

◆ kRchisqMax