| 1 | /* ======================================================================== *\
|
|---|
| 2 | !
|
|---|
| 3 | ! *
|
|---|
| 4 | ! * This file is part of MARS, the MAGIC Analysis and Reconstruction
|
|---|
| 5 | ! * Software. It is distributed to you in the hope that it can be a useful
|
|---|
| 6 | ! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
|
|---|
| 7 | ! * It is distributed WITHOUT ANY WARRANTY.
|
|---|
| 8 | ! *
|
|---|
| 9 | ! * Permission to use, copy, modify and distribute this software and its
|
|---|
| 10 | ! * documentation for any purpose is hereby granted without fee,
|
|---|
| 11 | ! * provided that the above copyright notice appear in all copies and
|
|---|
| 12 | ! * that both that copyright notice and this permission notice appear
|
|---|
| 13 | ! * in supporting documentation. It is provided "as is" without express
|
|---|
| 14 | ! * or implied warranty.
|
|---|
| 15 | ! *
|
|---|
| 16 | !
|
|---|
| 17 | !
|
|---|
| 18 | ! Author(s): Thomas Bretz, 6/2003 <mailto:tbretz@astro.uni-wuerzburg.de>
|
|---|
| 19 | !
|
|---|
| 20 | ! Copyright: MAGIC Software Development, 2000-2003
|
|---|
| 21 | !
|
|---|
| 22 | !
|
|---|
| 23 | \* ======================================================================== */
|
|---|
| 24 |
|
|---|
| 25 | //////////////////////////////////////////////////////////////////////////////
|
|---|
| 26 | //
|
|---|
| 27 | // MCamEvent
|
|---|
| 28 | //
|
|---|
| 29 | // A base class describing an event in the camera.
|
|---|
| 30 | //
|
|---|
| 31 | // If GetPixelContent returns kFALSE, it must not touch 'val'.
|
|---|
| 32 | //
|
|---|
| 33 | // You can derive a class in addition to TObject from MCamEvent, too.
|
|---|
| 34 | //
|
|---|
| 35 | // MCamEvent MUST be after TObject:
|
|---|
| 36 | // ALLOWED: class MyClass : public TObject, public MCamEvent
|
|---|
| 37 | // FORBIDDEN: class MyClass : public MCamEvent, public TObject
|
|---|
| 38 | //
|
|---|
| 39 | //////////////////////////////////////////////////////////////////////////////
|
|---|
| 40 | #include "MCamEvent.h"
|
|---|
| 41 |
|
|---|
| 42 | #include "MGeomCam.h"
|
|---|
| 43 |
|
|---|
| 44 | ClassImp(MCamEvent);
|
|---|
| 45 |
|
|---|
| 46 | // --------------------------------------------------------------------------
|
|---|
| 47 | //
|
|---|
| 48 | // You can overwrite this function if you want the container to be
|
|---|
| 49 | // initialized by MGeomApply with the geometry. If it is not overloaded
|
|---|
| 50 | // it calls InitSize with the corresponding pixel number. If this information
|
|---|
| 51 | // is enough for you it is enough to overload InitSize.
|
|---|
| 52 | //
|
|---|
| 53 | void MCamEvent::Init(const MGeomCam &geom)
|
|---|
| 54 | {
|
|---|
| 55 | InitSize(geom.GetNumPixels());
|
|---|
| 56 | }
|
|---|
| 57 |
|
|---|
| 58 | // --------------------------------------------------------------------------
|
|---|
| 59 | //
|
|---|
| 60 | // Return the mean of all corresponding GetPixelContent
|
|---|
| 61 | //
|
|---|
| 62 | Double_t MCamEvent::GetCameraMean(const MGeomCam &cam, Int_t type) const
|
|---|
| 63 | {
|
|---|
| 64 | Int_t num = 0;
|
|---|
| 65 | Double_t mean = 0;
|
|---|
| 66 | for (unsigned int i=0; i<cam.GetNumPixels(); i++)
|
|---|
| 67 | {
|
|---|
| 68 | Double_t val;
|
|---|
| 69 | if (!GetPixelContent(val, i, cam, type))
|
|---|
| 70 | continue;
|
|---|
| 71 |
|
|---|
| 72 | mean += val;
|
|---|
| 73 | num ++;
|
|---|
| 74 | }
|
|---|
| 75 | return num == 0 ? 0 : mean/num;
|
|---|
| 76 | }
|
|---|
