Skip to content

Commit 347dad9

Browse files
Merge pull request #208 from mach3-software/kskwarczynski/feat/BinningCI
feat: Expand Binning CI
2 parents 904a626 + 8d1c006 commit 347dad9

2 files changed

Lines changed: 743 additions & 55 deletions

File tree

CIValidations/BinningValidations.cpp

Lines changed: 30 additions & 28 deletions
Original file line numberDiff line numberDiff line change
@@ -26,44 +26,46 @@ XVarBins: [0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1]
2626

2727
std::vector<int> Samples = {0, 1};
2828
std::vector<double> XVars = {-1, 0, 0.5, 0.7, 1, 10};
29+
std::vector<double> YVars = {0.05, 0.15, 0.25, 0.35};
2930
std::vector<int> NomXBins = {0, 1, 2, 3};
3031
std::vector<int> NomYBins = {0, 1, 2, 3};
3132
// Loop over all combinations
3233
for (int sample : Samples) {
3334
for (double xvar : XVars) {
34-
for (int nomXBin : NomXBins) {
35-
for (int nomYBin : NomYBins) {
36-
std::vector<const double*> KinVar;
37-
std::vector<int> NomBin;
38-
39-
if(sample == 0){
40-
double yvar = 0;
41-
switch(nomYBin){
42-
case 0: yvar = 0.05; break;
43-
case 1: yvar = 0.15; break;
44-
case 2: yvar = 0.25; break;
45-
case 3: yvar = 0.35; break;
46-
default: yvar = 0; break;
35+
for (double yvar : YVars) {
36+
for (int nomXBin : NomXBins) {
37+
for (int nomYBin : NomYBins) {
38+
std::vector<const double*> KinVar;
39+
std::vector<int> NomBin;
40+
if(sample == 0) {
41+
KinVar = {&xvar, &yvar};
42+
NomBin = {nomXBin, nomYBin};
43+
} else if (sample == 1){
44+
KinVar = {&xvar};
45+
NomBin = {nomXBin};
46+
}
47+
48+
const int GlobalBin = Binning->FindGlobalBin(sample, KinVar, NomBin);
49+
outFile << "Sample " << sample
50+
<< ", XVar: " << xvar
51+
<< ", NomXBin: " << nomXBin
52+
<< ", NomYBin: " << nomYBin
53+
<< ", GlobalBin: " << GlobalBin
54+
<< std::endl;
4755
}
48-
49-
KinVar = {&xvar, &yvar};
50-
NomBin = {nomXBin, nomYBin};
51-
} else if (sample == 1){
52-
KinVar = {&xvar};
53-
NomBin = {nomXBin};
54-
}
55-
56-
const int GlobalBin = Binning->FindGlobalBin(sample, KinVar, NomBin);
57-
outFile << "Sample " << sample
58-
<< ", XVar: " << xvar
59-
<< ", NomXBin: " << nomXBin
60-
<< ", NomYBin: " << nomYBin
61-
<< ", GlobalBin: " << GlobalBin
62-
<< std::endl;
6356
}
6457
}
6558
}
6659
}
60+
61+
for(int iSam = 0; iSam < Samples.size(); iSam++) {
62+
for (int iBin = 0; iBin < Binning->GetNBins(iSam); iBin++ ) {
63+
outFile << "Sample " << iSam
64+
<< ", bin: " << iBin
65+
<< ", name: " << Binning->GetBinName(iSam, iBin)
66+
<< std::endl;
67+
}
68+
}
6769
}
6870

6971
int main(int argc, char *argv[])

0 commit comments

Comments
 (0)