//*-- AUTHOR : Olga Pechenova
//*-- Modified : 2006 by Olga Pechenova
//*-- Modified : 2013 by Olga Pechenova
//*-- Modified : 2015 by Vladimir Pechenov
//*-- Modified : 2019 by Olga Pechenova
//*-- Modified : 2023 by Olga Pechenova
//*-- Modified : 2024 by Olga Pechenova
//
//_HADES_CLASS_DESCRIPTION
/////////////////////////////////////////////////////////////////////////
//
//  HEnergyLossCorrPar
//
//  Container class keep parameters for energy loss correction
//
//
//  Double_t getDeltaMom(Int_t pId, Double_t mom, Double_t theta);
//     Function getDeltaMom returns the correction value (delP) for momentum.
//     Parameters: - pId is particle ID. If parameter pId is not correct
//  function getDeltaMom returns 0.
//  - mom is Runge-Kutta reconstructed momentum multiplied by particle charge (MeV).
//  - theta is theta angle in the interval [15;87] (degrees).
//  If parameter theta is out of this interval
//  function getDeltaMom returns 0.
//
//  Double_t getCorrMom(Int_t pId, Double_t mom, Double_t theta) {
//          return mom+getDeltaMom(pId,mom,theta);}
//     Function getCorrMom returns the corrected momentum.
//     Parameters have the same meanings as for function getDeltaMom.
//
//
//  Alternative method of getting parameters (without usage of runTimeDb):
//
//  Bool_t  setDefaultPar(const Char_t *rn);
//     Function setDefaultPar sets default values to parameters
//     for corresponding beam time.
//
//  Example:
//        HEnergyLossCorrPar enLossCorr;
//        enLossCorr.setDefaultPar("jan04"); // or "nov02", or "aug04", or "apr07", or "apr12", or "mar19", or "feb21", or "feb22", 
//        or "feb24_auau", or "feb24_cc", or
//        "jul14_w" - Wolfram target, "jul14_c3" - Carbon 3 segments, "jul14_pe" & "aug14_pe" - PE target, "aug14_c7" - Carbon 7 segments
//        Double_t delP;
//        ..............................
//        beginning of your loop;
//          ..............................................
//          delP = enLossCorr.getDeltaMom(14, mom, theta);
//          ..............................................
//        end of your loop;
/////////////////////////////////////////////////////////////////////////

#include "TMath.h"
#include "TString.h"
#include "henergylosscorrpar.h"
#include "hphysicsconstants.h"
#include "hparamlist.h"

#include <stdlib.h>
#include <cmath>

using namespace TMath;

HEnergyLossCorrPar* HEnergyLossCorrPar::gEnergyLossCorrPar = NULL;

ClassImp(HEnergyLossCorrPar)

HEnergyLossCorrPar::HEnergyLossCorrPar(const Char_t* name,const Char_t* title, const Char_t* context) {
  // constructor
  momInt1=720.;
  momInt2=2000.;
  nThetaReg=7;
  nParams=10;
  thetaStep=10.;
  thetaMiddle0=23.;
  for(Int_t i=0;i<nThetaReg*nParams;i++) parMomCorrH[i]=0.;
// for electrons
  nParamsElect=5;
  for(Int_t i=0;i<nParamsElect;i++) parMomCorrElect[i]=0.;
  typePar = 0;
  gEnergyLossCorrPar =  this;
}

void HEnergyLossCorrPar::clear(void) {
   for(Int_t i=0;i<nThetaReg*nParams;i++) parMomCorrH[i]=0.;
  for(Int_t i=0;i<nParamsElect;i++) parMomCorrElect[i]=0.;
  status=kFALSE;
  resetInputVersions();
  changed=kFALSE;
  typePar = 0;
}

void HEnergyLossCorrPar::putParams(HParamList* l) {
  // puts all parameters to the parameter list, which is used by the io
  if (!l) return;
  l->add("momInt1",        momInt1);
  l->add("momInt2",        momInt2);
  l->add("nThetaReg",      nThetaReg);
  l->add("nParams",    nParams);
  l->add("parMomCorrH", parMomCorrH,nThetaReg*nParams);
  l->add("thetaMiddle0",   thetaMiddle0);
  l->add("thetaStep",      thetaStep);

  l->add("nParamsElect",    nParamsElect);
  l->add("parMomCorrElect", parMomCorrElect,nParamsElect);
}


Bool_t HEnergyLossCorrPar::getParams(HParamList* l) {
  // gets all parameters from the parameter list, which is used by the io
  if (!l) return kFALSE;
  if (!(l->fill("momInt1",       &momInt1)))          return kFALSE;
  if (!(l->fill("momInt2",       &momInt2)))          return kFALSE;
  if (!(l->fill("nThetaReg",     &nThetaReg)))        return kFALSE;
  if (!(l->fill("nParams",   &nParams)))      return kFALSE;
  if (!(l->fill("parMomCorrH",parMomCorrH,nThetaReg*nParams)))   return kFALSE;
  if (!(l->fill("thetaMiddle0",  &thetaMiddle0)))     return kFALSE;
  if (!(l->fill("thetaStep",     &thetaStep)))        return kFALSE;

  if (!(l->fill("nParamsElect",       &nParamsElect)))     return kFALSE;
  if (!(l->fill("parMomCorrElect",parMomCorrElect,nParamsElect)))   return kFALSE;
  return kTRUE;
}

