Last change
on this file since 483 was
157,
checked in by campagne, 18 years ago
|
fix bug: 100km -> 130km factor (JEC)
|
File size:
727 bytes
|
Rev | Line | |
---|
[153] | 1 | { |
---|
[157] | 2 | |
---|
| 3 | ifstream fIn; |
---|
| 4 | fIn.open("spectre45minus.out"); |
---|
| 5 | if (!fIn) { |
---|
| 6 | std::cout << "Input file cannot be opened.\n"; |
---|
| 7 | return; |
---|
| 8 | } |
---|
[153] | 9 | |
---|
[157] | 10 | double bin,ene,numu,anumu,nue,anue,nutau,anutau = 0; |
---|
| 11 | double scale; // 100km -> 130km 1/L**2 law |
---|
[153] | 12 | |
---|
[157] | 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); |
---|
[153] | 25 | }//while |
---|
| 26 | |
---|
[157] | 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); |
---|
[153] | 35 | } |
---|
[157] | 36 | }//eo macro |
---|
Note: See
TracBrowser
for help on using the repository browser.