Changeset 778 for trunk/MagicSoft/Simulation/Corsika

Timestamp:

05/03/01 18:39:16 (24 years ago)

Author:

domingo

Message:

Version of GuiSimone for a "single file per run" output of Mmcs.
The "previous" event command is not operational at present. May need some
further thinking...

Location:

trunk/MagicSoft/Simulation/Corsika/GuiSimone

Files:

• MainFrameSimone.cxx (modified) (26 diffs)
• MainFrameSimone.hxx (modified) (1 diff)

trunk/MagicSoft/Simulation/Corsika/GuiSimone/MainFrameSimone.cxx

@@ -1 +1 @@
 //
-//
-//
-#define  START_DIRECTORY  "."
+////////////////////////////////////////////////////////////////////////////
+//
+// GuiSimone 
+//
+// Graphic User Interface for Simone (SIMultaed ONe Event) 
+//
+// Modification May 2001 by C.Domingo over the existing program:
+//   The header file MainFrameSimone.hxx has been changed: the ifsteam fCerFile
+//     and the string variable fCerFileName are now global  
+//   Some procedures in MainFrameSimone have been adapted to the "one
+//     file per run" Mmcs output
+//
+// GuiSimone reads Mmcs output files (one file per run version) and 
+// displays several histograms for each run and/or each event. Selection
+// of files is made through a window (frame) and graphics are displayed
+// in other windows (subframes)
+//
+// The ROOT package is needed
+//
+////////////////////////////////////////////////////////////////////////////
+//
+//
+
+#define  START_DIRECTORY  "/dat1"
 #define S_ISDIR(m) (((m)&(0170000)) == (0040000))  
 #define PI 3.1415926
@@ -10 +31 @@
 
 #include "MainFrameSimone.hxx"
-
 
 char NtupName[5][20] = {"WaveLen", "PosY:PosX", "Height", "Time", "DireY:DireX" } ;  
@@ -57 +77 @@
   fMenuBar = new TGMenuBar ( this, 1, 1, kHorizontalFrame ) ; 
   fMenuBar->AddPopup("File", fFileMenu, fLayMenuItem ) ; 
-  
+   
   AddFrame(fMenuBar, fLayMenuBar ) ; 
 
@@ -114 +134 @@
   
   fSelectButton = new TGTextButton(fTabF1b, "Select", M_BUTTON_SELECT );
-  fSelectButton->SetToolTipText("To select a path with Mmcs data files!");
+  fSelectButton->SetToolTipText("To select a Cerenkov data output file from Mmcs!");
   fSelectButton->SetState(kButtonDisabled); 
   fSelectButton->Associate(this) ; 
@@ -493 +513 @@
 // ======================================================================
 // ======================================================================
+//
+////////////////////////////////////////////////////////////////////////////
+//
+// Modification by C.Domingo May 2001
+//
+// Procedure CerFileExist NOT used at present
+//
+////////////////////////////////////////////////////////////////////////////
+//
 
 Bool_t MainFrameSimone::CerFileExist(Char_t *dir)
@@ -583 +612 @@
 // ======================================================================
 // ======================================================================        
+//
+////////////////////////////////////////////////////////////////////////////
+//
+// Modification by C.Domingo May 2001
+//
+// The fNtupFile is now the cerenkov FILE (not the directory where the
+// Cerenkov files are stored) which contains the information of a 
+// complete RUN
+//
+////////////////////////////////////////////////////////////////////////////
+//
 
 Bool_t MainFrameSimone::CheckNtup()
@@ -593 +633 @@
     new TGMsgBox(fClient->GetRoot(), this,
                  "Warning", 
-                 "No Input Path specified\n Do this before!!", 
+                 "No Cerenkov Input File specified\n Do this before!!", 
                  kMBIconExclamation, buttons, &retval);
     