Bool_t HEnergyLossCorrPar::setDefaultPar(TString run) {
  run.ToLower();

  if(run == "feb24_cc") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 86 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {

    303.006, -1.38858, 0.00240227, -1.64314e-06, 2.77135e-10, 88.039, -0.00363459, -0.000686509, 5.01122, 1.79959, -0.00065403, 2.96492,
    303.655, -1.44746, 0.00268667, -2.12323e-06, 5.49817e-10, 76.4089, -0.00347865, -0.000571898, 4.6041, 5.3473, -0.00108086, 2.86314,
    316.404, -1.6131, 0.00334704, -3.18153e-06, 1.14445e-09, 82.2106, -0.00363159, -0.000666929, 4.76375, 8.43975, -0.00128126, 2.84036,
    315.462, -1.58515, 0.00320647, -2.92685e-06, 9.89635e-10, 83.8184, -0.00365294, -0.000673311, 4.78907, 10.6355, -0.00137581, 2.85119,

    309.696, -1.47869, 0.00273315, -2.12665e-06, 5.25875e-10, 86.3822, -0.00366243, -0.000667636, 4.8427, 21.2177, -0.00154982, 2.87652,
    308.152, -1.46431, 0.00269889, -2.09942e-06, 5.21346e-10, 79.6432, -0.00352577, -0.000615128, 4.75629, 5.5796, -0.00104449, 2.87559,
    304.619, -1.44515, 0.002676, -2.11674e-06, 5.5255e-10, 85.7427, -0.00363281, -0.000689663, 4.9336, 4.39874, -0.000955355, 2.88182,
    309.073, -1.46652, 0.00271042, -2.13137e-06, 5.4713e-10, 86.147, -0.00362293, -0.00067359, 4.95785, 2.65796, -0.000757914, 2.90706,

    314.988, -1.53719, 0.00299159, -2.59004e-06, 8.11211e-10, 80.5933, -0.00350656, -0.000602258, 4.79388, 3.26996, -0.000845112, 2.9219,
    311.047, -1.54227, 0.00305142, -2.68747e-06, 8.57109e-10, 76.0539, -0.00350692, -0.000596977, 4.60262, 3.29042, -0.00085851, 2.77005,
    292.574, -1.55591, 0.00337729, -3.37381e-06, 1.28268e-09, 59.279, -0.00340819, -0.000505887, 3.90339, 3.45235, -0.0009639, 2.37661,
    296.278, -1.67995, 0.0039121, -4.21879e-06, 1.74126e-09, 58.4898, -0.00353243, -0.00048132, 3.53614, 3.87861, -0.00102926, 2.1173,

    289.429, -1.56918, 0.0034444, -3.45158e-06, 1.30542e-09, 64.4046, -0.00362938, -0.000495719, 3.66276, 3.10905, -0.000896729, 2.14958,
    298.526, -1.65956, 0.00380507, -4.05786e-06, 1.6638e-09, 64.8579, -0.00357006, -0.000505468, 3.78172, 13.8281, -0.00146572, 2.27977,
    309.069, -1.71917, 0.00394398, -4.20342e-06, 1.71963e-09, 65.9728, -0.00356041, -0.000546532, 3.98918, 1.2496, -0.000616838, 2.35983,
    296.173, -1.55166, 0.00330002, -3.20695e-06, 1.17634e-09, 62.7341, -0.00341173, -0.000496887, 3.99509, 4.73125, -0.0011588, 2.50118,

    305.678, -1.60798, 0.00346382, -3.44468e-06, 1.30925e-09, 87.1306, -0.00381197, -0.000679172, 4.57323, 3.67413, -0.000931057, 2.58745,
    298.311, -1.47303, 0.00297464, -2.77045e-06, 9.86393e-10, 84.93, -0.00359236, -0.00062634, 4.77203, 3.3842, -0.000867172, 2.84268,
    294.272, -1.32477, 0.00226923, -1.5562e-06, 2.77275e-10, 89.4942, -0.0035942, -0.00068521, 5.10163, 3.59946, -0.00090692, 3.04027,
    289.531, -1.21973, 0.00179578, -7.63219e-07, -1.7645e-10, 87.4161, -0.00354109, -0.000675191, 5.14265, 8.05539, -0.00120079, 3.11028,

    309.601, -1.41484, 0.00247939, -1.77495e-06, 3.59752e-10, 86.322, -0.00352193, -0.000674994, 5.17066, 5.58086, -0.00104167, 3.12005,
    304.577, -1.3973, 0.00245747, -1.76066e-06, 3.53114e-10, 86.1075, -0.00354863, -0.000717089, 5.22778, 6.42839, -0.00108417, 3.09305

    };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons pol4  pol2  pol4
    nParamsElect=13;

    parMomCorrElect[0] = 5.42717;
    parMomCorrElect[1] = -0.19857;
    parMomCorrElect[2] = 0.00401446;
    parMomCorrElect[3] = 1.5228e-05;
    parMomCorrElect[4] = -7.27793e-07;

    parMomCorrElect[5] = 10.1553;
    parMomCorrElect[6] = -0.155653;
    parMomCorrElect[7] = 0.000110523;

    parMomCorrElect[8] = -84.2622;
    parMomCorrElect[9] = 5.44626;
    parMomCorrElect[10] = -0.128162;
    parMomCorrElect[11] = 0.00132825;
    parMomCorrElect[12] = -5.07338e-06;

/*
    parMomCorrElect[0] = 5.95731;//-1.59259;//-9.13308;
    parMomCorrElect[1] = -0.257065;//0.707442;//1.40618;
    parMomCorrElect[2] = 0.00636724;//-0.0386977;//-0.0630954;
    parMomCorrElect[3] = -2.60238e-05;//0.000887115;//0.00123199;
    parMomCorrElect[4] = -4.60898e-07;//-7.23786e-06;//-0.00000884375;

    parMomCorrElect[5] = 7.80256;//90.1698;  //   114.165;
    parMomCorrElect[6] = -0.0635966;//-3.23091; //-4.14175;
    parMomCorrElect[7] = -0.000789466;//0.0294168;  //0.0380465;

    parMomCorrElect[8] = -101.73;//169.332;//-4.75336;
    parMomCorrElect[9] = 6.54043;//-9.95574;//0.389213;
    parMomCorrElect[10] = -0.153771;//0.22064;//-0.00910909;
    parMomCorrElect[11] = 0.00159367;//-0.00216254;//0.0000934938;
    parMomCorrElect[12] = -6.1012e-06;//7.95288e-06;//-0.000000348053;
*/
    typePar = 242;

    } else if(run == "feb24_auau") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 86 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {

    304.746, -1.64379, 0.00358106, -3.5583e-06, 1.33442e-09, 84.8596, -0.00408256, -0.000645663, 4.26532, 3.22121, -0.000930394, 2.43311,
    294.618, -1.54638, 0.00324714, -3.07213e-06, 1.07925e-09, 75.3093, -0.00385847, -0.000605914, 4.14173, 3.1804, -0.000934435, 2.37077,
    277.803, -1.36167, 0.00254398, -1.96417e-06, 4.61441e-10, 65.7023, -0.00365005, -0.000538012, 3.93448, 3.211, -0.000938208, 2.32759,
    293.994, -1.55937, 0.00331988, -3.19717e-06, 1.14956e-09, 69.4632, -0.00372962, -0.000570091, 3.94316, 3.1462, -0.000936711, 2.26858,

    297.558, -1.58989, 0.00341881, -3.33633e-06, 1.2208e-09, 67.6669, -0.00367021, -0.00055776, 3.89072, 3.06998, -0.000932255, 2.25319,
    292.855, -1.53267, 0.00320188, -3.00246e-06, 1.03919e-09, 66.5187, -0.00363682, -0.000552656, 3.88893, 3.13381, -0.000942005, 2.246,
    294.888, -1.55656, 0.00329931, -3.16449e-06, 1.13331e-09, 63.8593, -0.00358101, -0.000540152, 3.84645, 3.26386, -0.000957065, 2.25318,
    297.911, -1.57948, 0.00337396, -3.2752e-06, 1.19364e-09, 65.4094, -0.00360458, -0.000550764, 3.89243, 3.05379, -0.000934796, 2.27452,

    295.559, -1.54476, 0.0032341, -3.05547e-06, 1.07365e-09, 63.8856, -0.00355456, -0.00052988, 3.85437, 3.09311, -0.000939998, 2.27621,
    292.587, -1.55619, 0.00332653, -3.2261e-06, 1.17273e-09, 62.0811, -0.00356603, -0.000494075, 3.60556, 3.05252, -0.000959643, 2.13345,
    272.638, -1.59786, 0.00380569, -4.16958e-06, 1.74028e-09, 50.1228, -0.00361515, -0.000400746, 2.78222, 1.32278, -0.000701439, 1.59509,
    258.063, -1.60022, 0.004004, -4.58535e-06, 1.99167e-09, 31.8948, -0.0033026, -0.000274909, 2.11017, 39.3018, -0.00210209, 1.3332,

    259.232, -1.60419, 0.00400525, -4.5754e-06, 1.98181e-09, 33.9198, -0.00336783, -0.000281881, 2.0977, 2.33268, -0.00104033, 1.3287,
    258.003, -1.58737, 0.00394729, -4.49859e-06, 1.94679e-09, 35.5023, -0.00341252, -0.000290521, 2.21085, 36.0706, -0.00206541, 1.3647,
    257.823, -1.57278, 0.00386859, -4.35322e-06, 1.85782e-09, 34.9853, -0.00338252, -0.000288766, 2.22882, 7.82166, -0.00148414, 1.37755,
    261.252, -1.60105, 0.00396399, -4.49585e-06, 1.93561e-09, 35.6864, -0.00339117, -0.00029294, 2.22869, 2.44016, -0.00104133, 1.38704,

    261.442, -1.58339, 0.00386853, -4.32588e-06, 1.83546e-09, 35.0082, -0.00334965, -0.000293633, 2.26297, 34.1412, -0.00200452, 1.42761,
    263.758, -1.5951, 0.00389927, -4.36822e-06, 1.85845e-09, 37.2397, -0.00339793, -0.000308061, 2.33133, 2.4162, -0.0010316, 1.45084,
    265.287, -1.5956, 0.00388197, -4.3316e-06, 1.83705e-09, 37.8252, -0.00341342, -0.000322678, 2.37831, 2.53111, -0.00104784, 1.47767,
    266.437, -1.58887, 0.00383158, -4.23563e-06, 1.77878e-09, 38.9325, -0.00341692, -0.000326972, 2.45264, 58.1024, -0.00220484, 1.53649,

    267.186, -1.57694, 0.00376094, -4.10953e-06, 1.70537e-09, 41.7358, -0.00347897, -0.000349382, 2.53929, 2.84435, -0.00107643, 1.54575,
    276.75, -1.66725, 0.00408415, -4.60109e-06, 1.9725e-09, 44.4892, -0.00354458, -0.000369726, 2.65931, 35.9983, -0.00203962, 1.6179

    };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons pol4  (-pol2-)  pol4
    nParamsElect=13;
    parMomCorrElect[0] = -9.13308;
    parMomCorrElect[1] = 1.40618;
    parMomCorrElect[2] = -0.0630954;
    parMomCorrElect[3] = 0.00123199;
    parMomCorrElect[4] = -0.00000884375;
    parMomCorrElect[5] = 90.1698;  //   114.165;
    parMomCorrElect[6] = -3.23091; //-4.14175;
    parMomCorrElect[7] = 0.0294168;  //0.0380465;
    parMomCorrElect[8] = -4.75336;
    parMomCorrElect[9] = 0.389213;
    parMomCorrElect[10] = -0.00910909;
    parMomCorrElect[11] = 0.0000934938;
    parMomCorrElect[12] = -0.000000348053;

    typePar = 243;

    } else if(run == "feb22" || run == "feb21") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 87 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {

        359.248, -1.77623, 0.00363851, -3.46574e-06, 1.26721e-09, 95.9595, -0.00345436, -0.000762843, 6.05142, 4.80381, -0.00092181, 3.72064,
        353.04, -1.76783, 0.00366574, -3.53817e-06, 1.31346e-09, 91.3986, -0.00346478, -0.000731475, 5.74866, 4.43751, -0.000916614, 3.52476,
        349.678, -1.76931, 0.00370159, -3.6001e-06, 1.34492e-09, 87.9473, -0.00346051, -0.000695284, 5.51358, 4.31816, -0.00092009, 3.3923,
        349.159, -1.79593, 0.00382935, -3.80768e-06, 1.45902e-09, 84.835, -0.00345143, -0.000670064, 5.33114, 4.1564, -0.000921065, 3.28685,

        346.5, -1.78629, 0.00381137, -3.78552e-06, 1.44607e-09, 83.8763, -0.00346543, -0.000670669, 5.2715, 3.95739, -0.00090897, 3.23313,
        345.078, -1.78338, 0.00381562, -3.80223e-06, 1.4581e-09, 84.054, -0.00347904, -0.000667793, 5.22552, 4.16233, -0.000932374, 3.20043,
        345.597, -1.79644, 0.00386622, -3.87582e-06, 1.49549e-09, 82.8378, -0.00346924, -0.000647305, 5.14346, 3.9847, -0.000918991, 3.16682,
        343.929, -1.78345, 0.00382558, -3.81945e-06, 1.46655e-09, 86.1025, -0.00353439, -0.00068766, 5.19498, 3.90243, -0.000902841, 3.12157,

        343.513, -1.77554, 0.00378327, -3.74002e-06, 1.41742e-09, 89.5073, -0.0035892, -0.000674746, 5.07233, 4.87956, -0.00103212, 3.06659,
        334.984, -1.66747, 0.00332733, -2.97696e-06, 9.76991e-10, 83.8615, -0.00351493, -0.000643751, 5.01227, 6.20293, -0.00111573, 3.08511,
        335.114, -1.66738, 0.00334393, -3.02554e-06, 1.0122e-09, 86.798, -0.00357361, -0.00071215, 5.23325, 4.21407, -0.000935843, 3.10342,
        329.613, -1.60338, 0.00311138, -2.6834e-06, 8.3512e-10, 78.1484, -0.00339493, -0.000622829, 5.05874, 4.01635, -0.000930831, 3.14914,

        333.457, -1.64629, 0.0032956, -3.0063e-06, 1.03026e-09, 88.6664, -0.00358357, -0.000734817, 5.34336, 4.51024, -0.000968568, 3.17565,
        336.403, -1.65996, 0.00327639, -2.86973e-06, 9.00503e-10, 92.7996, -0.00363309, -0.000753519, 5.42233, 4.96443, -0.00100254, 3.2001,
        329.414, -1.56759, 0.00291584, -2.30682e-06, 5.9218e-10, 81.163, -0.00339407, -0.000625195, 5.18887, 3.98394, -0.000907139, 3.2494,
        328.766, -1.55178, 0.00290452, -2.39337e-06, 6.97968e-10, 91.9452, -0.00358296, -0.000743988, 5.48863, 4.80689, -0.000985434, 3.28096,

        332.567, -1.57675, 0.00295414, -2.41258e-06, 6.81701e-10, 89.5157, -0.00352186, -0.000739053, 5.53444, 3.36954, -0.000825273, 3.31992,
        332.362, -1.56082, 0.00291618, -2.40708e-06, 7.06648e-10, 89.2491, -0.00349655, -0.000722587, 5.55589, 3.56105, -0.000853404, 3.37547,
        335.395, -1.56875, 0.00286882, -2.22955e-06, 5.6045e-10, 92.1577, -0.00349204, -0.000714938, 5.60443, 5.63803, -0.0010531, 3.45033,
        344.479, -1.67277, 0.00328136, -2.8974e-06, 9.39316e-10, 106.595, -0.00370404, -0.000820909, 5.90645, 3.36002, -0.000834956, 3.49288,

        344.479, -1.67277, 0.00328136, -2.8974e-06, 9.39316e-10, 106.595, -0.00370404, -0.000820909, 5.90645, 3.36002, -0.000834956, 3.49288,
        344.479, -1.67277, 0.00328136, -2.8974e-06, 9.39316e-10, 106.595, -0.00370404, -0.000820909, 5.90645, 3.36002, -0.000834956, 3.49288

      };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons exp+pol1  exp  pol4
    nParamsElect=12;
    parMomCorrElect[0] = 0.880255;
    parMomCorrElect[1] = 0.0318177;
    parMomCorrElect[2] = -0.103786;
    parMomCorrElect[3] = 4.17457;
    parMomCorrElect[4] = 1.6266;
    parMomCorrElect[5] = -0.0473965;
    parMomCorrElect[6] = 2.93548;
    parMomCorrElect[7] = 0.847754;
    parMomCorrElect[8] = 0.186738;
    parMomCorrElect[9] = -0.00538821;
    parMomCorrElect[10] = 0.0000644817;
    parMomCorrElect[11] = -0.000000264453;

    typePar = 22;

    } else if(run == "mar19") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 86 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {

   331.127, -1.79734, 0.00405562, -4.27346e-06, 1.735e-09, 73.6073, -0.00345576, -0.000582559, 4.63143, 3.88641, -0.0009463, 2.85763,
   323.576, -1.79889, 0.00412357, -4.39188e-06, 1.79612e-09, 67.2306, -0.00350626, -0.000533308, 4.11317, 3.17558, -0.000925526, 2.50397,
   319.744, -1.7855, 0.00410667, -4.38536e-06, 1.79703e-09, 65.5867, -0.00350781, -0.000513687, 3.98569, 3.24107, -0.000941404, 2.43316,
   317.771, -1.78131, 0.00410794, -4.39378e-06, 1.80148e-09, 63.6918, -0.00350117, -0.000508197, 3.90963, 3.1139, -0.000937999, 2.38062,

   317.133, -1.78045, 0.00411114, -4.40261e-06, 1.80752e-09, 62.586, -0.00348684, -0.000496123, 3.86303, 3.16102, -0.000945441, 2.36431,
   316.082, -1.77001, 0.00407816, -4.36077e-06, 1.78908e-09, 63.4748, -0.00350647, -0.000501702, 3.87875, 3.12845, -0.0009425, 2.3655,
   318.695, -1.79788, 0.00417514, -4.50049e-06, 1.86117e-09, 63.9636, -0.00351873, -0.000503097, 3.87194, 3.06006, -0.000936491, 2.35765,
   317.629, -1.78201, 0.00411396, -4.40713e-06, 1.81089e-09, 64.1952, -0.00351285, -0.000501238, 3.88635, 3.0865, -0.000937421, 2.37409,

   317.305, -1.78027, 0.00410967, -4.40112e-06, 1.8073e-09, 64.5274, -0.00352598, -0.00050951, 3.90579, 3.15457, -0.000944896, 2.37474,
   307.19, -1.73526, 0.00402194, -4.31664e-06, 1.77418e-09, 58.3901, -0.00346572, -0.000454542, 3.61535, 3.04143, -0.000955256, 2.23254,
   269.605, -1.56635, 0.00371216, -4.05734e-06, 1.69292e-09, 47.1203, -0.00348349, -0.000365663, 2.89917, 2.50292, -0.000975642, 1.78891,
   252.504, -1.53331, 0.00376932, -4.25092e-06, 1.82238e-09, 37.7532, -0.00346699, -0.000292253, 2.33397, 2.38971, -0.00104616, 1.45284,

   253.178, -1.54433, 0.0038151, -4.32519e-06, 1.86432e-09, 38.5394, -0.00350577, -0.000296633, 2.33331, 2.31072, -0.00103731, 1.44315,
   253.705, -1.54277, 0.00379995, -4.29587e-06, 1.84681e-09, 38.8674, -0.00350881, -0.000302755, 2.36416, 2.23024, -0.00102253, 1.45752,
   255.392, -1.55332, 0.00382658, -4.32641e-06, 1.86004e-09, 39.5061, -0.00352028, -0.000308757, 2.39004, 2.3228, -0.00102958, 1.46958,
   257.281, -1.56217, 0.00384362, -4.34149e-06, 1.86503e-09, 40.0745, -0.00351683, -0.000309021, 2.41602, 2.30394, -0.00102783, 1.4921,

   259.739, -1.57502, 0.00387274, -4.37386e-06, 1.87949e-09, 40.735, -0.00351943, -0.000314866, 2.45859, 2.4122, -0.00103648, 1.51779,
   262.203, -1.58403, 0.00388233, -4.37164e-06, 1.87323e-09, 41.2454, -0.00350699, -0.00031877, 2.50649, 2.58409, -0.00105604, 1.55309,
   264.451, -1.59371, 0.00389996, -4.38753e-06, 1.87931e-09, 42.3023, -0.00351954, -0.000330261, 2.56556, 2.61165, -0.00105409, 1.58334,
   267.277, -1.6067, 0.00392585, -4.41284e-06, 1.88931e-09, 43.3361, -0.00351927, -0.000336194, 2.62386, 2.77194, -0.00107086, 1.62278,

   269.713, -1.6115, 0.00391444, -4.37416e-06, 1.86178e-09, 43.6594, -0.00348987, -0.000337886, 2.67942, 2.68782, -0.00105287, 1.66973,
   273.884, -1.62928, 0.00393977, -4.38161e-06, 1.85577e-09, 46.0191, -0.00353234, -0.00035993, 2.78447, 2.38215, -0.000999647, 1.7164

      };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons pol4  pol2  pol4
    nParamsElect=13;
    parMomCorrElect[0] = -9.13308;
    parMomCorrElect[1] = 1.40618;
    parMomCorrElect[2] = -0.0630954;
    parMomCorrElect[3] = 0.00123199;
    parMomCorrElect[4] = -0.00000884375;
    parMomCorrElect[5] = 114.165;
    parMomCorrElect[6] = -4.14175;
    parMomCorrElect[7] = 0.0380465;
    parMomCorrElect[8] = -4.75336;
    parMomCorrElect[9] = 0.389213;
    parMomCorrElect[10] = -0.00910909;
    parMomCorrElect[11] = 0.0000934938;
    parMomCorrElect[12] = -0.000000348053;

    typePar = 19;

    } else  if(run == "jul14_w") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 86 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {
      340.901,-1.50393,0.00278068,-0.00000238587,0.0000778523/100000.,119.233,-0.00334163,-0.00093208,8.25497,16.7969,-0.00141246,5.51394,
      359.447,-1.76346,0.00366044,-0.00000356899,0.000134229/100000.,104.042,-0.0033692,-0.000844377,7.23206,19.212,-0.00150658,4.77071,
      359.447,-1.76346,0.00366044,-0.00000356899,0.000134229/100000.,104.042,-0.0033692,-0.000844377,7.23206,19.212,-0.00150658,4.77071,
      359.447,-1.76346,0.00366044,-0.00000356899,0.000134229/100000.,104.042,-0.0033692,-0.000844377,7.23206,19.212,-0.00150658,4.77071,
      365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
      365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
      365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
      365.447,-1.79692,0.00377476,-0.00000375582,0.000145219/100000.,106.245,-0.00332119,-0.000823538,7.47061,22.6797,-0.00154827,5.03414,
      363.522,-1.69303,0.00336942,-0.00000317385,0.000116155/100000.,122.517,-0.00337433,-0.000959434,8.43587,29.2724,-0.00161355,5.63564,
      363.522,-1.69303,0.00336942,-0.00000317385,0.000116155/100000.,122.517,-0.00337433,-0.000959434,8.43587,29.2724,-0.00161355,5.63564,
      363.522,-1.69303,0.00336942,-0.00000317385,0.000116155/100000.,122.517,-0.00337433,-0.000959434,8.43587,29.2724,-0.00161355,5.63564,
      379.355,-1.74201,0.00348639,-0.00000335548,0.000127136/100000.,141.13,-0.00337822,-0.00111639,9.76923,25.7175,-0.00151109,6.5268,
      379.355,-1.74201,0.00348639,-0.00000335548,0.000127136/100000.,141.13,-0.00337822,-0.00111639,9.76923,25.7175,-0.00151109,6.5268,
      379.355,-1.74201,0.00348639,-0.00000335548,0.000127136/100000.,141.13,-0.00337822,-0.00111639,9.76923,25.7175,-0.00151109,6.5268,
      377.386,-1.61336,0.00297582,-0.00000259877,0.0000877779/100000.,147.441,-0.00325714,-0.00111654,10.7727,36.4515,-0.00158623,7.42574,
      368.669,-1.47622,0.00249971,-0.00000193881,0.0000549969/100000.,167.991,-0.0033598,-0.00137807,12.0193,28.8549,-0.00149969,8.04572,
      349.535,-1.24392,0.00174129,-0.00000092267,0.00000571075/100000.,166.974,-0.00325946,-0.00146554,12.9929,43.945,-0.00161646,8.77291,
      375.433,-1.41188,0.00225412,-0.00000162605,0.0000408913/100000.,185.174,-0.00326443,-0.00147595,13.9235,35.1007,-0.00147076,9.54521,
      367.516,-1.24152,0.00163618,-0.000000760779,-0.00000162891/100000.,202.942,-0.00326404,-0.00166685,15.5439,53.5905,-0.00159346,10.6089,
      367.516,-1.24152,0.00163618,-0.000000760779,-0.00000162891/100000.,202.942,-0.00326404,-0.00166685,15.5439,53.5905,-0.00159346,10.6089,
      410.091,-1.48901,0.00236275,-0.00000179646,0.0000538164/100000,211.017,-0.00321069,-0.00173945,16.8942,34.2646,-0.00141895,11.7299,
      338.12,-0.750049,0.0000830358,0.00000112099,-0.0000814397/100000.,249.746,-0.00326935,-0.00196529,18.2082,36.0361,-0.0014173,12.3444
      };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons pol6
    nParamsElect=7;
    parMomCorrElect[0] = 36.8903;
    parMomCorrElect[1] = -3.45134;
    parMomCorrElect[2] = 0.140583;
    parMomCorrElect[3] = -0.00259951;
    parMomCorrElect[4] = 0.0000198486;
    parMomCorrElect[5] = -0.000662151/100000.;
    parMomCorrElect[6] = -0.0000403811/100000.;

    typePar = 14;
  } else  if(run == "jul14_c3") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 86 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {
      363.813,-1.87234,0.00404871,-0.0000041178,0.00016226/100000.,90.8896,-0.00338883,-0.000731713,5.95128,5.35734,-0.000983874,3.73794,
      340.39,-1.78571,0.00390809,-0.00000400218,0.000158149/100000.,82.3474,-0.00344345,-0.000638642,5.1968,5.85834,-0.00108304,3.2622, 
      340.39,-1.78571,0.00390809,-0.00000400218,0.000158149/100000.,82.3474,-0.00344345,-0.000638642,5.1968,5.85834,-0.00108304,3.2622, 
      340.39,-1.78571,0.00390809,-0.00000400218,0.000158149/100000.,82.3474,-0.00344345,-0.000638642,5.1968,5.85834,-0.00108304,3.2622, 
      339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296, 
      339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296, 
      339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296, 
      339.486,-1.78318,0.00390811,-0.00000400864,0.000158707/100000.,81.6281,-0.00343837,-0.000641385,5.19119,5.62965,-0.00106968,3.25296, 
      345.581,-1.79395,0.00389273,-0.00000395957,0.000155694/100000.,84.0997,-0.00341154,-0.000663507,5.42551,5.65421,-0.00105822,3.41813,
      345.581,-1.79395,0.00389273,-0.00000395957,0.000155694/100000.,84.0997,-0.00341154,-0.000663507,5.42551,5.65421,-0.00105822,3.41813,
      345.581,-1.79395,0.00389273,-0.00000395957,0.000155694/100000.,84.0997,-0.00341154,-0.000663507,5.42551,5.65421,-0.00105822,3.41813,
      351.345,-1.81068,0.00390588,-0.00000395676,0.000155283/100000.,88.5021,-0.00344,-0.000700565,5.64477,5.92098,-0.00106608,3.53526, 
      351.345,-1.81068,0.00390588,-0.00000395676,0.000155283/100000.,88.5021,-0.00344,-0.000700565,5.64477,5.92098,-0.00106608,3.53526, 
      351.345,-1.81068,0.00390588,-0.00000395676,0.000155283/100000.,88.5021,-0.00344,-0.000700565,5.64477,5.92098,-0.00106608,3.53526, 
      355.35,-1.82723,0.00394577,-0.00000401542,0.000158839/100000.,92.513,-0.00347466,-0.00074445,5.83361,6.01486,-0.00106656,3.61997, 
      358.09,-1.82336,0.00386925,-0.00000383628,0.000146582/100000.,89.2857,-0.00340005,-0.000704394,5.80144,6.47893,-0.00108469,3.66993, 
      363.225,-1.86139,0.00399638,-0.00000402954,0.00015732/100000.,94.2199,-0.00345562,-0.000745102,5.96009,6.5169,-0.00108865,3.7278,
      367.319,-1.87481,0.00399567,-0.00000398468,0.000153338/100000.,92.8587,-0.00340985,-0.000728312,5.97965,6.31287,-0.00106943,3.78329, 
      365.256,-1.83977,0.00388251,-0.00000384856,0.000147824/100000.,96.1542,-0.00342691,-0.000763212,6.18846,6.87292,-0.00108557,3.89155, 
      365.256,-1.83977,0.00388251,-0.00000384856,0.000147824/100000.,96.1542,-0.00342691,-0.000763212,6.18846,6.87292,-0.00108557,3.89155, 
      377.967,-1.89942,0.0039926,-0.00000393773,0.000150343/100000.,99.3301,-0.00345412,-0.000839982,6.4699,6.84137,-0.00107652,3.98033,
      377.967,-1.89942,0.0039926,-0.00000393773,0.000150343/100000.,99.3301,-0.00345412,-0.000839982,6.4699,6.84137,-0.00107652,3.98033
      };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons pol6
    nParamsElect=7;
    parMomCorrElect[0] = 19.0501;
    parMomCorrElect[1] = -1.78541;
    parMomCorrElect[2] = 0.0799932;
    parMomCorrElect[3] = -0.00181062;
    parMomCorrElect[4] = 0.0000220232;
    parMomCorrElect[5] = -0.000000136325;
    parMomCorrElect[6] = 0.0000334713/100000.;

    typePar = 14;
  } else if(run == "jul14_pe" || run == "aug14_pe") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 86 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {
      395.219,-1.94427,0.00406475,-0.00000403356,0.000156339/100000.,113.782,-0.00340005,-0.000931158,7.49854,6.7844,-0.000981211,4.68952, 
      366.057,-1.81668,0.00377552,-0.00000368103,0.000138881/100000.,101.211,-0.00345251,-0.000818066,6.47664,6.09025,-0.00101831,4.02327, 
      366.057,-1.81668,0.00377552,-0.00000368103,0.000138881/100000.,101.211,-0.00345251,-0.000818066,6.47664,6.09025,-0.00101831,4.02327, 
      366.057,-1.81668,0.00377552,-0.00000368103,0.000138881/100000.,101.211,-0.00345251,-0.000818066,6.47664,6.09025,-0.00101831,4.02327, 
      352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111, 
      352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111, 
      352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111, 
      352.218,-1.79696,0.00382033,-0.00000380159,0.000146201/100000.,90.0212,-0.00345928,-0.00069451,5.62582,6.19342,-0.00108187,3.53111, 
      354.009,-1.86253,0.00408366,-0.00000419045,0.000166054/100000.,83.2562,-0.00343848,-0.000651417,5.28283,6.08154,-0.00110533,3.32809,
      354.009,-1.86253,0.00408366,-0.00000419045,0.000166054/100000.,83.2562,-0.00343848,-0.000651417,5.28283,6.08154,-0.00110533,3.32809,
      354.009,-1.86253,0.00408366,-0.00000419045,0.000166054/100000.,83.2562,-0.00343848,-0.000651417,5.28283,6.08154,-0.00110533,3.32809,
      348.382,-1.8238,0.00396499,-0.0000040237,0.000157397/100000.,82.5202,-0.00345177,-0.000645806,5.19877,6.15786,-0.00112109,3.26955, 
      348.382,-1.8238,0.00396499,-0.0000040237,0.000157397/100000.,82.5202,-0.00345177,-0.000645806,5.19877,6.15786,-0.00112109,3.26955, 
      348.382,-1.8238,0.00396499,-0.0000040237,0.000157397/100000.,82.5202,-0.00345177,-0.000645806,5.19877,6.15786,-0.00112109,3.26955, 
      353.894,-1.87342,0.00412526,-0.00000424506,0.000168473/100000.,83.1282,-0.00346394,-0.000648254,5.20019,6.28797,-0.00113295,3.27027, 
      353.85,-1.86659,0.00409331,-0.00000419413,0.00016577/100000.,83.926,-0.00347166,-0.000660473,5.2443,6.19997,-0.00112843,3.28492, 
      351.858,-1.84545,0.00402381,-0.00000409948,0.000161127/100000.,84.8005,-0.00347937,-0.000657066,5.26311,6.31226,-0.00113253,3.30581,
      356.568,-1.87939,0.00411911,-0.00000421744,0.000166505/100000.,86.025,-0.00349049,-0.00067037,5.32277,5.9717,-0.00111004,3.33111, 
      357.728,-1.87848,0.00410439,-0.00000419094,0.000165047/100000.,85.0026,-0.00344817,-0.000657249,5.3514,6.1102,-0.00111809,3.3878, 
      357.728,-1.87848,0.00410439,-0.00000419094,0.000165047/100000.,85.0026,-0.00344817,-0.000657249,5.3514,6.1102,-0.00111809,3.3878, 
      358.954,-1.877,0.00408277,-0.00000415016,0.000162742/100000.,78.3719,-0.0033126,-0.000622673,5.28907,7.45359,-0.00119627,3.41713,
      358.954,-1.877,0.00408277,-0.00000415016,0.000162742/100000.,78.3719,-0.0033126,-0.000622673,5.28907,7.45359,-0.00119627,3.41713
      };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons pol6
    nParamsElect=7;
    parMomCorrElect[0] = 15.1994;
    parMomCorrElect[1] = -1.15608;
    parMomCorrElect[2] = 0.0509546;
    parMomCorrElect[3] = -0.00122827;
    parMomCorrElect[4] = 0.0000165499;
    parMomCorrElect[5] = -0.000000116424;
    parMomCorrElect[6] = 0.0000333018/100000.;

    typePar = 14;
  } else if(run == "aug14_c7") {
        Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 12;
    nThetaReg    = 22; //from 20 grad till 86 grad
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    momInt1      = 700;
    momInt2      = 2500;
    Double_t par[12*22] = {
      377.29,-1.86591,0.00389833,-0.00000384956,0.000147982/100000.,108.716,-0.00346533,-0.000906983,7.00846,6.17732,-0.000989527,4.2911, 
      362.489,-1.8547,0.00397136,-0.00000398657,0.000154622/100000.,91.5262,-0.00341486,-0.000743405,5.95775,6.29919,-0.00105339,3.71595, 
      362.489,-1.8547,0.00397136,-0.00000398657,0.000154622/100000.,91.5262,-0.00341486,-0.000743405,5.95775,6.29919,-0.00105339,3.71595, 
      362.489,-1.8547,0.00397136,-0.00000398657,0.000154622/100000.,91.5262,-0.00341486,-0.000743405,5.95775,6.29919,-0.00105339,3.71595, 
      349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976, 
      349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976, 
      349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976, 
      349.768,-1.84078,0.00403105,-0.00000412447,0.000162726/100000.,81.3723,-0.00342136,-0.000642395,5.22003,5.82418,-0.00108067,3.27976, 
      338.529,-1.81398,0.00403322,-0.00000418336,0.000167129/100000.,76.5941,-0.00346187,-0.000600158,4.79564,5.78833,-0.00112573,3.00538,
      338.529,-1.81398,0.00403322,-0.00000418336,0.000167129/100000.,76.5941,-0.00346187,-0.000600158,4.79564,5.78833,-0.00112573,3.00538,
      338.529,-1.81398,0.00403322,-0.00000418336,0.000167129/100000.,76.5941,-0.00346187,-0.000600158,4.79564,5.78833,-0.00112573,3.00538,
      337.236,-1.79613,0.00397232,-0.00000409852,0.000162844/100000.,77.4287,-0.00345466,-0.000610697,4.88556,5.68106,-0.00111172,3.06302, 
      337.236,-1.79613,0.00397232,-0.00000409852,0.000162844/100000.,77.4287,-0.00345466,-0.000610697,4.88556,5.68106,-0.00111172,3.06302, 
      337.236,-1.79613,0.00397232,-0.00000409852,0.000162844/100000.,77.4287,-0.00345466,-0.000610697,4.88556,5.68106,-0.00111172,3.06302, 
      344.483,-1.82317,0.00401882,-0.00000414132,0.000164604/100000.,79.2906,-0.00340516,-0.000610728,5.07564,6.62067,-0.00114708,3.23332, 
      349.516,-1.83814,0.00402836,-0.00000412621,0.000162918/100000.,82.5493,-0.00342159,-0.00064238,5.28058,6.18647,-0.00111463,3.34096, 
      349.407,-1.80421,0.00387578,-0.00000388649,0.000150131/100000.,87.0329,-0.00345109,-0.000670834,5.47806,6.88032,-0.00113039,3.45412,
      356.488,-1.84309,0.00397902,-0.00000402071,0.000156785/100000.,88.8678,-0.00343584,-0.000695147,5.65209,6.5324,-0.00110471,3.56262, 
      362.666,-1.85719,0.00397536,-0.00000398646,0.000154446/100000.,93.5664,-0.0034514,-0.000736251,5.93534,6.36717,-0.00107681,3.72505, 
      362.666,-1.85719,0.00397536,-0.00000398646,0.000154446/100000.,93.5664,-0.0034514,-0.000736251,5.93534,6.36717,-0.00107681,3.72505, 
      377.254,-1.89525,0.00397359,-0.00000390074,0.00014792/100000.,95.1167,-0.00338462,-0.000754639,6.17869,7.24617,-0.00110272,3.90527,
      377.254,-1.89525,0.00397359,-0.00000390074,0.00014792/100000.,95.1167,-0.00338462,-0.000754639,6.17869,7.24617,-0.00110272,3.90527
      };
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons pol6
    nParamsElect=7;
    parMomCorrElect[0] = 22.9808;
    parMomCorrElect[1] = -2.53549;
    parMomCorrElect[2] = 0.139744;
    parMomCorrElect[3] = -0.00403039;
    parMomCorrElect[4] = 0.0000628111;
    parMomCorrElect[5] = -0.000000498781;
    parMomCorrElect[6] = 0.000158241/100000.;

    typePar = 14;
  } else if(run == "apr12") {
	Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    nParams      = 7;
    nThetaReg    = 22;
    thetaStep    = 3.;
    thetaMiddle0 = 21.5;
    Double_t par[7*22] = {
      -91.8651,0.0434063,678.858,-1120.65,992.81,-617.704,161.391,
      -98.9038,0.032164,717.974,-1187.02,1071.36,-686.091,185.454,
      -98.6181,0.0339081,716.166,-1183.66,1066.6,-681.832,184.06,
      -94.5398,0.0440444,696.528,-1155.27,1041.17,-663.52,178.263,
      -97.3125,0.0402414,713.659,-1185.34,1075.45,-691.278,187.469,
      -96.829,0.0388894,711.733,-1182.84,1073.93,-691.003,187.68,
      -99.0192,0.0378554,724.509,-1204.06,1095.97,-707.519,192.79,
      -98.4861,0.0402044,722.362,-1200.01,1089.7,-701.493,190.621,
      -100.162,0.0366257,731.766,-1215.2,1105.28,-713.395,194.44,
      -95.8465,0.0382859,708.559,-1184.38,1087.36,-706.755,193.663,
      -69.4872,0.0325539,550.01,-938.753,876.194,-573.516,157.748,
      -63.2317,0.0276368,511.495,-889.159,856.492,-575.452,161.753,
      -64.3656,0.0312026,520.12,-904.467,872.46,-587.034,165.186,
      -68.398,0.0260268,542.253,-938.799,904.784,-610.257,172.378,
      -69.6065,0.0253658,549.615,-949.945,913.837,-615.737,173.836,
      -71.149,0.0273896,558.977,-963.62,923.416,-620.118,174.521,
      -70.434,0.0328006,556.942,-959.752,915.907,-611.974,171.355,
      -73.0875,0.0344005,573.283,-985.193,938.032,-626.329,175.396,
      -73.3104,0.0371747,575.988,-988.623,937.965,-624.185,174.31,
      -74.2248,0.0409931,582.787,-998.393,943.378,-625.427,174.068,
      -78.2616,0.0396539,605.01,-1030.55,967.985,-639.483,177.565,
      -83.4854,0.0227933,629.292,-1062.32,989.642,-651.417,180.689 };
    fillParMomCorrH(nThetaReg*nParams,par);
    typePar = 12;
  } else if(run == "nov02") {
	Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    momInt1=600.;
    momInt2=2000.;
    nThetaReg=7;
    thetaStep=10.;
    thetaMiddle0=23.;
    nParams=10;
    Double_t par[70] = {
      0.0506414,-68.3320,30017.5,87.1936,-0.00405782,-0.000841998,4.52039,5.37248,-0.00121608,2.40080,   // 1st theta interval
      0.056843, -73.7507,30890.6,85.2889,-0.00406022,-0.000823641,4.41504,3.97163,-0.00105256,2.31750,   // 2nd theta interval
      0.0692018,-85.0326,33230.6,86.0512,-0.00408880,-0.000838878,4.39438,4.30892,-0.00110695,2.29046,   // 3rd theta interval
      0.05393,  -72.6496,31345.2,92.3181,-0.00414646,-0.000836083,4.41321,4.73116,-0.00115162,2.30843,   // 4th theta interval
      0.052317, -70.9782,31023.3,94.4471,-0.00416224,-0.000882047,4.56496,5.04969,-0.00117343,2.36712,   // 5th theta interval
      0.0524236,-71.4746,31505.2,97.8558,-0.00416946,-0.000919853,4.75295,5.53359,-0.00120439,2.46908,   // 6th theta interval
      0.0532758,-72.6315,32132.7,100.503,-0.00416881,-0.000962655,4.95385,5.78055,-0.00120932,2.57004 }; // 7th theta interval
    fillParMomCorrH(nThetaReg*nParams,par);
		
    //for electrons
    nParamsElect=5;
    parMomCorrElect[0]=2.29512;
    parMomCorrElect[1]=0.042546;
    parMomCorrElect[2]=-0.00175886;
    parMomCorrElect[3]=2.47346/100000;
    parMomCorrElect[4]=-1.0764/10000000;

    typePar = 1;
  } else if(run=="jan04" || run=="apr07") {
	Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    momInt1=720.;
    momInt2=2000.;
    nThetaReg=7;
    nParams=10;
    thetaStep=10.;
    thetaMiddle0=23.;
     
    Double_t par[70] = {
      0.0404892,-63.2212,31512.3,88.7387,-0.00383467,-0.000831376,4.83807,8.32566,-0.00134054,2.69398,  // 1st theta interval
      0.0392893,-61.7120,31187.0,67.0670,-0.00344743,-0.000770298,4.73321,4.08464,-0.00100119,2.70498,  // 2nd theta interval
      0.0402474,-63.1836,31819.4,95.3604,-0.00391875,-0.000930825,5.09812,4.57863,-0.00103624,2.68722,  // 3rd theta interval
      0.0421779,-64.9599,31849.2,85.6957,-0.00379927,-0.000812833,4.78868,9.17729,-0.00138972,2.64987,  // 4th theta interval
      0.0417347,-64.2374,31548.9,97.1928,-0.00399476,-0.000827159,4.80468,6.59782,-0.00125117,2.63429,  // 4th theta interval
      0.0404892,-63.2212,31512.3,88.7387,-0.00383467,-0.000831377,4.83807,8.32566,-0.00134054,2.69398,  // 6th theta interval
      0.0396831,-62.4906,31610.8,86.3941,-0.00374661,-0.000816029,4.88749,4.51922,-0.00104096,2.73135}; // 7th theta interval
    fillParMomCorrH(nThetaReg*nParams,par);

    //for electrons
    nParamsElect=5;
    parMomCorrElect[0]=2.21734;
    parMomCorrElect[1]=0.0794622;
    parMomCorrElect[2]=-0.00258538;
    parMomCorrElect[3]=3.03896/100000;
    parMomCorrElect[4]=-1.11898/10000000;

    typePar = 1;
  } else if(run=="aug04") {
	Info("setDefaultPar()","Setting dedx correction table for %s",run.Data());
    momInt1=800.;
    momInt2=2000.;
    nThetaReg=7;
    thetaStep=10.;
    thetaMiddle0=23.;
    nParams=10;
    Double_t par[70] = {
      0.0319386,-52.3932,28573.9,65.0824,-0.00340829,-0.000751249,4.91880,7.01293,-0.001202890,2.87892,   // 1st theta interval
      0.0354272,-57.1894,30552.1,73.0873,-0.00349743,-0.000814299,5.21396,5.66944,-0.001081390,2.99888,   // 2nd theta interval
      0.0369929,-59.2908,31560.5,74.1808,-0.00343821,-0.000746043,5.24870,5.80041,-0.001114980,3.16789,   // 3rd theta interval
      0.0381148,-62.3502,34058.2,74.2051,-0.00327892,-0.000714812,5.57025,7.58451,-0.001158990,3.49298,   // 4th theta interval
      0.0405264,-66.2530,36385.7,80.9644,-0.00328172,-0.000733450,6.02206,5.72358,-0.000985652,3.77823,   // 5th theta interval
      0.0442906,-72.6287,39923.8,112.382,-0.00364194,-0.001174410,7.26952,10.0128,-0.001196320,4.12544,   // 6th theta interval
      0.0419634,-70.9195,40627.6,122.213,-0.00366415,-0.001261510,7.75827,9.35643,-0.001134410,4.36657 }; // 7th theta interval
    fillParMomCorrH(nThetaReg*nParams,par);
		
    //for electrons
    nParamsElect=5;
    parMomCorrElect[0]=3.21963;
    parMomCorrElect[1]=-0.0257486;
    parMomCorrElect[2]=0.000407469;
    parMomCorrElect[3]=9.67088e-06;
    parMomCorrElect[4]=-9.93232e-08;
    typePar = 1;
  } else{
	Warning("setDefaultPar()","unknown run %s, no correction table loaded!",run.Data());
  }
  return typePar > 0;
}
void HEnergyLossCorrPar::fillParMomCorrH(Int_t size,Double_t *par) {
  for(Int_t i=0;i<size;i++) parMomCorrH[i] = par[i];
}

