Add momentum correction workaround
authoreserradi <eserradi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 27 Jun 2013 13:03:00 +0000 (13:03 +0000)
committereserradi <eserradi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 27 Jun 2013 13:03:00 +0000 (13:03 +0000)
PWGLF/SPECTRA/Nuclei/B2/AliAnalysisTaskB2.cxx
PWGLF/SPECTRA/Nuclei/B2/AliAnalysisTaskB2.h
PWGLF/SPECTRA/Nuclei/B2/macros/AddTaskB2.C
PWGLF/SPECTRA/Nuclei/B2/macros/MomentumCorrection.C [new file with mode: 0644]

index 66864e2..74b5ec8 100644 (file)
@@ -49,6 +49,7 @@
 #include <TH1D.h>
 #include <TH2D.h>
 #include <TList.h>
+#include <TProfile.h>
 
 #include "AliLnID.h"
 #include "AliLnHistoMap.h"
@@ -109,6 +110,8 @@ AliAnalysisTaskB2::AliAnalysisTaskB2()
 , fTOFmatch(0)
 , fMinM2(2.)
 , fMaxM2(6.)
+, fMomentumCorrection(0)
+, fMoCpfx(0)
 
 {
 //
@@ -172,6 +175,8 @@ AliAnalysisTaskB2::AliAnalysisTaskB2(const char* name)
 , fTOFmatch(0)
 , fMinM2(2.)
 , fMaxM2(6.)
+, fMomentumCorrection(0)
+, fMoCpfx(0)
 
 {
 //
@@ -274,6 +279,8 @@ AliAnalysisTaskB2::~AliAnalysisTaskB2()
        delete fTrigAna;
        
        if(fIsPidOwner) delete fESDpid;
+       
+       delete fMoCpfx;
 }
 
 void AliAnalysisTaskB2::SetParticleSpecies(const TString& species)
@@ -630,6 +637,7 @@ Int_t AliAnalysisTaskB2::GetTracks()
                
                // momentum
                
+               Double_t p       = iTrack->GetP()*z; // p at DCA
                Double_t pt      = iTrack->Pt()*z; // pt at DCA
                Double_t y       = this->GetRapidity(iTrack, fPartCode);
                Double_t pITS    = this->GetITSmomentum(iTrack);
@@ -640,6 +648,13 @@ Int_t AliAnalysisTaskB2::GetTracks()
                Int_t nPointsITS = this->GetITSnPointsPID(iTrack);
                Int_t nPointsTPC = iTrack->GetTPCsignalN();
                
+               if(fMomentumCorrection)
+               {
+                       pt += this->GetMomentumCorrection(pt);
+                       p   = TMath::Sqrt(pt*pt + iTrack->Pz()*iTrack->Pz());
+                       y   = this->GetRapidity(p, iTrack->Pz(), fPartCode);
+               }
+               
                Double_t beta = 0;
                Double_t mass = 0;
                Double_t m2   = 0;
@@ -1264,3 +1279,26 @@ Int_t AliAnalysisTaskB2::GetPidCode(const TString& species) const
        
        return -1;
 }
+
+Double_t AliAnalysisTaskB2::GetMomentumCorrection(Double_t ptrec) const
+{
+//
+// momentum correction for low pt
+//
+       if(fMoCpfx == 0) return 0;
+       
+       return fMoCpfx->Interpolate(ptrec);
+}
+
+Double_t AliAnalysisTaskB2::GetRapidity(Double_t p, Double_t pz, Int_t pid) const
+{
+//
+// Rapidity (for momentum correction)
+//
+       Double_t m  = AliPID::ParticleMass(pid);
+       Double_t e  = TMath::Sqrt(p*p + m*m);
+       
+       if(e <= pz) return 1.e+16;
+       
+       return 0.5*TMath::Log( (e+pz)/(e-pz) );
+}
index 3f6ab0d..7fd33b5 100644 (file)
@@ -19,6 +19,7 @@ class AliLnHistoMap;
 class AliLnID;
 class TParticle;
 class TList;
+class TProfile;
 
 class AliAnalysisTaskB2: public AliAnalysisTask
 {
@@ -55,6 +56,9 @@ class AliAnalysisTaskB2: public AliAnalysisTask
        void SetSimulation(Bool_t flag = kTRUE) { fSimulation = flag; }
        void SetHeavyIons(Bool_t flag = kTRUE) { fHeavyIons = flag; }
        
+       void SetMomentumCorrection(Bool_t flag = kTRUE) { fMomentumCorrection = flag; }
+       void SetMomentumCorrectionProfile(TProfile* pfx) { fMoCpfx = pfx; }
+       
        void SetHistogramMap(AliLnHistoMap* map) { fHistoMap = map; }
        
        void SetESDtrackCuts(AliESDtrackCuts* esdTrackCuts) {fESDtrackCuts = esdTrackCuts; }
@@ -101,6 +105,7 @@ class AliAnalysisTaskB2: public AliAnalysisTask
        Double_t GetPhi(const AliESDtrack* trk) const;
        Double_t GetTheta(const AliESDtrack* trk) const;
        Double_t GetRapidity(const AliESDtrack* trk, Int_t pid) const;
+       Double_t GetRapidity(Double_t p, Double_t pz, Int_t pid) const;
        Double_t GetITSmomentum(const AliESDtrack* trk) const;
        Double_t GetTOFmomentum(const AliESDtrack* trk) const;
        Double_t GetBeta(const AliESDtrack* trk) const;
@@ -112,6 +117,8 @@ class AliAnalysisTaskB2: public AliAnalysisTask
        
        Int_t GetPidCode(const TString& species) const;
        
+       Double_t GetMomentumCorrection(Double_t ptrec) const;
+       
        Double_t GetDiffM2(Double_t beta, Double_t p, Double_t m=1.875612793) const;
        
   private:
@@ -181,6 +188,9 @@ class AliAnalysisTaskB2: public AliAnalysisTask
        Double_t fMinM2; // minimum m2 for TPC+TOF pid
        Double_t fMaxM2; // maximum m2 for TPC+TOF pid
        
+       Bool_t fMomentumCorrection; // enable momentum correction
+       TProfile* fMoCpfx; // momentum correction from simulation
+       
        ClassDef(AliAnalysisTaskB2, 1)
 };
 
