boolean pick=false;
boolean picked=false; 
float m=0.25;
boolean showRefined=false;
Body B = new Body(5);

class Body {
   int maxnb=100;  
   int nb=0;  // number of balls
   int b=0; // index of picked ball
   float r[] = new float [maxnb];  // radii of balls
   pt C[] = new pt [maxnb];    // centers of balls 
   Body() {for (int i=0; i<maxnb; i++) C[i]=new pt();}
   Body(int n) {for (int i=0; i<maxnb; i++) C[i]=new pt(); nb=n; for (int i=0; i<nb; i++) {C[i].setTo(0,-i*dim/nb,0); r[i]=random(dim/10);};}
   void show() {pushMatrix(); for (float a=0; a<=TWO_PI; a+=PI/2) {showBalls(); showCaplets(24); rotateY(a);}; popMatrix();}
 
  void reset() {nb=0;}
  void add(pt P, float pr) {C[nb].setTo(P); r[nb++]=pr; }
  void changeRadius() {r[b]+=(mouseX-pmouseX)/10.; if(r[b]<1) r[b]=1;  }
 
  void showBalls() { for (int i=0; i<nb; i++)  C[i].show(r[i]); };    

  boolean pickBall (pt M) {boolean picked = false; float dm=10000.0; 
    for (int i=0; i<nb; i++)  {
      float dd = M.disTo(C[i]); 
      if((dd<r[i]*1.1)&&(dd<dm)) {picked=true; dm=M.disTo(C[i]); b=i; };  }; return picked;}
 
  void showCaplets(int s) {for (int i=0; i<nb-1; i++)  showCone(C[i],r[i],C[i+1],r[i+1],s);}
 
   }  // end Cube class

void showCone(pt A, float a, pt B, float b, int s) {
  vec T = A.vecTo(B); float d=T.norm(); T.makeUnit();
  vec U = new vec (1,0,0); vec N = cross(U,T);
  if (N.norm()<0.1) {U.setTo(0,1,0);  N.setTo(cross(U,T)); };
  N.makeUnit(); U.setTo(cross(N,T));
  float dd = A.disTo(B);
  float x=(a-b)/dd;
  float y=sqrt((dd-a+b)*(dd+a-b))/dd;
  pt PP = A.make();   pt QQ = B.make();
  PP.addScaledVec(a*x,T);   QQ.addScaledVec(b*x,T);
  pt P = new pt(); pt Q = new pt();
  beginShape(QUAD_STRIP);
  for (float t=0; t<TWO_PI; t+=TWO_PI/s) {
    P.setTo(PP);   Q.setTo(QQ);
    P.addScaledVec(a*y*cos(t),N);   Q.addScaledVec(b*y*cos(t),N); 
    P.addScaledVec(a*y*sin(t),U);   Q.addScaledVec(b*y*sin(t),U); 
    P.vert(); Q.vert();
    };
  endShape();
  }

void showInbetweenBall(pt A, float a, float s, pt B, float b) {
  if ((min(a,b)<s)&&(s<max(a,b))) {pt P = A.make();  P.moveTowards( (s-a)/(b-a),B);  P.show(a+s*(b-a)); }; }