Changeset 5269 for trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx

Timestamp:

10/13/04 18:05:05 (20 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

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

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

@@ -21 +21 @@
 //
 // $RCSfile: camera.cxx,v $
-// $Revision: 1.73 $
+// $Revision: 1.74 $
 // $Author: moralejo $ 
-// $Date: 2004-10-13 14:28:59 $
+// $Date: 2004-10-13 17:05:05 $
 //
 ////////////////////////////////////////////////////////////////////////
@@ -238 +238 @@
 static  float Spot_y=0.0;
 static  float Spotsigma=0.0;
-
-static  int CalibrationRun = 0;
 
   
@@ -717 +715 @@
   // get filenames
 
-  CalibrationRun = 1;
-  if (! CalibrationRun)
+  if (! is_calibration_run())
     for(int ict=0;ict<ct_Number;ict++)
       {
@@ -1682 +1679 @@
 
 
-  if (CalibrationRun)
+  if (is_calibration_run())
     {
       DoCalibration(Fadc_CT, Trigger_CT, camgeom, nsb_trigresp, nsb_fadcresp,
@@ -4455 +4452 @@
 {
 
-  int ntshow = 0;
+  int nevent = 0;
+  int maxevents;
   int lons_return;
 
@@ -4464 +4462 @@
   float timefirst = 0.;
   float timelast = 0.;
+  
+  float lambda, sigma_lambda, phot_per_pix, sigma_time;
+  UInt_t *first_pixel, *last_pixel;
+  int selected_pixel;
+
+  // Get parameters of the calibration:
+
+  get_calibration_properties( &lambda, &sigma_lambda, &phot_per_pix, 
+                              &sigma_time, &maxevents, &selected_pixel);
+
+  first_pixel = new UInt_t[ct_Number];
+  last_pixel  = new UInt_t[ct_Number];
+
+  for (int ict = 0; ict < ct_Number; ict++)
+    {
+      if (selected_pixel < 0)
+        {
+          first_pixel[ict] = 0;
+          last_pixel[ict]  = ((MGeomCam*)(camgeom.UncheckedAt(ict)))->GetNumPixels()-1;
+        }
+      else
+        first_pixel[ict] = last_pixel[ict] = selected_pixel;
+    }
+
 
   TArrayC *fadcValues;     //@< the analog Fadc High gain signal for pixels
@@ -4476 +4498 @@
   float *pheinpix = new float [ct_NPixels];
 
-  for (int calevent = 0; calevent < 1; calevent++) 
+  for (int calevent = 0; calevent < maxevents; calevent++) 
     { 
       //
@@ -4493 +4515 @@
         }
 
-      ntshow++;
-
-      if((ntshow+1)%100 == 1)
-        log(SIGNATURE, "Event %d\n", ntshow);
+      nevent++;
+
+      if((nevent+1)%100 == 1)
+        log(SIGNATURE, "Event %d\n", nevent);
         
       // Produce the photoelectrons
-      float phot_per_pix = 100.;
       for(int ict = 0; ict < ct_Number; ict++)
         produce_calib_phes( (MGeomCam*)(camgeom.UncheckedAt(ict)),
-                            phot_per_pix,
                             Trigger_CT[ict],
                             Fadc_CT[ict],
@@ -4508 +4528 @@
                             pheinpix,
                             &(ntotphot[ict]),
-                            ict);
+                            ict,
+                            lambda,
+                            sigma_lambda,
+                            phot_per_pix,
+                            sigma_time,
+                            selected_pixel
+                            );
 
       // NSB simulation
@@ -4518 +4544 @@
           for (int ict = 0; ict < ct_Number; ict++)
             {
-              for (UInt_t ui = 0;
-                  ui < ((MGeomCam*)(camgeom.UncheckedAt(ict)))->GetNumPixels();
-                  ui++)
+              for (UInt_t ui = first_pixel[ict]; ui <= last_pixel[ict]; ui++)
                 {
                   nphe_from_nsb[ict] += nsb_phepns[ict][ui];
@@ -4593 +4617 @@
           // Fill the FADC information
           //
-          for(UInt_t ui=0;
-              ui<((MGeomCam*)(camgeom.UncheckedAt(ict)))->GetNumPixels();
-              ui++)
+          for(UInt_t ui = first_pixel[ict]; ui <= last_pixel[ict]; ui++)
             {
               for (int jslice = 0; jslice < FADC_slices_written; jslice++)
@@ -4628 +4650 @@
 //------------------------------------------------------------
 
-int produce_calib_phes( MGeomCam *camgeom,      // The camera layout
-                        float     phot_per_pix, // Average number of photons per inner pixel
+int produce_calib_phes( MGeomCam *camgeom,       // The camera layout
                         MTrigger *trigger,
                         MFadc    *fadc,
-                        int      *itotnphe,     // total number of produced photoelectrons              
-                        float    *nphe,         // number of photoelectrons in each pixel
-                        int      *nphot,        // total number of photons in all pixels
-                        int       ict           // Telescope that is being analised to get the right QE.
+                        int      *itotnphe,      // total number of produced photoelectrons             
+                        float    *nphe,          // number of photoelectrons in each pixel
+                        int      *nphot,         // total number of photons in all pixels
+                        int       ict,           // Telescope that is being analised to get the right QE.
+                        float     lambda,        // Mean wavelength of light in nm
+                        float     sigma_lambda,  // Sigma of wavelengtgaussian spread
+                        float     phot_per_pix,  // Average # of photons per inner pixel
+                        float     sigma_time,    // Sigma of time spread of photons
+                        int       selected_pixel // if >= 0, only this pixel is used! 
                         )
 {
@@ -4643 +4669 @@
   float **qept;
   float radius_mm, focal_dist_mm;
+  int total_photons;
 
   // reset variables
@@ -4671 +4698 @@
 
   // Calculate total number of photons to be produced.
-  int total_photons = (int) (phot_per_pix * 3.14159265 * radius_mm * radius_mm / 
-                             (*camgeom)[0].GetA());
+  if (selected_pixel < 0)
+    total_photons = (int) (phot_per_pix * 3.14159265 * radius_mm * radius_mm / 
+                           (*camgeom)[0].GetA());
+  else
+    total_photons = (int) (phot_per_pix * 
+                           (*camgeom)[selected_pixel].GetA() / (*camgeom)[0].GetA());
 
   // loop over the photons
@@ -4679 +4710 @@
     {
 
-      time = 0.;
-
-      // Obtain photon coordinates on the camera. We assume a point source of light placed
-      // in the center of the mirror dish.
-
-      // polar angle 
-      float psi = RandomNumber * 2 * 3.14159265;
-      // angle between the telescope axis and the photon trajectory.
-      float epsilon = acos(1.-RandomNumber*(1.-cos_epsilon_max));
-      float tanepsilon = tan(epsilon);
-
-      cx = focal_dist_mm*tanepsilon*cos(psi); // mm
-      cy = focal_dist_mm*tanepsilon*sin(psi); // mm
-
-      // Angle between photon trajectory and camera plane:
-      phi = 3.14159265/2.-epsilon; // rad
+      // time
+      time = random.Gaus(4.*sigma_time, sigma_time);
+
+      cout << "### " << time << endl << endl;
 
       // wavelength
-          
-      wl = 500.; // nm
-
-
-      if( (wl > WAVEBANDBOUND6) || (wl < WAVEBANDBOUND1) || 
-          (sqrt(cx*cx + cy*cy) > radius_mm ) )
+      wl = lambda;
+
+      if( (wl > WAVEBANDBOUND6) || (wl < WAVEBANDBOUND1))
         continue;
 
-      //
-      // Pixelization
-      //
-      ipixnum = bpoint_is_in_pix(cx, cy, camgeom);
-
-      // -1 = the photon is in none of the pixels
-      //  0 = the phton is in the central pixel, which is not used
-      if (ipixnum==-1 || ipixnum==0)
-        continue;
+      if (selected_pixel < 0)
+        {
+          // Obtain photon coordinates on the camera. We assume a point source of light placed
+          // in the center of the mirror dish.
+
+          // polar angle 
+          float psi = RandomNumber * 2 * 3.14159265;
+          // angle between the telescope axis and the photon trajectory.
+          float epsilon = acos(1.-RandomNumber*(1.-cos_epsilon_max));
+          float tanepsilon = tan(epsilon);
+
+          cx = focal_dist_mm*tanepsilon*cos(psi); // mm
+          cy = focal_dist_mm*tanepsilon*sin(psi); // mm
+ 
+          if (sqrt(cx*cx + cy*cy) > radius_mm )
+            continue;
+
+          // Angle between photon trajectory and camera plane:
+          phi = 3.14159265/2.-epsilon; // rad
+
+          //
+          // Pixelization
+          //
+          ipixnum = bpoint_is_in_pix(cx, cy, camgeom);
+
+          // -1 = the photon is in none of the pixels
+          //  0 = the phton is in the central pixel, which is not used
+          if (ipixnum==-1 || ipixnum==0)
+            continue;
+        }
+      else
+        {
+          // Angle between photon trajectory and camera plane:
+          phi = atan2( focal_dist_mm, 
+                       sqrt( (*camgeom)[selected_pixel].GetX()*(*camgeom)[selected_pixel].GetX()+
+                             (*camgeom)[selected_pixel].GetY()*(*camgeom)[selected_pixel].GetY()));
+
+          ipixnum = selected_pixel;
+        }
+
 
       // increase number of photons within pixels
@@ -4793 +4840 @@
       *itotnphe += 1;
     }
-
+  cout << "$$" << endl;
   return(0);