]>
Commit | Line | Data |
---|---|---|
bf000c32 | 1 | //____________________________________________________________________ |
088f8e79 | 2 | // |
3 | // $Id$ | |
4 | // | |
5 | // Script that contains a class to draw eloss from hits, versus ADC | |
6 | // counts from digits, using the AliFMDInputHits class in the util library. | |
7 | // | |
8 | // It draws the energy loss versus the p/(mq^2). It can be overlayed | |
9 | // with the Bethe-Bloc curve to show how the simulation behaves | |
10 | // relative to the expected. | |
11 | // | |
12 | // Use the script `Compile.C' to compile this class using ACLic. | |
13 | // | |
14 | #include <TH2D.h> | |
15 | #include <AliFMDHit.h> | |
16 | #include <AliFMDDigit.h> | |
17 | #include <AliFMDInput.h> | |
18 | #include <AliFMDEdepMap.h> | |
19 | #include <iostream> | |
20 | #include <TStyle.h> | |
21 | #include <TArrayF.h> | |
8ec606c2 | 22 | #include <AliLog.h> |
088f8e79 | 23 | |
9b48326f | 24 | /** @class DrawHitsDigits |
25 | @brief Draw hit energy loss versus digit ADC | |
26 | @code | |
27 | Root> .L Compile.C | |
28 | Root> Compile("DrawHitsDigits.C") | |
29 | Root> DrawHitsDigits c | |
30 | Root> c.Run(); | |
31 | @endcode | |
32 | @ingroup FMD_script | |
33 | */ | |
8ec606c2 | 34 | class DrawHitsDigits : public AliFMDInput |
088f8e79 | 35 | { |
36 | private: | |
37 | TH2D* fElossVsAdc; // Histogram | |
38 | AliFMDEdepMap fMap; | |
39 | public: | |
bf000c32 | 40 | //__________________________________________________________________ |
088f8e79 | 41 | TArrayF MakeLogScale(Int_t n, Double_t min, Double_t max) |
42 | { | |
43 | TArrayF bins(n+1); | |
44 | Float_t dp = n / TMath::Log10(max / min); | |
45 | Float_t pmin = TMath::Log10(min); | |
46 | bins[0] = min; | |
47 | for (Int_t i = 1; i < n+1; i++) { | |
48 | Float_t p = pmin + i / dp; | |
49 | bins[i] = TMath::Power(10, p); | |
50 | } | |
51 | return bins; | |
52 | } | |
bf000c32 | 53 | //__________________________________________________________________ |
088f8e79 | 54 | DrawHitsDigits(Int_t n=900, Double_t emin=1e-3, Double_t emax=10, |
55 | Int_t m=1100, Double_t amin=-0.5, Double_t amax=1099.5) | |
56 | { | |
57 | AddLoad(kDigits); | |
8ec606c2 | 58 | AddLoad(kHits); |
088f8e79 | 59 | TArrayF eloss(MakeLogScale(n, emin, emax)); |
60 | TArrayF adcs(m+1); | |
61 | adcs[0] = amin; | |
62 | for (Int_t i = 1; i < m+1; i++) adcs[i] = adcs[i-1] + (amax-amin)/m; | |
63 | fElossVsAdc = new TH2D("bad", "#Delta E vs. ADC", | |
64 | eloss.fN-1, eloss.fArray, adcs.fN-1, adcs.fArray); | |
65 | fElossVsAdc->SetXTitle("#Delta E/#Delta x [MeV/cm]"); | |
66 | fElossVsAdc->SetYTitle("ADC value"); | |
67 | } | |
bf000c32 | 68 | //__________________________________________________________________ |
9b48326f | 69 | /** Begining of event |
70 | @param ev Event number | |
71 | @return @c false on error */ | |
088f8e79 | 72 | Bool_t Begin(Int_t ev) |
73 | { | |
74 | fMap.Reset(); | |
8ec606c2 | 75 | return AliFMDInput::Begin(ev); |
088f8e79 | 76 | } |
bf000c32 | 77 | //__________________________________________________________________ |
088f8e79 | 78 | Bool_t ProcessHit(AliFMDHit* hit, TParticle*) |
79 | { | |
80 | // Cache the energy loss | |
81 | if (!hit) return kFALSE; | |
82 | UShort_t det = hit->Detector(); | |
83 | Char_t rng = hit->Ring(); | |
84 | UShort_t sec = hit->Sector(); | |
85 | UShort_t str = hit->Strip(); | |
8f6ee336 | 86 | if (str > 511) { |
87 | AliWarning(Form("Bad strip number %d in hit", str)); | |
88 | return kTRUE; | |
89 | } | |
088f8e79 | 90 | fMap(det, rng, sec, str).fEdep += hit->Edep(); |
91 | fMap(det, rng, sec, str).fN++; | |
92 | return kTRUE; | |
93 | } | |
bf000c32 | 94 | //__________________________________________________________________ |
95 | Bool_t ProcessDigit(AliFMDDigit* digit) | |
088f8e79 | 96 | { |
8ec606c2 | 97 | if (!digit) return kFALSE; |
bf000c32 | 98 | UShort_t det = digit->Detector(); |
99 | Char_t rng = digit->Ring(); | |
100 | UShort_t sec = digit->Sector(); | |
101 | UShort_t str = digit->Strip(); | |
102 | if (str > 511) { | |
103 | AliWarning(Form("Bad strip number %d in digit", str)); | |
104 | continue; | |
088f8e79 | 105 | } |
bf000c32 | 106 | fElossVsAdc->Fill(fMap(det, rng, sec, str).fEdep, digit->Counts()); |
088f8e79 | 107 | return kTRUE; |
108 | } | |
bf000c32 | 109 | //__________________________________________________________________ |
088f8e79 | 110 | Bool_t Finish() |
111 | { | |
112 | gStyle->SetPalette(1); | |
113 | gStyle->SetOptTitle(0); | |
114 | gStyle->SetCanvasColor(0); | |
115 | gStyle->SetCanvasBorderSize(0); | |
116 | gStyle->SetPadColor(0); | |
117 | gStyle->SetPadBorderSize(0); | |
118 | fElossVsAdc->SetStats(kFALSE); | |
119 | fElossVsAdc->Draw("COLZ"); | |
120 | return kTRUE; | |
121 | } | |
8f6ee336 | 122 | |
123 | ClassDef(DrawHitsDigits,0); | |
088f8e79 | 124 | }; |
125 | ||
126 | //____________________________________________________________________ | |
127 | // | |
128 | // EOF | |
129 | // |