Changeset 366 for trunk/MagicSoft/Simulation

Timestamp:

02/18/00 17:42:39 (25 years ago)

Author:

petry

Message:

This version includes drastic changes and belongs to camera.cxx 1.5.
It is not yet finished and not immediately useful because the
trigger simulation is not yet re-implemented. I had to take it
out together with some other stuff in order to tidy the whole
program up. This is not meant as an insult to anyone. I needed
to do this in order to be able to work on it.

This version has been put in the repository in order to be
able to share the further development with others.

If you need something working, wait or take an earlier one.
See file README.

File:

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

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

@@ -19 +19 @@
 //=
 //= $RCSfile: camera.h,v $
-//= $Revision: 1.3 $
-//= $Author: harald $ 
-//= $Date: 1999-11-11 20:29:29 $
+//= $Revision: 1.4 $
+//= $Author: petry $ 
+//= $Date: 2000-02-18 17:42:39 $
 //=
 //=//////////////////////////////////////////////////////////////////////
@@ -61 +61 @@
 #include "moments.h"
 
-#include "lagrange.h"
-
 #include "MCEventHeader.hxx"
 #include "MCCphoton.hxx"
+#include "MTrigger.hxx"
 
 // command line options available
@@ -96 +95 @@
 //!@{
 
-#define SLICES            100
+#define SLICES            19
 #define WIDTH_TIMESLICE   3.3
  
+#define NOTIME 999999999
+
 #define SIN60   0.866025403784439
 #define COS60   0.500000000000000
@@ -104 +105 @@
 #define RandomNumber  drand48()
 
-#define PIX_ARRAY_SIDE       40
-#define PIX_ARRAY_HALF_SIDE  20
-#define PIXNUM               0
-#define PIXX                 1
-#define PIXY                 2
-
 #define iMAXNUMPIX  595 // total maximum possible number of pixels in the camera
 
-//@ the trigger threshold up to which the maximum passable threshold is tested 
-#define iMAX_THRESHOLD_PHE   50 
-
-//@ number of the 1st. pixel of a sector s in a ring r (central pixel: ring=0)
-#define FIRST_PIXEL(r,s)   ( ((r)>0) ? (3*(r)*((r)-1) + (r)*(s) + 1) : 0 )
-
-//@ number of the pixels include in a camera of r pixels
-#define NUMBER_PIXELS(r)   ( ((r)>0) ? FIRST_PIXEL((r)+1,0) : 1 )
+#define iMAXNUMPHE 50000 // maximum total number of photoelectrons in one event
+
+#define iNUMNSBPRODCALLS 1  // in order to average over the QE simulation, call the
+                            // phe production function for the NSB iNUMNSBPRODCALLS times
+
+#define iNUMWAVEBANDS 5  // number of wavebands for the NSB simulation
+
+#define WAVEBANDBOUND1 290. // iNUMWAVEBANDS+1 boundaries for waveband definitions (nm)
+#define WAVEBANDBOUND2 310. // for the NSB generation
+#define WAVEBANDBOUND3 400.
+#define WAVEBANDBOUND4 500.
+#define WAVEBANDBOUND5 600.
+#define WAVEBANDBOUND6 800.
+
+#define EXTWAVEBAND1 3.715  // iNUMWAVEBANDS extinction values in magnitudes per airmass
+#define EXTWAVEBAND2 0.642  // - these values were taken from D.L. King, Isaac Newton Group
+#define EXTWAVEBAND3 0.209  // Tech Note No. 31, they are for a clear night at the ORM,
+#define EXTWAVEBAND4 0.107  // La Palma. The values were averaged in the given wavebands
+#define EXTWAVEBAND5 0.053
+
+#define SIMTIMEOFFSET_NS 5  // determines how many ns before the first and after the last
+                            // shower photoelectron, there should be NSB photoelectrons 
 
 //@ now we define the list CT_ITEM_LIST of possible items in the CT def. file
@@ -225 +234 @@
   double dxpointcorr_deg; /* correction of the pixel coordinates; to be added to dxc[] to get correct value */
   double dypointcorr_deg; /* correction of the pixel coordinates; to be added to dxc[] to get correct value */
-  double di[iMAXNUMPIX]; /* i coordinate in JCs bi-axis hexagonal coordinate system */
-  double dj[iMAXNUMPIX]; /* j coordinate in JCs bi-axis hexagonal coordinate system */
  
+};
+
+class Photoelectron{
+ public:
+  int iarrtime_ns; // arrival time taken from the corresponding photon
+  int ipixnum;     // pixel in which the phe was evoked (starts at 0)
+  
+  Photoelectron(); // constructor (code in camera.cxx)
+
 };
 
@@ -268 +284 @@
                 float l, float m, float n,
                 float x, float y, float z);
