source: trunk/MagicSoft/Mars/macros/CT1collarea.C@ 1837

Last change on this file since 1837 was 1837, checked in by moralejo, 22 years ago
*** empty log message ***
File size: 2.5 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 CT1collarea(TString filename="MC_ON2.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 reader.DisableAutoScheme();
43 // reader.EnableBranch("fHadronness");
44
45 MMcCT1CollectionAreaCalc effi;
46 MHMcCT1CollectionArea* collarea = new MHMcCT1CollectionArea("MHMcCT1CollectionArea","",30,2.,5.);
47 parlist.AddToList(collarea);
48
49 tasklist.AddToList(&reader);
50
51 MF filterhadrons("MHadronness.fHadronness<0.25");
52 tasklist.AddToList(&filterhadrons);
53 effi.SetFilter(&filterhadrons);
54
55 tasklist.AddToList(&effi);
56
57 //
58 // set up the loop for the processing
59 //
60 MEvtLoop magic;
61 magic.SetParList(&parlist);
62
63 //
64 // Start to loop over all events
65 //
66 MProgressBar bar;
67 magic.SetProgressBar(&bar);
68 if (!magic.Eventloop())
69 return;
70
71 tasklist.PrintStatistics();
72
73 //
74 // Now the histogram we wanted to get out of the data is
75 // filled and can be displayed
76 //
77 parlist.FindObject("MHMcCT1CollectionArea")->DrawClone();
78
79 //
80 // Write histogram to a file in case an output
81 // filename has been supplied
82 //
83 if (outname.IsNull())
84 return;
85
86 TFile f("area.root","recreate");
87 collarea->GetHist()->Write();
88 collarea->GetHAll()->Write();
89 collarea->GetHSel()->Write();
90 f.Close();
91}
92
Note: See TracBrowser for help on using the repository browser.