#include <stdio.h>
#include <stdlib.h>
#include <math.h>

float PI=3.141592654;

float mt = 174.3;
float mb = 4.7;
float v = 246.0;

float MZ = 91.188;
float MW = 79.825;
float MH = 125.;
float s2tW;


float lambda=3;

float M_ini=300., M_fin=4000.;
float dtM=100.;

float s1, c1, cR, cRb;


float width_z(float M, float vert, float m);
float width_w(float M, float vert, float m);
float width_h(float M, float vert, float m);


int main () {



FILE * out_Ts;

float M,k=1;


float  Ts_Wb=0, Ts_Zt=0, Ts_Ht=0 ;
float  width_Ts;

float  MTs; /* masse bare */
float  MTsp; /* masse fisiche */

s2tW=1-pow(MW/MZ,2);


out_Ts =fopen("Width_Ts.txt","w");


fprintf(out_Ts,"M_Ts \t WIDTH (lambda 3) \t WIDTH (lambda 2) \t WIDTH (lambda 1) \n ");


for(M=M_ini; M<M_fin; M=M+dtM) {


MTsp=M;

 
 Ts_Wb = width_w(MTsp, lambda , mb);
 Ts_Zt = width_z(MTsp, lambda, mt);
 Ts_Ht = width_h(MTsp, lambda, mt);

 
 width_Ts = Ts_Wb + Ts_Zt + Ts_Ht;
 
 Ts_Wb = Ts_Wb/width_Ts;
 Ts_Zt = Ts_Zt/width_Ts;
 Ts_Ht = Ts_Ht/width_Ts; 
 

fprintf(out_Ts," %f \t %f \t %f \t %f\n", MTsp, width_Ts, width_Ts*4/9, width_Ts/9);


                            }

     fclose(out_Ts);

	return 0  ;
 }




float width_w(float M, float vert, float m){

float y;

y=pow(vert,2)/(32.*PI)*M*((1+pow(m/M,2)-pow(MW/M,2))*(1+pow(m/M,2)+2*pow(MW/M,2))-4*pow(m/M,2))*sqrt(1-2*(pow(m/M,2)+pow(MW/M,2))+pow((m*m-MW*MW)/(M*M),2));
return y;
}


float width_z(float M, float vert, float m){

float y;

y=1./2.*pow(vert,2)/(32.*PI)*M*((1+pow(m/M,2)-pow(MZ/M,2))*(1+pow(m/M,2)+2*pow(MZ/M,2))-4*pow(m/M,2))*sqrt(1-2*(pow(m/M,2)+pow(MZ/M,2))+pow((m*m-MZ*MZ)/(M*M),2));
return y;
}


float width_h(float M, float vert, float m){

float y;

y=1./2.*pow(vert,2)/(32.*PI)*M*(1-pow(MH/M,2)+pow(m/M,2))*sqrt(pow(1+MH*MH/(M*M)-m*m/(M*M),2)-4*pow(MH/M,2));
return y;
}