@@ -600 +640 @@
   
 
-  if ( ! strcmp ( fInputPath, fNtupPath ) ) {
+  if ( ! strcmp ( fCerFileName, fNtupFile ) ) {
     return (kTRUE) ; 
   }
@@ -615 +655 @@
 // ======================================================================
 // ======================================================================
-
+//
+////////////////////////////////////////////////////////////////////////////
+//
+// Modification by C.Domingo May 2001
+//
+// Major changes to procedure ReadDataToNtup
+//   One simone.root file per each RUN file is generated now
+//   One cerenkov file is read instead of all cerenkov files in a given directory
+//   All events in the cerenkov file selected are read, and all the photons of these
+//     events are skipped to proceed to next event
+//   
+//
+////////////////////////////////////////////////////////////////////////////
+//
 Bool_t MainFrameSimone::ReadDataToNtup()
 {
   FILE  *unit ; 
   Char_t  rootfile[256] ; 
-  Char_t  cername[256] ; 
   Char_t  datname[256] ;
   
-  ifstream   cerfile ; 
-  
-  COREventHeader Event ;  
+  Float_t lambda ;
+  Float_t energyofevent ;
+  Int_t iPhotonInShower ;
+  
+  COREventHeader Event ;
+  CORParticle Photon ;
+
+  Int_t buttons = 4, retval ;
   
   fNtup->Reset() ; 
 
-  sprintf (rootfile, "%s/simone.root", fInputPath ) ; 
-
-  //
-  //   check if there is a file simone.root in the directory fInputPath
+  sprintf (rootfile, "%s/simone%s.root", fInputPath, fRunSuffix ) ; 
+
+  //
+  //   check if there is a file simone######.root in the directory fInputPath
   //
   //   if the file is there open it an read the ntuple in, 
-  //   else open all cerenkov files and fill the histograms and store 
-  //   the ntuple as simone.root in the directory fInputPath
+  //   else open the selected cerenkov file and fill the histograms and store 
+  //   the ntuple as simone######.root in the directory fInputPath
   //
 
@@ -663 +720 @@
 
     //delete hfile ;
-
-    
 
     fNumEvents = (Int_t) fNtup->GetEntries()  ;
@@ -683 +738 @@
 
     
-    for (int i_cer = 1; i_cer <= 100000; i_cer++ ) {
-
-      //   inform about progress
-      
-      if (!( i_cer %1000) )
-        {
-          cout << i_cer << endl ; 
-          sprintf ( datname, "%d", fNumEvents ) ; 
+    fNumEvents = 0 ;
+      
+    fCerFile.open( fCerFileName );
+          
+    if ( fCerFile.bad() ) {         
+      new TGMsgBox(fClient->GetRoot(), this,
+                   "Warning", 
+                   "Could not open Cerenkov file\n file = bad (sataus)", 
+                   kMBIconExclamation, buttons, &retval);
+      return (kFALSE) ;
+    }
+      
+    // Dirty trick to read run header, as by chance it has the same structure 
+    // of an event header
+      
+    Event.read(fCerFile) ;
+      
+    // Now reading the events truly
+      
+    while ( ! (fCerFile.eof() || fCerFile.bad() ) ) { 
+      
+      Event.read( fCerFile );
+      energyofevent = Event.get_energy() ;
+      if (energyofevent == 0) break ;
+      fNumEvents++ ;
+
+      //
+      // Inform about progress
+      //
+
+      if (!( fNumEvents %1000) ) {
+        sprintf ( datname, "%d", fNumEvents ) ; 
           
-          fDirNum->RemoveEntry(1) ; 
-          fDirNum->AddEntry(datname,1) ;            
-          fDirNum->MapSubwindows() ; 
-          fDirNum->Layout() ;  
-        }
-      
-      
-      //   create the file names 
-      
-      sprintf ( cername, "%s/cer%06d", fInputPath, i_cer ) ; 
-      sprintf ( datname, "%s/dat%06d", fInputPath, i_cer ) ; 
-      
-      //   try to open the files
-      
-      cerfile.open( cername );
-    
-      if ( cerfile.bad() ) { 
-
-        fNumEvents = i_cer -1 ;
-        
-        sprintf ( datname, "%d", fNumEvents ) ; 
-        
-        fDirNum->RemoveEntry(1) ; 
-        fDirNum->AddEntry(datname,1) ;      
-        fDirNum->MapSubwindows() ; 
-        fDirNum->Layout() ;  
-        
-        break ; 
-      }
-      
-      Event.read( cerfile );
+        fDirNum->RemoveEntry(1) ; 
+        fDirNum->AddEntry(datname,1) ;      
+        fDirNum->MapSubwindows() ; 
+        fDirNum->Layout() ; 
+      } 
 
       //
@@ -729 +782 @@
       
       fNtup->Fill( Event.get_primary() ,
-                   Event.get_energy() , 
-                   Event.get_coreX(), 
-                   Event.get_coreY(), 
-                   Event.get_theta(),
-                   Event.get_phi(), 
-                   Event.get_firstInt() ) ;  
-
- 
-      cerfile.close();
-    
+                   Event.get_energy() , 
+                   Event.get_coreX(), 
+                   Event.get_coreY(), 
+                   Event.get_theta(),
+                   Event.get_phi(), 
+                   Event.get_firstInt() ) ;  
+
+      //
+      //   loop over the particles (cerenkov photons) in 
+      //   the file 
+      //
+    
+      iPhotonInShower = 0 ; 
+    
+      while( ! (fCerFile.eof() || fCerFile.bad() )) {
+
+        //
+        //   read in the particles
+        //
+      
+        Photon.read ( fCerFile ) ; 
+      
+        //
+        //   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. && ! (fCerFile.eof() || fCerFile.bad() ) ) { 
+          
+            Photon.read ( fCerFile ) ; 
+            lambda = Photon.get_wl() ;  
+            //    cout << lambda << endl ; 
+          } 
+        
+          //    go back in the file, because when the zeroes have finished
+          //    the following record has already been read as it were
+          //    a cerenkov foton. The pointer should be repositioned to
+          //    read the event end record
+        
+          fCerFile.seekg( -28, ios::cur  ) ; 
+        
+          //  read in the event end record. Again, it has the same structure (by
+          //  chance) thant the event header
+
+          Event.read( fCerFile );
+        
+          break ; 
+        } 
+
+        iPhotonInShower++ ; 
+      
+      }
+    
+      // cout << " Event " << fNumEvents << ": Photons " << iPhotonInShower << endl ;
+
     }
 
