myFunctions.cc File Reference

Go to the source code of this file.

Functions
double	crystalball (double *x, double *par)
double	ene1x1_xtal (double mymatrix[], int xtal)
double	ene3x3_xtal (double mymatrix[], int xtal)
double	ene5x5_xtal (double mymatrix[], int xtal)
double	ene7x7_xtal (double mymatrix[])
void	energy_ieta (double mymatrix[], double *energyieta)
void	energy_iphi (double mymatrix[], double *energyiphi)
int	maxAmplitInMatrix (double myMatrix[])
double	shapeFunction (double *x, double *par)

Function Documentation

double crystalball	(	double *	x,
		double *	par
	)

Definition at line 98 of file myFunctions.cc.

References JetChargeProducer_cfi::exp, dedxEstimators_cff::exponent, and funct::pow().

                                           {
  // par[0]:  mean
  // par[1]:  sigma
  // par[2]:  alpha, crossover point
  // par[3]:  n, length of tail
  // par[4]:  N, normalization
                                
  double cb = 0.0;
  double exponent = 0.0;
  double bla = 0.0;
  
  if (x[0] > par[0] - par[2]*par[1]) {
    exponent = (x[0] - par[0])/par[1];
    cb = exp(-exponent*exponent/2.);
  } else {
    double nenner  = pow(par[3]/par[2], par[3])*exp(-par[2]*par[2]/2.);
    double zaehler = (par[0] - x[0])/par[1] + par[3]/par[2] - par[2];
    zaehler = pow(zaehler, par[3]);
    cb = nenner/zaehler;
  }
  
  if (par[4] > 0.) {
    cb *= par[4];
  }
  return cb;
}

double ene1x1_xtal	(	double	mymatrix[],
		int	xtal
	)

Definition at line 21 of file myFunctions.cc.

{
  double E1x1 = 0.;
  if (mymatrix[xtal]<-50) { E1x1 = -1000.; }
  else { E1x1 = mymatrix[xtal]; }
  return E1x1;
}

double ene3x3_xtal	(	double	mymatrix[],
		int	xtal
	)

Definition at line 30 of file myFunctions.cc.

{
  double E3x3 = 0.;

  if ( (mymatrix[xtal-8]<-50) || (mymatrix[xtal-7]<-50)  || (mymatrix[xtal-6]<-50)  || (mymatrix[xtal-1]<-50) || (mymatrix[xtal]<-50)  || (mymatrix[xtal+1]<-50)  || (mymatrix[xtal+6]<-50) || (mymatrix[xtal+7]<-50)  || (mymatrix[xtal+8]<-50) ) 
    { E3x3 = -1000.; }
  else 
    { E3x3 = mymatrix[xtal-8] + mymatrix[xtal-7] + mymatrix[xtal-6] + mymatrix[xtal-1] + mymatrix[xtal] + mymatrix[xtal+1] + mymatrix[xtal+6] + mymatrix[xtal+7] + mymatrix[xtal+8]; }

  return E3x3;
}

double ene5x5_xtal	(	double	mymatrix[],
		int	xtal
	)

Definition at line 43 of file myFunctions.cc.

