#include "Bin.h"
#include <iostream>
using namespace std;
using namespace Bin;

//The latest binnings... 11-June-2002 JLK

void Bin::initPtAry(TArrayD* ary, int type)
{

  //Jon's binning as of 04-Dec-2002
  double ptAry0[]=
    {
    1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 
    2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7,
    4.0, 4.3, 4.7, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 12.5, 14.5 
    };
  int n0 = (int)sizeof(ptAry0)/sizeof(double);

  //Bum's binning
  double ptAry1[]=
    {
    1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4,
    2.6, 2.8, 3.0, 3.35, 3.8, 4.4, 5.1, 6.0, 7.0,
    8.0, 9.0, 10.0, 12.0, 14.0, 16.0
    };
  int n1 = (int)sizeof(ptAry1)/sizeof(double);

  //xT 130GeV Bum binning, where xT=2*pT/Sqrt(s)
  double ptAry2[]=
    {
    0.0246154, 0.0261538, 0.0276923, 0.0292308, 0.0307692,
    0.0323077, 0.0338462, 0.0353846, 0.0369231, 0.0400000,
    0.0430769, 0.0461538, 0.0515385, 0.0584615, 0.0676923,
    0.0784615, 0.0923077, 0.1076920, 0.1230770, 0.1384620,
    0.1538460, 0.1846150, 0.2153850, 0.2461540
    };
  int n2 = (int)sizeof(ptAry2)/sizeof(double);


  switch(type){
  case 0: ary->Set(n0,ptAry0); break;
  case 1: ary->Set(n1,ptAry1); break;
  case 2: ary->Set(n2,ptAry2); break;
  default:
    cout << "ERROR. type " << type << "out of range" << endl; exit(1);
  }
  
}