source: trunk/MagicSoft/Simulation/Corsika/Simone/simone.cxx@ 721

//
//
//
//
//
//
#include <stdlib.h>
#include <stdio.h>

#include "COREventHeader.hxx"
#include "CORParticle.hxx"
#include "CORStatfile.hxx"

int main(int argc, char **argv)
{
  char  path[100] = "" ; 

  char  cername[120] ; 
  char  datname[120] ; 
  
  ifstream   cerfile ; 
  
  COREventHeader Event ;   
  CORParticle    Photon ; 
  
  Int_t iPhotonInShower ; 
  
  Float_t lambda ; 

  cout << " ============================" << endl ; 
  cout << " SIMONE" << endl ; 
  cout << "                             " << endl ; 
  cout << " SImulated MONte carlo Events" << endl ; 
  cout << "                             " << endl ; 


  if ( argc <= 1 ) {
    cout << endl ; 
    cout << " INFO: You have to start the program with "<<endl <<endl ; 
    cout << "        -> simone DIRECTORY_WITH_CER_FILES" <<endl <<endl ; 
    cout << "        no SLASH at the end of the directory name"<<endl ; 
    cout << "        (example: -> simone /hd123/Protons   "<<endl ; 
    exit (-1) ; 
  }

  // copy the argument to filename 

  sprintf (cername , "%s", argv[1] ) ; 

    
  cout << cername  << endl ; 
      
  //
  //   try to open the files
  // 
    
  cerfile.open( cername );
    
  if ( cerfile.bad() ) { 
    cout << "Cannot open input file:  " << cername << endl ;
    exit (1) ; 
  }
    
  //   
  //    cout << " Read event " << endl ;
  //
  
  //   tricky stuff. We read in the RunHeader which has by chance the
  //   same sice than the EventHeader
  Event.read( cerfile ); 

  //    now we read in the event 
  
  Event.read( cerfile );
  
  Event.Print() ; 
  
  //
  //   loop over the particles (cerenkov photons) in 
  //   the file 
  //
  
  iPhotonInShower = 0 ; 
  
  while( ! (cerfile.eof() || cerfile.bad() )) {
    
    cout << "   testit " ; 
    //
    //   read in the particles
    //
    
    Photon.read ( cerfile ) ; 
    
    //
    //   only if the wavelength lambda is greater than 
    //   1.0 it is a real cerenkov photon 
    //
    lambda = Photon.get_wl() ; 

    Photon.print() ; 
    
    if ( lambda < 1.0 ) { 
      // here we got to the next event reading in all the ZERO 
      // from the file 
      while ( lambda < 1.) { 
        Photon.read ( cerfile ) ; 
        lambda = Photon.get_wl() ;      
      } 

      //      sprintf (" size of Photon: %d \n", (int) sizeof( CORParticle )) ; 
      fseek ( cerfile, -sizeof(CORParticle), SEEK_CUR(1) ); 

      cout << sizeof ( CORParticle ) ;
      cout << sizeof ( int ) ; 

      break ; 
      
      
      
      
    } 
    
    
    iPhotonInShower++ ; 
    
  }
  
  //
  //   close the file
  //
  
  cerfile.close();
  
  
  
}