Changeset 870 for trunk/MagicSoft/Simulation/Detector/ReflectorII

Timestamp:

07/18/01 09:37:32 (23 years ago)

Author:

bigongia

Message:

Some debugging printout added. They are activated using option verbose_level 4
C.Bigongiari
:wq

File:

• trunk/MagicSoft/Simulation/Detector/ReflectorII/reflector.c (modified) (10 diffs)

trunk/MagicSoft/Simulation/Detector/ReflectorII/reflector.c

@@ -33 +33 @@
 static long writep = 0L;        /*  write ptr (on rflfile)      */
 
-
 extern CerHeader *cheadp;       /*  var inited in header.c      */
 extern RflHeader *rheadp;       /*  var inited in header.c      */
@@ -43 +42 @@
 static int ProcessEvent(CerHeader *cheadp, FILE *cerfile, FILE *rflfile);
 
+FILE *chkf = NULL;     /**********************/
+long myloop;     /**********************/
+
 void main(void)
 {   long event = 0L;            /*  event counter       */
@@ -48 +50 @@
     /*  Read init & geometry parms, init files and vars */
     init(NULL);
+    chkf = fopen("check", "w");     /**********************/
 
     /*  Processing loop */
@@ -61 +64 @@
     Log(RFLF_CLOSE_LOG);
     fclose(rflfile);
+    fclose(chkf);     /**********************/
 
     /*  Clean memory and exit  */
@@ -78 +82 @@
     extern void makeOmegaI(float theta, float phi);
 
+
     /*  Various counters: phs = absphs + refphs[0..3] + cphs  */
     long phs,           /*  Number of incoming photons  */
@@ -134 +139 @@
                     1.f - Photons[ph].u*Photons[ph].u
                         - Photons[ph].v*Photons[ph].v));
+
             /*  Check absorption  */
-
             if (absorption(wlen, Photons[ph].h, theta))
                 absphs++;
@@ -141 +146 @@
             /*  Check reflection  */
             else if (0 != (ref_type =
-                           ph2cph(&Photons[ph], &CPhotons[cphs]))) 
+                           ph2cph(&Photons[ph], &CPhotons[cphs])))
                 refphs[ref_type-1]++; 
-
             else    /*  Photon passed   */
             {
+            Debug("Ph %d\t%f\t%f\t%f\t%f\t%f\n",ph,
+                   Photons[ph].x,Photons[ph].y,
+                   Photons[ph].u,Photons[ph].v,theta);
+            Debug("CPh %d\t%d\t%f\t%f\n\n",cphs,ph,CPhotons[cphs].x,
+                   CPhotons[cphs].y);
+
+
                 /*  Update first/last arrival times  */
                 if (first > CPhotons[cphs].t) first = CPhotons[cphs].t;
@@ -166 +177 @@
 
                 }   /*  if overflow  */
-            }  /*  else (=Photon passed)  */
+            }   /*  else (=Photon passed)  */
         }   /*  end of loop inside datablock    */
     }   /*  end of loop on datablocks           */
@@ -191 +202 @@
         rheadp->CPhotons     = (long) overflow * NR_OF_CPHOTONS + cphs;
         fwrite(rheadp, sizeof(RflHeader), 1, rflfile);
+        for (myloop=0; myloop<sizeof(RflHeader)/4; myloop++)
+            fprintf(chkf, "%e ", *((float *)rheadp+myloop)); 
+        fputc('\n', chkf);
 
         /*  If there was an overflow, append data from tempfile   */
@@ -252 +266 @@
             /*  Ok: set found at TRUE and exit loop  */
             if (strncmp(cheadp->EVTH, "EVTH", 4) == 0
-             && first_Event <= cheadp->EvtNumber
-             &&                cheadp->EvtNumber <= last_Event
+             && first_Event <= (long)cheadp->EvtNumber
+             &&                (long)cheadp->EvtNumber <= last_Event
              &&    low_Ecut <= cheadp->Etotal
              &&                cheadp->Etotal <= high_Ecut)