Last change
on this file since 429 was
157,
checked in by campagne, 18 years ago
|
fix bug: 100km -> 130km factor (JEC)
|
File size:
727 bytes
|
Line | |
---|
1 | { |
---|
2 | |
---|
3 | ifstream fIn; |
---|
4 | fIn.open("spectre45minus.out"); |
---|
5 | if (!fIn) { |
---|
6 | std::cout << "Input file cannot be opened.\n"; |
---|
7 | return; |
---|
8 | } |
---|
9 | |
---|
10 | double bin,ene,numu,anumu,nue,anue,nutau,anutau = 0; |
---|
11 | double scale; // 100km -> 130km 1/L**2 law |
---|
12 | |
---|
13 | scale = (100.0/130.0); |
---|
14 | scale *= scale; |
---|
15 | |
---|
16 | |
---|
17 | while (!fIn.eof()) { |
---|
18 | fIn >> bin |
---|
19 | >> ene |
---|
20 | >> numu |
---|
21 | >> anumu |
---|
22 | >> nue |
---|
23 | >> anue; |
---|
24 | printf("%.2f %.4e %.4e %.4e %.4e %.4e %.4e\n",ene,nue*scale,numu*scale,nutau,anue*scale,anumu*scale,anutau); |
---|
25 | }//while |
---|
26 | |
---|
27 | while (bin!=501) { |
---|
28 | bin++; |
---|
29 | ene +=0.02; |
---|
30 | nue=0.; |
---|
31 | numu=0.; |
---|
32 | anue=0.; |
---|
33 | anumu=0.; |
---|
34 | printf("%.2f %.4e %.4e %.4e %.4e %.4e %.4e\n",ene,nue*scale,numu*scale,nutau,anue*scale,anumu*scale,anutau); |
---|
35 | } |
---|
36 | }//eo macro |
---|
Note: See
TracBrowser
for help on using the repository browser.