Changeset 687 for trunk/MagicSoft/Simulation/Corsika/Mmcs/jcio.c

Timestamp:

03/05/01 16:20:45 (24 years ago)

Author:

harald

Message:

Changes of Ciro and Dennis to write all events of one run in only
one file. Be careful, from now on you need also a new reflector.

File:

• trunk/MagicSoft/Simulation/Corsika/Mmcs/jcio.c (modified) (6 diffs)

trunk/MagicSoft/Simulation/Corsika/Mmcs/jcio.c

@@ -1 +1 @@
 /*----------------------------------------------------------------------
-  jcio.c : 
-
-Rutinas para la gestion de ficheros en CORSIKA 5.20 Las posibilidades
-que contempla son: 
+  jcio.c :    
+
+Rutinas para la gestion de ficheros en CORSIKA 5.20 
+Author J.C.Gonzalez 
+
+Las posibilidades que contempla son: 
 
 - Creacion de un fichero de inicio de run: run###### 
@@ -14 +16 @@
 - La grabacion se realiza en C, por lo que no existen marcas de inicio
   ni final de bloque en los ficheros.
+
+Modified by C.Bigongiari 2000 Dec 19.
+No more RUN###### and END###### files created. 
+Three files, CER######, DAT###### and STA###### files, per run instead 
+of three files per shower.
+Now ###### is the RUN NUMBER not the shower number !!!
 ----------------------------------------------------------------------*/
 
@@ -30 +38 @@
 FILE *patape;
 FILE *sttape;
+/* CBC
 FILE *runtape;
 FILE *endtape;
+CBC*/
 char tpl[100];
 char certpl[100];
 char dattpl[100];
-char statpl[100];
+/* CBC char statpl[100]; */
 char cerfile[100];
 char datfile[100];
+/* CBC
 char stafile[100];
 char runfile[100];
 char endfile[100];
-
-/*--------------------------------------------------
-  jcinitio: 
-  graba un bloque al fichero de particulas 
+CBC */
+
+/*--------------------------------------------------
+  jcinitio:  
+
+Define file names 
 --------------------------------------------------*/
 #ifdef JC_UNDERSCORES
@@ -65 +78 @@
   strcpy (dattpl, tpl);
   strcat (dattpl, "dat%06d");
+  /* CBC
   strcpy (statpl, tpl);
   strcat (statpl, "sta%06d");
-
+  */
+  /* CBC
   sprintf (runfile, "%srun%06d", tpl, *runnum);
   sprintf (endfile, "%send%06d", tpl, *runnum);
+  CBC */
+
+  sprintf (cerfile, certpl,  *runnum);         /* Added by CB */
+  sprintf (datfile, dattpl,  *runnum);         /* Added by CB */
 }
 
@@ -107 +126 @@
   jcstartrun:
   abre el fichero run###### y graba el contenido
