#ifndef histMk_3D_setup_h #define histMk_3D_setup_h { // Setup file for histMK_3D_*.C //########################################################### // Change these values to meet your condition //########################################################### // Run1, Run2, Run3, Run4 or all Run1-4 if run==0 int run = 3; // Dlnu FF slope int newDlnuFF = 1; double rho_D_sq = 1.47292; //1.43262; TString DlnuFF("rho0000"); //"rho0000" // Use new parametrization or not int newSlope = 1; // D*lnu FF // use new FF parameters or not int newDstarFF = 1; TString DlnuFFMod(""); //"", "RhoMax", "RhoMin", "R1Max", etc... double R1_new = 1.396; //1.396 +- 0.112; double R2_new = 0.885; //0.885 +- 0.075; double rhoSq_new = 1.145; //1.145 +- 0.124; // D+ semilep BF int newDsemilepBF = 1; double DKlnu_new = 0.0871; double DKstarlnu_new = 0.0556; //############################################################ // tracking and PID correction int correctTrk = 1; // tracking efficiency and PID correction modification factors // modified correction = Corr * Modif double trkCorrModif = 1; double PidCorrModif = 1; // 3 variables TString varx("DlLeptMom"); TString vary("DMom"); TString varz("DLcosBY"); // binning enum { nbinsx = 4, nbinsy = 6, nbinsz = 5 }; Float_t xbins[5] = {0.8, 1.2, 1.6, 2.0, 2.4}; Float_t ybins[7] = {0.2, 0.6, 1.0, 1.4, 1.8, 2.2, 2.6}; Float_t zbins[6] = {-10, -2.5, -1.1, 0.0, 1.1, 5}; //Float_t xbins[6] = {1.0, 1.4, 1.6, 1.8, 2.0, 2.4}; //Float_t ybins[6] = {0.5, 0.9, 1.3, 1.7, 2.1, 2.5}; // D decay mode to select int DDecMod_D0KPi = 1; int DDecMod_DchKPiPi = 5; //############################################################ // D mass range float DMass_D0KPi1 = 1.816; // lower edge of sideband float DMass_D0KPi2 = 1.840; // lower edge of D mass peak float DMass_D0KPi3 = 1.888; // upper edge of D mass peak float DMass_D0KPi4 = 1.912; // upper edge of sideband float DMass_DchKPiPi1 = 1.821; // lower edge of sideband float DMass_DchKPiPi2 = 1.845; // lower edge of D mass peak float DMass_DchKPiPi3 = 1.893; // upper edge of D mass peak float DMass_DchKPiPi4 = 1.917; // upper edge of sideband // Input root-ntuple files TString baseDir_Run1("/home1x/OtherMounts/hep03/khamano/R16b/Run1"); TString inFile_Run1_OnPeak("BToDlnu-Run1-OnPeak-R16b"); int nStart_Run1_OnPeak = 1; int nStop_Run1_OnPeak = 440; TString inFile_Run1_OffPeak("BToDlnu-Run1-OffPeak-R16b"); int nStart_Run1_OffPeak = 1; int nStop_Run1_OffPeak = 26; TString inFile_Run1_BpBm("SP-1235-BToDlnu-Run1-R16b"); int nStart_Run1_BpBm = 1; int nStop_Run1_BpBm = 347; TString inFile_Run1_B0B0("SP-1237-BToDlnu-Run1-R16b"); int nStart_Run1_B0B0 = 1; int nStop_Run1_B0B0 = 372; TString inFile_Run1_ccbar("SP-1005-BToDlnu-Run1-R16b"); int nStart_Run1_ccbar = 1; int nStop_Run1_ccbar = 95; TString inFile_Run1_uds("SP-998-BToDlnu-Run1-R16b"); int nStart_Run1_uds = 1; int nStop_Run1_uds = 48; TString inFile_Run1_tautau("SP-3429-BToDlnu-Run1-R16b"); int nStart_Run1_tautau = 1; int nStop_Run1_tautau = 7; TString baseDir_Run2("/home1x/OtherMounts/hep03/khamano/R16b/Run2"); TString inFile_Run2_OnPeak("BToDlnu-Run2-OnPeak-R16b"); int nStart_Run2_OnPeak = 1; int nStop_Run2_OnPeak = 1569; //1569; TString inFile_Run2_OffPeak("BToDlnu-Run2-OffPeak-R16b"); int nStart_Run2_OffPeak = 1; int nStop_Run2_OffPeak = 84; //84; TString inFile_Run2_BpBm("SP-1235-BToDlnu-Run2-R16b"); int nStart_Run2_BpBm = 1; int nStop_Run2_BpBm = 1952; //1952; TString inFile_Run2_B0B0("SP-1237-BToDlnu-Run2-R16b"); int nStart_Run2_B0B0 = 1; int nStop_Run2_B0B0 = 1509; //1509; TString inFile_Run2_ccbar("SP-1005-BToDlnu-Run2-R16b"); int nStart_Run2_ccbar = 1; int nStop_Run2_ccbar = 492; //492; TString inFile_Run2_uds("SP-998-BToDlnu-Run2-R16b"); int nStart_Run2_uds = 1; int nStop_Run2_uds = 340; //340; TString inFile_Run2_tautau("SP-3429-BToDlnu-Run2-R16b"); int nStart_Run2_tautau = 1; int nStop_Run2_tautau = 35; //35; TString baseDir_Run3("/home1x/OtherMounts/hep03/khamano/R16b/Run3"); TString inFile_Run3_OnPeak("BToDlnu-Run3-OnPeak-R16b"); int nStart_Run3_OnPeak = 1; int nStop_Run3_OnPeak = 877; //877; //int nStop_Run3_OnPeak = 5; TString inFile_Run3_OffPeak("BToDlnu-Run3-OffPeak-R16b"); int nStart_Run3_OffPeak = 1; int nStop_Run3_OffPeak = 33; //33; //int nStop_Run3_OffPeak = 5; TString inFile_Run3_BpBm("SP-1235-BToDlnu-Run3-R16b"); int nStart_Run3_BpBm = 1; //1; int nStop_Run3_BpBm = 1193; //1193; //int nStop_Run3_BpBm = 5; TString inFile_Run3_B0B0("SP-1237-BToDlnu-Run3-R16b"); int nStart_Run3_B0B0 = 1; //1; int nStop_Run3_B0B0 = 855; //855; //int nStop_Run3_B0B0 = 5; TString inFile_Run3_ccbar("SP-1005-BToDlnu-Run3-R16b"); int nStart_Run3_ccbar = 1; int nStop_Run3_ccbar = 307; //307; //int nStop_Run3_ccbar = 5; TString inFile_Run3_uds("SP-998-BToDlnu-Run3-R16b"); int nStart_Run3_uds = 1; int nStop_Run3_uds = 282; //282; //int nStop_Run3_uds = 5; TString inFile_Run3_tautau("SP-3429-BToDlnu-Run3-R16b"); int nStart_Run3_tautau = 1; int nStop_Run3_tautau = 28; //28; //int nStop_Run3_tautau = 5; TString baseDir_Run4("/home1x/OtherMounts/hep03/khamano/R16b/Run4"); TString inFile_Run4_OnPeak("BToDlnu-Run4-OnPeak-R16b"); int nStart_Run4_OnPeak = 1; int nStop_Run4_OnPeak = 2898; TString inFile_Run4_OffPeak("BToDlnu-Run4-OffPeak-R16b"); int nStart_Run4_OffPeak = 1; int nStop_Run4_OffPeak = 144; TString inFile_Run4_BpBm("SP-1235-BToDlnu-Run4-R16b"); int nStart_Run4_BpBm = 1; int nStop_Run4_BpBm = 4488; TString inFile_Run4_B0B0("SP-1237-BToDlnu-Run4-R16b"); int nStart_Run4_B0B0 = 1; int nStop_Run4_B0B0 = 4505; TString inFile_Run4_ccbar("SP-1005-BToDlnu-Run4-R16b"); int nStart_Run4_ccbar = 1; int nStop_Run4_ccbar = 1298; TString inFile_Run4_uds("SP-998-BToDlnu-Run4-R16b"); int nStart_Run4_uds = 1; int nStop_Run4_uds = 763; TString inFile_Run4_tautau("SP-3429-BToDlnu-Run4-R16b"); int nStart_Run4_tautau = 1; int nStop_Run4_tautau = 74; // Name of the root tree in the input root files TString ntpName("ntp2"); // Luminosities double lumi_Run1_OnPeak = 17.801; double lumi_Run1_OffPeak = 2.330; double lumi_Run1_BpBm = 51.764; double lumi_Run1_B0B0 = 54.785; double lumi_Run1_ccbar = 16.420; double lumi_Run1_uds = 14.434; double lumi_Run1_tautau = 23.509; double lumi_Run2_OnPeak = 59.418; double lumi_Run2_OffPeak = 6.836; double lumi_Run2_BpBm = 279.12; double lumi_Run2_B0B0 = 213.19; double lumi_Run2_ccbar = 76.708; double lumi_Run2_uds = 84.326; double lumi_Run2_tautau = 106.18; double lumi_Run3_OnPeak = 30.630; double lumi_Run3_OffPeak = 2.394; double lumi_Run3_BpBm = 164.75; double lumi_Run3_B0B0 = 116.49; double lumi_Run3_ccbar = 44.385; double lumi_Run3_uds = 59.349; double lumi_Run3_tautau = 74.311; double lumi_Run4_OnPeak = 99.357; double lumi_Run4_OffPeak = 9.923; double lumi_Run4_BpBm = 616.80; double lumi_Run4_B0B0 = 611.71; double lumi_Run4_ccbar = 189.69; double lumi_Run4_uds = 165.81; double lumi_Run4_tautau = 164.79; double lumi_All_OnPeak = lumi_Run1_OnPeak + lumi_Run2_OnPeak + lumi_Run3_OnPeak + lumi_Run4_OnPeak; double lumi_All_OffPeak = lumi_Run1_OffPeak + lumi_Run2_OffPeak + lumi_Run3_OffPeak + lumi_Run4_OffPeak; double lumi_All_BpBm = lumi_Run1_BpBm + lumi_Run2_BpBm + lumi_Run3_BpBm + lumi_Run4_BpBm; double lumi_All_B0B0 = lumi_Run1_B0B0 + lumi_Run2_B0B0 + lumi_Run3_B0B0 + lumi_Run4_B0B0; double lumi_All_ccbar = lumi_Run1_ccbar + lumi_Run2_ccbar + lumi_Run3_ccbar + lumi_Run4_ccbar; double lumi_All_uds = lumi_Run1_uds + lumi_Run2_uds + lumi_Run3_uds + lumi_Run4_uds; double lumi_All_tautau = lumi_Run1_tautau + lumi_Run2_tautau + lumi_Run3_tautau + lumi_Run4_tautau; // Conversion factors double c_Run1_OffPeak = lumi_Run1_OnPeak / lumi_Run1_OffPeak; double c_Run1_BpBm = lumi_Run1_OnPeak / lumi_Run1_BpBm; double c_Run1_B0B0 = lumi_Run1_OnPeak / lumi_Run1_B0B0; double c_Run1_ccbar = lumi_Run1_OnPeak / lumi_Run1_ccbar; double c_Run1_uds = lumi_Run1_OnPeak / lumi_Run1_uds; double c_Run1_tautau = lumi_Run1_OnPeak / lumi_Run1_tautau; double c_Run2_OffPeak = lumi_Run2_OnPeak / lumi_Run2_OffPeak; double c_Run2_BpBm = lumi_Run2_OnPeak / lumi_Run2_BpBm; double c_Run2_B0B0 = lumi_Run2_OnPeak / lumi_Run2_B0B0; double c_Run2_ccbar = lumi_Run2_OnPeak / lumi_Run2_ccbar; double c_Run2_uds = lumi_Run2_OnPeak / lumi_Run2_uds; double c_Run2_tautau = lumi_Run2_OnPeak / lumi_Run2_tautau; double c_Run3_OffPeak = lumi_Run3_OnPeak / lumi_Run3_OffPeak; double c_Run3_BpBm = lumi_Run3_OnPeak / lumi_Run3_BpBm; double c_Run3_B0B0 = lumi_Run3_OnPeak / lumi_Run3_B0B0; double c_Run3_ccbar = lumi_Run3_OnPeak / lumi_Run3_ccbar; double c_Run3_uds = lumi_Run3_OnPeak / lumi_Run3_uds; double c_Run3_tautau = lumi_Run3_OnPeak / lumi_Run3_tautau; double c_Run4_OffPeak = lumi_Run4_OnPeak / lumi_Run4_OffPeak; double c_Run4_BpBm = lumi_Run4_OnPeak / lumi_Run4_BpBm; double c_Run4_B0B0 = lumi_Run4_OnPeak / lumi_Run4_B0B0; double c_Run4_ccbar = lumi_Run4_OnPeak / lumi_Run4_ccbar; double c_Run4_uds = lumi_Run4_OnPeak / lumi_Run4_uds; double c_Run4_tautau = lumi_Run4_OnPeak / lumi_Run4_tautau; double c_All_OffPeak = lumi_All_OnPeak / lumi_All_OffPeak; double c_All_BpBm = lumi_All_OnPeak / lumi_All_BpBm; double c_All_B0B0 = lumi_All_OnPeak / lumi_All_B0B0; double c_All_ccbar = lumi_All_OnPeak / lumi_All_ccbar; double c_All_uds = lumi_All_OnPeak / lumi_All_uds; double c_All_tautau = lumi_All_OnPeak / lumi_All_tautau; } #endif