Changeset 2517 in Sophya for trunk/SophyaPI/PIext/pigfd1.cc

Timestamp:

Mar 17, 2004, 7:04:57 PM (22 years ago)

Author:

cmv

Message:

GetDistanceToPoint cmv 17/03/04

File:

• trunk/SophyaPI/PIext/pigfd1.cc (modified) (3 diffs)

trunk/SophyaPI/PIext/pigfd1.cc

@@ -10 +10 @@
 SelectX(0);
 SelectErrBar(false,false);
+NptDraw = 0;
 }
 
@@ -60 +61 @@
 if(GetGraphicAtt().GetLineAtt()!=PI_NotDefLineAtt) SelectErrBar(true,true);
 double xp,yp,er;
+NptDraw = 0;
 for (int i=0; i<mGFD->NData(); i++) {
   xp = mGFD->Absc(VarX,i);
   yp = mGFD->Val(i);
   if(xp<xmin || xp>xmax || yp<ymin || yp>ymax)  continue;
+  NptDraw++;
   if(ErrX) {
     er = mGFD->EAbsc(VarX,i);
@@ -76 +79 @@
 return;
 }
+
+/* --Methode-- */
+double PIGenFitDat::GetDistanceToPoint(double x, double y)
+{
+ if(!mGFD) return 1.e+9;
+ if(mGFD->NData()<=0) return 1.e+9;
+ if(VarX<0) return 1.e+9;
+
+ const int nessai = 100;
+ long inc = (NptDraw>nessai) ? (long)(NptDraw/nessai)+1 : 1;
+
+ double dist = -1.e+18;
+ long n = 0;
+ for(long i=0; i<mGFD->NData(); i++) {
+   double xp=mGFD->Absc(VarX,i);
+   if(xp<XMin() || xp>XMax()) continue;
+   double yp=mGFD->Val(i);
+   if(yp<YMin() || yp>YMax()) continue;
+   if(n%inc==0) {
+     xp = (xp-x)/(XMax()-XMin())/0.5;
+     yp = (yp-y)/(YMax()-YMin())/0.5;
+     xp = xp*xp+yp*yp;
+     if(dist<0. || xp<dist) dist = xp;
+   }
+   n++;
+ }
+ dist=sqrt(fabs(dist));
+ //cout<<"PIGenFitDat: xlim="<<XMin()<<","<<XMax()<<" ylim="<<YMin()<<","<<YMax()
+ //    <<" NData="<<mGFD->NData()<<" inc="<<inc<<endl;
+ //cout<<"....d="<<dist<<" x="<<x<<" y="<<y<<" NptDraw="<<NptDraw<<endl;
+
+ return dist;
+}