7 #include "CAFAna/Cuts/QuantileCuts.h" 37 for (
int i = 0;
i < ndbins; ++
i) {
38 for (
int j = 0;
j < ndbins; ++
j) {
43 for (
int i = 0;
i < ndbins; ++
i) {
44 for (
int j = 0;
j < fdbins; ++
j) {
50 for (
int i = 0;
i < fdbins; ++
i) {
51 for (
int j = 0;
j < fdbins; ++
j) {
63 std::vector<double>
ret(ndbins);
64 for (
int i = 0;
i < ndbins; ++
i) {
81 for (
int i = 0;
i < ndbins; ++
i) {
82 for (
int j = 0;
j < ndbins; ++
j) {
87 for (
int i = 0;
i < ndbins; ++
i) {
88 for (
int j = 0;
j < fdbins; ++
j) {
94 for (
int i = 0;
i < fdbins; ++
i) {
95 for (
int j = 0;
j < fdbins; ++
j) {
106 std::vector<double>
ret(fdbins);
107 for (
int i = 0;
i < fdbins; ++
i) {
117 int ndbins,
int fdbins){
119 std::vector<double> decorrnd =
DecorrelateND(inmat, ndbins, fdbins);
120 std::vector<double> decorrfd =
DecorrelateFD(inmat, ndbins, fdbins);
122 std::vector<double> decorrerr;
123 for (
size_t i = 0;
i < decorrnd.size(); ++
i){
124 decorrerr.push_back(
std::sqrt(decorrnd.at(
i)));
126 for (
size_t i = 0;
i < decorrfd.size(); ++
i){
127 decorrerr.push_back(
std::sqrt(decorrfd.at(
i)));
std::vector< double > DecorrelateND(TMatrixD *inmat, int ndbins, int fdbins)
std::vector< double > DecorrelateFD(TMatrixD *inmat, int ndbins, int fdbins)
std::vector< double > GetDecorrelatedUncertainty(TMatrixD *inmat, int ndbins, int fdbins)