detectors/JFCalibration
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Consistent ratios.

Browse Source
This commit is contained in:
redford_s
2019-02-08 17:20:26 +01:00
parent 8e9dadc17a
commit 0de9d9366e
3 changed files with 51 additions and 40 deletions
Download Patch File Download Diff File Expand all files Collapse all files

46
CS_analysis.cpp
View File

@ -1394,28 +1394,31 @@ int main(int argc, char* argv[]) {
// get ratio measurements
if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) {
g0overg1hist->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1erhist->Fill(abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2)));
g0overg1map->Fill(i%NC,i/NC,fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1ermap->Fill(i%NC,i/NC,abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2)));
double this_g0overg1 = fit_g0->GetParameter(1)/fit_g1->GetParameter(1);
double this_g0overg1er = abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2));
g0overg1map->Fill(i%NC,i/NC,this_g0overg1);
g0overg1hist->Fill(this_g0overg1);
g0overg1ermap->Fill(i%NC,i/NC,this_g0overg1er);
g0overg1erhist->Fill(this_g0overg1er);
if (isEdge(i)) {
g0overg1_isEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1_isEdge->Fill(this_g0overg1);
}
if (isInnerEdge(i)) {
g0overg1_isInnerEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1_isInnerEdge->Fill(this_g0overg1);
}
if (isDouble(i)) {
g0overg1_isDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1_isDouble->Fill(this_g0overg1);
}
if (isNextToDouble(i)) {
g0overg1_isNextToDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1_isNextToDouble->Fill(this_g0overg1);
}
if (isQuad(i)) {
g0overg1_isQuad->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1_isQuad->Fill(this_g0overg1);
}
if (isBulk(i)) {
g0overg1_isBulk->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1_isBulk->Fill(this_g0overg1);
}
} else {
no_g0g1_calib_map->Fill(i%NC,i/NC,1);
@ -1423,28 +1426,31 @@ int main(int argc, char* argv[]) {
if (rangemax_g1 > rangemin_g1 && rangemax_g2 > rangemin_g2) {
g1overg2hist->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2erhist->Fill(abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2)));
g1overg2map->Fill(i%NC,i/NC,fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2ermap->Fill(i%NC,i/NC,abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2)));
double this_g1overg2 = fit_g1->GetParameter(1)/fit_g2->GetParameter(1);
double this_g1overg2er = abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2));
g1overg2map->Fill(i%NC,i/NC,this_g1overg2);
g1overg2hist->Fill(this_g1overg2);
g1overg2ermap->Fill(i%NC,i/NC,this_g1overg2er);
g1overg2erhist->Fill(this_g1overg2er);
if (isEdge(i)) {
g1overg2_isEdge->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2_isEdge->Fill(this_g1overg2);
}
if (isInnerEdge(i)) {
g1overg2_isInnerEdge->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2_isInnerEdge->Fill(this_g1overg2);
}
if (isDouble(i)) {
g1overg2_isDouble->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2_isDouble->Fill(this_g1overg2);
}
if (isNextToDouble(i)) {
g1overg2_isNextToDouble->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2_isNextToDouble->Fill(this_g1overg2);
}
if (isQuad(i)) {
g1overg2_isQuad->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2_isQuad->Fill(this_g1overg2);
}
if (isBulk(i)) {
g1overg2_isBulk->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
g1overg2_isBulk->Fill(this_g1overg2);
}
} else {
no_g1g2_calib_map->Fill(i%NC,i/NC,1);