x=7.5 //demi-largeur du scintillateur
y=14.25 //demi-longueur du scintillateur
z=20 //distance entre les 2pm
muon=0 //initialisation du nombre de muons détectés
thett=[] //initialisation 

for i=1:1000000
        //les 3 lignes suivantes me permettent dobtenir une loi de distribution pour mon angle zénithal en cos2
        u=grand(1,1,'unf',0,1)
        theta=acos(u^(0.25))
        thett=[thett,theta]
        //les deux lignes suivantes simulent les coordonnées d'arrivées d'un muon sur le scintillateur du haut (tirage uniforme) 
        a=grand(1,1,'unf',-x,x) 
        b=grand(1,1,'unf',-y,y)
        //tirage aléatoire uniforme de 0 ) pi sur l'angle longitudinale
        alpha=grand(1,1,'unf',0,2*%pi)
        r=z/cos(theta) //les contraintes imposent le paramètre r
        x2=r*sin(theta)*cos(alpha)+a+1 //position en x du muon sur le 2eme scintillateur
        y2=r*sin(theta)*sin(alpha)+b //position en y du muon
                     if abs(x2)<x & abs(y2)<y then muon=muon+1 // si le muon est dans l'intervalle définissant la position du 2eme scintillateur, j                     'incrémente ma variable muon
                     end
end

acc=(muon/size(thett,2))*100 // je fais le rapport muons détecté par les 2 scintillateurs sur le nombre total de muons qui a satisafait la condition de distribtion
