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