//*********************************************** //fonction créant une bille et deux rayons lasers function OT(x,y,z) { //newPhot(x-0.5,y-0.5,z,1,1,0); //newPhot(x+0.5,y-0.5,z,-1,1,0); newBead(x,y,z,0.06,1); } //*********************************************** //*************************************************** print(commands); setScalablepG(1); setMovable(0); setMovablepG(1); //Debut du programme principal setViewXYpG(0,-1); // Vue centrée sur la pince setViewX(0.5); setViewY(-0.5); //nouveau polymère newFJC; //Double pince optique OT(-0.6,0,0); OT(0.6,0,0); //Attache le polymère aux billes attachG(0); attachD(1); T=2; setT(T); // Graphique force extension Npt=1000; setWindow(0,0); newGraph(Npt); setBounds(0,-2,16,3.8); g = lastGraph; colorGraph(0,0,0); mean_window=10; setAxisX(1,0,2,0,0.1); setAxisY(1,0,0.5,0,0.1); setViewXYpG(0,-1); addLabel(1,1.5,'log(Force (~pN))'); addLabel(10,-0.5,'Length (~µm)'); tension=0; N=0; // pour le dragging nhold=0; hasdragged=0; distdragging=0.07; dd2=distdragging*distdragging; sync(0<1) { N++; if ((N%mean_window)==0) {addPtGraph((getBeadPosX(1)+0.6)*10,log(tension*50/20)/log(10));tension=0;} tension=tension+getTensionX; setBeadPosY(0,0); setBeadPosY(1,0); setBeadPosX(0,-0.6); // Reponse au dragging if (isMouseDragged==1) { if (hasdragged==1) { // continue to drag point kdrag=200; setBeadForce(nhold,kdrag*(getMouseX-getBeadPosX(nhold)),kdrag*(-getMouseY-getBeadPosY(nhold)),0); //setPtY(nhold,-getMouseY); } else { // look for a point to drag found=0;i=0;gN=getNBille;mx=getMouseX; my=-getMouseY; while ((found==0)&(i<gN)) { px=getBeadPosX(i);py=getBeadPosY(i); if (((mx-px)*(mx-px)+(my-py)*(my-py))<dd2) { found=1; nhold=i;} i=i+1;} if (found==1) { hasdragged=1; } } } else { if (hasdragged==1) { // free dragged point hasdragged=0; setBeadForce(nhold,0,0,0); } } if (kP==1) { a=rK; if(a=='p') { // increase NA T=T*1.01; setT(T); } if(a=='m') { // decrease NA T=T/1.01; setT(T); } } }