index c28e9de..6852bf5 100644 (file)
@@ -36,6 +36,7 @@ AliAnalysisTaskB2* AddTaskB2(  const TString& species
                              , Bool_t   V0AND            = kFALSE
                              , const TString& ztag       = ""
                              , Double_t maxVz            = 10
+                             , Bool_t momentumCorr       = kFALSE
                              , const TString& binSize    = ""
                              , Double_t minCentrality    = 0
                              , Double_t maxCentrality    = 20
@@ -111,6 +112,17 @@ AliAnalysisTaskB2* AddTaskB2(  const TString& species
                task->SetNoFastOnlyTrigger();
        }
        
+       // momentum correction
+       
+       if(momentumCorr && species=="Deuteron")
+       {
+               gROOT->LoadMacro("$ALICE_ROOT/PWGLF/SPECTRA/Nuclei/B2/macros/MomentumCorrection.C");
+               
+               TProfile* pfx = MomentumCorrection(species);
+               task->SetMomentumCorrectionProfile(pfx);
+               if(pfx != 0) task->SetMomentumCorrection();
+       }
+       
        // histograms
        
        gROOT->LoadMacro("$ALICE_ROOT/PWGLF/SPECTRA/Nuclei/B2/macros/CreateHistograms.C");
diff --git a/PWGLF/SPECTRA/Nuclei/B2/macros/MomentumCorrection.C b/PWGLF/SPECTRA/Nuclei/B2/macros/MomentumCorrection.C
new file mode 100644 (file)
index 0000000..4d15383
--- /dev/null
@@ -0,0 +1,1212 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+// momentum correction for deuteron and antideuterons
+// author: Eulogio Serradilla <eulogio.serradilla@cern.ch>
+
+TProfile* MomentumCorrection(const TString& species)
+{
+//
+// momentum correction for deuterons and antideuterons (LHC12a5b simulation)
+//
+       using namespace std;
+       
+       if(species != "Deuteron")
+       {
+               cerr << "MomentumCorrection: Unknown particle species " << species << endl;
+               return 0;
+       }
+       
+       TProfile* pfx = new TProfile("pfx","momentum correction",400,0,10,"");
+       pfx->SetBinEntries(10,6);
+       pfx->SetBinEntries(11,56);
+       pfx->SetBinEntries(12,99);
+       pfx->SetBinEntries(13,258);
+       pfx->SetBinEntries(14,518);
+       pfx->SetBinEntries(15,885);
+       pfx->SetBinEntries(16,1487);
+       pfx->SetBinEntries(17,2314);
+       pfx->SetBinEntries(18,2981);
+       pfx->SetBinEntries(19,4002);
+       pfx->SetBinEntries(20,5128);
+       pfx->SetBinEntries(21,6384);
+       pfx->SetBinEntries(22,7417);
+       pfx->SetBinEntries(23,8957);
+       pfx->SetBinEntries(24,10592);
+       pfx->SetBinEntries(25,11781);
+       pfx->SetBinEntries(26,13228);
+       pfx->SetBinEntries(27,14760);
+       pfx->SetBinEntries(28,16367);
+       pfx->SetBinEntries(29,17698);
+       pfx->SetBinEntries(30,18815);
+       pfx->SetBinEntries(31,20419);
+       pfx->SetBinEntries(32,21564);
+       pfx->SetBinEntries(33,22829);
+       pfx->SetBinEntries(34,23987);
+       pfx->SetBinEntries(35,25276);
+       pfx->SetBinEntries(36,26546);
+       pfx->SetBinEntries(37,27766);
+       pfx->SetBinEntries(38,28677);
+       pfx->SetBinEntries(39,29567);
+       pfx->SetBinEntries(40,30662);
+       pfx->SetBinEntries(41,31546);
+       pfx->SetBinEntries(42,32495);
+       pfx->SetBinEntries(43,33434);
+       pfx->SetBinEntries(44,33888);
+       pfx->SetBinEntries(45,34615);
+       pfx->SetBinEntries(46,35279);
+       pfx->SetBinEntries(47,35948);
+       pfx->SetBinEntries(48,36492);
+       pfx->SetBinEntries(49,37220);
+       pfx->SetBinEntries(50,38035);
+       pfx->SetBinEntries(51,38550);
+       pfx->SetBinEntries(52,39359);
+       pfx->SetBinEntries(53,40377);
+       pfx->SetBinEntries(54,40239);
+       pfx->SetBinEntries(55,40998);
+       pfx->SetBinEntries(56,41601);
+       pfx->SetBinEntries(57,41776);
+       pfx->SetBinEntries(58,42252);
+       pfx->SetBinEntries(59,42120);
+       pfx->SetBinEntries(60,42327);
+       pfx->SetBinEntries(61,42512);
+       pfx->SetBinEntries(62,42306);
+       pfx->SetBinEntries(63,42677);
+       pfx->SetBinEntries(64,42817);
+       pfx->SetBinEntries(65,42576);
+       pfx->SetBinEntries(66,42769);
+       pfx->SetBinEntries(67,42717);
+       pfx->SetBinEntries(68,43400);
+       pfx->SetBinEntries(69,43141);
+       pfx->SetBinEntries(70,43527);
+       pfx->SetBinEntries(71,43461);
+       pfx->SetBinEntries(72,43434);
+       pfx->SetBinEntries(73,43224);
+       pfx->SetBinEntries(74,43466);
+       pfx->SetBinEntries(75,43606);
+       pfx->SetBinEntries(76,43434);
+       pfx->SetBinEntries(77,43603);
+       pfx->SetBinEntries(78,43660);
+       pfx->SetBinEntries(79,43810);
+       pfx->SetBinEntries(80,43549);
+       pfx->SetBinEntries(81,43492);
+       pfx->SetBinEntries(82,43530);
+       pfx->SetBinEntries(83,43265);
+       pfx->SetBinEntries(84,43624);
+       pfx->SetBinEntries(85,43369);
+       pfx->SetBinEntries(86,43503);
+       pfx->SetBinEntries(87,43358);
+       pfx->SetBinEntries(88,43411);
+       pfx->SetBinEntries(89,43195);
+       pfx->SetBinEntries(90,43088);
+       pfx->SetBinEntries(91,43132);
+       pfx->SetBinEntries(92,42967);
+       pfx->SetBinEntries(93,43226);
+       pfx->SetBinEntries(94,43061);
+       pfx->SetBinEntries(95,43261);
+       pfx->SetBinEntries(96,42624);
+       pfx->SetBinEntries(97,42853);
+       pfx->SetBinEntries(98,42695);
+       pfx->SetBinEntries(99,42825);
+       pfx->SetBinEntries(100,42781);
+       pfx->SetBinEntries(101,42906);
+       pfx->SetBinEntries(102,42984);
+       pfx->SetBinEntries(103,42796);
+       pfx->SetBinEntries(104,42680);
+       pfx->SetBinEntries(105,42722);
+       pfx->SetBinEntries(106,42943);
+       pfx->SetBinEntries(107,42431);
+       pfx->SetBinEntries(108,42669);
+       pfx->SetBinEntries(109,42863);
+       pfx->SetBinEntries(110,42595);
+       pfx->SetBinEntries(111,42599);
+       pfx->SetBinEntries(112,42699);
+       pfx->SetBinEntries(113,42453);
+       pfx->SetBinEntries(114,42357);
+       pfx->SetBinEntries(115,42461);
+       pfx->SetBinEntries(116,42192);
+       pfx->SetBinEntries(117,42496);
+       pfx->SetBinEntries(118,42060);
+       pfx->SetBinEntries(119,42380);
+       pfx->SetBinEntries(120,42336);
+       pfx->SetBinEntries(121,41972);
+       pfx->SetBinEntries(122,41957);
+       pfx->SetBinEntries(123,42088);
+       pfx->SetBinEntries(124,42319);
+       pfx->SetBinEntries(125,42399);
+       pfx->SetBinEntries(126,42173);
+       pfx->SetBinEntries(127,41839);
+       pfx->SetBinEntries(128,42391);
+       pfx->SetBinEntries(129,42361);
+       pfx->SetBinEntries(130,42489);
+       pfx->SetBinEntries(131,42260);
+       pfx->SetBinEntries(132,41975);
+       pfx->SetBinEntries(133,41939);
+       pfx->SetBinEntries(134,42103);
+       pfx->SetBinEntries(135,41957);
+       pfx->SetBinEntries(136,41993);
+       pfx->SetBinEntries(137,41937);
+       pfx->SetBinEntries(138,42338);
+       pfx->SetBinEntries(139,42096);
+       pfx->SetBinEntries(140,42374);
+       pfx->SetBinEntries(141,42199);
+       pfx->SetBinEntries(142,42295);
+       pfx->SetBinEntries(143,42371);
+       pfx->SetBinEntries(144,42102);
+       pfx->SetBinEntries(145,42288);
+       pfx->SetBinEntries(146,41981);
+       pfx->SetBinEntries(147,42365);
+       pfx->SetBinEntries(148,42286);
+       pfx->SetBinEntries(149,42063);
+       pfx->SetBinEntries(150,41996);
+       pfx->SetBinEntries(151,42186);
+       pfx->SetBinEntries(152,42308);
+       pfx->SetBinEntries(153,42263);
+       pfx->SetBinEntries(154,42426);
+       pfx->SetBinEntries(155,42310);
+       pfx->SetBinEntries(156,41989);
+       pfx->SetBinEntries(157,42192);
+       pfx->SetBinEntries(158,41975);
+       pfx->SetBinEntries(159,42576);
+       pfx->SetBinEntries(160,42025);
+       pfx->SetBinEntries(161,42134);
+       pfx->SetBinEntries(162,42420);
+       pfx->SetBinEntries(163,42706);
+       pfx->SetBinEntries(164,42285);
+       pfx->SetBinEntries(165,42439);
+       pfx->SetBinEntries(166,42215);
+       pfx->SetBinEntries(167,42272);
+       pfx->SetBinEntries(168,42553);
+       pfx->SetBinEntries(169,42311);
+       pfx->SetBinEntries(170,41971);
+       pfx->SetBinEntries(171,42428);
+       pfx->SetBinEntries(172,42489);
+       pfx->SetBinEntries(173,42112);
+       pfx->SetBinEntries(174,42354);
+       pfx->SetBinEntries(175,42196);
+       pfx->SetBinEntries(176,41917);
+       pfx->SetBinEntries(177,41915);
+       pfx->SetBinEntries(178,42038);
+       pfx->SetBinEntries(179,42298);
+       pfx->SetBinEntries(180,42501);
+       pfx->SetBinEntries(181,42404);
+       pfx->SetBinEntries(182,42097);
+       pfx->SetBinEntries(183,42333);
+       pfx->SetBinEntries(184,42418);
+       pfx->SetBinEntries(185,42186);
+       pfx->SetBinEntries(186,42166);
+       pfx->SetBinEntries(187,42193);
+       pfx->SetBinEntries(188,42315);
+       pfx->SetBinEntries(189,42026);
+       pfx->SetBinEntries(190,42564);
+       pfx->SetBinEntries(191,42396);
+       pfx->SetBinEntries(192,42526);
+       pfx->SetBinEntries(193,41923);
+       pfx->SetBinEntries(194,42102);
+       pfx->SetBinEntries(195,42557);
+       pfx->SetBinEntries(196,42360);
+       pfx->SetBinEntries(197,42106);
+       pfx->SetBinEntries(198,42532);
+       pfx->SetBinEntries(199,42153);
+       pfx->SetBinEntries(200,42359);
+       pfx->SetBinEntries(201,42290);
+       pfx->SetBinEntries(202,42284);
+       pfx->SetBinEntries(203,42198);
+       pfx->SetBinEntries(204,42354);
+       pfx->SetBinEntries(205,42176);
+       pfx->SetBinEntries(206,42424);
+       pfx->SetBinEntries(207,41795);
+       pfx->SetBinEntries(208,42197);
+       pfx->SetBinEntries(209,42619);
+       pfx->SetBinEntries(210,42181);
+       pfx->SetBinEntries(211,42141);
+       pfx->SetBinEntries(212,42492);
+       pfx->SetBinEntries(213,42309);
+       pfx->SetBinEntries(214,42356);
+       pfx->SetBinEntries(215,41899);
+       pfx->SetBinEntries(216,42637);
+       pfx->SetBinEntries(217,42141);
+       pfx->SetBinEntries(218,42425);
+       pfx->SetBinEntries(219,42014);
+       pfx->SetBinEntries(220,42423);
+       pfx->SetBinEntries(221,42261);
+       pfx->SetBinEntries(222,42674);
+       pfx->SetBinEntries(223,42161);
+       pfx->SetBinEntries(224,42425);
+       pfx->SetBinEntries(225,42308);
+       pfx->SetBinEntries(226,42129);
+       pfx->SetBinEntries(227,42053);
+       pfx->SetBinEntries(228,41979);
+       pfx->SetBinEntries(229,42425);
+       pfx->SetBinEntries(230,42039);
+       pfx->SetBinEntries(231,42349);
+       pfx->SetBinEntries(232,42329);
+       pfx->SetBinEntries(233,42378);
+       pfx->SetBinEntries(234,42780);
+       pfx->SetBinEntries(235,42393);
+       pfx->SetBinEntries(236,42280);
+       pfx->SetBinEntries(237,42528);
+       pfx->SetBinEntries(238,42157);
+       pfx->SetBinEntries(239,42287);
+       pfx->SetBinEntries(240,42035);
+       pfx->SetBinEntries(241,42341);
+       pfx->SetBinEntries(242,42261);
+       pfx->SetBinEntries(243,42359);
+       pfx->SetBinEntries(244,42245);
+       pfx->SetBinEntries(245,42621);
+       pfx->SetBinEntries(246,42097);
+       pfx->SetBinEntries(247,42343);
+       pfx->SetBinEntries(248,42202);
+       pfx->SetBinEntries(249,42435);
+       pfx->SetBinEntries(250,41873);
+       pfx->SetBinEntries(251,41837);
+       pfx->SetBinEntries(252,42148);
+       pfx->SetBinEntries(253,42377);
+       pfx->SetBinEntries(254,42021);
+       pfx->SetBinEntries(255,42478);
+       pfx->SetBinEntries(256,42690);
+       pfx->SetBinEntries(257,42390);
+       pfx->SetBinEntries(258,42158);
+       pfx->SetBinEntries(259,42293);
+       pfx->SetBinEntries(260,42274);
+       pfx->SetBinEntries(261,42260);
+       pfx->SetBinEntries(262,42172);
+       pfx->SetBinEntries(263,42149);
+       pfx->SetBinEntries(264,42562);
+       pfx->SetBinEntries(265,42029);
+       pfx->SetBinEntries(266,42297);
+       pfx->SetBinEntries(267,42222);
+       pfx->SetBinEntries(268,42180);
+       pfx->SetBinEntries(269,42080);
+       pfx->SetBinEntries(270,42290);
+       pfx->SetBinEntries(271,42133);
+       pfx->SetBinEntries(272,42539);
+       pfx->SetBinEntries(273,42439);
+       pfx->SetBinEntries(274,42117);
+       pfx->SetBinEntries(275,42194);
+       pfx->SetBinEntries(276,42049);
+       pfx->SetBinEntries(277,41900);
+       pfx->SetBinEntries(278,42108);
+       pfx->SetBinEntries(279,42173);
+       pfx->SetBinEntries(280,42313);
+       pfx->SetBinEntries(281,42049);
+       pfx->SetBinEntries(282,42181);
+       pfx->SetBinEntries(283,42050);
+       pfx->SetBinEntries(284,41839);
+       pfx->SetBinEntries(285,41852);
+       pfx->SetBinEntries(286,41725);
+       pfx->SetBinEntries(287,41951);
+       pfx->SetBinEntries(288,42041);
+       pfx->SetBinEntries(289,41957);
+       pfx->SetBinEntries(290,42031);
+       pfx->SetBinEntries(291,42367);
+       pfx->SetBinEntries(292,42390);
+       pfx->SetBinEntries(293,42115);
+       pfx->SetBinEntries(294,41954);
+       pfx->SetBinEntries(295,41776);
+       pfx->SetBinEntries(296,41735);
+       pfx->SetBinEntries(297,42030);
+       pfx->SetBinEntries(298,41717);
+       pfx->SetBinEntries(299,41894);
+       pfx->SetBinEntries(300,41638);
+       pfx->SetBinEntries(301,41848);
+       pfx->SetBinEntries(302,41691);
+       pfx->SetBinEntries(303,41764);
+       pfx->SetBinEntries(304,42121);
+       pfx->SetBinEntries(305,41946);
+       pfx->SetBinEntries(306,41844);
+       pfx->SetBinEntries(307,41911);
+       pfx->SetBinEntries(308,42103);
+       pfx->SetBinEntries(309,41734);
+       pfx->SetBinEntries(310,41691);
+       pfx->SetBinEntries(311,41220);
+       pfx->SetBinEntries(312,41712);
+       pfx->SetBinEntries(313,41883);
+       pfx->SetBinEntries(314,41959);
+       pfx->SetBinEntries(315,41733);
+       pfx->SetBinEntries(316,41886);
+       pfx->SetBinEntries(317,41560);
+       pfx->SetBinEntries(318,41495);
+       pfx->SetBinEntries(319,41596);
+       pfx->SetBinEntries(320,41319);
+       pfx->SetBinEntries(321,41622);
+       pfx->SetBinEntries(322,41796);
+       pfx->SetBinEntries(323,41095);
+       pfx->SetBinEntries(324,41620);
+       pfx->SetBinEntries(325,41460);
+       pfx->SetBinEntries(326,41395);
+       pfx->SetBinEntries(327,41592);
+       pfx->SetBinEntries(328,41262);
+       pfx->SetBinEntries(329,41557);
+       pfx->SetBinEntries(330,41308);
+       pfx->SetBinEntries(331,41547);
+       pfx->SetBinEntries(332,41096);
+       pfx->SetBinEntries(333,41524);
+       pfx->SetBinEntries(334,41186);
+       pfx->SetBinEntries(335,41262);
+       pfx->SetBinEntries(336,41224);
+       pfx->SetBinEntries(337,41404);
+       pfx->SetBinEntries(338,41359);
+       pfx->SetBinEntries(339,41165);
+       pfx->SetBinEntries(340,40892);
+       pfx->SetBinEntries(341,40772);
+       pfx->SetBinEntries(342,41116);
+       pfx->SetBinEntries(343,41150);
+       pfx->SetBinEntries(344,40808);
+       pfx->SetBinEntries(345,41185);
+       pfx->SetBinEntries(346,41169);
+       pfx->SetBinEntries(347,40822);
+       pfx->SetBinEntries(348,40913);
+       pfx->SetBinEntries(349,40677);
+       pfx->SetBinEntries(350,41282);
+       pfx->SetBinEntries(351,40770);
+       pfx->SetBinEntries(352,40983);
+       pfx->SetBinEntries(353,40791);
+       pfx->SetBinEntries(354,40584);
+       pfx->SetBinEntries(355,40965);
+       pfx->SetBinEntries(356,41013);
+       pfx->SetBinEntries(357,40555);
+       pfx->SetBinEntries(358,40610);
+       pfx->SetBinEntries(359,40464);
+       pfx->SetBinEntries(360,40845);
+       pfx->SetBinEntries(361,40400);
+       pfx->SetBinEntries(362,40480);
+       pfx->SetBinEntries(363,40656);
+       pfx->SetBinEntries(364,40625);
+       pfx->SetBinEntries(365,40847);
+       pfx->SetBinEntries(366,40362);
+       pfx->SetBinEntries(367,40469);
+       pfx->SetBinEntries(368,40294);
+       pfx->SetBinEntries(369,40690);
+       pfx->SetBinEntries(370,40375);
+       pfx->SetBinEntries(371,40431);
+       pfx->SetBinEntries(372,40199);
+       pfx->SetBinEntries(373,39936);
+       pfx->SetBinEntries(374,40047);
+       pfx->SetBinEntries(375,40320);
+       pfx->SetBinEntries(376,40139);
+       pfx->SetBinEntries(377,40354);
+       pfx->SetBinEntries(378,39803);
+       pfx->SetBinEntries(379,40119);
+       pfx->SetBinEntries(380,40113);
+       pfx->SetBinEntries(381,40009);
+       pfx->SetBinEntries(382,39706);
+       pfx->SetBinEntries(383,38959);
+       pfx->SetBinEntries(384,38647);
+       pfx->SetBinEntries(385,37847);
+       pfx->SetBinEntries(386,37657);
+       pfx->SetBinEntries(387,36842);
+       pfx->SetBinEntries(388,35759);
+       pfx->SetBinEntries(389,34875);
+       pfx->SetBinEntries(390,34386);
+       pfx->SetBinEntries(391,33319);
+       pfx->SetBinEntries(392,32390);
+       pfx->SetBinEntries(393,30958);
+       pfx->SetBinEntries(394,29368);
+       pfx->SetBinEntries(395,28243);
+       pfx->SetBinEntries(396,26588);
+       pfx->SetBinEntries(397,24615);
+       pfx->SetBinEntries(398,23367);
+       pfx->SetBinEntries(399,21662);
+       pfx->SetBinEntries(400,20212);
+       pfx->SetBinEntries(401,137922);
+       pfx->SetBinContent(10,1.498);
+       pfx->SetBinContent(11,12.861);
+       pfx->SetBinContent(12,22.0525);
+       pfx->SetBinContent(13,49.161);
+       pfx->SetBinContent(14,91.483);
+       pfx->SetBinContent(15,147.269);
+       pfx->SetBinContent(16,229.526);
+       pfx->SetBinContent(17,332.995);
+       pfx->SetBinContent(18,392.96);
+       pfx->SetBinContent(19,486.707);
+       pfx->SetBinContent(20,575.162);
+       pfx->SetBinContent(21,656.894);
+       pfx->SetBinContent(22,700.745);
+       pfx->SetBinContent(23,779.266);
+       pfx->SetBinContent(24,844.906);
+       pfx->SetBinContent(25,867.478);
+       pfx->SetBinContent(26,899.238);
+       pfx->SetBinContent(27,933.53);
+       pfx->SetBinContent(28,951.69);
+       pfx->SetBinContent(29,958.502);
+       pfx->SetBinContent(30,944.578);
+       pfx->SetBinContent(31,955.648);
+       pfx->SetBinContent(32,936.921);
+       pfx->SetBinContent(33,930.725);
+       pfx->SetBinContent(34,906.534);
+       pfx->SetBinContent(35,896.805);
+       pfx->SetBinContent(36,877.532);
+       pfx->SetBinContent(37,863.957);
+       pfx->SetBinContent(38,835.466);
+       pfx->SetBinContent(39,816.927);
+       pfx->SetBinContent(40,796.281);
+       pfx->SetBinContent(41,767.347);
+       pfx->SetBinContent(42,747.322);
+       pfx->SetBinContent(43,728.795);
+       pfx->SetBinContent(44,696.885);
+       pfx->SetBinContent(45,669.992);
+       pfx->SetBinContent(46,644.381);
+       pfx->SetBinContent(47,610.65);
+       pfx->SetBinContent(48,580.651);
+       pfx->SetBinContent(49,553.92);
+       pfx->SetBinContent(50,525.714);
+       pfx->SetBinContent(51,496.879);
+       pfx->SetBinContent(52,471.746);
+       pfx->SetBinContent(53,444.109);
+       pfx->SetBinContent(54,408.403);
+       pfx->SetBinContent(55,381.361);
+       pfx->SetBinContent(56,359.956);
+       pfx->SetBinContent(57,338.895);
+       pfx->SetBinContent(58,325.912);
+       pfx->SetBinContent(59,309.628);
+       pfx->SetBinContent(60,288.36);
+       pfx->SetBinContent(61,280.942);
+       pfx->SetBinContent(62,258.548);
+       pfx->SetBinContent(63,253.744);
+       pfx->SetBinContent(64,246.012);
+       pfx->SetBinContent(65,236.367);
+       pfx->SetBinContent(66,220.55);
+       pfx->SetBinContent(67,218.771);
+       pfx->SetBinContent(68,217.38);
+       pfx->SetBinContent(69,216.244);
+       pfx->SetBinContent(70,208.884);
+       pfx->SetBinContent(71,200.755);
+       pfx->SetBinContent(72,201.114);
+       pfx->SetBinContent(73,204.406);
+       pfx->SetBinContent(74,206.551);
+       pfx->SetBinContent(75,200.242);
+       pfx->SetBinContent(76,206.853);
+       pfx->SetBinContent(77,204.993);
+       pfx->SetBinContent(78,205.45);
+       pfx->SetBinContent(79,217.247);
+       pfx->SetBinContent(80,225.895);
+       pfx->SetBinContent(81,213.368);
+       pfx->SetBinContent(82,216.671);
+       pfx->SetBinContent(83,217.432);
+       pfx->SetBinContent(84,221.402);
+       pfx->SetBinContent(85,220.779);
+       pfx->SetBinContent(86,227.347);
+       pfx->SetBinContent(87,207.753);
+       pfx->SetBinContent(88,222.917);
+       pfx->SetBinContent(89,213.499);
+       pfx->SetBinContent(90,221.03);
+       pfx->SetBinContent(91,215.225);
+       pfx->SetBinContent(92,212.267);
+       pfx->SetBinContent(93,208.294);
+       pfx->SetBinContent(94,205.216);
+       pfx->SetBinContent(95,203.488);
+       pfx->SetBinContent(96,211.815);
+       pfx->SetBinContent(97,211.315);
+       pfx->SetBinContent(98,200.546);
+       pfx->SetBinContent(99,209.305);
+       pfx->SetBinContent(100,204.662);
+       pfx->SetBinContent(101,214.492);
+       pfx->SetBinContent(102,203.016);
+       pfx->SetBinContent(103,214.964);
+       pfx->SetBinContent(104,212.629);
+       pfx->SetBinContent(105,198.984);
+       pfx->SetBinContent(106,207.573);
+       pfx->SetBinContent(107,201.142);
+       pfx->SetBinContent(108,192.879);
+       pfx->SetBinContent(109,186.784);
+       pfx->SetBinContent(110,211.004);
+       pfx->SetBinContent(111,201.427);
+       pfx->SetBinContent(112,203.052);
+       pfx->SetBinContent(113,203.737);
+       pfx->SetBinContent(114,197.814);
+       pfx->SetBinContent(115,190.139);
+       pfx->SetBinContent(116,180.696);
+       pfx->SetBinContent(117,181.102);
+       pfx->SetBinContent(118,180.142);
+       pfx->SetBinContent(119,178.683);
+       pfx->SetBinContent(120,213.262);
+       pfx->SetBinContent(121,188.73);
+       pfx->SetBinContent(122,183.369);
+       pfx->SetBinContent(123,196.616);
+       pfx->SetBinContent(124,184.075);
+       pfx->SetBinContent(125,182.267);
+       pfx->SetBinContent(126,180.315);
+       pfx->SetBinContent(127,170.849);
+       pfx->SetBinContent(128,192.712);
+       pfx->SetBinContent(129,162.001);
+       pfx->SetBinContent(130,171.485);
+       pfx->SetBinContent(131,176.516);
+       pfx->SetBinContent(132,158.295);
+       pfx->SetBinContent(133,168.941);
+       pfx->SetBinContent(134,184.208);
+       pfx->SetBinContent(135,178.957);
+       pfx->SetBinContent(136,182.447);
+       pfx->SetBinContent(137,166.983);
+       pfx->SetBinContent(138,158.138);
+       pfx->SetBinContent(139,185.531);
+       pfx->SetBinContent(140,171.153);
+       pfx->SetBinContent(141,156.523);
+       pfx->SetBinContent(142,144.669);
+       pfx->SetBinContent(143,164.95);
+       pfx->SetBinContent(144,156.538);
+       pfx->SetBinContent(145,151.055);
+       pfx->SetBinContent(146,163.336);
+       pfx->SetBinContent(147,157.911);
+       pfx->SetBinContent(148,163.984);
+       pfx->SetBinContent(149,168.494);
+       pfx->SetBinContent(150,149.863);
+       pfx->SetBinContent(151,173.995);
+       pfx->SetBinContent(152,180.54);
+       pfx->SetBinContent(153,185.124);
+       pfx->SetBinContent(154,175.516);
+       pfx->SetBinContent(155,147.426);
+       pfx->SetBinContent(156,170.625);
+       pfx->SetBinContent(157,155.471);
+       pfx->SetBinContent(158,169.775);
+       pfx->SetBinContent(159,158.756);
+       pfx->SetBinContent(160,167.447);
+       pfx->SetBinContent(161,137.689);
+       pfx->SetBinContent(162,177.396);
+       pfx->SetBinContent(163,173.447);
+       pfx->SetBinContent(164,153.374);
+       pfx->SetBinContent(165,177.654);
+       pfx->SetBinContent(166,181.88);
+       pfx->SetBinContent(167,187.271);
+       pfx->SetBinContent(168,194.266);
+       pfx->SetBinContent(169,194.034);
+       pfx->SetBinContent(170,192.259);
+       pfx->SetBinContent(171,199.7);
+       pfx->SetBinContent(172,182.962);
+       pfx->SetBinContent(173,184.878);
+       pfx->SetBinContent(174,185.94);
+       pfx->SetBinContent(175,213.635);
+       pfx->SetBinContent(176,209.018);
+       pfx->SetBinContent(177,202.12);
+       pfx->SetBinContent(178,187.68);
+       pfx->SetBinContent(179,201.394);
+       pfx->SetBinContent(180,195.191);
+       pfx->SetBinContent(181,189.723);
+       pfx->SetBinContent(182,204.585);
+       pfx->SetBinContent(183,198.534);
+       pfx->SetBinContent(184,216.789);
+       pfx->SetBinContent(185,224.281);
+       pfx->SetBinContent(186,176.761);
+       pfx->SetBinContent(187,208.037);
+       pfx->SetBinContent(188,195.28);
+       pfx->SetBinContent(189,203.074);
+       pfx->SetBinContent(190,201.993);
+       pfx->SetBinContent(191,202.994);
+       pfx->SetBinContent(192,223.693);
+       pfx->SetBinContent(193,220.655);
+       pfx->SetBinContent(194,203.101);
+       pfx->SetBinContent(195,249.928);
+       pfx->SetBinContent(196,254.327);
+       pfx->SetBinContent(197,246.391);
+       pfx->SetBinContent(198,216.177);
+       pfx->SetBinContent(199,204.509);
+       pfx->SetBinContent(200,236.13);
+       pfx->SetBinContent(201,190.965);
+       pfx->SetBinContent(202,260.77);
+       pfx->SetBinContent(203,285.472);
+       pfx->SetBinContent(204,223.67);
+       pfx->SetBinContent(205,297.955);
+       pfx->SetBinContent(206,256.319);
+       pfx->SetBinContent(207,234.07);
+       pfx->SetBinContent(208,218.88);
+       pfx->SetBinContent(209,303.794);
+       pfx->SetBinContent(210,265.703);
+       pfx->SetBinContent(211,265.09);
+       pfx->SetBinContent(212,268.681);
+       pfx->SetBinContent(213,268.929);
+       pfx->SetBinContent(214,249.607);
+       pfx->SetBinContent(215,325.919);
+       pfx->SetBinContent(216,297.889);
+       pfx->SetBinContent(217,264.758);
+       pfx->SetBinContent(218,277.494);
+       pfx->SetBinContent(219,257.925);
+       pfx->SetBinContent(220,261.516);
+       pfx->SetBinContent(221,258.697);
+       pfx->SetBinContent(222,310.868);
+       pfx->SetBinContent(223,293.226);
+       pfx->SetBinContent(224,274.617);
+       pfx->SetBinContent(225,298.093);
+       pfx->SetBinContent(226,323.255);
+       pfx->SetBinContent(227,336.509);
+       pfx->SetBinContent(228,325.168);
+       pfx->SetBinContent(229,315.734);
+       pfx->SetBinContent(230,283.758);
+       pfx->SetBinContent(231,342.337);
+       pfx->SetBinContent(232,306.045);
+       pfx->SetBinContent(233,359.416);
+       pfx->SetBinContent(234,325.759);
+       pfx->SetBinContent(235,323.507);
+       pfx->SetBinContent(236,346.311);
+       pfx->SetBinContent(237,358.701);
+       pfx->SetBinContent(238,358.237);
+       pfx->SetBinContent(239,400.736);
+       pfx->SetBinContent(240,324.672);
+       pfx->SetBinContent(241,407.269);
+       pfx->SetBinContent(242,347.431);
+       pfx->SetBinContent(243,345.247);
+       pfx->SetBinContent(244,338.534);
+       pfx->SetBinContent(245,351.188);
+       pfx->SetBinContent(246,373.245);
+       pfx->SetBinContent(247,358.937);
+       pfx->SetBinContent(248,380.312);
+       pfx->SetBinContent(249,395.081);
+       pfx->SetBinContent(250,365.726);
+       pfx->SetBinContent(251,339.133);
+       pfx->SetBinContent(252,366.415);
+       pfx->SetBinContent(253,399.1);
+       pfx->SetBinContent(254,376.422);
+       pfx->SetBinContent(255,379.349);
+       pfx->SetBinContent(256,371.868);
+       pfx->SetBinContent(257,382.217);
+       pfx->SetBinContent(258,357.484);
+       pfx->SetBinContent(259,399.403);
+       pfx->SetBinContent(260,414.294);
+       pfx->SetBinContent(261,402.987);
+       pfx->SetBinContent(262,411.965);
+       pfx->SetBinContent(263,429.635);
+       pfx->SetBinContent(264,362.686);
+       pfx->SetBinContent(265,386.67);
+       pfx->SetBinContent(266,409.149);
+       pfx->SetBinContent(267,364.283);
+       pfx->SetBinContent(268,353.16);
+       pfx->SetBinContent(269,379.556);
+       pfx->SetBinContent(270,381.357);
+       pfx->SetBinContent(271,381.931);
+       pfx->SetBinContent(272,456.608);
+       pfx->SetBinContent(273,413.225);
+       pfx->SetBinContent(274,395.008);
+       pfx->SetBinContent(275,375.988);
+       pfx->SetBinContent(276,394.091);
+       pfx->SetBinContent(277,351.908);
+       pfx->SetBinContent(278,408.176);
+       pfx->SetBinContent(279,428.575);
+       pfx->SetBinContent(280,388.168);
+       pfx->SetBinContent(281,421.85);
+       pfx->SetBinContent(282,404.47);
+       pfx->SetBinContent(283,418.938);
+       pfx->SetBinContent(284,430.91);
+       pfx->SetBinContent(285,393.423);
+       pfx->SetBinContent(286,454.77);
+       pfx->SetBinContent(287,388.166);
+       pfx->SetBinContent(288,496.85);
+       pfx->SetBinContent(289,381.939);
+       pfx->SetBinContent(290,444.934);
+       pfx->SetBinContent(291,430.5);
+       pfx->SetBinContent(292,453.678);
+       pfx->SetBinContent(293,468.948);
+       pfx->SetBinContent(294,370.429);
+       pfx->SetBinContent(295,431.494);
+       pfx->SetBinContent(296,426.921);
+       pfx->SetBinContent(297,436.992);
+       pfx->SetBinContent(298,445.874);
+       pfx->SetBinContent(299,465.897);
+       pfx->SetBinContent(300,462.106);
+       pfx->SetBinContent(301,442.342);
+       pfx->SetBinContent(302,431.233);
+       pfx->SetBinContent(303,465.8);
+       pfx->SetBinContent(304,443.678);
+       pfx->SetBinContent(305,440.162);
+       pfx->SetBinContent(306,344.091);
+       pfx->SetBinContent(307,428.903);
+       pfx->SetBinContent(308,423.999);
+       pfx->SetBinContent(309,453.669);
+       pfx->SetBinContent(310,442.489);
+       pfx->SetBinContent(311,382.607);
+       pfx->SetBinContent(312,478.946);
+       pfx->SetBinContent(313,454.324);
+       pfx->SetBinContent(314,485.697);
+       pfx->SetBinContent(315,421.75);
+       pfx->SetBinContent(316,480.34);
+       pfx->SetBinContent(317,422.854);
+       pfx->SetBinContent(318,465.354);
+       pfx->SetBinContent(319,414.979);
+       pfx->SetBinContent(320,401.623);
+       pfx->SetBinContent(321,470.957);
+       pfx->SetBinContent(322,421.762);
+       pfx->SetBinContent(323,385.705);
+       pfx->SetBinContent(324,418.958);
+       pfx->SetBinContent(325,453.342);
+       pfx->SetBinContent(326,384.963);
+       pfx->SetBinContent(327,417.203);
+       pfx->SetBinContent(328,419.263);
+       pfx->SetBinContent(329,428.683);
+       pfx->SetBinContent(330,435.04);
+       pfx->SetBinContent(331,500.481);
+       pfx->SetBinContent(332,446.147);
+       pfx->SetBinContent(333,428.009);
+       pfx->SetBinContent(334,438.342);
+       pfx->SetBinContent(335,492.926);
+       pfx->SetBinContent(336,403.515);
+       pfx->SetBinContent(337,516.049);
+       pfx->SetBinContent(338,465.155);
+       pfx->SetBinContent(339,447.05);
+       pfx->SetBinContent(340,419.668);
+       pfx->SetBinContent(341,437.139);
+       pfx->SetBinContent(342,433.226);
+       pfx->SetBinContent(343,456.919);
+       pfx->SetBinContent(344,410.913);
+       pfx->SetBinContent(345,395.724);
+       pfx->SetBinContent(346,428.227);
+       pfx->SetBinContent(347,363.568);
+       pfx->SetBinContent(348,411.832);
+       pfx->SetBinContent(349,381.018);
+       pfx->SetBinContent(350,396.683);
+       pfx->SetBinContent(351,456.052);
+       pfx->SetBinContent(352,401.231);
+       pfx->SetBinContent(353,346.811);
+       pfx->SetBinContent(354,420.841);
+       pfx->SetBinContent(355,425.153);
+       pfx->SetBinContent(356,426.565);
+       pfx->SetBinContent(357,341.458);
+       pfx->SetBinContent(358,430.449);
+       pfx->SetBinContent(359,452.143);
+       pfx->SetBinContent(360,366.352);
+       pfx->SetBinContent(361,461.106);
+       pfx->SetBinContent(362,476.188);
+       pfx->SetBinContent(363,349.512);
+       pfx->SetBinContent(364,394.65);
+       pfx->SetBinContent(365,482.426);
+       pfx->SetBinContent(366,347.851);
+       pfx->SetBinContent(367,334.394);
+       pfx->SetBinContent(368,400.737);
+       pfx->SetBinContent(369,411.661);
+       pfx->SetBinContent(370,386.772);
+       pfx->SetBinContent(371,417.292);
+       pfx->SetBinContent(372,486.134);
+       pfx->SetBinContent(373,413.773);
+       pfx->SetBinContent(374,422.605);
+       pfx->SetBinContent(375,365.114);
+       pfx->SetBinContent(376,358.082);
+       pfx->SetBinContent(377,374.774);
+       pfx->SetBinContent(378,345.616);
+       pfx->SetBinContent(379,325.011);
+       pfx->SetBinContent(380,409.335);
+       pfx->SetBinContent(381,284.952);
+       pfx->SetBinContent(382,172.566);
+       pfx->SetBinContent(383,77.5485);
+       pfx->SetBinContent(384,-139.827);
+       pfx->SetBinContent(385,-448.049);
+       pfx->SetBinContent(386,-619.486);
+       pfx->SetBinContent(387,-871.192);
+       pfx->SetBinContent(388,-1104.45);
+       pfx->SetBinContent(389,-1253.88);
+       pfx->SetBinContent(390,-1597.2);
+       pfx->SetBinContent(391,-1806.43);
+       pfx->SetBinContent(392,-2098.38);
+       pfx->SetBinContent(393,-2376.94);
+       pfx->SetBinContent(394,-2581.21);
+       pfx->SetBinContent(395,-2853.83);
+       pfx->SetBinContent(396,-2999.16);
+       pfx->SetBinContent(397,-3085.06);
+       pfx->SetBinContent(398,-3301.16);
+       pfx->SetBinContent(399,-3345.83);
+       pfx->SetBinContent(400,-3424.41);
+       pfx->SetBinContent(401,-36902.7);
+       pfx->SetBinError(10,0.616439);
+       pfx->SetBinError(11,1.72992);
+       pfx->SetBinError(12,2.23368);
+       pfx->SetBinError(13,3.10065);
+       pfx->SetBinError(14,4.07137);
+       pfx->SetBinError(15,5.01745);
+       pfx->SetBinError(16,6.05014);
+       pfx->SetBinError(17,7.05381);
+       pfx->SetBinError(18,7.34414);
+       pfx->SetBinError(19,7.87771);
+       pfx->SetBinError(20,8.24193);
+       pfx->SetBinError(21,8.46983);
+       pfx->SetBinError(22,8.38334);
+       pfx->SetBinError(23,8.51149);
+       pfx->SetBinError(24,8.51612);
+       pfx->SetBinError(25,8.31083);
+       pfx->SetBinError(26,8.17372);
+       pfx->SetBinError(27,8.03593);
+       pfx->SetBinError(28,7.82164);
+       pfx->SetBinError(29,7.59532);
+       pfx->SetBinError(30,7.28178);
+       pfx->SetBinError(31,7.10337);
+       pfx->SetBinError(32,6.82259);
+       pfx->SetBinError(33,6.59541);
+       pfx->SetBinError(34,6.33596);
+       pfx->SetBinError(35,6.12152);
+       pfx->SetBinError(36,5.88601);
+       pfx->SetBinError(37,5.7512);
+       pfx->SetBinError(38,5.51375);
+       pfx->SetBinError(39,5.33177);
+       pfx->SetBinError(40,5.16616);
+       pfx->SetBinError(41,5.00765);
+       pfx->SetBinError(42,4.83603);
+       pfx->SetBinError(43,4.73915);
+       pfx->SetBinError(44,4.59262);
+       pfx->SetBinError(45,4.45714);
+       pfx->SetBinError(46,4.32711);
+       pfx->SetBinError(47,4.19968);
+       pfx->SetBinError(48,4.09919);
+       pfx->SetBinError(49,4.00452);
+       pfx->SetBinError(50,3.92851);
+       pfx->SetBinError(51,3.80948);
+       pfx->SetBinError(52,3.75982);
+       pfx->SetBinError(53,3.71431);
+       pfx->SetBinError(54,3.68673);
+       pfx->SetBinError(55,3.68321);
+       pfx->SetBinError(56,3.57758);
+       pfx->SetBinError(57,3.56169);
+       pfx->SetBinError(58,3.64247);
+       pfx->SetBinError(59,3.59588);
+       pfx->SetBinError(60,3.6669);
+       pfx->SetBinError(61,3.67516);
+       pfx->SetBinError(62,3.72128);
+       pfx->SetBinError(63,3.77823);
+       pfx->SetBinError(64,3.85006);
+       pfx->SetBinError(65,3.90736);
+       pfx->SetBinError(66,3.95896);
+       pfx->SetBinError(67,3.90277);
+       pfx->SetBinError(68,4.10741);
+       pfx->SetBinError(69,4.15417);
+       pfx->SetBinError(70,4.15465);
+       pfx->SetBinError(71,4.30493);
+       pfx->SetBinError(72,4.45053);
+       pfx->SetBinError(73,4.4219);
+       pfx->SetBinError(74,4.68236);
+       pfx->SetBinError(75,4.56188);
+       pfx->SetBinError(76,4.66571);
+       pfx->SetBinError(77,4.87789);
+       pfx->SetBinError(78,4.95965);
+       pfx->SetBinError(79,5.00585);
+       pfx->SetBinError(80,5.10353);
+       pfx->SetBinError(81,5.10979);
+       pfx->SetBinError(82,5.15981);
+       pfx->SetBinError(83,5.35682);
+       pfx->SetBinError(84,5.41994);
+       pfx->SetBinError(85,5.47736);
+       pfx->SetBinError(86,5.50533);
+       pfx->SetBinError(87,5.65317);
+       pfx->SetBinError(88,5.81111);
+       pfx->SetBinError(89,5.87851);
+       pfx->SetBinError(90,6.01691);
+       pfx->SetBinError(91,6.01382);
+       pfx->SetBinError(92,6.24213);
+       pfx->SetBinError(93,6.30647);
+       pfx->SetBinError(94,6.45315);
+       pfx->SetBinError(95,6.41854);
+       pfx->SetBinError(96,6.53922);
+       pfx->SetBinError(97,6.7413);
+       pfx->SetBinError(98,6.82012);
+       pfx->SetBinError(99,6.89425);
+       pfx->SetBinError(100,6.94269);
+       pfx->SetBinError(101,7.09116);
+       pfx->SetBinError(102,7.21896);
+       pfx->SetBinError(103,7.35199);
+       pfx->SetBinError(104,7.40575);
+       pfx->SetBinError(105,7.60755);
+       pfx->SetBinError(106,7.65683);
+       pfx->SetBinError(107,7.73952);
+       pfx->SetBinError(108,8.0551);
+       pfx->SetBinError(109,8.10903);
+       pfx->SetBinError(110,8.19172);
+       pfx->SetBinError(111,8.26843);
+       pfx->SetBinError(112,8.45662);
+       pfx->SetBinError(113,8.61668);
+       pfx->SetBinError(114,8.69788);
+       pfx->SetBinError(115,8.82218);
+       pfx->SetBinError(116,8.90394);
+       pfx->SetBinError(117,9.00959);
+       pfx->SetBinError(118,9.21278);
+       pfx->SetBinError(119,9.30065);
+       pfx->SetBinError(120,9.46832);
+       pfx->SetBinError(121,9.54859);
+       pfx->SetBinError(122,9.59536);
+       pfx->SetBinError(123,9.79998);
+       pfx->SetBinError(124,9.86177);
+       pfx->SetBinError(125,10.0425);
+       pfx->SetBinError(126,10.1952);
+       pfx->SetBinError(127,10.3238);
+       pfx->SetBinError(128,10.335);
+       pfx->SetBinError(129,10.4642);
+       pfx->SetBinError(130,10.703);
+       pfx->SetBinError(131,10.7621);
+       pfx->SetBinError(132,10.8684);
+       pfx->SetBinError(133,11.0667);
+       pfx->SetBinError(134,11.2788);
+       pfx->SetBinError(135,11.3127);
+       pfx->SetBinError(136,11.4512);
+       pfx->SetBinError(137,11.5295);
+       pfx->SetBinError(138,11.8209);
+       pfx->SetBinError(139,11.8238);
+       pfx->SetBinError(140,11.9956);
+       pfx->SetBinError(141,12.0589);
+       pfx->SetBinError(142,12.2891);
+       pfx->SetBinError(143,12.5924);
+       pfx->SetBinError(144,12.4897);
+       pfx->SetBinError(145,12.8088);
+       pfx->SetBinError(146,12.8171);
+       pfx->SetBinError(147,13.072);
+       pfx->SetBinError(148,13.1604);
+       pfx->SetBinError(149,13.3261);
+       pfx->SetBinError(150,13.3748);
+       pfx->SetBinError(151,13.538);
+       pfx->SetBinError(152,13.6695);
+       pfx->SetBinError(153,13.9351);
+       pfx->SetBinError(154,13.91);
+       pfx->SetBinError(155,14.0188);
+       pfx->SetBinError(156,14.2745);
+       pfx->SetBinError(157,14.2934);
+       pfx->SetBinError(158,14.5653);
+       pfx->SetBinError(159,14.8644);
+       pfx->SetBinError(160,14.7273);
+       pfx->SetBinError(161,14.9128);
+       pfx->SetBinError(162,14.9543);
+       pfx->SetBinError(163,15.4038);
+       pfx->SetBinError(164,15.4801);
+       pfx->SetBinError(165,15.6339);
+       pfx->SetBinError(166,15.5855);
+       pfx->SetBinError(167,15.7744);
+       pfx->SetBinError(168,15.9024);
+       pfx->SetBinError(169,16.0523);
+       pfx->SetBinError(170,16.3169);
+       pfx->SetBinError(171,16.4003);
+       pfx->SetBinError(172,16.3849);
+       pfx->SetBinError(173,16.5948);
+       pfx->SetBinError(174,16.825);
+       pfx->SetBinError(175,16.9557);
+       pfx->SetBinError(176,17.0132);
+       pfx->SetBinError(177,17.3134);
+       pfx->SetBinError(178,17.2738);
+       pfx->SetBinError(179,17.5925);
+       pfx->SetBinError(180,17.8418);
+       pfx->SetBinError(181,17.9566);
+       pfx->SetBinError(182,17.9595);
+       pfx->SetBinError(183,18.1533);
+       pfx->SetBinError(184,18.2391);
+       pfx->SetBinError(185,18.3668);
+       pfx->SetBinError(186,18.6611);
+       pfx->SetBinError(187,18.8447);
+       pfx->SetBinError(188,18.762);
+       pfx->SetBinError(189,18.9102);
+       pfx->SetBinError(190,19.1766);
+       pfx->SetBinError(191,19.2297);
+       pfx->SetBinError(192,19.5992);
+       pfx->SetBinError(193,19.6522);
+       pfx->SetBinError(194,19.7484);
+       pfx->SetBinError(195,20.007);
+       pfx->SetBinError(196,20.2254);
+       pfx->SetBinError(197,20.1631);
+       pfx->SetBinError(198,20.3683);
+       pfx->SetBinError(199,20.5612);
+       pfx->SetBinError(200,20.6402);
+       pfx->SetBinError(201,20.9397);
+       pfx->SetBinError(202,20.8914);
+       pfx->SetBinError(203,21.3211);
+       pfx->SetBinError(204,21.1514);
+       pfx->SetBinError(205,21.6026);
+       pfx->SetBinError(206,21.5699);
+       pfx->SetBinError(207,21.517);
+       pfx->SetBinError(208,21.857);
+       pfx->SetBinError(209,22.2799);
+       pfx->SetBinError(210,22.0719);
+       pfx->SetBinError(211,22.3754);
+       pfx->SetBinError(212,22.5406);
+       pfx->SetBinError(213,22.6714);
+       pfx->SetBinError(214,22.8309);
+       pfx->SetBinError(215,22.8316);
+       pfx->SetBinError(216,23.1927);
+       pfx->SetBinError(217,23.1264);
+       pfx->SetBinError(218,23.2765);
+       pfx->SetBinError(219,23.5548);
+       pfx->SetBinError(220,23.6715);
+       pfx->SetBinError(221,23.6659);
+       pfx->SetBinError(222,24.2381);
+       pfx->SetBinError(223,24.0956);
+       pfx->SetBinError(224,24.3838);
+       pfx->SetBinError(225,24.2212);
+       pfx->SetBinError(226,24.516);
+       pfx->SetBinError(227,24.7666);
+       pfx->SetBinError(228,24.7664);
+       pfx->SetBinError(229,24.9714);
+       pfx->SetBinError(230,25.2049);
+       pfx->SetBinError(231,25.3975);
+       pfx->SetBinError(232,25.6472);
+       pfx->SetBinError(233,25.6502);
+       pfx->SetBinError(234,25.9516);
+       pfx->SetBinError(235,25.8526);
+       pfx->SetBinError(236,25.9609);
+       pfx->SetBinError(237,26.3826);
+       pfx->SetBinError(238,26.3914);
+       pfx->SetBinError(239,26.6302);
+       pfx->SetBinError(240,26.7064);
+       pfx->SetBinError(241,27.1028);
+       pfx->SetBinError(242,26.9378);
+       pfx->SetBinError(243,27.0565);
+       pfx->SetBinError(244,27.1023);
+       pfx->SetBinError(245,27.588);
+       pfx->SetBinError(246,27.5413);
+       pfx->SetBinError(247,27.7377);
+       pfx->SetBinError(248,27.9872);
+       pfx->SetBinError(249,27.9955);
+       pfx->SetBinError(250,27.9848);
+       pfx->SetBinError(251,28.1627);
+       pfx->SetBinError(252,28.3864);
+       pfx->SetBinError(253,28.7888);
+       pfx->SetBinError(254,28.7395);
+       pfx->SetBinError(255,29.0226);
+       pfx->SetBinError(256,29.1531);
+       pfx->SetBinError(257,29.3273);
+       pfx->SetBinError(258,29.392);
+       pfx->SetBinError(259,29.593);
+       pfx->SetBinError(260,29.8088);
+       pfx->SetBinError(261,29.8038);
+       pfx->SetBinError(262,30.019);
+       pfx->SetBinError(263,30.1702);
+       pfx->SetBinError(264,30.3958);
+       pfx->SetBinError(265,30.3525);
+       pfx->SetBinError(266,30.5933);
+       pfx->SetBinError(267,30.4665);
+       pfx->SetBinError(268,30.7536);
+       pfx->SetBinError(269,30.758);
+       pfx->SetBinError(270,31.0729);
+       pfx->SetBinError(271,31.0032);
+       pfx->SetBinError(272,31.7266);
+       pfx->SetBinError(273,31.876);
+       pfx->SetBinError(274,31.7338);
+       pfx->SetBinError(275,31.7029);
+       pfx->SetBinError(276,32.0293);
+       pfx->SetBinError(277,32.0787);
+       pfx->SetBinError(278,32.1201);
+       pfx->SetBinError(279,32.4714);
+       pfx->SetBinError(280,32.6217);
+       pfx->SetBinError(281,32.6574);
+       pfx->SetBinError(282,32.7444);
+       pfx->SetBinError(283,32.9243);
+       pfx->SetBinError(284,32.8649);
+       pfx->SetBinError(285,33.1107);
+       pfx->SetBinError(286,32.8976);
+       pfx->SetBinError(287,33.5005);
+       pfx->SetBinError(288,33.5264);
+       pfx->SetBinError(289,33.4835);
+       pfx->SetBinError(290,33.8739);
+       pfx->SetBinError(291,34.0327);
+       pfx->SetBinError(292,34.2425);
+       pfx->SetBinError(293,34.249);
+       pfx->SetBinError(294,34.3511);
+       pfx->SetBinError(295,34.4373);
+       pfx->SetBinError(296,34.4429);
+       pfx->SetBinError(297,34.5659);
+       pfx->SetBinError(298,34.805);
+       pfx->SetBinError(299,35.0643);
+       pfx->SetBinError(300,35.0076);
+       pfx->SetBinError(301,35.1115);
+       pfx->SetBinError(302,35.3025);
+       pfx->SetBinError(303,35.3357);
+       pfx->SetBinError(304,35.5786);
+       pfx->SetBinError(305,35.7108);
+       pfx->SetBinError(306,35.6469);
+       pfx->SetBinError(307,35.8856);
+       pfx->SetBinError(308,36.2562);
+       pfx->SetBinError(309,36.114);
+       pfx->SetBinError(310,36.2775);
+       pfx->SetBinError(311,36.0096);
+       pfx->SetBinError(312,36.6171);
+       pfx->SetBinError(313,36.5626);
+       pfx->SetBinError(314,36.9821);
+       pfx->SetBinError(315,36.8483);
+       pfx->SetBinError(316,36.9698);
+       pfx->SetBinError(317,36.9896);
+       pfx->SetBinError(318,37.1308);
+       pfx->SetBinError(319,37.1591);
+       pfx->SetBinError(320,37.2727);
+       pfx->SetBinError(321,37.4798);
+       pfx->SetBinError(322,37.5427);
+       pfx->SetBinError(323,37.5199);
+       pfx->SetBinError(324,37.9852);
+       pfx->SetBinError(325,37.8331);
+       pfx->SetBinError(326,38.0506);
+       pfx->SetBinError(327,38.2952);
+       pfx->SetBinError(328,38.0231);
+       pfx->SetBinError(329,38.4839);
+       pfx->SetBinError(330,38.3669);
+       pfx->SetBinError(331,38.6735);
+       pfx->SetBinError(332,38.3288);
+       pfx->SetBinError(333,38.9733);
+       pfx->SetBinError(334,38.7921);
+       pfx->SetBinError(335,38.8858);
+       pfx->SetBinError(336,39.0047);
+       pfx->SetBinError(337,39.2976);
+       pfx->SetBinError(338,39.0848);
+       pfx->SetBinError(339,39.359);
+       pfx->SetBinError(340,39.1442);
+       pfx->SetBinError(341,39.3762);
+       pfx->SetBinError(342,39.4192);
+       pfx->SetBinError(343,39.6347);
+       pfx->SetBinError(344,39.777);
+       pfx->SetBinError(345,39.7949);
+       pfx->SetBinError(346,40.0307);
+       pfx->SetBinError(347,40.1193);
+       pfx->SetBinError(348,40.193);
+       pfx->SetBinError(349,40.0144);
+       pfx->SetBinError(350,40.4891);
+       pfx->SetBinError(351,40.3722);
+       pfx->SetBinError(352,40.6616);
+       pfx->SetBinError(353,40.3099);
+       pfx->SetBinError(354,40.6565);
+       pfx->SetBinError(355,40.8992);
+       pfx->SetBinError(356,40.9322);
+       pfx->SetBinError(357,41.3121);
+       pfx->SetBinError(358,40.9736);
+       pfx->SetBinError(359,40.9073);
+       pfx->SetBinError(360,41.2808);
+       pfx->SetBinError(361,41.3006);
+       pfx->SetBinError(362,41.312);
+       pfx->SetBinError(363,41.2755);
+       pfx->SetBinError(364,41.4478);
+       pfx->SetBinError(365,41.9378);
+       pfx->SetBinError(366,41.4831);
+       pfx->SetBinError(367,41.606);
+       pfx->SetBinError(368,41.6426);
+       pfx->SetBinError(369,41.9574);
+       pfx->SetBinError(370,42.1263);
+       pfx->SetBinError(371,42.096);
+       pfx->SetBinError(372,42.1522);
+       pfx->SetBinError(373,41.8282);
+       pfx->SetBinError(374,42.2305);
+       pfx->SetBinError(375,42.3714);
+       pfx->SetBinError(376,42.6419);
+       pfx->SetBinError(377,42.4077);
+       pfx->SetBinError(378,42.2049);
+       pfx->SetBinError(379,42.5347);
+       pfx->SetBinError(380,42.8112);
+       pfx->SetBinError(381,42.2039);
+       pfx->SetBinError(382,41.5359);
+       pfx->SetBinError(383,40.6853);
+       pfx->SetBinError(384,39.8128);
+       pfx->SetBinError(385,38.9079);
+       pfx->SetBinError(386,38.0542);
+       pfx->SetBinError(387,37.0357);
+       pfx->SetBinError(388,35.9301);
+       pfx->SetBinError(389,34.838);
+       pfx->SetBinError(390,34.1853);
+       pfx->SetBinError(391,33.1597);
+       pfx->SetBinError(392,32.4606);
+       pfx->SetBinError(393,31.7274);
+       pfx->SetBinError(394,31.0704);
+       pfx->SetBinError(395,30.889);
+       pfx->SetBinError(396,30.418);
+       pfx->SetBinError(397,29.8384);
+       pfx->SetBinError(398,30.18);
+       pfx->SetBinError(399,29.8639);
+       pfx->SetBinError(400,30.1144);
+       pfx->SetBinError(401,108.979);
+       pfx->SetEntries(1.52688e+07);
+       
+       return pfx;
+}