detectors/JFCalibration
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity
Files
master
JFCalibration/CS_fit_M431.cpp
hinger_v 73938aad04 new file: BP_analysis_M431.cpp 
new file:   BP_analysis_thinSensors.cpp
	new file:   BP_fit.cpp
	new file:   BP_fit_M431.cpp
	new file:   BP_fit_thin.cpp
	new file:   BP_scan_both_speeds.sh
	new file:   BP_scan_thinSensor.sh
	new file:   CS_analysis_M431.cpp
	new file:   CS_analysis_M439.cpp
	new file:   CS_fit.cpp
	new file:   CS_fit_M431
	new file:   CS_fit_M431.cpp
	new file:   CS_scan_both_speeds.sh
	new file:   CuFluo_exposure_both_speeds.sh
	new file:   CuFluo_fit.cpp
	new file:   Default_pixels_arrays.cpp
2025-06-27 11:16:26 +02:00

1188 lines
37 KiB
C++
Raw Permalink Blame History

// to analyse the backplane pulsing data per module
// changes by VH 210906: to eliminate hardcoded absolute paths, uses location of the analysis root files as additional input argument (accordingly changed in filename_creator.sh)
#include "TApplication.h"
#include "sls_detector_calibration/jungfrauCommonHeader.h"
#include "sls_detector_calibration/jungfrauCommonFunctions.h"
#include "sls_detector_calibration/jungfrauFile.C"
#include "sls_detector_calibration/jungfrauPedestal.C"
#include "TGraph.h"
#include "TGraphErrors.h"
#include "TF1.h"
#include "TFile.h"
#include "TPaveStats.h"
#include "TLegend.h"
#include "TPaveText.h"
#include "TCanvas.h"
#include "TFitResultPtr.h"
#include "TFitResult.h"
#include <algorithm>
#include <sys/stat.h>
#include <sstream>
TGraphErrors *grap_g0;
TGraphErrors *grap_g1;
TGraphErrors *grap_g2;
TF1 *fit_g0=0;
TF1 *fit_g1=0;
TF1 *fit_g2=0;
TGraphErrors *norm_g0;
TGraphErrors *norm_g1;
TGraphErrors *norm_g2;
TF1 *flat_g0;
TF1 *flat_g2;
TF1 *lin_g1_p02pc;
TF1 *lin_g1_p05pc;
TF1 *lin_g1_m02pc;
TF1 *lin_g1_m05pc;
TF1 *lin_g2_p02pc;
TF1 *lin_g2_p05pc;
TF1 *lin_g2_m02pc;
TF1 *lin_g2_m05pc;
Double_t par1;
Double_t par2;
double checkRangeMinForRCEffect(double range_min) {
// check that the range minimum is at least 2 us, to avoid non-linearity of RC effect
if (range_min < 2.0) {
return 2.0;
} else {
return range_min;
}
}
bool isHGX=false;
int main(int argc, char* argv[]) {
//nonblock(NB_ENABLE);
cout <<"opening the rootapp" <<endl;
TApplication rootapp("example",&argc, argv);
jungfrauStyle();
//gROOT->SetBatch(1);
gStyle->SetOptFit(11);
/*
if (argc != 3) {
cout << "Correct usage:" << endl;
cout << "arg 1: specify module number" << endl;
cout << "arg 2: specify data location" << endl;
cout << "arg 3: specify column (x)" << endl;
cout << "arg 4: specify row (y)" << endl;
exit(1);
}
*/ //uncomment for SR
if (argc != 4) {
cout << "Correct usage:" << endl;
cout << "arg 1: specify module number" << endl;
cout << "arg 2: specify pixel x position from 0 to 1023" << endl;
cout << "arg 3: specify pixel y position from 0 to 511" << endl;
exit(1);
} //uncomment for VH 210906
string module_str = argv[1];
string str2 =("HG0G1G2");
string str3 =("HGOG1G2"); //filename creator had this bug
int column = atoi(argv[2]);
int row = atoi(argv[3]);
int pixel = column+row*1024;
char data_loc[256];
sprintf(data_loc,"/mnt/sls_det_storage/jungfrau_calib/jungfrau_ana_sophie/M%s_CalibAna/", module_str.c_str());
cout << data_loc << endl;
std::string folder_path(data_loc);
if (folder_path.find(str2) != string::npos) isHGX=true;
if (folder_path.find(str3) != string::npos) isHGX=true;
// cout << data_loc.find(str2)<<" " << string::npos << " " << str2 << " " << data_loc <<endl;
if (isHGX) {
cout << " HG0->HG1->HG2 sequence - dynamicHG0" <<endl;
// plotfolder_str="BackplanePulsing_HG0G1G2";
}
else {
cout << " G0->G1->G2 sequence - dynamicG0" <<endl;
// plotfolder_str="BackplanePulsing";
}
// Define vector with module name for the different exposure time scan
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;
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("065");
scan_v5.push_back("066");
scan_v5.push_back("068");
scan_v5.push_back("069");
scan_v5.push_back("088");
scan_v5.push_back("094");
vector<string> scan_v6;
scan_v6.push_back("008");
scan_v6.push_back("011");
scan_v6.push_back("021");
scan_v6.push_back("023");
scan_v6.push_back("025");
scan_v6.push_back("029");
scan_v6.push_back("035");
scan_v6.push_back("039");
scan_v6.push_back("043");
scan_v6.push_back("045");
scan_v6.push_back("046");
scan_v6.push_back("048");
scan_v6.push_back("049");
scan_v6.push_back("050");
scan_v6.push_back("051");
scan_v6.push_back("053");
scan_v6.push_back("054");
scan_v6.push_back("057");
scan_v6.push_back("058");
scan_v6.push_back("059");
scan_v6.push_back("060");
scan_v6.push_back("061");
scan_v6.push_back("062");
scan_v6.push_back("063");
scan_v6.push_back("064");
scan_v6.push_back("067");
scan_v6.push_back("070");
scan_v6.push_back("071");
scan_v6.push_back("072");
scan_v6.push_back("073");
scan_v6.push_back("074");
scan_v6.push_back("075");
scan_v6.push_back("076");
scan_v6.push_back("079");
scan_v6.push_back("080");
scan_v6.push_back("086");
scan_v6.push_back("089");
scan_v6.push_back("090");
scan_v6.push_back("091");
scan_v6.push_back("092");
scan_v6.push_back("093");
scan_v6.push_back("095");
scan_v6.push_back("097");
scan_v6.push_back("098");
scan_v6.push_back("099");
scan_v6.push_back("100");
scan_v6.push_back("102");
scan_v6.push_back("103");
scan_v6.push_back("104");
scan_v6.push_back("106");
scan_v6.push_back("107");
scan_v6.push_back("108");
scan_v6.push_back("109");
scan_v6.push_back("110");
scan_v6.push_back("111");
scan_v6.push_back("112");
scan_v6.push_back("113");
scan_v6.push_back("114");
scan_v6.push_back("115");
scan_v6.push_back("116");
scan_v6.push_back("117");
scan_v6.push_back("118");
scan_v6.push_back("119");
scan_v6.push_back("120");
scan_v6.push_back("121");
scan_v6.push_back("122");
scan_v6.push_back("123");
scan_v6.push_back("124");
scan_v6.push_back("125");
scan_v6.push_back("126");
scan_v6.push_back("129");
scan_v6.push_back("130");
scan_v6.push_back("131");
scan_v6.push_back("132");
scan_v6.push_back("133");
scan_v6.push_back("134");
scan_v6.push_back("135");
scan_v6.push_back("200");
scan_v6.push_back("201");
scan_v6.push_back("202");
scan_v6.push_back("203");
scan_v6.push_back("204");
scan_v6.push_back("205");
scan_v6.push_back("206");
scan_v6.push_back("208");
scan_v6.push_back("209");
scan_v6.push_back("210");
scan_v6.push_back("212");
scan_v6.push_back("214");
scan_v6.push_back("215");
scan_v6.push_back("216");
scan_v6.push_back("217");
scan_v6.push_back("218");
scan_v6.push_back("219");
scan_v6.push_back("221");
scan_v6.push_back("222");
scan_v6.push_back("223");
scan_v6.push_back("224");
scan_v6.push_back("225");
scan_v6.push_back("226");
scan_v6.push_back("227");
scan_v6.push_back("228");
scan_v6.push_back("229");
scan_v6.push_back("232");
scan_v6.push_back("233");
scan_v6.push_back("234");
scan_v6.push_back("235");
scan_v6.push_back("236");
scan_v6.push_back("242");
scan_v6.push_back("244");
scan_v6.push_back("245");
scan_v6.push_back("246");
scan_v6.push_back("248");
scan_v6.push_back("249");
scan_v6.push_back("250");
scan_v6.push_back("251");
scan_v6.push_back("252");
scan_v6.push_back("253");
scan_v6.push_back("254");
scan_v6.push_back("255");
scan_v6.push_back("260");
scan_v6.push_back("261");
scan_v6.push_back("263");
scan_v6.push_back("264");
scan_v6.push_back("265");
scan_v6.push_back("266");
scan_v6.push_back("267");
scan_v6.push_back("269");
scan_v6.push_back("270");
scan_v6.push_back("272");
scan_v6.push_back("273");
scan_v6.push_back("274");
scan_v6.push_back("275");
scan_v6.push_back("276");
scan_v6.push_back("277");
scan_v6.push_back("278");
scan_v6.push_back("279");
scan_v6.push_back("281");
scan_v6.push_back("282");
scan_v6.push_back("283");
scan_v6.push_back("285");
scan_v6.push_back("286");
scan_v6.push_back("287");
scan_v6.push_back("288");
scan_v6.push_back("289");
scan_v6.push_back("290");
scan_v6.push_back("291");
scan_v6.push_back("292");
scan_v6.push_back("294");
scan_v6.push_back("296");
scan_v6.push_back("298");
scan_v6.push_back("299");
scan_v6.push_back("301");
scan_v6.push_back("307");
scan_v6.push_back("312");
scan_v6.push_back("351");
scan_v6.push_back("352");
scan_v6.push_back("354");
scan_v6.push_back("360");
scan_v6.push_back("373");
scan_v6.push_back("380");
scan_v6.push_back("304");
scan_v6.push_back("321");
scan_v6.push_back("346");
scan_v6.push_back("375");
scan_v6.push_back("381");
scan_v6.push_back("383");
scan_v6.push_back("310");
scan_v6.push_back("313");
scan_v6.push_back("314");
scan_v6.push_back("334");
scan_v6.push_back("337");
scan_v6.push_back("347");
scan_v6.push_back("361");
scan_v6.push_back("362");
scan_v6.push_back("355");
scan_v6.push_back("364");
scan_v6.push_back("366");
scan_v6.push_back("243");
scan_v6.push_back("339");
scan_v6.push_back("341");
scan_v6.push_back("368");
scan_v6.push_back("376");
scan_v6.push_back("379");
scan_v6.push_back("359");
scan_v6.push_back("363");
scan_v6.push_back("365");
scan_v6.push_back("347");
scan_v6.push_back("366");
scan_v6.push_back("268");
scan_v6.push_back("300");
scan_v6.push_back("302");
scan_v6.push_back("303");
scan_v6.push_back("309");
scan_v6.push_back("326");
scan_v6.push_back("308");
scan_v6.push_back("322");
scan_v6.push_back("328");
scan_v6.push_back("332");
scan_v6.push_back("335");
scan_v6.push_back("353");
scan_v6.push_back("335_FW2_2");
scan_v6.push_back("353_FW2_2");
scan_v6.push_back("262_FWcheck_1_1_1");
scan_v6.push_back("262_FWcheck_1_1_1_test");
scan_v6.push_back("262_FWcheck_1_2_rc210831");
scan_v6.push_back("262_FWcheck_1_2_defdacs");
scan_v6.push_back("311");
scan_v6.push_back("316");
scan_v6.push_back("318");
scan_v6.push_back("319");
scan_v6.push_back("320");
scan_v6.push_back("324");
scan_v6.push_back("272_fullspeedtest");
scan_v6.push_back("272_halfspeedtest");
scan_v6.push_back("352_fullspeed"); //jf4mb2-00 (SLS)
scan_v6.push_back("307_fullspeed"); //jf4mb2-01 (SLS)
scan_v6.push_back("264_fullspeed"); //jf4mb2-02 (SLS)
scan_v6.push_back("253_fullspeed"); //jf4mb2-03 (SLS)
scan_v6.push_back("351_fullspeed"); //jf4mb2-04 (SLS)
scan_v6.push_back("312_fullspeed"); //jf4mb2-05 (SLS)
scan_v6.push_back("261_fullspeed"); //jf4mb2-06 (SLS)
scan_v6.push_back("310_fullspeed"); //jf4mb2-07 (SLS)
scan_v6.push_back("384"); //jfcalib-5
scan_v6.push_back("510"); //jfcalib-0
scan_v6.push_back("511"); //jfcalib-1
scan_v6.push_back("512"); //jfcalib-2
scan_v6.push_back("556"); //jfcalib-3
scan_v6.push_back("572"); //jfcalib-4
scan_v6.push_back("329"); //jfcalib-0
scan_v6.push_back("331"); //jfcalib-1
scan_v6.push_back("333"); //jfcalib-2
scan_v6.push_back("336"); //jfcalib-3
scan_v6.push_back("356"); //jfcalib-4
scan_v6.push_back("377"); //jfcalib-5
scan_v6.push_back("387"); //jfcalib-0
scan_v6.push_back("388"); //jfcalib-1
scan_v6.push_back("389"); //jfcalib-2
scan_v6.push_back("396"); //jfcalib-3
scan_v6.push_back("397"); //jfcalib-4
scan_v6.push_back("398"); //jfcalib-5
scan_v6.push_back("390"); //jfcalib-0
scan_v6.push_back("391"); //jfcalib-1
scan_v6.push_back("393"); //jfcalib-2
scan_v6.push_back("399"); //jfcalib-3
scan_v6.push_back("516"); //jfcalib-4
scan_v6.push_back("517"); //jfcalib-5
scan_v6.push_back("519"); //jfcalib-0
scan_v6.push_back("520"); //jfcalib-1
scan_v6.push_back("532"); //jfcalib-2
scan_v6.push_back("539"); //jfcalib-3
scan_v6.push_back("555"); //jfcalib-4
//Calib for VELA beamtime in March 2022
scan_v6.push_back("262"); //jfcalib-0
scan_v6.push_back("231"); //jfcalib-1
scan_v6.push_back("305"); //jfcalib-2
scan_v6.push_back("392"); //jfcalib-1
scan_v6.push_back("394"); //jfcalib-2
scan_v6.push_back("513"); //jfcalib-3
scan_v6.push_back("514"); //jfcalib-4
scan_v6.push_back("515"); //jfcalib-5
scan_v6.push_back("239"); //jfcalib-0 //bad
scan_v6.push_back("064"); //jfcalib-1 //bad //already added above
scan_v6.push_back("023"); //jfcalib-2 //bad //already added above
scan_v6.push_back("349"); //jfcalib-3 //bad
scan_v6.push_back("518"); //jfcalib-4
scan_v6.push_back("521"); //jfcalib-5
scan_v6.push_back("523"); //jfcalib-0
scan_v6.push_back("526"); //jfcalib-1
scan_v6.push_back("529"); //jfcalib-2
scan_v6.push_back("530"); //jfcalib-3
scan_v6.push_back("531"); //jfcalib-4
scan_v6.push_back("533"); //jfcalib-5
scan_v6.push_back("524"); //1M for ESRF
scan_v6.push_back("525"); //1M for ESRF
scan_v6.push_back("306"); //jfcalib-0 (strixel)
scan_v6.push_back("350"); //jfcalib-1 (strixel)
scan_v6.push_back("534"); //jfcalib-0
scan_v6.push_back("537"); //jfcalib-1
scan_v6.push_back("538"); //jfcalib-2
scan_v6.push_back("542"); //jfcalib-3
scan_v6.push_back("543"); //jfcalib-4
scan_v6.push_back("544"); //jfcalib-5
scan_v6.push_back("402"); //jfcalib-0
scan_v6.push_back("404"); //jfcalib-1
scan_v6.push_back("407"); //jfcalib-2
scan_v6.push_back("409"); //jfcalib-3
scan_v6.push_back("410"); //jfcalib-4
scan_v6.push_back("411"); //jfcalib-5
//scan_v6.push_back("008"); //jfcalib-0 //already in the list from previous calibration
scan_v6.push_back("433"); //jfcalib-1
scan_v6.push_back("550"); //jfcalib-2
scan_v6.push_back("545"); //jfcalib-3
scan_v6.push_back("546"); //jfcalib-4
scan_v6.push_back("547"); //jfcalib-5
//scan_v6.push_back("416"); //jfcalib-0 //destroyed, machanical sample
scan_v6.push_back("417"); //jfcalib-1
scan_v6.push_back("418"); //jfcalib-2
scan_v6.push_back("420"); //jfcalib-3
scan_v6.push_back("421"); //jfcalib-4
scan_v6.push_back("430"); //jfcalib-5
scan_v6.push_back("417_fullspeed"); //jfcalib-1
scan_v6.push_back("418_fullspeed"); //jfcalib-2
scan_v6.push_back("420_fullspeed"); //jfcalib-3
scan_v6.push_back("421_fullspeed"); //jfcalib-4
scan_v6.push_back("430_fullspeed"); //jfcalib-5
scan_v6.push_back("434"); //jfcalib-0
scan_v6.push_back("437"); //jfcalib-1
scan_v6.push_back("480"); //jfcalib-2
scan_v6.push_back("464"); //jfcalib-3
scan_v6.push_back("473"); //jfcalib-4
scan_v6.push_back("474"); //jfcalib-5
scan_v6.push_back("076"); //jfcalib-0
scan_v6.push_back("089"); //jfcalib-3
scan_v6.push_back("207"); //jfcalib-5
scan_v6.push_back("244"); //jfcalib-3
scan_v6.push_back("474_fullspeed"); //jfcalib-0
scan_v6.push_back("437_fullspeed"); //jfcalib-1
scan_v6.push_back("480_fullspeed"); //jfcalib-2
scan_v6.push_back("482_fullspeed"); //jfcalib-3
scan_v6.push_back("487_fullspeed"); //jfcalib-5
scan_v6.push_back("482"); //jfcalib-3
scan_v6.push_back("487"); //jfcalib-5
scan_v6.push_back("499_fullspeed"); //jfcalib-0
scan_v6.push_back("548_fullspeed"); //jfcalib-1
scan_v6.push_back("549_fullspeed"); //jfcalib-2
scan_v6.push_back("427_fullspeed"); //jfcalib-3
scan_v6.push_back("499"); //jfcalib-0
scan_v6.push_back("548"); //jfcalib-1
scan_v6.push_back("549"); //jfcalib-2
scan_v6.push_back("427"); //jfcalib-3
scan_v6.push_back("435"); //jfcalib-5
scan_v6.push_back("403"); //jfcalib-0
scan_v6.push_back("406"); //jfcalib-1
scan_v6.push_back("423"); //jfcalib-2
scan_v6.push_back("436"); //jfcalib-3
scan_v6.push_back("479"); //jfcalib-5
// Calibration 450 um sensor for ESRF
scan_v6.push_back("583"); //jfcalib-0
scan_v6.push_back("593"); //jfcalib-1
scan_v6.push_back("601"); //jfcalib-2
scan_v6.push_back("604"); //jfcalib-3
scan_v6.push_back("605"); //jfcalib-5
// Calibration 320 um sensors
scan_v6.push_back("441"); //jfcalib-0
scan_v6.push_back("460"); //jfcalib-1
scan_v6.push_back("463"); //jfcalib-2
scan_v6.push_back("476"); //jfcalib-3
scan_v6.push_back("371"); //jfcalib-4
scan_v6.push_back("429"); //jfcalib-0
scan_v6.push_back("442"); //jfcalib-1
scan_v6.push_back("444"); //jfcalib-2
scan_v6.push_back("446"); //jfcalib-3
scan_v6.push_back("447"); //jfcalib-4
scan_v6.push_back("428"); //jfcalib-0
scan_v6.push_back("439"); //jfcalib-1
scan_v6.push_back("443"); //jfcalib-2
scan_v6.push_back("445"); //jfcalib-3
scan_v6.push_back("448"); //jfcalib-0
scan_v6.push_back("454"); //jfcalib-1
scan_v6.push_back("455"); //jfcalib-2
scan_v6.push_back("456"); //jfcalib-3
scan_v6.push_back("466"); //jfcalib-4
scan_v6.push_back("471"); //jfcalib-5
scan_v6.push_back("481"); //jfcalib-0
scan_v6.push_back("488"); //jfcalib-1
scan_v6.push_back("490"); //jfcalib-2
scan_v6.push_back("491"); //jfcalib-3
scan_v6.push_back("492"); //jfcalib-4
scan_v6.push_back("493"); //jfcalib-5
scan_v6.push_back("494"); //jfcalib-0
scan_v6.push_back("495"); //jfcalib-1
scan_v6.push_back("497"); //jfcalib-2
scan_v6.push_back("498"); //jfcalib-3
scan_v6.push_back("501"); //jfcalib-4
scan_v6.push_back("502"); //jfcalib-5
scan_v6.push_back("036"); //jfcalib-0
scan_v6.push_back("560"); //jfcalib-1
scan_v6.push_back("536"); //jfcalib-2
scan_v6.push_back("535"); //jfcalib-3
scan_v6.push_back("527"); //jfcalib-4
scan_v6.push_back("522"); //jfcalib-5
scan_v6.push_back("571"); //jfcalib-0
scan_v6.push_back("571_fullspeed"); //jfcalib-0
scan_v6.push_back("573"); //jfcalib-1
scan_v6.push_back("573_fullspeed"); //jfcalib-1
scan_v6.push_back("574"); //jfcalib-2
scan_v6.push_back("574_fullspeed"); //jfcalib-2
scan_v6.push_back("575"); //jfcalib-3
scan_v6.push_back("575_fullspeed"); //jfcalib-3
scan_v6.push_back("576"); //jfcalib-4
scan_v6.push_back("576_fullspeed"); //jfcalib-4
scan_v6.push_back("577"); //jfcalib-5
scan_v6.push_back("577_fullspeed"); //jfcalib-5
scan_v6.push_back("563"); //jfcalib-0
scan_v6.push_back("563_fullspeed"); //jfcalib-0
scan_v6.push_back("564"); //jfcalib-1
scan_v6.push_back("564_fullspeed"); //jfcalib-1
scan_v6.push_back("566"); //jfcalib-2
scan_v6.push_back("566_fullspeed"); //jfcalib-2
scan_v6.push_back("567"); //jfcalib-3
scan_v6.push_back("567_fullspeed"); //jfcalib-3
scan_v6.push_back("568"); //jfcalib-4
scan_v6.push_back("568_fullspeed"); //jfcalib-4
scan_v6.push_back("569"); //jfcalib-5
scan_v6.push_back("569_fullspeed"); //jfcalib-5
scan_v6.push_back("579"); //QUAD thin
scan_v6.push_back("579_fullspeed"); //QUAD thin
scan_v6.push_back("454_fullspeed"); //jfcalib-0
scan_v6.push_back("466_fullspeed"); //jfcalib-1
scan_v6.push_back("536_fullspeed"); //jfcalib-2
scan_v6.push_back("561"); //jfcalib-4
scan_v6.push_back("561_fullspeed"); //jfcalib-3
scan_v6.push_back("562"); //jfcalib-4
scan_v6.push_back("562_fullspeed"); //jfcalib-4
scan_v6.push_back("038"); //jfcalib-5
scan_v6.push_back("038_fullspeed"); //jfcalib-5
scan_v6.push_back("431_JF12");
scan_v6.push_back("428_fullspeed"); //jfcalib-0
scan_v6.push_back("443_fullspeed"); //jfcalib-1
scan_v6.push_back("444_fullspeed"); //jfcalib-2
scan_v6.push_back("455_fullspeed"); //jfcalib-3
scan_v6.push_back("471_fullspeed"); //jfcalib-4
scan_v6.push_back("493_fullspeed"); //jfcalib-5
scan_v6.push_back("448_fullspeed"); //jfcalib-0
scan_v6.push_back("497_fullspeed"); //jfcalib-1
scan_v6.push_back("498_fullspeed"); //jfcalib-2
scan_v6.push_back("501_fullspeed"); //jfcalib-3
scan_v6.push_back("502_fullspeed"); //jfcalib-4
scan_v6.push_back("527_fullspeed"); //jfcalib-5
scan_v6.push_back("580_fullspeed"); //jfcalib-5
scan_v6.push_back("580"); //jfcalib-5
// filter out points at zero and points with 0 error (mean taken from one entry in hist)
int m_max = 0;
if (find(scan_v2.begin(), scan_v2.end(), module_str) != scan_v2.end()) {
m_max = 27;
} else if (find(scan_v4.begin(), scan_v4.end(), module_str) != scan_v4.end()) {
m_max = 31;
} else if (find(scan_v3.begin(), scan_v3.end(), module_str) != scan_v3.end()) {
m_max = 35;
} else if (find(scan_v5.begin(), scan_v5.end(), module_str) != scan_v5.end()) {
m_max = 51;
} else if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
m_max = 97;
} else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
m_max = 90;
}
std::vector<double> G0_pixel(97, 0);
std::vector<double> G0_err_pixel(97, 0);
std::vector<double> G1_pixel(97, 0);
std::vector<double> G1_err_pixel(97, 0);
std::vector<double> G2_pixel(97, 0);
std::vector<double> G2_err_pixel(97, 0);
char filename[256]; //uncomment for VH 210902
double filter[97];
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);
}
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 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);
}
for (int i = 0; i < 17; i++) {
filter[i+9] = 0.5+(i*0.25);
}
for (int i = 0; i < 17; i++) {
filter[i+9+17] = 5.+(i*2.5);
}
for (int i = 0; i < 8; i++) {
filter[i+9+17+17] = 50.+(i*10.);
}
} 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);
}
// loop 2a
for (int i = 0; i < 11; i++) {
filter[i+9] = 0.5+(i*0.05);
}
// loop 2b
for (int i = 0; i < 10; i++) {
filter[i+9+11] = 1.10+(i*0.10);
}
// loop 2c
for (int i = 0; i < 10; i++) {
filter[i+9+11+10] = 2.25+(i*0.25);
}
// loop 3a
for (int i = 0; i < 11; i++) {
filter[i+9+11+10+10] = 5.+(i*0.5);
}
// loop 3b
for (int i = 0; i < 10; i++) {
filter[i+9+11+10+10+11] = 11.+(i*1.0);
}
// loop 3c
for (int i = 0; i < 10; i++) {
filter[i+9+11+10+10+11+10] = 22.5+(i*2.5);
}
// loop 4a
for (int i = 0; i < 11; i++) {
filter[i+9+11+10+10+11+10+10] = 50.+(i*5.);
}
// loop 4b
for (int i = 0; i < 15; i++) {
filter[i+9+11+10+10+11+10+10+11] = 110.+(i*10.);
}
} 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);
}
for (int i = 0; i < 8; i++) {
filter[i+10] = 3+(i*1);
}
for (int i = 0; i < 25; i++) {
filter[i+10+8] = 12+(i*2);
}
for (int i = 0; i < 47; i++) {
filter[i+10+8+25] = 80+(i*20);
}
} 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);
}
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);
}
} 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);
}
for (int i = 0; i < 9; i++) {
filter[i+9] = 0.5+(i*0.25);
}
for (int i = 0; i < 9; i++) {
filter[i+9+9] = 5.+(i*2.5);
}
}
TCanvas *A2 = new TCanvas("A2","Plot scan",150,10,800,400);
TCanvas *A3 = new TCanvas("A3","Plot G1 fit",150,10,800,400);
TCanvas *A4 = new TCanvas("A4","Plot G1 residuals",150,10,800,400);
TCanvas *A5 = new TCanvas("A5","Plot G2 fit",150,10,800,400);
TCanvas *A6 = new TCanvas("A6","Plot G2 residuals",150,10,800,400);
TCanvas *A7 = new TCanvas("A7","Plot G0 fit",150,10,800,400);
TCanvas *A8 = new TCanvas("A8","Plot G0 residuals",150,10,800,400);
TCanvas *A9 = new TCanvas("A9","Plot G1 linearity",150,10,800,400);
//file name
sprintf(filename,"%sCS_histos_M%s.root",folder_path.c_str(), module_str.c_str());
cout << "Loading file " << filename << endl;
TFile* f = new TFile((const char *)(filename),"READ");
for (int j = 0; j < m_max; j++) {
TH2F* hist0=(TH2F*)f->Get(Form("adcmap_avg_g0_%d",j));
G0_pixel[j]=hist0->GetBinContent((column+1),(row+1));
TH2F* hist0er=(TH2F*)f->Get(Form("adcermap_avg_g0_%d",j));
G0_err_pixel[j]=hist0er->GetBinContent((column+1),(row+1));
TH2F* hist1=(TH2F*)f->Get(Form("adcmap_avg_g1_%d",j));
G1_pixel[j]=hist1->GetBinContent((column+1),(row+1));
TH2F* hist1er=(TH2F*)f->Get(Form("adcermap_avg_g1_%d",j));
G1_err_pixel[j]=hist1er->GetBinContent((column+1),(row+1));
TH2F* hist2=(TH2F*)f->Get(Form("adcmap_avg_g2_%d",j));
G2_pixel[j]=hist2->GetBinContent((column+1),(row+1));
TH2F* hist2er=(TH2F*)f->Get(Form("adcermap_avg_g2_%d",j));
G2_err_pixel[j]=hist2er->GetBinContent((column+1),(row+1));
//cout << "Data for pixel "<< pixel << "is loaded" << endl;
}
A2->SetLeftMargin(0.13);
A2->SetRightMargin(0.05);
vector<double> r0_adc;
vector<double> r0_filter;
vector<double> r0_adcerr;
vector<double> r0_ferr;
vector<double> r1_adc;
vector<double> r1_filter;
vector<double> r1_adcerr;
vector<double> r1_ferr;
vector<double> r2_adc;
vector<double> r2_filter;
vector<double> r2_adcerr;
vector<double> r2_ferr;
for (int j = 0; j < m_max; j++) {
if (G0_pixel[j] != 0) {
r0_filter.push_back(filter[j]);
r0_ferr.push_back(0.);
r0_adc.push_back(G0_pixel[j]);
r0_adcerr.push_back(G0_err_pixel[j]);
}
if (G1_pixel[j] != 0) {
r1_filter.push_back(filter[j]);
r1_ferr.push_back(0.);
r1_adc.push_back(G1_pixel[j]);
r1_adcerr.push_back(G1_err_pixel[j]);
}
if (G2_pixel[j] != 0) {
r2_filter.push_back(filter[j]);
r2_ferr.push_back(0.);
r2_adc.push_back(G2_pixel[j]);
r2_adcerr.push_back(G2_err_pixel[j]);
}
}
cout << "Filter array is created" << endl;
double rangemin_g0 = 0;
double rangemax_g0 = 0;
double rangemin_g1 = 0;
double rangemax_g1 = 0;
double rangemin_g2 = 0;
double rangemax_g2 = 0;
A2->cd();
// define graphs
if (r0_adc.size() > 1) {
grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0]));
grap_g0->SetMarkerStyle(20);
grap_g0->SetMarkerColor(kBlue);
grap_g0->SetLineColor(kBlue);
}
if (r1_adc.size() > 1) {
grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0]));
grap_g1->SetMarkerStyle(20);
grap_g1->SetMarkerColor(kGreen+2);
grap_g1->SetLineColor(kGreen+2);
}
if (r2_adc.size() > 1) {
grap_g2 = new TGraphErrors(r2_adc.size(),&(r2_filter[0]),&(r2_adc[0]),&(r2_ferr[0]),&(r2_adcerr[0]));
grap_g2->SetMarkerStyle(20);
grap_g2->SetMarkerColor(kRed);
grap_g2->SetLineColor(kRed);
}
// plot the datapoints
if (r1_adc.size() > 1) {
if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
grap_g1->GetXaxis()->SetLimits(0.009,5E2);
} else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
grap_g1->GetXaxis()->SetLimits(0.09,1.1E3);
} else {
grap_g1->GetXaxis()->SetLimits(0.009,1.1E2);
}
grap_g1->GetXaxis()->SetTitle("Integration time [#mus]");
grap_g1->GetYaxis()->SetTitle("ADC [ADU]");
grap_g1->GetYaxis()->SetTitleOffset(0.9);
grap_g1->SetMinimum(0);
grap_g1->SetMaximum(16000);
grap_g1->Draw("AP");
if (r2_adc.size() > 1) {
grap_g2->Draw("P");
}
if (r0_adc.size() > 1) {
grap_g0->Draw("P");
}
A2->SetLogx();
A2->Modified();
A2->Update();
}
// define fit ranges and fit
if (r0_adc.size() > 1) {
A7->cd();
rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end());
rangemin_g0 = checkRangeMinForRCEffect(rangemin_g0);
rangemax_g0 = highestPointBeforeSwitching(r0_filter, r1_filter);
if (rangemax_g0 > rangemin_g0) {
fit_g0 = new TF1("fit_g0","[0]+[1]*x",rangemin_g0, rangemax_g0);
fit_g0->SetParameter(0, 3000.);
fit_g0->SetParameter(1, 30000.);
fit_g0->SetLineColor(kBlue);
fit_g0->SetParName(0,"G0 const");
fit_g0->SetParName(1,"G0 grad");
grap_g0->Fit(fit_g0,"QR+","");
grap_g0->SetMinimum(0);
grap_g0->SetMaximum(18000);
grap_g0->GetXaxis()->SetTitle("Integration time [#mus]");
grap_g0->GetYaxis()->SetTitle("ADC [ADU]");
grap_g0->GetYaxis()->SetTitleOffset(0.9);
grap_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
grap_g0->Draw("AP");
fit_g0->Draw("same");
A7->Update();
if (TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats")) {
st0->SetX1NDC(0.2);
st0->SetX2NDC(0.54);
st0->SetY1NDC(0.71);
st0->SetY2NDC(0.90);
st0->SetBorderSize(0);
st0->SetTextColor(kBlue);
A7->SetLogx();
A7->Modified();
A7->Update();
//A7->cd();
}
vector<double> r0_adc_norm;
for (size_t j = 0; j < r0_adc.size(); j++) {
r0_adc_norm.push_back(r0_adc[j] - fit_g0->Eval(r0_filter[j]));
}
A8->cd();
TGraphErrors *norm_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_ferr[0]),&(r0_adcerr[0]));
norm_g0->SetMarkerColor(kBlue);
norm_g0->SetLineColor(kBlue);
TF1* flat_g0 = new TF1("flat_g0","0",rangemin_g0, rangemax_g0);
flat_g0->SetLineColor(kBlue);
norm_g0->GetXaxis()->SetTitle("Integration time [#mus]");
norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]");
norm_g0->GetYaxis()->SetTitleOffset(0.9);
norm_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
norm_g0->SetMinimum(-100);
norm_g0->SetMaximum(100);
norm_g0->Draw("AP");
flat_g0->Draw("same");
A8->Update();
//sprintf(savename,"plots/M%s/%s/pixel_%s_%d_g0norm_M%s.png", module_str.c_str(), plotfolder_str.c_str(),pixel_type.c_str(), i, module_str.c_str());
//mapcanvas->SaveAs((const char *)(savename));
//delete norm_g0;
}
}
A3->cd();
if (r1_adc.size() > 1) {
rangemin_g1 = lowestPointAfterSwitching(r1_filter, r0_filter);
//cout << rangemin_g1 << endl;
rangemin_g1 = checkRangeMinForRCEffect(rangemin_g1);
rangemax_g1 = highestPointBeforeSwitching(r1_filter, r2_filter);
//cout << rangemin_g1 << endl;
//cout << rangemax_g1 << endl;
if (rangemax_g1 > rangemin_g1) {
fit_g1 = new TF1("fit_g1","[0]+[1]*x",rangemin_g1,rangemax_g1);
fit_g1->SetParameter(0, 10000.);
fit_g1->SetParameter(1, -0.1);
fit_g1->SetLineColor(kGreen+2);
fit_g1->SetParName(0,"G1 const");
fit_g1->SetParName(1,"G1 grad");
grap_g1->Fit(fit_g1,"QR","");
TFitResultPtr r1 = grap_g1->Fit(fit_g1,"S");
par1 = r1->Parameter(1);
grap_g1->SetMinimum(0);
grap_g1->SetMaximum(16000);
grap_g1->GetXaxis()->SetTitle("Integration time [#mus]");
grap_g1->GetYaxis()->SetTitle("ADC [ADU]");
grap_g1->GetYaxis()->SetTitleOffset(0.9);
grap_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
grap_g1->Draw("AP");
fit_g1->Draw("same");
A3->Update();
TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats");
st0->SetX1NDC(0.2);
st0->SetX2NDC(0.54);
st0->SetY1NDC(0.17);
st0->SetY2NDC(0.36);
st0->SetBorderSize(0);
st0->SetTextColor(kGreen+2);
A3->SetLogx();
A3->Modified();
A3->Update();
vector<double> r1_adc_norm;
for (size_t j = 0; j < r1_adc.size(); j++) {
r1_adc_norm.push_back(r1_adc[j] - fit_g1->Eval(r1_filter[j]));
}
A4->cd();
TGraphErrors *norm_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcerr[0]));
norm_g1->SetMarkerColor(kGreen+2);
norm_g1->SetLineColor(kGreen+2);
TF1* flat_g1 = new TF1("flat_g1","0",rangemin_g1,rangemax_g1);
flat_g1->SetLineColor(kGreen+2);
TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_p02pc->SetParameter(0,fit_g1->GetParameter(0)/500.);
lin_g1_p02pc->SetParameter(1,fit_g1->GetParameter(1)/500.);
lin_g1_p02pc->SetLineColor(kOrange);
TF1* lin_g1_p05pc = new TF1("lin_g1_p05pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_p05pc->SetParameter(0,fit_g1->GetParameter(0)/200.);
lin_g1_p05pc->SetParameter(1,fit_g1->GetParameter(1)/200.);
lin_g1_p05pc->SetLineColor(kOrange+1);
TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_m02pc->SetParameter(0,fit_g1->GetParameter(0)/-500.);
lin_g1_m02pc->SetParameter(1,fit_g1->GetParameter(1)/-500.);
lin_g1_m02pc->SetLineColor(kOrange);
TF1* lin_g1_m05pc = new TF1("lin_g1_m05pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_m05pc->SetParameter(0,fit_g1->GetParameter(0)/-200.);
lin_g1_m05pc->SetParameter(1,fit_g1->GetParameter(1)/-200.);
lin_g1_m05pc->SetLineColor(kOrange+1);
norm_g1->GetXaxis()->SetTitle("Integration time [#mus]");
norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]");
norm_g1->GetYaxis()->SetTitleOffset(0.9);
norm_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
norm_g1->SetMinimum(1.5*lin_g1_m05pc->Eval(rangemin_g1));
norm_g1->SetMaximum(1.5*lin_g1_p05pc->Eval(rangemin_g1));
norm_g1->Draw("AP");
flat_g1->Draw("same");
lin_g1_p02pc->Draw("same");
lin_g1_p05pc->Draw("same");
lin_g1_m02pc->Draw("same");
lin_g1_m05pc->Draw("same");
norm_g1->Draw("P");
A4->Modified();
A4->Update();
}
}
if(fit_g1){
A9->cd();
TGraph* linearityGraph = new TGraph();
double x, y;
for(size_t i = 0; i < grap_g1->GetN(); i++){
grap_g1->GetPoint(i, x, y);
auto ADU = y;
auto fitADU = fit_g1->Eval(x);
linearityGraph->SetPoint(i, ADU, fitADU - ADU);
}
//linearityGraph->GetXaxis()->SetTitle("Injected signal [mV]");
linearityGraph->GetXaxis()->SetTitle("Injected signal [ADU]");
linearityGraph->GetYaxis()->SetTitle("fit - signal [ADU]");
linearityGraph->SetMarkerStyle(20);
linearityGraph->SetMarkerColor(kGreen+2);
linearityGraph->SetLineColor(kGreen+2);
linearityGraph->Draw("AP");
A9->Update();
}
A5->cd();
if (r2_adc.size() > 1) {
rangemin_g2 = lowestPointAfterSwitching(r2_filter, r1_filter);
rangemin_g2 = checkRangeMinForRCEffect(rangemin_g2);
rangemax_g2 = *max_element(r2_filter.begin(),r2_filter.end());
if (rangemax_g2 > rangemin_g2) {
fit_g2 = new TF1("fit_g2","[0]+[1]*x",rangemin_g2, rangemax_g2);
fit_g2->SetParameter(0, 10000.);
fit_g2->SetParameter(1, -0.01);
fit_g2->SetLineColor(kRed);
fit_g2->SetParName(0,"G2 const");
fit_g2->SetParName(1,"G2 grad");
grap_g2->Fit(fit_g2,"QR+","");
TFitResultPtr r = grap_g2->Fit(fit_g2,"S");
par2 = r->Parameter(1);
grap_g2->SetMinimum(0);
grap_g2->SetMaximum(16000);
grap_g2->GetXaxis()->SetTitle("Integration time [#mus]");
grap_g2->GetYaxis()->SetTitle("ADC [ADU]");
grap_g2->GetYaxis()->SetTitleOffset(0.9);
grap_g2->GetXaxis()->SetLimits(*min_element(r2_filter.begin(),r2_filter.end()),*max_element(r2_filter.begin(),r2_filter.end()));
grap_g2->Draw("AP");
fit_g2->Draw("same");
A5->Update();
TPaveStats *st0 = (TPaveStats*)grap_g2->FindObject("stats");
st0->SetX1NDC(0.2);
st0->SetX2NDC(0.54);
st0->SetY1NDC(0.17);
st0->SetY2NDC(0.36);
st0->SetBorderSize(0);
st0->SetTextColor(kRed);
A5->SetLogx();
A5->Modified();
A5->Update();
vector<double> r2_adc_norm;
for (size_t j = 0; j < r2_adc.size(); j++) {
r2_adc_norm.push_back(r2_adc[j] - fit_g2->Eval(r2_filter[j]));
}
A6->cd();
TGraphErrors *norm_g2 = new TGraphErrors(r2_adc.size(),&(r2_filter[0]),&(r2_adc_norm[0]),&(r2_ferr[0]),&(r2_adcerr[0]));
norm_g2->SetMarkerColor(kRed);
norm_g2->SetLineColor(kRed);
TF1* flat_g2 = new TF1("flat_g2","0",rangemin_g2, rangemax_g2);
flat_g2->SetLineColor(kRed);
TF1* lin_g2_p02pc = new TF1("lin_g2_p02pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
lin_g2_p02pc->SetParameter(0,fit_g2->GetParameter(0)/500.);
lin_g2_p02pc->SetParameter(1,fit_g2->GetParameter(1)/500.);
lin_g2_p02pc->SetLineColor(kOrange);
TF1* lin_g2_p05pc = new TF1("lin_g2_p05pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
lin_g2_p05pc->SetParameter(0,fit_g2->GetParameter(0)/200.);
lin_g2_p05pc->SetParameter(1,fit_g2->GetParameter(1)/200.);
lin_g2_p05pc->SetLineColor(kOrange+1);
TF1* lin_g2_m02pc = new TF1("lin_g2_m02pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
lin_g2_m02pc->SetParameter(0,fit_g2->GetParameter(0)/-500.);
lin_g2_m02pc->SetParameter(1,fit_g2->GetParameter(1)/-500.);
lin_g2_m02pc->SetLineColor(kOrange);
TF1* lin_g2_m05pc = new TF1("lin_g2_m05pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
lin_g2_m05pc->SetParameter(0,fit_g2->GetParameter(0)/-200.);
lin_g2_m05pc->SetParameter(1,fit_g2->GetParameter(1)/-200.);
lin_g2_m05pc->SetLineColor(kOrange+1);
norm_g2->GetXaxis()->SetTitle("Integration time [#mus]");
norm_g2->GetYaxis()->SetTitle("Normalised ADC [ADU]");
norm_g2->GetYaxis()->SetTitleOffset(0.9);
norm_g2->GetXaxis()->SetLimits(*min_element(r2_filter.begin(),r2_filter.end()),*max_element(r2_filter.begin(),r2_filter.end()));
norm_g2->SetMinimum(1.5*lin_g2_m05pc->Eval(rangemin_g2));
norm_g2->SetMaximum(1.5*lin_g2_p05pc->Eval(rangemin_g2));
norm_g2->Draw("AP");
flat_g2->Draw("same");
lin_g2_p02pc->Draw("same");
lin_g2_p05pc->Draw("same");
lin_g2_m02pc->Draw("same");
lin_g2_m05pc->Draw("same");
norm_g2->Draw("P");
A6->Modified();
A6->Update();
}
}
// cout << par1/par2 << endl;
rootapp.Run();
return 0;
}
Reference in New Issue View Git Blame Copy Permalink