CbmRoot
Loading...
Searching...
No Matches
CbmTrdRadiator.h
Go to the documentation of this file.
1/* Copyright (C) 2006-2021 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
2 SPDX-License-Identifier: GPL-3.0-only
3 Authors: Matus Kalisky [committer], Adrian Meyer-Ahrens, Cyrano Bergmann, Florian Uhlig */
4
5// -------------------------------------------------------------------------
6// ----- CbmTrdRadiator header file -----
7// ----- Created 10/11/04 by M.Kalisky -----
8// -------------------------------------------------------------------------
9
17#ifndef CBMTRDRADIATOR_H
18#define CBMTRDRADIATOR_H
19#define NCOMPONENTS 10
20
21#include <Rtypes.h> // for THashConsistencyHolder, ClassDef
22#include <RtypesCore.h> // for Float_t, Int_t, Bool_t, Double_t
23#include <TString.h> // for TString
24#include <TVector3.h> // for TVector3
25
26class CbmTrdPoint;
27class TH1D;
28typedef Float_t (*GetMuPtr)(Float_t);
29class CbmTrdRadiator {
30
31public:
32 struct CbmTrdEntranceWindow {
33 CbmTrdEntranceWindow(TString def = "");
34 Bool_t Init(TString def);
35
36 Int_t fN; // no of components
37 TString fMat[NCOMPONENTS]; // components' names
38 Float_t fThick[NCOMPONENTS]; // thickness [cm] of window components
39 Float_t fDens[NCOMPONENTS]; // density [g/cm^3] of window components
40 GetMuPtr GetMu[NCOMPONENTS]; // array of mu=f(E) for each component
41 };
42
44 CbmTrdRadiator(Bool_t SimpleTR = true, Int_t Nfoils = 337, Float_t FoilThick = 0.0012, Float_t GapThick = 0.09,
45 TString material = "pefoam20", TString window = "Kapton");
46
48 CbmTrdRadiator(Bool_t SimpleTR, TString prototype, TString window = "Kapton");
49 /* // implemented prototypes are:
50 "tdr18" 153x2mm PE foam foil layer (146 layers in box, 7 layers in carbon grid)
51 "B++" POKALON 24µm, 0.7mm gap, 350 foils
52 "K++" POKALON 24µm, 0.7mm gap, 350 foils micro-structured
53 "G30" 30xfiber layer
54
55
56 For the entrance window the general syntax would be c0;c1;c2;...
57 where ci are planar components. e.g. Al;C;Air;C;Al is the TRD-2D sandwich
58 There are two shortcuts for the Munster detector namely:
59 "Mylar" - simple mylar window
60 "Kapton" - Al(50nm)kapton(25µm )
61 */
62
64 virtual ~CbmTrdRadiator();
65
67 void CreateHistograms();
68
70 void Init();
71
73 void ProduceSpectra();
74
76 Int_t ELoss(Int_t index);
77
79 void ProcessTR();
80
82 Int_t TRspectrum();
83
85 Int_t WinTRspectrum();
86
88 Int_t DetTRspectrum();
89
91 Float_t GammaF();
92
94 Float_t Sigma(Float_t energykeV);
95
97 Float_t SigmaWin(Float_t energykeV);
98
100 Float_t SigmaDet(Float_t energykeV);
101
102 TString fWindowFoil;
103
105 static Float_t GetMuAl(Float_t energyMeV);
106 static Float_t GetMuPo(Float_t energyMeV);
107 static Float_t GetMuPok(Float_t energyMeV);
108 static Float_t GetMuKa(Float_t energyMeV);
109 static Float_t GetMuAir(Float_t energyMeV);
110 static Float_t GetMuXe(Float_t energyMeV);
111 static Float_t GetMuCO2(Float_t energyMeV);
112 static Float_t GetMuMy(Float_t energyMeV);
113 static Float_t GetMuC(Float_t energyMeV);
114
116 Float_t GetTR(TVector3 mom);
117
119 static Float_t Interpolate(Float_t energyMeV, const Float_t* en, const Float_t* mu, Int_t n);
120
122 static Int_t Locate(const Float_t* xv, Int_t n, Float_t xval, Int_t& kl, Float_t& dx);
123
124
126 void SetNFoils(Int_t n) { fNFoils = n; }
127 void SetFoilThick(Float_t t) { fFoilThick = t; }
128 void SetGapThick(Float_t t) { fGapThick = t; }
129 void SetSigma(Int_t SigmaT);
131 void SetEWwidths(Int_t n, Float_t* w);
132
133 Bool_t LatticeHit(const CbmTrdPoint* point);
134
135private:
136 CbmTrdRadiator& operator=(const CbmTrdRadiator&);
137 CbmTrdRadiator(const CbmTrdRadiator&);
138
139 /* Input parameters to be set */
140 TString fRadType;
141 Int_t fDetType; // 0: GSI 1: Muenster-Bucarest
142 Bool_t fFirstPass; // used for MB-Chamber with two gas gaps
143 Bool_t fSimpleTR; // mode of the TR production
144 Int_t fNFoils; // Number of foils in the radiator stack
145 Float_t fFoilThick; // Thickness of the foils (cm)
146 Float_t fGapThick; // Thickness of gaps between the foils (cm)
147 TString fFoilMaterial;
148 Float_t fGasThick; // Thickness of the active gas volume.
149
150
151 /* Parameters fixed in Init() */
152
153 Float_t fFoilDens; // Density of the radiator foils (g/cm^3)
154 Float_t fGapDens; // Dens. of gas in the radiator gaps (g/cm^3)
155 Float_t fFoilOmega; // Plasma frequency of the radiator foils
156 Float_t fGapOmega; // Plasma freq. of gas
157 Float_t fnPhotonCorr; // Correction of number of photon production efficiency
158
159 /* Parameters after correction of the angle of the particle */
160
161 Float_t fFoilThickCorr;
162 Float_t fGapThickCorr;
163 Float_t fGasThickCorr;
164
165
166 Float_t fnTRprod; // <nTR> produced
167 // Int_t fSigmaT; // set the absorption material
168
169 Float_t fWinDens;
170 Float_t fWinThick;
171 CbmTrdEntranceWindow WINDOW; // generic window structure
172
173 Float_t fCom1; // first component of the gas
174 Float_t fCom2; // second component of the gas
175
176 static const Int_t fSpNBins = 100; // Number of Bins for histos
177 static const Int_t fSpRange = 50; // Maximum (keV) of eloss spectrum
178 Float_t fSpBinWidth; // Bin width=fSpNBins/fSpRange
179
180 Float_t* fSigma;
181 Float_t* fSigmaWin;
182 Float_t* fSigmaDet;
183
184 TH1D* fSpectrum;
185 TH1D* fWinSpectrum;
186 TH1D* fDetSpectrumA;
187 TH1D* fDetSpectrum;
188
189 static const Int_t fNMom = 14; // number of momentum spectra
190 Double_t* fTrackMomentum;
191 // are available
192
193 TH1D* fFinal[fNMom];
194 Float_t fnTRabs[fNMom]; // <nTR> absorbed for differnt momenta
195 Float_t fnTRab; // <nTR> absorbed
196
197 Float_t fELoss; // The real result - in GeV
198
199 TVector3 fMom; // momentum of the electron
200
201 /* Photon absorption cross section for pokalon N470 */
202 static const Int_t fKN_pok = 46;
203 static constexpr Float_t fEn_pok[46] = {
204 1.000E-03, 1.500E-03, 2.000E-03, 3.000E-03, 4.000E-03, 5.000E-03, 6.000E-03, 8.000E-03, 1.000E-02, 1.500E-02,
205 2.000E-02, 3.000E-02, 4.000E-02, 5.000E-02, 6.000E-02, 8.000E-02, 1.000E-01, 1.500E-01, 2.000E-01, 3.000E-01,
206 4.000E-01, 5.000E-01, 6.000E-01, 8.000E-01, 1.000E+00, 1.022E+00, 1.250E+00, 1.500E+00, 2.000E+00, 2.044E+00,
207 3.000E+00, 4.000E+00, 5.000E+00, 6.000E+00, 7.000E+00, 8.000E+00, 9.000E+00, 1.000E+01, 1.100E+01, 1.200E+01,
208 1.300E+01, 1.400E+01, 1.500E+01, 1.600E+01, 1.800E+01, 2.000E+01};
209 static constexpr Float_t fMu_pok[46] = {
210 2.537E+03, 8.218E+02, 3.599E+02, 1.093E+02, 4.616E+01, 2.351E+01, 1.352E+01, 5.675E+00, 2.938E+00, 9.776E-01,
211 5.179E-01, 2.847E-01, 2.249E-01, 2.003E-01, 1.866E-01, 1.705E-01, 1.600E-01, 1.422E-01, 1.297E-01, 1.125E-01,
212 1.007E-01, 9.193E-02, 8.501E-02, 7.465E-02, 6.711E-02, 6.642E-02, 6.002E-02, 5.463E-02, 4.686E-02, 4.631E-02,
213 3.755E-02, 3.210E-02, 2.851E-02, 2.597E-02, 2.409E-02, 2.263E-02, 2.149E-02, 2.056E-02, 1.979E-02, 1.916E-02,
214 1.862E-02, 1.816E-02, 1.777E-02, 1.743E-02, 1.687E-02, 1.644E-02};
215
216 /* Photon absorption cross section for kapton */
217 static const Int_t fKN_ka = 46;
218 static constexpr Float_t fEn_ka[46] = {
219 1.000E-03, 1.500E-03, 2.000E-03, 3.000E-03, 4.000E-03, 5.000E-03, 6.000E-03, 8.000E-03, 1.000E-02, 1.500E-02,
220 2.000E-02, 3.000E-02, 4.000E-02, 5.000E-02, 6.000E-02, 8.000E-02, 1.000E-01, 1.500E-01, 2.000E-01, 3.000E-01,
221 4.000E-01, 5.000E-01, 6.000E-01, 8.000E-01, 1.000E+00, 1.022E+00, 1.250E+00, 1.500E+00, 2.000E+00, 2.044E+00,
222 3.000E+00, 4.000E+00, 5.000E+00, 6.000E+00, 7.000E+00, 8.000E+00, 9.000E+00, 1.000E+01, 1.100E+01, 1.200E+01,
223 1.300E+01, 1.400E+01, 1.500E+01, 1.600E+01, 1.800E+01, 2.000E+01};
224 static constexpr Float_t fMu_ka[46] = {
225 2.731E+03, 8.875E+02, 3.895E+02, 1.185E+02, 5.013E+01, 2.555E+01, 1.470E+01, 6.160E+00, 3.180E+00, 1.043E+00,
226 5.415E-01, 2.880E-01, 2.235E-01, 1.973E-01, 1.830E-01, 1.666E-01, 1.560E-01, 1.385E-01, 1.263E-01, 1.095E-01,
227 9.799E-02, 8.944E-02, 8.270E-02, 7.262E-02, 6.529E-02, 6.462E-02, 5.839E-02, 5.315E-02, 4.561E-02, 4.507E-02,
228 3.659E-02, 3.133E-02, 2.787E-02, 2.543E-02, 2.362E-02, 2.223E-02, 2.114E-02, 2.025E-02, 1.953E-02, 1.893E-02,
229 1.842E-02, 1.799E-02, 1.762E-02, 1.730E-02, 1.678E-02, 1.638E-02};
230
231 /* Photon absorption cross section for aluminum */
232 static const Int_t fKN_al = 48;
233 static constexpr Float_t fEn_al[48] = {
234 1.000E-03, 1.500E-03, 1.560E-03, 1.560E-03, 2.000E-03, 3.000E-03, 4.000E-03, 5.000E-03, 6.000E-03, 8.000E-03,
235 1.000E-02, 1.500E-02, 2.000E-02, 3.000E-02, 4.000E-02, 5.000E-02, 6.000E-02, 8.000E-02, 1.000E-01, 1.500E-01,
236 2.000E-01, 3.000E-01, 4.000E-01, 5.000E-01, 6.000E-01, 8.000E-01, 1.000E+00, 1.022E+00, 1.250E+00, 1.500E+00,
237 2.000E+00, 2.044E+00, 3.000E+00, 4.000E+00, 5.000E+00, 6.000E+00, 7.000E+00, 8.000E+00, 9.000E+00, 1.000E+01,
238 1.100E+01, 1.200E+01, 1.300E+01, 1.400E+01, 1.500E+01, 1.600E+01, 1.800E+01, 2.000E+01};
239 static constexpr Float_t fMu_al[48] = {
240 1.185E+03, 4.023E+02, 3.621E+02, 3.957E+03, 2.263E+03, 7.881E+02, 3.605E+02, 1.934E+02, 1.153E+02, 5.032E+01,
241 2.621E+01, 7.955E+00, 3.442E+00, 1.128E+00, 5.684E-01, 3.681E-01, 2.778E-01, 2.018E-01, 1.704E-01, 1.378E-01,
242 1.223E-01, 1.042E-01, 9.276E-02, 8.445E-02, 7.802E-02, 6.841E-02, 6.146E-02, 6.080E-02, 5.496E-02, 5.006E-02,
243 4.324E-02, 4.277E-02, 3.541E-02, 3.106E-02, 2.836E-02, 2.655E-02, 2.529E-02, 2.437E-02, 2.369E-02, 2.318E-02,
244 2.279E-02, 2.249E-02, 2.226E-02, 2.208E-02, 2.195E-02, 2.185E-02, 2.173E-02, 2.168E-02};
245
246 /* Photon absorption cross section for polypropylene/polyethylene */
247 static const Int_t fKN_po = 36;
248 static constexpr Float_t fEn_po[36] = {
249 1.000E-03, 1.500E-03, 2.000E-03, 3.000E-03, 4.000E-03, 5.000E-03, 6.000E-03, 8.000E-03, 1.000E-02,
250 1.500E-02, 2.000E-02, 3.000E-02, 4.000E-02, 5.000E-02, 6.000E-02, 8.000E-02, 1.000E-01, 1.500E-01,
251 2.000E-01, 3.000E-01, 4.000E-01, 5.000E-01, 6.000E-01, 8.000E-01, 1.000E+00, 1.250E+00, 1.500E+00,
252 2.000E+00, 3.000E+00, 4.000E+00, 5.000E+00, 6.000E+00, 8.000E+00, 1.000E+01, 1.500E+01, 2.000E+01};
253 static constexpr Float_t fMu_po[36] = {
254 1.894E+03, 5.999E+02, 2.593E+02, 7.743E+01, 3.242E+01, 1.643E+01, 9.432E+00, 3.975E+00, 2.088E+00,
255 7.452E-01, 4.315E-01, 2.706E-01, 2.275E-01, 2.084E-01, 1.970E-01, 1.823E-01, 1.719E-01, 1.534E-01,
256 1.402E-01, 1.217E-01, 1.089E-01, 9.947E-02, 9.198E-02, 8.078E-02, 7.262E-02, 6.495E-02, 5.910E-02,
257 5.064E-02, 4.045E-02, 3.444E-02, 3.045E-02, 2.760E-02, 2.383E-02, 2.145E-02, 1.819E-02, 1.658E-02};
258
259 /* Photon absorption cross section for carbon */
260 static const Int_t fKN_c = 46;
261 static constexpr Float_t fEn_c[46] = {
262 1.000E-03, 1.500E-03, 2.000E-03, 3.000E-03, 4.000E-03, 5.000E-03, 6.000E-03, 8.000E-03, 1.000E-02, 1.500E-02,
263 2.000E-02, 3.000E-02, 4.000E-02, 5.000E-02, 6.000E-02, 8.000E-02, 1.000E-01, 1.500E-01, 2.000E-01, 3.000E-01,
264 4.000E-01, 5.000E-01, 6.000E-01, 8.000E-01, 1.000E+00, 1.022E+00, 1.250E+00, 1.500E+00, 2.000E+00, 2.044E+00,
265 3.000E+00, 4.000E+00, 5.000E+00, 6.000E+00, 7.000E+00, 8.000E+00, 9.000E+00, 1.000E+01, 1.100E+01, 1.200E+01,
266 1.300E+01, 1.400E+01, 1.500E+01, 1.600E+01, 1.800E+01, 2.000E+01};
267 static constexpr Float_t fMu_c[46] = {
268 2.211e+03, 7.004e+02, 3.026e+02, 9.032e+01, 3.778e+01, 1.912e+01, 1.095e+01, 4.576e+00, 2.373e+00, 8.074e-01,
269 4.420e-01, 2.562e-01, 2.076e-01, 1.871e-01, 1.753e-01, 1.610e-01, 1.514e-01, 1.347e-01, 1.229e-01, 1.066e-01,
270 9.547e-02, 8.715e-02, 8.058e-02, 7.076e-02, 6.362e-02, 6.296e-02, 5.690e-02, 5.179e-02, 4.443e-02, 4.391e-02,
271 3.563e-02, 3.047e-02, 2.708e-02, 2.469e-02, 2.291e-02, 2.154e-02, 2.047e-02, 1.959e-02, 1.888e-02, 1.828e-02,
272 1.778e-02, 1.735e-02, 1.698e-02, 1.667e-02, 1.615e-02, 1.575e-02};
273
274 /* Photon absorption cross section for air */
275 static const Int_t fKN_air = 38;
276 static constexpr Float_t fEn_air[38] = {
277 0.10000E-02, 0.15000E-02, 0.20000E-02, 0.30000E-02, 0.32029E-02, 0.32029E-02, 0.40000E-02, 0.50000E-02,
278 0.60000E-02, 0.80000E-02, 0.10000E-01, 0.15000E-01, 0.20000E-01, 0.30000E-01, 0.40000E-01, 0.50000E-01,
279 0.60000E-01, 0.80000E-01, 0.10000E+00, 0.15000E+00, 0.20000E+00, 0.30000E+00, 0.40000E+00, 0.50000E+00,
280 0.60000E+00, 0.80000E+00, 0.10000E+01, 0.12500E+01, 0.15000E+01, 0.20000E+01, 0.30000E+01, 0.40000E+01,
281 0.50000E+01, 0.60000E+01, 0.80000E+01, 0.10000E+02, 0.15000E+02, 0.20000E+02};
282 static constexpr Float_t fMu_air[38] = {
283 0.35854E+04, 0.11841E+04, 0.52458E+03, 0.16143E+03, 0.15722E+03, 0.14250E+03, 0.77538E+02, 0.40099E+02,
284 0.23313E+02, 0.98816E+01, 0.51000E+01, 0.16079E+01, 0.77536E+00, 0.35282E+00, 0.24790E+00, 0.20750E+00,
285 0.18703E+00, 0.16589E+00, 0.15375E+00, 0.13530E+00, 0.12311E+00, 0.10654E+00, 0.95297E-01, 0.86939E-01,
286 0.80390E-01, 0.70596E-01, 0.63452E-01, 0.56754E-01, 0.51644E-01, 0.44382E-01, 0.35733E-01, 0.30721E-01,
287 0.27450E-01, 0.25171E-01, 0.22205E-01, 0.20399E-01, 0.18053E-01, 0.18057E-01};
288
289 /* Photon absorption cross section for xenon */
290 static const Int_t fKN_xe = 48;
291 static constexpr Float_t fEn_xe[48] = {
292 1.00000E-03, 1.07191E-03, 1.14900E-03, 1.14900E-03, 1.50000E-03, 2.00000E-03, 3.00000E-03, 4.00000E-03,
293 4.78220E-03, 4.78220E-03, 5.00000E-03, 5.10370E-03, 5.10370E-03, 5.27536E-03, 5.45280E-03, 5.45280E-03,
294 6.00000E-03, 8.00000E-03, 1.00000E-02, 1.50000E-02, 2.00000E-02, 3.00000E-02, 3.45614E-02, 3.45614E-02,
295 4.00000E-02, 5.00000E-02, 6.00000E-02, 8.00000E-02, 1.00000E-01, 1.50000E-01, 2.00000E-01, 3.00000E-01,
296 4.00000E-01, 5.00000E-01, 6.00000E-01, 8.00000E-01, 1.00000E+00, 1.25000E+00, 1.50000E+00, 2.00000E+00,
297 3.00000E+00, 4.00000E+00, 5.00000E+00, 6.00000E+00, 8.00000E+00, 1.00000E+01, 1.50000E+01, 2.00000E+01};
298 static constexpr Float_t fMu_xe[48] = {
299 9.413E+03, 8.151E+03, 7.035E+03, 7.338E+03, 4.085E+03, 2.088E+03, 7.780E+02, 3.787E+02, 2.408E+02, 6.941E+02,
300 6.392E+02, 6.044E+02, 8.181E+02, 7.579E+02, 6.991E+02, 8.064E+02, 6.376E+02, 3.032E+02, 1.690E+02, 5.743E+01,
301 2.652E+01, 8.930E+00, 6.129E+00, 3.316E+01, 2.270E+01, 1.272E+01, 7.825E+00, 3.633E+00, 2.011E+00, 7.202E-01,
302 3.760E-01, 1.797E-01, 1.223E-01, 9.699E-02, 8.281E-02, 6.696E-02, 5.785E-02, 5.054E-02, 4.594E-02, 4.078E-02,
303 3.681E-02, 3.577E-02, 3.583E-02, 3.634E-02, 3.797E-02, 3.987E-02, 4.445E-02, 4.815E-02};
304
305 /* Photon absorption cross section for CO2 */
306 static const Int_t fKN_co2 = 36;
307 static constexpr Float_t fEn_co2[36] = {0.10000E-02, 0.15000E-02, 0.20000E-02, 0.30000E-02, 0.40000E-02, 0.50000E-02,
308 0.60000E-02, 0.80000E-02, 0.10000E-01, 0.15000E-01, 0.20000E-01, 0.30000E-01,
309 0.40000E-01, 0.50000E-01, 0.60000E-01, 0.80000E-01, 0.10000E+00, 0.15000E+00,
310 0.20000E+00, 0.30000E+00, 0.40000E+00, 0.50000E+00, 0.60000E+00, 0.80000E+00,
311 0.10000E+01, 0.12500E+01, 0.15000E+01, 0.20000E+01, 0.30000E+01, 0.40000E+01,
312 0.50000E+01, 0.60000E+01, 0.80000E+01, 0.10000E+02, 0.15000E+02, 0.20000E+02};
313 static constexpr Float_t fMu_co2[36] = {0.39383E+04, 0.13166E+04, 0.58750E+03, 0.18240E+03, 0.77996E+02, 0.40024E+02,
314 0.23116E+02, 0.96997E+01, 0.49726E+01, 0.15543E+01, 0.74915E+00, 0.34442E+00,
315 0.24440E+00, 0.20589E+00, 0.18632E+00, 0.16578E+00, 0.15394E+00, 0.13558E+00,
316 0.12336E+00, 0.10678E+00, 0.95510E-01, 0.87165E-01, 0.80587E-01, 0.70769E-01,
317 0.63626E-01, 0.56894E-01, 0.51782E-01, 0.44499E-01, 0.35839E-01, 0.30825E-01,
318 0.27555E-01, 0.25269E-01, 0.22311E-01, 0.20516E-01, 0.18184E-01, 0.17152E-01};
319
320 /* Photon absorption cross section for mylar */
321 static const Int_t fKN_my = 36;
322 static constexpr Float_t fEn_my[36] = {1.00000E-03, 1.50000E-03, 2.00000E-03, 3.00000E-03, 4.00000E-03, 5.00000E-03,
323 6.00000E-03, 8.00000E-03, 1.00000E-02, 1.50000E-02, 2.00000E-02, 3.00000E-02,
324 4.00000E-02, 5.00000E-02, 6.00000E-02, 8.00000E-02, 1.00000E-01, 1.50000E-01,
325 2.00000E-01, 3.00000E-01, 4.00000E-01, 5.00000E-01, 6.00000E-01, 8.00000E-01,
326 1.00000E+00, 1.25000E+00, 1.50000E+00, 2.00000E+00, 3.00000E+00, 4.00000E+00,
327 5.00000E+00, 6.00000E+00, 8.00000E+00, 1.00000E+01, 1.50000E+01, 2.00000E+01};
328 static constexpr Float_t fMu_my[36] = {
329 2.911E+03, 9.536E+02, 4.206E+02, 1.288E+02, 5.466E+01, 2.792E+01, 1.608E+01, 6.750E+00, 3.481E+00,
330 1.132E+00, 5.798E-01, 3.009E-01, 2.304E-01, 2.020E-01, 1.868E-01, 1.695E-01, 1.586E-01, 1.406E-01,
331 1.282E-01, 1.111E-01, 9.947E-02, 9.079E-02, 8.395E-02, 7.372E-02, 6.628E-02, 5.927E-02, 5.395E-02,
332 4.630E-02, 3.715E-02, 3.181E-02, 2.829E-02, 2.582E-02, 2.257E-02, 2.057E-02, 1.789E-02, 1.664E-02};
333
334
335 ClassDef(CbmTrdRadiator, 3)
336};
337#endif
GetMuPtr
Float_t(* GetMuPtr)(Float_t)
Definition CbmTrdRadiator.h:28
NCOMPONENTS
#define NCOMPONENTS
Definition CbmTrdRadiator.h:19
CbmTrdRadiator::fMu_c
static constexpr Float_t fMu_c[46]
Definition CbmTrdRadiator.h:267
CbmTrdRadiator::fKN_xe
static const Int_t fKN_xe
Definition CbmTrdRadiator.h:290
CbmTrdRadiator::fMu_ka
static constexpr Float_t fMu_ka[46]
Definition CbmTrdRadiator.h:224
CbmTrdRadiator::SetFoilThick
void SetFoilThick(Float_t t)
Definition CbmTrdRadiator.h:127
CbmTrdRadiator::fEn_co2
static constexpr Float_t fEn_co2[36]
Definition CbmTrdRadiator.h:307
CbmTrdRadiator::Locate
static Int_t Locate(const Float_t *xv, Int_t n, Float_t xval, Int_t &kl, Float_t &dx)
Definition CbmTrdRadiator.cxx:1087
CbmTrdRadiator::SetGapThick
void SetGapThick(Float_t t)
Definition CbmTrdRadiator.h:128
CbmTrdRadiator::LatticeHit
Bool_t LatticeHit(const CbmTrdPoint *point)
Definition CbmTrdRadiator.cxx:435
CbmTrdRadiator::SetNFoils
void SetNFoils(Int_t n)
Definition CbmTrdRadiator.h:126
CbmTrdRadiator::CbmTrdRadiator
CbmTrdRadiator(const CbmTrdRadiator &)
CbmTrdRadiator::SigmaDet
Float_t SigmaDet(Float_t energykeV)
Definition CbmTrdRadiator.cxx:949
CbmTrdRadiator::fNMom
static const Int_t fNMom
TR passed through Detector.
Definition CbmTrdRadiator.h:189
CbmTrdRadiator::GetMuKa
static Float_t GetMuKa(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:982
CbmTrdRadiator::CbmTrdRadiator
CbmTrdRadiator(Bool_t SimpleTR=true, Int_t Nfoils=337, Float_t FoilThick=0.0012, Float_t GapThick=0.09, TString material="pefoam20", TString window="Kapton")
Definition CbmTrdRadiator.cxx:35
CbmTrdRadiator::Sigma
Float_t Sigma(Float_t energykeV)
Definition CbmTrdRadiator.cxx:889
CbmTrdRadiator::fEn_c
static constexpr Float_t fEn_c[46]
Definition CbmTrdRadiator.h:261
CbmTrdRadiator::GetMuC
static Float_t GetMuC(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:1002
CbmTrdRadiator::GetMuPo
static Float_t GetMuPo(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:1012
CbmTrdRadiator::SetSigma
void SetSigma(Int_t SigmaT)
Definition CbmTrdRadiator.cxx:852
CbmTrdRadiator::fEn_my
static constexpr Float_t fEn_my[36]
Definition CbmTrdRadiator.h:322
CbmTrdRadiator::fnTRabs
Float_t fnTRabs[fNMom]
Absorption spectra for different momenta.
Definition CbmTrdRadiator.h:194
CbmTrdRadiator::fSigmaDet
Float_t * fSigmaDet
[fSpNBins] Array of sigma values for the entrance window of detector
Definition CbmTrdRadiator.h:182
CbmTrdRadiator::fKN_pok
static const Int_t fKN_pok
Definition CbmTrdRadiator.h:202
CbmTrdRadiator::GetMuXe
static Float_t GetMuXe(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:1032
CbmTrdRadiator::fKN_co2
static const Int_t fKN_co2
Definition CbmTrdRadiator.h:306
CbmTrdRadiator::GetMuAir
static Float_t GetMuAir(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:1022
CbmTrdRadiator::fMu_my
static constexpr Float_t fMu_my[36]
Definition CbmTrdRadiator.h:328
CbmTrdRadiator::fEn_pok
static constexpr Float_t fEn_pok[46]
Definition CbmTrdRadiator.h:203
CbmTrdRadiator::WINDOW
CbmTrdEntranceWindow WINDOW
Definition CbmTrdRadiator.h:171
CbmTrdRadiator::fSigmaWin
Float_t * fSigmaWin
[fSpNBins] Array of sigma values for the foil of the radiator
Definition CbmTrdRadiator.h:181
CbmTrdRadiator::fKN_po
static const Int_t fKN_po
Definition CbmTrdRadiator.h:247
CbmTrdRadiator::fMu_al
static constexpr Float_t fMu_al[48]
Definition CbmTrdRadiator.h:239
CbmTrdRadiator::fMu_pok
static constexpr Float_t fMu_pok[46]
Definition CbmTrdRadiator.h:209
CbmTrdRadiator::Interpolate
static Float_t Interpolate(Float_t energyMeV, const Float_t *en, const Float_t *mu, Int_t n)
Definition CbmTrdRadiator.cxx:1062
CbmTrdRadiator::SigmaWin
Float_t SigmaWin(Float_t energykeV)
Definition CbmTrdRadiator.cxx:923
CbmTrdRadiator::fMu_co2
static constexpr Float_t fMu_co2[36]
Definition CbmTrdRadiator.h:313
CbmTrdRadiator::fDetSpectrumA
TH1D * fDetSpectrumA
TR spectra in gas-window foil.
Definition CbmTrdRadiator.h:186
CbmTrdRadiator::fDetSpectrum
TH1D * fDetSpectrum
TR absorbed in Detector.
Definition CbmTrdRadiator.h:187
CbmTrdRadiator::fKN_ka
static const Int_t fKN_ka
Definition CbmTrdRadiator.h:217
CbmTrdRadiator::GetTR
Float_t GetTR(TVector3 mom)
Definition CbmTrdRadiator.cxx:547
CbmTrdRadiator::fTrackMomentum
Double_t * fTrackMomentum
Definition CbmTrdRadiator.h:190
CbmTrdRadiator::fEn_ka
static constexpr Float_t fEn_ka[46]
Definition CbmTrdRadiator.h:218
CbmTrdRadiator::fSpNBins
static const Int_t fSpNBins
Definition CbmTrdRadiator.h:176
CbmTrdRadiator::fSpectrum
TH1D * fSpectrum
[fSpNBins] Array of sigma values for the active gas
Definition CbmTrdRadiator.h:184
CbmTrdRadiator::ELoss
Int_t ELoss(Int_t index)
Definition CbmTrdRadiator.cxx:629
CbmTrdRadiator::fWinSpectrum
TH1D * fWinSpectrum
TR photon energy spectrum.
Definition CbmTrdRadiator.h:185
CbmTrdRadiator::fMu_po
static constexpr Float_t fMu_po[36]
Definition CbmTrdRadiator.h:253
CbmTrdRadiator::fMu_air
static constexpr Float_t fMu_air[38]
Definition CbmTrdRadiator.h:282
CbmTrdRadiator::SetEWwidths
void SetEWwidths(Int_t n, Float_t *w)
Definition CbmTrdRadiator.cxx:424
CbmTrdRadiator::fSpRange
static const Int_t fSpRange
Definition CbmTrdRadiator.h:177
CbmTrdRadiator::GetMuAl
static Float_t GetMuAl(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:992
CbmTrdRadiator::operator=
CbmTrdRadiator & operator=(const CbmTrdRadiator &)
CbmTrdRadiator::fEn_al
static constexpr Float_t fEn_al[48]
Definition CbmTrdRadiator.h:233
CbmTrdRadiator::fKN_c
static const Int_t fKN_c
Definition CbmTrdRadiator.h:260
CbmTrdRadiator::fEn_po
static constexpr Float_t fEn_po[36]
Definition CbmTrdRadiator.h:248
CbmTrdRadiator::fMu_xe
static constexpr Float_t fMu_xe[48]
Definition CbmTrdRadiator.h:298
CbmTrdRadiator::fKN_air
static const Int_t fKN_air
Definition CbmTrdRadiator.h:275
CbmTrdRadiator::GetMuMy
static Float_t GetMuMy(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:1052
CbmTrdRadiator::~CbmTrdRadiator
virtual ~CbmTrdRadiator()
Definition CbmTrdRadiator.cxx:131
CbmTrdRadiator::fEn_xe
static constexpr Float_t fEn_xe[48]
Definition CbmTrdRadiator.h:291
CbmTrdRadiator::GetMuCO2
static Float_t GetMuCO2(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:1042
CbmTrdRadiator::GetMuPok
static Float_t GetMuPok(Float_t energyMeV)
Definition CbmTrdRadiator.cxx:973
CbmTrdRadiator::fKN_al
static const Int_t fKN_al
Definition CbmTrdRadiator.h:232
CbmTrdRadiator::fEn_air
static constexpr Float_t fEn_air[38]
Definition CbmTrdRadiator.h:276
CbmTrdRadiator::fKN_my
static const Int_t fKN_my
Definition CbmTrdRadiator.h:321
CbmTrdRadiator::fFinal
TH1D * fFinal[fNMom]
[fNMom] Track momenta for which spectra
Definition CbmTrdRadiator.h:193