source: trunk/MagicSoft/Mars/macros/collarea.C@ 1885

Last change on this file since 1885 was 1660, checked in by tbretz, 22 years ago
*** empty log message ***
File size: 2.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 12/2000 (tbretz@uni-sw.gwdg.de)
19!
20! Copyright: MAGIC Software Development, 2000-2001
21!
22!
23\* ======================================================================== */
24
25
26void collarea(TString filename="camera.root", TString outname="")
27{
28 //
29 // first we have to create our empty lists
30 //
31 MParList parlist;
32 MTaskList tasklist;
33
34 parlist.AddToList(&tasklist);
35
36 //
37 // Setup out tasks:
38 // - First we have to read the events
39 // - Then we can fill the efficiency histograms
40 //
41 MReadMarsFile reader("Events", filename);
42 MMcCollectionAreaCalc effi;
43
44 tasklist.AddToList(&reader);
45 tasklist.AddToList(&effi);
46
47 //
48 // set up the loop for the processing
49 //
50 MEvtLoop magic;
51 magic.SetParList(&parlist);
52
53 //
54 // Start to loop over all events
55 //
56 MProgressBar bar;
57 magic.SetProgressBar(&bar);
58 if (!magic.Eventloop())
59 return;
60
61 tasklist.PrintStatistics();
62
63 //
64 // Now the histogram we wanted to get out of the data is
65 // filled and can be displayed
66 //
67 parlist.FindObject("MHMcCollectionArea")->DrawClone();
68
69
70 //
71 // Write histogram to a file in case an output
72 // filename has been supplied
73 //
74 if (outname.IsNull())
75 return;
76
77 TFile f(outname, "recreate");
78 if (!f)
79 return;
80
81 parlist.FindObject("MHMcCollectionArea")->Write();
82}
Note: See TracBrowser for help on using the repository browser.