---------------------------------------------------*/
-#ifdef JC_UNDERSCORES
-void
-jcstartrun_ (float *runh)
-#else /* JC_NO_UNDERSCORES */
-void
-jcstartrun (float *runh)
-#endif /* JC_UNDERSCORES */
-{
+
+NOW open CER and DAT files  
+--------------------------------------------------*/
+#ifdef JC_UNDERSCORES
+void
+jcstartrun_ (void)
+#else /* JC_NO_UNDERSCORES */
+void
+jcstartrun (void)
+#endif /* JC_UNDERSCORES */
+{
+  /* CBC 
+
   if ((runtape = fopen (runfile, "wb")) == NULL) {
     printf ("JCIO:: Cannot open RUN file %s. Exiting.\n", runfile);
     exit (1);
   }
+  CBC */
+
   /* puts("JCIO:: saving runheader buffer..."); */
+
+  /* CBC
   fwrite (runh, sizeof (float) * MAXBUF, 1, runtape);
-}
-
-/*--------------------------------------------------
-  jcendrun:
-  abre el fichero run###### y graba el contenido
---------------------------------------------------*/
-#ifdef JC_UNDERSCORES
-void
-jcendrun_ (float *rune)
-#else /* JC_NO_UNDERSCORES */
-void
-jcendrun (float *rune)
-#endif /* JC_UNDERSCORES */
-{
-  fwrite (rune, sizeof (float) * MAXBUF, 1, runtape);
-  fclose (runtape);
-}
-
-/*--------------------------------------------------
-  jcnewcerfile:
-  abre un nuevo fichero Cherenkov
---------------------------------------------------*/
-#ifdef JC_UNDERSCORES
-void
-jcnewcerfile_ ()
-#else /* JC_NO_UNDERSCORES */
-void
-jcnewcerfile ()
-#endif /* JC_UNDERSCORES */
-{
-  sprintf (cerfile, certpl, nshow);
+
+  CBC */
+
+    /* Following lines up to CBC added by CB */ 
+
+
   if ((cetape = fopen (cerfile, "wb")) == NULL) {
     printf ("JCIO:: Cannot open CER file %s. Exiting.\n", cerfile);
     exit (1);
   }
-  timefirst = 9.0e10;
-  timelast = -9.0e10;
-}
-
-/*--------------------------------------------------
-  jcnewdatfile:
-  abre un nuevo fichero de particulas
---------------------------------------------------*/
-#ifdef JC_UNDERSCORES
-void
-jcnewdatfile_ (void)
-#else /* JC_NO_UNDERSCORES */
-void
-jcnewdatfile (void)
-#endif /* JC_UNDERSCORES */
-{
-  sprintf (datfile, dattpl, nshow);
   if ((patape = fopen (datfile, "wb")) == NULL) {
     printf ("JCIO:: Cannot open DAT file %s. Exiting.\n", datfile);
     exit (1);
   }
-}
-
-/*--------------------------------------------------
-  jcnewstafile:
-  abre un nuevo fichero de estadisticas
---------------------------------------------------*/
-#ifdef JC_UNDERSCORES
-void
-jcnewstafile_ (void)
-#else /* JC_NO_UNDERSCORES */
-void
-jcnewstafile (void)
-#endif /* JC_UNDERSCORES */
-{
-  sprintf (stafile, statpl, nshow);
-  if ((sttape = fopen (stafile, "wb")) == NULL) {
-    printf ("JCIO:: Cannot open STA file %s. Exiting.\n", stafile);
-    exit (1);
-  }
-}
-
-/*--------------------------------------------------
-  jcnewshower:
-  abre nuevos ficheros para la nueva cascada
---------------------------------------------------*/
-#ifdef JC_UNDERSCORES
-void
-jcnewshower_ (void)
-#else /* JC_NO_UNDERSCORES */
-void
-jcnewshower (void)
-#endif /* JC_UNDERSCORES */
-{
-  if (nshow > 0) {
-    fclose (patape);
-    fclose (cetape);
-    fclose (sttape);
-  }
-
-  nshow++;
-#ifdef JC_UNDERSCORES
-  jcnewdatfile_ ();
-  jcnewcerfile_ ();
-  jcnewstafile_ ();
-#else /* JC_NO_UNDERSCORES */
-  jcnewdatfile ();
-  jcnewcerfile ();
-  jcnewstafile ();
-#endif /* JC_UNDERSCORES */
-}
-
-/*--------------------------------------------------
-  jcenddata:
-  abre el fichero run###### y graba el contenido
---------------------------------------------------*/
-#ifdef JC_UNDERSCORES
-void
-jcenddata_ (float *runh, float *rune)
-#else /* JC_NO_UNDERSCORES */
-void
-jcenddata (float *runh, float *rune)
-#endif /* JC_UNDERSCORES */
-{
-  if ((endtape = fopen (endfile, "wb")) == NULL) {
-    printf ("JCIO:: Cannot open END file %s. Exiting.\n", endfile);
-    exit (1);
-  }
-  fwrite (runh, sizeof (float) * MAXBUF, 1, endtape);
-  fwrite (rune, sizeof (float) * MAXBUF, 1, endtape);
-  fclose (endtape);
-}
+
+  /* CBC */ 
+}
+
+/* WHAT FOLLOWS IS NO MORE USED !!!!!! CB */
+
+/*--------------------------------------------------
+  jcstartshower:
+  
+Function added by CB  
+Sets timefirst and timelast variables to their initial value at the beginnng 
+of new shower   
+--------------------------------------------------*/
+#ifdef JC_UNDERSCORES
+void
+jcstartshower_ (float *evth)
+#else /* JC_NO_UNDERSCORES */
+void
+jcstartshower (float *evth)
+#endif /* JC_UNDERSCORES */
+{
+  timefirst = 9.0e10;
+  timelast = -9.0e10;
+}
+
 
 /*--------------------------------------------------
@@ -357 +296 @@
 
 }
+
+/*--------------------------------------------------
+  jcendrun:
+  
+  Closes STA file 
+--------------------------------------------------*/
+#ifdef JC_UNDERSCORES
+void
+jcendrun_ (float *rune)
+#else /* JC_NO_UNDERSCORES */
+void
+jcendrun (float *rune)
+#endif /* JC_UNDERSCORES */
+{
+  /* fwrite (rune, sizeof (float) * MAXBUF, 1, sttape); */
+  fclose (sttape);
+}