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

//
//
//
//
//
//
#include <stdlib.h>
#include "TROOT.h"
#include "TApplication.h"

#include "TCanvas.h"
#include "TView.h"

#include "TSPHE.h"
#include "TBRIK.h"

#include "TNode.h"

#include "TFile.h"
#include "TNtuple.h"

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

#include "MSimone.hxx"

int main(int argc, char **argv)
{
  // initialise ROOT

  TROOT simple("simple", "SIMONE - statistic of MonteCarlo");

  TApplication theApp("App", &argc, argv);
  
  char  path[100] = "/dat1/Mmcs/Gamm/" ; 
  //    char  path[100] = "/hd02/Maggi/Data/prot_15/" ; 
  //    char  path[100] = "/data/mmcs/" ; 
  char  cername[120] ; 
  char  datname[120] ; 
  
  ifstream   cerfile ; 
  
  COREventHeader Event ;   
  CORParticle    Photon ; 
  
  MSimone    Infos ; 
    
  Int_t iPhotonInShower ; 
  
  Float_t lambda ; 
  
  cout << " SImulated MONte carlo Events" << endl ; 

  //
  //    now create a root-file for the ntuple output
  // 
  
  TFile *outfile = new TFile("simone.root","RECREATE");
  
  TNtuple *Ntup = new TNtuple("simone",
                              "Simone info of mmcs",
                              "fPartId:fEnergy:fTheta:fPhi:fCoreX:fCoreY:fCoreD:fNumCerPhot");
  
  
  for (int i_cer = 1; i_cer <= 100; i_cer++ ) {
    //
    //   info of progress
    //
    if (!( i_cer %10) )
      {
        cout << i_cer << endl ; 
      }

    //
    //   create the file names 
    //
    sprintf ( cername, "%s/cer%06d", path, i_cer ) ; 
    sprintf ( datname, "%s/dat%06d", path, i_cer ) ; 
    
    //      cout << cername  << endl ; 
    //      cout << datname  << endl ; 
    
    //
    //   try to open the files
    // 
    
    cerfile.open( cername );
    
    if ( cerfile.bad() ) { 
      cout << "Cannot open input file:  " << cername << endl ;
      continue ; 
    }
    
    //   
    //    cout << " Read event " << endl ;
    //
    
    Event.read( cerfile );
    Infos.Transfer ( &Event ) ; 
    
    //
    //   loop over the particles (cerenkov photons) in 
    //   the file 
    //
    
    iPhotonInShower = 0 ; 
    
    while( ! (cerfile.eof() || cerfile.bad() )) {
      
      //
      //   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() ; 
      
      if ( lambda < 1.0 ) 
        break ; 
      
      iPhotonInShower++ ; 
      
      //Photon.print() ; 
    }
    
    Infos.SetNumCerPhot( iPhotonInShower ) ; 
    
    //
    //   fill the ntuple
    //
    
    Infos.NtupFill ( Ntup ) ; 
    
    
    cerfile.close();
    
    
    
  } 
  
  
  //
  //   write all to file
  //
  outfile->Write() ; 

  
  //
  //
  //

  Ntup->Draw("fEnergy") ; 


  theApp.Run();
  return 0;                            
   
  
}