]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | function cumfit(x) |
2 | c fit for ALICE 94 RUN DATA | |
3 | c data p1,p2 /2.526, 1.104/ | |
4 | c94 data p1,p2 /2.300, 2.044/ | |
5 | c crystals for my geant | |
6 | c data p1,p2 /1.480, 1.600/ | |
7 | c 0rad | |
8 | data p1,p2 /1.700, 1.150/ | |
9 | c 0.5rad + | |
10 | c data p1,p2 /1.700, 0.800/ | |
11 | c 0.5rad - | |
12 | c data p1,p2 /5.500, 0.700/ | |
13 | c 0.3rad + | |
14 | c data p1,p2 /1.700, 0.900/ | |
15 | c 0.3rad - | |
16 | c data p1,p2 /4.00, 0.700/ | |
17 | c 0.2rad + | |
18 | c data p1,p2 /1.700, 1.000/ | |
19 | c 0.2rad - | |
20 | c data p1,p2 /3.50, 0.800/ | |
21 | c 0.1rad + | |
22 | c data p1,p2 /1.700, 1.100/ | |
23 | c 0.1rad - | |
24 | c data p1,p2 /2.50, 0.900/ | |
25 | c data angl /0./ | |
26 | data angl /-0.2/ | |
27 | ||
28 | aangl=abs(angl) | |
29 | if(angl.ge.0.) then | |
30 | pp1=p1 | |
31 | pp2=p2-0.7*aangl | |
32 | else | |
33 | pp1=p1+7.6*aangl | |
34 | pp2=p2-0.9*aangl | |
35 | endif | |
36 | x=abs(x) | |
37 | c lead glass for my geant | |
38 | c x=x*0.55 | |
39 | ||
40 | a=0.5*exp((pp1-sqrt(pp1**2+4.*pp2*x))/2./pp2) | |
41 | if(x.ge.0.) then | |
42 | cumfit=a | |
43 | else | |
44 | cumfit=1.-a | |
45 | endif | |
46 | end |