#ifndef plotMk_3D_setup_h #define plotMk_3D_setup_h { // Setup file for plotMK_3D_*.C //########################################################### // Change these values to meet your condition //########################################################### TString runNum("Run3"); //Run1,Run2, Run3, Run4 or RunAll TString corrEff("effCorr"); // "noCorr" or "effCorr" TString corrDstarFF("newFF"); // "newFF" or "oldFF" TString DlnuFF("rho147292"); // Dlnu FF slope of input histo int useOneHist = 1; // Use the above input hist or not double rho2 = 1.47292; //0.83456; TString DlnuFF2("rho147292"); // Dlnu FF slope for output file name TString Slope("New"); // Slope parametrization, "Old" or "New" TString corrBF("newBF"); // "oldBF" or "newBF" TString histDir("histogramsFit"); //"histogramsFit" or "valHisto1" // Fitted branching fractions double B_B1 = 0.0262713; //0.0210 double B_B2 = 0.0581354; //0.0560 double B_B31 = 0.00812992; //0.0020 //double B_B31 = B_B32*20/56; double B_B32 = 0.00; //0.0056 double B_B33 = 0.000647493; //0.0037 //double B_B33 = B_B32*37/56; double B_B34 = 0.00100122; //0.0037 //double B_B34 = B_B32*37/56; double B_B41 = 0.00; //0.0060 //double B_B41 = 0.00849; double B_B42 = 0.00816908; //0.0020 //double B_B42 = B_B41*67/54; // Fitted background double DiffB_D0 = 0.962483; double CascL_D0 = 0.930998; double LMisID_D0 = 0.228374; double RealD_D0 = 0.995462; double DiffB_Dch = 0.707311; double CascL_Dch = 0.753825; double LMisID_Dch = 0.910783; double RealD_Dch = 1.17802; // Fitted f_+0, t_0+ double f_p0 = 1.04069; double t_0p = 1; //0.92081; // = 1/1.086 double c_Dp0 = 0.971755; // 3 variables TString varx("LeptMom"); TString vary("DMom"); TString varz("cosBY"); // 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}; // Set max value of plot bool setMax = false; // for projection plots int xmax = 40000; int ymax = 40000; int zmax = 40000; int xmin = -100; int ymin = -100; int zmin = -100; bool requireSetMax = false; // for bin by bin plots int xmaxBin = 3000; int ymaxBin = 3000; int zmaxBin = 3000; int xminBin = -100; int yminBin = -100; int zminBin = -100; // Name base of input output files TString nameBase("3D_"+runNum+"_"+corrEff); TString oNameBase(nameBase+"_"+corrDstarFF+"_"+DlnuFF2+Slope); // name for output file // Old branching fractions double B_B1_SP6 = 0.0210; double B_B2_SP6 = 0.0560; double B_B31_SP6 = 0.0020; double B_B32_SP6 = 0.0056; double B_B33_SP6 = 0.0037; double B_B34_SP6 = 0.0037; double B_B41_SP6 = 0.0060; double B_B42_SP6 = 0.0020; // Input files enum { nFiles_D0KPi = 32, nFiles_DchKPiPi = 31 }; //TString inFile_OnPeak_data_D0KPi //("valHisto2/hist_AllMC_3D_Run3_noCorr_oldFF_rho0000Old_D0KPi.root"); TString inFile_OnPeak_data_D0KPi (histDir+"/histo_"+nameBase+"_OnPeak_data_D0KPi.root"); TString inFile_OffPeak_data_D0KPi (histDir+"/histo_"+nameBase+"_OffPeak_data_D0KPi.root"); TString inFile_BpBm_Dlnu_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_Dlnu_"+DlnuFF+Slope+"_D0KPi.root"); TString inFile_BpBm_Dstar_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_Dstar_"+corrDstarFF+"_D0KPi.root"); TString inFile_B0B0_Dstar_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_Dstar_"+corrDstarFF+"_D0KPi.root"); TString inFile_BpBm_D0star_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_D0star_D0KPi.root"); TString inFile_BpBm_D1_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_D1_D0KPi.root"); TString inFile_BpBm_D1prime_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_D1prime_D0KPi.root"); TString inFile_BpBm_D2star_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_D2star_D0KPi.root"); TString inFile_BpBm_DPi_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_DPi_D0KPi.root"); TString inFile_BpBm_DstarPi_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_DstarPi_D0KPi.root"); TString inFile_BpBm_DiffB_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_DiffB_D0KPi.root"); TString inFile_BpBm_CascL_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_CascL_"+corrBF+"_D0KPi.root"); TString inFile_BpBm_LMisID_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_LMisID_D0KPi.root"); TString inFile_BpBm_Comb_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_Comb_D0KPi.root"); TString inFile_B0B0_D0star_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_D0star_D0KPi.root"); TString inFile_B0B0_D1_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_D1_D0KPi.root"); TString inFile_B0B0_D1prime_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_D1prime_D0KPi.root"); TString inFile_B0B0_D2star_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_D2star_D0KPi.root"); TString inFile_B0B0_DPi_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_DPi_D0KPi.root"); TString inFile_B0B0_DstarPi_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_DstarPi_D0KPi.root"); TString inFile_B0B0_DiffB_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_DiffB_D0KPi.root"); TString inFile_B0B0_CascL_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_CascL_"+corrBF+"_D0KPi.root"); TString inFile_B0B0_LMisID_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_LMisID_D0KPi.root"); TString inFile_B0B0_Comb_D0KPi (histDir+"/histo_"+nameBase+"_B0B0_Comb_D0KPi.root"); TString inFile_ccbar_RealD_D0KPi (histDir+"/histo_"+nameBase+"_ccbar_RealD_"+corrBF+"_D0KPi.root"); TString inFile_ccbar_Cont_D0KPi (histDir+"/histo_"+nameBase+"_ccbar_Cont_D0KPi.root"); TString inFile_uds_Cont_D0KPi (histDir+"/histo_"+nameBase+"_uds_Cont_D0KPi.root"); TString inFile_tautau_Cont_D0KPi (histDir+"/histo_"+nameBase+"_tautau_Cont_D0KPi.root"); //TString inFile_OnPeak_data_DchKPiPi //("valHisto2/hist_AllMC_3D_Run3_noCorr_oldFF_rho0000Old_DchKPiPi.root"); TString inFile_OnPeak_data_DchKPiPi (histDir+"/histo_"+nameBase+"_OnPeak_data_DchKPiPi.root"); TString inFile_OffPeak_data_DchKPiPi (histDir+"/histo_"+nameBase+"_OffPeak_data_DchKPiPi.root"); TString inFile_BpBm_D0star_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_D0star_DchKPiPi.root"); TString inFile_BpBm_D1_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_D1_DchKPiPi.root"); TString inFile_BpBm_D1prime_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_D1prime_DchKPiPi.root"); TString inFile_BpBm_D2star_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_D2star_DchKPiPi.root"); TString inFile_BpBm_DPi_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_DPi_DchKPiPi.root"); TString inFile_BpBm_DstarPi_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_DstarPi_DchKPiPi.root"); TString inFile_BpBm_DiffB_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_DiffB_DchKPiPi.root"); TString inFile_BpBm_CascL_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_CascL_"+corrBF+"_DchKPiPi.root"); TString inFile_BpBm_LMisID_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_LMisID_DchKPiPi.root"); TString inFile_BpBm_Comb_DchKPiPi (histDir+"/histo_"+nameBase+"_BpBm_Comb_DchKPiPi.root"); TString inFile_B0B0_Dlnu_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_Dlnu_"+DlnuFF+Slope+"_DchKPiPi.root"); TString inFile_B0B0_Dstar_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_Dstar_"+corrDstarFF+"_DchKPiPi.root"); TString inFile_B0B0_D0star_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_D0star_DchKPiPi.root"); TString inFile_B0B0_D1_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_D1_DchKPiPi.root"); TString inFile_B0B0_D1prime_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_D1prime_DchKPiPi.root"); TString inFile_B0B0_D2star_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_D2star_DchKPiPi.root"); TString inFile_B0B0_DPi_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_DPi_DchKPiPi.root"); TString inFile_B0B0_DstarPi_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_DstarPi_DchKPiPi.root"); TString inFile_B0B0_DiffB_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_DiffB_DchKPiPi.root"); TString inFile_B0B0_CascL_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_CascL_"+corrBF+"_DchKPiPi.root"); TString inFile_B0B0_LMisID_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_LMisID_DchKPiPi.root"); TString inFile_B0B0_Comb_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_Comb_DchKPiPi.root"); TString inFile_ccbar_RealD_DchKPiPi (histDir+"/histo_"+nameBase+"_ccbar_RealD_"+corrBF+"_DchKPiPi.root"); TString inFile_ccbar_Cont_DchKPiPi (histDir+"/histo_"+nameBase+"_ccbar_Cont_DchKPiPi.root"); TString inFile_uds_Cont_DchKPiPi (histDir+"/histo_"+nameBase+"_uds_Cont_DchKPiPi.root"); TString inFile_tautau_Cont_DchKPiPi (histDir+"/histo_"+nameBase+"_tautau_Cont_DchKPiPi.root"); TString inFile_BpBm_Dlnu_0th_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_Dlnu_0th"+Slope+"_D0KPi.root"); TString inFile_BpBm_Dlnu_1st_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_Dlnu_1st"+Slope+"_D0KPi.root"); TString inFile_BpBm_Dlnu_2nd_D0KPi (histDir+"/histo_"+nameBase+"_BpBm_Dlnu_2nd"+Slope+"_D0KPi.root"); TString inFile_B0B0_Dlnu_0th_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_Dlnu_0th"+Slope+"_DchKPiPi.root"); TString inFile_B0B0_Dlnu_1st_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_Dlnu_1st"+Slope+"_DchKPiPi.root"); TString inFile_B0B0_Dlnu_2nd_DchKPiPi (histDir+"/histo_"+nameBase+"_B0B0_Dlnu_2nd"+Slope+"_DchKPiPi.root"); // Input files TString inFile_D0KPi[nFiles_D0KPi]; inFile_D0KPi[0] = inFile_OnPeak_data_D0KPi; inFile_D0KPi[1] = inFile_OffPeak_data_D0KPi; inFile_D0KPi[2] = inFile_BpBm_Dlnu_D0KPi; inFile_D0KPi[3] = inFile_BpBm_D0star_D0KPi; inFile_D0KPi[4] = inFile_B0B0_D0star_D0KPi; inFile_D0KPi[5] = inFile_BpBm_D1_D0KPi; inFile_D0KPi[6] = inFile_B0B0_D1_D0KPi; inFile_D0KPi[7] = inFile_BpBm_D1prime_D0KPi; inFile_D0KPi[8] = inFile_B0B0_D1prime_D0KPi; inFile_D0KPi[9] = inFile_BpBm_D2star_D0KPi; inFile_D0KPi[10] = inFile_B0B0_D2star_D0KPi; inFile_D0KPi[11] = inFile_BpBm_DPi_D0KPi; inFile_D0KPi[12] = inFile_B0B0_DPi_D0KPi; inFile_D0KPi[13] = inFile_BpBm_DstarPi_D0KPi; inFile_D0KPi[14] = inFile_B0B0_DstarPi_D0KPi; inFile_D0KPi[15] = inFile_BpBm_DiffB_D0KPi; inFile_D0KPi[16] = inFile_B0B0_DiffB_D0KPi; inFile_D0KPi[17] = inFile_BpBm_CascL_D0KPi; inFile_D0KPi[18] = inFile_B0B0_CascL_D0KPi; inFile_D0KPi[19] = inFile_BpBm_LMisID_D0KPi; inFile_D0KPi[20] = inFile_B0B0_LMisID_D0KPi; inFile_D0KPi[21] = inFile_BpBm_Comb_D0KPi; inFile_D0KPi[22] = inFile_B0B0_Comb_D0KPi; inFile_D0KPi[23] = inFile_ccbar_RealD_D0KPi; inFile_D0KPi[24] = inFile_ccbar_Cont_D0KPi; inFile_D0KPi[25] = inFile_uds_Cont_D0KPi; inFile_D0KPi[26] = inFile_tautau_Cont_D0KPi; inFile_D0KPi[27] = inFile_B0B0_Dstar_D0KPi; inFile_D0KPi[28] = inFile_BpBm_Dstar_D0KPi; inFile_D0KPi[29] = inFile_BpBm_Dlnu_0th_D0KPi; inFile_D0KPi[30] = inFile_BpBm_Dlnu_1st_D0KPi; inFile_D0KPi[31] = inFile_BpBm_Dlnu_2nd_D0KPi; TString inFile_DchKPiPi[nFiles_DchKPiPi]; inFile_DchKPiPi[0] = inFile_OnPeak_data_DchKPiPi; inFile_DchKPiPi[1] = inFile_OffPeak_data_DchKPiPi; inFile_DchKPiPi[2] = inFile_B0B0_Dlnu_DchKPiPi; inFile_DchKPiPi[3] = inFile_BpBm_D0star_DchKPiPi; inFile_DchKPiPi[4] = inFile_B0B0_D0star_DchKPiPi; inFile_DchKPiPi[5] = inFile_BpBm_D1_DchKPiPi; inFile_DchKPiPi[6] = inFile_B0B0_D1_DchKPiPi; inFile_DchKPiPi[7] = inFile_BpBm_D1prime_DchKPiPi; inFile_DchKPiPi[8] = inFile_B0B0_D1prime_DchKPiPi; inFile_DchKPiPi[9] = inFile_BpBm_D2star_DchKPiPi; inFile_DchKPiPi[10] = inFile_B0B0_D2star_DchKPiPi; inFile_DchKPiPi[11] = inFile_BpBm_DPi_DchKPiPi; inFile_DchKPiPi[12] = inFile_B0B0_DPi_DchKPiPi; inFile_DchKPiPi[13] = inFile_BpBm_DstarPi_DchKPiPi; inFile_DchKPiPi[14] = inFile_B0B0_DstarPi_DchKPiPi; inFile_DchKPiPi[15] = inFile_BpBm_DiffB_DchKPiPi; inFile_DchKPiPi[16] = inFile_B0B0_DiffB_DchKPiPi; inFile_DchKPiPi[17] = inFile_BpBm_CascL_DchKPiPi; inFile_DchKPiPi[18] = inFile_B0B0_CascL_DchKPiPi; inFile_DchKPiPi[19] = inFile_BpBm_LMisID_DchKPiPi; inFile_DchKPiPi[20] = inFile_B0B0_LMisID_DchKPiPi; inFile_DchKPiPi[21] = inFile_BpBm_Comb_DchKPiPi; inFile_DchKPiPi[22] = inFile_B0B0_Comb_DchKPiPi; inFile_DchKPiPi[23] = inFile_ccbar_RealD_DchKPiPi; inFile_DchKPiPi[24] = inFile_ccbar_Cont_DchKPiPi; inFile_DchKPiPi[25] = inFile_uds_Cont_DchKPiPi; inFile_DchKPiPi[26] = inFile_tautau_Cont_DchKPiPi; inFile_DchKPiPi[27] = inFile_B0B0_Dstar_DchKPiPi; inFile_DchKPiPi[28] = inFile_B0B0_Dlnu_0th_DchKPiPi; inFile_DchKPiPi[29] = inFile_B0B0_Dlnu_1st_DchKPiPi; inFile_DchKPiPi[30] = inFile_B0B0_Dlnu_2nd_DchKPiPi; // Coefficients double c_D0KPi[nFiles_D0KPi]; c_D0KPi[0] = 1; c_D0KPi[1] = 1; c_D0KPi[2] = f_p0 * B_B1 / B_B1_SP6; c_D0KPi[3] = f_p0 * B_B31 / B_B31_SP6; c_D0KPi[4] = t_0p * B_B31 / B_B31_SP6; c_D0KPi[5] = f_p0 * B_B32 / B_B32_SP6; c_D0KPi[6] = t_0p * B_B32 / B_B32_SP6; c_D0KPi[7] = f_p0 * B_B33 / B_B33_SP6; c_D0KPi[8] = t_0p * B_B33 / B_B33_SP6; c_D0KPi[9] = f_p0 * B_B34 / B_B34_SP6; c_D0KPi[10] = t_0p * B_B34 / B_B34_SP6; c_D0KPi[11] = f_p0 * B_B41 / B_B41_SP6; c_D0KPi[12] = t_0p * B_B41 / B_B41_SP6; c_D0KPi[13] = f_p0 * B_B42 / B_B42_SP6; c_D0KPi[14] = t_0p * B_B42 / B_B42_SP6; c_D0KPi[15] = f_p0 * DiffB_D0; c_D0KPi[16] = DiffB_D0; c_D0KPi[17] = f_p0 * CascL_D0; c_D0KPi[18] = CascL_D0; c_D0KPi[19] = f_p0 * LMisID_D0; c_D0KPi[20] = LMisID_D0; c_D0KPi[21] = f_p0; c_D0KPi[22] = 1; c_D0KPi[23] = RealD_D0; c_D0KPi[24] = 1; c_D0KPi[25] = 1; c_D0KPi[26] = 1; c_D0KPi[27] = t_0p * B_B2 / B_B2_SP6; c_D0KPi[28] = f_p0 * B_B2 / B_B2_SP6; c_D0KPi[29] = c_D0KPi[2]; c_D0KPi[30] = c_D0KPi[2]; c_D0KPi[31] = c_D0KPi[2]; double c_DchKPiPi[nFiles_DchKPiPi]; c_DchKPiPi[0] = 1; c_DchKPiPi[1] = 1; c_DchKPiPi[2] = c_Dp0 * t_0p * B_B1 / B_B1_SP6; c_DchKPiPi[3] = c_Dp0 * f_p0 * B_B31 / B_B31_SP6; c_DchKPiPi[4] = c_Dp0 * t_0p * B_B31 / B_B31_SP6; c_DchKPiPi[5] = c_Dp0 * f_p0 * B_B32 / B_B32_SP6; c_DchKPiPi[6] = c_Dp0 * t_0p * B_B32 / B_B32_SP6; c_DchKPiPi[7] = c_Dp0 * f_p0 * B_B33 / B_B33_SP6; c_DchKPiPi[8] = c_Dp0 * t_0p * B_B33 / B_B33_SP6; c_DchKPiPi[9] = c_Dp0 * f_p0 * B_B34 / B_B34_SP6; c_DchKPiPi[10] = c_Dp0 * t_0p * B_B34 / B_B34_SP6; c_DchKPiPi[11] = c_Dp0 * f_p0 * B_B41 / B_B41_SP6; c_DchKPiPi[12] = c_Dp0 * t_0p * B_B41 / B_B41_SP6; c_DchKPiPi[13] = c_Dp0 * f_p0 * B_B42 / B_B42_SP6; c_DchKPiPi[14] = c_Dp0 * t_0p * B_B42 / B_B42_SP6; c_DchKPiPi[15] = c_Dp0 * f_p0 * DiffB_Dch; c_DchKPiPi[16] = c_Dp0 * DiffB_Dch; c_DchKPiPi[17] = c_Dp0 * f_p0 * CascL_Dch; c_DchKPiPi[18] = c_Dp0 * CascL_Dch; c_DchKPiPi[19] = c_Dp0 * f_p0 * LMisID_Dch; c_DchKPiPi[20] = c_Dp0 * LMisID_Dch; c_DchKPiPi[21] = c_Dp0 * f_p0; c_DchKPiPi[22] = 1; c_DchKPiPi[23] = c_Dp0 * RealD_Dch; c_DchKPiPi[24] = 1; c_DchKPiPi[25] = 1; c_DchKPiPi[26] = 1; c_DchKPiPi[27] = c_Dp0 * t_0p * B_B2 / B_B2_SP6; c_DchKPiPi[28] = c_DchKPiPi[2]; c_DchKPiPi[29] = c_DchKPiPi[2]; c_DchKPiPi[30] = c_DchKPiPi[2]; } #endif