+    fCerFile.close();
+      
+    //
+    //  Refresh information on the screen window
+    //
+
+    sprintf ( datname, "%d", fNumEvents ) ; 
+          
+    fDirNum->RemoveEntry(1) ; 
+    fDirNum->AddEntry(datname,1) ;          
+    fDirNum->MapSubwindows() ; 
+    fDirNum->Layout() ;     
+
     //
     //  write Ntuple to the file
     //
+  
     fNtup->Write() ; 
 
@@ -749 +867 @@
     //  write the root file
     //
+  
     hfile->Write();
     
@@ -754 +873 @@
   }
 
-
   //
   //    set some value of this class
   //
-  sprintf ( fNtupPath, "%s", fInputPath) ;
+
+  sprintf ( fNtupFile, "%s", fCerFileName) ;
   
   return (kTRUE) ; 
@@ -766 +885 @@
 // ======================================================================
 // ======================================================================
+//
+////////////////////////////////////////////////////////////////////////////
+//
+// Modification by C.Domingo May 2001
+//
+// Major changes to procedure ReadInEvt
+//   Reads the photons from inside the selected cerenkov file (which contains 
+//     generated for a given run) instead of opening one file per event
+//   Checks if it must read the first event of a file, in such a case the 
+//     file is open
+//   
+////////////////////////////////////////////////////////////////////////////
+//
 
 Bool_t MainFrameSimone::ReadinEvt()
@@ -773 +905 @@
   Char_t  wort[256] ; 
   FILE  *unit ; 
-  Char_t  filename[256] ; 
-
-  ifstream   cerfile ; 
-  
+
   COREventHeader Event ;   
   CORParticle    Photon ; 