Double_t HEnergyLossCorrPar::getDeltaMom(Int_t pId, Double_t mom, Double_t theta) const {
  mom = TMath::Abs(mom);
  if(typePar == 1) return getDeltaMomT1(pId,mom,theta);
  if(typePar == 12) return getDeltaMomT12(pId,mom,theta);
  if(typePar==14 || typePar==19 || typePar==22 || typePar==242 || typePar==243) return getDeltaMomT14192224(pId,mom,theta);
  else return 0.;
}


Double_t HEnergyLossCorrPar::getDeltaMomT1(Int_t pId, Double_t mom, Double_t theta) const {

//  if(theta<15.||theta>87.) return 0.;
  if(theta<15.) theta = 15.;
  if(theta>87.) theta = 87.;

  if(pId==2||pId==3) {
//    if(mom<50.) return 0.;
    if(mom<40.) mom = 40.;
    // 4th order polynomial
    return parMomCorrElect[0]+theta*(parMomCorrElect[1]+theta*(parMomCorrElect[2]+theta*(parMomCorrElect[3]+parMomCorrElect[4]*theta)));
  }
  else {

    Double_t pMom, K[nParams];
    Double_t mPr   = HPhysicsConstants::mass(14);
    Double_t mPart = HPhysicsConstants::mass(pId);
    Int_t   chg   = HPhysicsConstants::charge(pId);
    if(mPart<=0.||(chg==0)) return 0.;

    if(pId==14) pMom=mom;
    else pMom = Sqrt(Power(mPr/mPart*(Sqrt(mom*mom+mPart*mPart)-mPart)+mPr,2)-Power(mPr,2));

//    if(pMom<170.) return 0.;
    if(pMom<170.) pMom = 170.;

    Int_t i=Int_t((theta-thetaMiddle0)/thetaStep+1);
    if(i==nThetaReg) i=nThetaReg-1;

    Int_t j;

    if(pMom<momInt1) j=0; else if(pMom<momInt2) j=3; else j=7;
    if(i==0) {
    // first theta interval
      Double_t thVal = (thetaMiddle0-theta)/thetaStep;
      K[j]   = parMomCorrH[j]   - (parMomCorrH[nParams+j]-parMomCorrH[j])*thVal;
      K[j+1] = parMomCorrH[j+1] - (parMomCorrH[nParams+j+1]-parMomCorrH[j+1])*thVal;
      K[j+2] = parMomCorrH[j+2] - (parMomCorrH[nParams+j+2]-parMomCorrH[j+2])*thVal;
      if(pMom>momInt1 && pMom<momInt2)
        K[j+3] = parMomCorrH[j+3] - (parMomCorrH[nParams+j+3]-parMomCorrH[j+3])*thVal;

//      K[j]=parMomCorrH[j]; K[j+1]=parMomCorrH[j+1]; K[j+2]=parMomCorrH[j+2];
//      if(pMom>momInt1&&pMom<momInt2) K[j+3]=parMomCorrH[j+3];
    }//i=0
    else {
      Int_t bound1 = (i-1)*nParams;
      Int_t bound2 = i*nParams;
      Double_t thInsideInt = (theta-(thetaMiddle0+thetaStep*(i-1)))/thetaStep;
      if(i==nThetaReg-1) {
      // last theta interval
        K[j]=parMomCorrH[bound2+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
        K[j+1]=parMomCorrH[bound2+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
        K[j+2]=parMomCorrH[bound2+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
        if(pMom>momInt1&&pMom<momInt2)
          K[j+3]=parMomCorrH[bound2+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
      }
      else {
      // theta intervals excluding first and last
        K[j]=parMomCorrH[bound1+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
        K[j+1]=parMomCorrH[bound1+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
        K[j+2]=parMomCorrH[bound1+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
        if(pMom>momInt1&&pMom<momInt2)
          K[j+3]=parMomCorrH[bound1+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
      } //else
    } //else

    if(pMom<momInt1) return (K[j]*pMom+K[j+1]+K[j+2]/pMom)*chg*chg;
    if(pMom<momInt2) return (K[j]*exp(K[j+1]*pMom)+K[j+2]*pMom+K[j+3])*chg*chg;
    return (K[j]*exp(K[j+1]*pMom)+K[j+2])*chg*chg;
  } // not leptons

  return 0.;
}


Double_t HEnergyLossCorrPar::binInter(Double_t b,Double_t rs, const Double_t *par) const {
  if(b>=1. || b<=0.) return 0.;

  Double_t b2    = b*b;
  Double_t bg2   = b/Sqrt(1.-b2);
  Double_t lnbg2 = Log(bg2)/b2;
  Double_t odb   = 1./b;
  Double_t osb   = Sqrt(odb);
  Double_t dp1 = par[0]*odb+par[1]*lnbg2+par[2]*osb + par[3]+par[4]*b+par[5]*b2+par[6]*b2*b;
  par += nParams;  // next theta region
  Double_t dp2 = par[0]*odb+par[1]*lnbg2+par[2]*osb + par[3]+par[4]*b+par[5]*b2+par[6]*b2*b;
  return  dp1 + (dp2-dp1)*rs;
}

Double_t HEnergyLossCorrPar::getDeltaMomT12(Int_t pId, Double_t mom, Double_t theta) const {

//  if(theta<15.||theta>87.) return 0.;
  if(theta<15.) theta = 15.;
  if(theta>87.) theta = 87.;

  if(pId==2||pId==3) {
//    if(mom<10.) return 0.;
    if(mom<40.) mom = 40.;

    if(theta<28.)   return 2.96406766784;
    if(theta<50.88) return -7.83473+1.31585*theta-0.0591126*Power(theta,2)+0.00115799*Power(theta,3)-0.00000833136*Power(theta,4);
    if(theta<51.81) return 27.1723-0.479349*theta;
    if(theta<80.28) return 2916.99-261.345*theta+9.7235*Power(theta,2)-0.192166*Power(theta,3)+0.00212771*Power(theta,4)
                               -0.0000125143*Power(theta,5)+0.0000000305468*Power(theta,6);
    return                 0.505533+0.0260891*theta;

  }  // leptons
  else {
  // beamtime = apr12
    Double_t mPart = HPhysicsConstants::mass(pId);
    Int_t    chg2  = HPhysicsConstants::charge(pId);
    if(mPart<=0.|| chg2==0) return 0.;
    chg2 = chg2*chg2;
    Double_t b = mom/Sqrt(mom*mom + mPart*mPart); // Beta
    if(b < 0.11) return 0; // b=0.11 => Pproton = 100 MeV/c
    Int_t i = Int_t((theta-thetaMiddle0)/thetaStep);
    if(i < 0) i = 0;
    if(i >= nThetaReg-1) i = nThetaReg-2;
    Double_t mixtCoef = (theta-(thetaMiddle0+thetaStep*i))/thetaStep;
    return chg2*binInter(b,mixtCoef, parMomCorrH+i*nParams);

  } // not leptons

  return 0.;
}


Double_t HEnergyLossCorrPar::getDeltaMomT14192224(Int_t pId, Double_t mom, Double_t theta) const {

// beamtime = jul14 and aug14 and mar19 and feb22 and feb24
//  if(theta<15.||theta>86.) return 0.;
  if(theta<15.) theta = 15.;
  if(theta>87.) theta = 87.;

  if(pId==2||pId==3) {
    if(mom<40.) mom = 40.;  //return 0.;
    if(typePar == 14)
    // pol6
    return parMomCorrElect[0]+theta*(parMomCorrElect[1]+theta*(parMomCorrElect[2]+theta*(parMomCorrElect[3]+theta*(parMomCorrElect[4]+theta*(parMomCorrElect[5]+parMomCorrElect[6]*theta)))));

    if(typePar == 19) {
    // pol4  pol2  pol4
    Double_t smCor = 0.3;

    if(theta<23.22) return 2.35+smCor;
    if(theta<50.4) return parMomCorrElect[0]+parMomCorrElect[1]*theta+parMomCorrElect[2]*Power(theta,2)+parMomCorrElect[3]*Power(theta,3)+parMomCorrElect[4]*Power(theta,4)+smCor;
    if(theta<53.5) return parMomCorrElect[5]+parMomCorrElect[6]*theta+parMomCorrElect[7]*Power(theta,2)+smCor;
    if(theta<83.3) return parMomCorrElect[8]+parMomCorrElect[9]*theta+parMomCorrElect[10]*Power(theta,2)+parMomCorrElect[11]*Power(theta,3)+parMomCorrElect[12]*Power(theta,4)+smCor;
    return         1.74+smCor;

    } // mar19

    if(typePar == 22) {
    // exp+pol1  exp  pol4

    if(theta<15.0) return 4.03;
    if(theta<31.4) return parMomCorrElect[0]*exp(parMomCorrElect[1]*theta)+parMomCorrElect[2]*theta+parMomCorrElect[3];
    if(theta<44.0) return parMomCorrElect[4]*exp(parMomCorrElect[5]*theta)+parMomCorrElect[6];
    if(theta<87.0) return parMomCorrElect[7]+parMomCorrElect[8]*theta+parMomCorrElect[9]*Power(theta,2)+parMomCorrElect[10]*Power(theta,3)+parMomCorrElect[11]*Power(theta,4);
    return 3.62;

    } // feb22

    if(typePar == 242) {
        // pol4  pol2  pol4
        Double_t smCor = -0.05;

        if(theta<21.) return 3.1+smCor;
        if(theta<49.) return parMomCorrElect[0]+parMomCorrElect[1]*theta+parMomCorrElect[2]*Power(theta,2)+parMomCorrElect[3]*Power(theta,3)+parMomCorrElect[4]*Power(theta,4)+smCor;
        if(theta<53.5) return parMomCorrElect[5]+parMomCorrElect[6]*theta+parMomCorrElect[7]*Power(theta,2)+smCor;
        if(theta<79.) return parMomCorrElect[8]+parMomCorrElect[9]*theta+parMomCorrElect[10]*Power(theta,2)+parMomCorrElect[11]*Power(theta,3)+parMomCorrElect[12]*Power(theta,4)+smCor;
        return         3.27+smCor;
    } //feb24_cc

    if(typePar == 243) {
        // pol4   pol4
        Double_t smCor = -0.05;

        if(theta<21.) return 2.3+smCor;
        if(theta<50.5) return parMomCorrElect[0]+parMomCorrElect[1]*theta+parMomCorrElect[2]*Power(theta,2)+parMomCorrElect[3]*Power(theta,3)+parMomCorrElect[4]*Power(theta,4)+smCor;
//        if(theta<53.5) return parMomCorrElect[5]+parMomCorrElect[6]*theta+parMomCorrElect[7]*Power(theta,2)+smCor;
        if(theta<79.) return parMomCorrElect[8]+parMomCorrElect[9]*theta+parMomCorrElect[10]*Power(theta,2)+parMomCorrElect[11]*Power(theta,3)+parMomCorrElect[12]*Power(theta,4)+smCor;
        return         1.69+smCor;
    } //feb24_auau

  } // leptons
  else {

    Double_t pMom, K[nParams];
    Double_t mPr   = HPhysicsConstants::mass(14);
    Double_t mPart = HPhysicsConstants::mass(pId);
    Int_t   chg   = HPhysicsConstants::charge(pId);
    if(mPart<=0.||(chg==0)) return 0.;

    if(pId==14) pMom=mom;
    else        pMom = Sqrt(Power(mPr/mPart*(Sqrt(mom*mom+mPart*mPart)-mPart)+mPr,2)-Power(mPr,2));

    if(typePar == 14) if(pMom<170.) pMom = 170.;
    if(typePar == 19) if(pMom<100.) pMom = 100.;
    if(typePar == 22) if(pMom<70.) pMom = 70.;
    if(typePar==242||typePar==243) if(pMom<50.) pMom = 50.;

    Int_t i=Int_t((theta-thetaMiddle0)/thetaStep+1);
    if(i==nThetaReg) i=nThetaReg-1;

    Int_t j;

    if(pMom<momInt1) j=0; else if(pMom<momInt2) j=5; else j=9;
    if(i==0) {
    // first theta interval
      Double_t thVal = (thetaMiddle0-theta)/thetaStep;
      K[j]   = parMomCorrH[j]   - (parMomCorrH[nParams+j]-parMomCorrH[j])*thVal;
      K[j+1] = parMomCorrH[j+1] - (parMomCorrH[nParams+j+1]-parMomCorrH[j+1])*thVal;
      K[j+2] = parMomCorrH[j+2] - (parMomCorrH[nParams+j+2]-parMomCorrH[j+2])*thVal;
      if(pMom<momInt2)
        K[j+3] = parMomCorrH[j+3] - (parMomCorrH[nParams+j+3]-parMomCorrH[j+3])*thVal;
      if(pMom<momInt1)
        K[j+4] = parMomCorrH[j+4] - (parMomCorrH[nParams+j+4]-parMomCorrH[j+4])*thVal;

    }//i=0
    else {
      Int_t bound1 = (i-1)*nParams;
      Int_t bound2 = i*nParams;
      Double_t thInsideInt = (theta-(thetaMiddle0+thetaStep*(i-1)))/thetaStep;
      if(i==nThetaReg-1) {
      // last theta interval
        K[j]=parMomCorrH[bound2+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
        K[j+1]=parMomCorrH[bound2+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
        K[j+2]=parMomCorrH[bound2+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
        if(pMom<momInt2)
          K[j+3]=parMomCorrH[bound2+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
        if(pMom<momInt1)
          K[j+4]=parMomCorrH[bound2+j+4] + (parMomCorrH[bound2+j+4]-parMomCorrH[bound1+j+4])*thInsideInt;
      }
      else {
      // theta intervals excluding first and last
        K[j]=parMomCorrH[bound1+j] + (parMomCorrH[bound2+j]-parMomCorrH[bound1+j])*thInsideInt;
        K[j+1]=parMomCorrH[bound1+j+1] + (parMomCorrH[bound2+j+1]-parMomCorrH[bound1+j+1])*thInsideInt;
        K[j+2]=parMomCorrH[bound1+j+2] + (parMomCorrH[bound2+j+2]-parMomCorrH[bound1+j+2])*thInsideInt;
        if(pMom<momInt2)
          K[j+3]=parMomCorrH[bound1+j+3] + (parMomCorrH[bound2+j+3]-parMomCorrH[bound1+j+3])*thInsideInt;
        if(pMom<momInt1)
          K[j+4]=parMomCorrH[bound1+j+4] + (parMomCorrH[bound2+j+4]-parMomCorrH[bound1+j+4])*thInsideInt;
      } //else
    } // i!=0

    if(pMom<momInt1) return (K[j]+pMom*(K[j+1]+pMom*(K[j+2]+pMom*(K[j+3]+K[j+4]*pMom))))*chg*chg;
    if(pMom<momInt2) return (K[j]*exp(K[j+1]*pMom)+K[j+2]*pMom+K[j+3])*chg*chg;
    return (K[j]*exp(K[j+1]*pMom)+K[j+2])*chg*chg;
  } // not leptons

  return 0.;
}