{
  double E5x5 = 0.;

  if( (mymatrix[xtal-16]<-50) || (mymatrix[xtal-15]<-50) || (mymatrix[xtal-14]<-50) || (mymatrix[xtal-13]<-50) || (mymatrix[xtal-12]<-50) || (mymatrix[xtal-9]<-50) || (mymatrix[xtal-8]<-50) || (mymatrix[xtal-7]<-50) || (mymatrix[xtal-6]<-50) || (mymatrix[xtal-5]<-50) || (mymatrix[xtal-2]<-50) || (mymatrix[xtal-1]<-50) || (mymatrix[xtal]<-50) || (mymatrix[xtal+1]<-50) || (mymatrix[xtal+2]<-50) || (mymatrix[xtal+5]<-50) || (mymatrix[xtal+6]<-50) || (mymatrix[xtal+7]<-50) || (mymatrix[xtal+8]<-50) || (mymatrix[xtal+9]<-50) || (mymatrix[xtal+12]<-50) || (mymatrix[xtal+13]<-50) || (mymatrix[xtal+14]<-50) || (mymatrix[xtal+15]<-50) || (mymatrix[xtal+16]<-50) )
    { E5x5 = -1000.; }
  else
    { E5x5 = mymatrix[xtal-16] + mymatrix[xtal-15] + mymatrix[xtal-14] + mymatrix[xtal-13] + mymatrix[xtal-12] + mymatrix[xtal-9] + mymatrix[xtal-8] + mymatrix[xtal-7] + mymatrix[xtal-6] + mymatrix[xtal-5] + mymatrix[xtal-2] + mymatrix[xtal-1] + mymatrix[xtal] + mymatrix[xtal+1] + mymatrix[xtal+2] + mymatrix[xtal+5] + mymatrix[xtal+6] + mymatrix[xtal+7] + mymatrix[xtal+8] + mymatrix[xtal+9] + mymatrix[xtal+12] + mymatrix[xtal+13] + mymatrix[xtal+14] + mymatrix[xtal+15] + mymatrix[xtal+16]; }

  return E5x5;
}

double ene7x7_xtal

(

double

mymatrix[]

)

Definition at line 77 of file myFunctions.cc.

References cuy::ii.

{
  double E7x7 = 0.;
  
  for (int ii=0;ii<49;ii++)
    {
      if (mymatrix[ii]<-50)
    E7x7 = -1000.;
      else
    E7x7 += mymatrix[ii];
    }

  return E7x7;
}

void energy_ieta	(	double	mymatrix[],
		double *	energyieta
	)

Definition at line 56 of file myFunctions.cc.

References cuy::ii, and findQualityFiles::jj.

{
  for(int jj = 0 ; jj < 5 ; jj++) energyieta[jj] = 0.0; 

  for(int jj = 0 ; jj < 5 ; jj++){
    for (int ii = 0; ii < 5 ; ii++) energyieta[jj] += mymatrix[(8+jj)+7*ii]; 
  }
}

void energy_iphi	(	double	mymatrix[],
		double *	energyiphi
	)

Definition at line 66 of file myFunctions.cc.

References cuy::ii, and findQualityFiles::jj.

{
  for(int jj = 0 ; jj < 5 ; jj++) energyiphi[jj] = 0.0; 

  for(int jj = 0 ; jj < 5 ; jj++){
    for (int ii = 0; ii < 5 ; ii++) energyiphi[jj] += mymatrix[(8+7*jj)+ii]; 
  }
}

int maxAmplitInMatrix

(

double

myMatrix[]

)

Definition at line 6 of file myFunctions.cc.

{
  int maxXtal   = 999;
  double maxADC = -999.;
  
  for (int icry=0; icry<49; icry++) {
    if (myMatrix[icry] > maxADC){
      maxADC  = myMatrix[icry];
      maxXtal = icry;
    }}
  
  return maxXtal;
}

double shapeFunction	(	double *	x,
		double *	par
	)

Definition at line 125 of file myFunctions.cc.

                                             {

  // par[0] = constant value
  // par[1], par[2], par[3]: a3, a2 x>0 pol3  
  // par[4], par[5], par[6]: a3, a2 x<=0 pol3  
  
  double ret_val;
  if(x[0]>0.)
    ret_val = par[0] + par[1]*x[0]*x[0] + par[2]*x[0]*x[0]*x[0] + par[3]*x[0]*x[0]*x[0]*x[0];
  else
    ret_val = par[0] + par[4]*x[0]*x[0] + par[5]*x[0]*x[0]*x[0] + par[6]*x[0]*x[0]*x[0]*x[0];

  return ret_val;
}