@@ -791 +920 @@
     new TGMsgBox(fClient->GetRoot(), this,
                  "Warning", 
-                 "No Input Path specified\n Do this before!!", 
+                 "No Cerenkov Input File specified\n Do this before!!", 
                  kMBIconExclamation, buttons, &retval);
     
@@ -802 +931 @@
   }
 
-  sprintf ( filename, "%s/cer%06d", fInputPath, fEvtNb ) ;
-
   fNtupEvt->Reset() ; 
   
   //
-  //   open the file
-  // 
-  
-  cerfile.open( filename );
-    
-  if ( cerfile.bad() ) { 
-
+  //   if first event, open the file and read run header
+  //  
+  
+  if ( fEvtNb == 1) {
+    fCerFile.open( fCerFileName );    
+    if ( fCerFile.bad() ) { 
+
+      new TGMsgBox(fClient->GetRoot(), this,
+                   "ERROR", 
+                   "Can't open the cer file!!", 
+                   kMBIconExclamation, buttons, &retval);
+    
+      return (kFALSE) ; 
+    }
+
+    // Dirty trick to read run header (see comment in RaedDataToNtup for
+    // further details)
+      
+    Event.read(fCerFile) ;
+      
+  }
+  
+  //
+  //  For all events, read in event header
+  //
+      
+  Event.read( fCerFile );
+  
+  //
+  //  Check for RUN end (an ENDRUN record will return a 0 value of the energy)
+  //
+
+  if (Event.get_energy() == 0) {
     new TGMsgBox(fClient->GetRoot(), this,
-                 "ERROR", 
-                 "Can't open the cer file!!", 
-                 kMBIconExclamation, buttons, &retval);
-    
-    return (kFALSE) ; 
-  }
-
-  //
-  //   read in the header
-  //
-      
-  Event.read( cerfile );
-  
+                "Information", 
+                "End of file detected",
+                kMBIconAsterisk, buttons, &retval);
+    fCerFile.close () ;
+    return (kFALSE) ;
+  }
 
   //
@@ -860 +1006 @@
   //   loop over the particles
   //
-  
     
   iPhotonInShower = 0 ; 
-    
-  while( ! (cerfile.eof() || cerfile.bad() )) {
-      
+      
+  while( ! (fCerFile.eof() || fCerFile.bad() )) {
+
     //
     //   read in the particles
     //
       
-    Photon.read ( cerfile ) ; 
-          //
+    Photon.read ( fCerFile ) ; 
+      
+    //
     //   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++ ; 
-    
-    fNtupEvt->Fill(Photon.get_wl(), 
-                   Photon.get_x(), 
-                   Photon.get_y(), 
-                   Photon.get_h(), 
-                   Photon.get_t(), 
-                   Photon.get_u(), 
-                   Photon.get_v()  ) ; 
-
+ 
+    if ( lambda >= 1.0 ) {     
+      iPhotonInShower++ ; 
+      fNtupEvt->Fill(Photon.get_wl(), 
+                     Photon.get_x(), 
+                     Photon.get_y(), 
+                     Photon.get_h(), 
+                     Photon.get_t(), 
+                     Photon.get_u(), 
+                     Photon.get_v()  ) ;     
     }
-  
-  
+
+    // here we got to the next event reading in all the ZERO 
+    // from the file 
+
+    else {
+      while ( lambda < 1. && ! (fCerFile.eof() || fCerFile.bad() ) ) {     
+        Photon.read ( fCerFile ) ; 
+        lambda = Photon.get_wl() ;      
+      } 
+        
+      //    go back in the file, because when the zeroes have finished
+      //    the following record has already been read as it were
+      //    a cerenkov foton. The pointer should be repositioned to
+      //    read the event end record
+
+      fCerFile.seekg( -28, ios::cur  ) ; 
+        
+      //  read in the event end record. Again, it has the same structure (by
+      //  chance) thant the event header
+
+      Event.read( fCerFile );
+      break ;
+    }
+
+  }
+    
   //
   //   close the file
   //
-    
-  cerfile.close();
- 
+
   sprintf( wort," %10d", iPhotonInShower ) ; 
-  fInfoNbPhot->SetText( wort ) ; 
-
+  fInfoNbPhot->SetText( wort ) ;      
+
+  //  fCerFile.close();
+  
+  //  sprintf( wort," %10d", iPhotonInShower ) ; 
+  //  fInfoNbPhot->SetText( wort ) ; 
+  
   return (kTRUE) ; 
-
 }
