HYDRA_development_version
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros
henergylosscorrpar.cc
Go to the documentation of this file.
1 //*-- AUTHOR : Olga Pechenova
2 //*-- Modified : 2006 by Olga Pechenova
3 //*-- Modified : 2013 by Olga Pechenova
4 //*-- Modified : 2015 by Vladimir Pechenov
5 //*-- Modified : 2019 by Olga Pechenova
6 //
7 //_HADES_CLASS_DESCRIPTION
8 /////////////////////////////////////////////////////////////////////////
9 //
10 // HEnergyLossCorrPar
11 //
12 // Container class keep parameters for energy loss correction
13 //
14 //
15 // Double_t getDeltaMom(Int_t pId, Double_t mom, Double_t theta);
16 // Function getDeltaMom returns the correction value (delP) for momentum.
17 // Parameters: - pId is particle ID. If parameter pId is not correct
18 // function getDeltaMom returns 0.
19 // - mom is Runge-Kutta reconstructed momentum multiplied by particle charge (MeV).
20 // - theta is theta angle in the interval [15;87] (degrees).
21 // If parameter theta is out of this interval
22 // function getDeltaMom returns 0.
23 //
24 // Double_t getCorrMom(Int_t pId, Double_t mom, Double_t theta) {
25 // return mom+getDeltaMom(pId,mom,theta);}
26 // Function getCorrMom returns the corrected momentum.
27 // Parameters have the same meanings as for function getDeltaMom.
28 //
29 //
30 // Alternative method of getting parameters (without usage of runTimeDb):
31 //
32 // Bool_t setDefaultPar(const Char_t *rn);
33 // Function setDefaultPar sets default values to parameters
34 // for corresponding beam time.
35 //
36 // Example:
37 // HEnergyLossCorrPar enLossCorr;
38 // enLossCorr.setDefaultPar("jan04"); // or "nov02", or "aug04", or "apr07", or "apr12", or "mar19", or
39 // "jul14_w" - Wolfram target, "jul14_c3" - Carbon 3 segments, "jul14_pe" & "aug14_pe" - PE target, "aug14_c7" - Carbon 7 segments
40 // Double_t delP;
41 // ..............................
42 // beginning of your loop;
43 // ..............................................
44 // delP = enLossCorr.getDeltaMom(14, mom, theta);
45 // ..............................................
46 // end of your loop;
47 /////////////////////////////////////////////////////////////////////////
48 
49 #include "TMath.h"
50 #include "TString.h"
51 #include "henergylosscorrpar.h"
52 #include "hphysicsconstants.h"
53 #include "hparamlist.h"
54 
55 #include <stdlib.h>
56 #include <cmath>
57 
58 using namespace TMath;
59 
60 HEnergyLossCorrPar* HEnergyLossCorrPar::gEnergyLossCorrPar = NULL;
61 
62 ClassImp(HEnergyLossCorrPar)
63 
64 HEnergyLossCorrPar::HEnergyLossCorrPar(const Char_t* name,const Char_t* title, const Char_t* context) {
65  // constructor
66  momInt1=720.;
67  momInt2=2000.;
68  nThetaReg=7;
69  nParams=10;
70  thetaStep=10.;
71  thetaMiddle0=23.;
72  for(Int_t i=0;i<nThetaReg*nParams;i++) parMomCorrH[i]=0.;
73 // for electrons
74  nParamsElect=5;
75  for(Int_t i=0;i<nParamsElect;i++) parMomCorrElect[i]=0.;
76  typePar = 0;
77  gEnergyLossCorrPar = this;
78 }
79 
80 void HEnergyLossCorrPar::clear(void) {
81  for(Int_t i=0;i<nThetaReg*nParams;i++) parMomCorrH[i]=0.;
82  for(Int_t i=0;i<nParamsElect;i++) parMomCorrElect[i]=0.;
83  status=kFALSE;
84  resetInputVersions();
85  changed=kFALSE;
86  typePar = 0;
87 }
88 
89 void HEnergyLossCorrPar::putParams(HParamList* l) {
90  // puts all parameters to the parameter list, which is used by the io
91  if (!l) return;
92  l->add("momInt1", momInt1);
93  l->add("momInt2", momInt2);
94  l->add("nThetaReg", nThetaReg);
95  l->add("nParams", nParams);
96  l->add("parMomCorrH", parMomCorrH,nThetaReg*nParams);
97  l->add("thetaMiddle0", thetaMiddle0);
98  l->add("thetaStep", thetaStep);
99 
100  l->add("nParamsElect", nParamsElect);
101  l->add("parMomCorrElect", parMomCorrElect,nParamsElect);
102 }
103 
104 
105 Bool_t HEnergyLossCorrPar::getParams(HParamList* l) {
106  // gets all parameters from the parameter list, which is used by the io
107  if (!l) return kFALSE;
108  if (!(l->fill("momInt1", &momInt1))) return kFALSE;
109  if (!(l->fill("momInt2", &momInt2))) return kFALSE;
110  if (!(l->fill("nThetaReg", &nThetaReg))) return kFALSE;
111  if (!(l->fill("nParams", &nParams))) return kFALSE;
112  if (!(l->fill("parMomCorrH",parMomCorrH,nThetaReg*nParams))) return kFALSE;
113  if (!(l->fill("thetaMiddle0", &thetaMiddle0))) return kFALSE;
114  if (!(l->fill("thetaStep", &thetaStep))) return kFALSE;
115 
116  if (!(l->fill("nParamsElect", &nParamsElect))) return kFALSE;
117  if (!(l->fill("parMomCorrElect",parMomCorrElect,nParamsElect))) return kFALSE;
118  return kTRUE;
119 }
120 
121 Bool_t HEnergyLossCorrPar::setDefaultPar(TString run) {
122  run.ToLower();
123  if(run == "mar19") {
124  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
125  nParams = 12;
126  nThetaReg = 22; //from 20 grad till 86 grad
127  thetaStep = 3.;
128  thetaMiddle0 = 21.5;
129  momInt1 = 700;
130  momInt2 = 2500;
131  Double_t par[12*22] = {
132 
133  331.127, -1.79734, 0.00405562, -4.27346e-06, 1.735e-09, 73.6073, -0.00345576, -0.000582559, 4.63143, 3.88641, -0.0009463, 2.85763,
134  323.576, -1.79889, 0.00412357, -4.39188e-06, 1.79612e-09, 67.2306, -0.00350626, -0.000533308, 4.11317, 3.17558, -0.000925526, 2.50397,
135  319.744, -1.7855, 0.00410667, -4.38536e-06, 1.79703e-09, 65.5867, -0.00350781, -0.000513687, 3.98569, 3.24107, -0.000941404, 2.43316,
136  317.771, -1.78131, 0.00410794, -4.39378e-06, 1.80148e-09, 63.6918, -0.00350117, -0.000508197, 3.90963, 3.1139, -0.000937999, 2.38062,
137 
138  317.133, -1.78045, 0.00411114, -4.40261e-06, 1.80752e-09, 62.586, -0.00348684, -0.000496123, 3.86303, 3.16102, -0.000945441, 2.36431,
139  316.082, -1.77001, 0.00407816, -4.36077e-06, 1.78908e-09, 63.4748, -0.00350647, -0.000501702, 3.87875, 3.12845, -0.0009425, 2.3655,
140  318.695, -1.79788, 0.00417514, -4.50049e-06, 1.86117e-09, 63.9636, -0.00351873, -0.000503097, 3.87194, 3.06006, -0.000936491, 2.35765,
141  317.629, -1.78201, 0.00411396, -4.40713e-06, 1.81089e-09, 64.1952, -0.00351285, -0.000501238, 3.88635, 3.0865, -0.000937421, 2.37409,
142 
143  317.305, -1.78027, 0.00410967, -4.40112e-06, 1.8073e-09, 64.5274, -0.00352598, -0.00050951, 3.90579, 3.15457, -0.000944896, 2.37474,
144  307.19, -1.73526, 0.00402194, -4.31664e-06, 1.77418e-09, 58.3901, -0.00346572, -0.000454542, 3.61535, 3.04143, -0.000955256, 2.23254,
145  269.605, -1.56635, 0.00371216, -4.05734e-06, 1.69292e-09, 47.1203, -0.00348349, -0.000365663, 2.89917, 2.50292, -0.000975642, 1.78891,
146  252.504, -1.53331, 0.00376932, -4.25092e-06, 1.82238e-09, 37.7532, -0.00346699, -0.000292253, 2.33397, 2.38971, -0.00104616, 1.45284,
147 
148  253.178, -1.54433, 0.0038151, -4.32519e-06, 1.86432e-09, 38.5394, -0.00350577, -0.000296633, 2.33331, 2.31072, -0.00103731, 1.44315,
149  253.705, -1.54277, 0.00379995, -4.29587e-06, 1.84681e-09, 38.8674, -0.00350881, -0.000302755, 2.36416, 2.23024, -0.00102253, 1.45752,
150  255.392, -1.55332, 0.00382658, -4.32641e-06, 1.86004e-09, 39.5061, -0.00352028, -0.000308757, 2.39004, 2.3228, -0.00102958, 1.46958,
151  257.281, -1.56217, 0.00384362, -4.34149e-06, 1.86503e-09, 40.0745, -0.00351683, -0.000309021, 2.41602, 2.30394, -0.00102783, 1.4921,
152 
153  259.739, -1.57502, 0.00387274, -4.37386e-06, 1.87949e-09, 40.735, -0.00351943, -0.000314866, 2.45859, 2.4122, -0.00103648, 1.51779,
154  262.203, -1.58403, 0.00388233, -4.37164e-06, 1.87323e-09, 41.2454, -0.00350699, -0.00031877, 2.50649, 2.58409, -0.00105604, 1.55309,
155  264.451, -1.59371, 0.00389996, -4.38753e-06, 1.87931e-09, 42.3023, -0.00351954, -0.000330261, 2.56556, 2.61165, -0.00105409, 1.58334,
156  267.277, -1.6067, 0.00392585, -4.41284e-06, 1.88931e-09, 43.3361, -0.00351927, -0.000336194, 2.62386, 2.77194, -0.00107086, 1.62278,
157 
158  269.713, -1.6115, 0.00391444, -4.37416e-06, 1.86178e-09, 43.6594, -0.00348987, -0.000337886, 2.67942, 2.68782, -0.00105287, 1.66973,
159  273.884, -1.62928, 0.00393977, -4.38161e-06, 1.85577e-09, 46.0191, -0.00353234, -0.00035993, 2.78447, 2.38215, -0.000999647, 1.7164
160 
161  };
162  fillParMomCorrH(nThetaReg*nParams,par);
163 
164  //for electrons pol4 pol2 pol4
165  nParamsElect=13;
166  parMomCorrElect[0] = -9.13308;
167  parMomCorrElect[1] = 1.40618;
168  parMomCorrElect[2] = -0.0630954;
169  parMomCorrElect[3] = 0.00123199;
170  parMomCorrElect[4] = -0.00000884375;
171  parMomCorrElect[5] = 114.165;
172  parMomCorrElect[6] = -4.14175;
173  parMomCorrElect[7] = 0.0380465;
174  parMomCorrElect[8] = -4.75336;
175  parMomCorrElect[9] = 0.389213;
176  parMomCorrElect[10] = -0.00910909;
177  parMomCorrElect[11] = 0.0000934938;
178  parMomCorrElect[12] = -0.000000348053;
179 
180  typePar = 19;
181 
182  } else if(run == "jul14_w") {
183  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
184  nParams = 12;
185  nThetaReg = 22; //from 20 grad till 86 grad
186  thetaStep = 3.;
187  thetaMiddle0 = 21.5;
188  momInt1 = 700;
189  momInt2 = 2500;
190  Double_t par[12*22] = {
191  340.901,-1.50393,0.00278068,-0.00000238587,0.0000778523/100000.,119.233,-0.00334163,-0.00093208,8.25497,16.7969,-0.00141246,5.51394,
192  359.447,-1.76346,0.00366044,-0.00000356899,0.000134229/100000.,104.042,-0.0033692,-0.000844377,7.23206,19.212,-0.00150658,4.77071,
193  359.447,-1.76346,0.00366044,-0.00000356899,0.000134229/100000.,104.042,-0.0033692,-0.000844377,7.23206,19.212,-0.00150658,4.77071,
194  359.447,-1.76346,0.00366044,-0.00000356899,0.000134229/100000.,104.042,-0.0033692,-0.000844377,7.23206,19.212,-0.00150658,4.77071,
195  365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
196  365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
197  365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
198  365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
199  363.522,-1.69303,0.00336942,-0.00000317385,0.000116155/100000.,122.517,-0.00337433,-0.000959434,8.43587,29.2724,-0.00161355,5.63564,
200  363.522,-1.69303,0.00336942,-0.00000317385,0.000116155/100000.,122.517,-0.00337433,-0.000959434,8.43587,29.2724,-0.00161355,5.63564,
201  363.522,-1.69303,0.00336942,-0.00000317385,0.000116155/100000.,122.517,-0.00337433,-0.000959434,8.43587,29.2724,-0.00161355,5.63564,
202  379.355,-1.74201,0.00348639,-0.00000335548,0.000127136/100000.,141.13,-0.00337822,-0.00111639,9.76923,25.7175,-0.00151109,6.5268,
203  379.355,-1.74201,0.00348639,-0.00000335548,0.000127136/100000.,141.13,-0.00337822,-0.00111639,9.76923,25.7175,-0.00151109,6.5268,
204  379.355,-1.74201,0.00348639,-0.00000335548,0.000127136/100000.,141.13,-0.00337822,-0.00111639,9.76923,25.7175,-0.00151109,6.5268,
205  377.386,-1.61336,0.00297582,-0.00000259877,0.0000877779/100000.,147.441,-0.00325714,-0.00111654,10.7727,36.4515,-0.00158623,7.42574,
206  368.669,-1.47622,0.00249971,-0.00000193881,0.0000549969/100000.,167.991,-0.0033598,-0.00137807,12.0193,28.8549,-0.00149969,8.04572,
207  349.535,-1.24392,0.00174129,-0.00000092267,0.00000571075/100000.,166.974,-0.00325946,-0.00146554,12.9929,43.945,-0.00161646,8.77291,
208  375.433,-1.41188,0.00225412,-0.00000162605,0.0000408913/100000.,185.174,-0.00326443,-0.00147595,13.9235,35.1007,-0.00147076,9.54521,
209  367.516,-1.24152,0.00163618,-0.000000760779,-0.00000162891/100000.,202.942,-0.00326404,-0.00166685,15.5439,53.5905,-0.00159346,10.6089,
210  367.516,-1.24152,0.00163618,-0.000000760779,-0.00000162891/100000.,202.942,-0.00326404,-0.00166685,15.5439,53.5905,-0.00159346,10.6089,
211  410.091,-1.48901,0.00236275,-0.00000179646,0.0000538164/100000,211.017,-0.00321069,-0.00173945,16.8942,34.2646,-0.00141895,11.7299,
212  338.12,-0.750049,0.0000830358,0.00000112099,-0.0000814397/100000.,249.746,-0.00326935,-0.00196529,18.2082,36.0361,-0.0014173,12.3444
213  };
214  fillParMomCorrH(nThetaReg*nParams,par);
215 
216  //for electrons pol6
217  nParamsElect=7;
218  parMomCorrElect[0] = 36.8903;
219  parMomCorrElect[1] = -3.45134;
220  parMomCorrElect[2] = 0.140583;
221  parMomCorrElect[3] = -0.00259951;
222  parMomCorrElect[4] = 0.0000198486;
223  parMomCorrElect[5] = -0.000662151/100000.;
224  parMomCorrElect[6] = -0.0000403811/100000.;
225 
226  typePar = 14;
227  } else if(run == "jul14_c3") {
228  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
229  nParams = 12;
230  nThetaReg = 22; //from 20 grad till 86 grad
231  thetaStep = 3.;
232  thetaMiddle0 = 21.5;
233  momInt1 = 700;
234  momInt2 = 2500;
235  Double_t par[12*22] = {
236  363.813,-1.87234,0.00404871,-0.0000041178,0.00016226/100000.,90.8896,-0.00338883,-0.000731713,5.95128,5.35734,-0.000983874,3.73794,
237  340.39,-1.78571,0.00390809,-0.00000400218,0.000158149/100000.,82.3474,-0.00344345,-0.000638642,5.1968,5.85834,-0.00108304,3.2622,
238  340.39,-1.78571,0.00390809,-0.00000400218,0.000158149/100000.,82.3474,-0.00344345,-0.000638642,5.1968,5.85834,-0.00108304,3.2622,
239  340.39,-1.78571,0.00390809,-0.00000400218,0.000158149/100000.,82.3474,-0.00344345,-0.000638642,5.1968,5.85834,-0.00108304,3.2622,
240  339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296,
241  339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296,
242  339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296,
243  339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296,
244  345.581,-1.79395,0.00389273,-0.00000395957,0.000155694/100000.,84.0997,-0.00341154,-0.000663507,5.42551,5.65421,-0.00105822,3.41813,
245  345.581,-1.79395,0.00389273,-0.00000395957,0.000155694/100000.,84.0997,-0.00341154,-0.000663507,5.42551,5.65421,-0.00105822,3.41813,
246  345.581,-1.79395,0.00389273,-0.00000395957,0.000155694/100000.,84.0997,-0.00341154,-0.000663507,5.42551,5.65421,-0.00105822,3.41813,
247  351.345,-1.81068,0.00390588,-0.00000395676,0.000155283/100000.,88.5021,-0.00344,-0.000700565,5.64477,5.92098,-0.00106608,3.53526,
248  351.345,-1.81068,0.00390588,-0.00000395676,0.000155283/100000.,88.5021,-0.00344,-0.000700565,5.64477,5.92098,-0.00106608,3.53526,
249  351.345,-1.81068,0.00390588,-0.00000395676,0.000155283/100000.,88.5021,-0.00344,-0.000700565,5.64477,5.92098,-0.00106608,3.53526,
250  355.35,-1.82723,0.00394577,-0.00000401542,0.000158839/100000.,92.513,-0.00347466,-0.00074445,5.83361,6.01486,-0.00106656,3.61997,
251  358.09,-1.82336,0.00386925,-0.00000383628,0.000146582/100000.,89.2857,-0.00340005,-0.000704394,5.80144,6.47893,-0.00108469,3.66993,
252  363.225,-1.86139,0.00399638,-0.00000402954,0.00015732/100000.,94.2199,-0.00345562,-0.000745102,5.96009,6.5169,-0.00108865,3.7278,
253  367.319,-1.87481,0.00399567,-0.00000398468,0.000153338/100000.,92.8587,-0.00340985,-0.000728312,5.97965,6.31287,-0.00106943,3.78329,
254  365.256,-1.83977,0.00388251,-0.00000384856,0.000147824/100000.,96.1542,-0.00342691,-0.000763212,6.18846,6.87292,-0.00108557,3.89155,
255  365.256,-1.83977,0.00388251,-0.00000384856,0.000147824/100000.,96.1542,-0.00342691,-0.000763212,6.18846,6.87292,-0.00108557,3.89155,
256  377.967,-1.89942,0.0039926,-0.00000393773,0.000150343/100000.,99.3301,-0.00345412,-0.000839982,6.4699,6.84137,-0.00107652,3.98033,
257  377.967,-1.89942,0.0039926,-0.00000393773,0.000150343/100000.,99.3301,-0.00345412,-0.000839982,6.4699,6.84137,-0.00107652,3.98033
258  };
259  fillParMomCorrH(nThetaReg*nParams,par);
260 
261  //for electrons pol6
262  nParamsElect=7;
263  parMomCorrElect[0] = 19.0501;
264  parMomCorrElect[1] = -1.78541;
265  parMomCorrElect[2] = 0.0799932;
266  parMomCorrElect[3] = -0.00181062;
267  parMomCorrElect[4] = 0.0000220232;
268  parMomCorrElect[5] = -0.000000136325;
269  parMomCorrElect[6] = 0.0000334713/100000.;
270 
271  typePar = 14;
272  } else if(run == "jul14_pe" || run == "aug14_pe") {
273  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
274  nParams = 12;
275  nThetaReg = 22; //from 20 grad till 86 grad
276  thetaStep = 3.;
277  thetaMiddle0 = 21.5;
278  momInt1 = 700;
279  momInt2 = 2500;
280  Double_t par[12*22] = {
281  395.219,-1.94427,0.00406475,-0.00000403356,0.000156339/100000.,113.782,-0.00340005,-0.000931158,7.49854,6.7844,-0.000981211,4.68952,
282  366.057,-1.81668,0.00377552,-0.00000368103,0.000138881/100000.,101.211,-0.00345251,-0.000818066,6.47664,6.09025,-0.00101831,4.02327,
283  366.057,-1.81668,0.00377552,-0.00000368103,0.000138881/100000.,101.211,-0.00345251,-0.000818066,6.47664,6.09025,-0.00101831,4.02327,
284  366.057,-1.81668,0.00377552,-0.00000368103,0.000138881/100000.,101.211,-0.00345251,-0.000818066,6.47664,6.09025,-0.00101831,4.02327,
285  352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111,
286  352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111,
287  352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111,
288  352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111,
289  354.009,-1.86253,0.00408366,-0.00000419045,0.000166054/100000.,83.2562,-0.00343848,-0.000651417,5.28283,6.08154,-0.00110533,3.32809,
290  354.009,-1.86253,0.00408366,-0.00000419045,0.000166054/100000.,83.2562,-0.00343848,-0.000651417,5.28283,6.08154,-0.00110533,3.32809,
291  354.009,-1.86253,0.00408366,-0.00000419045,0.000166054/100000.,83.2562,-0.00343848,-0.000651417,5.28283,6.08154,-0.00110533,3.32809,
292  348.382,-1.8238,0.00396499,-0.0000040237,0.000157397/100000.,82.5202,-0.00345177,-0.000645806,5.19877,6.15786,-0.00112109,3.26955,
293  348.382,-1.8238,0.00396499,-0.0000040237,0.000157397/100000.,82.5202,-0.00345177,-0.000645806,5.19877,6.15786,-0.00112109,3.26955,
294  348.382,-1.8238,0.00396499,-0.0000040237,0.000157397/100000.,82.5202,-0.00345177,-0.000645806,5.19877,6.15786,-0.00112109,3.26955,
295  353.894,-1.87342,0.00412526,-0.00000424506,0.000168473/100000.,83.1282,-0.00346394,-0.000648254,5.20019,6.28797,-0.00113295,3.27027,
296  353.85,-1.86659,0.00409331,-0.00000419413,0.00016577/100000.,83.926,-0.00347166,-0.000660473,5.2443,6.19997,-0.00112843,3.28492,
297  351.858,-1.84545,0.00402381,-0.00000409948,0.000161127/100000.,84.8005,-0.00347937,-0.000657066,5.26311,6.31226,-0.00113253,3.30581,
298  356.568,-1.87939,0.00411911,-0.00000421744,0.000166505/100000.,86.025,-0.00349049,-0.00067037,5.32277,5.9717,-0.00111004,3.33111,
299  357.728,-1.87848,0.00410439,-0.00000419094,0.000165047/100000.,85.0026,-0.00344817,-0.000657249,5.3514,6.1102,-0.00111809,3.3878,
300  357.728,-1.87848,0.00410439,-0.00000419094,0.000165047/100000.,85.0026,-0.00344817,-0.000657249,5.3514,6.1102,-0.00111809,3.3878,
301  358.954,-1.877,0.00408277,-0.00000415016,0.000162742/100000.,78.3719,-0.0033126,-0.000622673,5.28907,7.45359,-0.00119627,3.41713,
302  358.954,-1.877,0.00408277,-0.00000415016,0.000162742/100000.,78.3719,-0.0033126,-0.000622673,5.28907,7.45359,-0.00119627,3.41713
303  };
304  fillParMomCorrH(nThetaReg*nParams,par);
305 
306  //for electrons pol6
307  nParamsElect=7;
308  parMomCorrElect[0] = 15.1994;
309  parMomCorrElect[1] = -1.15608;
310  parMomCorrElect[2] = 0.0509546;
311  parMomCorrElect[3] = -0.00122827;
312  parMomCorrElect[4] = 0.0000165499;
313  parMomCorrElect[5] = -0.000000116424;
314  parMomCorrElect[6] = 0.0000333018/100000.;
315 
316  typePar = 14;
317  } else if(run == "aug14_c7") {
318  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
319  nParams = 12;
320  nThetaReg = 22; //from 20 grad till 86 grad
321  thetaStep = 3.;
322  thetaMiddle0 = 21.5;
323  momInt1 = 700;
324  momInt2 = 2500;
325  Double_t par[12*22] = {
326  377.29,-1.86591,0.00389833,-0.00000384956,0.000147982/100000.,108.716,-0.00346533,-0.000906983,7.00846,6.17732,-0.000989527,4.2911,
327  362.489,-1.8547,0.00397136,-0.00000398657,0.000154622/100000.,91.5262,-0.00341486,-0.000743405,5.95775,6.29919,-0.00105339,3.71595,
328  362.489,-1.8547,0.00397136,-0.00000398657,0.000154622/100000.,91.5262,-0.00341486,-0.000743405,5.95775,6.29919,-0.00105339,3.71595,
329  362.489,-1.8547,0.00397136,-0.00000398657,0.000154622/100000.,91.5262,-0.00341486,-0.000743405,5.95775,6.29919,-0.00105339,3.71595,
330  349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976,
331  349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976,
332  349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976,
333  349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976,
334  338.529,-1.81398,0.00403322,-0.00000418336,0.000167129/100000.,76.5941,-0.00346187,-0.000600158,4.79564,5.78833,-0.00112573,3.00538,
335  338.529,-1.81398,0.00403322,-0.00000418336,0.000167129/100000.,76.5941,-0.00346187,-0.000600158,4.79564,5.78833,-0.00112573,3.00538,
336  338.529,-1.81398,0.00403322,-0.00000418336,0.000167129/100000.,76.5941,-0.00346187,-0.000600158,4.79564,5.78833,-0.00112573,3.00538,
337  337.236,-1.79613,0.00397232,-0.00000409852,0.000162844/100000.,77.4287,-0.00345466,-0.000610697,4.88556,5.68106,-0.00111172,3.06302,
338  337.236,-1.79613,0.00397232,-0.00000409852,0.000162844/100000.,77.4287,-0.00345466,-0.000610697,4.88556,5.68106,-0.00111172,3.06302,
339  337.236,-1.79613,0.00397232,-0.00000409852,0.000162844/100000.,77.4287,-0.00345466,-0.000610697,4.88556,5.68106,-0.00111172,3.06302,
340  344.483,-1.82317,0.00401882,-0.00000414132,0.000164604/100000.,79.2906,-0.00340516,-0.000610728,5.07564,6.62067,-0.00114708,3.23332,
341  349.516,-1.83814,0.00402836,-0.00000412621,0.000162918/100000.,82.5493,-0.00342159,-0.00064238,5.28058,6.18647,-0.00111463,3.34096,
342  349.407,-1.80421,0.00387578,-0.00000388649,0.000150131/100000.,87.0329,-0.00345109,-0.000670834,5.47806,6.88032,-0.00113039,3.45412,
343  356.488,-1.84309,0.00397902,-0.00000402071,0.000156785/100000.,88.8678,-0.00343584,-0.000695147,5.65209,6.5324,-0.00110471,3.56262,
344  362.666,-1.85719,0.00397536,-0.00000398646,0.000154446/100000.,93.5664,-0.0034514,-0.000736251,5.93534,6.36717,-0.00107681,3.72505,
345  362.666,-1.85719,0.00397536,-0.00000398646,0.000154446/100000.,93.5664,-0.0034514,-0.000736251,5.93534,6.36717,-0.00107681,3.72505,
346  377.254,-1.89525,0.00397359,-0.00000390074,0.00014792/100000.,95.1167,-0.00338462,-0.000754639,6.17869,7.24617,-0.00110272,3.90527,
347  377.254,-1.89525,0.00397359,-0.00000390074,0.00014792/100000.,95.1167,-0.00338462,-0.000754639,6.17869,7.24617,-0.00110272,3.90527
348  };
349  fillParMomCorrH(nThetaReg*nParams,par);
350 
351  //for electrons pol6
352  nParamsElect=7;
353  parMomCorrElect[0] = 22.9808;
354  parMomCorrElect[1] = -2.53549;
355  parMomCorrElect[2] = 0.139744;
356  parMomCorrElect[3] = -0.00403039;
357  parMomCorrElect[4] = 0.0000628111;
358  parMomCorrElect[5] = -0.000000498781;
359  parMomCorrElect[6] = 0.000158241/100000.;
360 
361  typePar = 14;
362  } else if(run == "apr12") {
363  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
364  nParams = 7;
365  nThetaReg = 22;
366  thetaStep = 3.;
367  thetaMiddle0 = 21.5;
368  Double_t par[7*22] = {
369  -91.8651,0.0434063,678.858,-1120.65,992.81,-617.704,161.391,
370  -98.9038,0.032164,717.974,-1187.02,1071.36,-686.091,185.454,
371  -98.6181,0.0339081,716.166,-1183.66,1066.6,-681.832,184.06,
372  -94.5398,0.0440444,696.528,-1155.27,1041.17,-663.52,178.263,
373  -97.3125,0.0402414,713.659,-1185.34,1075.45,-691.278,187.469,
374  -96.829,0.0388894,711.733,-1182.84,1073.93,-691.003,187.68,
375  -99.0192,0.0378554,724.509,-1204.06,1095.97,-707.519,192.79,
376  -98.4861,0.0402044,722.362,-1200.01,1089.7,-701.493,190.621,
377  -100.162,0.0366257,731.766,-1215.2,1105.28,-713.395,194.44,
378  -95.8465,0.0382859,708.559,-1184.38,1087.36,-706.755,193.663,
379  -69.4872,0.0325539,550.01,-938.753,876.194,-573.516,157.748,
380  -63.2317,0.0276368,511.495,-889.159,856.492,-575.452,161.753,
381  -64.3656,0.0312026,520.12,-904.467,872.46,-587.034,165.186,
382  -68.398,0.0260268,542.253,-938.799,904.784,-610.257,172.378,
383  -69.6065,0.0253658,549.615,-949.945,913.837,-615.737,173.836,
384  -71.149,0.0273896,558.977,-963.62,923.416,-620.118,174.521,
385  -70.434,0.0328006,556.942,-959.752,915.907,-611.974,171.355,
386  -73.0875,0.0344005,573.283,-985.193,938.032,-626.329,175.396,
387  -73.3104,0.0371747,575.988,-988.623,937.965,-624.185,174.31,
388  -74.2248,0.0409931,582.787,-998.393,943.378,-625.427,174.068,
389  -78.2616,0.0396539,605.01,-1030.55,967.985,-639.483,177.565,
390  -83.4854,0.0227933,629.292,-1062.32,989.642,-651.417,180.689 };
391  fillParMomCorrH(nThetaReg*nParams,par);
392  typePar = 12;
393  } else if(run == "nov02") {
394  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
395  momInt1=600.;
396  momInt2=2000.;
397  nThetaReg=7;
398  thetaStep=10.;
399  thetaMiddle0=23.;
400  nParams=10;
401  Double_t par[70] = {
402  0.0506414,-68.3320,30017.5,87.1936,-0.00405782,-0.000841998,4.52039,5.37248,-0.00121608,2.40080, // 1st theta interval
403  0.056843, -73.7507,30890.6,85.2889,-0.00406022,-0.000823641,4.41504,3.97163,-0.00105256,2.31750, // 2nd theta interval
404  0.0692018,-85.0326,33230.6,86.0512,-0.00408880,-0.000838878,4.39438,4.30892,-0.00110695,2.29046, // 3rd theta interval
405  0.05393, -72.6496,31345.2,92.3181,-0.00414646,-0.000836083,4.41321,4.73116,-0.00115162,2.30843, // 4th theta interval
406  0.052317, -70.9782,31023.3,94.4471,-0.00416224,-0.000882047,4.56496,5.04969,-0.00117343,2.36712, // 5th theta interval
407  0.0524236,-71.4746,31505.2,97.8558,-0.00416946,-0.000919853,4.75295,5.53359,-0.00120439,2.46908, // 6th theta interval
408  0.0532758,-72.6315,32132.7,100.503,-0.00416881,-0.000962655,4.95385,5.78055,-0.00120932,2.57004 }; // 7th theta interval
409  fillParMomCorrH(nThetaReg*nParams,par);
410 
411  //for electrons
412  nParamsElect=5;
413  parMomCorrElect[0]=2.29512;
414  parMomCorrElect[1]=0.042546;
415  parMomCorrElect[2]=-0.00175886;
416  parMomCorrElect[3]=2.47346/100000;
417  parMomCorrElect[4]=-1.0764/10000000;
418 
419  typePar = 1;
420  } else if(run=="jan04" || run=="apr07") {
421  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
422  momInt1=720.;
423  momInt2=2000.;
424  nThetaReg=7;
425  nParams=10;
426  thetaStep=10.;
427  thetaMiddle0=23.;
428 
429  Double_t par[70] = {
430  0.0404892,-63.2212,31512.3,88.7387,-0.00383467,-0.000831376,4.83807,8.32566,-0.00134054,2.69398, // 1st theta interval
431  0.0392893,-61.7120,31187.0,67.0670,-0.00344743,-0.000770298,4.73321,4.08464,-0.00100119,2.70498, // 2nd theta interval
432  0.0402474,-63.1836,31819.4,95.3604,-0.00391875,-0.000930825,5.09812,4.57863,-0.00103624,2.68722, // 3rd theta interval
433  0.0421779,-64.9599,31849.2,85.6957,-0.00379927,-0.000812833,4.78868,9.17729,-0.00138972,2.64987, // 4th theta interval
434  0.0417347,-64.2374,31548.9,97.1928,-0.00399476,-0.000827159,4.80468,6.59782,-0.00125117,2.63429, // 4th theta interval
435  0.0404892,-63.2212,31512.3,88.7387,-0.00383467,-0.000831377,4.83807,8.32566,-0.00134054,2.69398, // 6th theta interval
436  0.0396831,-62.4906,31610.8,86.3941,-0.00374661,-0.000816029,4.88749,4.51922,-0.00104096,2.73135}; // 7th theta interval
437  fillParMomCorrH(nThetaReg*nParams,par);
438 
439  //for electrons
440  nParamsElect=5;
441  parMomCorrElect[0]=2.21734;
442  parMomCorrElect[1]=0.0794622;
443  parMomCorrElect[2]=-0.00258538;
444  parMomCorrElect[3]=3.03896/100000;
445  parMomCorrElect[4]=-1.11898/10000000;
446 
447  typePar = 1;
448  } else if(run=="aug04") {
449  Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
450  momInt1=800.;
451  momInt2=2000.;
452  nThetaReg=7;
453  thetaStep=10.;
454  thetaMiddle0=23.;
455  nParams=10;
456  Double_t par[70] = {
457  0.0319386,-52.3932,28573.9,65.0824,-0.00340829,-0.000751249,4.91880,7.01293,-0.001202890,2.87892, // 1st theta interval
458  0.0354272,-57.1894,30552.1,73.0873,-0.00349743,-0.000814299,5.21396,5.66944,-0.001081390,2.99888, // 2nd theta interval
459  0.0369929,-59.2908,31560.5,74.1808,-0.00343821,-0.000746043,5.24870,5.80041,-0.001114980,3.16789, // 3rd theta interval
460  0.0381148,-62.3502,34058.2,74.2051,-0.00327892,-0.000714812,5.57025,7.58451,-0.001158990,3.49298, // 4th theta interval
461  0.0405264,-66.2530,36385.7,80.9644,-0.00328172,-0.000733450,6.02206,5.72358,-0.000985652,3.77823, // 5th theta interval
462  0.0442906,-72.6287,39923.8,112.382,-0.00364194,-0.001174410,7.26952,10.0128,-0.001196320,4.12544, // 6th theta interval
463  0.0419634,-70.9195,40627.6,122.213,-0.00366415,-0.001261510,7.75827,9.35643,-0.001134410,4.36657 }; // 7th theta interval
464  fillParMomCorrH(nThetaReg*nParams,par);
465 
466  //for electrons
467  nParamsElect=5;
468  parMomCorrElect[0]=3.21963;
469  parMomCorrElect[1]=-0.0257486;
470  parMomCorrElect[2]=0.000407469;
471  parMomCorrElect[3]=9.67088e-06;
472  parMomCorrElect[4]=-9.93232e-08;
473  typePar = 1;
474  } else{
475  Warning("setDefaultPar()","unknown run %s, no correction table loaded!",run.Data());
476  }
477  return typePar > 0;
478 }
479 void HEnergyLossCorrPar::fillParMomCorrH(Int_t size,Double_t *par) {
480  for(Int_t i=0;i<size;i++) parMomCorrH[i] = par[i];
481 }
482 
483 Double_t HEnergyLossCorrPar::getDeltaMom(Int_t pId, Double_t mom, Double_t theta) const {
484  mom = TMath::Abs(mom);
485  if(typePar == 1) return getDeltaMomT1(pId,mom,theta);
486  if(typePar == 12) return getDeltaMomT12(pId,mom,theta);
487  if(typePar == 14 || typePar == 19) return getDeltaMomT1419(pId,mom,theta);
488  else return 0.;
489 }
490 
491 
492 Double_t HEnergyLossCorrPar::getDeltaMomT1(Int_t pId, Double_t mom, Double_t theta) const {
493 
494 // if(theta<15.||theta>87.) return 0.;
495  if(theta<15.) theta = 15.;
496  if(theta>87.) theta = 87.;
497 
498  if(pId==2||pId==3) {
499 // if(mom<50.) return 0.;
500  if(mom<40.) mom = 40.;
501  // 4th order polynomial
502  return parMomCorrElect[0]+theta*(parMomCorrElect[1]+theta*(parMomCorrElect[2]+theta*(parMomCorrElect[3]+parMomCorrElect[4]*theta)));
503  }
504  else {
505 
506  Double_t pMom, K[nParams];
507  Double_t mPr = HPhysicsConstants::mass(14);
508  Double_t mPart = HPhysicsConstants::mass(pId);
509  Int_t chg = HPhysicsConstants::charge(pId);
510  if(mPart<=0.||(chg==0)) return 0.;
511 
512  if(pId==14) pMom=mom;
513  else pMom = Sqrt(Power(mPr/mPart*(Sqrt(mom*mom+mPart*mPart)-mPart)+mPr,2)-Power(mPr,2));
514 
515 // if(pMom<170.) return 0.;
516  if(pMom<170.) pMom = 170.;
517 
518  Int_t i=Int_t((theta-thetaMiddle0)/thetaStep+1);
519  if(i==nThetaReg) i=nThetaReg-1;
520 
521  Int_t j;
522 
523  if(pMom<momInt1) j=0; else if(pMom<momInt2) j=3; else j=7;
524  if(i==0) {
525  // first theta interval
526  Double_t thVal = (thetaMiddle0-theta)/thetaStep;
527  K[j] = parMomCorrH[j] - (parMomCorrH[nParams+j]-parMomCorrH[j])*thVal;
528  K[j+1] = parMomCorrH[j+1] - (parMomCorrH[nParams+j+1]-parMomCorrH[j+1])*thVal;
529  K[j+2] = parMomCorrH[j+2] - (parMomCorrH[nParams+j+2]-parMomCorrH[j+2])*thVal;
530  if(pMom>momInt1 && pMom<momInt2)
531  K[j+3] = parMomCorrH[j+3] - (parMomCorrH[nParams+j+3]-parMomCorrH[j+3])*thVal;
532 
533 // K[j]=parMomCorrH[j]; K[j+1]=parMomCorrH[j+1]; K[j+2]=parMomCorrH[j+2];
534 // if(pMom>momInt1&&pMom<momInt2) K[j+3]=parMomCorrH[j+3];
535  }//i=0
536  else {
537  Int_t bound1 = (i-1)*nParams;
538  Int_t bound2 = i*nParams;
539  Double_t thInsideInt = (theta-(thetaMiddle0+thetaStep*(i-1)))/thetaStep;
540  if(i==nThetaReg-1) {
541  // last theta interval
542  K[j]=parMomCorrH[bound2+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
543  K[j+1]=parMomCorrH[bound2+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
544  K[j+2]=parMomCorrH[bound2+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
545  if(pMom>momInt1&&pMom<momInt2)
546  K[j+3]=parMomCorrH[bound2+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
547  }
548  else {
549  // theta intervals excluding first and last
550  K[j]=parMomCorrH[bound1+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
551  K[j+1]=parMomCorrH[bound1+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
552  K[j+2]=parMomCorrH[bound1+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
553  if(pMom>momInt1&&pMom<momInt2)
554  K[j+3]=parMomCorrH[bound1+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
555  } //else
556  } //else
557 
558  if(pMom<momInt1) return (K[j]*pMom+K[j+1]+K[j+2]/pMom)*chg*chg;
559  if(pMom<momInt2) return (K[j]*exp(K[j+1]*pMom)+K[j+2]*pMom+K[j+3])*chg*chg;
560  return (K[j]*exp(K[j+1]*pMom)+K[j+2])*chg*chg;
561  } // not leptons
562 
563  return 0.;
564 }
565 
566 
567 Double_t HEnergyLossCorrPar::binInter(Double_t b,Double_t rs, const Double_t *par) const {
568  if(b>=1. || b<=0.) return 0.;
569 
570  Double_t b2 = b*b;
571  Double_t bg2 = b/Sqrt(1.-b2);
572  Double_t lnbg2 = Log(bg2)/b2;
573  Double_t odb = 1./b;
574  Double_t osb = Sqrt(odb);
575  Double_t dp1 = par[0]*odb+par[1]*lnbg2+par[2]*osb + par[3]+par[4]*b+par[5]*b2+par[6]*b2*b;
576  par += nParams; // next theta region
577  Double_t dp2 = par[0]*odb+par[1]*lnbg2+par[2]*osb + par[3]+par[4]*b+par[5]*b2+par[6]*b2*b;
578  return dp1 + (dp2-dp1)*rs;
579 }
580 
581 Double_t HEnergyLossCorrPar::getDeltaMomT12(Int_t pId, Double_t mom, Double_t theta) const {
582 
583 // if(theta<15.||theta>87.) return 0.;
584  if(theta<15.) theta = 15.;
585  if(theta>87.) theta = 87.;
586 
587  if(pId==2||pId==3) {
588 // if(mom<10.) return 0.;
589  if(mom<40.) mom = 40.;
590 
591  if(theta<28.) return 2.96406766784;
592  if(theta<50.88) return -7.83473+1.31585*theta-0.0591126*Power(theta,2)+0.00115799*Power(theta,3)-0.00000833136*Power(theta,4);
593  if(theta<51.81) return 27.1723-0.479349*theta;
594  if(theta<80.28) return 2916.99-261.345*theta+9.7235*Power(theta,2)-0.192166*Power(theta,3)+0.00212771*Power(theta,4)
595  -0.0000125143*Power(theta,5)+0.0000000305468*Power(theta,6);
596  return 0.505533+0.0260891*theta;
597 
598  } // leptons
599  else {
600  // beamtime = apr12
601  Double_t mPart = HPhysicsConstants::mass(pId);
602  Int_t chg2 = HPhysicsConstants::charge(pId);
603  if(mPart<=0.|| chg2==0) return 0.;
604  chg2 = chg2*chg2;
605  Double_t b = mom/Sqrt(mom*mom + mPart*mPart); // Beta
606  if(b < 0.11) return 0; // b=0.11 => Pproton = 100 MeV/c
607  Int_t i = Int_t((theta-thetaMiddle0)/thetaStep);
608  if(i < 0) i = 0;
609  if(i >= nThetaReg-1) i = nThetaReg-2;
610  Double_t mixtCoef = (theta-(thetaMiddle0+thetaStep*i))/thetaStep;
611  return chg2*binInter(b,mixtCoef, parMomCorrH+i*nParams);
612 
613  } // not leptons
614 
615  return 0.;
616 }
617 
618 
619 Double_t HEnergyLossCorrPar::getDeltaMomT1419(Int_t pId, Double_t mom, Double_t theta) const {
620 
621 // beamtime = jul14 and aug14 and mar19
622 // if(theta<15.||theta>86.) return 0.;
623  if(theta<15.) theta = 15.;
624  if(theta>87.) theta = 87.;
625 
626  if(pId==2||pId==3) {
627  if(mom<40.) mom = 40.; //return 0.;
628  if(typePar == 14)
629  // pol6
630  return parMomCorrElect[0]+theta*(parMomCorrElect[1]+theta*(parMomCorrElect[2]+theta*(parMomCorrElect[3]+theta*(parMomCorrElect[4]+theta*(parMomCorrElect[5]+parMomCorrElect[6]*theta)))));
631 
632  if(typePar == 19) {
633  // pol4 pol2 pol4
634  Double_t smCor = 0.3;
635 
636  if(theta<23.22) return 2.35+smCor;
637  if(theta<50.4) return parMomCorrElect[0]+parMomCorrElect[1]*theta+parMomCorrElect[2]*Power(theta,2)+parMomCorrElect[3]*Power(theta,3)+parMomCorrElect[4]*Power(theta,4)+smCor;
638  if(theta<53.5) return parMomCorrElect[5]+parMomCorrElect[6]*theta+parMomCorrElect[7]*Power(theta,2)+smCor;
639  if(theta<83.3) return parMomCorrElect[8]+parMomCorrElect[9]*theta+parMomCorrElect[10]*Power(theta,2)+parMomCorrElect[11]*Power(theta,3)+parMomCorrElect[12]*Power(theta,4)+smCor;
640  return 1.74+smCor;
641 
642 
643  } // mar19
644  } // leptons
645  else {
646 
647  Double_t pMom, K[nParams];
648  Double_t mPr = HPhysicsConstants::mass(14);
649  Double_t mPart = HPhysicsConstants::mass(pId);
650  Int_t chg = HPhysicsConstants::charge(pId);
651  if(mPart<=0.||(chg==0)) return 0.;
652 
653  if(pId==14) pMom=mom;
654  else pMom = Sqrt(Power(mPr/mPart*(Sqrt(mom*mom+mPart*mPart)-mPart)+mPr,2)-Power(mPr,2));
655 
656  if(typePar == 14) if(pMom<170.) pMom = 170.;
657  if(typePar == 19) if(pMom<100.) pMom = 100.;
658 
659  Int_t i=Int_t((theta-thetaMiddle0)/thetaStep+1);
660  if(i==nThetaReg) i=nThetaReg-1;
661 
662  Int_t j;
663 
664  if(pMom<momInt1) j=0; else if(pMom<momInt2) j=5; else j=9;
665  if(i==0) {
666  // first theta interval
667  Double_t thVal = (thetaMiddle0-theta)/thetaStep;
668  K[j] = parMomCorrH[j] - (parMomCorrH[nParams+j]-parMomCorrH[j])*thVal;
669  K[j+1] = parMomCorrH[j+1] - (parMomCorrH[nParams+j+1]-parMomCorrH[j+1])*thVal;
670  K[j+2] = parMomCorrH[j+2] - (parMomCorrH[nParams+j+2]-parMomCorrH[j+2])*thVal;
671  if(pMom<momInt2)
672  K[j+3] = parMomCorrH[j+3] - (parMomCorrH[nParams+j+3]-parMomCorrH[j+3])*thVal;
673  if(pMom<momInt1)
674  K[j+4] = parMomCorrH[j+4] - (parMomCorrH[nParams+j+4]-parMomCorrH[j+4])*thVal;
675 
676  }//i=0
677  else {
678  Int_t bound1 = (i-1)*nParams;
679  Int_t bound2 = i*nParams;
680  Double_t thInsideInt = (theta-(thetaMiddle0+thetaStep*(i-1)))/thetaStep;
681  if(i==nThetaReg-1) {
682  // last theta interval
683  K[j]=parMomCorrH[bound2+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
684  K[j+1]=parMomCorrH[bound2+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
685  K[j+2]=parMomCorrH[bound2+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
686  if(pMom<momInt2)
687  K[j+3]=parMomCorrH[bound2+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
688  if(pMom<momInt1)
689  K[j+4]=parMomCorrH[bound2+j+4] + (parMomCorrH[bound2+j+4]-parMomCorrH[bound1+j+4])*thInsideInt;
690  }
691  else {
692  // theta intervals excluding first and last
693  K[j]=parMomCorrH[bound1+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
694  K[j+1]=parMomCorrH[bound1+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
695  K[j+2]=parMomCorrH[bound1+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
696  if(pMom<momInt2)
697  K[j+3]=parMomCorrH[bound1+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
698  if(pMom<momInt1)
699  K[j+4]=parMomCorrH[bound1+j+4] + (parMomCorrH[bound2+j+4]-parMomCorrH[bound1+j+4])*thInsideInt;
700  } //else
701  } // i!=0
702 
703  if(pMom<momInt1) return (K[j]+pMom*(K[j+1]+pMom*(K[j+2]+pMom*(K[j+3]+K[j+4]*pMom))))*chg*chg;
704  if(pMom<momInt2) return (K[j]*exp(K[j+1]*pMom)+K[j+2]*pMom+K[j+3])*chg*chg;
705  return (K[j]*exp(K[j+1]*pMom)+K[j+2])*chg*chg;
706  } // not leptons
707 
708  return 0.;
709 }
710 
711 
HEnergyLossCorrPar::getDeltaMomT1
Double_t getDeltaMomT1(Int_t pId, Double_t mom, Double_t theta) const
Definition: henergylosscorrpar.cc:492
HPhysicsConstants::charge
static Int_t charge(const Int_t id)
Definition: hphysicsconstants.cc:491
HEnergyLossCorrPar::binInter
Double_t binInter(Double_t b, Double_t rs, const Double_t *par) const
Definition: henergylosscorrpar.cc:567
HEnergyLossCorrPar::setDefaultPar
Bool_t setDefaultPar(TString run)
Definition: henergylosscorrpar.cc:121
HEnergyLossCorrPar::getParams
Bool_t getParams(HParamList *)
Definition: henergylosscorrpar.cc:105
HEnergyLossCorrPar::getDeltaMomT1419
Double_t getDeltaMomT1419(Int_t pId, Double_t mom, Double_t theta) const
Definition: henergylosscorrpar.cc:619
theta
Double_t theta
Definition: drawResolution.C:101
HParamList::add
void add(HParamObj &)
Definition: hparamlist.cc:415
HPhysicsConstants::mass
static Float_t mass(const Int_t id)
Definition: hphysicsconstants.cc:558
ClassImp
ClassImp(HEnergyLossCorrPar) HEnergyLossCorrPar
Definition: henergylosscorrpar.cc:62
HEnergyLossCorrPar::gEnergyLossCorrPar
static HEnergyLossCorrPar * gEnergyLossCorrPar
Definition: henergylosscorrpar.h:22
HEnergyLossCorrPar::getDeltaMom
Double_t getDeltaMom(Int_t pId, Double_t mom, Double_t theta) const
Definition: henergylosscorrpar.cc:483
HEnergyLossCorrPar::getDeltaMomT12
Double_t getDeltaMomT12(Int_t pId, Double_t mom, Double_t theta) const
Definition: henergylosscorrpar.cc:581
HParamList::fill
Bool_t fill(const Text_t *, Text_t *, const Int_t)
Definition: hparamlist.cc:561
HEnergyLossCorrPar::fillParMomCorrH
void fillParMomCorrH(Int_t size, Double_t *par)
Definition: henergylosscorrpar.cc:479
HEnergyLossCorrPar::putParams
void putParams(HParamList *)
Definition: henergylosscorrpar.cc:89