Changeset 1431 for trunk/MagicSoft/Simulation/Detector/ReflectorII/geometry.c

Timestamp:

07/24/02 15:36:10 (22 years ago)

Author:

bigongia

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Simulation/Detector/ReflectorII/geometry.c (modified) (7 diffs)

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

@@ -8 +8 @@
 
 extern  char line[];            /*  parsing buf. (init) */
+extern  char axisdev_filename[256], reflectivity_filename[256];
+
+float   mean_refl;              /*  Mirror mean reflectivity 270-610 nm. 
+                                 *  AM June 2002.
+                                 */
 
 float   ct_Focal_mean;          /*  focal dist. (mean) (cm)     */
-float   ct_Focal_std;           /*  focal dist. (std) (cm)      */
 float   ct_PSpread_mean;        /*  pt. spread fn. (mean) (cm)  */
-float   ct_PSpread_std;         /*  pt. spread fn. (std) (cm)   */
-float   ct_Adjustment_std;      /*  adjustment dev. (std) (cm)  */
 float   ct_BlackSpot_rad;       /*  black spot radius (cm)      */
 float   ct_RMirror;             /*  rad. of single mirror (cm)  */
@@ -43 +45 @@
     Log(MIRR_ALLOC_LOG, ct_NMirrors);
     Log(MIRR_TABLE_LOG);
-    if (ct_BinaryData)
-    {   Log(BINF_OPEN__LOG, ct_BinaryName);
-        fread(ct_data, sizeof(mirror), ct_NMirrors, ct_BinaryData);
-        fclose(ct_BinaryData); }
-    else
-    {   /*  read ASCII data  */
-        Log(READ_ASCII_LOG);
-        for (i=0; i<ct_NMirrors; i++)
-        {   if (12 != fscanf(geofile, "%d %f %f %f %f %f %f %f %f %f %f %f",
-                &ct_data[i].i,      &ct_data[i].f,
-                &ct_data[i].sx,     &ct_data[i].sy,
-                &ct_data[i].x,      &ct_data[i].y,      &ct_data[i].z,
-                &ct_data[i].theta,  &ct_data[i].phi,
-                &ct_data[i].xn,     &ct_data[i].yn,     &ct_data[i].zn))
-                break;
-            Log("[%d]", ct_data[i].i);  }
-        Log("%c", '\n');
-        if (i < ct_NMirrors)
-            FatalError(MIRR_FEW___FTL, i);
-        /*  Data Ok: save binary data for next time     */
-        if ((ct_BinaryData=fopen(ct_BinaryName, "w"))==NULL)
-            Log(BINF_ERROR_LOG, ct_BinaryName);
-        else
-        {   Log(BINF_WRITE_LOG, ct_BinaryName);
-            fwrite(ct_data, sizeof(mirror), ct_NMirrors, ct_BinaryData);
-            fclose(ct_BinaryData);      }
-    }   /*  end of if: reading ASCII data  */
+
+    Log(READ_ASCII_LOG);
+    for (i=0; i<ct_NMirrors; i++)
+      {   
+        if (12 != fscanf(geofile, "%d %f %f %f %f %f %f %f %f %f %f %f",
+                         &ct_data[i].i,     &ct_data[i].f,
+                         &ct_data[i].sx,            &ct_data[i].sy,
+                         &ct_data[i].x,     &ct_data[i].y,      &ct_data[i].z,
+                         &ct_data[i].theta,  &ct_data[i].phi,
+                         &ct_data[i].xn,            &ct_data[i].yn,     &ct_data[i].zn))
+          break;
+      }
+    if (i < ct_NMirrors)
+      FatalError(MIRR_FEW___FTL, i);
+
 }   /*  end of ReadMirrorTable  */
 
 static void ReadReflectivity(char *datname)
