Index: trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx =================================================================== --- trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx (revision 2351) +++ trunk/MagicSoft/Simulation/Detector/Camera/camera.cxx (revision 2352) @@ -21,7 +21,7 @@ // // $RCSfile: camera.cxx,v $ -// $Revision: 1.59 $ +// $Revision: 1.60 $ // $Author: blanch $ -// $Date: 2003-09-22 11:38:47 $ +// $Date: 2003-09-23 16:50:55 $ // //////////////////////////////////////////////////////////////////////// @@ -134,65 +134,16 @@ @"*/ +/*!@" + + And this is the information about the whole telescope. + + @"*/ + //!@{ -static int ct_Type; //@< Type of telescope: 0:CT1, 1:MAGIC -//!@} - -/*!@" - - And this is the information about the whole telescope. - - @"*/ - -//!@{ // parameters of the CT (from the CT definition file) - -////@: Focal distances [cm] -//static float *ct_Focal; - -//@: Mean Focal distances [cm] -static float ct_Focal_mean; - -//@: STDev. Focal distances [cm] -static float ct_Focal_std; - -//@: Mean Point Spread function [cm] -static float ct_PSpread_mean; - -//@: STDev. Point Spread function [cm] -static float ct_PSpread_std; - -//@: STDev. Adjustmente deviation [cm] -static float ct_Adjustment_std; - -//@: Radius of the Black Spot in mirror [cm] -static float ct_BlackSpot_rad; - -//@: Radius of one mirror [cm] -static float ct_RMirror; - -//@: Camera width [cm] -static float ct_CameraWidth; - -//@: Pixel width [cm] -static float ct_PixelWidth; - -//@: ct_PixelWidth_corner_2_corner = ct_PixelWidth / cos(60) -static float ct_PixelWidth_corner_2_corner; - -//@: Number of mirrors -static int ct_NMirrors = 0; //@: Number of pixels static int ct_NPixels; - -//@: Number of pixels -static int ct_NCentralPixels; - -//@: Number of pixels -static int ct_NGapPixels; - -//@: name of the CT definition file to use -static char ct_filename[256]; //@: Number of CT @@ -286,7 +237,4 @@ -//@: coordinates x,y for each pixel -static float **pixary; - //@: contents of the pixels (ph.e.) static float *fnpix; @@ -331,44 +279,4 @@ @"*/ - -//!@{ -// Pointer to a table with the following info.: - -static float **ct_data; - -/* - * TYPE=0 (CT1) - * i s rho theta x y z thetan phin xn yn zn - * - * i : number of the mirror - * s : arc length [cm] - * rho : polar rho of the position of the center of the mirror [cm] - * theta : polar angle of the position of the center of the mirror [cm] - * x : x coordinate of the center of the mirror [cm] - * y : y coordinate of the center of the mirror [cm] - * z : z coordinate of the center of the mirror [cm] - * thetan : polar theta angle of the direction where the mirror points to - * phin : polar phi angle of the direction where the mirror points to - * xn : xn coordinate of the normal vector in the center (normalized) - * yn : yn coordinate of the normal vector in the center (normalized) - * zn : zn coordinate of the normal vector in the center (normalized) - * - * TYPE=1 (MAGIC) - * i f sx sy x y z thetan phin - * - * i : number of the mirror - * f : focal distance of that mirror - * sx : curvilinear coordinate of mirror's center in X[cm] - * sy : curvilinear coordinate of mirror's center in X[cm] - * x : x coordinate of the center of the mirror [cm] - * y : y coordinate of the center of the mirror [cm] - * z : z coordinate of the center of the mirror [cm] - * thetan : polar theta angle of the direction where the mirror points to - * phin : polar phi angle of the direction where the mirror points to - * xn : xn coordinate of the normal vector in the center (normalized) - * yn : yn coordinate of the normal vector in the center (normalized) - * zn : zn coordinate of the normal vector in the center (normalized) - */ -//!@} /*!@" @@ -737,5 +645,4 @@ strcpy( rootname, get_root_filename() ); strcpy( rootname_loop, get_loop_filename() ); - strcpy( ct_filename, get_ct_filename() ); strcpy( nsbpathname, get_nsb_directory() ); strcpy( nsbpath_outer, get_nsb_directory_outer() ); @@ -753,59 +660,66 @@ Starfield_rotate = get_starfield_rotate(); - log(SIGNATURE, - "%s:\n\t%20s: %s\n", - "Starfield Rotate", - "Rotate Starfield", ONoff(Starfield_rotate) - ); // log filenames information - + for(Int_t ict=0;ict (0:CT1 ¦ 1:MAGIC) - - // get focal distance - - sscanf(line, "%s %d", token, &ct_Type); - - log( "read_ct_file", ": %s\n", - ((ct_Type==0) ? "CT1" : "MAGIC") ); - - break; - - case focal_distance: // [cm] - - // get focal distance - - sscanf(line, "%s %f", token, &ct_Focal_mean); - - log( "read_ct_file", ": %f cm\n", ct_Focal_mean ); - - break; - - case focal_std: // s(focal distance) [cm] - - // get focal distance - - sscanf(line, "%s %f", token, &ct_Focal_std); - - log( "read_ct_file", "s(Focal distance): %f cm\n", ct_Focal_std ); - - break; - - case point_spread: // [cm] - - // get point spread - - sscanf(line, "%s %f", token, &ct_PSpread_mean); - - log( "read_ct_file", ": %f cm\n", ct_PSpread_mean ); - - break; - - case point_std: // s(point spread) [cm] - - // get point spread - - sscanf(line, "%s %f", token, &ct_PSpread_std); - - log( "read_ct_file", "s(Point spread): %f cm\n", ct_PSpread_std ); - - break; - - case adjustment_dev: // s(adjustment_dev) [cm] - - // get point spread - - sscanf(line, "%s %f", token, &ct_Adjustment_std); - - log( "read_ct_file", "s(Adjustment): %f cm\n", ct_Adjustment_std ); - - break; - - case black_spot: // radius of the black spot in the center [cm] - - // get black spot radius - - sscanf(line, "%s %f", token, &ct_BlackSpot_rad); - - log( "read_ct_file", "Radius of the black spots: %f cm\n", - ct_BlackSpot_rad); - - break; - - case r_mirror: // radius of the mirrors [cm] - - // get radius of mirror - - sscanf(line, "%s %f", token, &ct_RMirror); - - log( "read_ct_file", "Radii of the mirrors: %f cm\n", ct_RMirror ); - - break; - - case n_mirrors: // number of mirrors - - // get the name of the output_file from the line - - sscanf(line, "%s %d", token, &ct_NMirrors); - - log( "read_ct_file", "Number of mirrors: %d\n", ct_NMirrors ); - - break; - - case camera_width: // camera width [cm] - - // get the name of the ct_file from the line - - sscanf(line, "%s %f", token, &ct_CameraWidth); - - log( "read_ct_file", "Camera width: %f cm\n", ct_CameraWidth ); - - break; - - case n_pixels: // number of pixels - - // get the name of the output_file from the line - - sscanf(line, "%s %d", token, &ct_NPixels); - - log( "read_ct_file", "Number of pixels: %d\n", ct_NPixels ); - - break; - - case n_centralpixels: // number of central pixels - - // get the name of the output_file from the line - - sscanf(line, "%s %d", token, &ct_NCentralPixels); - - log( "read_ct_file", "Number of central pixels: %d\n", ct_NCentralPixels ); - - break; - - case n_gappixels: // number of gap pixels - - // get the name of the output_file from the line - - sscanf(line, "%s %d", token, &ct_NGapPixels); - - log( "read_ct_file", "Number of gap pixels: %d\n", ct_NGapPixels ); - - break; - - case pixel_width: // pixel width [cm] - - // get the name of the ct_file from the line - - sscanf(line, "%s %f", token, &ct_PixelWidth); - - ct_PixelWidth_corner_2_corner = ct_PixelWidth / cos(RAD(30.0)); - - log( "read_ct_file", "Pixel width: %f cm\n", ct_PixelWidth ); - - break; - - case define_mirrors: // read table with the parameters of the mirrors - - log( "read_ct_file", "Table of mirrors data:\n" ); - - // check whether the number of mirrors was already set - - if ( ct_NMirrors == 0 ) - error( "read_ct_file", "NMirrors was not set.\n" ); - - // allocate memory for paths list - - log( "read_ct_file", "Allocating memory for ct_data\n" ); - - ct_data = new float*[ct_NMirrors]; - - for (i=0; i> ct_data[i][j]; - - break; - - } // switch ( i ) - - } // end while - - // end - - log( "read_ct_file", "done.\n" ); - - return; -} -//!@} - - -//!----------------------------------------------------------- // @name read_QE // @@ -3209,195 +2892,4 @@ //!@} -//!----------------------------------------------------------- -// @name read_pixels -// -// @desc read pixels data -// -// @date Fri Mar 12 16:33:34 MET 1999 -//------------------------------------------------------------ -// @function - -//!@{ -void -read_pixels(struct camera *pcam) -{ - /* ifstream qefile; - char line[LINE_MAX_LENGTH]; - int n, i, j, icount; - float qe; - - //------------------------------------------------------------ - // first, pixels' coordinates - - pcam->inumpixels = ct_NPixels; - pcam->inumcentralpixels = ct_NCentralPixels; - pcam->inumgappixels = ct_NGapPixels; - pcam->inumbigpixels = ct_NPixels - ct_NCentralPixels - ct_NGapPixels; - pcam->dpixdiameter_cm = ct_PixelWidth; - - // initialize pixel numbers - - pixary = new float* [2*ct_NCentralPixels]; - for ( i=0; i<2*ct_NCentralPixels; ++i ) - pixary[i] = new float[2]; - - pixneig = new int* [ct_NCentralPixels]; - for ( i=0; i -1) && - ((j-1) < pointsQE) && ((j-1) > -1) ) { - QE[i-1][0][j-1] = QElambda[j-1]; - QE[i-1][1][j-1] = qe; - } - - if ( i > ct_NPixels) break; - - icount++; - - } - - if(icount/pointsQE < ct_NPixels){ - error( "read_pixels", "The quantum efficiency file is faulty\n (found only %d pixels instead of %d).\n", - icount/pointsQE, ct_NPixels ); - } - - // close file - - qefile.close(); - - // test QE - - for(icount=0; icount< ct_NPixels; icount++){ - for(i=0; i 1000. || - QE[icount][1][i] < 0. || QE[icount][1][i] > 100.){ - error( "read_pixels", "The quantum efficiency file is faulty\n pixel %d, point %d is % f, %f\n", - icount, i, QE[icount][0][i], QE[icount][1][i] ); - } - } - } - - // end - - log("read_pixels", "Done.\n"); - */ -} -//!@} //!----------------------------------------------------------- @@ -3534,29 +3026,4 @@ } -//!----------------------------------------------------------- -// @name pixels_are_neig -// -// @desc check whether two pixels are neighbours -// -// @var pix1 Number of the first pixel -// @var pix2 Number of the second pixel -// @return TRUE: both pixels are neighbours; FALSE: oth. -// -// @date Wed Sep 9 17:58:37 MET DST 1998 -//------------------------------------------------------------ -// @function - -//!@{ -int -pixels_are_neig(int pix1, int pix2) -{ - if ( sqrt(SQR( pixary[pix1][0] - pixary[pix2][0] ) + - SQR( pixary[pix1][1] - pixary[pix2][1] ) ) - > ct_PixelWidth_corner_2_corner ) - return ( FALSE ); - else - return ( TRUE ); -} -//!@} //!-----------------------------------------------------------