SiStripPulseShape.cc File Reference

#include "DQM/SiStripCommissioningAnalysis/interface/SiStripPulseShape.h"
#include <TF1.h>
#include <TMath.h>

Go to the source code of this file.

Functions
double	fdecay (double *x, double *par)
double	fdeconv (double *x, double *par)
double	fdeconv_convoluted (double *x, double *par)
double	fpeak (double *x, double *par)
double	fpeak_convoluted (double *x, double *par)
double	fturnOn (double *x, double *par)
double	pulse (double x, double y, double z, double t)
double	pulse_raw (double x, double y, double z, double t)
double	pulse_x0 (double y, double z, double t)
double	pulse_x0_yz (double z, double t)
double	pulse_yz (double x, double z, double t)

Function Documentation

fdecay()

double fdecay	(	double *	x,
		double *	par
	)

Definition at line 91 of file SiStripPulseShape.cc.

References alignCSCRings::s, and x.

Referenced by CalibrationScanAlgorithm::analyse().

                                      {
  double s = par[0];
  double c_exp = par[1];
  double c_pow = par[2];

  return s * TMath::Exp(x[0] * c_exp) * (1 + x[0] * c_pow);
}

fdeconv()

double fdeconv	(	double *	x,
		double *	par
	)

Definition at line 101 of file SiStripPulseShape.cc.

References fpeak(), x, and geometryCSVtoXML::xz.

Referenced by CalibrationAlgorithm::analyse(), CalibrationScanAlgorithm::analyse(), and SiStripPulseShape::getNormalizedValue().

                                       {
  double xm = par[6] * (x[0] - 25);
  double xp = par[6] * (x[0] + 25);
  double xz = par[6] * x[0];
  return 1.2131 * fpeak(&xp, par) - 1.4715 * fpeak(&xz, par) + 0.4463 * fpeak(&xm, par);
}

fdeconv_convoluted()

double fdeconv_convoluted	(	double *	x,
		double *	par
	)

Definition at line 113 of file SiStripPulseShape.cc.

References fpeak_convoluted(), x, and geometryCSVtoXML::xz.

Referenced by SamplingAlgorithm::SamplingAlgorithm().

                                                  {
  double xm = (x[0] - 25);
  double xp = (x[0] + 25);
  double xz = x[0];
  return 1.2131 * fpeak_convoluted(&xp, par) - 1.4715 * fpeak_convoluted(&xz, par) +
         0.4463 * fpeak_convoluted(&xm, par);
}

fpeak()

double fpeak	(	double *	x,
		double *	par
	)

Definition at line 66 of file SiStripPulseShape.cc.

References pulse(), alignCSCRings::s, submitPVValidationJobs::t, x, geometryCSVtoXML::xx, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by CalibrationAlgorithm::analyse(), fdeconv(), fpeak_convoluted(), and SiStripPulseShape::getNormalizedValue().

                                     {
  double xx = par[0];
  double y = par[1];
  double z = par[2];
  double a_0 = par[3];
  double s = par[4];
  double t_0 = par[5];
  double t = x[0] - t_0;

  // below turn-on time return just a constant
  if (x[0] < t_0)
    return a_0;
  // elswhere return the pulse
  return a_0 + s * pulse(xx, y, z, t);
}

fpeak_convoluted()

double fpeak_convoluted	(	double *	x,
		double *	par
	)

Definition at line 108 of file SiStripPulseShape.cc.

References f, fpeak(), and x.

Referenced by fdeconv_convoluted(), and SamplingAlgorithm::SamplingAlgorithm().

                                                {
  TF1 f("peak_convoluted", fpeak, 0, 250, 4);
  return f.IntegralError(x[0] - par[4] / 2., x[0] + par[4] / 2., par, nullptr, 1.) / (par[4]);
}

fturnOn()

double fturnOn	(	double *	x,
		double *	par
	)

Definition at line 82 of file SiStripPulseShape.cc.

References alignCSCRings::s, ApeEstimator_cff::width, and x.

Referenced by CalibrationScanAlgorithm::analyse().

                                       {
  double a_0 = par[0];
  double s = par[1];
  double t_0 = par[2];
  double width = par[3];

  return a_0 + s * TMath::Erf((x[0] - t_0) / width);
}

pulse()

double pulse	(	double	x,
		double	y,
		double	z,
		double	t
	)

Definition at line 49 of file SiStripPulseShape.cc.

References pulse_raw(), pulse_x0(), pulse_x0_yz(), pulse_yz(), submitPVValidationJobs::t, x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by EcalUncalibRecHitTimingCCAlgo::computeTimeCC(), fpeak(), pos::PixelFEDTestDAC::levelEncoder(), HcalSiPMHitResponse::makeSiPMSignal(), FitterFuncs::PulseShapeFunctor::PulseShapeFunctor(), AbsElectronicODERHS::setInputPulse(), HcalQie::shape(), and EcalUncalibRecHitTimeWeightsAlgo< EBDataFrame >::time().

                                                     {
  if (x > y) {
    double pivot = x;
    x = y;
    y = pivot;
  }

  if ((x == 0) && (y == z))
    return pulse_x0_yz(z, t);
  else if (x == 0)
    return pulse_x0(y, z, t);
  else if (y == z)
    return pulse_yz(x, z, t);
  else
    return pulse_raw(x, y, z, t);
}

pulse_raw()

double pulse_raw	(	double	x,
		double	y,
		double	z,
		double	t
	)

Definition at line 19 of file SiStripPulseShape.cc.

References JetChargeProducer_cfi::exp, funct::pow(), submitPVValidationJobs::t, x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by pulse().

                                                         {
  double result1, result2, result3;

  result1 = z * y * exp(-t / y);
  result1 /= pow(y, 2) - (x + z) * y + z * x;

  result2 = z * x * exp(-t / x);
  result2 /= pow(x, 2) - (x - z) * y - z * x;

  result3 = pow(z, 2) * exp(-t / z);
  result3 /= pow(z, 2) + (x - z) * y - z * x;

  return result1 + result2 + result3;
}

pulse_x0()

double pulse_x0	(	double	y,
		double	z,
		double	t
	)

Definition at line 34 of file SiStripPulseShape.cc.

References JetChargeProducer_cfi::exp, submitPVValidationJobs::t, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by pulse().

{ return z / (y - z) * (exp(-t / y) - exp(-t / z)); }

pulse_x0_yz()

double pulse_x0_yz	(	double	z,
		double	t
	)

Definition at line 47 of file SiStripPulseShape.cc.

References JetChargeProducer_cfi::exp, submitPVValidationJobs::t, and detailsBasic3DVector::z.

Referenced by pulse().

{ return t / z * exp(-t / z); }

pulse_yz()

double pulse_yz	(	double	x,
		double	z,
		double	t
	)

Definition at line 36 of file SiStripPulseShape.cc.

References JetChargeProducer_cfi::exp, submitPVValidationJobs::t, x, and detailsBasic3DVector::z.

Referenced by pulse().

                                              {
  double result1, result2;

  result1 = exp(-t / x) - exp(-t / z);
  result1 *= z * x / (z - x);

  result2 = t * exp(-t / z);

  return (result1 + result2) / (z - x);
}