Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/inclxx/incl_physics/src/G4INCLNKbElasticChannel.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/models/inclxx/incl_physics/src/G4INCLNKbElasticChannel.cc (Version 11.3.0) and /processes/hadronic/models/inclxx/incl_physics/src/G4INCLNKbElasticChannel.cc (Version 11.0.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 // INCL++ intra-nuclear cascade model              26 // INCL++ intra-nuclear cascade model
 27 // Alain Boudard, CEA-Saclay, France               27 // Alain Boudard, CEA-Saclay, France
 28 // Joseph Cugnon, University of Liege, Belgium     28 // Joseph Cugnon, University of Liege, Belgium
 29 // Jean-Christophe David, CEA-Saclay, France       29 // Jean-Christophe David, CEA-Saclay, France
 30 // Pekka Kaitaniemi, CEA-Saclay, France, and H     30 // Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
 31 // Sylvie Leray, CEA-Saclay, France                31 // Sylvie Leray, CEA-Saclay, France
 32 // Davide Mancusi, CEA-Saclay, France              32 // Davide Mancusi, CEA-Saclay, France
 33 //                                                 33 //
 34 #define INCLXX_IN_GEANT4_MODE 1                    34 #define INCLXX_IN_GEANT4_MODE 1
 35                                                    35 
 36 #include "globals.hh"                              36 #include "globals.hh"
 37                                                    37 
 38 #include "G4INCLNKbElasticChannel.hh"              38 #include "G4INCLNKbElasticChannel.hh"
 39 #include "G4INCLKinematicsUtils.hh"                39 #include "G4INCLKinematicsUtils.hh"
 40 #include "G4INCLBinaryCollisionAvatar.hh"          40 #include "G4INCLBinaryCollisionAvatar.hh"
 41 #include "G4INCLRandom.hh"                         41 #include "G4INCLRandom.hh"
 42 #include "G4INCLGlobals.hh"                        42 #include "G4INCLGlobals.hh"
 43 #include "G4INCLLogger.hh"                         43 #include "G4INCLLogger.hh"
 44 #include <algorithm>                               44 #include <algorithm>
 45 #include "G4INCLPhaseSpaceGenerator.hh"            45 #include "G4INCLPhaseSpaceGenerator.hh"
 46                                                    46 
 47 namespace G4INCL {                                 47 namespace G4INCL {
 48                                                    48   
 49   NKbElasticChannel::NKbElasticChannel(Particl     49   NKbElasticChannel::NKbElasticChannel(Particle *p1, Particle *p2)
 50     : particle1(p1), particle2(p2)                 50     : particle1(p1), particle2(p2)
 51     {}                                             51     {}
 52                                                    52   
 53   NKbElasticChannel::~NKbElasticChannel(){}        53   NKbElasticChannel::~NKbElasticChannel(){}
 54                                                    54   
 55   void NKbElasticChannel::fillFinalState(Final     55   void NKbElasticChannel::fillFinalState(FinalState *fs){
 56                                                    56     
 57     Particle *nucleon;                             57     Particle *nucleon;
 58     Particle *kaon;                                58     Particle *kaon;
 59                                                    59     
 60     if(particle1->isNucleon()){                    60     if(particle1->isNucleon()){
 61       nucleon = particle1;                         61       nucleon = particle1;
 62       kaon = particle2;                            62       kaon = particle2;
 63     }                                              63     }
 64     else{                                          64     else{
 65       nucleon = particle2;                         65       nucleon = particle2;
 66       kaon = particle1;                            66       kaon = particle1;
 67     }                                              67     }
 68                                                    68     
 69     ThreeVector mom_kaon = KaonMomentum(kaon,n     69     ThreeVector mom_kaon = KaonMomentum(kaon,nucleon); // Caution: angular distribution comes from K-p which is "Lambda resonant"
 70                                                    70 
 71     kaon->setMomentum(mom_kaon);                   71     kaon->setMomentum(mom_kaon);
 72     nucleon->setMomentum(-mom_kaon);               72     nucleon->setMomentum(-mom_kaon);
 73                                                    73     
 74     nucleon->adjustEnergyFromMomentum();           74     nucleon->adjustEnergyFromMomentum();
 75     kaon->adjustEnergyFromMomentum();              75     kaon->adjustEnergyFromMomentum();
 76                                                    76     
 77     fs->addModifiedParticle(nucleon);              77     fs->addModifiedParticle(nucleon);
 78     fs->addModifiedParticle(kaon);                 78     fs->addModifiedParticle(kaon);
 79                                                    79         
 80   }                                                80   }
 81                                                    81   
 82   ThreeVector NKbElasticChannel::KaonMomentum(     82   ThreeVector NKbElasticChannel::KaonMomentum(Particle const * const kaon, Particle const * const nucleon){
 83                                                    83     
 84     const G4double norm = KinematicsUtils::mom     84     const G4double norm = KinematicsUtils::momentumInCM(kaon,nucleon);
 85     const G4double pLab = KinematicsUtils::mom     85     const G4double pLab = KinematicsUtils::momentumInLab(kaon,nucleon);
 86                                                    86     
 87     if(pLab < 225.) return Random::normVector(     87     if(pLab < 225.) return Random::normVector()*norm; // isotropic
 88                                                    88     
 89     G4double cos_theta = 1.;                       89     G4double cos_theta = 1.;
 90     G4double sin_theta = 0.;                       90     G4double sin_theta = 0.;
 91     const G4double cos_phi = std::cos(Random::     91     const G4double cos_phi = std::cos(Random::shoot()*Math::twoPi);
 92     const G4double sin_phi = std::sqrt(1-cos_p     92     const G4double sin_phi = std::sqrt(1-cos_phi*cos_phi);
 93                                                    93     
 94     const G4double x = kaon->getMomentum().get     94     const G4double x = kaon->getMomentum().getX();
 95     const G4double y = kaon->getMomentum().get     95     const G4double y = kaon->getMomentum().getY();
 96     const G4double z = kaon->getMomentum().get     96     const G4double z = kaon->getMomentum().getZ();
 97                                                    97     
 98     const G4double r = std::sqrt(x*x+y*y+z*z);     98     const G4double r = std::sqrt(x*x+y*y+z*z);
 99     const G4double rho = std::sqrt(x*x+y*y);       99     const G4double rho = std::sqrt(x*x+y*y);
100                                                   100     
101     if(pLab >= 2375.){                            101     if(pLab >= 2375.){
102       const G4double b = 12. * pLab/2375.; //     102       const G4double b = 12. * pLab/2375.; // correspond to the forward slope description at 2375 MeV/c
103       cos_theta = std::log(Random::shoot()*(st    103       cos_theta = std::log(Random::shoot()*(std::exp(b)-std::exp(-b))+std::exp(-b))/b;
104       sin_theta = std::sqrt(1-cos_theta*cos_th    104       sin_theta = std::sqrt(1-cos_theta*cos_theta);
105                                                   105       
106     }                                             106     }
107     else{                                         107     else{
108       const G4double Legendre_coef[431][11] =     108       const G4double Legendre_coef[431][11] = {
109         {225,0.14007,-0.1329,0.00044,-0.00469,    109         {225,0.14007,-0.1329,0.00044,-0.00469,-0.00254,-0.00011,0,0,0,0},
110         {230,0.14143,-0.1335,0.00058,-0.00417,    110         {230,0.14143,-0.1335,0.00058,-0.00417,-0.00248,-4E-05,0,0,0,0},
111         {235,0.14283,-0.1179,0.00077,-0.00364,    111         {235,0.14283,-0.1179,0.00077,-0.00364,-0.00243,2E-05,0,0,0,0},
112         {240,0.14428,-0.08345,0.00102,-0.00312    112         {240,0.14428,-0.08345,0.00102,-0.00312,-0.00238,9E-05,0,0,0,0},
113         {245,0.14582,-0.04111,0.00133,-0.00259    113         {245,0.14582,-0.04111,0.00133,-0.00259,-0.00232,0.00015,0,0,0,0},
114         {250,0.14746,-0.00771,0.00173,-0.00206    114         {250,0.14746,-0.00771,0.00173,-0.00206,-0.00226,0.00022,0,0,0,0},
115         {255,0.14924,0.00768,0.00222,-0.00153,    115         {255,0.14924,0.00768,0.00222,-0.00153,-0.0022,0.00028,0,0,0,0},
116         {260,0.15117,0.00618,0.00283,-0.00098,    116         {260,0.15117,0.00618,0.00283,-0.00098,-0.00214,0.00035,0,0,0,0},
117         {265,0.15327,-0.00651,0.00358,-0.00043    117         {265,0.15327,-0.00651,0.00358,-0.00043,-0.00206,0.00041,0,0,0,0},
118         {270,0.15552,-0.02353,0.00449,0.00014,    118         {270,0.15552,-0.02353,0.00449,0.00014,-0.00198,0.00048,0,0,0,0},
119         {275,0.15793,-0.03863,0.00558,0.00072,    119         {275,0.15793,-0.03863,0.00558,0.00072,-0.0019,0.00054,0,0,0,0},
120         {280,0.16051,-0.04941,0.00685,0.00131,    120         {280,0.16051,-0.04941,0.00685,0.00131,-0.0018,0.0006,0,0,0,0},
121         {285,0.16325,-0.05717,0.00834,0.00192,    121         {285,0.16325,-0.05717,0.00834,0.00192,-0.00168,0.00067,0,0,0,0},
122         {290,0.16617,-0.06091,0.01006,0.00255,    122         {290,0.16617,-0.06091,0.01006,0.00255,-0.00155,0.00073,0,0,0,0},
123         {295,0.16926,-0.05462,0.01201,0.0032,-    123         {295,0.16926,-0.05462,0.01201,0.0032,-0.00141,0.00079,0,0,0,0},
124         {300,0.17254,-0.03426,0.0142,0.00387,-    124         {300,0.17254,-0.03426,0.0142,0.00387,-0.00124,0.00085,0,0,0,0},
125         {305,0.17599,-0.00416,0.01663,0.00457,    125         {305,0.17599,-0.00416,0.01663,0.00457,-0.00105,0.00091,0,0,0,0},
126         {310,0.17961,0.02807,0.01931,0.00529,-    126         {310,0.17961,0.02807,0.01931,0.00529,-0.00084,0.00097,0,0,0,0},
127         {315,0.18341,0.05924,0.02221,0.00603,-    127         {315,0.18341,0.05924,0.02221,0.00603,-0.0006,0.00102,0,0,0,0},
128         {320,0.18737,0.09129,0.02534,0.00681,-    128         {320,0.18737,0.09129,0.02534,0.00681,-0.00033,0.00107,0,0,0,0},
129         {325,0.19147,0.12688,0.02868,0.0076,-2    129         {325,0.19147,0.12688,0.02868,0.0076,-2E-05,0.00112,0,0,0,0},
130         {330,0.19568,0.16704,0.03221,0.00841,0    130         {330,0.19568,0.16704,0.03221,0.00841,0.00032,0.00116,0,0,0,0},
131         {335,0.20001,0.21052,0.03591,0.00927,0    131         {335,0.20001,0.21052,0.03591,0.00927,0.0007,0.0012,0,0,0,0},
132         {340,0.20444,0.25417,0.03977,0.01019,0    132         {340,0.20444,0.25417,0.03977,0.01019,0.00113,0.00123,0,0,0,0},
133         {345,0.20897,0.29637,0.04374,0.01121,0    133         {345,0.20897,0.29637,0.04374,0.01121,0.0016,0.00125,0,0,0,0},
134         {350,0.21361,0.34404,0.04781,0.01235,0    134         {350,0.21361,0.34404,0.04781,0.01235,0.00211,0.00126,0,0,0,0},
135         {355,0.21836,0.41178,0.05192,0.01363,0    135         {355,0.21836,0.41178,0.05192,0.01363,0.00268,0.00126,0,0,0,0},
136         {360,0.22323,0.51001,0.05605,0.01505,0    136         {360,0.22323,0.51001,0.05605,0.01505,0.00329,0.00125,0,0,0,0},
137         {365,0.22824,0.63548,0.06014,0.01663,0    137         {365,0.22824,0.63548,0.06014,0.01663,0.00396,0.00123,0,0,0,0},
138         {370,0.2334,0.7775,0.06415,0.01839,0.0    138         {370,0.2334,0.7775,0.06415,0.01839,0.00469,0.0012,0,0,0,0},
139         {375,0.23875,0.92674,0.06804,0.02031,0    139         {375,0.23875,0.92674,0.06804,0.02031,0.00548,0.00114,0,0,0,0},
140         {380,0.24435,1.07258,0.07178,0.02238,0    140         {380,0.24435,1.07258,0.07178,0.02238,0.00632,0.00107,0,0,0,0},
141         {385,0.25024,1.19622,0.07534,0.0246,0.    141         {385,0.25024,1.19622,0.07534,0.0246,0.00723,0.00099,1E-05,0,0,0},
142         {390,0.25648,1.27148,0.07871,0.02695,0    142         {390,0.25648,1.27148,0.07871,0.02695,0.00819,0.00088,1E-05,0,0,0},
143         {395,0.26314,1.27773,0.08191,0.02941,0    143         {395,0.26314,1.27773,0.08191,0.02941,0.00922,0.00074,1E-05,0,0,0},
144         {400,0.2703,1.21317,0.08493,0.03197,0.    144         {400,0.2703,1.21317,0.08493,0.03197,0.0103,0.00058,1E-05,0,0,0},
145         {405,0.27805,1.09842,0.08781,0.03461,0    145         {405,0.27805,1.09842,0.08781,0.03461,0.01144,0.0004,1E-05,0,0,0},
146         {410,0.28645,0.96106,0.09058,0.03732,0    146         {410,0.28645,0.96106,0.09058,0.03732,0.01264,0.00018,1E-05,0,0,0},
147         {415,0.29558,0.82322,0.09326,0.04009,0    147         {415,0.29558,0.82322,0.09326,0.04009,0.01389,-7E-05,1E-05,0,0,0},
148         {420,0.30549,0.70361,0.09589,0.04293,0    148         {420,0.30549,0.70361,0.09589,0.04293,0.01518,-0.00035,2E-05,0,0,0},
149         {425,0.31624,0.61949,0.09849,0.04579,0    149         {425,0.31624,0.61949,0.09849,0.04579,0.01652,-0.00067,2E-05,0,0,0},
150         {430,0.32784,0.57051,0.10106,0.04865,0    150         {430,0.32784,0.57051,0.10106,0.04865,0.0179,-0.00102,2E-05,0,0,0},
151         {435,0.34034,0.53542,0.10361,0.05148,0    151         {435,0.34034,0.53542,0.10361,0.05148,0.0193,-0.00142,2E-05,0,0,0},
152         {440,0.35375,0.49078,0.10615,0.05423,0    152         {440,0.35375,0.49078,0.10615,0.05423,0.02072,-0.00186,2E-05,0,0,0},
153         {445,0.36809,0.42743,0.10865,0.05688,0    153         {445,0.36809,0.42743,0.10865,0.05688,0.02216,-0.00234,2E-05,0,0,0},
154         {450,0.38335,0.35685,0.11111,0.05938,0    154         {450,0.38335,0.35685,0.11111,0.05938,0.02361,-0.00287,3E-05,0,0,0},
155         {455,0.39954,0.29473,0.11352,0.0617,0.    155         {455,0.39954,0.29473,0.11352,0.0617,0.02506,-0.00345,3E-05,0,0,0},
156         {460,0.41664,0.24563,0.11587,0.0638,0.    156         {460,0.41664,0.24563,0.11587,0.0638,0.02653,-0.00407,3E-05,0,0,0},
157         {465,0.43459,0.21189,0.11819,0.06569,0    157         {465,0.43459,0.21189,0.11819,0.06569,0.02798,-0.00474,3E-05,0,0,0},
158         {470,0.45332,0.19982,0.12049,0.06737,0    158         {470,0.45332,0.19982,0.12049,0.06737,0.0294,-0.00546,3E-05,0,0,0},
159         {475,0.47275,0.20898,0.12283,0.06886,0    159         {475,0.47275,0.20898,0.12283,0.06886,0.03078,-0.00623,4E-05,0,0,0},
160         {480,0.49282,0.22837,0.12528,0.07016,0    160         {480,0.49282,0.22837,0.12528,0.07016,0.03211,-0.00705,4E-05,0,0,0},
161         {485,0.51344,0.24427,0.12795,0.0713,0.    161         {485,0.51344,0.24427,0.12795,0.0713,0.03337,-0.00791,4E-05,0,0,0},
162         {490,0.53457,0.25093,0.13095,0.07231,0    162         {490,0.53457,0.25093,0.13095,0.07231,0.03456,-0.00882,4E-05,0,0,0},
163         {495,0.55613,0.25424,0.13442,0.07322,0    163         {495,0.55613,0.25424,0.13442,0.07322,0.03567,-0.00977,4E-05,0,0,0},
164         {500,0.57806,0.26108,0.13848,0.07406,0    164         {500,0.57806,0.26108,0.13848,0.07406,0.03671,-0.01076,4E-05,0,0,0},
165         {505,0.60029,0.27305,0.14325,0.07486,0    165         {505,0.60029,0.27305,0.14325,0.07486,0.03765,-0.01179,5E-05,0,0,0},
166         {510,0.62275,0.2865,0.14878,0.07566,0.    166         {510,0.62275,0.2865,0.14878,0.07566,0.03851,-0.01285,5E-05,0,0,0},
167         {515,0.64538,0.30095,0.15509,0.07648,0    167         {515,0.64538,0.30095,0.15509,0.07648,0.03926,-0.01394,5E-05,0,0,0},
168         {520,0.66812,0.32366,0.16211,0.07734,0    168         {520,0.66812,0.32366,0.16211,0.07734,0.03992,-0.01506,5E-05,0,0,0},
169         {525,0.69093,0.35704,0.16973,0.07824,0    169         {525,0.69093,0.35704,0.16973,0.07824,0.04047,-0.0162,5E-05,0,0,0},
170         {530,0.71377,0.36827,0.17776,0.07915,0    170         {530,0.71377,0.36827,0.17776,0.07915,0.04093,-0.01737,4E-05,0,0,0},
171         {535,0.7366,0.38051,0.18596,0.08008,0.    171         {535,0.7366,0.38051,0.18596,0.08008,0.0413,-0.01855,4E-05,0,0,0},
172         {540,0.75939,0.39356,0.19409,0.08101,0    172         {540,0.75939,0.39356,0.19409,0.08101,0.04158,-0.01975,4E-05,-1E-05,0,0},
173         {545,0.78211,0.40692,0.20194,0.08198,0    173         {545,0.78211,0.40692,0.20194,0.08198,0.04178,-0.02094,4E-05,-1E-05,0,0},
174         {550,0.80474,0.42044,0.20931,0.083,0.0    174         {550,0.80474,0.42044,0.20931,0.083,0.04194,-0.02212,3E-05,-1E-05,0,0},
175         {555,0.82725,0.43387,0.21611,0.08409,0    175         {555,0.82725,0.43387,0.21611,0.08409,0.04205,-0.02328,3E-05,-1E-05,0,0},
176         {560,0.84963,0.447,0.2222,0.08529,0.04    176         {560,0.84963,0.447,0.2222,0.08529,0.04216,-0.02441,2E-05,-1E-05,0,0},
177         {565,0.87184,0.45974,0.22759,0.08661,0    177         {565,0.87184,0.45974,0.22759,0.08661,0.04226,-0.0255,2E-05,-1E-05,0,0},
178         {570,0.89386,0.47209,0.23233,0.08808,0    178         {570,0.89386,0.47209,0.23233,0.08808,0.04239,-0.02654,1E-05,-1E-05,0,0},
179         {575,0.91567,0.4841,0.23658,0.08972,0.    179         {575,0.91567,0.4841,0.23658,0.08972,0.04255,-0.02752,0,-1E-05,0,0},
180         {580,0.93726,0.4957,0.24036,0.09155,0.    180         {580,0.93726,0.4957,0.24036,0.09155,0.04277,-0.02843,-1E-05,-1E-05,0,0},
181         {585,0.95862,0.50708,0.24394,0.09359,0    181         {585,0.95862,0.50708,0.24394,0.09359,0.04306,-0.02924,-2E-05,-1E-05,0,0},
182         {590,0.97974,0.51822,0.24738,0.09584,0    182         {590,0.97974,0.51822,0.24738,0.09584,0.04344,-0.02996,-3E-05,-1E-05,0,0},
183         {595,1.00062,0.52933,0.25094,0.09833,0    183         {595,1.00062,0.52933,0.25094,0.09833,0.04391,-0.03055,-4E-05,-1E-05,0,0},
184         {600,1.02125,0.54045,0.2547,0.10105,0.    184         {600,1.02125,0.54045,0.2547,0.10105,0.0445,-0.03102,-6E-05,-1E-05,0,0},
185         {605,1.04163,0.55179,0.25885,0.10403,0    185         {605,1.04163,0.55179,0.25885,0.10403,0.04522,-0.03134,-7E-05,-1E-05,0,0},
186         {610,1.06175,0.56353,0.2635,0.10727,0.    186         {610,1.06175,0.56353,0.2635,0.10727,0.04608,-0.03151,-9E-05,-1E-05,0,0},
187         {615,1.08161,0.57592,0.26878,0.1108,0.    187         {615,1.08161,0.57592,0.26878,0.1108,0.0471,-0.0315,-0.00011,-1E-05,0,0},
188         {620,1.10122,0.58905,0.27473,0.1146,0.    188         {620,1.10122,0.58905,0.27473,0.1146,0.0483,-0.03132,-0.00013,-1E-05,0,0},
189         {625,1.12057,0.6033,0.28148,0.11868,0.    189         {625,1.12057,0.6033,0.28148,0.11868,0.0497,-0.03096,-0.00015,-1E-05,0,0},
190         {630,1.13968,0.61885,0.28905,0.12305,0    190         {630,1.13968,0.61885,0.28905,0.12305,0.05131,-0.03043,-0.00017,0,0,0},
191         {635,1.15854,0.63595,0.2975,0.1277,0.0    191         {635,1.15854,0.63595,0.2975,0.1277,0.05316,-0.02972,-0.0002,0,0,0},
192         {640,1.17716,0.65493,0.30687,0.13264,0    192         {640,1.17716,0.65493,0.30687,0.13264,0.05526,-0.02883,-0.00022,0,0,0},
193         {645,1.19553,0.67602,0.31718,0.13783,0    193         {645,1.19553,0.67602,0.31718,0.13783,0.05761,-0.02778,-0.00025,0,0,0},
194         {650,1.21364,0.69939,0.32841,0.14324,0    194         {650,1.21364,0.69939,0.32841,0.14324,0.06018,-0.02657,-0.00028,0,0,0},
195         {655,1.23147,0.72528,0.34057,0.14882,0    195         {655,1.23147,0.72528,0.34057,0.14882,0.06296,-0.02522,-0.00031,0,0,0},
196         {660,1.249,0.7539,0.35367,0.15454,0.06    196         {660,1.249,0.7539,0.35367,0.15454,0.06592,-0.02375,-0.00034,0,0,0},
197         {665,1.26621,0.78533,0.36764,0.16034,0    197         {665,1.26621,0.78533,0.36764,0.16034,0.06903,-0.02217,-0.00038,0,0,0},
198         {670,1.28308,0.8195,0.3824,0.1662,0.07    198         {670,1.28308,0.8195,0.3824,0.1662,0.07229,-0.0205,-0.00041,0,0,0},
199         {675,1.29958,0.85642,0.39788,0.17209,0    199         {675,1.29958,0.85642,0.39788,0.17209,0.07568,-0.01875,-0.00045,0,1E-05,0},
200         {680,1.31567,0.89591,0.41398,0.17799,0    200         {680,1.31567,0.89591,0.41398,0.17799,0.07917,-0.01692,-0.00049,1E-05,1E-05,0},
201         {685,1.3313,0.93777,0.43055,0.1839,0.0    201         {685,1.3313,0.93777,0.43055,0.1839,0.08278,-0.01503,-0.00052,1E-05,1E-05,0},
202         {690,1.3464,0.98156,0.44745,0.1898,0.0    202         {690,1.3464,0.98156,0.44745,0.1898,0.08648,-0.01309,-0.00056,1E-05,1E-05,0},
203         {695,1.36092,1.02685,0.46446,0.19569,0    203         {695,1.36092,1.02685,0.46446,0.19569,0.09029,-0.01109,-0.0006,1E-05,1E-05,0},
204         {700,1.3748,1.07309,0.48139,0.20155,0.    204         {700,1.3748,1.07309,0.48139,0.20155,0.09418,-0.00906,-0.00064,2E-05,1E-05,0},
205         {705,1.38798,1.11974,0.49805,0.20738,0    205         {705,1.38798,1.11974,0.49805,0.20738,0.09815,-0.00698,-0.00068,2E-05,1E-05,0},
206         {710,1.4004,1.16608,0.51421,0.21317,0.    206         {710,1.4004,1.16608,0.51421,0.21317,0.10219,-0.00489,-0.00072,2E-05,1E-05,0},
207         {715,1.41199,1.2115,0.52966,0.21892,0.    207         {715,1.41199,1.2115,0.52966,0.21892,0.10629,-0.00279,-0.00076,2E-05,1E-05,0},
208         {720,1.4227,1.25536,0.54424,0.22461,0.    208         {720,1.4227,1.25536,0.54424,0.22461,0.11042,-0.00069,-0.0008,3E-05,2E-05,0},
209         {725,1.43247,1.29714,0.55781,0.23025,0    209         {725,1.43247,1.29714,0.55781,0.23025,0.11459,0.00139,-0.00083,3E-05,2E-05,0},
210         {730,1.44128,1.33634,0.57029,0.23586,0    210         {730,1.44128,1.33634,0.57029,0.23586,0.11881,0.00347,-0.00087,4E-05,2E-05,0},
211         {735,1.4491,1.37258,0.58161,0.24145,0.    211         {735,1.4491,1.37258,0.58161,0.24145,0.12307,0.00557,-0.0009,4E-05,2E-05,0},
212         {740,1.4559,1.40563,0.59178,0.24704,0.    212         {740,1.4559,1.40563,0.59178,0.24704,0.1274,0.0077,-0.00092,4E-05,2E-05,0},
213         {745,1.46166,1.43537,0.60085,0.25266,0    213         {745,1.46166,1.43537,0.60085,0.25266,0.1318,0.00989,-0.00094,5E-05,2E-05,0},
214         {750,1.4664,1.46184,0.60893,0.25834,0.    214         {750,1.4664,1.46184,0.60893,0.25834,0.1363,0.01213,-0.00096,5E-05,3E-05,0},
215         {755,1.47012,1.48507,0.61613,0.2641,0.    215         {755,1.47012,1.48507,0.61613,0.2641,0.14092,0.01446,-0.00097,6E-05,3E-05,0},
216         {760,1.47283,1.50532,0.62261,0.27,0.14    216         {760,1.47283,1.50532,0.62261,0.27,0.14569,0.01687,-0.00098,6E-05,3E-05,0},
217         {765,1.47454,1.52274,0.62849,0.27606,0    217         {765,1.47454,1.52274,0.62849,0.27606,0.15063,0.01938,-0.00097,7E-05,3E-05,0},
218         {770,1.47527,1.5377,0.634,0.2823,0.155    218         {770,1.47527,1.5377,0.634,0.2823,0.15572,0.02196,-0.00096,7E-05,3E-05,0},
219         {775,1.47506,1.55045,0.63926,0.28875,0    219         {775,1.47506,1.55045,0.63926,0.28875,0.16095,0.02458,-0.00094,8E-05,3E-05,0},
220         {780,1.47395,1.56132,0.64445,0.29542,0    220         {780,1.47395,1.56132,0.64445,0.29542,0.16628,0.02721,-0.0009,8E-05,3E-05,0},
221         {785,1.47197,1.57062,0.64973,0.30235,0    221         {785,1.47197,1.57062,0.64973,0.30235,0.17174,0.02983,-0.00086,9E-05,4E-05,0},
222         {790,1.46919,1.5787,0.65524,0.30957,0.    222         {790,1.46919,1.5787,0.65524,0.30957,0.1773,0.03242,-0.0008,9E-05,4E-05,0},
223         {795,1.46566,1.58579,0.66108,0.31709,0    223         {795,1.46566,1.58579,0.66108,0.31709,0.18297,0.03496,-0.00072,0.0001,4E-05,0},
224         {800,1.46143,1.59213,0.66731,0.32494,0    224         {800,1.46143,1.59213,0.66731,0.32494,0.18871,0.03743,-0.00063,0.0001,4E-05,0},
225         {805,1.45658,1.59801,0.67402,0.33312,0    225         {805,1.45658,1.59801,0.67402,0.33312,0.19452,0.03981,-0.00052,0.00011,4E-05,0},
226         {810,1.45114,1.60371,0.68124,0.34165,0    226         {810,1.45114,1.60371,0.68124,0.34165,0.20038,0.04209,-0.00039,0.00011,4E-05,0},
227         {815,1.44519,1.60936,0.68893,0.35054,0    227         {815,1.44519,1.60936,0.68893,0.35054,0.20629,0.04424,-0.00023,0.00011,4E-05,0},
228         {820,1.43878,1.61515,0.69704,0.35982,0    228         {820,1.43878,1.61515,0.69704,0.35982,0.21224,0.04625,-6E-05,0.00012,4E-05,0},
229         {825,1.43195,1.62131,0.70553,0.36949,0    229         {825,1.43195,1.62131,0.70553,0.36949,0.21824,0.04813,0.00015,0.00012,4E-05,0},
230         {830,1.42477,1.62791,0.71427,0.37956,0    230         {830,1.42477,1.62791,0.71427,0.37956,0.22429,0.04984,0.00038,0.00012,3E-05,0},
231         {835,1.41729,1.63504,0.72318,0.39006,0    231         {835,1.41729,1.63504,0.72318,0.39006,0.23038,0.05138,0.00064,0.00012,3E-05,0},
232         {840,1.40956,1.64282,0.73218,0.401,0.2    232         {840,1.40956,1.64282,0.73218,0.401,0.23651,0.05274,0.00093,0.00012,3E-05,0},
233         {845,1.40161,1.65132,0.74124,0.41239,0    233         {845,1.40161,1.65132,0.74124,0.41239,0.2427,0.05393,0.00126,0.00012,2E-05,0},
234         {850,1.3935,1.66056,0.7503,0.42424,0.2    234         {850,1.3935,1.66056,0.7503,0.42424,0.24895,0.05496,0.00162,0.00011,2E-05,0},
235         {855,1.38526,1.67054,0.75937,0.43656,0    235         {855,1.38526,1.67054,0.75937,0.43656,0.25527,0.05582,0.00201,0.00011,1E-05,0},
236         {860,1.37693,1.68128,0.76851,0.44938,0    236         {860,1.37693,1.68128,0.76851,0.44938,0.26166,0.05654,0.00245,0.0001,1E-05,0},
237         {865,1.36855,1.69275,0.77776,0.46269,0    237         {865,1.36855,1.69275,0.77776,0.46269,0.26815,0.05712,0.00292,9E-05,0,0},
238         {870,1.36017,1.70495,0.78724,0.47652,0    238         {870,1.36017,1.70495,0.78724,0.47652,0.27478,0.05759,0.00343,8E-05,-1E-05,0},
239         {875,1.35182,1.71783,0.79705,0.49087,0    239         {875,1.35182,1.71783,0.79705,0.49087,0.28156,0.05797,0.00399,7E-05,-2E-05,0},
240         {880,1.34355,1.73135,0.80728,0.50576,0    240         {880,1.34355,1.73135,0.80728,0.50576,0.28853,0.05828,0.00458,5E-05,-3E-05,0},
241         {885,1.33538,1.74545,0.81803,0.52118,0    241         {885,1.33538,1.74545,0.81803,0.52118,0.29569,0.05854,0.00522,4E-05,-4E-05,0},
242         {890,1.32737,1.76009,0.82938,0.53712,0    242         {890,1.32737,1.76009,0.82938,0.53712,0.30306,0.05876,0.00589,1E-05,-5E-05,0},
243         {895,1.31955,1.77516,0.8414,0.55357,0.    243         {895,1.31955,1.77516,0.8414,0.55357,0.31064,0.05896,0.00661,-1E-05,-7E-05,0},
244         {900,1.31195,1.79066,0.85422,0.57052,0    244         {900,1.31195,1.79066,0.85422,0.57052,0.31846,0.05916,0.00737,-4E-05,-8E-05,0},
245         {905,1.30462,1.80645,0.86774,0.58794,0    245         {905,1.30462,1.80645,0.86774,0.58794,0.32651,0.05937,0.00817,-7E-05,-0.0001,0},
246         {910,1.29758,1.82248,0.88204,0.60577,0    246         {910,1.29758,1.82248,0.88204,0.60577,0.33479,0.05957,0.009,-0.0001,-0.00012,0},
247         {915,1.29087,1.83865,0.89711,0.62396,0    247         {915,1.29087,1.83865,0.89711,0.62396,0.34327,0.05975,0.00987,-0.00014,-0.00014,0},
248         {920,1.28451,1.85489,0.91289,0.64244,0    248         {920,1.28451,1.85489,0.91289,0.64244,0.35194,0.05992,0.01077,-0.00019,-0.00016,0},
249         {925,1.27854,1.87113,0.92941,0.66113,0    249         {925,1.27854,1.87113,0.92941,0.66113,0.36077,0.06006,0.01169,-0.00024,-0.00018,0},
250         {930,1.27299,1.88729,0.94663,0.67995,0    250         {930,1.27299,1.88729,0.94663,0.67995,0.36971,0.06017,0.01263,-0.00029,-0.00021,0},
251         {935,1.26789,1.90333,0.9646,0.69881,0.    251         {935,1.26789,1.90333,0.9646,0.69881,0.37872,0.06024,0.01358,-0.00034,-0.00024,0},
252         {940,1.26327,1.91917,0.98329,0.71765,0    252         {940,1.26327,1.91917,0.98329,0.71765,0.38775,0.06027,0.01453,-0.00041,-0.00026,0},
253         {945,1.25914,1.93478,1.00276,0.73636,0    253         {945,1.25914,1.93478,1.00276,0.73636,0.39678,0.06025,0.01548,-0.00047,-0.00029,0},
254         {950,1.25553,1.95014,1.02309,0.75486,0    254         {950,1.25553,1.95014,1.02309,0.75486,0.40576,0.06017,0.01642,-0.00054,-0.00032,0},
255         {955,1.25247,1.96526,1.04444,0.77305,0    255         {955,1.25247,1.96526,1.04444,0.77305,0.41466,0.06003,0.01733,-0.00062,-0.00036,0},
256         {960,1.24996,1.98014,1.06695,0.79084,0    256         {960,1.24996,1.98014,1.06695,0.79084,0.42344,0.05981,0.0182,-0.0007,-0.00039,0},
257         {965,1.24804,1.99479,1.09078,0.80815,0    257         {965,1.24804,1.99479,1.09078,0.80815,0.43208,0.0595,0.01902,-0.00079,-0.00042,0},
258         {970,1.24672,2.00926,1.11616,0.82491,0    258         {970,1.24672,2.00926,1.11616,0.82491,0.44055,0.0591,0.01979,-0.00088,-0.00046,0},
259         {975,1.24601,2.02357,1.14324,0.84102,0    259         {975,1.24601,2.02357,1.14324,0.84102,0.44882,0.05858,0.02051,-0.00098,-0.00049,0},
260         {980,1.24594,2.03775,1.17221,0.8564,0.    260         {980,1.24594,2.03775,1.17221,0.8564,0.45685,0.05794,0.02117,-0.00108,-0.00053,0},
261         {985,1.2465,2.0518,1.20312,0.87098,0.4    261         {985,1.2465,2.0518,1.20312,0.87098,0.46461,0.05716,0.02178,-0.00118,-0.00057,0},
262         {990,1.24771,2.06568,1.23594,0.88468,0    262         {990,1.24771,2.06568,1.23594,0.88468,0.47206,0.05623,0.02233,-0.00129,-0.0006,0},
263         {995,1.24956,2.0793,1.27044,0.89743,0.    263         {995,1.24956,2.0793,1.27044,0.89743,0.47913,0.05514,0.02282,-0.0014,-0.00064,0},
264         {1000,1.25203,2.09256,1.30631,0.90917,    264         {1000,1.25203,2.09256,1.30631,0.90917,0.48577,0.05388,0.02322,-0.00151,-0.00067,0},
265         {1005,1.25511,2.10532,1.34305,0.91985,    265         {1005,1.25511,2.10532,1.34305,0.91985,0.49192,0.05246,0.02355,-0.00162,-0.0007,0},
266         {1010,1.2588,2.11731,1.37985,0.92939,0    266         {1010,1.2588,2.11731,1.37985,0.92939,0.49755,0.05085,0.02379,-0.00174,-0.00073,0},
267         {1015,1.26308,2.12833,1.41594,0.93779,    267         {1015,1.26308,2.12833,1.41594,0.93779,0.5026,0.04907,0.02393,-0.00185,-0.00076,0},
268         {1020,1.26795,2.13814,1.45041,0.94503,    268         {1020,1.26795,2.13814,1.45041,0.94503,0.50705,0.0471,0.02399,-0.00197,-0.00078,0},
269         {1025,1.27338,2.14656,1.48242,0.95109,    269         {1025,1.27338,2.14656,1.48242,0.95109,0.51088,0.04496,0.02395,-0.00208,-0.0008,0},
270         {1030,1.27938,2.15331,1.51104,0.95597,    270         {1030,1.27938,2.15331,1.51104,0.95597,0.51404,0.04264,0.02382,-0.00219,-0.00082,0},
271         {1035,1.28592,2.15832,1.53575,0.95964,    271         {1035,1.28592,2.15832,1.53575,0.95964,0.51651,0.04015,0.0236,-0.00229,-0.00083,0},
272         {1040,1.293,2.16155,1.55614,0.96211,0.    272         {1040,1.293,2.16155,1.55614,0.96211,0.51826,0.03747,0.02327,-0.00239,-0.00083,0},
273         {1045,1.30061,2.163,1.57204,0.96334,0.    273         {1045,1.30061,2.163,1.57204,0.96334,0.51926,0.03463,0.02286,-0.00249,-0.00083,0},
274         {1050,1.30872,2.16272,1.58343,0.96336,    274         {1050,1.30872,2.16272,1.58343,0.96336,0.51949,0.03161,0.02235,-0.00257,-0.00082,0},
275         {1055,1.31732,2.16081,1.5905,0.96221,0    275         {1055,1.31732,2.16081,1.5905,0.96221,0.51898,0.02843,0.02175,-0.00264,-0.00081,0},
276         {1060,1.32638,2.15741,1.59355,0.95996,    276         {1060,1.32638,2.15741,1.59355,0.95996,0.51774,0.0251,0.02108,-0.00271,-0.00078,0},
277         {1065,1.33587,2.15274,1.59306,0.95665,    277         {1065,1.33587,2.15274,1.59306,0.95665,0.51581,0.02165,0.02035,-0.00275,-0.00074,0},
278         {1070,1.34578,2.14696,1.58939,0.95236,    278         {1070,1.34578,2.14696,1.58939,0.95236,0.51321,0.01808,0.01956,-0.00278,-0.00069,0},
279         {1075,1.35609,2.14027,1.58309,0.94716,    279         {1075,1.35609,2.14027,1.58309,0.94716,0.50999,0.01443,0.01874,-0.0028,-0.00063,0},
280         {1080,1.36676,2.13286,1.57454,0.94113,    280         {1080,1.36676,2.13286,1.57454,0.94113,0.50619,0.01071,0.01787,-0.00279,-0.00056,0},
281         {1085,1.37779,2.12492,1.56429,0.93433,    281         {1085,1.37779,2.12492,1.56429,0.93433,0.50184,0.00695,0.01698,-0.00275,-0.00047,0},
282         {1090,1.38914,2.11665,1.55271,0.92686,    282         {1090,1.38914,2.11665,1.55271,0.92686,0.49697,0.00317,0.01608,-0.00269,-0.00037,0},
283         {1095,1.40076,2.10821,1.54028,0.91878,    283         {1095,1.40076,2.10821,1.54028,0.91878,0.49165,-0.00062,0.01518,-0.0026,-0.00025,0},
284         {1100,1.41261,2.09974,1.52734,0.9102,0    284         {1100,1.41261,2.09974,1.52734,0.9102,0.4859,-0.00439,0.01429,-0.00247,-0.00011,0},
285         {1105,1.42462,2.09138,1.51429,0.90121,    285         {1105,1.42462,2.09138,1.51429,0.90121,0.47978,-0.00813,0.01344,-0.00231,5E-05,0},
286         {1110,1.43677,2.08324,1.50142,0.89189,    286         {1110,1.43677,2.08324,1.50142,0.89189,0.47332,-0.01183,0.01264,-0.0021,0.00023,0},
287         {1115,1.44898,2.07537,1.48901,0.88234,    287         {1115,1.44898,2.07537,1.48901,0.88234,0.46657,-0.01546,0.0119,-0.00185,0.00044,0},
288         {1120,1.46122,2.06783,1.4773,0.87264,0    288         {1120,1.46122,2.06783,1.4773,0.87264,0.45958,-0.019,0.01126,-0.00154,0.00066,0},
289         {1125,1.47343,2.06063,1.46643,0.86287,    289         {1125,1.47343,2.06063,1.46643,0.86287,0.45238,-0.02245,0.01072,-0.00118,0.00091,0},
290         {1130,1.48556,2.05378,1.45653,0.85308,    290         {1130,1.48556,2.05378,1.45653,0.85308,0.44503,-0.02577,0.01029,-0.00076,0.00119,0},
291         {1135,1.49758,2.04726,1.44767,0.84333,    291         {1135,1.49758,2.04726,1.44767,0.84333,0.43758,-0.02897,0.00999,-0.00028,0.00149,0},
292         {1140,1.50943,2.04106,1.43987,0.83366,    292         {1140,1.50943,2.04106,1.43987,0.83366,0.43008,-0.03202,0.00984,0.00027,0.00183,0},
293         {1145,1.52109,2.03518,1.43314,0.82413,    293         {1145,1.52109,2.03518,1.43314,0.82413,0.4226,-0.03489,0.00984,0.00089,0.00219,0},
294         {1150,1.53252,2.0296,1.42743,0.81479,0    294         {1150,1.53252,2.0296,1.42743,0.81479,0.4152,-0.03756,0.01,0.00159,0.00258,0},
295         {1155,1.54367,2.02432,1.42269,0.8057,0    295         {1155,1.54367,2.02432,1.42269,0.8057,0.4079,-0.04002,0.01034,0.00238,0.003,0},
296         {1160,1.5545,2.01935,1.41883,0.79694,0    296         {1160,1.5545,2.01935,1.41883,0.79694,0.40076,-0.04224,0.01087,0.00325,0.00345,0},
297         {1165,1.56496,2.01471,1.41576,0.78856,    297         {1165,1.56496,2.01471,1.41576,0.78856,0.3938,-0.04419,0.01158,0.00421,0.00393,0},
298         {1170,1.57503,2.01042,1.41337,0.78064,    298         {1170,1.57503,2.01042,1.41337,0.78064,0.38707,-0.04587,0.01249,0.00526,0.00444,0},
299         {1175,1.58469,2.1965,1.41155,0.77321,0    299         {1175,1.58469,2.1965,1.41155,0.77321,0.38061,-0.04724,0.01361,0.00642,0.00498,0},
300         {1180,1.5939,2.00299,1.41016,0.76633,0    300         {1180,1.5939,2.00299,1.41016,0.76633,0.37444,-0.04829,0.01496,0.00767,0.00554,0},
301         {1185,1.60265,1.9999,1.40908,0.76005,0    301         {1185,1.60265,1.9999,1.40908,0.76005,0.3686,-0.049,0.01653,0.00902,0.00613,0},
302         {1190,1.61092,1.99724,1.40818,0.7544,0    302         {1190,1.61092,1.99724,1.40818,0.7544,0.36312,-0.04934,0.01834,0.01048,0.00674,1E-05},
303         {1195,1.6187,1.99499,1.40729,0.7494,0.    303         {1195,1.6187,1.99499,1.40729,0.7494,0.35805,-0.0493,0.02039,0.01204,0.00737,1E-05},
304         {1200,1.62597,1.99317,1.40635,0.74509,    304         {1200,1.62597,1.99317,1.40635,0.74509,0.35343,-0.04887,0.02268,0.01371,0.00801,1E-05},
305         {1205,1.63272,1.99173,1.40521,0.74151,    305         {1205,1.63272,1.99173,1.40521,0.74151,0.34931,-0.04801,0.02523,0.01547,0.00866,1E-05},
306         {1210,1.63895,1.99063,1.40379,0.73867,    306         {1210,1.63895,1.99063,1.40379,0.73867,0.34575,-0.0467,0.02803,0.01733,0.00932,1E-05},
307         {1215,1.64465,1.98982,1.40203,0.73662,    307         {1215,1.64465,1.98982,1.40203,0.73662,0.34282,-0.0449,0.03109,0.01929,0.00997,1E-05},
308         {1220,1.64984,1.98925,1.39991,0.73538,    308         {1220,1.64984,1.98925,1.39991,0.73538,0.34058,-0.04259,0.03443,0.02134,0.01062,2E-05},
309         {1225,1.65456,1.98885,1.39743,0.73499,    309         {1225,1.65456,1.98885,1.39743,0.73499,0.33906,-0.03974,0.03803,0.02346,0.01124,2E-05},
310         {1230,1.65885,1.98857,1.39463,0.73546,    310         {1230,1.65885,1.98857,1.39463,0.73546,0.33833,-0.03633,0.0419,0.02566,0.01185,2E-05},
311         {1235,1.66276,1.9884,1.39163,0.73683,0    311         {1235,1.66276,1.9884,1.39163,0.73683,0.33842,-0.03235,0.04604,0.02791,0.01241,3E-05},
312         {1240,1.66633,1.98832,1.38853,0.7391,0    312         {1240,1.66633,1.98832,1.38853,0.7391,0.33937,-0.0278,0.05045,0.03022,0.01294,3E-05},
313         {1245,1.66959,1.98831,1.38547,0.74229,    313         {1245,1.66959,1.98831,1.38547,0.74229,0.34116,-0.02272,0.05512,0.03256,0.01341,3E-05},
314         {1250,1.67261,1.98841,1.38262,0.74638,    314         {1250,1.67261,1.98841,1.38262,0.74638,0.34381,-0.01713,0.06004,0.03492,0.01381,4E-05},
315         {1255,1.67544,1.98864,1.38016,0.75139,    315         {1255,1.67544,1.98864,1.38016,0.75139,0.34735,-0.01103,0.06523,0.03731,0.01416,4E-05},
316         {1260,1.67811,1.98907,1.37827,0.75735,    316         {1260,1.67811,1.98907,1.37827,0.75735,0.35181,-0.00444,0.07068,0.03972,0.01444,4E-05},
317         {1265,1.6807,1.98975,1.37713,0.76427,0    317         {1265,1.6807,1.98975,1.37713,0.76427,0.35722,0.00267,0.07642,0.04215,0.01467,5E-05},
318         {1270,1.68325,1.99077,1.37693,0.77218,    318         {1270,1.68325,1.99077,1.37693,0.77218,0.36363,0.0103,0.08243,0.04461,0.01486,5E-05},
319         {1275,1.68581,1.99224,1.37785,0.7811,0    319         {1275,1.68581,1.99224,1.37785,0.7811,0.37107,0.01845,0.08874,0.0471,0.01499,5E-05},
320         {1280,1.68844,1.99422,1.38006,0.79104,    320         {1280,1.68844,1.99422,1.38006,0.79104,0.37957,0.02713,0.09535,0.04962,0.01508,6E-05},
321         {1285,1.69117,1.99685,1.38376,0.80203,    321         {1285,1.69117,1.99685,1.38376,0.80203,0.38917,0.03633,0.10226,0.05216,0.01514,6E-05},
322         {1290,1.69404,2.00021,1.38909,0.81409,    322         {1290,1.69404,2.00021,1.38909,0.81409,0.39987,0.04604,0.10947,0.05472,0.01516,7E-05},
323         {1295,1.69705,2.00443,1.39628,0.8272,0    323         {1295,1.69705,2.00443,1.39628,0.8272,0.41163,0.05625,0.11694,0.05729,0.01514,7E-05},
324         {1300,1.70021,2.00965,1.4055,0.84135,0    324         {1300,1.70021,2.00965,1.4055,0.84135,0.42441,0.06693,0.12467,0.05985,0.01508,8E-05},
325         {1305,1.70353,2.01597,1.41698,0.85652,    325         {1305,1.70353,2.01597,1.41698,0.85652,0.43817,0.07806,0.13261,0.0624,0.015,8E-05},
326         {1310,1.70703,2.02352,1.43096,0.87269,    326         {1310,1.70703,2.02352,1.43096,0.87269,0.45284,0.08962,0.14074,0.06492,0.01488,8E-05},
327         {1315,1.7107,2.0324,1.44765,0.88983,0.    327         {1315,1.7107,2.0324,1.44765,0.88983,0.46839,0.10157,0.14903,0.0674,0.01473,9E-05},
328         {1320,1.71455,2.04268,1.46721,0.90794,    328         {1320,1.71455,2.04268,1.46721,0.90794,0.48476,0.11389,0.15746,0.06983,0.01456,9E-05},
329         {1325,1.7186,2.0544,1.48977,0.92699,0.    329         {1325,1.7186,2.0544,1.48977,0.92699,0.50191,0.12657,0.166,0.0722,0.01435,9E-05},
330         {1330,1.72283,2.06752,1.51535,0.94694,    330         {1330,1.72283,2.06752,1.51535,0.94694,0.5198,0.13956,0.17463,0.07448,0.01412,0.0001},
331         {1335,1.72726,2.082,1.54394,0.96778,0.    331         {1335,1.72726,2.082,1.54394,0.96778,0.53838,0.15284,0.18333,0.07669,0.01386,0.0001},
332         {1340,1.73188,2.09761,1.57517,0.98947,    332         {1340,1.73188,2.09761,1.57517,0.98947,0.55762,0.16638,0.1921,0.07879,0.01357,0.0001},
333         {1345,1.73669,2.11416,1.60872,1.01198,    333         {1345,1.73669,2.11416,1.60872,1.01198,0.57747,0.18014,0.2009,0.08078,0.01327,0.0001},
334         {1350,1.74168,2.13141,1.64411,1.03528,    334         {1350,1.74168,2.13141,1.64411,1.03528,0.5979,0.19409,0.20973,0.08266,0.01294,0.0001},
335         {1355,1.74687,2.14911,1.68074,1.05933,    335         {1355,1.74687,2.14911,1.68074,1.05933,0.61886,0.2082,0.21857,0.0844,0.01259,0.0001},
336         {1360,1.75224,2.16699,1.71798,1.08412,    336         {1360,1.75224,2.16699,1.71798,1.08412,0.64031,0.22244,0.22739,0.086,0.01222,0.0001},
337         {1365,1.7578,2.1849,1.75528,1.1096,0.6    337         {1365,1.7578,2.1849,1.75528,1.1096,0.66221,0.23678,0.23619,0.08744,0.01183,0.0001},
338         {1370,1.76355,2.20267,1.79212,1.13575,    338         {1370,1.76355,2.20267,1.79212,1.13575,0.68453,0.25118,0.24495,0.08872,0.01142,0.0001},
339         {1375,1.76948,2.22025,1.82818,1.16253,    339         {1375,1.76948,2.22025,1.82818,1.16253,0.70722,0.26561,0.25364,0.08982,0.011,9E-05},
340         {1380,1.7756,2.23755,1.863,1.18992,0.7    340         {1380,1.7756,2.23755,1.863,1.18992,0.73024,0.28004,0.26226,0.09073,0.01056,9E-05},
341         {1385,1.78191,2.25468,1.89666,1.21789,    341         {1385,1.78191,2.25468,1.89666,1.21789,0.75355,0.29444,0.27078,0.09145,0.0101,8E-05},
342         {1390,1.78839,2.27153,1.92896,1.2464,0    342         {1390,1.78839,2.27153,1.92896,1.2464,0.77711,0.30879,0.2792,0.09196,0.00964,7E-05},
343         {1395,1.79506,2.28807,1.95988,1.2754,0    343         {1395,1.79506,2.28807,1.95988,1.2754,0.80087,0.32308,0.28749,0.0923,0.00916,6E-05},
344         {1400,1.80191,2.3043,1.9895,1.30486,0.    344         {1400,1.80191,2.3043,1.9895,1.30486,0.8248,0.33732,0.29565,0.09246,0.00868,5E-05},
345         {1405,1.80892,2.32031,2.01814,1.33474,    345         {1405,1.80892,2.32031,2.01814,1.33474,0.84883,0.35149,0.30367,0.09248,0.00819,4E-05},
346         {1410,1.81611,2.33611,2.04597,1.36499,    346         {1410,1.81611,2.33611,2.04597,1.36499,0.87294,0.36561,0.31154,0.09237,0.00769,3E-05},
347         {1415,1.82347,2.35183,2.0734,1.39557,0    347         {1415,1.82347,2.35183,2.0734,1.39557,0.89707,0.37966,0.31924,0.09214,0.0072,1E-05},
348         {1420,1.83098,2.36761,2.10085,1.42644,    348         {1420,1.83098,2.36761,2.10085,1.42644,0.92119,0.39364,0.32677,0.09182,0.0067,-1E-05},
349         {1425,1.83866,2.3837,2.12881,1.45756,0    349         {1425,1.83866,2.3837,2.12881,1.45756,0.94524,0.40755,0.33411,0.09141,0.00621,-3E-05},
350         {1430,1.8465,2.40026,2.15761,1.48889,0    350         {1430,1.8465,2.40026,2.15761,1.48889,0.96918,0.42139,0.34125,0.09094,0.00572,-5E-05},
351         {1435,1.85448,2.41772,2.18799,1.52038,    351         {1435,1.85448,2.41772,2.18799,1.52038,0.99296,0.43515,0.34818,0.09043,0.00524,-7E-05},
352         {1440,1.86261,2.43632,2.22012,1.55198,    352         {1440,1.86261,2.43632,2.22012,1.55198,1.01656,0.44883,0.35489,0.08987,0.00476,-0.0001},
353         {1445,1.87088,2.4564,2.25442,1.58366,1    353         {1445,1.87088,2.4564,2.25442,1.58366,1.03994,0.46243,0.36137,0.08928,0.00429,-0.00013},
354         {1450,1.87926,2.47821,2.29095,1.61534,    354         {1450,1.87926,2.47821,2.29095,1.61534,1.06307,0.47594,0.36761,0.08864,0.00384,-0.00016},
355         {1455,1.88775,2.50203,2.32979,1.64699,    355         {1455,1.88775,2.50203,2.32979,1.64699,1.08592,0.48935,0.37359,0.08797,0.00338,-0.0002},
356         {1460,1.89633,2.52791,2.37059,1.67854,    356         {1460,1.89633,2.52791,2.37059,1.67854,1.10846,0.50267,0.37931,0.08724,0.00294,-0.00023},
357         {1465,1.90498,2.55575,2.41275,1.70994,    357         {1465,1.90498,2.55575,2.41275,1.70994,1.13066,0.51587,0.38474,0.08647,0.0025,-0.00027},
358         {1470,1.9137,2.58528,2.45552,1.74115,1    358         {1470,1.9137,2.58528,2.45552,1.74115,1.1525,0.52897,0.38989,0.08565,0.00207,-0.00031},
359         {1475,1.92247,2.61583,2.49765,1.77211,    359         {1475,1.92247,2.61583,2.49765,1.77211,1.17394,0.54196,0.39473,0.08478,0.00164,-0.00036},
360         {1480,1.93128,2.64682,2.53822,1.80277,    360         {1480,1.93128,2.64682,2.53822,1.80277,1.19496,0.55482,0.39926,0.08386,0.00122,-0.00041},
361         {1485,1.9401,2.67749,2.57631,1.83306,1    361         {1485,1.9401,2.67749,2.57631,1.83306,1.21553,0.56756,0.40347,0.08288,0.00081,-0.00046},
362         {1490,1.94894,2.70719,2.61131,1.86296,    362         {1490,1.94894,2.70719,2.61131,1.86296,1.23562,0.58015,0.40734,0.08184,0.0004,-0.00051},
363         {1495,1.95779,2.73518,2.64271,1.89239,    363         {1495,1.95779,2.73518,2.64271,1.89239,1.2552,0.59258,0.41088,0.08076,0,-0.00057},
364         {1500,1.96663,2.76107,2.67053,1.92133,    364         {1500,1.96663,2.76107,2.67053,1.92133,1.27426,0.60481,0.4141,0.07962,-0.0004,-0.00063},
365         {1505,1.97546,2.78466,2.6951,1.94971,1    365         {1505,1.97546,2.78466,2.6951,1.94971,1.29278,0.61681,0.41699,0.07844,-0.0008,-0.00069},
366         {1510,1.98427,2.80574,2.71662,1.9775,1    366         {1510,1.98427,2.80574,2.71662,1.9775,1.31074,0.62856,0.41955,0.07722,-0.00119,-0.00076},
367         {1515,1.99305,2.82446,2.73579,2.00463,    367         {1515,1.99305,2.82446,2.73579,2.00463,1.3281,0.64003,0.42179,0.07595,-0.00158,-0.00083},
368         {1520,2.00181,2.84098,2.7532,2.03108,1    368         {1520,2.00181,2.84098,2.7532,2.03108,1.34486,0.65118,0.42371,0.07465,-0.00196,-0.0009},
369         {1525,2.01052,2.85555,2.76943,2.05678,    369         {1525,2.01052,2.85555,2.76943,2.05678,1.36098,0.66199,0.42531,0.07332,-0.00235,-0.00098},
370         {1530,2.01919,2.86842,2.78501,2.08169,    370         {1530,2.01919,2.86842,2.78501,2.08169,1.37646,0.67243,0.4266,0.07196,-0.00273,-0.00106},
371         {1535,2.0278,2.87996,2.80048,2.10576,1    371         {1535,2.0278,2.87996,2.80048,2.10576,1.39126,0.68247,0.42757,0.07057,-0.0031,-0.00114},
372         {1540,2.03635,2.8903,2.8161,2.12896,1.    372         {1540,2.03635,2.8903,2.8161,2.12896,1.40537,0.6921,0.42823,0.06915,-0.00348,-0.00123},
373         {1545,2.04483,2.89977,2.83218,2.15129,    373         {1545,2.04483,2.89977,2.83218,2.15129,1.41881,0.7013,0.42858,0.06771,-0.00385,-0.00132},
374         {1550,2.05324,2.90864,2.84896,2.17277,    374         {1550,2.05324,2.90864,2.84896,2.17277,1.43157,0.71009,0.42864,0.06625,-0.00422,-0.00141},
375         {1555,2.06156,2.91713,2.86654,2.1934,1    375         {1555,2.06156,2.91713,2.86654,2.1934,1.44367,0.71845,0.42841,0.06476,-0.00459,-0.00151},
376         {1560,2.06979,2.92546,2.88488,2.2132,1    376         {1560,2.06979,2.92546,2.88488,2.2132,1.45512,0.7264,0.4279,0.06325,-0.00494,-0.00161},
377         {1565,2.07792,2.93381,2.90387,2.23216,    377         {1565,2.07792,2.93381,2.90387,2.23216,1.46592,0.73394,0.42712,0.06171,-0.00529,-0.00171},
378         {1570,2.08595,2.94231,2.92325,2.25032,    378         {1570,2.08595,2.94231,2.92325,2.25032,1.47609,0.74106,0.42607,0.06015,-0.00563,-0.00181},
379         {1575,2.09387,2.95113,2.94275,2.26766,    379         {1575,2.09387,2.95113,2.94275,2.26766,1.48563,0.74777,0.42477,0.05856,-0.00596,-0.00191},
380         {1580,2.10166,2.96017,2.96173,2.28421,    380         {1580,2.10166,2.96017,2.96173,2.28421,1.49455,0.75406,0.42322,0.05695,-0.00628,-0.00201},
381         {1585,2.10933,2.96958,2.97993,2.29998,    381         {1585,2.10933,2.96958,2.97993,2.29998,1.50287,0.75995,0.42143,0.05532,-0.00659,-0.00211},
382         {1590,2.11686,2.97912,2.99666,2.31496,    382         {1590,2.11686,2.97912,2.99666,2.31496,1.51058,0.76542,0.41941,0.05366,-0.00689,-0.00221},
383         {1595,2.12426,2.98872,3.01163,2.32918,    383         {1595,2.12426,2.98872,3.01163,2.32918,1.5177,0.77048,0.41717,0.05198,-0.00716,-0.0023},
384         {1600,2.1315,2.99818,3.0245,2.34264,1.    384         {1600,2.1315,2.99818,3.0245,2.34264,1.52424,0.77514,0.41471,0.05028,-0.00743,-0.0024},
385         {1605,2.13859,3.00725,3.03509,2.35536,    385         {1605,2.13859,3.00725,3.03509,2.35536,1.53021,0.77939,0.41204,0.04855,-0.00768,-0.00249},
386         {1610,2.14551,3.01569,3.04332,2.36733,    386         {1610,2.14551,3.01569,3.04332,2.36733,1.53561,0.78324,0.40917,0.0468,-0.00791,-0.00258},
387         {1615,2.15227,3.02325,3.04931,2.37858,    387         {1615,2.15227,3.02325,3.04931,2.37858,1.54046,0.78668,0.40611,0.04503,-0.00811,-0.00267},
388         {1620,2.15884,3.02974,3.05324,2.38911,    388         {1620,2.15884,3.02974,3.05324,2.38911,1.54476,0.78972,0.40287,0.04324,-0.0083,-0.00275},
389         {1625,2.16523,3.03506,3.0554,2.39893,1    389         {1625,2.16523,3.03506,3.0554,2.39893,1.54852,0.79235,0.39945,0.04142,-0.00847,-0.00282},
390         {1630,2.17143,3.03914,3.05609,2.40805,    390         {1630,2.17143,3.03914,3.05609,2.40805,1.55175,0.79459,0.39587,0.03958,-0.00862,-0.0029},
391         {1635,2.17743,3.04206,3.05566,2.41649,    391         {1635,2.17743,3.04206,3.05566,2.41649,1.55446,0.79642,0.39212,0.03772,-0.00874,-0.00296},
392         {1640,2.18323,3.0438,3.05428,2.42425,1    392         {1640,2.18323,3.0438,3.05428,2.42425,1.55667,0.79786,0.38823,0.03583,-0.00884,-0.00302},
393         {1645,2.18883,3.04456,3.05224,2.43135,    393         {1645,2.18883,3.04456,3.05224,2.43135,1.5584,0.79893,0.38421,0.03394,-0.0089,-0.00308},
394         {1650,2.19425,3.04449,3.04972,2.43784,    394         {1650,2.19425,3.04449,3.04972,2.43784,1.55969,0.79964,0.38008,0.03203,-0.00894,-0.00312},
395         {1655,2.19949,3.04374,3.04683,2.44372,    395         {1655,2.19949,3.04374,3.04683,2.44372,1.56058,0.80001,0.37586,0.03012,-0.00894,-0.00316},
396         {1660,2.20457,3.04245,3.04363,2.44904,    396         {1660,2.20457,3.04245,3.04363,2.44904,1.5611,0.80007,0.37157,0.02822,-0.0089,-0.00318},
397         {1665,2.20949,3.04076,3.04018,2.4538,1    397         {1665,2.20949,3.04076,3.04018,2.4538,1.56127,0.79983,0.36722,0.02633,-0.00882,-0.00319},
398         {1670,2.21426,3.03876,3.03649,2.45804,    398         {1670,2.21426,3.03876,3.03649,2.45804,1.56113,0.79931,0.36285,0.02445,-0.00869,-0.00319},
399         {1675,2.21889,3.03657,3.03256,2.46178,    399         {1675,2.21889,3.03657,3.03256,2.46178,1.56073,0.79853,0.35847,0.0226,-0.00851,-0.00317},
400         {1680,2.22339,3.03424,3.02844,2.46505,    400         {1680,2.22339,3.03424,3.02844,2.46505,1.56008,0.79752,0.35409,0.02078,-0.00828,-0.00314},
401         {1685,2.22778,3.03182,3.02406,2.46788,    401         {1685,2.22778,3.03182,3.02406,2.46788,1.55922,0.79629,0.34975,0.019,-0.00799,-0.00309},
402         {1690,2.23204,3.02936,3.01948,2.47029,    402         {1690,2.23204,3.02936,3.01948,2.47029,1.5582,0.79487,0.34545,0.01726,-0.00764,-0.00302},
403         {1695,2.2362,3.02689,3.01476,2.47235,1    403         {1695,2.2362,3.02689,3.01476,2.47235,1.55707,0.79332,0.34124,0.0156,-0.00723,-0.00294},
404         {1700,2.24024,3.02446,3.00998,2.47413,    404         {1700,2.24024,3.02446,3.00998,2.47413,1.55589,0.79171,0.33714,0.01403,-0.00674,-0.00283},
405         {1705,2.24417,3.0221,3.00522,2.47568,1    405         {1705,2.24417,3.0221,3.00522,2.47568,1.55472,0.79009,0.33318,0.01259,-0.00618,-0.00269},
406         {1710,2.24799,3.01984,3.00059,2.47707,    406         {1710,2.24799,3.01984,3.00059,2.47707,1.55362,0.78852,0.3294,0.01129,-0.00554,-0.00253},
407         {1715,2.2517,3.01771,2.99623,2.47837,1    407         {1715,2.2517,3.01771,2.99623,2.47837,1.55264,0.78708,0.32582,0.01016,-0.00482,-0.00234},
408         {1720,2.2553,3.01574,2.99226,2.47962,1    408         {1720,2.2553,3.01574,2.99226,2.47962,1.55186,0.78581,0.32248,0.00922,-0.004,-0.00212},
409         {1725,2.25879,3.01395,2.98884,2.48091,    409         {1725,2.25879,3.01395,2.98884,2.48091,1.55133,0.78478,0.31941,0.00851,-0.00308,-0.00186},
410         {1730,2.26217,3.01237,2.98611,2.48228,    410         {1730,2.26217,3.01237,2.98611,2.48228,1.55111,0.78405,0.31663,0.00804,-0.00207,-0.00158},
411         {1735,2.26543,3.01101,2.98419,2.48381,    411         {1735,2.26543,3.01101,2.98419,2.48381,1.55125,0.78368,0.31418,0.00784,-0.00095,-0.00125},
412         {1740,2.2686,3.00988,2.98319,2.48555,1    412         {1740,2.2686,3.00988,2.98319,2.48555,1.55181,0.78373,0.31209,0.00793,0.00028,-0.00089},
413         {1745,2.27167,3.00897,2.98323,2.48753,    413         {1745,2.27167,3.00897,2.98323,2.48753,1.55279,0.7842,0.31036,0.00833,0.00164,-0.00048},
414         {1750,2.27466,3.00828,2.98436,2.48976,    414         {1750,2.27466,3.00828,2.98436,2.48976,1.5542,0.78513,0.309,0.00904,0.00312,-4E-05},
415         {1755,2.27759,3.00777,2.9866,2.49229,1    415         {1755,2.27759,3.00777,2.9866,2.49229,1.55605,0.78653,0.30801,0.01008,0.00473,0.00046},
416         {1760,2.28046,3.00743,2.98993,2.49512,    416         {1760,2.28046,3.00743,2.98993,2.49512,1.55834,0.78842,0.30742,0.01146,0.00649,0.00101},
417         {1765,2.28329,3.00722,2.99427,2.49831,    417         {1765,2.28329,3.00722,2.99427,2.49831,1.56109,0.7908,0.30721,0.0132,0.0084,0.00161},
418         {1770,2.2861,3.00712,2.99948,2.50185,1    418         {1770,2.2861,3.00712,2.99948,2.50185,1.56429,0.79371,0.30741,0.0153,0.01047,0.00226},
419         {1775,2.28889,3.00709,3.00539,2.5058,1    419         {1775,2.28889,3.00709,3.00539,2.5058,1.56796,0.79716,0.30802,0.01778,0.0127,0.00297},
420         {1780,2.29169,3.1971,3.01177,2.51016,1    420         {1780,2.29169,3.1971,3.01177,2.51016,1.5721,0.80117,0.30905,0.02065,0.01511,0.00375},
421         {1785,2.2945,3.00713,3.01836,2.51497,1    421         {1785,2.2945,3.00713,3.01836,2.51497,1.57672,0.80575,0.3105,0.02392,0.01769,0.00458},
422         {1790,2.29733,3.00727,3.02501,2.52024,    422         {1790,2.29733,3.00727,3.02501,2.52024,1.58182,0.81092,0.31238,0.0276,0.02047,0.00549},
423         {1795,2.30019,3.00743,3.03129,2.52598,    423         {1795,2.30019,3.00743,3.03129,2.52598,1.58741,0.81667,0.3147,0.03168,0.02343,0.00646},
424         {1800,2.30308,3.00769,3.03701,2.53218,    424         {1800,2.30308,3.00769,3.03701,2.53218,1.5935,0.82299,0.31748,0.03616,0.02657,0.0075},
425         {1805,2.30599,3.1981,3.04192,2.53883,1    425         {1805,2.30599,3.1981,3.04192,2.53883,1.60008,0.82989,0.3207,0.04101,0.0299,0.00861},
426         {1810,2.30894,3.00875,3.04579,2.54593,    426         {1810,2.30894,3.00875,3.04579,2.54593,1.60716,0.83734,0.3244,0.04625,0.03341,0.00979},
427         {1815,2.31191,3.00975,3.04841,2.55348,    427         {1815,2.31191,3.00975,3.04841,2.55348,1.61474,0.84535,0.32856,0.05185,0.0371,0.01104},
428         {1820,2.31492,3.01122,3.04963,2.56148,    428         {1820,2.31492,3.01122,3.04963,2.56148,1.62284,0.85391,0.3332,0.05781,0.04098,0.01236},
429         {1825,2.31795,3.01333,3.04935,2.56991,    429         {1825,2.31795,3.01333,3.04935,2.56991,1.63144,0.86302,0.33833,0.06413,0.04503,0.01375},
430         {1830,2.32102,3.01624,3.04762,2.57878,    430         {1830,2.32102,3.01624,3.04762,2.57878,1.64056,0.87266,0.34395,0.07078,0.04926,0.01521},
431         {1835,2.32413,3.02014,3.04457,2.58808,    431         {1835,2.32413,3.02014,3.04457,2.58808,1.6502,0.88283,0.35007,0.07777,0.05367,0.01675},
432         {1840,2.32727,3.02528,3.04081,2.59781,    432         {1840,2.32727,3.02528,3.04081,2.59781,1.66037,0.89353,0.35671,0.08509,0.05826,0.01836},
433         {1845,2.33044,3.03168,3.03673,2.60797,    433         {1845,2.33044,3.03168,3.03673,2.60797,1.67106,0.90475,0.36386,0.09272,0.06301,0.02005},
434         {1850,2.33366,3.03949,3.03328,2.61854,    434         {1850,2.33366,3.03949,3.03328,2.61854,1.68229,0.91647,0.37154,0.10066,0.06795,0.02181},
435         {1855,2.33691,3.04874,3.03143,2.62953,    435         {1855,2.33691,3.04874,3.03143,2.62953,1.69405,0.9287,0.37975,0.1089,0.07305,0.02365},
436         {1860,2.3402,3.0594,3.03217,2.64094,1.    436         {1860,2.3402,3.0594,3.03217,2.64094,1.70635,0.94143,0.3885,0.11743,0.07833,0.02557},
437         {1865,2.34353,3.07137,3.03635,2.65275,    437         {1865,2.34353,3.07137,3.03635,2.65275,1.71919,0.95466,0.3978,0.12624,0.08378,0.02756},
438         {1870,2.3469,3.08454,3.04459,2.66497,1    438         {1870,2.3469,3.08454,3.04459,2.66497,1.73259,0.96837,0.40765,0.13533,0.0894,0.02963},
439         {1875,2.35032,3.09878,3.05724,2.67759,    439         {1875,2.35032,3.09878,3.05724,2.67759,1.74653,0.98256,0.41807,0.14468,0.09518,0.03179},
440         {1880,2.35378,3.11396,3.07431,2.69061,    440         {1880,2.35378,3.11396,3.07431,2.69061,1.76103,0.99722,0.42906,0.15429,0.10114,0.03402},
441         {1885,2.35728,3.12996,3.09565,2.70402,    441         {1885,2.35728,3.12996,3.09565,2.70402,1.77609,1.01235,0.44063,0.16415,0.10726,0.03633},
442         {1890,2.36083,3.1466,3.12083,2.71781,1    442         {1890,2.36083,3.1466,3.12083,2.71781,1.7917,1.02793,0.45277,0.17425,0.11354,0.03872},
443         {1895,2.3644,3.16367,3.14918,2.73195,1    443         {1895,2.3644,3.16367,3.14918,2.73195,1.80781,1.04392,0.46544,0.18455,0.11996,0.04119},
444         {1900,2.36801,3.181,3.18008,2.74642,1.    444         {1900,2.36801,3.181,3.18008,2.74642,1.82439,1.06029,0.4786,0.19504,0.1265,0.04371},
445         {1905,2.37163,3.19834,3.21275,2.76118,    445         {1905,2.37163,3.19834,3.21275,2.76118,1.84138,1.07698,0.4922,0.20568,0.13315,0.0463},
446         {1910,2.37525,3.2154,3.24639,2.77619,1    446         {1910,2.37525,3.2154,3.24639,2.77619,1.85874,1.09395,0.50619,0.21644,0.13988,0.04894},
447         {1915,2.37887,3.23188,3.28016,2.79144,    447         {1915,2.37887,3.23188,3.28016,2.79144,1.87643,1.11116,0.52052,0.2273,0.14668,0.05162},
448         {1920,2.38247,3.24747,3.31324,2.80688,    448         {1920,2.38247,3.24747,3.31324,2.80688,1.89439,1.12858,0.53516,0.23822,0.15352,0.05433},
449         {1925,2.38605,3.26188,3.34488,2.82248,    449         {1925,2.38605,3.26188,3.34488,2.82248,1.91258,1.14615,0.55006,0.24919,0.16039,0.05707},
450         {1930,2.3896,3.27487,3.37443,2.83822,1    450         {1930,2.3896,3.27487,3.37443,2.83822,1.93096,1.16384,0.56516,0.26017,0.16726,0.05983},
451         {1935,2.3931,3.28629,3.40141,2.85406,1    451         {1935,2.3931,3.28629,3.40141,2.85406,1.94947,1.18161,0.58043,0.27113,0.17413,0.0626},
452         {1940,2.39655,3.29596,3.42532,2.86997,    452         {1940,2.39655,3.29596,3.42532,2.86997,1.96808,1.1994,0.59582,0.28205,0.18096,0.06538},
453         {1945,2.39994,3.30388,3.44597,2.88593,    453         {1945,2.39994,3.30388,3.44597,2.88593,1.98673,1.21718,0.61129,0.2929,0.18774,0.06815},
454         {1950,2.40326,3.31012,3.46328,2.90189,    454         {1950,2.40326,3.31012,3.46328,2.90189,2.00539,1.23491,0.62678,0.30364,0.19445,0.07092},
455         {1955,2.4065,3.31482,3.47731,2.91783,2    455         {1955,2.4065,3.31482,3.47731,2.91783,2.02399,1.25255,0.64225,0.31426,0.20108,0.07366},
456         {1960,2.40965,3.31816,3.48822,2.93371,    456         {1960,2.40965,3.31816,3.48822,2.93371,2.04251,1.27005,0.65766,0.32472,0.20759,0.07638},
457         {1965,2.41269,3.32044,3.49632,2.94951,    457         {1965,2.41269,3.32044,3.49632,2.94951,2.06089,1.28736,0.67296,0.335,0.21398,0.07907},
458         {1970,2.41563,3.32179,3.50179,2.96519,    458         {1970,2.41563,3.32179,3.50179,2.96519,2.07908,1.30446,0.68811,0.34507,0.22022,0.08172},
459         {1975,2.41845,3.32254,3.50506,2.98073,    459         {1975,2.41845,3.32254,3.50506,2.98073,2.09704,1.32129,0.70305,0.35489,0.22629,0.08432},
460         {1980,2.42114,3.32293,3.50649,2.99609,    460         {1980,2.42114,3.32293,3.50649,2.99609,2.11473,1.33781,0.71775,0.36445,0.23218,0.08686},
461         {1985,2.42369,3.32318,3.50645,3.01123,    461         {1985,2.42369,3.32318,3.50645,3.01123,2.13209,1.35398,0.73216,0.3737,0.23786,0.08934},
462         {1990,2.4261,3.32346,3.50529,3.02615,2    462         {1990,2.4261,3.32346,3.50529,3.02615,2.1491,1.36979,0.74625,0.38265,0.24332,0.09175},
463         {1995,2.42837,3.3239,3.50335,3.04085,2    463         {1995,2.42837,3.3239,3.50335,3.04085,2.16577,1.38523,0.76001,0.3913,0.24858,0.0941},
464         {2000,2.43052,3.3246,3.50095,3.05532,2    464         {2000,2.43052,3.3246,3.50095,3.05532,2.18212,1.40032,0.77347,0.39966,0.25362,0.09637},
465         {2005,2.43254,3.32561,3.49838,3.0696,2    465         {2005,2.43254,3.32561,3.49838,3.0696,2.19817,1.4151,0.78662,0.40776,0.25847,0.09857},
466         {2010,2.43445,3.32696,3.49592,3.08368,    466         {2010,2.43445,3.32696,3.49592,3.08368,2.21393,1.42957,0.79948,0.41561,0.26313,0.10071},
467         {2015,2.43626,3.32865,3.49379,3.09757,    467         {2015,2.43626,3.32865,3.49379,3.09757,2.22941,1.44376,0.81205,0.42322,0.2676,0.10277},
468         {2020,2.43797,3.33067,3.4922,3.11129,2    468         {2020,2.43797,3.33067,3.4922,3.11129,2.24465,1.45768,0.82434,0.43061,0.2719,0.10477},
469         {2025,2.43959,3.33296,3.49131,3.12484,    469         {2025,2.43959,3.33296,3.49131,3.12484,2.25964,1.47135,0.83637,0.4378,0.27602,0.10669},
470         {2030,2.44113,3.3355,3.49125,3.13823,2    470         {2030,2.44113,3.3355,3.49125,3.13823,2.27442,1.48479,0.84813,0.4448,0.27998,0.10855},
471         {2035,2.4426,3.33826,3.49211,3.15148,2    471         {2035,2.4426,3.33826,3.49211,3.15148,2.28899,1.49802,0.85963,0.45163,0.28378,0.11033},
472         {2040,2.444,3.34111,3.49372,3.16458,2.    472         {2040,2.444,3.34111,3.49372,3.16458,2.30338,1.51106,0.87089,0.45831,0.28743,0.11205},
473         {2045,2.44535,3.34405,3.49616,3.17756,    473         {2045,2.44535,3.34405,3.49616,3.17756,2.31759,1.52392,0.88191,0.46485,0.29094,0.11369},
474         {2050,2.44665,3.34702,3.49935,3.19042,    474         {2050,2.44665,3.34702,3.49935,3.19042,2.33166,1.53664,0.8927,0.47126,0.29431,0.11527},
475         {2055,2.44791,3.34996,3.50318,3.20317,    475         {2055,2.44791,3.34996,3.50318,3.20317,2.34559,1.54921,0.90327,0.47758,0.29755,0.11678},
476         {2060,2.44914,3.35283,3.50751,3.21581,    476         {2060,2.44914,3.35283,3.50751,3.21581,2.3594,1.56167,0.91363,0.4838,0.30066,0.11821},
477         {2065,2.45034,3.35558,3.51224,3.22837,    477         {2065,2.45034,3.35558,3.51224,3.22837,2.37312,1.57403,0.92377,0.48995,0.30366,0.11957},
478         {2070,2.45152,3.35816,3.51725,3.24084,    478         {2070,2.45152,3.35816,3.51725,3.24084,2.38675,1.58631,0.93373,0.49604,0.30655,0.12087},
479         {2075,2.4527,3.36053,3.52245,3.25325,2    479         {2075,2.4527,3.36053,3.52245,3.25325,2.40031,1.59853,0.94349,0.50209,0.30933,0.12209},
480         {2080,2.45388,3.36265,3.52782,3.26559,    480         {2080,2.45388,3.36265,3.52782,3.26559,2.41382,1.61071,0.95306,0.50811,0.31202,0.12324},
481         {2085,2.45506,3.36448,3.53335,3.27788,    481         {2085,2.45506,3.36448,3.53335,3.27788,2.4273,1.62286,0.96247,0.51413,0.31462,0.12433},
482         {2090,2.45626,3.36599,3.53911,3.29012,    482         {2090,2.45626,3.36599,3.53911,3.29012,2.44076,1.63501,0.97171,0.52015,0.31713,0.12534},
483         {2095,2.45747,3.36713,3.54518,3.30232,    483         {2095,2.45747,3.36713,3.54518,3.30232,2.4542,1.64714,0.98078,0.52617,0.31957,0.12627},
484         {2100,2.4587,3.36785,3.5517,3.31448,2.    484         {2100,2.4587,3.36785,3.5517,3.31448,2.46763,1.65927,0.98971,0.53219,0.32192,0.12714},
485         {2105,2.45995,3.3681,3.55882,3.32661,2    485         {2105,2.45995,3.3681,3.55882,3.32661,2.48105,1.67138,0.99848,0.5382,0.32419,0.12793},
486         {2110,2.46121,3.36785,3.56668,3.3387,2    486         {2110,2.46121,3.36785,3.56668,3.3387,2.49446,1.68349,1.1971,0.5442,0.32639,0.12864},
487         {2115,2.46248,3.36709,3.57531,3.35077,    487         {2115,2.46248,3.36709,3.57531,3.35077,2.50787,1.69558,1.01559,0.55019,0.32851,0.12928},
488         {2120,2.46378,3.36582,3.58467,3.3628,2    488         {2120,2.46378,3.36582,3.58467,3.3628,2.52126,1.70766,1.02394,0.55616,0.33056,0.12983},
489         {2125,2.4651,3.36413,3.59451,3.37482,2    489         {2125,2.4651,3.36413,3.59451,3.37482,2.53466,1.71973,1.03217,0.56211,0.33254,0.13031},
490         {2130,2.46643,3.36214,3.60448,3.38682,    490         {2130,2.46643,3.36214,3.60448,3.38682,2.54806,1.73178,1.04027,0.56805,0.33445,0.1307},
491         {2135,2.46779,3.36001,3.61413,3.39879,    491         {2135,2.46779,3.36001,3.61413,3.39879,2.56146,1.74382,1.04825,0.57395,0.33629,0.13101},
492         {2140,2.46917,3.35797,3.62309,3.41076,    492         {2140,2.46917,3.35797,3.62309,3.41076,2.57486,1.75585,1.05612,0.57983,0.33807,0.13124},
493         {2145,2.47056,3.35605,3.63094,3.42271,    493         {2145,2.47056,3.35605,3.63094,3.42271,2.58827,1.76786,1.06388,0.58568,0.33978,0.13138},
494         {2150,2.47199,3.35441,3.63757,3.43466,    494         {2150,2.47199,3.35441,3.63757,3.43466,2.60169,1.77986,1.07154,0.59149,0.34143,0.13143},
495         {2155,2.47343,3.3531,3.64301,3.4466,2.    495         {2155,2.47343,3.3531,3.64301,3.4466,2.61512,1.79184,1.0791,0.59726,0.34301,0.1314},
496         {2160,2.4749,3.35214,3.64737,3.45854,2    496         {2160,2.4749,3.35214,3.64737,3.45854,2.62857,1.8038,1.08657,0.60299,0.34454,0.13127},
497         {2165,2.47639,3.35154,3.65085,3.47048,    497         {2165,2.47639,3.35154,3.65085,3.47048,2.64203,1.81575,1.09396,0.60868,0.34602,0.13106},
498         {2170,2.47791,3.3513,3.65368,3.48243,2    498         {2170,2.47791,3.3513,3.65368,3.48243,2.65551,1.82768,1.10126,0.61432,0.34743,0.13075},
499         {2175,2.47945,3.35147,3.65608,3.49439,    499         {2175,2.47945,3.35147,3.65608,3.49439,2.66901,1.83959,1.10849,0.61991,0.3488,0.13035},
500         {2180,2.48102,3.35223,3.65849,3.50635,    500         {2180,2.48102,3.35223,3.65849,3.50635,2.68254,1.85148,1.11565,0.62544,0.35011,0.12985},
501         {2185,2.48262,3.35353,3.66096,3.51833,    501         {2185,2.48262,3.35353,3.66096,3.51833,2.69609,1.86335,1.12274,0.63092,0.35137,0.12925},
502         {2190,2.48425,3.35566,3.66378,3.53033,    502         {2190,2.48425,3.35566,3.66378,3.53033,2.70967,1.8752,1.12977,0.63634,0.35259,0.12856},
503         {2195,2.4859,3.3589,3.66741,3.54235,2.    503         {2195,2.4859,3.3589,3.66741,3.54235,2.72328,1.88704,1.13675,0.64169,0.35375,0.12777},
504         {2200,2.48758,3.36365,3.6723,3.55439,2    504         {2200,2.48758,3.36365,3.6723,3.55439,2.73693,1.89885,1.14368,0.64697,0.35488,0.12688},
505         {2205,2.4893,3.3704,3.67895,3.56646,2.    505         {2205,2.4893,3.3704,3.67895,3.56646,2.7506,1.91064,1.15057,0.65219,0.35596,0.12589},
506         {2210,2.49104,3.37963,3.68794,3.57856,    506         {2210,2.49104,3.37963,3.68794,3.57856,2.76432,1.9224,1.15742,0.65733,0.357,0.12479},
507         {2215,2.49282,3.39172,3.69969,3.59069,    507         {2215,2.49282,3.39172,3.69969,3.59069,2.77807,1.93415,1.16424,0.6624,0.358,0.1236},
508         {2220,2.49463,3.40674,3.7144,3.60285,2    508         {2220,2.49463,3.40674,3.7144,3.60285,2.79187,1.94587,1.17103,0.66738,0.35896,0.12229},
509         {2225,2.49647,3.42421,3.73174,3.61506,    509         {2225,2.49647,3.42421,3.73174,3.61506,2.80571,1.95756,1.17779,0.67228,0.35989,0.12088},
510         {2230,2.49835,3.44326,3.75104,3.62731,    510         {2230,2.49835,3.44326,3.75104,3.62731,2.8196,1.96923,1.18454,0.6771,0.36078,0.11936},
511         {2235,2.50025,3.46231,3.77088,3.6396,2    511         {2235,2.50025,3.46231,3.77088,3.6396,2.83354,1.98088,1.19128,0.68183,0.36165,0.11773},
512         {2240,2.5022,3.48011,3.78998,3.65194,2    512         {2240,2.5022,3.48011,3.78998,3.65194,2.84753,1.9925,1.19802,0.68646,0.36248,0.11598},
513         {2245,2.50418,3.49581,3.80762,3.66433,    513         {2245,2.50418,3.49581,3.80762,3.66433,2.86157,2.1941,1.20475,0.691,0.36328,0.11413},
514         {2250,2.5062,3.50913,3.82344,3.67678,2    514         {2250,2.5062,3.50913,3.82344,3.67678,2.87567,2.01566,1.21149,0.69544,0.36406,0.11216},
515         {2255,2.50825,3.52022,3.83754,3.68929,    515         {2255,2.50825,3.52022,3.83754,3.68929,2.88983,2.0272,1.21823,0.69977,0.36481,0.11008},
516         {2260,2.51034,3.52949,3.85027,3.70185,    516         {2260,2.51034,3.52949,3.85027,3.70185,2.90405,2.03872,1.225,0.704,0.36554,0.10788},
517         {2265,2.51247,3.53741,3.86203,3.71448,    517         {2265,2.51247,3.53741,3.86203,3.71448,2.91833,2.0502,1.23178,0.70812,0.36624,0.10556},
518         {2270,2.51464,3.54443,3.87322,3.72718,    518         {2270,2.51464,3.54443,3.87322,3.72718,2.93268,2.06165,1.23859,0.71213,0.36693,0.10313},
519         {2275,2.51685,3.5509,3.88415,3.73994,2    519         {2275,2.51685,3.5509,3.88415,3.73994,2.94709,2.07308,1.24543,0.71603,0.3676,0.10057},
520         {2280,2.5191,3.55707,3.89503,3.75278,2    520         {2280,2.5191,3.55707,3.89503,3.75278,2.96158,2.08447,1.25231,0.71981,0.36826,0.09789},
521         {2285,2.52139,3.56312,3.90596,3.7657,2    521         {2285,2.52139,3.56312,3.90596,3.7657,2.97614,2.09584,1.25923,0.72346,0.3689,0.09509},
522         {2290,2.52372,3.56913,3.91699,3.77869,    522         {2290,2.52372,3.56913,3.91699,3.77869,2.99077,2.10717,1.26619,0.72699,0.36952,0.09217},
523         {2295,2.52609,3.57517,3.92809,3.79175,    523         {2295,2.52609,3.57517,3.92809,3.79175,3.00546,2.11846,1.27319,0.73041,0.37014,0.08914},
524         {2300,2.5285,3.58122,3.93914,3.80486,3    524         {2300,2.5285,3.58122,3.93914,3.80486,3.0202,2.12972,1.28021,0.73371,0.37074,0.08599},
525         {2305,2.53093,3.58725,3.95002,3.81803,    525         {2305,2.53093,3.58725,3.95002,3.81803,3.03497,2.14093,1.28726,0.73689,0.37132,0.08276},
526         {2310,2.5334,3.5932,3.96055,3.83123,3.    526         {2310,2.5334,3.5932,3.96055,3.83123,3.04978,2.15209,1.29431,0.73998,0.37189,0.07943},
527         {2315,2.53589,3.599,3.97054,3.84446,3.    527         {2315,2.53589,3.599,3.97054,3.84446,3.06461,2.16319,1.30137,0.74296,0.37245,0.07603},
528         {2320,2.5384,3.60457,3.97979,3.85771,3    528         {2320,2.5384,3.60457,3.97979,3.85771,3.07944,2.17424,1.30842,0.74584,0.37299,0.07255},
529         {2325,2.54093,3.60985,3.98817,3.87097,    529         {2325,2.54093,3.60985,3.98817,3.87097,3.09426,2.18523,1.31545,0.74863,0.37351,0.06901},
530         {2330,2.54348,3.61467,3.99541,3.88423,    530         {2330,2.54348,3.61467,3.99541,3.88423,3.10907,2.19615,1.32246,0.75133,0.37402,0.06542},
531         {2335,2.54603,3.61904,4.00148,3.89749,    531         {2335,2.54603,3.61904,4.00148,3.89749,3.12385,2.20701,1.32944,0.75395,0.3745,0.06178},
532         {2340,2.5486,3.6229,4.00632,3.91072,3.    532         {2340,2.5486,3.6229,4.00632,3.91072,3.1386,2.21778,1.33638,0.75649,0.37498,0.0581},
533         {2345,2.55117,3.62622,4.00994,3.92392,    533         {2345,2.55117,3.62622,4.00994,3.92392,3.15329,2.22848,1.34326,0.75895,0.37543,0.0544},
534         {2350,2.55375,3.62901,4.01239,3.93709,    534         {2350,2.55375,3.62901,4.01239,3.93709,3.16793,2.2391,1.35009,0.76134,0.37586,0.05068},
535         {2355,2.55632,3.6313,4.01379,3.9502,3.    535         {2355,2.55632,3.6313,4.01379,3.9502,3.1825,2.24963,1.35684,0.76366,0.37628,0.04694},
536         {2360,2.5589,3.6331,4.01428,3.96326,3.    536         {2360,2.5589,3.6331,4.01428,3.96326,3.19698,2.26006,1.36353,0.76592,0.37668,0.0432},
537         {2365,2.56146,3.63449,4.01403,3.97626,    537         {2365,2.56146,3.63449,4.01403,3.97626,3.21137,2.2704,1.37012,0.76813,0.37705,0.03947},
538         {2370,2.56401,3.63552,4.0132,3.98917,3    538         {2370,2.56401,3.63552,4.0132,3.98917,3.22566,2.28065,1.37663,0.77027,0.37741,0.03576},
539         {2375,2.56655,3.63624,4.01196,4.002,3.    539         {2375,2.56655,3.63624,4.01196,4.002,3.23983,2.29079,1.38303,0.77237,0.37774,0.03206}};
540                                                   540       
541       const G4int coef_ener = G4int((pLab-Lege    541       const G4int coef_ener = G4int((pLab-Legendre_coef[0][0])/5);
542       const G4double sup_ener = pLab/5. - coef    542       const G4double sup_ener = pLab/5. - coef_ener -Legendre_coef[0][0]/5;
543       //std::cout << "sup_ener\t" << sup_ener     543       //std::cout << "sup_ener\t" << sup_ener << std::endl;
544                                                   544       
545 // assert(pLab >= Legendre_coef[coef_ener][0]     545 // assert(pLab >= Legendre_coef[coef_ener][0] && pLab < Legendre_coef[coef_ener+1][0]);
546                                                   546       
547       // Legendre coefficient normalized          547       // Legendre coefficient normalized
548       const G4double A0 = 1.;                     548       const G4double A0 = 1.;
549       const G4double A1 = (1-sup_ener)*Legendr    549       const G4double A1 = (1-sup_ener)*Legendre_coef[coef_ener][1] + sup_ener*Legendre_coef[coef_ener+1][1];
550       const G4double A2 = (1-sup_ener)*Legendr    550       const G4double A2 = (1-sup_ener)*Legendre_coef[coef_ener][2] + sup_ener*Legendre_coef[coef_ener+1][2];
551       const G4double A3 = (1-sup_ener)*Legendr    551       const G4double A3 = (1-sup_ener)*Legendre_coef[coef_ener][3] + sup_ener*Legendre_coef[coef_ener+1][3];
552       const G4double A4 = (1-sup_ener)*Legendr    552       const G4double A4 = (1-sup_ener)*Legendre_coef[coef_ener][4] + sup_ener*Legendre_coef[coef_ener+1][4];
553       const G4double A5 = (1-sup_ener)*Legendr    553       const G4double A5 = (1-sup_ener)*Legendre_coef[coef_ener][5] + sup_ener*Legendre_coef[coef_ener+1][5];
554       const G4double A6 = (1-sup_ener)*Legendr    554       const G4double A6 = (1-sup_ener)*Legendre_coef[coef_ener][6] + sup_ener*Legendre_coef[coef_ener+1][6];
555       const G4double A7 = (1-sup_ener)*Legendr    555       const G4double A7 = (1-sup_ener)*Legendre_coef[coef_ener][7] + sup_ener*Legendre_coef[coef_ener+1][7];
556       const G4double A8 = (1-sup_ener)*Legendr    556       const G4double A8 = (1-sup_ener)*Legendre_coef[coef_ener][8] + sup_ener*Legendre_coef[coef_ener+1][8];
557       const G4double A9 = (1-sup_ener)*Legendr    557       const G4double A9 = (1-sup_ener)*Legendre_coef[coef_ener][9] + sup_ener*Legendre_coef[coef_ener+1][9];
558       const G4double A10 = (1-sup_ener)*Legend    558       const G4double A10 = (1-sup_ener)*Legendre_coef[coef_ener][10] + sup_ener*Legendre_coef[coef_ener+1][10];
559                                                   559       
560       // Theoritical max if all Ai > 0 (often     560       // Theoritical max if all Ai > 0 (often the case)
561       const G4double A = std::fabs(A0) + std::    561       const G4double A = std::fabs(A0) + std::fabs(A1) + std::fabs(A2) + std::fabs(A3) + std::fabs(A4) + std::fabs(A5) + std::fabs(A6) + std::fabs(A7) + std::fabs(A8) + std::fabs(A9) + std::fabs(A10);
562                                                   562       
563       G4bool success = false;                     563       G4bool success = false;
564       G4int maxloop = 0;                          564       G4int maxloop = 0;
565                                                   565       
566       while(!success && maxloop < 1000){          566       while(!success && maxloop < 1000){
567                                                   567         
568         cos_theta = Random::shoot()*2-1.; // n    568         cos_theta = Random::shoot()*2-1.; // not optimized
569                                                   569         
570         // Legendre Polynomial                    570         // Legendre Polynomial
571         G4double P0 = A0;                         571         G4double P0 = A0;
572         G4double P1 = A1*cos_theta;               572         G4double P1 = A1*cos_theta;
573         G4double P2 = A2/2.*(3*std::pow(cos_th    573         G4double P2 = A2/2.*(3*std::pow(cos_theta,2)-1);
574         G4double P3 = A3/2.*(5*std::pow(cos_th    574         G4double P3 = A3/2.*(5*std::pow(cos_theta,3)-3*cos_theta);
575         G4double P4 = A4/8.*(35*std::pow(cos_t    575         G4double P4 = A4/8.*(35*std::pow(cos_theta,4)-30*std::pow(cos_theta,2)+3);
576         G4double P5 = A5/8.*(63*std::pow(cos_t    576         G4double P5 = A5/8.*(63*std::pow(cos_theta,5)-70*std::pow(cos_theta,3)+15*cos_theta);
577         G4double P6 = A6/16.*(231*std::pow(cos    577         G4double P6 = A6/16.*(231*std::pow(cos_theta,6)-315*std::pow(cos_theta,4)+105*std::pow(cos_theta,2)-5);
578         G4double P7 = A7/16.*(429*std::pow(cos    578         G4double P7 = A7/16.*(429*std::pow(cos_theta,7)-693*std::pow(cos_theta,5)+315*std::pow(cos_theta,3)-35*cos_theta);
579         G4double P8 = A8/128.*(6435*std::pow(c    579         G4double P8 = A8/128.*(6435*std::pow(cos_theta,8)-12012*std::pow(cos_theta,6)+6930*std::pow(cos_theta,4)-1260*std::pow(cos_theta,2)+35);
580         G4double P9 = A9/128.*(12155*std::pow(    580         G4double P9 = A9/128.*(12155*std::pow(cos_theta,9)-25740*std::pow(cos_theta,7)+18018*std::pow(cos_theta,5)-4620*std::pow(cos_theta,3)+315*cos_theta);
581         G4double P10 = A10/256.*(46189*std::po    581         G4double P10 = A10/256.*(46189*std::pow(cos_theta,10)-109395*std::pow(cos_theta,8)+90090*std::pow(cos_theta,6)-30030*std::pow(cos_theta,4)+3465*std::pow(cos_theta,2)-63);
582                                                   582         
583         G4double P = (P0 + P1 + P2 + P3 + P4 +    583         G4double P = (P0 + P1 + P2 + P3 + P4 + P5 + P6 + P7 + P8 + P9 + P10);
584                                                   584         
585         if(Random::shoot()*A < P) success = tr    585         if(Random::shoot()*A < P) success = true;
586         maxloop +=1 ;                             586         maxloop +=1 ;
587         if(maxloop==1000) cos_theta = std::log    587         if(maxloop==1000) cos_theta = std::log(Random::shoot()*(std::exp(10.)-std::exp(-10.))+std::exp(-10.))/10.; // if no success in 1E4 shoot, probably angulard distribution piked very forward 
588       }                                           588       }
589       sin_theta = std::sqrt(1-cos_theta*cos_th    589       sin_theta = std::sqrt(1-cos_theta*cos_theta);
590     }                                             590     }
591                                                   591     
592     if(rho == 0) return ThreeVector(sin_theta*    592     if(rho == 0) return ThreeVector(sin_theta*cos_phi,sin_theta*sin_phi,cos_theta)*norm;
593     // Rotation in the direction of the incide    593     // Rotation in the direction of the incident kaon
594     const G4double px = x/r*cos_theta - y/rho*    594     const G4double px = x/r*cos_theta - y/rho*sin_theta*cos_phi + z/r*x/rho*sin_theta*sin_phi;
595     const G4double py = y/r*cos_theta + x/rho*    595     const G4double py = y/r*cos_theta + x/rho*sin_theta*cos_phi + z/r*y/rho*sin_theta*sin_phi;
596     const G4double pz = z/r*cos_theta - rho/r*    596     const G4double pz = z/r*cos_theta - rho/r*sin_theta*sin_phi;
597                                                   597     
598                                                   598     
599     return ThreeVector(px,py,pz)*norm;            599     return ThreeVector(px,py,pz)*norm;
600   }                                               600   }
601 }                                                 601 }
602                                                   602