-     
+int check_reflector_file(FILE *infile);
+float lin_interpol(float x1, float y1, float x2, float y2, float x);
+int produce_phes( FILE *sp, // the input file
+                  struct camera *cam, // the camera layout
+                  float minwl_nm, // the minimum accepted wavelength
+                  float maxwl_nm, // the maximum accepted wavelength
+                  class Photoelectron phe[iMAXNUMPHE], // the generated phes
+                  int *itotnphe, // total number of produced photoelectrons
+                  float *nphe, // number of photoelectrons in each pixel
+                  int *incph,    // total number of cph read
+                  float *tmin_ns,   // minimum arrival time of all phes
+                  float *tmax_ns    // maximum arrival time of all phes
+                  );
+
+int produce_nsbrates( char *inname,
+                      struct camera *cam,
+                      class Photoelectron phe[iMAXNUMPHE],
+                      float nsbrate[][iNUMWAVEBANDS]
+                      );
+
+int produce_nsb_phes( float *atmin_ns,
+                      float *atmax_ns,
+                      float theta_rad,
+                      struct camera *cam,
+                      float nsbr_phepns[iMAXNUMPIX][iNUMWAVEBANDS],
+                      float dnsb_phepns[iMAXNUMPIX],
+                      float extinction[iNUMWAVEBANDS],
+                      float fnpx[iMAXNUMPIX],
+                      Photoelectron photo[iMAXNUMPHE],
+                      int *inphe,
+                      float base_mv[iMAXNUMPIX]); 
+
 //!@}
 
