diff --git a/CIValidations/SamplePDFValidations.cpp b/CIValidations/SamplePDFValidations.cpp index ec7d322..abe72af 100755 --- a/CIValidations/SamplePDFValidations.cpp +++ b/CIValidations/SamplePDFValidations.cpp @@ -78,10 +78,16 @@ void SampleLLHValidation(std::ostream& outFile, const std::string& OriginalSampl << std::endl; } - Sample->PrintIntegral(iSample); + TString Name = "Rates_Debug_" + Sample->GetSampleTitle(iSample)+".txt"; + Sample->PrintIntegral(iSample, Name); + std::ifstream RatesFile(Name.Data()); + outFile << "Contents of " << Name << ":" << std::endl; + std::string line; + while (std::getline(RatesFile, line)) { + outFile << line << std::endl; + } + std::remove(Name.Data()); } - - std::remove(tempConfigPath.c_str()); } diff --git a/CIValidations/TestOutputs/SampleOut.txt b/CIValidations/TestOutputs/SampleOut.txt index 589346b..0abbed6 100644 --- a/CIValidations/TestOutputs/SampleOut.txt +++ b/CIValidations/TestOutputs/SampleOut.txt @@ -210057,11 +210057,53 @@ Sample: Tutorial update w2: true LLH is: 0.828340 using: Barlow-Beeston test sta Sample: Tutorial update w2: true LLH is: 85.485128 using: IceCube test stat Sample: Tutorial update w2: true LLH is: 1.773274 using: Pearson test stat Sample: Tutorial update w2: true LLH is: 0.828206 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_Tutorial.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: Tutorial} +\label{Tutorial-EventRate} +\begin{tabular}{|l|c|c|c|c|c|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} &#nu_{e} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{#mu} &#bar{#nu}_{e} #rightarrow #bar{#nu}_{e} &#nu_{#mu} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e} & Total:\\ \hline + CCQE& 5941.22340 & 15563.94656 & 10003.75090 & 13724.01899 & 71.76749 & 38.30346 & 45343.01081 \\ \hline + CC2p2h& 560.20464 & 1415.53714 & 1529.86208 & 2194.67300 & 1.69351 & 3.74035 & 5705.71070 \\ \hline + CC1pipm& 6671.96406 & 19565.76533 & 11114.37085 & 16700.44934 & 22.62902 & 16.89931 & 54092.07791 \\ \hline + CCMPi& 1938.69011 & 5453.50187 & 2583.18829 & 4716.75100 & 2.08748 & 3.44328 & 14697.66203 \\ \hline + CCDIS& 0.00000 & 0.00000 & 8.86516 & 20.62613 & 0.00000 & 0.00404 & 29.49533 \\ \hline + NCQE& 7279.00000 & 4807.00000 & 7401.00000 & 4779.00000 & 0.00000 & 0.00000 & 24266.00000 \\ \hline + NC1pipm& 1308.00000 & 1583.00000 & 1561.00000 & 1618.00000 & 0.00000 & 0.00000 & 6070.00000 \\ \hline + NC1pi0& 1906.00000 & 2075.00000 & 2201.00000 & 2005.00000 & 0.00000 & 0.00000 & 8187.00000 \\ \hline + Total:& 25605.08220 & 50463.75090 & 36403.03728 & 45758.51846 & 98.17750 & 62.39043 & 158390.95678 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: Tutorial update w2: false LLH is: 1.772403 using: Poisson test stat Sample: Tutorial update w2: false LLH is: 0.829403 using: Barlow-Beeston test stat Sample: Tutorial update w2: false LLH is: 85.469129 using: IceCube test stat Sample: Tutorial update w2: false LLH is: 1.773274 using: Pearson test stat Sample: Tutorial update w2: false LLH is: 0.829263 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_Tutorial.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: Tutorial} +\label{Tutorial-EventRate} +\begin{tabular}{|l|c|c|c|c|c|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} &#nu_{e} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{#mu} &#bar{#nu}_{e} #rightarrow #bar{#nu}_{e} &#nu_{#mu} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e} & Total:\\ \hline + CCQE& 5941.22340 & 15563.94656 & 10003.75090 & 13724.01899 & 71.76749 & 38.30346 & 45343.01081 \\ \hline + CC2p2h& 560.20464 & 1415.53714 & 1529.86208 & 2194.67300 & 1.69351 & 3.74035 & 5705.71070 \\ \hline + CC1pipm& 6671.96406 & 19565.76533 & 11114.37085 & 16700.44934 & 22.62902 & 16.89931 & 54092.07791 \\ \hline + CCMPi& 1938.69011 & 5453.50187 & 2583.18829 & 4716.75100 & 2.08748 & 3.44328 & 14697.66203 \\ \hline + CCDIS& 0.00000 & 0.00000 & 8.86516 & 20.62613 & 0.00000 & 0.00404 & 29.49533 \\ \hline + NCQE& 7279.00000 & 4807.00000 & 7401.00000 & 4779.00000 & 0.00000 & 0.00000 & 24266.00000 \\ \hline + NC1pipm& 1308.00000 & 1583.00000 & 1561.00000 & 1618.00000 & 0.00000 & 0.00000 & 6070.00000 \\ \hline + NC1pi0& 1906.00000 & 2075.00000 & 2201.00000 & 2005.00000 & 0.00000 & 0.00000 & 8187.00000 \\ \hline + Total:& 25605.08220 & 50463.75090 & 36403.03728 & 45758.51846 & 98.17750 & 62.39043 & 158390.95678 \\ \hline +\end{tabular} +\end{center} +\end{table} + Info for sample: TutorialConfigs/Samples/SampleHandler_Tutorial.yaml Rates Prior: 165774.581297 Likelihood: 1.772403 @@ -410078,11 +410120,53 @@ Sample: Tutorial ATM update w2: true LLH is: 0.080879 using: Barlow-Beeston test Sample: Tutorial ATM update w2: true LLH is: 74.416429 using: IceCube test stat Sample: Tutorial ATM update w2: true LLH is: 0.203848 using: Pearson test stat Sample: Tutorial ATM update w2: true LLH is: 0.080988 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_Tutorial ATM.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: Tutorial ATM} +\label{Tutorial ATM-EventRate} +\begin{tabular}{|l|c|c|c|c|c|c|c|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} &#nu_{e} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{#mu} &#bar{#nu}_{e} #rightarrow #bar{#nu}_{e} &#nu_{#mu} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e} &#nu_{e} #rightarrow #nu_{#mu} &#bar{#nu}_{e} #rightarrow #bar{#nu}_{#mu} & Total:\\ \hline + CCQE& 3618.85920 & 2995.28102 & 1483.38237 & 941.18400 & 13.40642 & 3.23346 & 11.73559 & 1.63455 & 9068.71660 \\ \hline + CC2p2h& 549.56867 & 480.64439 & 264.62893 & 179.16742 & 2.14132 & 0.36551 & 1.41395 & 0.10672 & 1478.03690 \\ \hline + CC1pipm& 9819.45550 & 7365.09731 & 2241.62249 & 1523.00494 & 44.31114 & 3.73757 & 19.81892 & 2.02600 & 21019.07388 \\ \hline + CCMPi& 5148.29219 & 3465.74575 & 823.93656 & 526.10507 & 20.62713 & 1.36887 & 16.54939 & 0.51682 & 10003.14178 \\ \hline + CCDIS& 0.00000 & 0.00000 & 3.34241 & 2.00633 & 0.00000 & 0.00725 & 0.00000 & 0.01025 & 5.36624 \\ \hline + NCQE& 1292.00000 & 869.00000 & 628.00000 & 302.00000 & 0.00000 & 0.00000 & 0.00000 & 0.00000 & 3091.00000 \\ \hline + NC1pipm& 1034.00000 & 555.00000 & 290.00000 & 143.00000 & 0.00000 & 0.00000 & 0.00000 & 0.00000 & 2022.00000 \\ \hline + NC1pi0& 1035.00000 & 654.00000 & 354.00000 & 163.00000 & 0.00000 & 0.00000 & 0.00000 & 0.00000 & 2206.00000 \\ \hline + Total:& 22497.17557 & 16384.76846 & 6088.91277 & 3779.46775 & 80.48601 & 8.71266 & 49.51786 & 4.29433 & 48893.33540 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: Tutorial ATM update w2: false LLH is: 0.204242 using: Poisson test stat Sample: Tutorial ATM update w2: false LLH is: 0.081416 using: Barlow-Beeston test stat Sample: Tutorial ATM update w2: false LLH is: 74.386708 using: IceCube test stat Sample: Tutorial ATM update w2: false LLH is: 0.203848 using: Pearson test stat Sample: Tutorial ATM update w2: false LLH is: 0.081524 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_Tutorial ATM.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: Tutorial ATM} +\label{Tutorial ATM-EventRate} +\begin{tabular}{|l|c|c|c|c|c|c|c|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} &#nu_{e} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{#mu} &#bar{#nu}_{e} #rightarrow #bar{#nu}_{e} &#nu_{#mu} #rightarrow #nu_{e} &#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e} &#nu_{e} #rightarrow #nu_{#mu} &#bar{#nu}_{e} #rightarrow #bar{#nu}_{#mu} & Total:\\ \hline + CCQE& 3618.85920 & 2995.28102 & 1483.38237 & 941.18400 & 13.40642 & 3.23346 & 11.73559 & 1.63455 & 9068.71660 \\ \hline + CC2p2h& 549.56867 & 480.64439 & 264.62893 & 179.16742 & 2.14132 & 0.36551 & 1.41395 & 0.10672 & 1478.03690 \\ \hline + CC1pipm& 9819.45550 & 7365.09731 & 2241.62249 & 1523.00494 & 44.31114 & 3.73757 & 19.81892 & 2.02600 & 21019.07388 \\ \hline + CCMPi& 5148.29219 & 3465.74575 & 823.93656 & 526.10507 & 20.62713 & 1.36887 & 16.54939 & 0.51682 & 10003.14178 \\ \hline + CCDIS& 0.00000 & 0.00000 & 3.34241 & 2.00633 & 0.00000 & 0.00725 & 0.00000 & 0.01025 & 5.36624 \\ \hline + NCQE& 1292.00000 & 869.00000 & 628.00000 & 302.00000 & 0.00000 & 0.00000 & 0.00000 & 0.00000 & 3091.00000 \\ \hline + NC1pipm& 1034.00000 & 555.00000 & 290.00000 & 143.00000 & 0.00000 & 0.00000 & 0.00000 & 0.00000 & 2022.00000 \\ \hline + NC1pi0& 1035.00000 & 654.00000 & 354.00000 & 163.00000 & 0.00000 & 0.00000 & 0.00000 & 0.00000 & 2206.00000 \\ \hline + Total:& 22497.17557 & 16384.76846 & 6088.91277 & 3779.46775 & 80.48601 & 8.71266 & 49.51786 & 4.29433 & 48893.33540 \\ \hline +\end{tabular} +\end{center} +\end{table} + Info for sample: TutorialConfigs/Samples/SampleHandler_Tutorial_ATM.yaml Rates Prior: 52470.872110 Likelihood: 0.204242 @@ -443809,51 +443893,261 @@ Sample: ND_CC0pi update w2: true LLH is: 0.000000 using: Barlow-Beeston test sta Sample: ND_CC0pi update w2: true LLH is: 308.034777 using: IceCube test stat Sample: ND_CC0pi update w2: true LLH is: 0.000000 using: Pearson test stat Sample: ND_CC0pi update w2: true LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_CC0pi.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_CC0pi} +\label{ND_CC0pi-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 17933.28372 & 17933.28372 \\ \hline + CC2p2h& 726.95591 & 726.95591 \\ \hline + CC1pipm& 1363.85053 & 1363.85053 \\ \hline + CCMPi& 101.10184 & 101.10184 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 0.00000 & 0.00000 \\ \hline + NC1pipm& 0.00000 & 0.00000 \\ \hline + NC1pi0& 0.00000 & 0.00000 \\ \hline + Total:& 20125.19200 & 20125.19200 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_CC1pip update w2: true LLH is: 0.000000 using: Poisson test stat Sample: ND_CC1pip update w2: true LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_CC1pip update w2: true LLH is: 325.823261 using: IceCube test stat Sample: ND_CC1pip update w2: true LLH is: 0.000000 using: Pearson test stat Sample: ND_CC1pip update w2: true LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_CC1pip.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_CC1pip} +\label{ND_CC1pip-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 120.64101 & 120.64101 \\ \hline + CC2p2h& 7.48513 & 7.48513 \\ \hline + CC1pipm& 6690.18955 & 6690.18955 \\ \hline + CCMPi& 104.48775 & 104.48775 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 0.00000 & 0.00000 \\ \hline + NC1pipm& 0.00000 & 0.00000 \\ \hline + NC1pi0& 0.00000 & 0.00000 \\ \hline + Total:& 6922.80343 & 6922.80343 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_CCPi0 update w2: true LLH is: 0.000000 using: Poisson test stat Sample: ND_CCPi0 update w2: true LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_CCPi0 update w2: true LLH is: 272.910885 using: IceCube test stat Sample: ND_CCPi0 update w2: true LLH is: 0.000000 using: Pearson test stat Sample: ND_CCPi0 update w2: true LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_CCPi0.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_CCPi0} +\label{ND_CCPi0-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 59.18189 & 59.18189 \\ \hline + CC2p2h& 2.37297 & 2.37297 \\ \hline + CC1pipm& 2263.05571 & 2263.05571 \\ \hline + CCMPi& 60.06771 & 60.06771 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 0.00000 & 0.00000 \\ \hline + NC1pipm& 0.00000 & 0.00000 \\ \hline + NC1pi0& 0.00000 & 0.00000 \\ \hline + Total:& 2384.67829 & 2384.67829 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_NC0pi update w2: true LLH is: 0.000000 using: Poisson test stat Sample: ND_NC0pi update w2: true LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_NC0pi update w2: true LLH is: 144.116287 using: IceCube test stat Sample: ND_NC0pi update w2: true LLH is: 0.000000 using: Pearson test stat Sample: ND_NC0pi update w2: true LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_NC0pi.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_NC0pi} +\label{ND_NC0pi-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 0.00000 & 0.00000 \\ \hline + CC2p2h& 0.00000 & 0.00000 \\ \hline + CC1pipm& 0.00000 & 0.00000 \\ \hline + CCMPi& 0.00000 & 0.00000 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 5799.00000 & 5799.00000 \\ \hline + NC1pipm& 161.00000 & 161.00000 \\ \hline + NC1pi0& 243.00000 & 243.00000 \\ \hline + Total:& 6203.00000 & 6203.00000 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_NC1pi0 update w2: true LLH is: 0.000000 using: Poisson test stat Sample: ND_NC1pi0 update w2: true LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_NC1pi0 update w2: true LLH is: 55.386077 using: IceCube test stat Sample: ND_NC1pi0 update w2: true LLH is: 0.000000 using: Pearson test stat Sample: ND_NC1pi0 update w2: true LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_NC1pi0.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_NC1pi0} +\label{ND_NC1pi0-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 0.00000 & 0.00000 \\ \hline + CC2p2h& 0.00000 & 0.00000 \\ \hline + CC1pipm& 0.00000 & 0.00000 \\ \hline + CCMPi& 0.00000 & 0.00000 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 9.00000 & 9.00000 \\ \hline + NC1pipm& 40.00000 & 40.00000 \\ \hline + NC1pi0& 1294.00000 & 1294.00000 \\ \hline + Total:& 1343.00000 & 1343.00000 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_CC0pi update w2: false LLH is: 0.000000 using: Poisson test stat Sample: ND_CC0pi update w2: false LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_CC0pi update w2: false LLH is: 308.034777 using: IceCube test stat Sample: ND_CC0pi update w2: false LLH is: 0.000000 using: Pearson test stat Sample: ND_CC0pi update w2: false LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_CC0pi.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_CC0pi} +\label{ND_CC0pi-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 17933.28372 & 17933.28372 \\ \hline + CC2p2h& 726.95591 & 726.95591 \\ \hline + CC1pipm& 1363.85053 & 1363.85053 \\ \hline + CCMPi& 101.10184 & 101.10184 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 0.00000 & 0.00000 \\ \hline + NC1pipm& 0.00000 & 0.00000 \\ \hline + NC1pi0& 0.00000 & 0.00000 \\ \hline + Total:& 20125.19200 & 20125.19200 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_CC1pip update w2: false LLH is: 0.000000 using: Poisson test stat Sample: ND_CC1pip update w2: false LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_CC1pip update w2: false LLH is: 325.823261 using: IceCube test stat Sample: ND_CC1pip update w2: false LLH is: 0.000000 using: Pearson test stat Sample: ND_CC1pip update w2: false LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_CC1pip.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_CC1pip} +\label{ND_CC1pip-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 120.64101 & 120.64101 \\ \hline + CC2p2h& 7.48513 & 7.48513 \\ \hline + CC1pipm& 6690.18955 & 6690.18955 \\ \hline + CCMPi& 104.48775 & 104.48775 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 0.00000 & 0.00000 \\ \hline + NC1pipm& 0.00000 & 0.00000 \\ \hline + NC1pi0& 0.00000 & 0.00000 \\ \hline + Total:& 6922.80343 & 6922.80343 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_CCPi0 update w2: false LLH is: 0.000000 using: Poisson test stat Sample: ND_CCPi0 update w2: false LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_CCPi0 update w2: false LLH is: 272.910885 using: IceCube test stat Sample: ND_CCPi0 update w2: false LLH is: 0.000000 using: Pearson test stat Sample: ND_CCPi0 update w2: false LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_CCPi0.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_CCPi0} +\label{ND_CCPi0-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 59.18189 & 59.18189 \\ \hline + CC2p2h& 2.37297 & 2.37297 \\ \hline + CC1pipm& 2263.05571 & 2263.05571 \\ \hline + CCMPi& 60.06771 & 60.06771 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 0.00000 & 0.00000 \\ \hline + NC1pipm& 0.00000 & 0.00000 \\ \hline + NC1pi0& 0.00000 & 0.00000 \\ \hline + Total:& 2384.67829 & 2384.67829 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_NC0pi update w2: false LLH is: 0.000000 using: Poisson test stat Sample: ND_NC0pi update w2: false LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_NC0pi update w2: false LLH is: 144.116287 using: IceCube test stat Sample: ND_NC0pi update w2: false LLH is: 0.000000 using: Pearson test stat Sample: ND_NC0pi update w2: false LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_NC0pi.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_NC0pi} +\label{ND_NC0pi-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 0.00000 & 0.00000 \\ \hline + CC2p2h& 0.00000 & 0.00000 \\ \hline + CC1pipm& 0.00000 & 0.00000 \\ \hline + CCMPi& 0.00000 & 0.00000 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 5799.00000 & 5799.00000 \\ \hline + NC1pipm& 161.00000 & 161.00000 \\ \hline + NC1pi0& 243.00000 & 243.00000 \\ \hline + Total:& 6203.00000 & 6203.00000 \\ \hline +\end{tabular} +\end{center} +\end{table} + Sample: ND_NC1pi0 update w2: false LLH is: 0.000000 using: Poisson test stat Sample: ND_NC1pi0 update w2: false LLH is: 0.000000 using: Barlow-Beeston test stat Sample: ND_NC1pi0 update w2: false LLH is: 55.386077 using: IceCube test stat Sample: ND_NC1pi0 update w2: false LLH is: 0.000000 using: Pearson test stat Sample: ND_NC1pi0 update w2: false LLH is: 0.000000 using: Dembinski-Abdelmotteleb test stat +Contents of Rates_Debug_ND_NC1pi0.txt: +\begin{table}[ht] +\begin{center} +\caption{Integral breakdown for sample: ND_NC1pi0} +\label{ND_NC1pi0-EventRate} +\begin{tabular}{|l|c|c|} +\hline + Mode:&#nu_{#mu} #rightarrow #nu_{#mu} & Total:\\ \hline + CCQE& 0.00000 & 0.00000 \\ \hline + CC2p2h& 0.00000 & 0.00000 \\ \hline + CC1pipm& 0.00000 & 0.00000 \\ \hline + CCMPi& 0.00000 & 0.00000 \\ \hline + CCDIS& 0.00000 & 0.00000 \\ \hline + NCQE& 9.00000 & 9.00000 \\ \hline + NC1pipm& 40.00000 & 40.00000 \\ \hline + NC1pi0& 1294.00000 & 1294.00000 \\ \hline + Total:& 1343.00000 & 1343.00000 \\ \hline +\end{tabular} +\end{center} +\end{table} + Info for sample: TutorialConfigs/Samples/SampleHandler_Tutorial_ND.yaml Rates Prior: 20125.192000 Likelihood: 17190.607414 diff --git a/TutorialConfigs/Samples/SampleHandler_Tutorial_ND.yaml b/TutorialConfigs/Samples/SampleHandler_Tutorial_ND.yaml index a0315bc..dfe0cbc 100755 --- a/TutorialConfigs/Samples/SampleHandler_Tutorial_ND.yaml +++ b/TutorialConfigs/Samples/SampleHandler_Tutorial_ND.yaml @@ -4,9 +4,6 @@ SampleHandlerName: "ND2137" POT: 2.1428e+21 MaCh3ModeConfig: "TutorialConfigs/MaCh3Modes.yaml" -NuOsc: - NuOscConfigFile: "TutorialConfigs/NuOscillator/NuFASTLinear.yaml" - EqualBinningPerOscChannel: true LikelihoodOptions: #False means you calculate W2 histogram only during first reweight, advisable for Barlow-Beeston