Index: trunk/MagicSoft/Mars/mgeom/GeomIncl.h =================================================================== --- trunk/MagicSoft/Mars/mgeom/GeomIncl.h (revision 1436) +++ trunk/MagicSoft/Mars/mgeom/GeomIncl.h (revision 1436) @@ -0,0 +1,3 @@ +#ifndef __CINT__ + +#endif // __CINT__ Index: trunk/MagicSoft/Mars/mgeom/GeomLinkDef.h =================================================================== --- trunk/MagicSoft/Mars/mgeom/GeomLinkDef.h (revision 1436) +++ trunk/MagicSoft/Mars/mgeom/GeomLinkDef.h (revision 1436) @@ -0,0 +1,12 @@ +#ifdef __CINT__ + +#pragma link off all globals; +#pragma link off all classes; +#pragma link off all functions; + +#pragma link C++ class MGeomPix+; +#pragma link C++ class MGeomCam+; +#pragma link C++ class MGeomCamCT1+; +#pragma link C++ class MGeomCamMagic+; + +#endif Index: trunk/MagicSoft/Mars/mgeom/MGeomCam.cc =================================================================== --- trunk/MagicSoft/Mars/mgeom/MGeomCam.cc (revision 1436) +++ trunk/MagicSoft/Mars/mgeom/MGeomCam.cc (revision 1436) @@ -0,0 +1,119 @@ +/* ======================================================================== *\ +! +! * +! * This file is part of MARS, the MAGIC Analysis and Reconstruction +! * Software. It is distributed to you in the hope that it can be a useful +! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. +! * It is distributed WITHOUT ANY WARRANTY. +! * +! * Permission to use, copy, modify and distribute this software and its +! * documentation for any purpose is hereby granted without fee, +! * provided that the above copyright notice appear in all copies and +! * that both that copyright notice and this permission notice appear +! * in supporting documentation. It is provided "as is" without express +! * or implied warranty. +! * +! +! +! Author(s): Thomas Bretz 12/2000 +! Author(s): Harald Kornmayer 1/2001 +! +! Copyright: MAGIC Software Development, 2000-2002 +! +! +\* ======================================================================== */ + +/////////////////////////////////////////////////////////////////////// +// +// MGeomCam +// +// This is the base class of different camera geometries. It creates +// a pixel object for a given number of pixels and defines the +// interface of how to acccess the geometry information. +// +/////////////////////////////////////////////////////////////////////// + +#include "MGeomCam.h" + +#include "MLog.h" +#include "MLogManip.h" + +#include "MGeomPix.h" +#include "MHexagon.h" + +ClassImp(MGeomCam); + +// -------------------------------------------------------------------------- +// +// Initializes a Camera Geometry with npix pixels. All pixels +// are deleted when the corresponding array is deleted. +// +MGeomCam::MGeomCam(UInt_t npix, Float_t dist, const char *name, const char *title) + : fNumPixels(npix), fCamDist(dist), fConvMm2Deg(kRad2Deg/(dist*1000)) +{ + fName = name ? name : "MGeomCam"; + fTitle = title ? title : "Storage container for a camera geometry"; + + fPixels = new TObjArray(npix); + + // + // make sure that the destructor delete all contained objects + // + fPixels->SetOwner(); + + for (UInt_t i=0; ifMaxRadius) + fMaxRadius = maxr; + } +} + +// -------------------------------------------------------------------------- +// +// Prints the Geometry information of all pixels in the camera +// +void MGeomCam::Print(Option_t *) const +{ + // + // Print Information about the Geometry of the camera + // + *fLog << all << " Number of Pixels (" << GetTitle() << "): " << fNumPixels << endl; + + for (UInt_t i=0; i +#endif + +class MGeomPix; + +class MGeomCam : public MParContainer +{ +private: + UInt_t fNumPixels; // Number of pixels in this camera + Float_t fMaxRadius; // maximum radius of the camera (eg. for GUI layout) + + Float_t fCamDist; // [m] Average distance of the camera from the mirror + Float_t fConvMm2Deg; // conversion factor to convert mm in the camera plain into degrees + + TObjArray *fPixels; // Array of singel pixels storing the geometry + +protected: + void CalcMaxRadius(); + +public: + + MGeomCam(UInt_t npix, Float_t dist, const char *name=NULL, const char *title=NULL); + + virtual ~MGeomCam(); + + Float_t GetCameraDist() const { return fCamDist; } + Float_t GetConvMm2Deg() const { return fConvMm2Deg; } + + UInt_t GetNumPixels() const { return fNumPixels; } + Float_t GetMaxRadius() const { return fMaxRadius; } + + MGeomPix &operator[](Int_t i) { return *(MGeomPix*)fPixels->UncheckedAt(i); } + MGeomPix &operator[](Int_t i) const { return *(MGeomPix*)fPixels->UncheckedAt(i); } + + virtual void Print(Option_t *opt=NULL) const; + + ClassDef(MGeomCam, 1) // Geometry base class for the camera +}; + +#endif + Index: trunk/MagicSoft/Mars/mgeom/MGeomCamCT1.cc =================================================================== --- trunk/MagicSoft/Mars/mgeom/MGeomCamCT1.cc (revision 1436) +++ trunk/MagicSoft/Mars/mgeom/MGeomCamCT1.cc (revision 1436) @@ -0,0 +1,267 @@ +/* ======================================================================== *\ +! +! * +! * This file is part of MARS, the MAGIC Analysis and Reconstruction +! * Software. It is distributed to you in the hope that it can be a useful +! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. +! * It is distributed WITHOUT ANY WARRANTY. +! * +! * Permission to use, copy, modify and distribute this software and its +! * documentation for any purpose is hereby granted without fee, +! * provided that the above copyright notice appear in all copies and +! * that both that copyright notice and this permission notice appear +! * in supporting documentation. It is provided "as is" without express +! * or implied warranty. +! * +! +! +! Author(s): Thomas Bretz 12/2000 +! Author(s): Harald Kornmayer 1/2001 +! +! Copyright: MAGIC Software Development, 2000-2002 +! +! +\* ======================================================================== */ + +///////////////////////////////////////////////////////////////////////////// +// +// MGeomCamCT1 +// +// This class stores the geometrz information of the CT1 camera. +// The next neighbor information comes from a table, the geometry layout +// is calculated (for Algorithm see CreateCam +// +//////////////////////////////////////////////////////////////////////////// + +#include "MGeomCamCT1.h" + + +#include // floor + +#include "TCanvas.h" + +#include "MLog.h" +#include "MLogManip.h" + +#include "MGeomPix.h" + +ClassImp(MGeomCamCT1); + +// -------------------------------------------------------------------------- +// +// CT1 camera has 127 pixels. For geometry and Next Neighbor info see +// CreateCam and CreateNN +// +MGeomCamCT1::MGeomCamCT1(const char *name) + : MGeomCam(127, 4.88, name, "Geometry information of CT1 camera") +{ + CreateCam(); + CreateNN(); + CalcMaxRadius(); +} + +// -------------------------------------------------------------------------- +// +// Create Next Neighbors: Fill the NN Info into the pixel objects. +// +void MGeomCamCT1::CreateNN() +{ + const Short_t nn[127][6] = { // Neighbors of # + { 1, 2, 3, 4, 5, 6}, // 0 + { 0, 2, 6, 7, 8, 18}, + { 0, 1, 3, 8, 9, 10}, + { 0, 2, 4, 10, 11, 12}, + { 0, 3, 5, 12, 13, 14}, + { 0, 4, 6, 14, 15, 16}, + { 0, 1, 5, 16, 17, 18}, + { 1, 8, 18, 19, 20, 36}, + { 1, 2, 7, 9, 20, 21}, + { 2, 8, 10, 21, 22, 23}, + { 2, 3, 9, 11, 23, 24}, // 10 + { 3, 10, 12, 24, 25, 26}, + { 3, 4, 11, 13, 26, 27}, + { 4, 12, 14, 27, 28, 29}, + { 4, 5, 13, 15, 29, 30}, + { 5, 14, 16, 30, 31, 32}, + { 5, 6, 15, 17, 32, 33}, + { 6, 16, 18, 33, 34, 35}, + { 1, 6, 7, 17, 35, 36}, + { 7, 20, 36, 37, 38, 60}, + { 7, 8, 19, 21, 38, 39}, // 20 + { 8, 9, 20, 22, 39, 40}, + { 9, 21, 23, 40, 41, 42}, + { 9, 10, 22, 24, 42, 43}, + { 10, 11, 23, 25, 43, 44}, + { 11, 24, 26, 44, 45, 46}, + { 11, 12, 25, 27, 46, 47}, + { 12, 13, 26, 28, 47, 48}, + { 13, 27, 29, 48, 49, 50}, + { 13, 14, 28, 30, 50, 51}, + { 14, 15, 29, 31, 51, 52}, // 30 + { 15, 30, 32, 52, 53, 54}, + { 15, 16, 31, 33, 54, 55}, + { 16, 17, 32, 34, 55, 56}, + { 17, 33, 35, 56, 57, 58}, + { 17, 18, 34, 36, 58, 59}, + { 7, 18, 19, 35, 59, 60}, + { 19, 38, 60, 61, 62, 90}, + { 19, 20, 37, 39, 62, 63}, + { 20, 21, 38, 40, 63, 64}, + { 21, 22, 39, 41, 64, 65}, // 40 + { 22, 40, 42, 65, 66, 67}, + { 22, 23, 41, 43, 67, 68}, + { 23, 24, 42, 44, 68, 69}, + { 24, 25, 43, 45, 69, 70}, + { 25, 44, 46, 70, 71, 72}, + { 25, 26, 45, 47, 72, 73}, + { 26, 27, 46, 48, 73, 74}, + { 27, 28, 47, 49, 74, 75}, + { 28, 48, 50, 75, 76, 77}, + { 28, 29, 49, 51, 77, 78}, // 50 + { 29, 30, 50, 52, 78, 79}, + { 30, 31, 51, 53, 79, 80}, + { 31, 52, 54, 80, 81, 82}, + { 31, 32, 53, 55, 82, 83}, + { 32, 33, 54, 56, 83, 84}, + { 33, 34, 55, 57, 84, 85}, + { 34, 56, 58, 85, 86, 87}, + { 34, 35, 57, 59, 87, 88}, + { 35, 36, 58, 60, 88, 89}, + { 19, 36, 37, 59, 89, 90}, // 60 + { 37, 62, 90, 91, 92, 126}, + { 37, 38, 61, 63, 92, 93}, + { 38, 39, 62, 64, 93, 94}, + { 39, 40, 63, 65, 94, 95}, + { 40, 41, 64, 66, 95, 96}, + { 41, 65, 67, 96, 97, 98}, + { 41, 42, 66, 68, 98, 99}, + { 42, 43, 67, 69, 99, 100}, + { 43, 44, 68, 70, 100, 101}, + { 44, 45, 69, 71, 101, 102}, // 70 + { 45, 70, 72, 102, 103, 104}, + { 45, 46, 71, 73, 104, 105}, + { 46, 47, 72, 74, 105, 106}, + { 47, 48, 73, 75, 106, 107}, + { 48, 49, 74, 76, 107, 108}, + { 49, 75, 77, 108, 109, 110}, + { 49, 50, 76, 78, 110, 111}, + { 50, 51, 77, 79, 111, 112}, + { 51, 52, 78, 80, 112, 113}, + { 52, 53, 79, 81, 113, 114}, // 80 + { 53, 80, 82, 114, 115, 116}, + { 53, 54, 81, 83, 116, 117}, + { 54, 55, 82, 84, 117, 118}, + { 55, 56, 83, 85, 118, 119}, + { 56, 57, 84, 86, 119, 120}, + { 57, 85, 87, 120, 121, 122}, + { 57, 58, 86, 88, 122, 123}, + { 58, 59, 87, 89, 123, 124}, + { 59, 60, 88, 90, 124, 125}, + { 37, 60, 61, 89, 125, 126}, // 90 + { 61, 92, 126, -1, -1, -1}, + { 61, 62, 91, 93, -1, -1}, + { 62, 63, 92, 94, -1, -1}, + { 63, 64, 93, 95, -1, -1}, + { 64, 65, 94, 96, -1, -1}, + { 65, 66, 95, 97, -1, -1}, + { 66, 96, 98, -1, -1, -1}, + { 66, 67, 97, 99, -1, -1}, + { 67, 68, 98, 100, -1, -1}, + { 68, 69, 99, 101, -1, -1}, // 100 + { 69, 70, 100, 102, -1, -1}, + { 70, 71, 101, 103, -1, -1}, + { 71, 102, 104, -1, -1, -1}, + { 71, 72, 103, 105, -1, -1}, + { 72, 73, 104, 106, -1, -1}, + { 73, 74, 105, 107, -1, -1}, + { 74, 75, 106, 108, -1, -1}, + { 75, 76, 107, 109, -1, -1}, + { 76, 108, 110, -1, -1, -1}, + { 76, 77, 109, 111, -1, -1}, // 110 + { 77, 78, 110, 112, -1, -1}, + { 78, 79, 111, 113, -1, -1}, + { 79, 80, 112, 114, -1, -1}, + { 80, 81, 113, 115, -1, -1}, + { 81, 114, 116, -1, -1, -1}, + { 81, 82, 115, 117, -1, -1}, + { 82, 83, 116, 118, -1, -1}, + { 83, 84, 117, 119, -1, -1}, + { 84, 85, 118, 120, -1, -1}, + { 85, 86, 119, 121, -1, -1}, // 120 + { 86, 120, 122, -1, -1, -1}, + { 86, 87, 121, 123, -1, -1}, + { 87, 88, 122, 124, -1, -1}, + { 88, 89, 123, 125, -1, -1}, + { 89, 90, 124, 126, -1, -1}, + { 61, 90, 91, 125, -1, -1} // 126 + }; + + for (Int_t i=0; i<127; i++) + (*this)[i].SetNeighbors(nn[i][0], nn[i][1], nn[i][2], + nn[i][3], nn[i][4], nn[i][5]); +} + +// -------------------------------------------------------------------------- +// +// Calculate the geometry information of CT1 and fill this information +// into the pixel objects. +// +void MGeomCamCT1::CreateCam() +{ + // + // fill the geometry class with the coordinates of the CT1 camera + // + *fLog << inf << " Create CT1 geometry " << endl; + + // + // this algorithm is from Martin Kestel originally + // it was punt into a root/C++ context by Harald Kornmayer and Thomas Bretz + + const Float_t diameter = 21; // units are mm + const Float_t kS32 = sqrt(3)/2; + + // + // add the first pixel to the list + // + Int_t pixnum = 0; + + (*this)[pixnum++].Set(0, 0, diameter); + + for (Int_t ring=1; ring<7; ring++) + { + // + // calc. coords for this ring counting from the + // starting number to the ending number + // + for (int i=0; i +! Author(s): Harald Kornmayer 1/2001 +! +! Copyright: MAGIC Software Development, 2000-2002 +! +! +\* ======================================================================== */ + +///////////////////////////////////////////////////////////////////////////// +// +// MGeomCamMagic +// +// This class stores the geometry information of the Magic camera. +// All information are copied from tables, see source code. +// +//////////////////////////////////////////////////////////////////////////// + +#include "MGeomCamMagic.h" + +#include "TCanvas.h" + +#include "MLog.h" +#include "MLogManip.h" + +#include "MGeomPix.h" + +ClassImp(MGeomCamMagic); + +// -------------------------------------------------------------------------- +// +// Magic camera has 577 pixels. For geometry and Next Neighbor info see +// CreateCam and CreateNN +// +MGeomCamMagic::MGeomCamMagic(const char *name) + : MGeomCam(577, 17, name, "Geometry information of Magic Camera") +{ + CreateCam(); + CreateNN(); + CalcMaxRadius(); +} + +// -------------------------------------------------------------------------- +// +// This fills the geometry information from a table into the pixel objects. +// +void MGeomCamMagic::CreateCam() +{ + // + // fill the geometry class with the coordinates of the MAGIC camera + // + *fLog << inf << " Creating Magic geometry " << endl ; + + // + // here define the hardwire things of the magic telescope + // + const Float_t xtemp[577] = { + +000.000, 30.000, 15.000, -15.000, -30.000, -15.000, 15.000, 60.000, // 0 + +045.000, 30.000, 0.000, -30.000, -45.000, -60.000, -45.000, -30.000, // 8 + +000.000, 30.000, 45.000, 90.000, 75.000, 60.000, 45.000, 15.000, // 16 + -015.000, -45.000, -60.000, -75.000, -90.000, -75.000, -60.000, -45.000, // 24 + -015.000, 15.000, 45.000, 60.000, 75.000, 120.000, 105.000, 90.000, // 32 + +075.000, 60.000, 30.000, 0.000, -30.000, -60.000, -75.000, -90.000, // 40 + -105.000, -120.000, -105.000, -90.000, -75.000, -60.000, -30.000, 0.000, // 48 + +030.000, 60.000, 75.000, 90.000, 105.000, 150.000, 135.000, 120.000, // 56 + +105.000, 90.000, 75.000, 45.000, 15.000, -15.000, -45.000, -75.000, // 64 + -090.000, -105.000, -120.000, -135.000, -150.000, -135.000, -120.000, -105.000, // 72 + -090.000, -75.000, -45.000, -15.000, 15.000, 45.000, 75.000, 90.000, // 80 + +105.000, 120.000, 135.000, 180.000, 165.000, 150.000, 135.000, 120.000, // 88 + +105.000, 90.000, 60.000, 30.000, 0.000, -30.000, -60.000, -90.000, // 96 + -105.000, -120.000, -135.000, -150.000, -165.000, -180.000, -165.000, -150.000, // 104 + -135.000, -120.000, -105.000, -90.000, -60.000, -30.000, 0.000, 30.000, // 112 + +060.000, 90.000, 105.000, 120.000, 135.000, 150.000, 165.000, 210.000, // 120 + +195.000, 180.000, 165.000, 150.000, 135.000, 120.000, 105.000, 75.000, // 128 + +045.000, 15.000, -15.000, -45.000, -75.000, -105.000, -120.000, -135.000, // 136 + -150.000, -165.000, -180.000, -195.000, -210.000, -195.000, -180.000, -165.000, // 144 + -150.000, -135.000, -120.000, -105.000, -75.000, -45.000, -15.000, 15.000, // 152 + +045.000, 75.000, 105.000, 120.000, 135.000, 150.000, 165.000, 180.000, // 160 + +195.000, 240.000, 225.000, 210.000, 195.000, 180.000, 165.000, 150.000, // 168 + +135.000, 120.000, 90.000, 60.000, 30.000, 0.000, -30.000, -60.000, // 176 + -090.000, -120.000, -135.000, -150.000, -165.000, -180.000, -195.000, -210.000, // 184 + -225.000, -240.000, -225.000, -210.000, -195.000, -180.000, -165.000, -150.000, // 192 + -135.000, -120.000, -90.000, -60.000, -30.000, 0.000, 30.000, 60.000, // 200 + +090.000, 120.000, 135.000, 150.000, 165.000, 180.000, 195.000, 210.000, // 208 + +225.000, 270.000, 255.000, 240.000, 225.000, 210.000, 195.000, 180.000, // 216 + +165.000, 150.000, 135.000, 105.000, 75.000, 45.000, 15.000, -15.000, // 224 + -045.000, -75.000, -105.000, -135.000, -150.000, -165.000, -180.000, -195.000, // 232 + -210.000, -225.000, -240.000, -255.000, -270.000, -255.000, -240.000, -225.000, // 240 + -210.000, -195.000, -180.000, -165.000, -150.000, -135.000, -105.000, -75.000, // 248 + -045.000, -15.000, 15.000, 45.000, 75.000, 105.000, 135.000, 150.000, // 256 + +165.000, 180.000, 195.000, 210.000, 225.000, 240.000, 255.000, 300.000, // 264 + +285.000, 270.000, 255.000, 240.000, 225.000, 210.000, 195.000, 180.000, // 272 + +165.000, 150.000, 120.000, 90.000, 60.000, 30.000, 0.000, -30.000, // 280 + -060.000, -90.000, -120.000, -150.000, -165.000, -180.000, -195.000, -210.000, // 288 + -225.000, -240.000, -255.000, -270.000, -285.000, -300.000, -285.000, -270.000, // 296 + -255.000, -240.000, -225.000, -210.000, -195.000, -180.000, -165.000, -150.000, // 304 + -120.000, -90.000, -60.000, -30.000, 0.000, 30.000, 60.000, 90.000, // 312 + +120.000, 150.000, 165.000, 180.000, 195.000, 210.000, 225.000, 240.000, // 320 + +255.000, 270.000, 285.000, 330.000, 315.000, 300.000, 285.000, 270.000, // 328 + +255.000, 240.000, 225.000, 210.000, 195.000, 180.000, 165.000, 135.000, // 336 + +105.000, 75.000, 45.000, 15.000, -15.000, -45.000, -75.000, -105.000, // 344 + -135.000, -165.000, -180.000, -195.000, -210.000, -225.000, -240.000, -255.000, // 352 + -270.000, -285.000, -300.000, -315.000, -330.000, -315.000, -300.000, -285.000, // 360 + -270.000, -255.000, -240.000, -225.000, -210.000, -195.000, -180.000, -165.000, // 368 + -135.000, -105.000, -75.000, -45.000, -15.000, 15.000, 45.000, 75.000, // 376 + +105.000, 135.000, 165.000, 180.000, 195.000, 210.000, 225.000, 240.000, // 384 + +255.000, 270.000, 285.000, 300.000, 315.000, 360.000, 330.000, 300.000, // 392 + +270.000, 240.000, 210.000, 150.000, 90.000, 30.000, -30.000, -90.000, // 400 + -150.000, -210.000, -240.000, -270.000, -300.000, -330.000, -360.000, -360.000, // 408 + -330.000, -300.000, -270.000, -240.000, -210.000, -150.000, -90.000, -30.000, // 416 + +030.000, 90.000, 150.000, 210.000, 240.000, 270.000, 300.000, 330.000, // 424 + +360.000, 420.000, 390.000, 360.000, 330.000, 300.000, 270.000, 240.000, // 432 + +180.000, 120.000, 60.000, 0.000, -60.000, -120.000, -180.000, -240.000, // 440 + -270.000, -300.000, -330.000, -360.000, -390.000, -420.000, -420.000, -390.000, // 448 + -360.000, -330.000, -300.000, -270.000, -240.000, -180.000, -120.000, -60.000, // 456 + +000.000, 60.000, 120.000, 180.000, 240.000, 270.000, 300.000, 330.000, // 464 + +360.000, 390.000, 420.000, 480.000, 450.000, 420.000, 390.000, 360.000, // 472 + +330.000, 300.000, 270.000, 210.000, 150.000, 90.000, 30.000, -30.000, // 480 + -090.000, -150.000, -210.000, -270.000, -300.000, -330.000, -360.000, -390.000, // 488 + -420.000, -450.000, -480.000, -480.000, -450.000, -420.000, -390.000, -360.000, // 496 + -330.000, -300.000, -270.000, -210.000, -150.000, -90.000, -30.000, 30.000, // 504 + +090.000, 150.000, 210.000, 270.000, 300.000, 330.000, 360.000, 390.000, // 512 + +420.000, 450.000, 480.000, 540.000, 510.000, 480.000, 450.000, 420.000, // 520 + +390.000, 360.000, 330.000, 300.000, 240.000, 180.000, 120.000, 60.000, // 528 + +000.000, -60.000, -120.000, -180.000, -240.000, -300.000, -330.000, -360.000, // 536 + -390.000, -420.000, -450.000, -480.000, -510.000, -540.000, -540.000, -510.000, // 544 + -480.000, -450.000, -420.000, -390.000, -360.000, -330.000, -300.000, -240.000, // 552 + -180.000, -120.000, -60.000, 0.000, 60.000, 120.000, 180.000, 240.000, // 560 + +300.000, 330.000, 360.000, 390.000, 420.000, 450.000, 480.000, 510.000, // 568 + +540.000 // 576 + }; + + const Float_t ytemp[577] = { + +000.000, 0.000, 25.981, 25.981, 0.000, -25.981, -25.981, 0.000, // 0 + +025.981, 51.961, 51.961, 51.961, 25.981, 0.000, -25.981, -51.961, // 8 + -051.961, -51.961, -25.981, 0.000, 25.981, 51.961, 77.942, 77.942, // 16 + +077.942, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 24 + -077.942, -77.942, -77.942, -51.961, -25.981, 0.000, 25.981, 51.961, // 32 + +077.942, 103.923, 103.923, 103.923, 103.923, 103.923, 77.942, 51.961, // 40 + +025.981, 0.000, -25.981, -51.961, -77.942, -103.923, -103.923, -103.923, // 48 + -103.923, -103.923, -77.942, -51.961, -25.981, 0.000, 25.981, 51.961, // 56 + +077.942, 103.923, 129.904, 129.904, 129.904, 129.904, 129.904, 129.904, // 64 + +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 72 + -103.923, -129.904, -129.904, -129.904, -129.904, -129.904, -129.904, -103.923, // 80 + -077.942, -51.961, -25.981, 0.000, 25.981, 51.961, 77.942, 103.923, // 88 + +129.904, 155.885, 155.885, 155.885, 155.885, 155.885, 155.885, 155.885, // 96 + +129.904, 103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, // 104 + -077.942, -103.923, -129.904, -155.885, -155.885, -155.885, -155.885, -155.885, // 112 + -155.885, -155.885, -129.904, -103.923, -77.942, -51.961, -25.981, 0.000, // 120 + +025.981, 51.961, 77.942, 103.923, 129.904, 155.885, 181.865, 181.865, // 128 + +181.865, 181.865, 181.865, 181.865, 181.865, 181.865, 155.885, 129.904, // 136 + +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 144 + -103.923, -129.904, -155.885, -181.865, -181.865, -181.865, -181.865, -181.865, // 152 + -181.865, -181.865, -181.865, -155.885, -129.904, -103.923, -77.942, -51.961, // 160 + -025.981, 0.000, 25.981, 51.961, 77.942, 103.923, 129.904, 155.885, // 168 + +181.865, 207.846, 207.846, 207.846, 207.846, 207.846, 207.846, 207.846, // 176 + +207.846, 207.846, 181.865, 155.885, 129.904, 103.923, 77.942, 51.961, // 184 + +025.981, 0.000, -25.981, -51.961, -77.942, -103.923, -129.904, -155.885, // 192 + -181.865, -207.846, -207.846, -207.846, -207.846, -207.846, -207.846, -207.846, // 200 + -207.846, -207.846, -181.865, -155.885, -129.904, -103.923, -77.942, -51.961, // 208 + -025.981, 0.000, 25.981, 51.961, 77.942, 103.923, 129.904, 155.885, // 216 + +181.865, 207.846, 233.827, 233.827, 233.827, 233.827, 233.827, 233.827, // 224 + +233.827, 233.827, 233.827, 233.827, 207.846, 181.865, 155.885, 129.904, // 232 + +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 240 + -103.923, -129.904, -155.885, -181.865, -207.846, -233.827, -233.827, -233.827, // 248 + -233.827, -233.827, -233.827, -233.827, -233.827, -233.827, -233.827, -207.846, // 256 + -181.865, -155.885, -129.904, -103.923, -77.942, -51.961, -25.981, 0.000, // 264 + +025.981, 51.961, 77.942, 103.923, 129.904, 155.885, 181.865, 207.846, // 272 + +233.827, 259.808, 259.808, 259.808, 259.808, 259.808, 259.808, 259.808, // 280 + +259.808, 259.808, 259.808, 259.808, 233.827, 207.846, 181.865, 155.885, // 288 + +129.904, 103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, // 296 + -077.942, -103.923, -129.904, -155.885, -181.865, -207.846, -233.827, -259.808, // 304 + -259.808, -259.808, -259.808, -259.808, -259.808, -259.808, -259.808, -259.808, // 312 + -259.808, -259.808, -233.827, -207.846, -181.865, -155.885, -129.904, -103.923, // 320 + -077.942, -51.961, -25.981, 0.000, 25.981, 51.961, 77.942, 103.923, // 328 + +129.904, 155.885, 181.865, 207.846, 233.827, 259.808, 285.788, 285.788, // 336 + +285.788, 285.788, 285.788, 285.788, 285.788, 285.788, 285.788, 285.788, // 344 + +285.788, 285.788, 259.808, 233.827, 207.846, 181.865, 155.885, 129.904, // 352 + +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 360 + -103.923, -129.904, -155.885, -181.865, -207.846, -233.827, -259.808, -285.788, // 368 + -285.788, -285.788, -285.788, -285.788, -285.788, -285.788, -285.788, -285.788, // 376 + -285.788, -285.788, -285.788, -259.808, -233.827, -207.846, -181.865, -155.885, // 384 + -129.904, -103.923, -77.942, -51.961, -25.981, 34.641, 86.603, 138.564, // 392 + +190.526, 242.487, 294.449, 329.090, 329.090, 329.090, 329.090, 329.090, // 400 + +329.090, 294.449, 242.487, 190.526, 138.564, 86.603, 34.641, -34.641, // 408 + -086.603, -138.564, -190.526, -242.487, -294.449, -329.090, -329.090, -329.090, // 416 + -329.090, -329.090, -329.090, -294.449, -242.487, -190.526, -138.564, -86.603, // 424 + -034.641, 34.641, 86.603, 138.564, 190.526, 242.487, 294.449, 346.410, // 432 + +381.051, 381.051, 381.051, 381.051, 381.051, 381.051, 381.051, 346.410, // 440 + +294.449, 242.487, 190.526, 138.564, 86.603, 34.641, -34.641, -86.603, // 448 + -138.564, -190.526, -242.487, -294.449, -346.410, -381.051, -381.051, -381.051, // 456 + -381.051, -381.051, -381.051, -381.051, -346.410, -294.449, -242.487, -190.526, // 464 + -138.564, -86.603, -34.641, 34.641, 86.603, 138.564, 190.526, 242.487, // 472 + +294.449, 346.410, 398.372, 433.013, 433.013, 433.013, 433.013, 433.013, // 480 + +433.013, 433.013, 433.013, 398.372, 346.410, 294.449, 242.487, 190.526, // 488 + +138.564, 86.603, 34.641, -34.641, -86.603, -138.564, -190.526, -242.487, // 496 + -294.449, -346.410, -398.372, -433.013, -433.013, -433.013, -433.013, -433.013, // 504 + -433.013, -433.013, -433.013, -398.372, -346.410, -294.449, -242.487, -190.526, // 512 + -138.564, -86.603, -34.641, 34.641, 86.603, 138.564, 190.526, 242.487, // 520 + +294.449, 346.410, 398.372, 450.333, 484.974, 484.974, 484.974, 484.974, // 528 + +484.974, 484.974, 484.974, 484.974, 484.974, 450.333, 398.372, 346.410, // 536 + +294.449, 242.487, 190.526, 138.564, 86.603, 34.641, -34.641, -86.603, // 544 + -138.564, -190.526, -242.487, -294.449, -346.410, -398.372, -450.333, -484.974, // 552 + -484.974, -484.974, -484.974, -484.974, -484.974, -484.974, -484.974, -484.974, // 560 + -450.333, -398.372, -346.410, -294.449, -242.487, -190.526, -138.564, -86.603, // 568 + -034.641 // 576 + }; + + const Float_t rtemp[577] = { + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 0 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 8 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 16 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 24 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 32 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 40 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 48 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 56 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 64 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 72 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 80 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 88 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 96 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 104 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 112 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 120 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 128 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 136 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 144 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 152 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 160 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 168 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 176 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 184 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 192 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 200 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 208 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 216 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 224 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 232 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 240 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 248 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 256 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 264 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 272 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 280 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 288 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 296 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 304 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 312 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 320 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 328 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 336 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 344 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 352 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 360 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 368 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 376 + 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 384 + 30.00, 30.00, 30.00, 30.00, 30.00, 60.00, 60.00, 60.00, // 392 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 400 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 408 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 416 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 424 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 432 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 440 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 448 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 456 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 464 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 472 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 480 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 488 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 496 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 504 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 512 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 520 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 528 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 536 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 544 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 552 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 560 + 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 568 + 60.00 }; // 576 + + // + // fill the pixels list with this data + // + + for (UInt_t i=0; i +! Author(s): Harald Kornmayer 1/2001 +! +! Copyright: MAGIC Software Development, 2000-2002 +! +! +\* ======================================================================== */ + +////////////////////////////////////////////////////////////////////////////// +// +// MGeomPix +// +// This container stores the geometry (position) information of +// a single pixel together with the information about next neighbors. +// +//////////////////////////////////////////////////////////////////////////// + +#include "MGeomPix.h" + +#include + +#include "MLog.h" +#include "MLogManip.h" + +ClassImp(MGeomPix); + +// -------------------------------------------------------------------------- +// +// Initialiyes one pixel +// +MGeomPix::MGeomPix(Float_t x, Float_t y, Float_t r) : fX(x), fY(y), fR(r) +{ + // default constructor +} + +// -------------------------------------------------------------------------- +// +// Return the area of the pixel. A hegagonal shape is assumed. +// +Float_t MGeomPix::GetA() const +{ + return fR*fR*tan(60/kRad2Deg); +} + +// -------------------------------------------------------------------------- +// +// Initialiyes Next Neighbors. +// +// WARNING: This function is public, but it is not ment for user access. +// It should only be used from geometry classes (like MGeomCam) +// +void MGeomPix::SetNeighbors(Short_t i0, Short_t i1, Short_t i2, + Short_t i3, Short_t i4, Short_t i5) +{ + fNeighbors[0] = i0; + fNeighbors[1] = i1; + fNeighbors[2] = i2; + fNeighbors[3] = i3; + fNeighbors[4] = i4; + fNeighbors[5] = i5; + + int i; + for (i=0; i<6; i++) + if (fNeighbors[i]<0) + break; + + fNumNeighbors = i; +} + +// -------------------------------------------------------------------------- +// +// Print the geometry information of one pixel. +// +void MGeomPix::Print(Option_t *opt) const +{ + // information about a pixel + *fLog << all << "MPixGeom: x= " << fX + << " y= " << fY + << " r= " << fR + << endl ; +} + Index: trunk/MagicSoft/Mars/mgeom/MGeomPix.h =================================================================== --- trunk/MagicSoft/Mars/mgeom/MGeomPix.h (revision 1436) +++ trunk/MagicSoft/Mars/mgeom/MGeomPix.h (revision 1436) @@ -0,0 +1,46 @@ +#ifndef MARS_MGeomPix +#define MARS_MGeomPix + +#ifndef MARS_MParContainer +#include "MParContainer.h" +#endif + +class MGeomPix : public MParContainer +{ +private: + Float_t fX; // [mm] the x coordinate of the center + Float_t fY; // [mm] the y coordinate of the center + Float_t fR; // [mm] the r coordinate of the pixel (dist between two parallel sides) + + Byte_t fNumNeighbors; // number of valid neighbors + Short_t fNeighbors[6]; // the IDs of the pixel next to it + // we are assuming an hexagonal geometry + +public: + MGeomPix(Float_t x=0, Float_t y=0, Float_t r=0); + + void Print(Option_t *opt=NULL) const; + + void Set(Float_t x, Float_t y, Float_t r) { fX=x; fY=y; fR=r; } + + void SetNeighbors(Short_t i0=-1, Short_t i1=-1, Short_t i2=-1, + Short_t i3=-1, Short_t i4=-1, Short_t i5=-1); + + void SetX(Float_t x) { fX = x; } + void SetY(Float_t y) { fY = y; } + void SetR(Float_t r) { fR = r; } + + Float_t GetX() const { return fX; } + Float_t GetY() const { return fY; } + Float_t GetR() const { return fR; } + + Float_t GetA() const; + + Byte_t GetNumNeighbors() const { return fNumNeighbors; } + Short_t GetNeighbor(Byte_t i) const { return fNeighbors[i]; } + + ClassDef(MGeomPix, 1) // Geometry class describing the geometry of one pixel +}; + +#endif + Index: trunk/MagicSoft/Mars/mgeom/Makefile =================================================================== --- trunk/MagicSoft/Mars/mgeom/Makefile (revision 1436) +++ trunk/MagicSoft/Mars/mgeom/Makefile (revision 1436) @@ -0,0 +1,50 @@ +################################################################## +# +# makefile +# +# for the MARS software +# +################################################################## +include ../Makefile.conf.$(OSTYPE) +include ../Makefile.conf.general + +# +# Handling name of the Root Dictionary Files +# +CINT = Geom + +# +# Library name to creatre +# +LIB = mgeom.a + +# +# connect the include files defined in the config.mk file +# +INCLUDES = -I. -I../mbase + +#------------------------------------------------------------------------------ + +.SUFFIXES: .c .cc .cxx .h .hxx .o + +SRCFILES = MGeomPix.cc \ + MGeomCam.cc \ + MGeomCamCT1.cc \ + MGeomCamMagic.cc + +SRCS = $(SRCFILES) +HEADERS = $(SRCFILES:.cc=.h) +OBJS = $(SRCFILES:.cc=.o) + +############################################################ + +all: $(LIB) + +include ../Makefile.rules + +clean: rmcint rmobjs rmcore rmlib + +mrproper: clean rmbak + +# @endcode + Index: trunk/MagicSoft/Mars/mgui/GuiLinkDef.h =================================================================== --- trunk/MagicSoft/Mars/mgui/GuiLinkDef.h (revision 1435) +++ trunk/MagicSoft/Mars/mgui/GuiLinkDef.h (revision 1436) @@ -6,12 +6,5 @@ #pragma link C++ class MHexagon+; - -#pragma link C++ class MGeomPix+; -#pragma link C++ class MGeomCam+; -#pragma link C++ class MGeomCamCT1+; -#pragma link C++ class MGeomCamMagic+; - #pragma link C++ class MCamDisplay+; -#pragma link C++ class MineSweeper+; #endif Index: trunk/MagicSoft/Mars/mgui/MGeomCam.cc =================================================================== --- trunk/MagicSoft/Mars/mgui/MGeomCam.cc (revision 1435) +++ (revision ) @@ -1,119 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz 12/2000 -! Author(s): Harald Kornmayer 1/2001 -! -! Copyright: MAGIC Software Development, 2000-2002 -! -! -\* ======================================================================== */ - -/////////////////////////////////////////////////////////////////////// -// -// MGeomCam -// -// This is the base class of different camera geometries. It creates -// a pixel object for a given number of pixels and defines the -// interface of how to acccess the geometry information. -// -/////////////////////////////////////////////////////////////////////// - -#include "MGeomCam.h" - -#include "MLog.h" -#include "MLogManip.h" - -#include "MGeomPix.h" -#include "MHexagon.h" - -ClassImp(MGeomCam); - -// -------------------------------------------------------------------------- -// -// Initializes a Camera Geometry with npix pixels. All pixels -// are deleted when the corresponding array is deleted. -// -MGeomCam::MGeomCam(UInt_t npix, Float_t dist, const char *name, const char *title) - : fNumPixels(npix), fCamDist(dist), fConvMm2Deg(kRad2Deg/(dist*1000)) -{ - fName = name ? name : "MGeomCam"; - fTitle = title ? title : "Storage container for a camera geometry"; - - fPixels = new TObjArray(npix); - - // - // make sure that the destructor delete all contained objects - // - fPixels->SetOwner(); - - for (UInt_t i=0; ifMaxRadius) - fMaxRadius = maxr; - } -} - -// -------------------------------------------------------------------------- -// -// Prints the Geometry information of all pixels in the camera -// -void MGeomCam::Print(Option_t *) const -{ - // - // Print Information about the Geometry of the camera - // - *fLog << all << " Number of Pixels (" << GetTitle() << "): " << fNumPixels << endl; - - for (UInt_t i=0; i -#endif - -class MGeomPix; - -class MGeomCam : public MParContainer -{ -private: - UInt_t fNumPixels; // Number of pixels in this camera - Float_t fMaxRadius; // maximum radius of the camera (eg. for GUI layout) - - Float_t fCamDist; // [m] Average distance of the camera from the mirror - Float_t fConvMm2Deg; // conversion factor to convert mm in the camera plain into degrees - - TObjArray *fPixels; // Array of singel pixels storing the geometry - -protected: - void CalcMaxRadius(); - -public: - - MGeomCam(UInt_t npix, Float_t dist, const char *name=NULL, const char *title=NULL); - - virtual ~MGeomCam(); - - Float_t GetCameraDist() const { return fCamDist; } - Float_t GetConvMm2Deg() const { return fConvMm2Deg; } - - UInt_t GetNumPixels() const { return fNumPixels; } - Float_t GetMaxRadius() const { return fMaxRadius; } - - MGeomPix &operator[](Int_t i) { return *(MGeomPix*)fPixels->UncheckedAt(i); } - MGeomPix &operator[](Int_t i) const { return *(MGeomPix*)fPixels->UncheckedAt(i); } - - virtual void Print(Option_t *opt=NULL) const; - - ClassDef(MGeomCam, 1) // Geometry base class for the camera -}; - -#endif - Index: trunk/MagicSoft/Mars/mgui/MGeomCamCT1.cc =================================================================== --- trunk/MagicSoft/Mars/mgui/MGeomCamCT1.cc (revision 1435) +++ (revision ) @@ -1,267 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * This file is part of MARS, the MAGIC Analysis and Reconstruction -! * Software. It is distributed to you in the hope that it can be a useful -! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. -! * It is distributed WITHOUT ANY WARRANTY. -! * -! * Permission to use, copy, modify and distribute this software and its -! * documentation for any purpose is hereby granted without fee, -! * provided that the above copyright notice appear in all copies and -! * that both that copyright notice and this permission notice appear -! * in supporting documentation. It is provided "as is" without express -! * or implied warranty. -! * -! -! -! Author(s): Thomas Bretz 12/2000 -! Author(s): Harald Kornmayer 1/2001 -! -! Copyright: MAGIC Software Development, 2000-2002 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// MGeomCamCT1 -// -// This class stores the geometrz information of the CT1 camera. -// The next neighbor information comes from a table, the geometry layout -// is calculated (for Algorithm see CreateCam -// -//////////////////////////////////////////////////////////////////////////// - -#include "MGeomCamCT1.h" - - -#include // floor - -#include "TCanvas.h" - -#include "MLog.h" -#include "MLogManip.h" - -#include "MGeomPix.h" - -ClassImp(MGeomCamCT1); - -// -------------------------------------------------------------------------- -// -// CT1 camera has 127 pixels. For geometry and Next Neighbor info see -// CreateCam and CreateNN -// -MGeomCamCT1::MGeomCamCT1(const char *name) - : MGeomCam(127, 4.88, name, "Geometry information of CT1 camera") -{ - CreateCam(); - CreateNN(); - CalcMaxRadius(); -} - -// -------------------------------------------------------------------------- -// -// Create Next Neighbors: Fill the NN Info into the pixel objects. -// -void MGeomCamCT1::CreateNN() -{ - const Short_t nn[127][6] = { // Neighbors of # - { 1, 2, 3, 4, 5, 6}, // 0 - { 0, 2, 6, 7, 8, 18}, - { 0, 1, 3, 8, 9, 10}, - { 0, 2, 4, 10, 11, 12}, - { 0, 3, 5, 12, 13, 14}, - { 0, 4, 6, 14, 15, 16}, - { 0, 1, 5, 16, 17, 18}, - { 1, 8, 18, 19, 20, 36}, - { 1, 2, 7, 9, 20, 21}, - { 2, 8, 10, 21, 22, 23}, - { 2, 3, 9, 11, 23, 24}, // 10 - { 3, 10, 12, 24, 25, 26}, - { 3, 4, 11, 13, 26, 27}, - { 4, 12, 14, 27, 28, 29}, - { 4, 5, 13, 15, 29, 30}, - { 5, 14, 16, 30, 31, 32}, - { 5, 6, 15, 17, 32, 33}, - { 6, 16, 18, 33, 34, 35}, - { 1, 6, 7, 17, 35, 36}, - { 7, 20, 36, 37, 38, 60}, - { 7, 8, 19, 21, 38, 39}, // 20 - { 8, 9, 20, 22, 39, 40}, - { 9, 21, 23, 40, 41, 42}, - { 9, 10, 22, 24, 42, 43}, - { 10, 11, 23, 25, 43, 44}, - { 11, 24, 26, 44, 45, 46}, - { 11, 12, 25, 27, 46, 47}, - { 12, 13, 26, 28, 47, 48}, - { 13, 27, 29, 48, 49, 50}, - { 13, 14, 28, 30, 50, 51}, - { 14, 15, 29, 31, 51, 52}, // 30 - { 15, 30, 32, 52, 53, 54}, - { 15, 16, 31, 33, 54, 55}, - { 16, 17, 32, 34, 55, 56}, - { 17, 33, 35, 56, 57, 58}, - { 17, 18, 34, 36, 58, 59}, - { 7, 18, 19, 35, 59, 60}, - { 19, 38, 60, 61, 62, 90}, - { 19, 20, 37, 39, 62, 63}, - { 20, 21, 38, 40, 63, 64}, - { 21, 22, 39, 41, 64, 65}, // 40 - { 22, 40, 42, 65, 66, 67}, - { 22, 23, 41, 43, 67, 68}, - { 23, 24, 42, 44, 68, 69}, - { 24, 25, 43, 45, 69, 70}, - { 25, 44, 46, 70, 71, 72}, - { 25, 26, 45, 47, 72, 73}, - { 26, 27, 46, 48, 73, 74}, - { 27, 28, 47, 49, 74, 75}, - { 28, 48, 50, 75, 76, 77}, - { 28, 29, 49, 51, 77, 78}, // 50 - { 29, 30, 50, 52, 78, 79}, - { 30, 31, 51, 53, 79, 80}, - { 31, 52, 54, 80, 81, 82}, - { 31, 32, 53, 55, 82, 83}, - { 32, 33, 54, 56, 83, 84}, - { 33, 34, 55, 57, 84, 85}, - { 34, 56, 58, 85, 86, 87}, - { 34, 35, 57, 59, 87, 88}, - { 35, 36, 58, 60, 88, 89}, - { 19, 36, 37, 59, 89, 90}, // 60 - { 37, 62, 90, 91, 92, 126}, - { 37, 38, 61, 63, 92, 93}, - { 38, 39, 62, 64, 93, 94}, - { 39, 40, 63, 65, 94, 95}, - { 40, 41, 64, 66, 95, 96}, - { 41, 65, 67, 96, 97, 98}, - { 41, 42, 66, 68, 98, 99}, - { 42, 43, 67, 69, 99, 100}, - { 43, 44, 68, 70, 100, 101}, - { 44, 45, 69, 71, 101, 102}, // 70 - { 45, 70, 72, 102, 103, 104}, - { 45, 46, 71, 73, 104, 105}, - { 46, 47, 72, 74, 105, 106}, - { 47, 48, 73, 75, 106, 107}, - { 48, 49, 74, 76, 107, 108}, - { 49, 75, 77, 108, 109, 110}, - { 49, 50, 76, 78, 110, 111}, - { 50, 51, 77, 79, 111, 112}, - { 51, 52, 78, 80, 112, 113}, - { 52, 53, 79, 81, 113, 114}, // 80 - { 53, 80, 82, 114, 115, 116}, - { 53, 54, 81, 83, 116, 117}, - { 54, 55, 82, 84, 117, 118}, - { 55, 56, 83, 85, 118, 119}, - { 56, 57, 84, 86, 119, 120}, - { 57, 85, 87, 120, 121, 122}, - { 57, 58, 86, 88, 122, 123}, - { 58, 59, 87, 89, 123, 124}, - { 59, 60, 88, 90, 124, 125}, - { 37, 60, 61, 89, 125, 126}, // 90 - { 61, 92, 126, -1, -1, -1}, - { 61, 62, 91, 93, -1, -1}, - { 62, 63, 92, 94, -1, -1}, - { 63, 64, 93, 95, -1, -1}, - { 64, 65, 94, 96, -1, -1}, - { 65, 66, 95, 97, -1, -1}, - { 66, 96, 98, -1, -1, -1}, - { 66, 67, 97, 99, -1, -1}, - { 67, 68, 98, 100, -1, -1}, - { 68, 69, 99, 101, -1, -1}, // 100 - { 69, 70, 100, 102, -1, -1}, - { 70, 71, 101, 103, -1, -1}, - { 71, 102, 104, -1, -1, -1}, - { 71, 72, 103, 105, -1, -1}, - { 72, 73, 104, 106, -1, -1}, - { 73, 74, 105, 107, -1, -1}, - { 74, 75, 106, 108, -1, -1}, - { 75, 76, 107, 109, -1, -1}, - { 76, 108, 110, -1, -1, -1}, - { 76, 77, 109, 111, -1, -1}, // 110 - { 77, 78, 110, 112, -1, -1}, - { 78, 79, 111, 113, -1, -1}, - { 79, 80, 112, 114, -1, -1}, - { 80, 81, 113, 115, -1, -1}, - { 81, 114, 116, -1, -1, -1}, - { 81, 82, 115, 117, -1, -1}, - { 82, 83, 116, 118, -1, -1}, - { 83, 84, 117, 119, -1, -1}, - { 84, 85, 118, 120, -1, -1}, - { 85, 86, 119, 121, -1, -1}, // 120 - { 86, 120, 122, -1, -1, -1}, - { 86, 87, 121, 123, -1, -1}, - { 87, 88, 122, 124, -1, -1}, - { 88, 89, 123, 125, -1, -1}, - { 89, 90, 124, 126, -1, -1}, - { 61, 90, 91, 125, -1, -1} // 126 - }; - - for (Int_t i=0; i<127; i++) - (*this)[i].SetNeighbors(nn[i][0], nn[i][1], nn[i][2], - nn[i][3], nn[i][4], nn[i][5]); -} - -// -------------------------------------------------------------------------- -// -// Calculate the geometry information of CT1 and fill this information -// into the pixel objects. -// -void MGeomCamCT1::CreateCam() -{ - // - // fill the geometry class with the coordinates of the CT1 camera - // - *fLog << inf << " Create CT1 geometry " << endl; - - // - // this algorithm is from Martin Kestel originally - // it was punt into a root/C++ context by Harald Kornmayer and Thomas Bretz - - const Float_t diameter = 21; // units are mm - const Float_t kS32 = sqrt(3)/2; - - // - // add the first pixel to the list - // - Int_t pixnum = 0; - - (*this)[pixnum++].Set(0, 0, diameter); - - for (Int_t ring=1; ring<7; ring++) - { - // - // calc. coords for this ring counting from the - // starting number to the ending number - // - for (int i=0; i -! Author(s): Harald Kornmayer 1/2001 -! -! Copyright: MAGIC Software Development, 2000-2002 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// MGeomCamMagic -// -// This class stores the geometry information of the Magic camera. -// All information are copied from tables, see source code. -// -//////////////////////////////////////////////////////////////////////////// - -#include "MGeomCamMagic.h" - -#include "TCanvas.h" - -#include "MLog.h" -#include "MLogManip.h" - -#include "MGeomPix.h" - -ClassImp(MGeomCamMagic); - -// -------------------------------------------------------------------------- -// -// Magic camera has 577 pixels. For geometry and Next Neighbor info see -// CreateCam and CreateNN -// -MGeomCamMagic::MGeomCamMagic(const char *name) - : MGeomCam(577, 17, name, "Geometry information of Magic Camera") -{ - CreateCam(); - CreateNN(); - CalcMaxRadius(); -} - -// -------------------------------------------------------------------------- -// -// This fills the geometry information from a table into the pixel objects. -// -void MGeomCamMagic::CreateCam() -{ - // - // fill the geometry class with the coordinates of the MAGIC camera - // - *fLog << inf << " Creating Magic geometry " << endl ; - - // - // here define the hardwire things of the magic telescope - // - const Float_t xtemp[577] = { - +000.000, 30.000, 15.000, -15.000, -30.000, -15.000, 15.000, 60.000, // 0 - +045.000, 30.000, 0.000, -30.000, -45.000, -60.000, -45.000, -30.000, // 8 - +000.000, 30.000, 45.000, 90.000, 75.000, 60.000, 45.000, 15.000, // 16 - -015.000, -45.000, -60.000, -75.000, -90.000, -75.000, -60.000, -45.000, // 24 - -015.000, 15.000, 45.000, 60.000, 75.000, 120.000, 105.000, 90.000, // 32 - +075.000, 60.000, 30.000, 0.000, -30.000, -60.000, -75.000, -90.000, // 40 - -105.000, -120.000, -105.000, -90.000, -75.000, -60.000, -30.000, 0.000, // 48 - +030.000, 60.000, 75.000, 90.000, 105.000, 150.000, 135.000, 120.000, // 56 - +105.000, 90.000, 75.000, 45.000, 15.000, -15.000, -45.000, -75.000, // 64 - -090.000, -105.000, -120.000, -135.000, -150.000, -135.000, -120.000, -105.000, // 72 - -090.000, -75.000, -45.000, -15.000, 15.000, 45.000, 75.000, 90.000, // 80 - +105.000, 120.000, 135.000, 180.000, 165.000, 150.000, 135.000, 120.000, // 88 - +105.000, 90.000, 60.000, 30.000, 0.000, -30.000, -60.000, -90.000, // 96 - -105.000, -120.000, -135.000, -150.000, -165.000, -180.000, -165.000, -150.000, // 104 - -135.000, -120.000, -105.000, -90.000, -60.000, -30.000, 0.000, 30.000, // 112 - +060.000, 90.000, 105.000, 120.000, 135.000, 150.000, 165.000, 210.000, // 120 - +195.000, 180.000, 165.000, 150.000, 135.000, 120.000, 105.000, 75.000, // 128 - +045.000, 15.000, -15.000, -45.000, -75.000, -105.000, -120.000, -135.000, // 136 - -150.000, -165.000, -180.000, -195.000, -210.000, -195.000, -180.000, -165.000, // 144 - -150.000, -135.000, -120.000, -105.000, -75.000, -45.000, -15.000, 15.000, // 152 - +045.000, 75.000, 105.000, 120.000, 135.000, 150.000, 165.000, 180.000, // 160 - +195.000, 240.000, 225.000, 210.000, 195.000, 180.000, 165.000, 150.000, // 168 - +135.000, 120.000, 90.000, 60.000, 30.000, 0.000, -30.000, -60.000, // 176 - -090.000, -120.000, -135.000, -150.000, -165.000, -180.000, -195.000, -210.000, // 184 - -225.000, -240.000, -225.000, -210.000, -195.000, -180.000, -165.000, -150.000, // 192 - -135.000, -120.000, -90.000, -60.000, -30.000, 0.000, 30.000, 60.000, // 200 - +090.000, 120.000, 135.000, 150.000, 165.000, 180.000, 195.000, 210.000, // 208 - +225.000, 270.000, 255.000, 240.000, 225.000, 210.000, 195.000, 180.000, // 216 - +165.000, 150.000, 135.000, 105.000, 75.000, 45.000, 15.000, -15.000, // 224 - -045.000, -75.000, -105.000, -135.000, -150.000, -165.000, -180.000, -195.000, // 232 - -210.000, -225.000, -240.000, -255.000, -270.000, -255.000, -240.000, -225.000, // 240 - -210.000, -195.000, -180.000, -165.000, -150.000, -135.000, -105.000, -75.000, // 248 - -045.000, -15.000, 15.000, 45.000, 75.000, 105.000, 135.000, 150.000, // 256 - +165.000, 180.000, 195.000, 210.000, 225.000, 240.000, 255.000, 300.000, // 264 - +285.000, 270.000, 255.000, 240.000, 225.000, 210.000, 195.000, 180.000, // 272 - +165.000, 150.000, 120.000, 90.000, 60.000, 30.000, 0.000, -30.000, // 280 - -060.000, -90.000, -120.000, -150.000, -165.000, -180.000, -195.000, -210.000, // 288 - -225.000, -240.000, -255.000, -270.000, -285.000, -300.000, -285.000, -270.000, // 296 - -255.000, -240.000, -225.000, -210.000, -195.000, -180.000, -165.000, -150.000, // 304 - -120.000, -90.000, -60.000, -30.000, 0.000, 30.000, 60.000, 90.000, // 312 - +120.000, 150.000, 165.000, 180.000, 195.000, 210.000, 225.000, 240.000, // 320 - +255.000, 270.000, 285.000, 330.000, 315.000, 300.000, 285.000, 270.000, // 328 - +255.000, 240.000, 225.000, 210.000, 195.000, 180.000, 165.000, 135.000, // 336 - +105.000, 75.000, 45.000, 15.000, -15.000, -45.000, -75.000, -105.000, // 344 - -135.000, -165.000, -180.000, -195.000, -210.000, -225.000, -240.000, -255.000, // 352 - -270.000, -285.000, -300.000, -315.000, -330.000, -315.000, -300.000, -285.000, // 360 - -270.000, -255.000, -240.000, -225.000, -210.000, -195.000, -180.000, -165.000, // 368 - -135.000, -105.000, -75.000, -45.000, -15.000, 15.000, 45.000, 75.000, // 376 - +105.000, 135.000, 165.000, 180.000, 195.000, 210.000, 225.000, 240.000, // 384 - +255.000, 270.000, 285.000, 300.000, 315.000, 360.000, 330.000, 300.000, // 392 - +270.000, 240.000, 210.000, 150.000, 90.000, 30.000, -30.000, -90.000, // 400 - -150.000, -210.000, -240.000, -270.000, -300.000, -330.000, -360.000, -360.000, // 408 - -330.000, -300.000, -270.000, -240.000, -210.000, -150.000, -90.000, -30.000, // 416 - +030.000, 90.000, 150.000, 210.000, 240.000, 270.000, 300.000, 330.000, // 424 - +360.000, 420.000, 390.000, 360.000, 330.000, 300.000, 270.000, 240.000, // 432 - +180.000, 120.000, 60.000, 0.000, -60.000, -120.000, -180.000, -240.000, // 440 - -270.000, -300.000, -330.000, -360.000, -390.000, -420.000, -420.000, -390.000, // 448 - -360.000, -330.000, -300.000, -270.000, -240.000, -180.000, -120.000, -60.000, // 456 - +000.000, 60.000, 120.000, 180.000, 240.000, 270.000, 300.000, 330.000, // 464 - +360.000, 390.000, 420.000, 480.000, 450.000, 420.000, 390.000, 360.000, // 472 - +330.000, 300.000, 270.000, 210.000, 150.000, 90.000, 30.000, -30.000, // 480 - -090.000, -150.000, -210.000, -270.000, -300.000, -330.000, -360.000, -390.000, // 488 - -420.000, -450.000, -480.000, -480.000, -450.000, -420.000, -390.000, -360.000, // 496 - -330.000, -300.000, -270.000, -210.000, -150.000, -90.000, -30.000, 30.000, // 504 - +090.000, 150.000, 210.000, 270.000, 300.000, 330.000, 360.000, 390.000, // 512 - +420.000, 450.000, 480.000, 540.000, 510.000, 480.000, 450.000, 420.000, // 520 - +390.000, 360.000, 330.000, 300.000, 240.000, 180.000, 120.000, 60.000, // 528 - +000.000, -60.000, -120.000, -180.000, -240.000, -300.000, -330.000, -360.000, // 536 - -390.000, -420.000, -450.000, -480.000, -510.000, -540.000, -540.000, -510.000, // 544 - -480.000, -450.000, -420.000, -390.000, -360.000, -330.000, -300.000, -240.000, // 552 - -180.000, -120.000, -60.000, 0.000, 60.000, 120.000, 180.000, 240.000, // 560 - +300.000, 330.000, 360.000, 390.000, 420.000, 450.000, 480.000, 510.000, // 568 - +540.000 // 576 - }; - - const Float_t ytemp[577] = { - +000.000, 0.000, 25.981, 25.981, 0.000, -25.981, -25.981, 0.000, // 0 - +025.981, 51.961, 51.961, 51.961, 25.981, 0.000, -25.981, -51.961, // 8 - -051.961, -51.961, -25.981, 0.000, 25.981, 51.961, 77.942, 77.942, // 16 - +077.942, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 24 - -077.942, -77.942, -77.942, -51.961, -25.981, 0.000, 25.981, 51.961, // 32 - +077.942, 103.923, 103.923, 103.923, 103.923, 103.923, 77.942, 51.961, // 40 - +025.981, 0.000, -25.981, -51.961, -77.942, -103.923, -103.923, -103.923, // 48 - -103.923, -103.923, -77.942, -51.961, -25.981, 0.000, 25.981, 51.961, // 56 - +077.942, 103.923, 129.904, 129.904, 129.904, 129.904, 129.904, 129.904, // 64 - +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 72 - -103.923, -129.904, -129.904, -129.904, -129.904, -129.904, -129.904, -103.923, // 80 - -077.942, -51.961, -25.981, 0.000, 25.981, 51.961, 77.942, 103.923, // 88 - +129.904, 155.885, 155.885, 155.885, 155.885, 155.885, 155.885, 155.885, // 96 - +129.904, 103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, // 104 - -077.942, -103.923, -129.904, -155.885, -155.885, -155.885, -155.885, -155.885, // 112 - -155.885, -155.885, -129.904, -103.923, -77.942, -51.961, -25.981, 0.000, // 120 - +025.981, 51.961, 77.942, 103.923, 129.904, 155.885, 181.865, 181.865, // 128 - +181.865, 181.865, 181.865, 181.865, 181.865, 181.865, 155.885, 129.904, // 136 - +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 144 - -103.923, -129.904, -155.885, -181.865, -181.865, -181.865, -181.865, -181.865, // 152 - -181.865, -181.865, -181.865, -155.885, -129.904, -103.923, -77.942, -51.961, // 160 - -025.981, 0.000, 25.981, 51.961, 77.942, 103.923, 129.904, 155.885, // 168 - +181.865, 207.846, 207.846, 207.846, 207.846, 207.846, 207.846, 207.846, // 176 - +207.846, 207.846, 181.865, 155.885, 129.904, 103.923, 77.942, 51.961, // 184 - +025.981, 0.000, -25.981, -51.961, -77.942, -103.923, -129.904, -155.885, // 192 - -181.865, -207.846, -207.846, -207.846, -207.846, -207.846, -207.846, -207.846, // 200 - -207.846, -207.846, -181.865, -155.885, -129.904, -103.923, -77.942, -51.961, // 208 - -025.981, 0.000, 25.981, 51.961, 77.942, 103.923, 129.904, 155.885, // 216 - +181.865, 207.846, 233.827, 233.827, 233.827, 233.827, 233.827, 233.827, // 224 - +233.827, 233.827, 233.827, 233.827, 207.846, 181.865, 155.885, 129.904, // 232 - +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 240 - -103.923, -129.904, -155.885, -181.865, -207.846, -233.827, -233.827, -233.827, // 248 - -233.827, -233.827, -233.827, -233.827, -233.827, -233.827, -233.827, -207.846, // 256 - -181.865, -155.885, -129.904, -103.923, -77.942, -51.961, -25.981, 0.000, // 264 - +025.981, 51.961, 77.942, 103.923, 129.904, 155.885, 181.865, 207.846, // 272 - +233.827, 259.808, 259.808, 259.808, 259.808, 259.808, 259.808, 259.808, // 280 - +259.808, 259.808, 259.808, 259.808, 233.827, 207.846, 181.865, 155.885, // 288 - +129.904, 103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, // 296 - -077.942, -103.923, -129.904, -155.885, -181.865, -207.846, -233.827, -259.808, // 304 - -259.808, -259.808, -259.808, -259.808, -259.808, -259.808, -259.808, -259.808, // 312 - -259.808, -259.808, -233.827, -207.846, -181.865, -155.885, -129.904, -103.923, // 320 - -077.942, -51.961, -25.981, 0.000, 25.981, 51.961, 77.942, 103.923, // 328 - +129.904, 155.885, 181.865, 207.846, 233.827, 259.808, 285.788, 285.788, // 336 - +285.788, 285.788, 285.788, 285.788, 285.788, 285.788, 285.788, 285.788, // 344 - +285.788, 285.788, 259.808, 233.827, 207.846, 181.865, 155.885, 129.904, // 352 - +103.923, 77.942, 51.961, 25.981, 0.000, -25.981, -51.961, -77.942, // 360 - -103.923, -129.904, -155.885, -181.865, -207.846, -233.827, -259.808, -285.788, // 368 - -285.788, -285.788, -285.788, -285.788, -285.788, -285.788, -285.788, -285.788, // 376 - -285.788, -285.788, -285.788, -259.808, -233.827, -207.846, -181.865, -155.885, // 384 - -129.904, -103.923, -77.942, -51.961, -25.981, 34.641, 86.603, 138.564, // 392 - +190.526, 242.487, 294.449, 329.090, 329.090, 329.090, 329.090, 329.090, // 400 - +329.090, 294.449, 242.487, 190.526, 138.564, 86.603, 34.641, -34.641, // 408 - -086.603, -138.564, -190.526, -242.487, -294.449, -329.090, -329.090, -329.090, // 416 - -329.090, -329.090, -329.090, -294.449, -242.487, -190.526, -138.564, -86.603, // 424 - -034.641, 34.641, 86.603, 138.564, 190.526, 242.487, 294.449, 346.410, // 432 - +381.051, 381.051, 381.051, 381.051, 381.051, 381.051, 381.051, 346.410, // 440 - +294.449, 242.487, 190.526, 138.564, 86.603, 34.641, -34.641, -86.603, // 448 - -138.564, -190.526, -242.487, -294.449, -346.410, -381.051, -381.051, -381.051, // 456 - -381.051, -381.051, -381.051, -381.051, -346.410, -294.449, -242.487, -190.526, // 464 - -138.564, -86.603, -34.641, 34.641, 86.603, 138.564, 190.526, 242.487, // 472 - +294.449, 346.410, 398.372, 433.013, 433.013, 433.013, 433.013, 433.013, // 480 - +433.013, 433.013, 433.013, 398.372, 346.410, 294.449, 242.487, 190.526, // 488 - +138.564, 86.603, 34.641, -34.641, -86.603, -138.564, -190.526, -242.487, // 496 - -294.449, -346.410, -398.372, -433.013, -433.013, -433.013, -433.013, -433.013, // 504 - -433.013, -433.013, -433.013, -398.372, -346.410, -294.449, -242.487, -190.526, // 512 - -138.564, -86.603, -34.641, 34.641, 86.603, 138.564, 190.526, 242.487, // 520 - +294.449, 346.410, 398.372, 450.333, 484.974, 484.974, 484.974, 484.974, // 528 - +484.974, 484.974, 484.974, 484.974, 484.974, 450.333, 398.372, 346.410, // 536 - +294.449, 242.487, 190.526, 138.564, 86.603, 34.641, -34.641, -86.603, // 544 - -138.564, -190.526, -242.487, -294.449, -346.410, -398.372, -450.333, -484.974, // 552 - -484.974, -484.974, -484.974, -484.974, -484.974, -484.974, -484.974, -484.974, // 560 - -450.333, -398.372, -346.410, -294.449, -242.487, -190.526, -138.564, -86.603, // 568 - -034.641 // 576 - }; - - const Float_t rtemp[577] = { - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 0 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 8 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 16 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 24 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 32 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 40 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 48 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 56 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 64 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 72 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 80 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 88 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 96 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 104 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 112 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 120 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 128 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 136 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 144 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 152 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 160 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 168 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 176 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 184 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 192 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 200 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 208 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 216 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 224 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 232 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 240 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 248 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 256 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 264 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 272 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 280 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 288 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 296 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 304 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 312 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 320 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 328 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 336 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 344 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 352 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 360 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 368 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 376 - 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, 30.00, // 384 - 30.00, 30.00, 30.00, 30.00, 30.00, 60.00, 60.00, 60.00, // 392 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 400 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 408 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 416 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 424 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 432 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 440 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 448 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 456 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 464 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 472 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 480 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 488 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 496 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 504 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 512 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 520 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 528 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 536 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 544 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 552 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 560 - 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, 60.00, // 568 - 60.00 }; // 576 - - // - // fill the pixels list with this data - // - - for (UInt_t i=0; i -! Author(s): Harald Kornmayer 1/2001 -! -! Copyright: MAGIC Software Development, 2000-2002 -! -! -\* ======================================================================== */ - -////////////////////////////////////////////////////////////////////////////// -// -// MGeomPix -// -// This container stores the geometry (position) information of -// a single pixel together with the information about next neighbors. -// -//////////////////////////////////////////////////////////////////////////// - -#include "MGeomPix.h" - -#include - -#include "MLog.h" -#include "MLogManip.h" - -ClassImp(MGeomPix); - -// -------------------------------------------------------------------------- -// -// Initialiyes one pixel -// -MGeomPix::MGeomPix(Float_t x, Float_t y, Float_t r) : fX(x), fY(y), fR(r) -{ - // default constructor -} - -// -------------------------------------------------------------------------- -// -// Return the area of the pixel. A hegagonal shape is assumed. -// -Float_t MGeomPix::GetA() const -{ - return fR*fR*tan(60/kRad2Deg); -} - -// -------------------------------------------------------------------------- -// -// Initialiyes Next Neighbors. -// -// WARNING: This function is public, but it is not ment for user access. -// It should only be used from geometry classes (like MGeomCam) -// -void MGeomPix::SetNeighbors(Short_t i0, Short_t i1, Short_t i2, - Short_t i3, Short_t i4, Short_t i5) -{ - fNeighbors[0] = i0; - fNeighbors[1] = i1; - fNeighbors[2] = i2; - fNeighbors[3] = i3; - fNeighbors[4] = i4; - fNeighbors[5] = i5; - - int i; - for (i=0; i<6; i++) - if (fNeighbors[i]<0) - break; - - fNumNeighbors = i; -} - -// -------------------------------------------------------------------------- -// -// Print the geometry information of one pixel. -// -void MGeomPix::Print(Option_t *opt) const -{ - // information about a pixel - *fLog << all << "MPixGeom: x= " << fX - << " y= " << fY - << " r= " << fR - << endl ; -} - Index: trunk/MagicSoft/Mars/mgui/MGeomPix.h =================================================================== --- trunk/MagicSoft/Mars/mgui/MGeomPix.h (revision 1435) +++ (revision ) @@ -1,46 +1,0 @@ -#ifndef MARS_MGeomPix -#define MARS_MGeomPix - -#ifndef MARS_MParContainer -#include "MParContainer.h" -#endif - -class MGeomPix : public MParContainer -{ -private: - Float_t fX; // [mm] the x coordinate of the center - Float_t fY; // [mm] the y coordinate of the center - Float_t fR; // [mm] the r coordinate of the pixel (dist between two parallel sides) - - Byte_t fNumNeighbors; // number of valid neighbors - Short_t fNeighbors[6]; // the IDs of the pixel next to it - // we are assuming an hexagonal geometry - -public: - MGeomPix(Float_t x=0, Float_t y=0, Float_t r=0); - - void Print(Option_t *opt=NULL) const; - - void Set(Float_t x, Float_t y, Float_t r) { fX=x; fY=y; fR=r; } - - void SetNeighbors(Short_t i0=-1, Short_t i1=-1, Short_t i2=-1, - Short_t i3=-1, Short_t i4=-1, Short_t i5=-1); - - void SetX(Float_t x) { fX = x; } - void SetY(Float_t y) { fY = y; } - void SetR(Float_t r) { fR = r; } - - Float_t GetX() const { return fX; } - Float_t GetY() const { return fY; } - Float_t GetR() const { return fR; } - - Float_t GetA() const; - - Byte_t GetNumNeighbors() const { return fNumNeighbors; } - Short_t GetNeighbor(Byte_t i) const { return fNeighbors[i]; } - - ClassDef(MGeomPix, 1) // Geometry class describing the geometry of one pixel -}; - -#endif - Index: trunk/MagicSoft/Mars/mgui/Makefile =================================================================== --- trunk/MagicSoft/Mars/mgui/Makefile (revision 1435) +++ trunk/MagicSoft/Mars/mgui/Makefile (revision 1436) @@ -22,5 +22,5 @@ # connect the include files defined in the config.mk file # -INCLUDES = -I. -I../mbase -I../manalysis -I../mdatacheck -I../meventdisp +INCLUDES = -I. -I../mbase -I../mgeom -I../manalysis #------------------------------------------------------------------------------ @@ -29,10 +29,5 @@ SRCFILES = MHexagon.cc \ - MGeomPix.cc \ - MGeomCam.cc \ - MCamDisplay.cc \ - MineSweeper.cc \ - MGeomCamCT1.cc \ - MGeomCamMagic.cc + MCamDisplay.cc SRCS = $(SRCFILES)