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

Only one place to insert the new module number now. Then it knows the scan parameters everywhere.

Browse Source
This commit is contained in:
redford_s
2018-04-27 16:18:41 +02:00
parent c4480848c4
commit d31a5afdb2
1 changed files with 76 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

97
CS_analysis.cpp
View File

@ -37,6 +37,61 @@ int main(int argc, char* argv[]) {
string data_loc = argv[2];
string data_file = argv[3];
vector<string> scan_v1;
scan_v1.push_back("002");
scan_v1.push_back("003");
vector<string> scan_v2;
scan_v2.push_back("032");
vector<string> scan_v3;
scan_v3.push_back("008");
vector<string> scan_v4;
scan_v4.push_back("040");
scan_v4.push_back("044");
vector<string> scan_v5;
scan_v5.push_back("006");
scan_v5.push_back("022");
scan_v5.push_back("055");
scan_v5.push_back("059");
scan_v5.push_back("065");
scan_v5.push_back("066");
scan_v5.push_back("068");
scan_v5.push_back("069");
scan_v5.push_back("079");
scan_v5.push_back("088");
scan_v5.push_back("094");
vector<string> scan_v6;
scan_v6.push_back("021");
scan_v6.push_back("029");
scan_v6.push_back("035");
scan_v6.push_back("039");
scan_v6.push_back("045");
scan_v6.push_back("048");
scan_v6.push_back("050");
scan_v6.push_back("057");
scan_v6.push_back("062");
scan_v6.push_back("064");
scan_v6.push_back("067");
scan_v6.push_back("080");
scan_v6.push_back("086");
scan_v6.push_back("100");
scan_v6.push_back("102");
scan_v6.push_back("103");
scan_v6.push_back("104");
scan_v6.push_back("108");
scan_v6.push_back("109");
scan_v6.push_back("110");
scan_v6.push_back("113");
scan_v6.push_back("119");
scan_v6.push_back("121");
scan_v6.push_back("124");
scan_v6.push_back("126");
scan_v6.push_back("226");
jungfrauFile *thisfile = new jungfrauFile();
jungfrauPixelMask *pixelMaskObject = new jungfrauPixelMask();
@ -293,15 +348,15 @@ int main(int argc, char* argv[]) {
int j = 0;
int nfiles = 0;
if (module_str == "032") {
if (find(scan_v2.begin(), scan_v2.end(), module_str) != scan_v2.end()) {
nfiles = 2;
} else if (module_str == "008" || module_str == "040" || module_str == "044") {
} else if (find(scan_v3.begin(), scan_v3.end(), module_str) != scan_v3.end() || find(scan_v4.begin(), scan_v4.end(), module_str) != scan_v4.end()) {
nfiles = 3;
} else if (module_str == "006" || module_str == "022" || module_str == "065" || module_str == "069" || module_str == "079" || module_str == "068" || module_str == "088" || module_str == "055" || module_str == "059" || module_str == "094" || module_str == "066") {
} else if (find(scan_v5.begin(), scan_v5.end(), module_str) != scan_v5.end()) {
nfiles = 4;
} else if (module_str == "021" || module_str == "029" || module_str == "035" || module_str == "039" || module_str == "045" || module_str == "048" || module_str == "050" || module_str == "057" || module_str == "062" || module_str == "064" || module_str == "067" || module_str == "080" || module_str == "086" || module_str == "100" || module_str == "102" || module_str == "103" || module_str == "104" || module_str == "108" || module_str == "109" || module_str == "110" || module_str == "113" || module_str == "119" || module_str == "121" || module_str == "124" || module_str == "126" || module_str == "226") {
} else if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
nfiles = 7;
} else if (module_str == "002" || module_str == "003") {
} else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
nfiles = 6;
}
@ -394,13 +449,13 @@ int main(int argc, char* argv[]) {
// catch the last frame (bad ctrl-c exit)
// new receiver doesn't have this problem
int last_frame = 0;
if (module_str == "032") {
if (find(scan_v2.begin(), scan_v2.end(), module_str) != scan_v2.end()) {
last_frame = 19199-1920;
} else if (module_str == "008") {
} else if (find(scan_v3.begin(), scan_v3.end(), module_str) != scan_v3.end()) {
last_frame = 24319-1920;
} else if (module_str == "040" || module_str == "044") {
} else if (find(scan_v4.begin(), scan_v4.end(), module_str) != scan_v4.end()) {
last_frame = 21759-1920;
} else if (module_str == "002" || module_str == "003") {
} else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
last_frame = 59519-1920;
}
@ -634,7 +689,7 @@ int main(int argc, char* argv[]) {
TH1F *g1overg2_isBulk = new TH1F("g1overg2_isBulk","",100,11,15);
double filter[97];
if (module_str == "040" || module_str == "044") {
if (find(scan_v4.begin(), scan_v4.end(), module_str) != scan_v4.end()) {
for (int i = 0; i < 5; i++) {
filter[i] = 0.05+(i*0.05);
}
@ -644,7 +699,7 @@ int main(int argc, char* argv[]) {
for (int i = 0; i < 17; i++) {
filter[i+5+9] = 5.+(i*2.5);
}
} else if (module_str == "006" || module_str == "022" || module_str == "065" || module_str == "069" || module_str == "079" || module_str == "068" || module_str == "088" || module_str == "055" || module_str == "059" || module_str == "094" || module_str == "066") {
} else if (find(scan_v5.begin(), scan_v5.end(), module_str) != scan_v5.end()) {
for (int i = 0; i < 9; i++) {
filter[i] = 0.05+(i*0.05);
}
@ -657,7 +712,7 @@ int main(int argc, char* argv[]) {
for (int i = 0; i < 8; i++) {
filter[i+9+17+17] = 50.+(i*10.);
}
} else if (module_str == "021" || module_str == "029" || module_str == "035" || module_str == "039" || module_str == "045" || module_str == "048" || module_str == "050" || module_str == "057" || module_str == "062" || module_str == "064" || module_str == "067" || module_str == "080" || module_str == "086" || module_str == "100" || module_str == "102" || module_str == "103" || module_str == "104" || module_str == "108" || module_str == "109" || module_str == "110" || module_str == "113" || module_str == "119" || module_str == "121" || module_str == "124" || module_str == "126" || module_str == "226") {
} else if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
// loop 1
for (int i = 0; i < 9; i++) {
filter[i] = 0.05+(i*0.05);
@ -694,7 +749,7 @@ int main(int argc, char* argv[]) {
for (int i = 0; i < 15; i++) {
filter[i+9+11+10+10+11+10+10+11] = 110.+(i*10.);
}
} else if (module_str == "002" || module_str == "003") {
} else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
for (int i = 0; i < 10; i++) {
filter[i] = 0.2+(i*0.2);
}
@ -707,7 +762,7 @@ int main(int argc, char* argv[]) {
for (int i = 0; i < 47; i++) {
filter[i+10+8+25] = 80+(i*20);
}
} else if (module_str == "008") {
} else if (find(scan_v3.begin(), scan_v3.end(), module_str) != scan_v3.end()) {
for (int i = 0; i < 9; i++) {
filter[i] = 0.05+(i*0.025);
}
@ -717,7 +772,7 @@ int main(int argc, char* argv[]) {
for (int i = 0; i < 17; i++) {
filter[i+9+9] = 5.+(i*2.5);
}
} else if (module_str == "032") {
} else if (find(scan_v2.begin(), scan_v2.end(), module_str) != scan_v2.end()) {
for (int i = 0; i < 9; i++) {
filter[i] = 0.05+(i*0.025);
}
@ -754,17 +809,17 @@ int main(int argc, char* argv[]) {
// filter out points at zero and points with 0 error (mean taken from one entry in hist)
int m_max = 0;
if (module_str == "032") {
if (find(scan_v2.begin(), scan_v2.end(), module_str) != scan_v2.end()) {
m_max = 27;
} else if (module_str == "040" || module_str == "044") {
} else if (find(scan_v4.begin(), scan_v4.end(), module_str) != scan_v4.end()) {
m_max = 31;
} else if (module_str == "008") {
} else if (find(scan_v3.begin(), scan_v3.end(), module_str) != scan_v3.end()) {
m_max = 35;
} else if (module_str == "006" || module_str == "022" || module_str == "065" || module_str == "069" || module_str == "079" || module_str == "068" || module_str == "088" || module_str == "055" || module_str == "059" || module_str == "094" || module_str == "066") {
} else if (find(scan_v5.begin(), scan_v5.end(), module_str) != scan_v5.end()) {
m_max = 51;
} else if (module_str == "021" || module_str == "029" || module_str == "035" || module_str == "039" || module_str == "045" || module_str == "048" || module_str == "050" || module_str == "057" || module_str == "062" || module_str == "064" || module_str == "067" || module_str == "080" || module_str == "086" || module_str == "100" || module_str == "102" || module_str == "103" || module_str == "104" || module_str == "108" || module_str == "109" || module_str == "110" || module_str == "113" || module_str == "119" || module_str == "121" || module_str == "124" || module_str == "126" || module_str == "226") {
} else if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
m_max = 97;
} else if (module_str == "002" || module_str == "003") {
} else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
m_max = 90;
}