-{   FILE *datfile = fopen(datname, "r");
+{
+    FILE *datfile = fopen(datname, "r");
     int current = 0;
+
+    mean_refl = 0.;
 
     if (datfile == NULL)
         FatalError(RFLF_ERROR_FTL, datname);
+    else
+      printf("Reading file %s\n", datname);
+
     while (fgets(line, LINE_MAX_LENGTH, datfile))
-    {   if (line[0] == '#') continue;
-        if (nReflectivity == 0)
-        {   nReflectivity = atoi(line);
-            if (nReflectivity)
-            {   if ((Reflectivity[0] =
-                    (float *) malloc(sizeof(float) * nReflectivity)) == NULL
-                 || (Reflectivity[1] =
-                    (float *) malloc(sizeof(float) * nReflectivity)) == NULL)
-                    FatalError(REFL_ALLOC_FTL, nReflectivity);  }}
-        else if (2 == sscanf(line, "%f %f",
-                &Reflectivity[0][current], &Reflectivity[1][current]));
-        {   current++;
-            if (current >= nReflectivity) break; }}
+    {   
+      if (line[0] == '#') continue;
+
+      if (nReflectivity == 0)
+        {   
+          nReflectivity = atoi(line);
+          if (nReflectivity)
+            {
+              if ((Reflectivity[0] =
+                   (float *) malloc(sizeof(float) * nReflectivity)) == NULL ||
+                  (Reflectivity[1] =
+                   (float *) malloc(sizeof(float) * nReflectivity)) == NULL)
+                FatalError(REFL_ALLOC_FTL, nReflectivity);
+            }
+        }
+      else if (2 == sscanf(line, "%f %f", &Reflectivity[0][current], 
+                             &Reflectivity[1][current]))
+      { 
+        // Added June 2002, AM:
+        mean_refl += Reflectivity[1][current];
+        current++;
+        if (current >= nReflectivity) break; 
+      }
+    }
     fclose(datfile);
 
     nReflectivity = current;
+    if (current > 0)
+      mean_refl /= (float) current;
+
 }   /*  end of ReadReflectivity  */
 
@@ -103 +113 @@
     if (datfile == NULL)
         FatalError(AXIS_ERROR_FTL, datname);
+    else
+      printf("Reading file %s\n", axisdev_filename);
+
     if ((AxisDeviation[0]=
         (float *) malloc(sizeof(float) * ct_NMirrors)) == NULL
@@ -149 +162 @@
                  break;
             case focal_distance:
-                 Log(LOG__FLOAT_LOG, "focal distance (average)",
+                 Log(LOG__FLOAT_LOG, "focal distance (average, cm)",
                      ct_Focal_mean = (float) atof(value_ptr));
                  break;
-            case focal_std:
-                 Log(LOG__FLOAT_LOG, "focal distance (std. dev.)",
-                     ct_Focal_std = (float) atof(value_ptr));
+            case focal_std:  /* not implemented. */
                  break;
             case point_spread:
-                 Log(LOG__FLOAT_LOG, "point spread fn. (average)",
+                 Log(LOG__FLOAT_LOG, "point spread fn. sigma (average, cm)",
                      ct_PSpread_mean = (float) atof(value_ptr));
                  break;
-            case point_std:
-                 Log(LOG__FLOAT_LOG, "point spread fn. (std. dev.)",
-                     ct_PSpread_std = (float) atof(value_ptr));
-                 break;
-            case adjustment_dev:
-                 Log(LOG__FLOAT_LOG, "adjustment dev. (std. dev.)",
-                     ct_Adjustment_std = (float) atof(value_ptr));
+            case point_std:  /* not implemented */
+                 break;
+            case adjustment_dev:  /* not implemented */
                  break;
             case black_spot:
@@ -172 +179 @@
                      ct_BlackSpot_rad = (float) atof(value_ptr));
                  break;
+            case n_mirrors:
+                 Log(LOG__INT___LOG, "number of mirrors",
+                     ct_NMirrors = atoi(value_ptr));
+                 break;
             case r_mirror:
                  Log(LOG__FLOAT_LOG, "single mirror radius (cm)",
                      ct_RMirror = (float) atof(value_ptr));
                  break;
-            case n_mirrors:
-                 Log(LOG__INT___LOG, "number of mirrors",
-                     ct_NMirrors = atoi(value_ptr));
-                 break;
             case camera_width:
-                 Log(LOG__FLOAT_LOG, "camera width (cm)",
+                 Log(LOG__FLOAT_LOG, "camera radius (cm)",
                      ct_CameraWidth = (float) atof(value_ptr));
                  break;
@@ -191 +198 @@
                  Log(LOG__FLOAT_LOG, "pixel width (cm)",
                      ct_PixelWidth = (float) atof(value_ptr));
+                 break;
+            case n_centralpixels:
+                 /* this parameter is for camera, not for reflector */
+                 break;
+            case n_gappixels:
+                 /* this parameter is for camera, not for reflector */
                  break;
             case refl_file:
@@ -208 +221 @@
     fclose(geofile);
 
-    if (Reflectivity[0] == NULL) ReadReflectivity(REFLECTIVITY_FILE);
-    if (AxisDeviation[0]== NULL) ReadAxisDev(AXISDEVIATION_FILE);
+    if (strlen(reflectivity_filename) == 0)
+      strcpy(reflectivity_filename, REFLECTIVITY_FILE);
+    if (Reflectivity[0] == NULL) ReadReflectivity(reflectivity_filename);
+
+    if (strlen(axisdev_filename) == 0)
+      strcpy(axisdev_filename, AXISDEVIATION_FILE);
+    if (AxisDeviation[0]== NULL) ReadAxisDev(axisdev_filename);
+
+
     ReadFocals();
 }   /*  end of GeometrySwitch  */
+