+
+
 // ======================================================================
 // ======================================================================
+//
+////////////////////////////////////////////////////////////////////////////
+//
+// Modification by C.Domingo May 2001
+//
+// Minor changes to procedure ProcessMessage
+//   case M_BUTTON_SELECT takes into account that we no longer select 
+//     directories, but files.
+//   case M_BUTTON_PREV functionality has (temporarily) been removed
+//   case kCT_ITEMDBLCLICK takes into account that we deal with are going to
+//     select cerenkov files, not directories containing them. Test for detection
+//     of possible cerenkov files is added.
+//   
+////////////////////////////////////////////////////////////////////////////
+//
 
 Bool_t MainFrameSimone::ProcessMessage(Long_t msg, Long_t parm1, Long_t)
@@ -916 +1100 @@
   Char_t  testfile[100] ;
   Char_t  extens[5] ;
+  Char_t  idcer[5] ;
   Char_t  command[110] ;
  
@@ -985 +1170 @@
           //      
           item = (TGFileItem *) fFileCont->GetNextSelected(&np) ; 
-
-          if ( S_ISDIR(item->GetType()) ) {   // file is directory
-            
-            sprintf (fInputPath, "%s/%s",
-                     gSystem->WorkingDirectory(), 
-                     item->GetItemName()->GetString()) ;
-
-            fDir2->RemoveEntry(1) ; 
-            fDir2->AddEntry(fInputPath,1) ;         
-            fDir2->MapSubwindows() ; 
-            fDir2->Layout() ;   
+           
+          // cout << "fCerFileName =" << fCerFileName ;
+ 
+          //
+          // Check if fCerFile is already declared and open
+          // If so, close it before selecting next file
+          //
+
+          if (fCerFileName != "\0" &&  fCerFile ) 
+              fCerFile.close();
+
+          //
+          //  Compose the InputPath name and the Cerenkov filename
+          //
+
+          sprintf (fInputPath, "%s",
+                   gSystem->WorkingDirectory() ) ;
+
+          sprintf (fCerFileName, "%s/%s" ,
+                   gSystem->WorkingDirectory(), 
+                   item->GetItemName()->GetString()) ;           
+            
+          //cout << fInputPath << endl;
+          //cout << fCerFileName << endl ;
+
+          // 
+          //  Refresh values in the frames
+          //         
+
+          fDir2->RemoveEntry(1) ;     
+          fDir2->AddEntry(fCerFileName,1) ;
+          fDir2->MapSubwindows() ; 
+          fDir2->Layout() ;   
               
-            fDir3->RemoveEntry(1) ; 
-            fDir3->AddEntry(fInputPath,1) ;         
-            fDir3->MapSubwindows() ; 
-            fDir3->Layout() ;                 
-
-            fEvtNb = 0 ; 
-            sprintf (wort2, "%d", fEvtNb ) ; 
-            fDirEvtNb->RemoveEntry(1) ;
-            fDirEvtNb->AddEntry( wort2, 1 ) ; 
-            fDirEvtNb->MapSubwindows() ; 
-            fDirEvtNb->Layout() ; 
-
-          }
+          fDir3->RemoveEntry(1) ;     
+          fDir3->AddEntry(fCerFileName,1) ;         
+          fDir3->MapSubwindows() ; 
+          fDir3->Layout() ;           
+
+          fEvtNb = 0 ; 
+          sprintf (wort2, "%d", fEvtNb ) ; 
+          fDirEvtNb->RemoveEntry(1) ;
+          fDirEvtNb->AddEntry( wort2, 1 ) ; 
+          fDirEvtNb->MapSubwindows() ; 
+          fDirEvtNb->Layout() ; 
+
         }
         break;
