Plus compliqué: gaussienne 2D

// Fonction de fit Gaussienne 2D

Double_t Gaus2D(Double_t *x, Double_t *par)

{

if(par[2] > 0 && par[4] > 0)

{

double rx=(x[0]-par[1])/par[2];

double ry=(x[1]-par[3])/par[4];

return par[0]*TMath::Exp(-(rx*rx+ry*ry)/2.);

}

else

return 0.;

}

Tableau des paramètres

x[0] : x par[0] : Constante

x[1] : y par[1] : X0

par[2] : sx

par[3] : Y0

par[4] : sy

Tableau d'arguments

cp /home/mra1/Jour3/Gaus2D.C ./