source: GLBFrejus/HEAD/data/4.5GeV/convert.C@ 689

{

  ifstream fIn;
  fIn.open("spectre45minus.out"); 
  if (!fIn) {
    std::cout << "Input file cannot be opened.\n";
    return;
  }
  
  double bin,ene,numu,anumu,nue,anue,nutau,anutau = 0;
  double scale; // 100km -> 130km 1/L**2 law

  scale = (100.0/130.0);
  scale *= scale;


  while (!fIn.eof()) {
    fIn >> bin 
        >> ene 
        >> numu
        >> anumu
        >> nue
        >> anue;
    printf("%.2f %.4e %.4e %.4e %.4e %.4e %.4e\n",ene,nue*scale,numu*scale,nutau,anue*scale,anumu*scale,anutau);
  }//while
  
  while (bin!=501) {
    bin++;
    ene +=0.02;
    nue=0.;
    numu=0.;
    anue=0.;
    anumu=0.;
    printf("%.2f %.4e %.4e %.4e %.4e %.4e %.4e\n",ene,nue*scale,numu*scale,nutau,anue*scale,anumu*scale,anutau);  
  }
}//eo macro