@@ -1044 +1250 @@
 
 
-
       case M_BUTTON_PREV:
-        fEvtNb-- ; 
-        
+//      fEvtNb-- ; 
+//      
         if (fEvtNb < 1 ) {
           fEvtNb++ ; 
           
         }
-          
-        sprintf (wort2, "%d", fEvtNb ) ; 
-        fDirEvtNb->RemoveEntry(1) ;
-        fDirEvtNb->AddEntry( wort2, 1 ) ; 
-        fDirEvtNb->MapSubwindows() ; 
-        fDirEvtNb->Layout() ; 
-
-        ReadinEvt() ; 
-
-        DrawEvtParam() ; 
-
+//        
+//      sprintf (wort2, "%d", fEvtNb ) ; 
+//      fDirEvtNb->RemoveEntry(1) ;
+//      fDirEvtNb->AddEntry( wort2, 1 ) ; 
+//      fDirEvtNb->MapSubwindows() ; 
+//      fDirEvtNb->Layout() ; 
+//
+//      ReadinEvt() ; 
+//
+//      DrawEvtParam() ; 
+//
+//
+//      Do nothing !!
+//
+        new TGMsgBox(fClient->GetRoot(), this,
+                     "Warning", 
+                     "This button is not operational at present!!", 
+                     kMBIconExclamation, buttons, &retval);
+ 
         break;  
-
-
+//
+//
       case M_BUTTON_NEXT:
         fEvtNb++ ; 
@@ -1080 +1293 @@
 
         break;  
-
-
 
       case M_CTRL_EXIT:
@@ -1153 +1364 @@
             sprintf (wort, "%s", item->GetItemName()->GetString()) ;
             
+            sprintf (fInputPath,"\0") ; 
+            fSelectButton->SetState(kButtonDisabled) ;
+            fDir2->RemoveEntry(1) ; 
+            fDir2->AddEntry(fInputPath,1) ;
+            fDir2->MapSubwindows() ; 
+            fDir2->Layout() ;
+              
+            fDir3->RemoveEntry(1) ; 
+            fDir3->AddEntry(fInputPath,1) ;
+            fDir3->MapSubwindows() ; 
+            fDir3->Layout() ;
+
+            sprintf (wort2, "\0" ) ; 
+            fDirEvtNb->RemoveEntry(1) ;
+            fDirEvtNb->AddEntry( wort2, 1 ) ; 
+            fDirEvtNb->MapSubwindows() ; 
+            fDirEvtNb->Layout() ; 
+ 
+            //
+            //  goto directory 
             //
-            //  check if in this directory the file cer000001 exits
-            //  if it exists, do not open the stuff
-            //
-
-            if ( CerFileExist( wort ) ) {
-
-              fSelectButton->SetState(kButtonUp) ; 
+            fFileCont->ChangeDirectory ( wort ) ; 
+            gSystem->ChangeDirectory( wort ) ; 
               
-              new TGMsgBox(fClient->GetRoot(), this,
-                           "Information", 
-                           "This is a directory containing cer files\n click on select to chose it!",
-                           kMBIconAsterisk, buttons, &retval);
-
-            }  
-            else {
-              sprintf (fInputPath,"\0") ; 
-              fSelectButton->SetState(kButtonDisabled) ;
-              fDir2->RemoveEntry(1) ; 
-              fDir2->AddEntry(fInputPath,1) ;
-              fDir2->MapSubwindows() ; 
-              fDir2->Layout() ;
-              
-              fDir3->RemoveEntry(1) ; 
-              fDir3->AddEntry(fInputPath,1) ;
-              fDir3->MapSubwindows() ; 
-              fDir3->Layout() ;
-
-              sprintf (wort2, "\0" ) ; 
-              fDirEvtNb->RemoveEntry(1) ;
-              fDirEvtNb->AddEntry( wort2, 1 ) ; 
-              fDirEvtNb->MapSubwindows() ; 
-              fDirEvtNb->Layout() ; 
- 
-              //
-              //  goto directory 
-              //
-              fFileCont->ChangeDirectory ( wort ) ; 
-              gSystem->ChangeDirectory( wort ) ; 
-              
-              sprintf (wort, "%s", gSystem->WorkingDirectory() ) ; 
-              fDir->RemoveEntry(1) ; 
-              fDir->AddEntry( wort, 1 ) ; 
-              fDir->MapSubwindows() ; 
-              fDir->Layout() ; 
-            }
+            sprintf (wort, "%s", gSystem->WorkingDirectory() ) ; 
+            fDir->RemoveEntry(1) ; 
+            fDir->AddEntry( wort, 1 ) ; 
+            fDir->MapSubwindows() ; 
+            fDir->Layout() ; 
             
           }
