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	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	get_compensation (double x)
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

double fdeconv	(	double *	x,
		double *	par
	)

Definition at line 102 of file SiStripPulseShape.cc.

References fpeak(), and geometryCSVtoXML::xz.

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

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

double fdeconv_convoluted	(	double *	x,
		double *	par
	)

Definition at line 116 of file SiStripPulseShape.cc.

References fpeak_convoluted(), 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);
}

double fpeak	(	double *	x,
		double *	par
	)

Definition at line 87 of file SiStripPulseShape.cc.

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

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

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

  double y = get_compensation(xx);

  if(x[0] < t_0) return a_0;
  return a_0 + s * pulse(xx, y, z, t);
}

double fpeak_convoluted	(	double *	x,
		double *	par
	)

Definition at line 110 of file SiStripPulseShape.cc.

References f, and fpeak().

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

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

double get_compensation

(

double

x

)

Definition at line 76 of file SiStripPulseShape.cc.

References funct::pow().

Referenced by fpeak().

{
       return 49.9581
              - 1.7941      * x
     - 0.110089    * pow(x,2)
     + 0.0113809   * pow(x,3)
     - 0.000388111 * pow(x,4)
     + 5.9761e-6   * pow(x,5)
     - 3.51805e-8  * pow(x,6);
}

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

Definition at line 58 of file SiStripPulseShape.cc.

References pulse_raw(), pulse_x0(), pulse_x0_yz(), pulse_yz(), vertices_cff::x, and detailsBasic3DVector::y.

Referenced by CSCUpgradeCathodeLCTProcessor::findLCTs(), CSCCathodeLCTProcessor::findLCTs(), fpeak(), pos::PixelFEDTestDAC::levelEncoder(), HcalSiPMHitResponse::makeSiPMSignal(), CSCCathodeLCTProcessor::patternFinding(), CSCUpgradeCathodeLCTProcessor::preTrigger(), CSCCathodeLCTProcessor::preTrigger(), CSCCathodeLCTProcessor::pulseExtension(), 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);
}

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

Definition at line 20 of file SiStripPulseShape.cc.

References JetChargeProducer_cfi::exp, funct::pow(), 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;
}

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

Definition at line 36 of file SiStripPulseShape.cc.

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

Referenced by pulse().

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

double pulse_x0_yz	(	double	z,
		double	t
	)

Definition at line 53 of file SiStripPulseShape.cc.

References JetChargeProducer_cfi::exp.

Referenced by pulse().

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

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

Definition at line 41 of file SiStripPulseShape.cc.

References JetChargeProducer_cfi::exp, and vertices_cff::x.

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);
}