#ifndef CONSTANTS_HH
#define CONSTANTS_HH

#include "TCut.h"
#include "RooRealVar.h"

using namespace RooFit;

//-----Define some constants

#define DSMASS 2.01029
#define PIONMASS 0.13957018
#define D0MASS 1.86486
#define DSmD0MASS 0.145421
#define D1MASS 2.4214
#define D1WIDTH 0.0274
#define D2sMASS 2.4626
#define D2sWIDTH 0.049
#define D1brMASS 2.427
#define D1brWIDTH 0.384
#define DJs_2650_MASS 2.4604
#define DJs_2650_WIDTH 0.14
#define DJs_2760_MASS 2.7611
#define DJs_2760_WIDTH 0.0744
#define DJ_2580_MASS 2.5795
#define DJ_2580_WIDTH 0.1775
#define DJ_2740_MASS 2.737
#define DJ_2740_WIDTH 0.0732
#define DJ_3000_MASS 2.9718
#define DJ_3000_WIDTH 0.1881

//------Define fiducial volumes for fitting and selecting data

// m(D*) - m(D0)
int bin_DeltaM = 210;
double min_DeltaM = 0.1405;
double max_DeltaM = 0.160;

double min_DeltaM_fit = 0.1405;
double max_DeltaM_fit = 0.160;

// LogIP
int bin_PiMipBV = 110;
double min_PiMipBV = -8.0;
double max_PiMipBV = 3.0;

double min_PiMipBV_fit = -8.0;
double max_PiMipBV_fit = 3.0;

// m(D*pi)
int bin_DsPi = 130;
double min_DsPi = 2.2;
double max_DsPi = 3.5;

double min_DsPi_fit= 2.2;
double max_DsPi_fit= 3.5;

// m(D*pi) resolution
int bin_res = 80;
double min_res = -0.02;
double max_res = 0.02;

double min_res_fit = -0.02;
double max_res_fit = 0.02;

// D1(2420)^0
int bin_D1 = 150;
double min_D1 = 2.2;
double max_D1 = 2.8;

// D2*(2460)^0
int bin_D2s = 150;
double min_D2s = 2.2;
double max_D2s = 2.8;

// D1(2430)^0
int bin_D1br = 150;
double min_D1br = 2.2;
double max_D1br = 2.8;

//-----Create constant roorealvars
RooRealVar DsMass("DsMass","M(D^{*})",DSMASS,"GeV/c^{2}");
RooRealVar PionMass("PionMass","M(#pi)",PIONMASS,"GeV/c^{2}");
RooRealVar spin0("spin0","D^{**} decay wave",0,"");
RooRealVar spin2("spin2","D^{**} decay wave",2,"");
RooRealVar mg("mg","#mu",0.0,"GeV/c^{2}");
RooRealVar r("r","D^{**} radius",4.0,"(GeV/c)^{-1}");

//-----Define monk pion preselection (before TMVA training)
TCut presel = "isPion==1 && isTagPWrongCharge==0 && nTagParticles>0 && PiMipPU>-99998 && PiMipPV>-99998 && PiMipBV>-99998 && PiMipBV<=3.0 && PIDNNPion>0.1 && DsPiMass>=2.2 && DsPiMass<=3.5";

//---Define BDTA cuts
double BDTA_VeryLooseCuts = -999.0;
double BDTA_LooseCuts_PiMBdphi = -999.0;
double BDTA_LooseCuts_PiMBdeta = -999.0;
double BDTA_LooseCuts_PiMBdQ = -999.0;
double BDTA_LooseCuts_PiMDsdphi = -999.0;
double BDTA_LooseCuts_PiMDsdeta = -999.0;
double BDTA_AllCuts = 0.1021;

#endif