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

double fdecay	(	double *	x,
		double *	par
	)

Definition at line 95 of file SiStripPulseShape.cc.

References alignCSCRings::s, and vertices_cff::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);
}

double fdeconv	(	double *	x,
		double *	par
	)

Definition at line 107 of file SiStripPulseShape.cc.

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

double fdeconv_convoluted	(	double *	x,
		double *	par
	)

Definition at line 119 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 69 of file SiStripPulseShape.cc.

References pulse(), alignCSCRings::s, lumiQTWidget::t, 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);
}

double fpeak_convoluted	(	double *	x,
		double *	par
	)

Definition at line 114 of file SiStripPulseShape.cc.

References f, and fpeak().

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

double fturnOn	(	double *	x,
		double *	par
	)

Definition at line 85 of file SiStripPulseShape.cc.

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

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

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

Definition at line 55 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 51 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 40 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);
}