void f2gauss5(double x,double a[],double *y,double dyda[],int na)
{
- //Function describing a sum of 2D gaussians with 5 parameters.
- //There number of gaussians is na/5.
+ /*Function describing a sum of 2D gaussians with 5 parameters.
+ There number of gaussians is na/5. */
int i,index;
double fac,fac1,fac2,ex,ex1,ex2,arg1,arg2,u,v;
- //printf("fitting na %d, with pad %f time %f amplitude %f padwidth %f timewidth %f\n",na,a[2],a[4],a[1],a[3],a[5]);
+ /*printf("fitting na %d, with pad %f time %f amplitude %f padwidth %f timewidth %f\n",na,a[2],a[4],a[1],a[3],a[5]);*/
index = nint(x);
if( index < 0 || index >=FIT_PTS )
{
}
u = plane[index].u;
v = plane[index].v;
- //printf("u %f v %f\n",u,v);
+ /*printf("u %f v %f\n",u,v);*/
*y=0.0;
for (i=1;i<=na-1;i+=5)
{
void nrerror(char error_text[])
/* Numerical Recipes standard error handler */
{
- //printf("%s\n",error_text);
- //exit(1);
+ /*printf("%s\n",error_text);
+ exit(1);*/
}
void free_vector(double *v, long nl, long nh)
return 0;
}
else {
- //if( control_g.print_fit_errors==2 )
+ /*if( control_g.print_fit_errors==2 )*/
fprintf( stderr, " runtime error\n" );
free_matrix(alpha,1,MA,1,MA);
-// @(#) $Id$
+/*@(#) $Id$*/
#ifndef AliL3FitUtilities
#define AliL3FitUtilities
-//This we do because this file is read both with c and c++ compiler,
-//and extern "C" is needed only in case of c++.
+/*This we do because this file is read both with c and c++ compiler,
+ and extern "C" is needed only in case of c++. */
#ifdef __cplusplus
extern "C"
#endif
#define NR_END 1
#define FREE_ARG char*
#define EPSILON 1.0E-12
-//#define TRUE 1
-//#define FALSE 1
#define FIT_PTS 2000
#define FIT_MAXPAR 41
#define NUM_PARS 5
};
typedef struct DPOINT DPOINT;
-extern DPOINT *plane; //!
+extern DPOINT *plane;
typedef struct {
long float_size;