@@ -1204 +1398 @@
             sprintf (testfile, "%s", item->GetItemName()->GetString() ) ; 
             
-            //
-            //    determine the file type by extensions
-            //
+            //
+            //   Is it a Cerenkov file?
+            //
+            
+            for (Int_t i = 0 ; i<5; i++) 
+                idcer[i] = '\0' ;
+            for (Int_t i = 0 ; i<10; i++)
+                fRunSuffix[i]= '\0' ;
+
+            strncpy (idcer, &testfile[0], 3 ) ;
+            if ( ! strcmp (idcer , "cer" ) ) {
+              for ( Int_t i=3; 15; i++) {
+                if ( testfile[i] == '\0') {
+                  strncpy (fRunSuffix, &testfile[3], i-3 ) ; 
+                  break ; 
+                }
+              }
+
+              // cout << "fRunSuffix = "<< fRunSuffix << endl ;
+
+              fSelectButton->SetState(kButtonUp) ; 
+              
+              new TGMsgBox(fClient->GetRoot(), this,
+                           "Information", 
+                           "This file probably contains CERENKOV events\n click on the SELECT button to chose it!",
+                           kMBIconAsterisk, buttons, &retval);                                    
+            }
+            
+            else {
+
+              fSelectButton->SetState(kButtonDisabled) ;
+              
+              //
+              //    determine the file type by extensions
+              //
             
-            for (Int_t i = 0 ; i<5; i++) 
-              extens[i] = '\0' ; 
+              for (Int_t i = 0 ; i<5; i++) 
+                extens[i] = '\0' ; 
             
-            for ( Int_t i=0; testfile[i] != '\0'; i++) {
-              if ( testfile[i] == '.') {
-                strncpy (extens, &testfile[i+1], 4 ) ; 
-                break ; 
+              for ( Int_t i=0; testfile[i] != '\0'; i++) {
+                if ( testfile[i] == '.') {
+                  strncpy (extens, &testfile[i+1], 4 ) ; 
+                  break ; 
+                }
               }
-            }
-
-            //
-            //     postscript file
-            //
+
+              //
+              //     postscript file
+              //
             
-            if ( ! strcmp(extens, "ps") ) {
-              sprintf ( command, "gv %s/%s &", 
-                        gSystem->WorkingDirectory(), 
-                        testfile ) ; 
-              gSystem->Exec(command) ; 
-              break ; 
-            }
+              if ( ! strcmp(extens, "ps") ) {
+                sprintf ( command, "gv %s/%s &", 
+                          gSystem->WorkingDirectory(), 
+                          testfile ) ; 
+                gSystem->Exec(command) ; 
+                break ; 
+              }
             
-                
+            }   
             
               
@@ -1248 +1474 @@
   default:
     break;
-  }
+
+  }
+
   return kTRUE;
+
 }
 

trunk/MagicSoft/Simulation/Corsika/GuiSimone/MainFrameSimone.hxx

@@ -56 +56 @@
   Char_t            fInputPath[200] ; 
   Int_t             fNumEvents ; 
-  Char_t            fNtupPath[200] ; 
+//  Char_t            fNtupPath[200] ;
+  Char_t            fNtupFile[256] ;
+  Char_t            fCerFileName[256] ;
+  Char_t            fRunSuffix[10] ;
 
+  ifstream          fCerFile ; 
   Int_t             fEvtNb ;