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

Added calibration for module 044. Some small changes in data paths.

Browse Source
This commit is contained in:
redford_s
2017-02-22 16:17:33 +01:00
parent 6ca856c5bb
commit ae9ebeaad0
5 changed files with 53 additions and 37 deletions
Download Patch File Download Diff File Expand all files Collapse all files
JFMC_CurrentSourceScan.cpp
+30 -16
View File
@@ -84,14 +84,26 @@ int main(int argc, char* argv[]) {
}
double filter[35];
for (int i = 0; i < 9; i++) {
filter[i] = 0.05+(i*0.025);
}
for (int i = 0; i < 9; i++) {
filter[i+9] = 0.5+(i*0.25);
}
for (int i = 0; i < 17; i++) {
filter[i+9+9] = 5.+(i*2.5);
if (module_str == "044") {
for (int i = 0; i < 5; i++) {
filter[i] = 0.05+(i*0.05);
}
for (int i = 0; i < 9; i++) {
filter[i+5] = 0.5+(i*0.25);
}
for (int i = 0; i < 17; i++) {
filter[i+5+9] = 5.+(i*2.5);
}
} else {
for (int i = 0; i < 9; i++) {
filter[i] = 0.05+(i*0.025);
}
for (int i = 0; i < 9; i++) {
filter[i+9] = 0.5+(i*0.25);
}
for (int i = 0; i < 17; i++) {
filter[i+9+9] = 5.+(i*2.5);
}
}
static double adcmeans_g1[NCH][35];
@@ -151,7 +163,7 @@ int main(int argc, char* argv[]) {
int f_max = 0;
if (module_str == "032" || module_str == "021" || module_str == "022") {
f_max = 3;
} else if (module_str == "006" || module_str == "008") {
} else if (module_str == "006" || module_str == "008" || module_str == "044") {
f_max = 4;
}
@@ -234,6 +246,8 @@ int main(int argc, char* argv[]) {
last_frame = 19199;
} else if (module_str == "006" || module_str == "008") {
last_frame = 24319;
} else if (module_str == "044") {
last_frame = 24319;
}
if (j%640 == 0 || j == last_frame) {
@@ -339,13 +353,13 @@ int main(int argc, char* argv[]) {
TH1F *g1_avg;
if (module_str == "032") {
g1_avg = new TH1F("g1_avg","",100,-12000,-4000);
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008") {
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") {
g1_avg = new TH1F("g1_avg","",100,-3000,0);
}
TH1F *g2_avg;
if (module_str == "032") {
g2_avg = new TH1F("g2_avg","",100,-1000,-200);
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008") {
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") {
g2_avg = new TH1F("g2_avg","",100,-300,0);
}
@@ -389,7 +403,7 @@ int main(int argc, char* argv[]) {
// filter out points at zero
int m_max = 0;
if (module_str == "032" || module_str == "021" || module_str == "022") {
if (module_str == "032" || module_str == "021" || module_str == "022" || module_str == "044") {
m_max = 27;
} else if (module_str == "006" || module_str == "008") {
m_max = 35;
@@ -428,7 +442,7 @@ int main(int argc, char* argv[]) {
if (module_str == "032") {
rangemin2 = 0.75;
rangemax2 = *max_element(r1_filter.begin(),r1_filter.end());
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008") {
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") {
rangemin2 = 0.075;
rangemax2 = 2.5;
}
@@ -461,7 +475,7 @@ int main(int argc, char* argv[]) {
if (module_str == "032" || module_str == "022" || module_str == "021") {
rangemin3 = 5;
rangemax3 = *max_element(r3_filter.begin(),r3_filter.end());
} else if (module_str == "006" || module_str == "008") {
} else if (module_str == "006" || module_str == "008" || module_str == "044") {
rangemin3 = 7.5;
rangemax3 = *max_element(r3_filter.begin(),r3_filter.end());
}
@@ -642,7 +656,7 @@ int main(int argc, char* argv[]) {
g1map_avg->GetYaxis()->SetTitleOffset(0.7);
if (module_str == "032") {
g1map_avg->GetZaxis()->SetRangeUser(-12000,-4000);
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008") {
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") {
g1map_avg->GetZaxis()->SetRangeUser(-3000,0);
}
g1map_avg->Draw("colz");
@@ -654,7 +668,7 @@ int main(int argc, char* argv[]) {
g2map_avg->GetYaxis()->SetTitleOffset(0.7);
if (module_str == "032") {
g2map_avg->GetZaxis()->SetRangeUser(-800,-400);
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008") {
} else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") {
g2map_avg->GetZaxis()->SetRangeUser(-300,0);
}
g2map_avg->Draw("colz");