ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Utils/src/PhotonFix.cc
Revision: 1.2
Committed: Thu Sep 8 15:51:24 2011 UTC (13 years, 7 months ago) by bendavid
Content type: text/plain
Branch: MAIN
CVS Tags: Mit_025d, Mit_025c, Mit_025b, Mit_025a, Mit_025, Mit_025pre2, Mit_025pre1
Changes since 1.1: +2147 -1672 lines
Log Message: 
Add tool and required changes for photon energy regression

File Contents

# User Rev Content
1 bendavid 1.2 #include <cmath>
2     #include <cassert>
3     #include <fstream>
4     #include <iomanip>
5    
6     // ensure that this include points to the appropriate location for PhotonFix.h
7     #include "../interface/PhotonFix.h"
8    
9     /*PhotonFix::PhotonFix(double e, double eta, double phi, double r9) :
10     _e(e), _eta(eta), _phi(phi), _r9(r9) {
11    
12     setup();
13     }*/
14    
15     void PhotonFix::setup(double e, double eta, double phi, double r9){
16     // Check constants have been set up
17     assert(_initialised);
18    
19     _e = e;
20     _eta = eta;
21     _phi = phi;
22     _r9 = r9;
23    
24     // Determine if EB or EE
25     _be=(fabs(_eta)<1.48?0:1);
26    
27     // Determine if high or low R9
28     if(_be==0) _hl=(_r9>=0.94?0:1);
29     else _hl=(_r9>=0.95?0:1);
30    
31     // Coordinates relative to cracks
32     double r2Min;
33     if(_be==0) {
34    
35     r2Min=1.0e6;
36     for(unsigned i(0);i<169;i++) {
37     for(unsigned j(0);j<360;j++) {
38     double de(_eta-_barrelCGap[i][j][0]);
39     double df(dPhi(_phi,_barrelCGap[i][j][1]));
40     double r2(de*de+df*df);
41    
42     if(r2<r2Min) {
43     r2Min=r2;
44     if(i>=84) {
45     _aC= de;
46     _bC=-df;
47     } else {
48     _aC=-de;
49     _bC= df;
50     }
51     }
52     }
53     }
54    
55     r2Min=1.0e6;
56     for(unsigned i(0);i<33;i++) {
57     for(unsigned j(0);j<180;j++) {
58     double de(_eta-_barrelSGap[i][j][0]);
59     double df(dPhi(_phi,_barrelSGap[i][j][1]));
60     double r2(de*de+df*df);
61    
62     if(r2<r2Min) {
63     r2Min=r2;
64     if(i>=16) {
65     _aS= de;
66     _bS=-df;
67     } else {
68     _aS=-de;
69     _bS= df;
70     }
71     }
72     }
73     }
74    
75     r2Min=1.0e6;
76     for(unsigned i(0);i<7;i++) {
77     for(unsigned j(0);j<18;j++) {
78     double de(_eta-_barrelMGap[i][j][0]);
79     double df(dPhi(_phi,_barrelMGap[i][j][1]));
80     double r2(de*de+df*df);
81    
82     if(r2<r2Min) {
83     r2Min=r2;
84     if(i>=3) {
85     _aM= de;
86     _bM=-df;
87     } else {
88     _aM=-de;
89     _bM= df;
90     }
91     }
92     }
93     }
94    
95     } else {
96     unsigned iz(_eta>=0.0?0:1);
97     double r[2]={xZ(),yZ()};
98    
99     r2Min=1.0e6;
100     for(unsigned i(0);i<7080;i++) {
101     double dx(r[0]-_endcapCGap[iz][i][0]);
102     double dy(r[1]-_endcapCGap[iz][i][1]);
103     double r2(dx*dx+dy*dy);
104    
105     if(r2<r2Min) {
106     r2Min=r2;
107     if(r[0]>0.0) _aC= dx;
108     else _aC=-dx;
109     if(r[1]>0.0) _bC= dy;
110     else _bC=-dy;
111     }
112     }
113    
114     r2Min=1.0e6;
115     for(unsigned i(0);i<264;i++) {
116     double dx(r[0]-_endcapSGap[iz][i][0]);
117     double dy(r[1]-_endcapSGap[iz][i][1]);
118     double r2(dx*dx+dy*dy);
119    
120     if(r2<r2Min) {
121     r2Min=r2;
122     if(r[0]>0.0) _aS= dx;
123     else _aS=-dx;
124     if(r[1]>0.0) _bS= dy;
125     else _bS=-dy;
126     }
127     }
128    
129     r2Min=1.0e6;
130     for(unsigned i(0);i<1;i++) {
131     double dx(r[0]-_endcapMGap[iz][i][0]);
132     double dy(r[1]-_endcapMGap[iz][i][1]);
133     double r2(dx*dx+dy*dy);
134    
135     if(r2<r2Min) {
136     r2Min=r2;
137     if(iz==0) {_aM= dx;_bM= dy;}
138     else {_aM=-dx;_bM=-dy;}
139     }
140     }
141     }
142     }
143    
144     double PhotonFix::fixedEnergy() const {
145     double f(0.0);
146    
147     // Overall scale and energy(T) dependence
148     f =_meanScale[_be][_hl][0];
149     f+=_meanScale[_be][_hl][1]*_e;
150     f+=_meanScale[_be][_hl][2]*_e/cosh(_eta);
151     f+=_meanScale[_be][_hl][3]*cosh(_eta)/_e;
152    
153     // General eta or zeta dependence
154     if(_be==0) {
155     f+=_meanAT[_be][_hl][0]*_eta*_eta;
156     f+=expCorrection(_eta,_meanBT[_be][_hl]);
157     } else {
158     f+=_meanAT[_be][_hl][0]*xZ()*xZ();
159     f+=_meanBT[_be][_hl][0]*yZ()*yZ();
160     }
161    
162     // Eta or x crystal, submodule and module dependence
163     f+=expCorrection(_aC,_meanAC[_be][_hl]);
164     f+=expCorrection(_aS,_meanAS[_be][_hl]);
165     f+=expCorrection(_aM,_meanAM[_be][_hl]);
166    
167     // Phi or y crystal, submodule and module dependence
168     f+=expCorrection(_bC,_meanBC[_be][_hl]);
169     f+=expCorrection(_bS,_meanBS[_be][_hl]);
170     f+=expCorrection(_bM,_meanBM[_be][_hl]);
171    
172     // R9 dependence
173     if(_hl==0) {
174     f+=_meanR9[_be][_hl][1]*(_r9-_meanR9[_be][_hl][0])*(_r9-_meanR9[_be][_hl][0])
175     +_meanR9[_be][_hl][2]*(_r9-_meanR9[_be][_hl][0])*(_r9-_meanR9[_be][_hl][0])*(_r9-_meanR9[_be][_hl][0]);
176     } else {
177     f+=_meanR9[_be][_hl][0]*_r9+_meanR9[_be][_hl][1]*_r9*_r9+_meanR9[_be][_hl][2]*_r9*_r9*_r9;
178     }
179    
180     return _e*f;
181     }
182    
183     double PhotonFix::sigmaEnergy() const {
184    
185     // Overall resolution scale vs energy
186     double sigma;
187     if(_be==0) {
188     sigma =_sigmaScale[_be][_hl][0]*_sigmaScale[_be][_hl][0];
189     //std::cout << "PhotonFix::sigmaEnergy 1 sigma = " << sigma << std::endl;
190     sigma+=_sigmaScale[_be][_hl][1]*_sigmaScale[_be][_hl][1]*_e;
191     //std::cout << "PhotonFix::sigmaEnergy 2 sigma = " << sigma << std::endl;
192     sigma+=_sigmaScale[_be][_hl][2]*_sigmaScale[_be][_hl][2]*_e*_e;
193     //std::cout << "PhotonFix::sigmaEnergy 3 sigma = " << sigma << std::endl;
194     } else {
195     sigma =_sigmaScale[_be][_hl][0]*_sigmaScale[_be][_hl][0]*cosh(_eta)*cosh(_eta);
196     sigma+=_sigmaScale[_be][_hl][1]*_sigmaScale[_be][_hl][1]*_e;
197     sigma+=_sigmaScale[_be][_hl][2]*_sigmaScale[_be][_hl][2]*_e*_e;
198     }
199     sigma=sqrt(sigma);
200    
201     double f(1.0);
202    
203     // General eta or zeta dependence
204     if(_be==0) {
205     f+=_sigmaAT[_be][_hl][0]*_eta*_eta;
206     //std::cout << "PhotonFix::sigmaEnergy 4 f = " << f << std::endl;
207     f+=expCorrection(_eta,_sigmaBT[_be][_hl]);
208     //std::cout << "PhotonFix::sigmaEnergy 5 f = " << f << std::endl;
209     } else {
210     f+=_sigmaAT[_be][_hl][0]*xZ()*xZ();
211     f+=_sigmaBT[_be][_hl][0]*yZ()*yZ();
212     }
213    
214     // Eta or x crystal, submodule and module dependence
215     f+=expCorrection(_aC,_sigmaAC[_be][_hl]);
216     //std::cout << "PhotonFix::sigmaEnergy 6 f = " << f << std::endl;
217     f+=expCorrection(_aS,_sigmaAS[_be][_hl]);
218     //std::cout << "PhotonFix::sigmaEnergy 7 f = " << f << std::endl;
219     f+=expCorrection(_aM,_sigmaAM[_be][_hl]);
220     //std::cout << "PhotonFix::sigmaEnergy 8 f = " << f << std::endl;
221    
222     // Phi or y crystal, submodule and module dependence
223     f+=expCorrection(_bC,_sigmaBC[_be][_hl]);
224     //std::cout << "PhotonFix::sigmaEnergy 9 f = " << f << std::endl;
225     f+=expCorrection(_bS,_sigmaBS[_be][_hl]);
226     //std::cout << "PhotonFix::sigmaEnergy 10 f = " << f << std::endl;
227     f+=expCorrection(_bM,_sigmaBM[_be][_hl]);
228     //std::cout << "PhotonFix::sigmaEnergy 11 f = " << f << std::endl;
229    
230     // R9 dependence
231     if(_hl==0) {
232     f+=_sigmaR9[_be][_hl][1]*(_r9-_sigmaR9[_be][_hl][0])*(_r9-_sigmaR9[_be][_hl][0])
233     +_sigmaR9[_be][_hl][2]*(_r9-_sigmaR9[_be][_hl][0])*(_r9-_sigmaR9[_be][_hl][0])*(_r9-_sigmaR9[_be][_hl][0]);
234     //std::cout << "PhotonFix::sigmaEnergy 12 f = " << f << std::endl;
235     } else {
236     f+=_sigmaR9[_be][_hl][0]*_r9+_sigmaR9[_be][_hl][1]*_r9*_r9+_sigmaR9[_be][_hl][2]*_r9*_r9*_r9;
237     //std::cout << "PhotonFix::sigmaEnergy 13 f = " << f << std::endl;
238     }
239    
240     return sigma*f;
241     }
242    
243     double PhotonFix::rawEnergy() const {
244     return _e;
245     }
246    
247     double PhotonFix::eta() const {
248     return _eta;
249     }
250    
251     double PhotonFix::phi() const {
252     return _phi;
253     }
254    
255     double PhotonFix::r9() const {
256     return _r9;
257     }
258    
259     double PhotonFix::etaC() const {
260     assert(_be==0);
261     return _aC;
262     }
263    
264     double PhotonFix::etaS() const {
265     assert(_be==0);
266     return _aS;
267     }
268    
269     double PhotonFix::etaM() const {
270     assert(_be==0);
271     return _aM;
272     }
273    
274     double PhotonFix::phiC() const {
275     assert(_be==0);
276     return _bC;
277     }
278    
279     double PhotonFix::phiS() const {
280     assert(_be==0);
281     return _bS;
282     }
283    
284     double PhotonFix::phiM() const {
285     assert(_be==0);
286     return _bM;
287     }
288    
289     double PhotonFix::xZ() const {
290     assert(_be==1);
291     return asinh(cos(_phi)/sinh(_eta));
292     }
293    
294     double PhotonFix::xC() const {
295     assert(_be==1);
296     return _aC;
297     }
298    
299     double PhotonFix::xS() const {
300     assert(_be==1);
301     return _aS;
302     }
303    
304     double PhotonFix::xM() const {
305     assert(_be==1);
306     return _aM;
307     }
308    
309     double PhotonFix::yZ() const {
310     assert(_be==1);
311     return asinh(sin(_phi)/sinh(_eta));
312     }
313    
314     double PhotonFix::yC() const {
315     assert(_be==1);
316     return _bC;
317     }
318    
319     double PhotonFix::yS() const {
320     assert(_be==1);
321     return _bS;
322     }
323    
324     double PhotonFix::yM() const {
325     assert(_be==1);
326     return _bM;
327     }
328    
329     double PhotonFix::GetaPhi(double f0, double f1) const {
330     return aPhi(f0,f1);
331     }
332    
333     void PhotonFix::barrelCGap(unsigned i, unsigned j, unsigned k, double c){
334     _barrelCGap[i][j][k] = c;
335     }
336     void PhotonFix::barrelSGap(unsigned i, unsigned j, unsigned k, double c){
337     _barrelSGap[i][j][k] = c;
338     }
339     void PhotonFix::barrelMGap(unsigned i, unsigned j, unsigned k, double c){
340     _barrelMGap[i][j][k] = c;
341     }
342     void PhotonFix::endcapCrystal(unsigned i, unsigned j, bool c){
343     _endcapCrystal[i][j] = c;
344     }
345     void PhotonFix::endcapCGap(unsigned i, unsigned j, unsigned k, double c){
346     _endcapCGap[i][j][k] = c;
347     }
348     void PhotonFix::endcapSGap(unsigned i, unsigned j, unsigned k, double c){
349     _endcapSGap[i][j][k] = c;
350     }
351     void PhotonFix::endcapMGap(unsigned i, unsigned j, unsigned k, double c){
352     _endcapMGap[i][j][k] = c;
353     }
354    
355    
356     void PhotonFix::print() const {
357     std::cout << "PhotonFix: e,eta,phi,r9 = " << _e << ", " << _eta << ", " << _phi << ", " << _r9 << ", gaps "
358     << _aC << ", " << _aS << ", " << _aM << ", "
359     << _bC << ", " << _bS << ", " << _bM << std::endl;
360     }
361    
362     void PhotonFix::setParameters(unsigned be, unsigned hl, const double *p) {
363     for(unsigned i(0);i<4;i++) {
364     _meanScale[be][hl][i] =p[i+ 0*4];
365     _meanAT[be][hl][i] =p[i+ 1*4];
366     _meanAC[be][hl][i] =p[i+ 2*4];
367     _meanAS[be][hl][i] =p[i+ 3*4];
368     _meanAM[be][hl][i] =p[i+ 4*4];
369     _meanBT[be][hl][i] =p[i+ 5*4];
370     _meanBC[be][hl][i] =p[i+ 6*4];
371     _meanBS[be][hl][i] =p[i+ 7*4];
372     _meanBM[be][hl][i] =p[i+ 8*4];
373     _meanR9[be][hl][i] =p[i+ 9*4];
374    
375     _sigmaScale[be][hl][i]=p[i+10*4];
376     _sigmaAT[be][hl][i] =p[i+11*4];
377     _sigmaAC[be][hl][i] =p[i+12*4];
378     _sigmaAS[be][hl][i] =p[i+13*4];
379     _sigmaAM[be][hl][i] =p[i+14*4];
380     _sigmaBT[be][hl][i] =p[i+15*4];
381     _sigmaBC[be][hl][i] =p[i+16*4];
382     _sigmaBS[be][hl][i] =p[i+17*4];
383     _sigmaBM[be][hl][i] =p[i+18*4];
384     _sigmaR9[be][hl][i] =p[i+19*4];
385     }
386     }
387    
388     void PhotonFix::getParameters(unsigned be, unsigned hl, double *p) {
389     for(unsigned i(0);i<4;i++) {
390     p[i+ 0*4]=_meanScale[be][hl][i];
391     p[i+ 1*4]=_meanAT[be][hl][i];
392     p[i+ 2*4]=_meanAC[be][hl][i];
393     p[i+ 3*4]=_meanAS[be][hl][i];
394     p[i+ 4*4]=_meanAM[be][hl][i];
395     p[i+ 5*4]=_meanBT[be][hl][i];
396     p[i+ 6*4]=_meanBC[be][hl][i];
397     p[i+ 7*4]=_meanBS[be][hl][i];
398     p[i+ 8*4]=_meanBM[be][hl][i];
399     p[i+ 9*4]=_meanR9[be][hl][i];
400    
401     p[i+10*4]=_sigmaScale[be][hl][i];
402     p[i+11*4]=_sigmaAT[be][hl][i];
403     p[i+12*4]=_sigmaAC[be][hl][i];
404     p[i+13*4]=_sigmaAS[be][hl][i];
405     p[i+14*4]=_sigmaAM[be][hl][i];
406     p[i+15*4]=_sigmaBT[be][hl][i];
407     p[i+16*4]=_sigmaBC[be][hl][i];
408     p[i+17*4]=_sigmaBS[be][hl][i];
409     p[i+18*4]=_sigmaBM[be][hl][i];
410     p[i+19*4]=_sigmaR9[be][hl][i];
411     }
412     }
413    
414     void PhotonFix::dumpParameters(std::ostream &o) {
415     o << std::setprecision(9);
416    
417     for(unsigned be(0);be<2;be++) {
418     for(unsigned hl(0);hl<2;hl++) {
419     for(unsigned i(0);i<4;i++) {
420     o << " _meanScale[" << be << "][" << hl << "][" << i << "]=" << _meanScale[be][hl][i] << ";" << std::endl;
421     }
422     for(unsigned i(0);i<4;i++) {
423     o << " _meanAT[" << be << "][" << hl << "][" << i << "]=" << _meanAT[be][hl][i] << ";" << std::endl;
424     }
425     for(unsigned i(0);i<4;i++) {
426     o << " _meanAC[" << be << "][" << hl << "][" << i << "]=" << _meanAC[be][hl][i] << ";" << std::endl;
427     }
428     for(unsigned i(0);i<4;i++) {
429     o << " _meanAS[" << be << "][" << hl << "][" << i << "]=" << _meanAS[be][hl][i] << ";" << std::endl;
430     }
431     for(unsigned i(0);i<4;i++) {
432     o << " _meanAM[" << be << "][" << hl << "][" << i << "]=" << _meanAM[be][hl][i] << ";" << std::endl;
433     }
434     for(unsigned i(0);i<4;i++) {
435     o << " _meanBT[" << be << "][" << hl << "][" << i << "]=" << _meanBT[be][hl][i] << ";" << std::endl;
436     }
437     for(unsigned i(0);i<4;i++) {
438     o << " _meanBC[" << be << "][" << hl << "][" << i << "]=" << _meanBC[be][hl][i] << ";" << std::endl;
439     }
440     for(unsigned i(0);i<4;i++) {
441     o << " _meanBS[" << be << "][" << hl << "][" << i << "]=" << _meanBS[be][hl][i] << ";" << std::endl;
442     }
443     for(unsigned i(0);i<4;i++) {
444     o << " _meanBM[" << be << "][" << hl << "][" << i << "]=" << _meanBM[be][hl][i] << ";" << std::endl;
445     }
446     for(unsigned i(0);i<4;i++) {
447     o << " _meanR9[" << be << "][" << hl << "][" << i << "]=" << _meanR9[be][hl][i] << ";" << std::endl;
448     }
449     o << std::endl;
450    
451     for(unsigned i(0);i<4;i++) {
452     o << " _sigmaScale[" << be << "][" << hl << "][" << i << "]=" << _sigmaScale[be][hl][i] << ";" << std::endl;
453     }
454     for(unsigned i(0);i<4;i++) {
455     o << " _sigmaAT[" << be << "][" << hl << "][" << i << "]=" << _sigmaAT[be][hl][i] << ";" << std::endl;
456     }
457     for(unsigned i(0);i<4;i++) {
458     o << " _sigmaAC[" << be << "][" << hl << "][" << i << "]=" << _sigmaAC[be][hl][i] << ";" << std::endl;
459     }
460     for(unsigned i(0);i<4;i++) {
461     o << " _sigmaAS[" << be << "][" << hl << "][" << i << "]=" << _sigmaAS[be][hl][i] << ";" << std::endl;
462     }
463     for(unsigned i(0);i<4;i++) {
464     o << " _sigmaAM[" << be << "][" << hl << "][" << i << "]=" << _sigmaAM[be][hl][i] << ";" << std::endl;
465     }
466     for(unsigned i(0);i<4;i++) {
467     o << " _sigmaBT[" << be << "][" << hl << "][" << i << "]=" << _sigmaBT[be][hl][i] << ";" << std::endl;
468     }
469     for(unsigned i(0);i<4;i++) {
470     o << " _sigmaBC[" << be << "][" << hl << "][" << i << "]=" << _sigmaBC[be][hl][i] << ";" << std::endl;
471     }
472     for(unsigned i(0);i<4;i++) {
473     o << " _sigmaBS[" << be << "][" << hl << "][" << i << "]=" << _sigmaBS[be][hl][i] << ";" << std::endl;
474     }
475     for(unsigned i(0);i<4;i++) {
476     o << " _sigmaBM[" << be << "][" << hl << "][" << i << "]=" << _sigmaBM[be][hl][i] << ";" << std::endl;
477     }
478     for(unsigned i(0);i<4;i++) {
479     o << " _sigmaR9[" << be << "][" << hl << "][" << i << "]=" << _sigmaR9[be][hl][i] << ";" << std::endl;
480     }
481     o << std::endl;
482     }
483     }
484     }
485    
486     void PhotonFix::printParameters(std::ostream &o) {
487     o << "PhotonFix::printParameters()" << std::endl;
488    
489     for(unsigned be(0);be<2;be++) {
490     for(unsigned hl(0);hl<2;hl++) {
491     o << " Parameters for " << (be==0?"barrel":"endcap")
492     << ", " << (hl==0?"high":"low") << " R9" << std::endl;
493    
494     o << " Mean scaling ";
495     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanScale[be][hl][i];
496     o << std::endl;
497     o << " Mean " << (be==0?"Eta ":"ZetaX") << " total ";
498     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAT[be][hl][i];
499     o << std::endl;
500     o << " Mean " << (be==0?"Eta ":"ZetaX") << " crystal ";
501     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAC[be][hl][i];
502     o << std::endl;
503     o << " Mean " << (be==0?"Eta ":"ZetaX") << " submodule";
504     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAS[be][hl][i];
505     o << std::endl;
506     o << " Mean " << (be==0?"Eta ":"ZetaX") << " module ";
507     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAM[be][hl][i];
508     o << std::endl;
509     o << " Mean " << (be==0?"Eta zero ":"ZetaY total ");
510     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBT[be][hl][i];
511     o << std::endl;
512     o << " Mean " << (be==0?"Phi ":"ZetaY") << " crystal ";
513     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBC[be][hl][i];
514     o << std::endl;
515     o << " Mean " << (be==0?"Phi ":"ZetaY") << " submodule";
516     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBS[be][hl][i];
517     o << std::endl;
518     o << " Mean " << (be==0?"Phi ":"ZetaY") << " module ";
519     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBM[be][hl][i];
520     o << std::endl;
521     o << " Mean R9 ";
522     for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanR9[be][hl][i];
523     o << std::endl;
524    
525     o << " Sigma scaling ";
526     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaScale[be][hl][i];
527     o << std::endl;
528     o << " Sigma " << (be==0?"Eta ":"ZetaX") << " total ";
529     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAT[be][hl][i];
530     o << std::endl;
531     o << " Sigma " << (be==0?"Eta ":"ZetaX") << " crystal ";
532     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAC[be][hl][i];
533     o << std::endl;
534     o << " Sigma " << (be==0?"Eta ":"ZetaX") << " submodule";
535     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAS[be][hl][i];
536     o << std::endl;
537     o << " Sigma " << (be==0?"Eta ":"ZetaX") << " module ";
538     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAM[be][hl][i];
539     o << std::endl;
540     o << " Sigma " << (be==0?"Eta ":"ZetaY") << " total ";
541     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBT[be][hl][i];
542     o << std::endl;
543     o << " Sigma " << (be==0?"Eta ":"ZetaY") << " crystal ";
544     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBC[be][hl][i];
545     o << std::endl;
546     o << " Sigma " << (be==0?"Phi ":"ZetaY") << " submodule";
547     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBS[be][hl][i];
548     o << std::endl;
549     o << " Sigma " << (be==0?"Phi ":"ZetaY") << " module ";
550     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBM[be][hl][i];
551     o << std::endl;
552     o << " Sigma R9 ";
553     for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaR9[be][hl][i];
554     o << std::endl;
555     }
556     }
557     }
558    
559     double PhotonFix::asinh(double s) const {
560     if(s>=0.0) return log(sqrt(s*s+1.0)+s);
561     else return -log(sqrt(s*s+1.0)-s);
562     }
563    
564     void PhotonFix::dumpGaps(std::ostream &o) {
565     o << std::setprecision(15);
566    
567     for(unsigned i(0);i<169;i++) {
568     for(unsigned j(0);j<360;j++) {
569     for(unsigned k(0);k<2;k++) {
570     o << _barrelCGap[i][j][k] << std::endl;
571     }
572     }
573     }
574    
575     for(unsigned i(0);i<33;i++) {
576     for(unsigned j(0);j<180;j++) {
577     for(unsigned k(0);k<2;k++) {
578     o << _barrelSGap[i][j][k] << std::endl;
579     }
580     }
581     }
582    
583     for(unsigned i(0);i<7;i++) {
584     for(unsigned j(0);j<18;j++) {
585     for(unsigned k(0);k<2;k++) {
586     o << _barrelMGap[i][j][k] << std::endl;
587     }
588     }
589     }
590    
591     for(unsigned i(0);i<100;i++) {
592     for(unsigned j(0);j<100;j++) {
593     if(_endcapCrystal[i][j]) o << 0 << std::endl;
594     else o << 1 << std::endl;
595     }
596     }
597    
598     for(unsigned i(0);i<2;i++) {
599     for(unsigned j(0);j<7080;j++) {
600     for(unsigned k(0);k<2;k++) {
601     o << _endcapCGap[i][j][k] << std::endl;
602     }
603     }
604     }
605    
606     for(unsigned i(0);i<2;i++) {
607     for(unsigned j(0);j<264;j++) {
608     for(unsigned k(0);k<2;k++) {
609     o << _endcapSGap[i][j][k] << std::endl;
610     }
611     }
612     }
613    
614     for(unsigned i(0);i<2;i++) {
615     for(unsigned j(0);j<1;j++) {
616     for(unsigned k(0);k<2;k++) {
617     o << _endcapMGap[i][j][k] << std::endl;
618     }
619     }
620     }
621     }
622    
623     double PhotonFix::dPhi(double f0, double f1) const {
624     double df(f0-f1);
625     if(df> _onePi) df-=_twoPi;
626     if(df<-_onePi) df+=_twoPi;
627     return df;
628     }
629    
630     double PhotonFix::aPhi(double f0, double f1) const {
631     double af(0.5*(f0+f1));
632     if(fabs(dPhi(af,f0))>0.5*_onePi) {
633     if(af>=0.0) af-=_onePi;
634     else af+=_onePi;
635     }
636    
637     assert(fabs(dPhi(af,f0))<0.5*_onePi);
638     assert(fabs(dPhi(af,f1))<0.5*_onePi);
639    
640     return af;
641     }
642    
643     double PhotonFix::expCorrection(double a, const double *p) const {
644     if(p[1]==0.0 || p[2]==0.0 || p[3]==0.0) return 0.0;
645    
646     double b(a-p[0]);
647     if(b>=0.0) return p[1]*exp(-fabs(p[2])*b);
648     else return p[1]*exp( fabs(p[3])*b);
649     }
650    
651     double PhotonFix::gausCorrection(double a, const double *p) const {
652     if(p[1]==0.0 || p[2]==0.0 || p[3]==0.0) return 0.0;
653    
654     double b(a-p[0]);
655     if(b>=0.0) return p[1]*exp(-0.5*p[2]*p[2]*b*b);
656     else return p[1]*exp(-0.5*p[3]*p[3]*b*b);
657     }
658     bool PhotonFix::initialised() {
659     return _initialised;
660     }
661     bool PhotonFix::initialise(const std::string &s, const std::string &infile) {
662     if(_initialised) return false;
663    
664    
665     initialiseParameters(s);
666     initialiseGeometry(s,infile);
667     return true;
668     }
669    
670     bool PhotonFix::initialiseParameters(const std::string &s) {
671     _initialised=false;
672    
673     if(s=="Nominal") {
674     for(unsigned be(0);be<2;be++) {
675     for(unsigned hl(0);hl<2;hl++) {
676     for(unsigned i(0);i<4;i++) {
677     _meanScale[be][hl][i]=0;
678     _meanAT[be][hl][i]=0;
679     _meanAC[be][hl][i]=0;
680     _meanAS[be][hl][i]=0;
681     _meanAM[be][hl][i]=0;
682     _meanBT[be][hl][i]=0;
683     _meanBC[be][hl][i]=0;
684     _meanBS[be][hl][i]=0;
685     _meanBM[be][hl][i]=0;
686     _meanR9[be][hl][i]=0;
687    
688     _sigmaScale[be][hl][i]=0;
689     _sigmaAT[be][hl][i]=0;
690     _sigmaAC[be][hl][i]=0;
691     _sigmaAS[be][hl][i]=0;
692     _sigmaAM[be][hl][i]=0;
693     _sigmaBT[be][hl][i]=0;
694     _sigmaBC[be][hl][i]=0;
695     _sigmaBS[be][hl][i]=0;
696     _sigmaBM[be][hl][i]=0;
697     _sigmaR9[be][hl][i]=0;
698     }
699    
700     _meanScale[be][hl][0]=1.0;
701     if(be==0) {
702     _sigmaScale[be][hl][0]=0.2;
703     _sigmaScale[be][hl][1]=0.03;
704     _sigmaScale[be][hl][2]=0.006;
705     } else {
706     _sigmaScale[be][hl][0]=0.25;
707     _sigmaScale[be][hl][1]=0.05;
708     _sigmaScale[be][hl][2]=0.010;
709     }
710     }
711     }
712    
713     _initialised=true;
714     }
715    
716     if(s=="3_8") {
717     _meanScale[0][0][0]=0.994724;
718     _meanScale[0][0][1]=1.98102e-06;
719     _meanScale[0][0][2]=1.43015e-05;
720     _meanScale[0][0][3]=-0.0908525;
721     _meanAT[0][0][0]=0.0;
722     _meanAT[0][0][1]=0.0;
723     _meanAT[0][0][2]=0.0;
724     _meanAT[0][0][3]=0.0;
725     _meanAC[0][0][0]=-0.00352041;
726     _meanAC[0][0][1]=0.00982015;
727     _meanAC[0][0][2]=434.32;
728     _meanAC[0][0][3]=529.508;
729     _meanAS[0][0][0]=-1.1;
730     _meanAS[0][0][1]=0.00135995;
731     _meanAS[0][0][2]=295.712;
732     _meanAS[0][0][3]=5.13202e+07;
733     _meanAM[0][0][0]=-0.00140562;
734     _meanAM[0][0][1]=0.156322;
735     _meanAM[0][0][2]=263.097;
736     _meanAM[0][0][3]=222.294;
737     _meanBT[0][0][0]=0.0;
738     _meanBT[0][0][1]=0.0;
739     _meanBT[0][0][2]=0.0;
740     _meanBT[0][0][3]=0.0;
741     _meanBC[0][0][0]=-0.00294295;
742     _meanBC[0][0][1]=0.011533;
743     _meanBC[0][0][2]=562.905;
744     _meanBC[0][0][3]=421.097;
745     _meanBS[0][0][0]=-0.00204373;
746     _meanBS[0][0][1]=0.00347592;
747     _meanBS[0][0][2]=36.5614;
748     _meanBS[0][0][3]=1265.25;
749     _meanBM[0][0][0]=-0.00275381;
750     _meanBM[0][0][1]=0.0812447;
751     _meanBM[0][0][2]=216.885;
752     _meanBM[0][0][3]=264.754;
753     _meanR9[0][0][0]=0.952584;
754     _meanR9[0][0][1]=22.7119;
755     _meanR9[0][0][2]=402.816;
756     _meanR9[0][0][3]=0;
757    
758     _sigmaScale[0][0][0]=0.167184;
759     _sigmaScale[0][0][1]=6.14323e-11;
760     _sigmaScale[0][0][2]=0.00769693;
761     _sigmaScale[0][0][3]=0;
762     _sigmaAT[0][0][0]=0.228255;
763     _sigmaAT[0][0][1]=0;
764     _sigmaAT[0][0][2]=0;
765     _sigmaAT[0][0][3]=0;
766     _sigmaAC[0][0][0]=-0.00411906;
767     _sigmaAC[0][0][1]=0.077799;
768     _sigmaAC[0][0][2]=23.1033;
769     _sigmaAC[0][0][3]=-3e+17;
770     _sigmaAS[0][0][0]=0;
771     _sigmaAS[0][0][1]=0;
772     _sigmaAS[0][0][2]=0;
773     _sigmaAS[0][0][3]=0;
774     _sigmaAM[0][0][0]=-0.000130695;
775     _sigmaAM[0][0][1]=11.2121;
776     _sigmaAM[0][0][2]=468.535;
777     _sigmaAM[0][0][3]=407.652;
778     _sigmaBT[0][0][0]=1.33384e-05;
779     _sigmaBT[0][0][1]=8.77098;
780     _sigmaBT[0][0][2]=324.048;
781     _sigmaBT[0][0][3]=239.868;
782     _sigmaBC[0][0][0]=-0.00281964;
783     _sigmaBC[0][0][1]=0.125811;
784     _sigmaBC[0][0][2]=538.949;
785     _sigmaBC[0][0][3]=1358.76;
786     _sigmaBS[0][0][0]=0;
787     _sigmaBS[0][0][1]=0;
788     _sigmaBS[0][0][2]=0;
789     _sigmaBS[0][0][3]=0;
790     _sigmaBM[0][0][0]=-0.00293676;
791     _sigmaBM[0][0][1]=8.88276;
792     _sigmaBM[0][0][2]=350.032;
793     _sigmaBM[0][0][3]=580.354;
794     _sigmaR9[0][0][0]=0.955876;
795     _sigmaR9[0][0][1]=2254.5;
796     _sigmaR9[0][0][2]=14627;
797     _sigmaR9[0][0][3]=0;
798    
799     _meanScale[0][1][0]=0.888348;
800     _meanScale[0][1][1]=1.20452e-05;
801     _meanScale[0][1][2]=-1.04458e-05;
802     _meanScale[0][1][3]=-0.542383;
803     _meanAT[0][1][0]=0.0;
804     _meanAT[0][1][1]=0.0;
805     _meanAT[0][1][2]=0.0;
806     _meanAT[0][1][3]=0.0;
807     _meanAC[0][1][0]=-0.00320856;
808     _meanAC[0][1][1]=0.0240109;
809     _meanAC[0][1][2]=115.145;
810     _meanAC[0][1][3]=205.859;
811     _meanAS[0][1][0]=0.0349736;
812     _meanAS[0][1][1]=-0.00232864;
813     _meanAS[0][1][2]=318.584;
814     _meanAS[0][1][3]=1.4e+09;
815     _meanAM[0][1][0]=-0.00104798;
816     _meanAM[0][1][1]=0.208249;
817     _meanAM[0][1][2]=297.049;
818     _meanAM[0][1][3]=220.609;
819     _meanBT[0][1][0]=0.0;
820     _meanBT[0][1][1]=0.0;
821     _meanBT[0][1][2]=0.0;
822     _meanBT[0][1][3]=0.0;
823     _meanBC[0][1][0]=-0.00420429;
824     _meanBC[0][1][1]=0.00203991;
825     _meanBC[0][1][2]=172.278;
826     _meanBC[0][1][3]=410.677;
827     _meanBS[0][1][0]=-0.0430854;
828     _meanBS[0][1][1]=0.0961883;
829     _meanBS[0][1][2]=0.196958;
830     _meanBS[0][1][3]=11442.2;
831     _meanBM[0][1][0]=-0.00389457;
832     _meanBM[0][1][1]=0.0449086;
833     _meanBM[0][1][2]=78.9252;
834     _meanBM[0][1][3]=103.237;
835     _meanR9[0][1][0]=0.0182102;
836     _meanR9[0][1][1]=-0.03752;
837     _meanR9[0][1][2]=0.0198881;
838     _meanR9[0][1][3]=0;
839    
840     _sigmaScale[0][1][0]=0.386681;
841     _sigmaScale[0][1][1]=0.0913412;
842     _sigmaScale[0][1][2]=0.00119232;
843     _sigmaScale[0][1][3]=0;
844     _sigmaAT[0][1][0]=1.36562;
845     _sigmaAT[0][1][1]=0;
846     _sigmaAT[0][1][2]=0;
847     _sigmaAT[0][1][3]=0;
848     _sigmaAC[0][1][0]=-0.00504613;
849     _sigmaAC[0][1][1]=-1.09115;
850     _sigmaAC[0][1][2]=8.57406;
851     _sigmaAC[0][1][3]=57.1351;
852     _sigmaAS[0][1][0]=0;
853     _sigmaAS[0][1][1]=0;
854     _sigmaAS[0][1][2]=0;
855     _sigmaAS[0][1][3]=0;
856     _sigmaAM[0][1][0]=-0.00014319;
857     _sigmaAM[0][1][1]=5.39527;
858     _sigmaAM[0][1][2]=432.566;
859     _sigmaAM[0][1][3]=265.165;
860     _sigmaBT[0][1][0]=-0.040161;
861     _sigmaBT[0][1][1]=2.65711;
862     _sigmaBT[0][1][2]=-0.398357;
863     _sigmaBT[0][1][3]=-0.440649;
864     _sigmaBC[0][1][0]=0.00580015;
865     _sigmaBC[0][1][1]=-0.631833;
866     _sigmaBC[0][1][2]=18594.3;
867     _sigmaBC[0][1][3]=4.00955e+08;
868     _sigmaBS[0][1][0]=0;
869     _sigmaBS[0][1][1]=0;
870     _sigmaBS[0][1][2]=0;
871     _sigmaBS[0][1][3]=0;
872     _sigmaBM[0][1][0]=-0.00376665;
873     _sigmaBM[0][1][1]=3.74316;
874     _sigmaBM[0][1][2]=102.72;
875     _sigmaBM[0][1][3]=157.396;
876     _sigmaR9[0][1][0]=-3.12696;
877     _sigmaR9[0][1][1]=1.75114;
878     _sigmaR9[0][1][2]=0;
879     _sigmaR9[0][1][3]=0;
880    
881     _meanScale[1][0][0]=0.999461;
882     _meanScale[1][0][1]=4.37414e-06;
883     _meanScale[1][0][2]=4.92078e-06;
884     _meanScale[1][0][3]=-0.121609;
885     _meanAT[1][0][0]=0.0;
886     _meanAT[1][0][1]=0.0;
887     _meanAT[1][0][2]=0.0;
888     _meanAT[1][0][3]=0.0;
889     _meanAC[1][0][0]=-0.000396058;
890     _meanAC[1][0][1]=0.0144837;
891     _meanAC[1][0][2]=1374.93;
892     _meanAC[1][0][3]=945.634;
893     _meanAS[1][0][0]=-0.000871036;
894     _meanAS[1][0][1]=0.0442747;
895     _meanAS[1][0][2]=645.709;
896     _meanAS[1][0][3]=962.845;
897     _meanAM[1][0][0]=0.000434298;
898     _meanAM[1][0][1]=0.0658628;
899     _meanAM[1][0][2]=1928.49;
900     _meanAM[1][0][3]=728.522;
901     _meanBT[1][0][0]=0.0;
902     _meanBT[1][0][1]=0.0;
903     _meanBT[1][0][2]=0.0;
904     _meanBT[1][0][3]=0.0;
905     _meanBC[1][0][0]=-0.000452212;
906     _meanBC[1][0][1]=0.0129968;
907     _meanBC[1][0][2]=1056.08;
908     _meanBC[1][0][3]=759.102;
909     _meanBS[1][0][0]=-0.000786157;
910     _meanBS[1][0][1]=0.0346555;
911     _meanBS[1][0][2]=592.239;
912     _meanBS[1][0][3]=854.285;
913     _meanBM[1][0][0]=-0.0665038;
914     _meanBM[1][0][1]=-0.00211713;
915     _meanBM[1][0][2]=4.84395;
916     _meanBM[1][0][3]=11.6644;
917     _meanR9[1][0][0]=0.971355;
918     _meanR9[1][0][1]=47.2751;
919     _meanR9[1][0][2]=536.907;
920     _meanR9[1][0][3]=0;
921    
922     _sigmaScale[1][0][0]=0.254641;
923     _sigmaScale[1][0][1]=0.00264818;
924     _sigmaScale[1][0][2]=0.0114953;
925     _sigmaScale[1][0][3]=0;
926     _sigmaAT[1][0][0]=0.935839;
927     _sigmaAT[1][0][1]=0;
928     _sigmaAT[1][0][2]=0;
929     _sigmaAT[1][0][3]=0;
930     _sigmaAC[1][0][0]=-0.00476475;
931     _sigmaAC[1][0][1]=2.14548;
932     _sigmaAC[1][0][2]=29937;
933     _sigmaAC[1][0][3]=2.6e+11;
934     _sigmaAS[1][0][0]=-8.17285e-05;
935     _sigmaAS[1][0][1]=1.5821;
936     _sigmaAS[1][0][2]=1928.83;
937     _sigmaAS[1][0][3]=902.519;
938     _sigmaAM[1][0][0]=0.0278577;
939     _sigmaAM[1][0][1]=0.58439;
940     _sigmaAM[1][0][2]=43.3575;
941     _sigmaAM[1][0][3]=19.7836;
942     _sigmaBT[1][0][0]=-0.456051;
943     _sigmaBT[1][0][1]=0;
944     _sigmaBT[1][0][2]=0;
945     _sigmaBT[1][0][3]=0;
946     _sigmaBC[1][0][0]=-0.00264527;
947     _sigmaBC[1][0][1]=0.696043;
948     _sigmaBC[1][0][2]=7.49509e+12;
949     _sigmaBC[1][0][3]=96843;
950     _sigmaBS[1][0][0]=0.000258933;
951     _sigmaBS[1][0][1]=1.28387;
952     _sigmaBS[1][0][2]=1668.71;
953     _sigmaBS[1][0][3]=730.716;
954     _sigmaBM[1][0][0]=0.00121506;
955     _sigmaBM[1][0][1]=0.938541;
956     _sigmaBM[1][0][2]=9003.57;
957     _sigmaBM[1][0][3]=288.897;
958     _sigmaR9[1][0][0]=1.01207;
959     _sigmaR9[1][0][1]=-816.244;
960     _sigmaR9[1][0][2]=-16283.8;
961     _sigmaR9[1][0][3]=0;
962    
963     _meanScale[1][1][0]=0.324634;
964     _meanScale[1][1][1]=9.48206e-05;
965     _meanScale[1][1][2]=1.0e-12;
966     _meanScale[1][1][3]=1.0e-12;
967     _meanAT[1][1][0]=0.0;
968     _meanAT[1][1][1]=0.0;
969     _meanAT[1][1][2]=0.0;
970     _meanAT[1][1][3]=0.0;
971     _meanAC[1][1][0]=-0.00158311;
972     _meanAC[1][1][1]=0.0106161;
973     _meanAC[1][1][2]=338.964;
974     _meanAC[1][1][3]=797.172;
975     _meanAS[1][1][0]=-0.00960269;
976     _meanAS[1][1][1]=-0.00496491;
977     _meanAS[1][1][2]=934.472;
978     _meanAS[1][1][3]=8.32667e-16;
979     _meanAM[1][1][0]=-0.00219814;
980     _meanAM[1][1][1]=0.653906;
981     _meanAM[1][1][2]=0.0949848;
982     _meanAM[1][1][3]=0.0977831;
983     _meanBT[1][1][0]=0.0;
984     _meanBT[1][1][1]=0.0;
985     _meanBT[1][1][2]=0.0;
986     _meanBT[1][1][3]=0.0;
987     _meanBC[1][1][0]=-0.00423472;
988     _meanBC[1][1][1]=0.0279695;
989     _meanBC[1][1][2]=28073.7;
990     _meanBC[1][1][3]=118612;
991     _meanBS[1][1][0]=-0.0012476;
992     _meanBS[1][1][1]=0.02744;
993     _meanBS[1][1][2]=390.697;
994     _meanBS[1][1][3]=727.861;
995     _meanBM[1][1][0]=-1.36573e-05;
996     _meanBM[1][1][1]=0.0667504;
997     _meanBM[1][1][2]=-80154.4;
998     _meanBM[1][1][3]=576.637;
999     _meanR9[1][1][0]=0.113317;
1000     _meanR9[1][1][1]=0.0142669;
1001     _meanR9[1][1][2]=-0.125721;
1002     _meanR9[1][1][3]=0;
1003    
1004     _sigmaScale[1][1][0]=0.471767;
1005     _sigmaScale[1][1][1]=0.211196;
1006     _sigmaScale[1][1][2]=0.0240124;
1007     _sigmaScale[1][1][3]=0;
1008     _sigmaAT[1][1][0]=0.404395;
1009     _sigmaAT[1][1][1]=0;
1010     _sigmaAT[1][1][2]=0;
1011     _sigmaAT[1][1][3]=0;
1012     _sigmaAC[1][1][0]=0.00173151;
1013     _sigmaAC[1][1][1]=-0.479291;
1014     _sigmaAC[1][1][2]=11583.5;
1015     _sigmaAC[1][1][3]=-7e+09;
1016     _sigmaAS[1][1][0]=0.000450387;
1017     _sigmaAS[1][1][1]=0.662978;
1018     _sigmaAS[1][1][2]=924.051;
1019     _sigmaAS[1][1][3]=448.417;
1020     _sigmaAM[1][1][0]=0.00335603;
1021     _sigmaAM[1][1][1]=0.648407;
1022     _sigmaAM[1][1][2]=134.672;
1023     _sigmaAM[1][1][3]=27.4139;
1024     _sigmaBT[1][1][0]=0.602402;
1025     _sigmaBT[1][1][1]=0;
1026     _sigmaBT[1][1][2]=0;
1027     _sigmaBT[1][1][3]=0;
1028     _sigmaBC[1][1][0]=-0.00256192;
1029     _sigmaBC[1][1][1]=2.01276;
1030     _sigmaBC[1][1][2]=114558;
1031     _sigmaBC[1][1][3]=2.15421e+06;
1032     _sigmaBS[1][1][0]=0.00151576;
1033     _sigmaBS[1][1][1]=0.359084;
1034     _sigmaBS[1][1][2]=329.414;
1035     _sigmaBS[1][1][3]=154.509;
1036     _sigmaBM[1][1][0]=-0.0452587;
1037     _sigmaBM[1][1][1]=1.26253;
1038     _sigmaBM[1][1][2]=1.9e+09;
1039     _sigmaBM[1][1][3]=1058.76;
1040     _sigmaR9[1][1][0]=4.59667;
1041     _sigmaR9[1][1][1]=-5.14404;
1042     _sigmaR9[1][1][2]=0;
1043     _sigmaR9[1][1][3]=0;
1044    
1045     _initialised=true;
1046     }
1047    
1048     if(s=="3_11") {
1049     _meanScale[0][0][0]=0.994363;
1050     _meanScale[0][0][1]=4.84904e-07;
1051     _meanScale[0][0][2]=1.54475e-05;
1052     _meanScale[0][0][3]=-0.103309;
1053     _meanAT[0][0][0]=0.0;
1054     _meanAT[0][0][1]=0.0;
1055     _meanAT[0][0][2]=0.0;
1056     _meanAT[0][0][3]=0.0;
1057     _meanAC[0][0][0]=-0.00360057;
1058     _meanAC[0][0][1]=0.00970858;
1059     _meanAC[0][0][2]=409.406;
1060     _meanAC[0][0][3]=527.952;
1061     _meanAS[0][0][0]=-1.1;
1062     _meanAS[0][0][1]=0.00135995;
1063     _meanAS[0][0][2]=295.712;
1064     _meanAS[0][0][3]=5.13202e+07;
1065     _meanAM[0][0][0]=-0.00129854;
1066     _meanAM[0][0][1]=0.151466;
1067     _meanAM[0][0][2]=261.828;
1068     _meanAM[0][0][3]=214.662;
1069     _meanBT[0][0][0]=0.0;
1070     _meanBT[0][0][1]=0.0;
1071     _meanBT[0][0][2]=0.0;
1072     _meanBT[0][0][3]=0.0;
1073     _meanBC[0][0][0]=-0.00286864;
1074     _meanBC[0][0][1]=0.0114118;
1075     _meanBC[0][0][2]=563.962;
1076     _meanBC[0][0][3]=412.922;
1077     _meanBS[0][0][0]=-0.00210996;
1078     _meanBS[0][0][1]=0.00327867;
1079     _meanBS[0][0][2]=23.617;
1080     _meanBS[0][0][3]=1018.45;
1081     _meanBM[0][0][0]=-0.002287;
1082     _meanBM[0][0][1]=0.0848984;
1083     _meanBM[0][0][2]=235.575;
1084     _meanBM[0][0][3]=260.773;
1085     _meanR9[0][0][0]=0.951724;
1086     _meanR9[0][0][1]=23.7181;
1087     _meanR9[0][0][2]=177.34;
1088     _meanR9[0][0][3]=0;
1089    
1090     _sigmaScale[0][0][0]=0.187578;
1091     _sigmaScale[0][0][1]=-0.000901045;
1092     _sigmaScale[0][0][2]=0.00673186;
1093     _sigmaScale[0][0][3]=0;
1094     _sigmaAT[0][0][0]=0.183777;
1095     _sigmaAT[0][0][1]=0;
1096     _sigmaAT[0][0][2]=0;
1097     _sigmaAT[0][0][3]=0;
1098     _sigmaAC[0][0][0]=-0.00430202;
1099     _sigmaAC[0][0][1]=0.122501;
1100     _sigmaAC[0][0][2]=51.9772;
1101     _sigmaAC[0][0][3]=-3e+17;
1102     _sigmaAS[0][0][0]=0;
1103     _sigmaAS[0][0][1]=0;
1104     _sigmaAS[0][0][2]=0;
1105     _sigmaAS[0][0][3]=0;
1106     _sigmaAM[0][0][0]=0.00101883;
1107     _sigmaAM[0][0][1]=11.2009;
1108     _sigmaAM[0][0][2]=593.111;
1109     _sigmaAM[0][0][3]=345.433;
1110     _sigmaBT[0][0][0]=-6.02356e-05;
1111     _sigmaBT[0][0][1]=6.99896;
1112     _sigmaBT[0][0][2]=235.996;
1113     _sigmaBT[0][0][3]=196;
1114     _sigmaBC[0][0][0]=-0.00282254;
1115     _sigmaBC[0][0][1]=0.18764;
1116     _sigmaBC[0][0][2]=509.825;
1117     _sigmaBC[0][0][3]=1400.14;
1118     _sigmaBS[0][0][0]=0;
1119     _sigmaBS[0][0][1]=0;
1120     _sigmaBS[0][0][2]=0;
1121     _sigmaBS[0][0][3]=0;
1122     _sigmaBM[0][0][0]=-0.00252199;
1123     _sigmaBM[0][0][1]=39.1544;
1124     _sigmaBM[0][0][2]=612.481;
1125     _sigmaBM[0][0][3]=905.994;
1126     _sigmaR9[0][0][0]=0.95608;
1127     _sigmaR9[0][0][1]=2203.31;
1128     _sigmaR9[0][0][2]=-22454.2;
1129     _sigmaR9[0][0][3]=0;
1130    
1131     _meanScale[0][1][0]=0.889415;
1132     _meanScale[0][1][1]=1.21788e-05;
1133     _meanScale[0][1][2]=-4.3438e-06;
1134     _meanScale[0][1][3]=-0.629968;
1135     _meanAT[0][1][0]=0.0;
1136     _meanAT[0][1][1]=0.0;
1137     _meanAT[0][1][2]=0.0;
1138     _meanAT[0][1][3]=0.0;
1139     _meanAC[0][1][0]=-0.00313701;
1140     _meanAC[0][1][1]=0.0227998;
1141     _meanAC[0][1][2]=128.653;
1142     _meanAC[0][1][3]=234.333;
1143     _meanAS[0][1][0]=0.0346198;
1144     _meanAS[0][1][1]=-0.00261336;
1145     _meanAS[0][1][2]=177.983;
1146     _meanAS[0][1][3]=1.19839e+14;
1147     _meanAM[0][1][0]=-0.00100745;
1148     _meanAM[0][1][1]=0.264247;
1149     _meanAM[0][1][2]=337.255;
1150     _meanAM[0][1][3]=251.454;
1151     _meanBT[0][1][0]=0.0;
1152     _meanBT[0][1][1]=0.0;
1153     _meanBT[0][1][2]=0.0;
1154     _meanBT[0][1][3]=0.0;
1155     _meanBC[0][1][0]=-0.00397794;
1156     _meanBC[0][1][1]=0.00219079;
1157     _meanBC[0][1][2]=176.842;
1158     _meanBC[0][1][3]=450.29;
1159     _meanBS[0][1][0]=-2e+07;
1160     _meanBS[0][1][1]=0.0957598;
1161     _meanBS[0][1][2]=-8.88573e-27;
1162     _meanBS[0][1][3]=11442.2;
1163     _meanBM[0][1][0]=-0.00366315;
1164     _meanBM[0][1][1]=0.0622186;
1165     _meanBM[0][1][2]=94.5155;
1166     _meanBM[0][1][3]=126.404;
1167     _meanR9[0][1][0]=0.00636789;
1168     _meanR9[0][1][1]=0.000336062;
1169     _meanR9[0][1][2]=-0.0092699;
1170     _meanR9[0][1][3]=0;
1171    
1172     _sigmaScale[0][1][0]=0.685096;
1173     _sigmaScale[0][1][1]=0.129065;
1174     _sigmaScale[0][1][2]=-0.00212486;
1175     _sigmaScale[0][1][3]=0;
1176     _sigmaAT[0][1][0]=0.898865;
1177     _sigmaAT[0][1][1]=0;
1178     _sigmaAT[0][1][2]=0;
1179     _sigmaAT[0][1][3]=0;
1180     _sigmaAC[0][1][0]=-0.00492979;
1181     _sigmaAC[0][1][1]=-1.20123;
1182     _sigmaAC[0][1][2]=2.89231;
1183     _sigmaAC[0][1][3]=18.2059;
1184     _sigmaAS[0][1][0]=0;
1185     _sigmaAS[0][1][1]=0;
1186     _sigmaAS[0][1][2]=0;
1187     _sigmaAS[0][1][3]=0;
1188     _sigmaAM[0][1][0]=-0.000727825;
1189     _sigmaAM[0][1][1]=8.42395;
1190     _sigmaAM[0][1][2]=512.032;
1191     _sigmaAM[0][1][3]=415.962;
1192     _sigmaBT[0][1][0]=-0.0336364;
1193     _sigmaBT[0][1][1]=2.45182;
1194     _sigmaBT[0][1][2]=-0.284353;
1195     _sigmaBT[0][1][3]=-0.31679;
1196     _sigmaBC[0][1][0]=0.00510553;
1197     _sigmaBC[0][1][1]=-0.953869;
1198     _sigmaBC[0][1][2]=113872;
1199     _sigmaBC[0][1][3]=1.35966e+09;
1200     _sigmaBS[0][1][0]=0;
1201     _sigmaBS[0][1][1]=0;
1202     _sigmaBS[0][1][2]=0;
1203     _sigmaBS[0][1][3]=0;
1204     _sigmaBM[0][1][0]=-0.0034071;
1205     _sigmaBM[0][1][1]=4.19719;
1206     _sigmaBM[0][1][2]=128.952;
1207     _sigmaBM[0][1][3]=180.604;
1208     _sigmaR9[0][1][0]=-3.38988;
1209     _sigmaR9[0][1][1]=2.0714;
1210     _sigmaR9[0][1][2]=0;
1211     _sigmaR9[0][1][3]=0;
1212    
1213     _meanScale[1][0][0]=1.0009;
1214     _meanScale[1][0][1]=-4.79805e-06;
1215     _meanScale[1][0][2]=3.34625e-05;
1216     _meanScale[1][0][3]=-0.194267;
1217     _meanAT[1][0][0]=0.0;
1218     _meanAT[1][0][1]=0.0;
1219     _meanAT[1][0][2]=0.0;
1220     _meanAT[1][0][3]=0.0;
1221     _meanAC[1][0][0]=-0.000177563;
1222     _meanAC[1][0][1]=0.0122839;
1223     _meanAC[1][0][2]=1798.92;
1224     _meanAC[1][0][3]=776.856;
1225     _meanAS[1][0][0]=-0.000533039;
1226     _meanAS[1][0][1]=0.0642604;
1227     _meanAS[1][0][2]=969.596;
1228     _meanAS[1][0][3]=1004.15;
1229     _meanAM[1][0][0]=0.000163185;
1230     _meanAM[1][0][1]=0.085936;
1231     _meanAM[1][0][2]=1593.17;
1232     _meanAM[1][0][3]=681.623;
1233     _meanBT[1][0][0]=0.0;
1234     _meanBT[1][0][1]=0.0;
1235     _meanBT[1][0][2]=0.0;
1236     _meanBT[1][0][3]=0.0;
1237     _meanBC[1][0][0]=-0.000518186;
1238     _meanBC[1][0][1]=0.0121868;
1239     _meanBC[1][0][2]=1112.53;
1240     _meanBC[1][0][3]=933.281;
1241     _meanBS[1][0][0]=-0.000750734;
1242     _meanBS[1][0][1]=0.03859;
1243     _meanBS[1][0][2]=547.579;
1244     _meanBS[1][0][3]=775.887;
1245     _meanBM[1][0][0]=-0.190395;
1246     _meanBM[1][0][1]=-0.00362647;
1247     _meanBM[1][0][2]=5.25687;
1248     _meanBM[1][0][3]=-2.8e+08;
1249     _meanR9[1][0][0]=0.972346;
1250     _meanR9[1][0][1]=53.9185;
1251     _meanR9[1][0][2]=1354.5;
1252     _meanR9[1][0][3]=0;
1253    
1254     _sigmaScale[1][0][0]=0.348019;
1255     _sigmaScale[1][0][1]=-6.43731e-11;
1256     _sigmaScale[1][0][2]=0.0158647;
1257     _sigmaScale[1][0][3]=0;
1258     _sigmaAT[1][0][0]=0.215239;
1259     _sigmaAT[1][0][1]=0;
1260     _sigmaAT[1][0][2]=0;
1261     _sigmaAT[1][0][3]=0;
1262     _sigmaAC[1][0][0]=-0.00492298;
1263     _sigmaAC[1][0][1]=-3.40058;
1264     _sigmaAC[1][0][2]=17263.9;
1265     _sigmaAC[1][0][3]=2.6e+11;
1266     _sigmaAS[1][0][0]=-0.000237998;
1267     _sigmaAS[1][0][1]=3.0258;
1268     _sigmaAS[1][0][2]=1811.25;
1269     _sigmaAS[1][0][3]=1846.79;
1270     _sigmaAM[1][0][0]=0.0210134;
1271     _sigmaAM[1][0][1]=0.328359;
1272     _sigmaAM[1][0][2]=22.49;
1273     _sigmaAM[1][0][3]=14.5021;
1274     _sigmaBT[1][0][0]=-0.495072;
1275     _sigmaBT[1][0][1]=0;
1276     _sigmaBT[1][0][2]=0;
1277     _sigmaBT[1][0][3]=0;
1278     _sigmaBC[1][0][0]=-0.00265007;
1279     _sigmaBC[1][0][1]=0.970549;
1280     _sigmaBC[1][0][2]=-6.89119e+07;
1281     _sigmaBC[1][0][3]=180110;
1282     _sigmaBS[1][0][0]=0.00045833;
1283     _sigmaBS[1][0][1]=2.16342;
1284     _sigmaBS[1][0][2]=3582.4;
1285     _sigmaBS[1][0][3]=1100.36;
1286     _sigmaBM[1][0][0]=0.00188871;
1287     _sigmaBM[1][0][1]=1.66177;
1288     _sigmaBM[1][0][2]=3.2e+08;
1289     _sigmaBM[1][0][3]=2163.81;
1290     _sigmaR9[1][0][0]=-220.415;
1291     _sigmaR9[1][0][1]=5.19136e-08;
1292     _sigmaR9[1][0][2]=3.04028e-10;
1293     _sigmaR9[1][0][3]=0;
1294    
1295     _meanScale[1][1][0]=0.338011;
1296     _meanScale[1][1][1]=9.47815e-05;
1297     _meanScale[1][1][2]=-0.000238735;
1298     _meanScale[1][1][3]=-0.846414;
1299     _meanAT[1][1][0]=0.0;
1300     _meanAT[1][1][1]=0.0;
1301     _meanAT[1][1][2]=0.0;
1302     _meanAT[1][1][3]=0.0;
1303     _meanAC[1][1][0]=-0.00125367;
1304     _meanAC[1][1][1]=0.013324;
1305     _meanAC[1][1][2]=203.988;
1306     _meanAC[1][1][3]=431.951;
1307     _meanAS[1][1][0]=0.000282607;
1308     _meanAS[1][1][1]=0.0307431;
1309     _meanAS[1][1][2]=343.509;
1310     _meanAS[1][1][3]=274.957;
1311     _meanAM[1][1][0]=0.0020258;
1312     _meanAM[1][1][1]=0.643913;
1313     _meanAM[1][1][2]=0.0693877;
1314     _meanAM[1][1][3]=0.0816029;
1315     _meanBT[1][1][0]=0.0;
1316     _meanBT[1][1][1]=0.0;
1317     _meanBT[1][1][2]=0.0;
1318     _meanBT[1][1][3]=0.0;
1319     _meanBC[1][1][0]=-0.00513833;
1320     _meanBC[1][1][1]=5.94424e+08;
1321     _meanBC[1][1][2]=-62814.9;
1322     _meanBC[1][1][3]=118612;
1323     _meanBS[1][1][0]=-0.00152129;
1324     _meanBS[1][1][1]=0.0234694;
1325     _meanBS[1][1][2]=186.483;
1326     _meanBS[1][1][3]=754.201;
1327     _meanBM[1][1][0]=-0.000404987;
1328     _meanBM[1][1][1]=0.156384;
1329     _meanBM[1][1][2]=-1.7e+08;
1330     _meanBM[1][1][3]=1793.83;
1331     _meanR9[1][1][0]=0.0645278;
1332     _meanR9[1][1][1]=0.161614;
1333     _meanR9[1][1][2]=-0.215822;
1334     _meanR9[1][1][3]=0;
1335    
1336     _sigmaScale[1][1][0]=1.07376;
1337     _sigmaScale[1][1][1]=7.47238e-13;
1338     _sigmaScale[1][1][2]=0.0289594;
1339     _sigmaScale[1][1][3]=0;
1340     _sigmaAT[1][1][0]=-0.520907;
1341     _sigmaAT[1][1][1]=0;
1342     _sigmaAT[1][1][2]=0;
1343     _sigmaAT[1][1][3]=0;
1344     _sigmaAC[1][1][0]=0.00165941;
1345     _sigmaAC[1][1][1]=-0.351422;
1346     _sigmaAC[1][1][2]=8968.94;
1347     _sigmaAC[1][1][3]=-7e+09;
1348     _sigmaAS[1][1][0]=0.000490279;
1349     _sigmaAS[1][1][1]=0.554531;
1350     _sigmaAS[1][1][2]=469.111;
1351     _sigmaAS[1][1][3]=457.541;
1352     _sigmaAM[1][1][0]=0.00102079;
1353     _sigmaAM[1][1][1]=0.628055;
1354     _sigmaAM[1][1][2]=53.9452;
1355     _sigmaAM[1][1][3]=72.911;
1356     _sigmaBT[1][1][0]=-0.461542;
1357     _sigmaBT[1][1][1]=0;
1358     _sigmaBT[1][1][2]=0;
1359     _sigmaBT[1][1][3]=0;
1360     _sigmaBC[1][1][0]=-0.00219303;
1361     _sigmaBC[1][1][1]=0.874327;
1362     _sigmaBC[1][1][2]=71353.2;
1363     _sigmaBC[1][1][3]=2.09924e+08;
1364     _sigmaBS[1][1][0]=0.00104021;
1365     _sigmaBS[1][1][1]=0.236098;
1366     _sigmaBS[1][1][2]=482.954;
1367     _sigmaBS[1][1][3]=191.984;
1368     _sigmaBM[1][1][0]=-0.000116086;
1369     _sigmaBM[1][1][1]=2.4438;
1370     _sigmaBM[1][1][2]=1.9e+09;
1371     _sigmaBM[1][1][3]=-700.271;
1372     _sigmaR9[1][1][0]=4.59374;
1373     _sigmaR9[1][1][1]=-5.06202;
1374     _sigmaR9[1][1][2]=0;
1375     _sigmaR9[1][1][3]=0;
1376    
1377     _initialised=true;
1378     }
1379    
1380     if(s=="4_2") {
1381     _meanScale[0][0][0]=0.995941423;
1382     _meanScale[0][0][1]=-1.41986304e-05;
1383     _meanScale[0][0][2]=3.66129541e-05;
1384     _meanScale[0][0][3]=-0.0774047233;
1385     _meanAT[0][0][0]=0.000720281545;
1386     _meanAT[0][0][1]=0;
1387     _meanAT[0][0][2]=0;
1388     _meanAT[0][0][3]=0;
1389     _meanAC[0][0][0]=-0.00344862444;
1390     _meanAC[0][0][1]=0.0101395802;
1391     _meanAC[0][0][2]=466.112225;
1392     _meanAC[0][0][3]=507.628173;
1393     _meanAS[0][0][0]=0;
1394     _meanAS[0][0][1]=0;
1395     _meanAS[0][0][2]=0;
1396     _meanAS[0][0][3]=0;
1397     _meanAM[0][0][0]=-0.000871553792;
1398     _meanAM[0][0][1]=0.141419889;
1399     _meanAM[0][0][2]=281.104504;
1400     _meanAM[0][0][3]=195.875679;
1401     _meanBT[0][0][0]=0;
1402     _meanBT[0][0][1]=0.026344491;
1403     _meanBT[0][0][2]=-104.20518;
1404     _meanBT[0][0][3]=-176099;
1405     _meanBC[0][0][0]=-0.00272095949;
1406     _meanBC[0][0][1]=0.012411788;
1407     _meanBC[0][0][2]=587.318903;
1408     _meanBC[0][0][3]=381.415059;
1409     _meanBS[0][0][0]=-0.00201265145;
1410     _meanBS[0][0][1]=0.00372948657;
1411     _meanBS[0][0][2]=41.2773112;
1412     _meanBS[0][0][3]=748.890936;
1413     _meanBM[0][0][0]=-0.00168471013;
1414     _meanBM[0][0][1]=0.0685484442;
1415     _meanBM[0][0][2]=217.983503;
1416     _meanBM[0][0][3]=207.660928;
1417     _meanR9[0][0][0]=0.946581139;
1418     _meanR9[0][0][1]=20.6034189;
1419     _meanR9[0][0][2]=187.28856;
1420     _meanR9[0][0][3]=0;
1421    
1422     _sigmaScale[0][0][0]=0.206349443;
1423     _sigmaScale[0][0][1]=0.0206592338;
1424     _sigmaScale[0][0][2]=0.00653752299;
1425     _sigmaScale[0][0][3]=0;
1426     _sigmaAT[0][0][0]=0.178629422;
1427     _sigmaAT[0][0][1]=0;
1428     _sigmaAT[0][0][2]=0;
1429     _sigmaAT[0][0][3]=0;
1430     _sigmaAC[0][0][0]=-0.00335501889;
1431     _sigmaAC[0][0][1]=0.0997921532;
1432     _sigmaAC[0][0][2]=93.6397821;
1433     _sigmaAC[0][0][3]=1519.43272;
1434     _sigmaAS[0][0][0]=0;
1435     _sigmaAS[0][0][1]=0;
1436     _sigmaAS[0][0][2]=0;
1437     _sigmaAS[0][0][3]=0;
1438     _sigmaAM[0][0][0]=0.000927325527;
1439     _sigmaAM[0][0][1]=10.2678389;
1440     _sigmaAM[0][0][2]=619.975988;
1441     _sigmaAM[0][0][3]=285.190815;
1442     _sigmaBT[0][0][0]=0;
1443     _sigmaBT[0][0][1]=0.895041707;
1444     _sigmaBT[0][0][2]=94.6834192;
1445     _sigmaBT[0][0][3]=62.3012502;
1446     _sigmaBC[0][0][0]=-0.00169896783;
1447     _sigmaBC[0][0][1]=0.323973706;
1448     _sigmaBC[0][0][2]=1234.03309;
1449     _sigmaBC[0][0][3]=907.352988;
1450     _sigmaBS[0][0][0]=0;
1451     _sigmaBS[0][0][1]=0;
1452     _sigmaBS[0][0][2]=0;
1453     _sigmaBS[0][0][3]=0;
1454     _sigmaBM[0][0][0]=-0.00249508825;
1455     _sigmaBM[0][0][1]=57.8982306;
1456     _sigmaBM[0][0][2]=665.068952;
1457     _sigmaBM[0][0][3]=1075.1094;
1458     _sigmaR9[0][0][0]=0.952890416;
1459     _sigmaR9[0][0][1]=1958.37946;
1460     _sigmaR9[0][0][2]=21612.0219;
1461     _sigmaR9[0][0][3]=0;
1462    
1463     _meanScale[0][1][0]=0.982680412;
1464     _meanScale[0][1][1]=3.13860176e-05;
1465     _meanScale[0][1][2]=-2.89107109e-05;
1466     _meanScale[0][1][3]=-0.458678502;
1467     _meanAT[0][1][0]=-0.00204222443;
1468     _meanAT[0][1][1]=0;
1469     _meanAT[0][1][2]=0;
1470     _meanAT[0][1][3]=0;
1471     _meanAC[0][1][0]=-0.00329797061;
1472     _meanAC[0][1][1]=0.0212879256;
1473     _meanAC[0][1][2]=135.879912;
1474     _meanAC[0][1][3]=238.247576;
1475     _meanAS[0][1][0]=0;
1476     _meanAS[0][1][1]=0;
1477     _meanAS[0][1][2]=0;
1478     _meanAS[0][1][3]=0;
1479     _meanAM[0][1][0]=-0.000512006976;
1480     _meanAM[0][1][1]=0.124281288;
1481     _meanAM[0][1][2]=480.326634;
1482     _meanAM[0][1][3]=286.165783;
1483     _meanBT[0][1][0]=0;
1484     _meanBT[0][1][1]=0.204384889;
1485     _meanBT[0][1][2]=303.764745;
1486     _meanBT[0][1][3]=408.14741;
1487     _meanBC[0][1][0]=-0.0035698745;
1488     _meanBC[0][1][1]=0.00402323151;
1489     _meanBC[0][1][2]=980.296598;
1490     _meanBC[0][1][3]=869.711616;
1491     _meanBS[0][1][0]=0;
1492     _meanBS[0][1][1]=0;
1493     _meanBS[0][1][2]=0;
1494     _meanBS[0][1][3]=0;
1495     _meanBM[0][1][0]=-0.00321305828;
1496     _meanBM[0][1][1]=0.0454848819;
1497     _meanBM[0][1][2]=147.827487;
1498     _meanBM[0][1][3]=227.625382;
1499     _meanR9[0][1][0]=0.0253777359;
1500     _meanR9[0][1][1]=-0.0420810898;
1501     _meanR9[0][1][2]=0.0181966013;
1502     _meanR9[0][1][3]=0;
1503    
1504     _sigmaScale[0][1][0]=1.53707929;
1505     _sigmaScale[0][1][1]=0.0946423194;
1506     _sigmaScale[0][1][2]=-0.00765920151;
1507     _sigmaScale[0][1][3]=0;
1508     _sigmaAT[0][1][0]=0.808880052;
1509     _sigmaAT[0][1][1]=0;
1510     _sigmaAT[0][1][2]=0;
1511     _sigmaAT[0][1][3]=0;
1512     _sigmaAC[0][1][0]=-0.00195542375;
1513     _sigmaAC[0][1][1]=-2.09949949;
1514     _sigmaAC[0][1][2]=4.30292193;
1515     _sigmaAC[0][1][3]=5.09475964;
1516     _sigmaAS[0][1][0]=0;
1517     _sigmaAS[0][1][1]=0;
1518     _sigmaAS[0][1][2]=0;
1519     _sigmaAS[0][1][3]=0;
1520     _sigmaAM[0][1][0]=-0.00105652021;
1521     _sigmaAM[0][1][1]=5.83420851;
1522     _sigmaAM[0][1][2]=506.986527;
1523     _sigmaAM[0][1][3]=468.330744;
1524     _sigmaBT[0][1][0]=0;
1525     _sigmaBT[0][1][1]=2.83411417;
1526     _sigmaBT[0][1][2]=-0.211242292;
1527     _sigmaBT[0][1][3]=-0.198231087;
1528     _sigmaBC[0][1][0]=0.00580038243;
1529     _sigmaBC[0][1][1]=0.165505659;
1530     _sigmaBC[0][1][2]=4133.45418;
1531     _sigmaBC[0][1][3]=375000000;
1532     _sigmaBS[0][1][0]=0;
1533     _sigmaBS[0][1][1]=0;
1534     _sigmaBS[0][1][2]=0;
1535     _sigmaBS[0][1][3]=0;
1536     _sigmaBM[0][1][0]=-0.00269993666;
1537     _sigmaBM[0][1][1]=3.42390459;
1538     _sigmaBM[0][1][2]=171.300481;
1539     _sigmaBM[0][1][3]=284.718025;
1540     _sigmaR9[0][1][0]=-3.75255938;
1541     _sigmaR9[0][1][1]=4.3849733;
1542     _sigmaR9[0][1][2]=-1.81745726;
1543     _sigmaR9[0][1][3]=0;
1544    
1545     _meanScale[1][0][0]=0.990082016;
1546     _meanScale[1][0][1]=-3.75802712e-06;
1547     _meanScale[1][0][2]=2.56693516e-05;
1548     _meanScale[1][0][3]=-0.0492813428;
1549     _meanAT[1][0][0]=0.072352478;
1550     _meanAT[1][0][1]=0;
1551     _meanAT[1][0][2]=0;
1552     _meanAT[1][0][3]=0;
1553     _meanAC[1][0][0]=-0.0002936899;
1554     _meanAC[1][0][1]=0.0160546814;
1555     _meanAC[1][0][2]=1183.48593;
1556     _meanAC[1][0][3]=761.29774;
1557     _meanAS[1][0][0]=-0.000462243216;
1558     _meanAS[1][0][1]=0.0795658256;
1559     _meanAS[1][0][2]=887.080242;
1560     _meanAS[1][0][3]=1067.72442;
1561     _meanAM[1][0][0]=0.000354495505;
1562     _meanAM[1][0][1]=0.516700576;
1563     _meanAM[1][0][2]=4376.14811;
1564     _meanAM[1][0][3]=2093.33478;
1565     _meanBT[1][0][0]=0.077752944;
1566     _meanBT[1][0][1]=0;
1567     _meanBT[1][0][2]=0;
1568     _meanBT[1][0][3]=0;
1569     _meanBC[1][0][0]=-0.000411367107;
1570     _meanBC[1][0][1]=0.0161135906;
1571     _meanBC[1][0][2]=1414.07982;
1572     _meanBC[1][0][3]=951.556042;
1573     _meanBS[1][0][0]=8.51070829e-05;
1574     _meanBS[1][0][1]=0.0699037982;
1575     _meanBS[1][0][2]=1565.72963;
1576     _meanBS[1][0][3]=841.509573;
1577     _meanBM[1][0][0]=-0.00252281385;
1578     _meanBM[1][0][1]=0.00600665031;
1579     _meanBM[1][0][2]=268.761304;
1580     _meanBM[1][0][3]=46.5945865;
1581     _meanR9[1][0][0]=0.964231565;
1582     _meanR9[1][0][1]=30.1631606;
1583     _meanR9[1][0][2]=414.510458;
1584     _meanR9[1][0][3]=0;
1585    
1586     _sigmaScale[1][0][0]=0.218991853;
1587     _sigmaScale[1][0][1]=6.93889e-18;
1588     _sigmaScale[1][0][2]=0.00939222285;
1589     _sigmaScale[1][0][3]=0;
1590     _sigmaAT[1][0][0]=1.61339852;
1591     _sigmaAT[1][0][1]=0;
1592     _sigmaAT[1][0][2]=0;
1593     _sigmaAT[1][0][3]=0;
1594     _sigmaAC[1][0][0]=0.00019476922;
1595     _sigmaAC[1][0][1]=0.697650974;
1596     _sigmaAC[1][0][2]=-0.000125668382;
1597     _sigmaAC[1][0][3]=12.8659982;
1598     _sigmaAS[1][0][0]=-1.68218147e-05;
1599     _sigmaAS[1][0][1]=6.57794255;
1600     _sigmaAS[1][0][2]=1555.93015;
1601     _sigmaAS[1][0][3]=1401.542;
1602     _sigmaAM[1][0][0]=0.0570038229;
1603     _sigmaAM[1][0][1]=0.633551691;
1604     _sigmaAM[1][0][2]=9.59639e+11;
1605     _sigmaAM[1][0][3]=16.4637695;
1606     _sigmaBT[1][0][0]=-0.0591443023;
1607     _sigmaBT[1][0][1]=0;
1608     _sigmaBT[1][0][2]=0;
1609     _sigmaBT[1][0][3]=0;
1610     _sigmaBC[1][0][0]=-0.00320070019;
1611     _sigmaBC[1][0][1]=25.5502578;
1612     _sigmaBC[1][0][2]=7.49509e+12;
1613     _sigmaBC[1][0][3]=3798165.72;
1614     _sigmaBS[1][0][0]=9.63685051e-05;
1615     _sigmaBS[1][0][1]=6.91673581;
1616     _sigmaBS[1][0][2]=2447.68053;
1617     _sigmaBS[1][0][3]=1721.11327;
1618     _sigmaBM[1][0][0]=0.00148006;
1619     _sigmaBM[1][0][1]=28;
1620     _sigmaBM[1][0][2]=5400000;
1621     _sigmaBM[1][0][3]=-9000000;
1622     _sigmaR9[1][0][0]=187.987786;
1623     _sigmaR9[1][0][1]=-1.91777372e-07;
1624     _sigmaR9[1][0][2]=8.29820105e-09;
1625     _sigmaR9[1][0][3]=0;
1626    
1627     _meanScale[1][1][0]=0.331585644;
1628     _meanScale[1][1][1]=-4.97323079e-05;
1629     _meanScale[1][1][2]=0.000208912195;
1630     _meanScale[1][1][3]=-1.36032052;
1631     _meanAT[1][1][0]=-0.0640673292;
1632     _meanAT[1][1][1]=0;
1633     _meanAT[1][1][2]=0;
1634     _meanAT[1][1][3]=0;
1635     _meanAC[1][1][0]=-0.00129027954;
1636     _meanAC[1][1][1]=0.00733510902;
1637     _meanAC[1][1][2]=182.714706;
1638     _meanAC[1][1][3]=621.652554;
1639     _meanAS[1][1][0]=-0.000490574173;
1640     _meanAS[1][1][1]=0.0308208884;
1641     _meanAS[1][1][2]=385.372647;
1642     _meanAS[1][1][3]=492.313289;
1643     _meanAM[1][1][0]=-0.0064828927;
1644     _meanAM[1][1][1]=0.649443452;
1645     _meanAM[1][1][2]=0.0573092773;
1646     _meanAM[1][1][3]=0.0743069;
1647     _meanBT[1][1][0]=-0.147343956;
1648     _meanBT[1][1][1]=0;
1649     _meanBT[1][1][2]=0;
1650     _meanBT[1][1][3]=0;
1651     _meanBC[1][1][0]=-0.00503351921;
1652     _meanBC[1][1][1]=-57691.5085;
1653     _meanBC[1][1][2]=46202.9758;
1654     _meanBC[1][1][3]=118612;
1655     _meanBS[1][1][0]=-0.000793147706;
1656     _meanBS[1][1][1]=0.0238305184;
1657     _meanBS[1][1][2]=402.215233;
1658     _meanBS[1][1][3]=455.848092;
1659     _meanBM[1][1][0]=0.000434549102;
1660     _meanBM[1][1][1]=0.0443539812;
1661     _meanBM[1][1][2]=-39970930.5;
1662     _meanBM[1][1][3]=-635.815445;
1663     _meanR9[1][1][0]=-0.411370898;
1664     _meanR9[1][1][1]=1.30133082;
1665     _meanR9[1][1][2]=-0.890618718;
1666     _meanR9[1][1][3]=0;
1667    
1668     _sigmaScale[1][1][0]=1.49352299;
1669     _sigmaScale[1][1][1]=1.38778e-17;
1670     _sigmaScale[1][1][2]=0.0248352105;
1671     _sigmaScale[1][1][3]=0;
1672     _sigmaAT[1][1][0]=-1.18239629;
1673     _sigmaAT[1][1][1]=0;
1674     _sigmaAT[1][1][2]=0;
1675     _sigmaAT[1][1][3]=0;
1676     _sigmaAC[1][1][0]=0.00155030534;
1677     _sigmaAC[1][1][1]=-0.673931391;
1678     _sigmaAC[1][1][2]=134075.829;
1679     _sigmaAC[1][1][3]=-7e+09;
1680     _sigmaAS[1][1][0]=6.95848091e-05;
1681     _sigmaAS[1][1][1]=0.522471203;
1682     _sigmaAS[1][1][2]=463.305497;
1683     _sigmaAS[1][1][3]=1159.49992;
1684     _sigmaAM[1][1][0]=-0.00509006951;
1685     _sigmaAM[1][1][1]=0.945276887;
1686     _sigmaAM[1][1][2]=46.4072512;
1687     _sigmaAM[1][1][3]=7.11474e+12;
1688     _sigmaBT[1][1][0]=-1.59480683;
1689     _sigmaBT[1][1][1]=0;
1690     _sigmaBT[1][1][2]=0;
1691     _sigmaBT[1][1][3]=0;
1692     _sigmaBC[1][1][0]=-0.00202302997;
1693     _sigmaBC[1][1][1]=15.4301057;
1694     _sigmaBC[1][1][2]=-33315545.5;
1695     _sigmaBC[1][1][3]=-6e+09;
1696     _sigmaBS[1][1][0]=0.00271126099;
1697     _sigmaBS[1][1][1]=0.325669289;
1698     _sigmaBS[1][1][2]=2322.66097;
1699     _sigmaBS[1][1][3]=298.692034;
1700     _sigmaBM[1][1][0]=-0.0454765849;
1701     _sigmaBM[1][1][1]=6.81541098;
1702     _sigmaBM[1][1][2]=1.9e+09;
1703     _sigmaBM[1][1][3]=-26353.4449;
1704     _sigmaR9[1][1][0]=41.1074567;
1705     _sigmaR9[1][1][1]=-86.9595346;
1706     _sigmaR9[1][1][2]=45.7818889;
1707     _sigmaR9[1][1][3]=0;
1708    
1709     _initialised=true;
1710     }
1711    
1712     if(s=="4_2e") {
1713     _meanScale[0][0][0]=1.03294629;
1714     _meanScale[0][0][1]=-0.000210626517;
1715     _meanScale[0][0][2]=0.000268568795;
1716     _meanScale[0][0][3]=0.338053561;
1717     _meanAT[0][0][0]=0.0200811135;
1718     _meanAT[0][0][1]=0;
1719     _meanAT[0][0][2]=0;
1720     _meanAT[0][0][3]=0;
1721     _meanAC[0][0][0]=-0.00326696352;
1722     _meanAC[0][0][1]=0.010765809;
1723     _meanAC[0][0][2]=513.763513;
1724     _meanAC[0][0][3]=546.438243;
1725     _meanAS[0][0][0]=0;
1726     _meanAS[0][0][1]=0;
1727     _meanAS[0][0][2]=0;
1728     _meanAS[0][0][3]=0;
1729     _meanAM[0][0][0]=-0.00135522301;
1730     _meanAM[0][0][1]=0.166490439;
1731     _meanAM[0][0][2]=278.324187;
1732     _meanAM[0][0][3]=245.998361;
1733     _meanBT[0][0][0]=0;
1734     _meanBT[0][0][1]=0;
1735     _meanBT[0][0][2]=0;
1736     _meanBT[0][0][3]=0;
1737     _meanBC[0][0][0]=-0.00332906015;
1738     _meanBC[0][0][1]=0.00792585358;
1739     _meanBC[0][0][2]=514.766605;
1740     _meanBC[0][0][3]=488.870257;
1741     _meanBS[0][0][0]=-0.00199241828;
1742     _meanBS[0][0][1]=0.0037942702;
1743     _meanBS[0][0][2]=29.9438726;
1744     _meanBS[0][0][3]=1077.1644;
1745     _meanBM[0][0][0]=-0.00159080193;
1746     _meanBM[0][0][1]=0.107998922;
1747     _meanBM[0][0][2]=229.934523;
1748     _meanBM[0][0][3]=231.786153;
1749     _meanR9[0][0][0]=0.857844414;
1750     _meanR9[0][0][1]=-16.8494499;
1751     _meanR9[0][0][2]=125.493331;
1752     _meanR9[0][0][3]=0;
1753    
1754     _sigmaScale[0][0][0]=0.392737806;
1755     _sigmaScale[0][0][1]=0.0353140568;
1756     _sigmaScale[0][0][2]=-0.00613223131;
1757     _sigmaScale[0][0][3]=0;
1758     _sigmaAT[0][0][0]=1.02977565;
1759     _sigmaAT[0][0][1]=0;
1760     _sigmaAT[0][0][2]=0;
1761     _sigmaAT[0][0][3]=0;
1762     _sigmaAC[0][0][0]=-0.00350109526;
1763     _sigmaAC[0][0][1]=-0.951103069;
1764     _sigmaAC[0][0][2]=-54434.4267;
1765     _sigmaAC[0][0][3]=-3e+17;
1766     _sigmaAS[0][0][0]=0;
1767     _sigmaAS[0][0][1]=0;
1768     _sigmaAS[0][0][2]=0;
1769     _sigmaAS[0][0][3]=0;
1770     _sigmaAM[0][0][0]=0.00127749544;
1771     _sigmaAM[0][0][1]=5.03867192;
1772     _sigmaAM[0][0][2]=563.047721;
1773     _sigmaAM[0][0][3]=272.293234;
1774     _sigmaBT[0][0][0]=0.00480679465;
1775     _sigmaBT[0][0][1]=7.56230742;
1776     _sigmaBT[0][0][2]=-33600000;
1777     _sigmaBT[0][0][3]=-257.677353;
1778     _sigmaBC[0][0][0]=-0.00169935002;
1779     _sigmaBC[0][0][1]=2790083.26;
1780     _sigmaBC[0][0][2]=-97275416.4;
1781     _sigmaBC[0][0][3]=23710676.7;
1782     _sigmaBS[0][0][0]=0;
1783     _sigmaBS[0][0][1]=0;
1784     _sigmaBS[0][0][2]=0;
1785     _sigmaBS[0][0][3]=0;
1786     _sigmaBM[0][0][0]=-0.00194553738;
1787     _sigmaBM[0][0][1]=7.77713222;
1788     _sigmaBM[0][0][2]=264.960159;
1789     _sigmaBM[0][0][3]=363.487107;
1790     _sigmaR9[0][0][0]=0.952571;
1791     _sigmaR9[0][0][1]=0;
1792     _sigmaR9[0][0][2]=0;
1793     _sigmaR9[0][0][3]=0;
1794    
1795     _meanScale[0][1][0]=0.86164193;
1796     _meanScale[0][1][1]=-0.0001184458;
1797     _meanScale[0][1][2]=0.000232979403;
1798     _meanScale[0][1][3]=0.310305987;
1799     _meanAT[0][1][0]=0.0103409006;
1800     _meanAT[0][1][1]=0;
1801     _meanAT[0][1][2]=0;
1802     _meanAT[0][1][3]=0;
1803     _meanAC[0][1][0]=-0.00325081301;
1804     _meanAC[0][1][1]=0.0208748426;
1805     _meanAC[0][1][2]=165.245698;
1806     _meanAC[0][1][3]=292.03632;
1807     _meanAS[0][1][0]=0.0330004;
1808     _meanAS[0][1][1]=-148569.764;
1809     _meanAS[0][1][2]=87999432.1;
1810     _meanAS[0][1][3]=7787218.96;
1811     _meanAM[0][1][0]=-0.000867413605;
1812     _meanAM[0][1][1]=0.10580464;
1813     _meanAM[0][1][2]=396.92529;
1814     _meanAM[0][1][3]=263.112883;
1815     _meanBT[0][1][0]=0;
1816     _meanBT[0][1][1]=0.216283067;
1817     _meanBT[0][1][2]=312.543466;
1818     _meanBT[0][1][3]=463.601293;
1819     _meanBC[0][1][0]=-0.00505883024;
1820     _meanBC[0][1][1]=0.00182528255;
1821     _meanBC[0][1][2]=507.478054;
1822     _meanBC[0][1][3]=-6837.26736;
1823     _meanBS[0][1][0]=-166707004;
1824     _meanBS[0][1][1]=0.0928055999;
1825     _meanBS[0][1][2]=-5.30004162e-11;
1826     _meanBS[0][1][3]=11442.2;
1827     _meanBM[0][1][0]=-5.93998135e-05;
1828     _meanBM[0][1][1]=0.0096852184;
1829     _meanBM[0][1][2]=59.8040186;
1830     _meanBM[0][1][3]=-440000000;
1831     _meanR9[0][1][0]=0.0716647946;
1832     _meanR9[0][1][1]=-0.204241803;
1833     _meanR9[0][1][2]=0.154962477;
1834     _meanR9[0][1][3]=0;
1835    
1836     _sigmaScale[0][1][0]=0.469123815;
1837     _sigmaScale[0][1][1]=-0.090283052;
1838     _sigmaScale[0][1][2]=0.000469934719;
1839     _sigmaScale[0][1][3]=0;
1840     _sigmaAT[0][1][0]=1.77629522;
1841     _sigmaAT[0][1][1]=0;
1842     _sigmaAT[0][1][2]=0;
1843     _sigmaAT[0][1][3]=0;
1844     _sigmaAC[0][1][0]=-0.00636220086;
1845     _sigmaAC[0][1][1]=-0.781271127;
1846     _sigmaAC[0][1][2]=4.90734224;
1847     _sigmaAC[0][1][3]=65.6835127;
1848     _sigmaAS[0][1][0]=0;
1849     _sigmaAS[0][1][1]=0;
1850     _sigmaAS[0][1][2]=0;
1851     _sigmaAS[0][1][3]=0;
1852     _sigmaAM[0][1][0]=0.000179292631;
1853     _sigmaAM[0][1][1]=7.62815501;
1854     _sigmaAM[0][1][2]=743.55507;
1855     _sigmaAM[0][1][3]=354.656661;
1856     _sigmaBT[0][1][0]=-0.0507778073;
1857     _sigmaBT[0][1][1]=3.00903133;
1858     _sigmaBT[0][1][2]=-0.526032834;
1859     _sigmaBT[0][1][3]=-0.630748789;
1860     _sigmaBC[0][1][0]=0.00490009575;
1861     _sigmaBC[0][1][1]=-1.53772346;
1862     _sigmaBC[0][1][2]=553415.545;
1863     _sigmaBC[0][1][3]=2.36808e+19;
1864     _sigmaBS[0][1][0]=0;
1865     _sigmaBS[0][1][1]=0;
1866     _sigmaBS[0][1][2]=0;
1867     _sigmaBS[0][1][3]=0;
1868     _sigmaBM[0][1][0]=-0.00113947453;
1869     _sigmaBM[0][1][1]=3.74348887;
1870     _sigmaBM[0][1][2]=91.9478901;
1871     _sigmaBM[0][1][3]=101.304882;
1872     _sigmaR9[0][1][0]=-0.261512815;
1873     _sigmaR9[0][1][1]=-1.69974425;
1874     _sigmaR9[0][1][2]=0;
1875     _sigmaR9[0][1][3]=0;
1876    
1877     _meanScale[1][0][0]=0.961072344;
1878     _meanScale[1][0][1]=8.81367775e-05;
1879     _meanScale[1][0][2]=-0.000270690177;
1880     _meanScale[1][0][3]=0.745461418;
1881     _meanAT[1][0][0]=0.532495533;
1882     _meanAT[1][0][1]=0;
1883     _meanAT[1][0][2]=0;
1884     _meanAT[1][0][3]=0;
1885     _meanAC[1][0][0]=-0.000539999855;
1886     _meanAC[1][0][1]=0.0100918811;
1887     _meanAC[1][0][2]=953.905309;
1888     _meanAC[1][0][3]=808.944612;
1889     _meanAS[1][0][0]=-0.000597157153;
1890     _meanAS[1][0][1]=0.0571921693;
1891     _meanAS[1][0][2]=700.692431;
1892     _meanAS[1][0][3]=924.653733;
1893     _meanAM[1][0][0]=0.000230736156;
1894     _meanAM[1][0][1]=1.77368196;
1895     _meanAM[1][0][2]=4461.03178;
1896     _meanAM[1][0][3]=3300.73792;
1897     _meanBT[1][0][0]=0.483274186;
1898     _meanBT[1][0][1]=0;
1899     _meanBT[1][0][2]=0;
1900     _meanBT[1][0][3]=0;
1901     _meanBC[1][0][0]=-0.000651403853;
1902     _meanBC[1][0][1]=0.0111101805;
1903     _meanBC[1][0][2]=1276.07724;
1904     _meanBC[1][0][3]=1489.51887;
1905     _meanBS[1][0][0]=-0.000251246189;
1906     _meanBS[1][0][1]=0.0530409004;
1907     _meanBS[1][0][2]=767.699586;
1908     _meanBS[1][0][3]=835.195311;
1909     _meanBM[1][0][0]=-0.187856578;
1910     _meanBM[1][0][1]=-0.00821848896;
1911     _meanBM[1][0][2]=0.891813494;
1912     _meanBM[1][0][3]=-580000000;
1913     _meanR9[1][0][0]=0.96358076;
1914     _meanR9[1][0][1]=28.7116938;
1915     _meanR9[1][0][2]=697.709731;
1916     _meanR9[1][0][3]=0;
1917    
1918     _sigmaScale[1][0][0]=0.46256953;
1919     _sigmaScale[1][0][1]=-2.50963561e-08;
1920     _sigmaScale[1][0][2]=0.0139636379;
1921     _sigmaScale[1][0][3]=0;
1922     _sigmaAT[1][0][0]=6.47165025;
1923     _sigmaAT[1][0][1]=0;
1924     _sigmaAT[1][0][2]=0;
1925     _sigmaAT[1][0][3]=0;
1926     _sigmaAC[1][0][0]=48.1275;
1927     _sigmaAC[1][0][1]=150005000;
1928     _sigmaAC[1][0][2]=21231.6;
1929     _sigmaAC[1][0][3]=2.6e+11;
1930     _sigmaAS[1][0][0]=0.000209127817;
1931     _sigmaAS[1][0][1]=2.19868731;
1932     _sigmaAS[1][0][2]=1695.98579;
1933     _sigmaAS[1][0][3]=967.250228;
1934     _sigmaAM[1][0][0]=0.0217972665;
1935     _sigmaAM[1][0][1]=1.26317651;
1936     _sigmaAM[1][0][2]=34.0924905;
1937     _sigmaAM[1][0][3]=55.1895282;
1938     _sigmaBT[1][0][0]=5.21983754;
1939     _sigmaBT[1][0][1]=0;
1940     _sigmaBT[1][0][2]=0;
1941     _sigmaBT[1][0][3]=0;
1942     _sigmaBC[1][0][0]=-0.004;
1943     _sigmaBC[1][0][1]=-120000;
1944     _sigmaBC[1][0][2]=7.49509e+12;
1945     _sigmaBC[1][0][3]=36643600;
1946     _sigmaBS[1][0][0]=0.000250338051;
1947     _sigmaBS[1][0][1]=1.98819262;
1948     _sigmaBS[1][0][2]=1967.55308;
1949     _sigmaBS[1][0][3]=1098.23855;
1950     _sigmaBM[1][0][0]=0.00101799874;
1951     _sigmaBM[1][0][1]=88.0546723;
1952     _sigmaBM[1][0][2]=8.47552e+10;
1953     _sigmaBM[1][0][3]=-132255.757;
1954     _sigmaR9[1][0][0]=144.031062;
1955     _sigmaR9[1][0][1]=-6.11507616e-07;
1956     _sigmaR9[1][0][2]=1.18181734e-08;
1957     _sigmaR9[1][0][3]=0;
1958    
1959     _meanScale[1][1][0]=0.288888347;
1960     _meanScale[1][1][1]=6.52038486e-06;
1961     _meanScale[1][1][2]=0.000173654897;
1962     _meanScale[1][1][3]=0.422671325;
1963     _meanAT[1][1][0]=0.0614964598;
1964     _meanAT[1][1][1]=0;
1965     _meanAT[1][1][2]=0;
1966     _meanAT[1][1][3]=0;
1967     _meanAC[1][1][0]=-0.00123181641;
1968     _meanAC[1][1][1]=0.0133568947;
1969     _meanAC[1][1][2]=165.847556;
1970     _meanAC[1][1][3]=332.705784;
1971     _meanAS[1][1][0]=-0.00088161986;
1972     _meanAS[1][1][1]=0.0304986746;
1973     _meanAS[1][1][2]=382.755876;
1974     _meanAS[1][1][3]=616.470187;
1975     _meanAM[1][1][0]=0.000980695422;
1976     _meanAM[1][1][1]=0.63575757;
1977     _meanAM[1][1][2]=0.0336097848;
1978     _meanAM[1][1][3]=0.043315868;
1979     _meanBT[1][1][0]=0.11623414;
1980     _meanBT[1][1][1]=0;
1981     _meanBT[1][1][2]=0;
1982     _meanBT[1][1][3]=0;
1983     _meanBC[1][1][0]=-0.00716072255;
1984     _meanBC[1][1][1]=-0.440696266;
1985     _meanBC[1][1][2]=1887.74154;
1986     _meanBC[1][1][3]=118612;
1987     _meanBS[1][1][0]=-0.000492035977;
1988     _meanBS[1][1][1]=0.0292167014;
1989     _meanBS[1][1][2]=433.232787;
1990     _meanBS[1][1][3]=484.310448;
1991     _meanBM[1][1][0]=0.00299476541;
1992     _meanBM[1][1][1]=0.0149328977;
1993     _meanBM[1][1][2]=-48728700;
1994     _meanBM[1][1][3]=37.0041547;
1995     _meanR9[1][1][0]=0.19617696;
1996     _meanR9[1][1][1]=-0.350976375;
1997     _meanR9[1][1][2]=0.181094838;
1998     _meanR9[1][1][3]=0;
1999    
2000     _sigmaScale[1][1][0]=1.26164895;
2001     _sigmaScale[1][1][1]=-6.61150347e-07;
2002     _sigmaScale[1][1][2]=0.0280532297;
2003     _sigmaScale[1][1][3]=0;
2004     _sigmaAT[1][1][0]=-0.232612761;
2005     _sigmaAT[1][1][1]=0;
2006     _sigmaAT[1][1][2]=0;
2007     _sigmaAT[1][1][3]=0;
2008     _sigmaAC[1][1][0]=0.00137406444;
2009     _sigmaAC[1][1][1]=-0.377659364;
2010     _sigmaAC[1][1][2]=27171.5802;
2011     _sigmaAC[1][1][3]=-560000000;
2012     _sigmaAS[1][1][0]=0.00022943714;
2013     _sigmaAS[1][1][1]=0.335082568;
2014     _sigmaAS[1][1][2]=590.511812;
2015     _sigmaAS[1][1][3]=387.352521;
2016     _sigmaAM[1][1][0]=-0.000780390674;
2017     _sigmaAM[1][1][1]=1.05127796;
2018     _sigmaAM[1][1][2]=33.7378914;
2019     _sigmaAM[1][1][3]=61.3730807;
2020     _sigmaBT[1][1][0]=0.529507693;
2021     _sigmaBT[1][1][1]=0;
2022     _sigmaBT[1][1][2]=0;
2023     _sigmaBT[1][1][3]=0;
2024     _sigmaBC[1][1][0]=-0.00203996;
2025     _sigmaBC[1][1][1]=93000;
2026     _sigmaBC[1][1][2]=61225800;
2027     _sigmaBC[1][1][3]=-4.43323e+17;
2028     _sigmaBS[1][1][0]=0.00125939613;
2029     _sigmaBS[1][1][1]=0.31048111;
2030     _sigmaBS[1][1][2]=295.258764;
2031     _sigmaBS[1][1][3]=263.974257;
2032     _sigmaBM[1][1][0]=-0.046100748;
2033     _sigmaBM[1][1][1]=1.22348596;
2034     _sigmaBM[1][1][2]=1.9e+09;
2035     _sigmaBM[1][1][3]=1254.99;
2036     _sigmaR9[1][1][0]=9.09347838;
2037     _sigmaR9[1][1][1]=-10.0390435;
2038     _sigmaR9[1][1][2]=0;
2039     _sigmaR9[1][1][3]=0;
2040    
2041     _initialised=true;
2042     }
2043    
2044     assert(_initialised);
2045     return true;
2046     }
2047    
2048     // Get the geometry of cracks and gaps from file
2049     bool PhotonFix::initialiseGeometry(const std::string &s,const std::string &infile) {
2050    
2051     std::ifstream fin(infile.c_str());
2052     assert(fin);
2053    
2054     std::cout << "Reading in here" << std::endl;
2055     for(unsigned i(0);i<169;i++) {
2056     for(unsigned j(0);j<360;j++) {
2057     for(unsigned k(0);k<2;k++) {
2058     fin >> _barrelCGap[i][j][k];
2059     }
2060     }
2061     }
2062    
2063     for(unsigned i(0);i<33;i++) {
2064     for(unsigned j(0);j<180;j++) {
2065     for(unsigned k(0);k<2;k++) {
2066     fin >> _barrelSGap[i][j][k];
2067     }
2068     }
2069     }
2070    
2071     for(unsigned i(0);i<7;i++) {
2072     for(unsigned j(0);j<18;j++) {
2073     for(unsigned k(0);k<2;k++) {
2074     fin >> _barrelMGap[i][j][k];
2075     }
2076     }
2077     }
2078     for(unsigned i(0);i<100;i++) {
2079     for(unsigned j(0);j<100;j++) {
2080     unsigned k;
2081     fin >> k;
2082     _endcapCrystal[i][j]=(k==0);
2083     }
2084     }
2085    
2086     for(unsigned i(0);i<2;i++) {
2087     for(unsigned j(0);j<7080;j++) {
2088     for(unsigned k(0);k<2;k++) {
2089     fin >> _endcapCGap[i][j][k];
2090     }
2091     }
2092     }
2093    
2094     for(unsigned i(0);i<2;i++) {
2095     for(unsigned j(0);j<264;j++) {
2096     for(unsigned k(0);k<2;k++) {
2097     fin >> _endcapSGap[i][j][k];
2098     }
2099     }
2100     }
2101    
2102     for(unsigned i(0);i<2;i++) {
2103     for(unsigned j(0);j<1;j++) {
2104     for(unsigned k(0);k<2;k++) {
2105     fin >> _endcapMGap[i][j][k];
2106     }
2107     }
2108     }
2109    
2110     assert(fin);
2111    
2112     return true;
2113     }
2114    
2115    
2116     // bool PhotonFix::_initialised=false;
2117     //
2118     // double PhotonFix::_meanScale[2][2][4];
2119     // double PhotonFix::_meanAT[2][2][4];
2120     // double PhotonFix::_meanAC[2][2][4];
2121     // double PhotonFix::_meanAS[2][2][4];
2122     // double PhotonFix::_meanAM[2][2][4];
2123     // double PhotonFix::_meanBT[2][2][4];
2124     // double PhotonFix::_meanBC[2][2][4];
2125     // double PhotonFix::_meanBS[2][2][4];
2126     // double PhotonFix::_meanBM[2][2][4];
2127     // double PhotonFix::_meanR9[2][2][4];
2128     //
2129     // double PhotonFix::_sigmaScale[2][2][4];
2130     // double PhotonFix::_sigmaAT[2][2][4];
2131     // double PhotonFix::_sigmaAC[2][2][4];
2132     // double PhotonFix::_sigmaAS[2][2][4];
2133     // double PhotonFix::_sigmaAM[2][2][4];
2134     // double PhotonFix::_sigmaBT[2][2][4];
2135     // double PhotonFix::_sigmaBC[2][2][4];
2136     // double PhotonFix::_sigmaBS[2][2][4];
2137     // double PhotonFix::_sigmaBM[2][2][4];
2138     // double PhotonFix::_sigmaR9[2][2][4];
2139     //
2140     // double PhotonFix::_barrelCGap[169][360][2];
2141     // double PhotonFix::_barrelSGap[33][180][2];
2142     // double PhotonFix::_barrelMGap[7][18][2];
2143     //
2144     // bool PhotonFix::_endcapCrystal[100][100];
2145     // double PhotonFix::_endcapCGap[2][7080][2];
2146     // double PhotonFix::_endcapSGap[2][264][2];
2147     // double PhotonFix::_endcapMGap[2][1][2];