source: branches/MarsISDCBranchBasedOn17887/mimage/MHHillas.cc@ 19921

Last change on this file since 19921 was 12866, checked in by tbretz, 13 years ago
Removed the camera display, because the display is uusally in degrees and not in millimeters.
File size: 9.1 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 2001 <mailto:tbretz@astro.uni-wuerzburg.de>
19! Author(s): Wolfgang Wittek 2002 <mailto:wittek@mppmu.mpg.de>
20!
21! Copyright: MAGIC Software Development, 2000-2009
22!
23!
24\* ======================================================================== */
25
26/////////////////////////////////////////////////////////////////////////////
27//
28// MHHillas
29//
30// This class contains histograms for the source independent image parameters
31//
32// ClassVersion 2:
33// ---------------
34// - fMm2Deg
35// - fUseMmScale
36//
37/////////////////////////////////////////////////////////////////////////////
38#include "MHHillas.h"
39
40#include <math.h>
41
42#include <TH1.h>
43#include <TH2.h>
44#include <TPad.h>
45#include <TStyle.h>
46#include <TCanvas.h>
47
48#include "MLog.h"
49#include "MLogManip.h"
50
51#include "MParList.h"
52
53#include "MHillas.h"
54#include "MGeomCam.h"
55#include "MBinning.h"
56
57#include "MHCamera.h"
58
59ClassImp(MHHillas);
60
61using namespace std;
62
63// --------------------------------------------------------------------------
64//
65// Setup four histograms for Width, Length
66//
67MHHillas::MHHillas(const char *name, const char *title)
68 : fGeomCam(0)
69{
70 //
71 // set the name and title of this object
72 //
73 fName = name ? name : "MHHillas";
74 fTitle = title ? title : "Source independent image parameters";
75
76 fLength = new TH1F("Length", "Length of Ellipse", 100, 0, 1.0);
77 fWidth = new TH1F("Width", "Width of Ellipse", 100, 0, 1.0);
78 fDistC = new TH1F("DistC", "Distance from center of camera", 100, 0, 1.5);
79 fDelta = new TH1F("Delta", "Angle (Main axis - x-axis)", 101, -90, 90);
80
81 fLength->SetLineColor(kBlue);
82
83 fLength->SetDirectory(NULL);
84 fWidth->SetDirectory(NULL);
85 fDistC->SetDirectory(NULL);
86 fDelta->SetDirectory(NULL);
87
88 fLength->SetXTitle("Length [\\circ]");
89 fWidth->SetXTitle("Width [\\circ]");
90 fDistC->SetXTitle("Distance [\\circ]");
91 fDelta->SetXTitle("Delta [\\circ]");
92
93 fLength->SetYTitle("Counts");
94 fWidth->SetYTitle("Counts");
95 fDistC->SetYTitle("Counts");
96 fDelta->SetYTitle("Counts");
97
98 fDelta->SetMinimum(0);
99
100 MBinning bins;
101 bins.SetEdgesLog(50, 1, 1e7);
102
103 fSize = new TH1F;
104 fSize->SetName("Size");
105 fSize->SetTitle("Number of Photons");
106 fSize->SetDirectory(NULL);
107 fSize->SetXTitle("Size");
108 fSize->SetYTitle("Counts");
109 fSize->GetXaxis()->SetTitleOffset(1.2);
110 fSize->GetXaxis()->SetLabelOffset(-0.015);
111 fSize->SetFillStyle(4000);
112 fSize->UseCurrentStyle();
113
114 bins.Apply(*fSize);
115
116 fCenter = new TH2F("Center", "Center of gravity", 51, -1.5, 1.5, 51, -1.5, 1.5);
117 fCenter->SetDirectory(NULL);
118 fCenter->SetXTitle("x [mm]");
119 fCenter->SetYTitle("y [mm]");
120 fCenter->SetZTitle("Counts");
121}
122
123// --------------------------------------------------------------------------
124//
125// Delete the histograms
126//
127MHHillas::~MHHillas()
128{
129 delete fLength;
130 delete fWidth;
131
132 delete fDistC;
133 delete fDelta;
134
135 delete fSize;
136 delete fCenter;
137}
138
139// --------------------------------------------------------------------------
140//
141// Setup the Binning for the histograms automatically if the correct
142// instances of MBinning (with the names 'BinningWidth' and 'BinningLength')
143// are found in the parameter list
144// Use this function if you want to set the conversion factor which
145// is used to convert the mm-scale in the camera plain into the deg-scale
146// used for histogram presentations. The conversion factor is part of
147// the camera geometry. Please create a corresponding MGeomCam container.
148//
149Bool_t MHHillas::SetupFill(const MParList *plist)
150{
151 fGeomCam = (MGeomCam*)plist->FindObject("MGeomCam");
152 if (!fGeomCam)
153 {
154 *fLog << err << "MGeomCam not found... abort." << endl;
155 return kFALSE;
156 }
157
158
159 ApplyBinning(*plist, "Width", *fWidth);
160 ApplyBinning(*plist, "Length", *fLength);
161 ApplyBinning(*plist, "Dist", *fDistC);
162 ApplyBinning(*plist, "Delta", *fDelta);
163 ApplyBinning(*plist, "Size", *fSize);
164
165 const MBinning *bins = (MBinning*)plist->FindObject("BinningCamera");
166 if (!bins)
167 {
168 const Float_t r = fGeomCam->GetMaxRadius()*fGeomCam->GetConvMm2Deg();
169
170 const MBinning b(61, -r, r);
171 SetBinning(*fCenter, b, b);
172 }
173 else
174 SetBinning(*fCenter, *bins, *bins);
175
176
177 return kTRUE;
178}
179
180// --------------------------------------------------------------------------
181//
182// Fill the histograms with data from a MHillas-Container.
183// Be careful: Only call this with an object of type MHillas
184//
185Int_t MHHillas::Fill(const MParContainer *par, const Stat_t w)
186{
187 if (!par)
188 {
189 *fLog << err << "MHHillas::Fill: Pointer (!=NULL) expected." << endl;
190 return kERROR;
191 }
192
193 const MHillas &h = *(MHillas*)par;
194
195 const Double_t scale = fGeomCam->GetConvMm2Deg();
196
197 fLength->Fill(scale*h.GetLength(), w);
198 fWidth ->Fill(scale*h.GetWidth(), w);
199 fDistC ->Fill(scale*h.GetDist0(), w);
200 fCenter->Fill(scale*h.GetMeanX(), scale*h.GetMeanY(), w);
201 fDelta ->Fill(kRad2Deg*h.GetDelta(), w);
202 fSize ->Fill(h.GetSize(), w);
203
204 return kTRUE;
205}
206
207// --------------------------------------------------------------------------
208//
209// Creates a new canvas and draws the four histograms into it.
210// Be careful: The histograms belongs to this object and won't get deleted
211// together with the canvas.
212//
213void MHHillas::Draw(Option_t *o)
214{
215 TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this);
216 pad->SetBorderMode(0);
217
218 AppendPad("");
219
220 TString opt(o);
221 opt.ToLower();
222
223 // FIXME: If same-option given make two independant y-axis!
224 const Bool_t same = opt.Contains("same");
225
226 if (!same)
227 pad->Divide(2,3);
228 else
229 {
230 fLength->SetName("LengthSame");
231 fWidth->SetName("WidthSame");
232 fDistC->SetName("DistCSame");
233 fDelta->SetName("DeltaSame");
234 fSize->SetName("SizeSame");
235 fCenter->SetName("CenterSame");
236
237 fLength->SetDirectory(0);
238 fWidth->SetDirectory(0);
239 fDistC->SetDirectory(0);
240 fDelta->SetDirectory(0);
241 fSize->SetDirectory(0);
242 fCenter->SetDirectory(0);
243
244 fDistC->SetLineColor(kBlue);
245 fSize->SetLineColor(kBlue);
246 fDelta->SetLineColor(kBlue);
247 fWidth->SetLineColor(kMagenta);
248 fLength->SetLineColor(kCyan);
249 }
250
251 pad->cd(1);
252 gPad->SetBorderMode(0);
253 gPad->SetGridx();
254 gPad->SetGridy();
255 RemoveFromPad("WidthSame");
256 RemoveFromPad("LengthSame");
257 MH::DrawSame(*fWidth, *fLength, "Width'n'Length", same);
258
259 pad->cd(2);
260 gPad->SetBorderMode(0);
261 gPad->SetGridx();
262 gPad->SetGridy();
263 RemoveFromPad("DistCSame");
264 fDistC->Draw(same?"same":"");
265
266 pad->cd(3);
267 gPad->SetGridx();
268 gPad->SetGridy();
269 gPad->SetBorderMode(0);
270 gPad->SetLogx();
271 gPad->SetLogy();
272 RemoveFromPad("SizeSame");
273 fSize->Draw(same?"same":"");
274
275 pad->cd(4);
276 gPad->SetBorderMode(0);
277 gPad->SetGridx();
278 gPad->SetGridy();
279 gPad->SetPad(0.51, 0.01, 0.99, 0.65);
280 if (same)
281 {
282 TH2 *h=dynamic_cast<TH2*>(gPad->FindObject("Center"));
283 if (h)
284 {
285 // This causes crashes in THistPainter::PaintTable
286 // if the z-axis is not kept. No idea why...
287 h->SetDrawOption("z");
288 h->SetMarkerColor(kBlack);
289 }
290
291 RemoveFromPad("CenterSame");
292 fCenter->SetMarkerColor(kBlue);
293 fCenter->Draw("same");
294 }
295 else
296 fCenter->Draw("colz");
297
298 if (fGeomCam)
299 {
300// MHCamera *cam = new MHCamera(*fGeomCam);
301// cam->Draw("same");
302// cam->SetBit(kCanDelete);
303 }
304
305 pad->cd(5);
306 gPad->SetBorderMode(0);
307 gPad->SetGridx();
308 gPad->SetGridy();
309 RemoveFromPad("DeltaSame");
310 fDelta->Draw(same?"same":"");
311
312 pad->cd(6);
313 if (gPad && !same)
314 delete gPad;
315}
316
317TH1 *MHHillas::GetHistByName(const TString name) const
318{
319 if (name.Contains("Width", TString::kIgnoreCase))
320 return fWidth;
321 if (name.Contains("Length", TString::kIgnoreCase))
322 return fLength;
323 if (name.Contains("Size", TString::kIgnoreCase))
324 return fSize;
325 if (name.Contains("Delta", TString::kIgnoreCase))
326 return fDelta;
327 if (name.Contains("DistC", TString::kIgnoreCase))
328 return fDistC;
329 if (name.Contains("Center", TString::kIgnoreCase))
330 return fCenter;
331
332 return NULL;
333}
334
335void MHHillas::Paint(Option_t *opt)
336{
337 MH::SetPalette("pretty");
338 MH::Paint();
339}
Note: See TracBrowser for help on using the repository browser.