source: trunk/MagicSoft/Mars/mtemp/mucm/Flux6.C@ 6724

Last change on this file since 6724 was 5456, checked in by marcos, 20 years ago
*** empty log message ***
File size: 3.2 KB
Line 
1/*############################################################################
2 #
3 # Input:
4 # - A 3D hist of Alpha vs Energy vs Theta
5 # - A 2D hist of DiffTime vs Time
6 # - A 2D with Effective Collection Area
7 #
8 # Then calculate:
9 # - Effecive on Time
10 # - Nex for the Alpha plot foreach bin (energy,theta)
11 #
12 #
13 # Fixed a bug: the loop for extracing the signal started in bin=0 instead in
14 # bin 1.
15 #
16 # Fixed a bug: the above change introduced a extremely important bug, since
17 # the loop started from 1 but we have in the arrays:
18 # signal[i] = fNex;
19 # errorsignal[i] = fdNex;
20 # so the signal[0] was actually the signal for the last bin. For that we
21 # had a very low point in the first bin and all the spectrum was displaced.
22 # Now it is fixed with:
23 # signal[i-1] = fNex;
24 # errorsignal[i-1] = fdNex;
25 #
26 # authors: Marcos Lopez
27 # Rober Wagner
28 # email: marcos@gae.ucm.es
29 # date: 31-11-2004
30 # last revision: 31-11-2004
31 ###########################################################################*/
32
33void Flux6()
34{
35
36 //TString datafile = "flux/Alpha10_H3.New.Resized.root";
37 //TString datafile = "flux/Alpha10_HadCut_Variable.Resized.root";
38 TString datafile = "flux/Alpha10_HadCut_Variable.SizeAbove250.Resized.root";
39 //TString datafile = "flux/Alpha_E_50_2000_10bins_HadCut_Variable.SizeAbove250.Resized.root";
40
41
42 //TString areafile = "flux/area_HadCut_03_and_AlphaCut_20.Resized.root";
43 //TString areafile = "../mcdata/collarea.root";
44 //TString areafile = "area_HadCut03_and_AlphaCut20.root";
45 //TString areafile = "area_HacCut_Variable_and_AlphaCut_20.Resized.root";
46 //TString areafile = "area_HacCut_Variable.Resized.root";
47 TString areafile = "flux/area_HadCut_Variable.SizeAbove250.Resized.root";
48
49 // ------------------------------------------------------------------------
50 //
51 // Load histograms from file
52 //
53 TFile* file = new TFile(datafile);
54
55 //
56 // Read Hist of Alpha Vs. Energy and Theta
57 //
58 MHAlphaEnergyTheta hAlpha;
59 hAlpha.Read("MHAlphaEnergyTheta");
60 hAlpha.DrawClone();
61
62 //
63 // Read Hist of EffectiveTime vs. Theta and Time
64 //
65 MHEffectiveOnTime hEffTime;
66 hEffTime.Read("MHEffectiveOnTime");
67 hEffTime.DrawClone();
68
69 //
70 // Read CollectionArea
71 //
72 TFile* file3 = new TFile(areafile);
73 MHMcCollectionArea area;
74 area.Read("MHMcCollectionArea");
75 area.DrawClone();
76
77
78 // -----------------------------------------------------------------------
79 //
80 // Calculate # Excess events vs. Energy and Theta
81 //
82 MHExcessEnergyTheta *hex = new MHExcessEnergyTheta;
83 hex->Calc(&hAlpha);
84 hex->Draw();
85
86 //
87 // Calculate diferential Flux vs. Enregy and Theta
88 //
89 MHFlux* hFlux = new MHFlux;
90 hFlux->Calc(hex, &area, &hEffTime);
91 hFlux->Draw();
92
93
94
95
96 //
97 // Write flux into root file
98 //
99 TFile outfile("flux/Flux.root","RECREATE");
100 hFlux->Write();
101 outfile.Close();
102
103}
104
Note: See TracBrowser for help on using the repository browser.