@@ -278 +325 @@
 /*
  *$Log: not supported by cvs2svn $
+ *Revision 1.3  1999/11/11 20:29:29  harald
+ *Small changes to run the new version on a linux machine.
+ *
  *Revision 1.2  1999/11/10 07:42:41  harald
  *Small change to read the right data files in.
@@ -290 +340 @@
  *The "rootification" was done by Dirk Petry and Harald Kornmayer. 
  *
- *In the following you can see the README file of that version:
- *
- *==================================================
- *
- *Fri Oct 22  1999   D.P.
- *
- *The MAGIC Monte Carlo System
- *
- *Camera Simulation Programme
- *---------------------------
- *
- *1) Description
- *
- *This version is the result of the fusion of H.K.'s
- *root_camera which is described below (section 2)
- *and another version by D.P. which had a few additional
- *useful features.
- *
- *The version compiles under Linux with ROOT 2.22 installed
- *(variable ROOTSYS has to be set).
- *
- *Compile as before simply using "make" in the root_camera
- *directory.
- *
- *All features of H.K.'s root_camera were retained.
- *
- *Additional features of this version are:
- *
- *  a) HBOOK is no longer used and all references are removed.
- *
- *  b) Instead of HBOOK, the user is given now the possibility of 
- *     having Diagnostic data in ROOT format as a complement
- *     to the ROOT Raw data.
- *
- *     This data is written to the file which is determined by
- *     the new input parameter "diag_file" in the camera parameter
- *     file.
- *
- *     All source code file belonging to this part have filenames
- *     starting with "MDiag".
- *
- *     The user can read the output file using the following commands
- *     in an interactive ROOT session:
- *
- *              root [0] .L MDiag.so
- *      root [1] new TFile("diag.root");
- *      root [2] new TTreeViewer("T");
- *      
- *     This brings up a viewer from which all variables of the
- *     TTree can be accessed and histogrammed. This example
- *     assumes that you have named the file "diag.root", that
- *     you are using ROOT version 2.22 or later and that you have
- *     the shared object library "MDiag.so" which is produced
- *     by the Makefile along with the executable "camera".
- *       
- * !   The contents of the so-called diag file is not yet fixed.
- * !   At the moment it is what J.C.G. used to put into the HBOOK
- * !   ntuple. In future versions the moments calculation can be
- * !   removed and the parameter list be modified correspondingly.
- *
- *  c) Now concatenated reflector files can be read. This is useful
- *     if you have run the reflector with different parameters but
- *     you want to continue the analysis with all reflector data
- *     going into ONE ROOT outputfile.
- *
- *     The previous camera version contained a bug which made reading 
- *     of two or more concatenated reflector files impossible.
- *
- *  d) The reflector output format was changed. It is now version
- *     0.4 .
- *     The change solely consists in a shortening of the flag
- *     definition in the file 
- *
- *           include-MC/MCCphoton.hxx  
- *
- * !   IF YOU WANT TO READ REFLECTOR FORMAT 0.3, you can easily
- * !   do so by recompiling camera with the previous version of
- * !   include-MC/MCCphoton.hxx.
- *
- *     The change was necessary for saving space and better
- *     debugging. From now on, this format can be frozen.
- *
- * !   For producing reflector output in the new format, you
- * !   of course have to recompile your reflector with the
- * !   new include-MC/MCCphoton.hxx .
- *
- *  e) A first version of the pixelization with the larger
- *     outer pixels is implemented. THIS IS NOT YET FULLY
- *     TESTED, but first rough tests show that it works
- *     at least to a good approximation.
- *
- *     The present version implements the camera outline
- *     with 18 "gap-pixels" and 595 pixels in total as
- *     shown in 
- *
- *        http://sarastro.ifae.es/internal/home/hardware/camera/numbering.ps
- *
- *     This change involved 
- *
- *      (i) The file pixels.dat is no longer needed. Instead
- *          the coordinates are generated by the program itself
- *          (takes maybe 1 second). In the file 
- *
- *              pixel-coords.txt
- *
- *        in the same directory as this README, you find a list
- *          of the coordinates generated by this new routine. It
- *          has the format
- *
- *              number   i   j   x  y  size-factor
- *
- *          where i and j are J.C.G.'s so called biaxis hexagonal
- *          coordinates (for internal use) and x and y are the
- *          coordinates of the pixel centers in the standard camera
- *          coordinate system in units of centimeters. The value
- *          of "size-factor" determines the linear size of the pixel
- *          relative to the central pixels. 
- *
- *        (ii) The magic.def file has two additional parameters
- *          which give the number of central pixels and the
- *          number of gap pixels
- *
- *        (iii) In camera.h and camera.cxx several changes were 
- *          necessary, among them the introduction of several
- *          new functions 
- *
- *     The newly suggested outline with asymmetric Winston cones
- *     will be implemented in a later version.
- *
- *  f) phe files can no longer be read since this contradicts
- *     our philosophy that the analysis should be done with other
- *     programs like e.g. EVITA and not with "camera" itself.
- *     This possibility was removed. 
- *
- *  g) ROOT is no longer invoked with an interactive interface.
- *     In this way, camera can better be run as a batch program and
- *     it uses less memory.
- *
- *  h) small changes concerning the variable "t_chan" were necessary in
- *     order to avoid segmentation faults: The variable is used as an
- *     index and it went sometimes outside the limits when camera
- *     was reading proton data. This is because the reflector files
- *     don't contain the photons in a chronological order and also
- *     the timespread can be considerably longer that the foreseen
- *     digitisation timespan. Please see the source code of camera.cxx
- *     round about line 1090.
- *
- *  j) several unused variables were removed, a few warning messages
- *     occur when you compile camera.cxx but these can be ignored at
- *     the moment.
- *
- *In general the program is of course not finished. It still needs
- *debugging, proper trigger simulation, simulation of the asymmetric
- *version of the outer pixels, proper NSB simulation, adaption of
- *the diag "ntuple" contents to our need and others small improvements.
- *
- *In the directory rfl-files there is now a file in reflector format 0.4
- *containing a single event produced by the starfiled adder. It has
- *a duration of 30 ns and represents the region around the Crab Nebula.
- *Using the enclosed input parameter file, camera should process this
- *file without problems.
- *
- *2) The README for the previous version of root_camera
- *
- *README for a preliminary version of the 
- *root_camera program. 
- *
- *root_camera is based on the program "camera"of Jose Carlos
- *Gonzalez. It was changed in the way that only the pixelisation 
- *and the distibution of the phe to the FADCs works in a 
- *first version. 
- *
- *Using the #undef command most possibilities of the orignal 
- *program are switched of. 
- *
- *The new parts are signed by 
- *
- *- ROOT or __ROOT__ 
- *  nearly all  important codelines for ROOT output are enclosed 
- *  in structures like 
- *  #ifdef __ROOT__ 
- *  
- *    code 
- *
- *  #endif __ROOT__ 
- *
- *  In same case the new lines are signed by a comment with the word 
- *  ROOT in it. 
- *
- *  For timing of the pulse some variable names are changed. 
- *  (t0, t1, t  -->  t_ini, t_fin, t_1st, t_chan,...) 
- *  Look also for this changes. 
- *
- *  For the new root-file is also a change in readparm-files
- *
- *
- *- __DETAIL_TRIGGER__
- *
- *  This is for the implementation of the current work on trigger 
- *  studies. Because the class MTrigger is not well documented it 
- *  isn´t a part of this tar file. Only a dummy File exists. 
- *
- *
- *
- *With all files in the archive, the root_camera program should run. 
- *
- *A reflector file is in the directory rfl-files
- *
- *==================================================
- *
- *From now on, use CVS for development!!!!
- *
- *
- *
  *Revision 1.3  1999/10/22 15:32:56  petry
  *tidied-up version, really sent to H.K. and N.M., 22-10-99