source: trunk/MagicSoft/Mars/mgui/MCamEvent.cc@ 9450

Last change on this file since 9450 was 8907, checked in by tbretz, 16 years ago
*** empty log message ***
File size: 3.8 KB
Line 
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 <TMath.h>
43
44#include "MGeomCam.h"
45
46ClassImp(MCamEvent);
47
48// --------------------------------------------------------------------------
49//
50// You can overwrite this function if you want the container to be
51// initialized by MGeomApply with the geometry. If it is not overloaded
52// it calls InitSize with the corresponding pixel number. If this information
53// is enough for you it is enough to overload InitSize.
54//
55void MCamEvent::Init(const MGeomCam &geom)
56{
57 InitSize(geom.GetNumPixels());
58}
59
60// --------------------------------------------------------------------------
61//
62// Return the mean of all corresponding GetPixelContent
63//
64Double_t MCamEvent::GetCameraMean(const MGeomCam &cam, Int_t type) const
65{
66 return GetCameraStat(cam, type)[0];
67 /*
68 Int_t num = 0;
69 Double_t mean = 0;
70 for (unsigned int i=0; i<cam.GetNumPixels(); i++)
71 {
72 Double_t val;
73 if (!GetPixelContent(val, i, cam, type))
74 continue;
75
76 mean += val;
77 num ++;
78 }
79 return num == 0 ? 0 : mean/num;
80 */
81}
82
83// --------------------------------------------------------------------------
84//
85// Return the rms of all corresponding GetPixelContent
86//
87Double_t MCamEvent::GetCameraRMS(const MGeomCam &cam, Int_t type) const
88{
89 return GetCameraStat(cam, type)[1];
90 /*
91 Int_t num = 0;
92 Double_t sum = 0;
93 Double_t sq = 0;
94
95 for (unsigned int i=0; i<cam.GetNumPixels(); i++)
96 {
97 Double_t val;
98 if (!GetPixelContent(val, i, cam, type))
99 continue;
100
101 sum += val;
102 sq += val;
103 num ++;
104 }
105
106 if (num==0)
107 return 0;
108
109 sum /= num;
110 sq /= num;
111
112 return TMath::Sqrt(sq-sum*sum);
113 */
114}
115
116// --------------------------------------------------------------------------
117//
118// Return the rms of all corresponding GetPixelContent
119//
120TArrayD MCamEvent::GetCameraStat(const MGeomCam &cam, Int_t type) const
121{
122 Int_t num = 0;
123 Double_t sum = 0;
124 Double_t sq = 0;
125
126 for (unsigned int i=0; i<cam.GetNumPixels(); i++)
127 {
128 Double_t val;
129 if (!GetPixelContent(val, i, cam, type))
130 continue;
131
132 sum += val;
133 sq += val;
134 num ++;
135 }
136
137 if (num==0)
138 return 0;
139
140 sum /= num;
141 sq /= num;
142
143 TArrayD stat(2);
144 stat[0] = sum;
145 stat[1] = TMath::Sqrt(sq-sum*sum);
146 return stat;
147}
Note: See TracBrowser for help on using the repository browser.