#include "diabolo.h"
#include "bolo.h"
#include "arcunit.h"

//*********  coefficients  pour  les  mesures  bolo    ***************************************************
//													//
//	-1-	loi de reponse thermique des bolos  avec R en ohms et T en Kelvin  coef 0,1,2		//
//													//
// 	T =  coef2  *   (   ln (  R / coef1) **  ( -1 / coef0 )						//
//													//
//	-2-	fuite thermique du bolo   coef 3,4							//
//													//
//	Ptot =  coef3  *  ( Tb ** coef4    -    Tcryo  **  coef4 )					//
//													//
//	-3-	calcul empirique de Pciel et de  tau	  coef  5,6					//
//													//
//		Pciel	=	V * I -  coef5			coef5= I * Ai (tables xavier)		//
//		tau	=	- ln ( 1 +  Pciel / coef6 )	coef6= I * Bi (tables xavier)		//
//													//
//******************************************************************************************************//

#define	c(i)		(1e-4*(double)parametr.nom_coef[parametr.bolo[fen-1].numero_nom_coef].coef[i])	


void mesures_bolo(int fen,int carr,int tri,int flag)
{
double I,V,R,T,Pelec,Tcryo,Ptot,Pciel,tau;
double a;
def_gains;

if(parametr.bolo[fen-1].bolo_diviseur)		//  BEBO normale ou MLPA
	{
	I = (double)tri / pt_micA(fen);						//	I en µA
	V= (double) carr / pt_micV(fen);					//	V en µVolts
	if(I>0.0000001)		R=V/I;	else	R=0;				//	R  en ‡
	if(R>=1e6) ecritD(fen,b_res,"%7.3fM‡ %6.3fnA %7.3fmV",R*1e-6,I*1e3,V*1e-3);// R  en M‡ I en nA et V en mV
	else	 ecritD(fen,b_res,"%6.2fK‡ %6.3fnA %7.3fmV",R*1e-3,I*1e3,V*1e-3);// R  en K‡   I en nA et V en mV
	V=V-xbolbrut(fen-1);	//  corrigée du déséquilibre  en  µV
	}
else						//  pour carte BEBO modifiée (mesure temperature)
	{
	I = 1e-3 * (double)carr * 2441. / parametr.bolo[fen-1].bolo_capa;	//	I en µA
	V=xbolbrut(fen-1);							//	V en µVolts
	if(I>0.0000001)		R=V/I;	else	R=0;				//	R  en ‡
	if(R>=1e3) ecritD(fen,b_res,"%7.3fK‡ %6.3fnA %7.3fmV",R*1e-3,I*1e3,V*1e-3);// R  en K‡ I en nA et V en mV
	else	 ecritD(fen,b_res,"%6.2f‡ %6.3fnA %7.3fmV",R,I*1e3,V*1e-3);	// R  en ‡   I en nA et V en mV
	}


if(!flag) return;

if(!fenetre(fenetre_mesures_bolo)) nouveauT(fenetre_mesures_bolo,mesures_bolo_id,"mesures bolos");

	//****************  tension corrigée du déséquilibre     *******************

if(I<0.000001) return;
if(V<-1000000) return;
if(V>1000000) return;


Pelec=V*I;		// pW

//	-1-	loi de reponse thermique des bolos  avec R en ohms et T en Kelvin  coef 0,1,2		//
//													//
// 	T =  coef2  *   (   ln (  R / coef1) **  ( -1 / coef0 )						//

a=1;	if( (R>0) && (c(1) >0.01) )		a= log ( R / c(1) );
T=0;	if( (a>0) && (c(0)>0.01) )		T= c(2) * pow( a , -1 / c(0) );

//	-2-	fuite thermique du bolo   coef 3,4							//
//													//
//	Ptot =  coef3  *  ( (10*Tb) ** coef4    -    (10*Tcryo)  **  coef4 )				//
/* Desormais 30/06/1999, Tcryo est calculee avec la germanium */
//Tcryo=0.1;
Tcryo= gg->temperature_cryo ;
Ptot=0;	if( (c(4)>0.01) && (T>0.01) )	Ptot = c(3) * ( pow(10.*T,c(4)) - pow(10.*Tcryo,c(4)) );

//	-3-	calcul empirique de Pciel et de  tau	  coef  5,6					//
//													//
//		Pciel	=	coef5  - V * I   		coef5= I * Ai (tables xavier)		//
//		tau	=	- ln ( 1 +  Pciel / coef6 )	coef6= I * Bi (tables xavier)		//

//Pciel = 0;  if( c(5)>Pelec)	Pciel = c(5) - Pelec;
Pciel = 0;  if( Ptot>Pelec)	Pciel = Ptot - Pelec;  /* Pour avoir la puissance de rayonnement absorbé */
//a=1;	if( c(6) >0.1 )		a =  1 - Pciel / c(6) ;
//tau=0;	if( a>0 )		tau =  - log ( a );
tau=Tcryo ;  /* FXD pour avoir Tcryo 16/07/99 */


ecritT(fenetre_mesures_bolo,fin_f,"%d %s :",fen,parametr.nom_coef[parametr.bolo[fen-1].numero_nom_coef].bolo_nom);

//ecritT(fenetre_mesures_bolo,fin_f,"%s (b%d) :"/* bebo=%d  num=%d  gain_pa=%d   coef0=%g  coef1=%g  coef2=%g */,
//		parametr.bolo_nom[fen-1],fen,parametr.bolo[fen-1].bolo_bebo,parametr.bolo[fen-1].bolo_num,parametr.bolo_gain[fen-1]
//		,parametr.coef[fen-1][0],parametr.coef[fen-1][1],parametr.coef[fen-1][2]);


ecritT(fenetre_mesures_bolo,fin_f," R=%7.3fM‡  I=%6.3fnA  V=%7.3fmV   Pelec=%6.2fpW",R*1e-6,I*1e3,V*1e-3,Pelec);

if(T>0.01) ecritT(fenetre_mesures_bolo,fin_f,"  T=%6.1fmK  ", T*1e3);

if( Ptot>0.01) ecritT(fenetre_mesures_bolo,fin_f," Ptot=%6.2fpW   Pciel=%6.2fpW   Tcryo=%6.1f  ",Ptot,Pciel,tau*1e3);


//ecritT(fenetre_mesures_bolo,fin_f,": C=%g  ptmic=%g\n",capa(fen),pt_micA(fen));

ecritT(fenetre_mesures_bolo,fin_f,"\n");
}