Changeset 1685 for trunk/MagicSoft/Simulation/Detector

Timestamp:

12/13/02 10:12:01 (22 years ago)

Author:

blanch

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx (modified) (8 diffs)

trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx

@@ -21 +21 @@
 //
 // $RCSfile: camera.cxx,v $
-// $Revision: 1.48 $
+// $Revision: 1.49 $
 // $Author: blanch $ 
-// $Date: 2002-12-12 17:40:50 $
+// $Date: 2002-12-13 10:04:07 $
 //
 ////////////////////////////////////////////////////////////////////////
@@ -2030 +2030 @@
                     0);
 
+  //  Store qe for each PMT in output file
+  TArrayF qe_pmt;
+  TArrayF wav_pmt;
+  McConfigRunHeader->InitSizePMTs(ct_NPixels);
+  for(int i=0; i<ct_NPixels;i++){
+    McConfigRunHeader->AddPMT(i);
+    MGeomPMT &pmt  = McConfigRunHeader->GetPMT(i);
+    qe_pmt.Set(pointsQE,QE[i][1]);
+    wav_pmt.Set(pointsQE,QElambda);
+    pmt.SetArraySize(pointsQE);
+    pmt.SetPMTContent(i,wav_pmt,qe_pmt);
+  }
+
+  // Store Light Guides factors in the output file
+  TArrayF theta_lg;
+  TArrayF factor_lg;
+  theta_lg.Set(pointsWC,WC[0]);
+  factor_lg.Set(pointsWC,WC[1]);
+  McConfigRunHeader->SetLightGuides(theta_lg,factor_lg);
+
   //  Fill the Header Tree with the current leaves of each branch
   HeaderTree.Fill() ;
@@ -2605 +2625 @@
   
   if ( wcfile.bad() )
-    error( "read_WC", "Cannot open \"%s\". Exiting.\n", QE_FILE );
+    error( "read_WC", "Cannot open \"%s\". Exiting.\n", WC_FILE );
   
   // read file
@@ -2854 +2874 @@
   Float_t camwidth = -1.0;
 
+  Int_t imir;
+  Float_t f,sx,sy,x,y,z,thetan,phin,xn,yn,zn;
+  Float_t dx,dy;
+  Int_t nummir, numref;
+  Float_t wav,ref;
+
   Char_t token[40];
   Char_t line[511];
   Char_t flag;
-  cout<<"To be implemented..."<<endl;
+
   while(1){
     if((flag=fgetc(sp))==EOF)
@@ -2863 +2889 @@
     fgets(&line[1],500,sp);
     line[0]=flag;
+    if (line[0]== '#' && line[2]=='n' && line[3]=='u' && line[4]=='m' &&
+        line[5]== 'b' && line[12]=='d' && line[13]=='a' && line[14]=='t'){
+      fgets(line,500,sp);
+      sscanf(line,"%i",&numref);
+      for(int i=0; i<numref;i++){
+        fgets(line,500,sp);
+        sscanf(line,"%f %f",&wav,&ref);
+      }
+      continue;
+    }
     if (line[0]== '#')
       continue;
     if (strstr(line, "type") == line)
       continue;
-    if (strstr(line, "focal_distance") == line)
+    if (strstr(line, "focal_distance") == line){
       sscanf(line,"%s %f",token,&focal);
       continue;
-    if (strstr(line, "focal_std") == line)
+    }
+    if (strstr(line, "focal_std") == line){
+      sscanf(line,"%s %f",token,&stdfocal);
       continue;
-    if (strstr(line, "point_spread") == line)
+    }
+    if (strstr(line, "point_spread") == line){
+      sscanf(line,"%s %f",token,&point);
       continue;
-    if (strstr(line, "point_std") == line)
+    }
+    if (strstr(line, "point_std") == line){
+      sscanf(line,"%s %f",token,&stdpoint);
       continue;
-    if (strstr(line, "adjustment_dev") == line)
+    }
+    if (strstr(line, "adjustment_dev") == line){
+      sscanf(line,"%s %f",token,&adjust);
       continue;
-    if (strstr(line, "black_spot") == line)
+    }
+    if (strstr(line, "black_spot") == line){
+      sscanf(line,"%s %f",token,&spot);
       continue;
-    if (strstr(line, "camera_width") == line)
+    }
+    if (strstr(line, "camera_width") == line){
+      sscanf(line,"%s %f",token,&camwidth);
       continue;
+    }
     if (strstr(line, "n_pixels") == line)
       continue;
@@ -2890 +2939 @@
     if (strstr(line, "n_gappixels") == line)
       continue;
-    if (strstr(line, "n_mirrors") == line)
+    if (strstr(line, "n_mirrors") == line){
+      sscanf(line,"%s %i",token,&nummir);
+      config->InitSizeMirror(nummir);
       continue;
-    if (strstr(line, "r_mirror") == line)
+    }
+    if (strstr(line, "r_mirror") == line){
+      sscanf(line,"%s %f",token,&radius);
       continue;
-    if (strstr(line, "define_mirrors") == line)
-      break;
-  }
+    }
+    if (strstr(line, "define_mirrors") == line){
+      for(int i=0;i<nummir;i++){
+        fgets(line,500,sp);
+        sscanf(line,"%i %f %f %f %f %f %f %f %f %f %f %f",
+               &imir,&f,&sx,&sy,&x,&y,&z,&thetan,&phin,&xn,&yn,&zn);
+        config->AddMirror(i);
+        MGeomMirror &mirror = config->GetMirror(i);
+        mirror.SetMirrorContent(imir,f,sx,sy,x,y,z,thetan,phin,xn,yn,zn);
+      }
+      fgets(line,500,sp);
+      while(line[2]!='w' || line[3]!='i' || line[4]!='t' || line[5]!='h'){
+        fgets(line,500,sp);
+      }
+      fgets(line,500,sp);
+      for(int i=0;i<nummir;i++){
+      fgets(line,500,sp);
+        sscanf(line,"%f %f",&dx,&dy);
+        MGeomMirror &mirror = config->GetMirror(i);
+        mirror.SetMirrorDeviations(dx,dy);
+      }      
+      continue;
+    }
+  }
+  config->SetMagicDef(radius, focal, stdfocal, point,
+                      stdpoint, adjust, spot, camwidth);
 }
 
@@ -3857 +3933 @@
 //
 // $Log: not supported by cvs2svn $
+// Revision 1.48  2002/12/12 17:40:50  blanch
+// *** empty log message ***
+//
 // Revision 1.47  2002/12/10 17:19:31  blanch
 // *** empty log message ***
@@ -4002 +4081 @@
 //
 // $Log: not supported by cvs2svn $
+// Revision 1.48  2002/12/12 17:40:50  blanch
+// *** empty log message ***
+//
 // Revision 1.47  2002/12/10 17:19:31  blanch
 // *** empty log message ***