-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathEnergyLoss.java~
48 lines (40 loc) · 1.3 KB
/
EnergyLoss.java~
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
import java.io.*;
class EnergyLoss
{
static PrintWriter screen=new PrintWriter (System.out, true);
protected String material;
protected double Z, A, Ro, E, b, energyLoss, p, y, I;
protected final double mu= 106;
protected final double me= 0.511;
public double getEnergyL( double E)
{
p = Math.sqrt(E*E-mu*mu);
b = (p/(Math.sqrt((p*p)+(mu*mu))));
y = (1/(Math.sqrt(1-(b*b))));
I= (1.35e-5)*Z;
//energyLoss = 0.307*Ro*(Z/A)*(1/(b*b))*((Math.log((2*me*b*b*y*y)/I))-(b*b));
return 0.307*Ro*(Z/A)*(1/(b*b))*((Math.log((2*me*b*b*y*y)/I))-(b*b));
}
public EnergyLoss(String material, double z, double a, double ro, double e)
{
Z=z;
A=a;
Ro=ro;
E=e;
/*
double b = (p/(Math.sqrt((p*p)+(mu*mu))));
double y = (1/(Math.sqrt(1-(b*b))));
double I= (1.35e-5)*z;
*/
//double energyLoss = 0.307*ro*(z/a)*(1/(b*b))*((Math.log((2*me*b*b*y*y)/I))-(b*b));
double EL=getEnergyL(e);
screen.println(" Energy loose "+EL);
screen.println(" momentum "+p);
screen.println(" beta "+b);
screen.println(" gamma "+y);
}
public double getBeta()
{
return b;
}
}