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

Clean up comments in SC scripts

Browse Source
This commit is contained in:
hinger_v
2022-05-19 10:49:02 +02:00
parent ce27817a9d
commit 4364efe87b
1 changed files with 8 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

29
CuFluo_analysis_sc_singlethread_data.cpp
View File

@ -20,7 +20,6 @@
#include "TPaveStats.h" #include "TPaveStats.h"
#include "TLegend.h" #include "TLegend.h"
#include "TPaveText.h" #include "TPaveText.h"
//#include "TObjArray.h"
#include <sys/stat.h> #include <sys/stat.h>
#include <vector> #include <vector>
@ -54,16 +53,13 @@ int main(int argc, char* argv[]) {
string data_loc = argv[3]; string data_loc = argv[3];
string pede_file = argv[4]; string pede_file = argv[4];
string data_file = argv[5]; string data_file = argv[5];
string anadata_loc = argv[6]; //uncomment for VH 210906 string anadata_loc = argv[6];
bool save_to_afs = false; bool save_to_afs = false;
int createHistoFile = 1; int createHistoFile = 1;
//char histoname[128]; char savename[256];
//char savename[128];
//char histoname[256]; // VH 210902
char savename[256]; // VH 210902
char plotfolder[256]; char plotfolder[256];
char rootdatafolder[256]; char rootdatafolder[256];
int filen = 352; int filen = 352;
@ -86,9 +82,7 @@ int main(int argc, char* argv[]) {
// plots/Mxxx/CuFluo/HG0 // plots/Mxxx/CuFluo/HG0
sprintf(savename,"plots/M%s/CuFluo/HG0", module_str.c_str()); sprintf(savename,"plots/M%s/CuFluo/HG0", module_str.c_str());
mkdir(savename, S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH); mkdir(savename, S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH);
// /mnt/pcmoench_jungfrau_data/jungfrau_ana_sophie/Mxxx_CalibAna sprintf(savename,"%s", anadata_loc.c_str());
//sprintf(savename,"/mnt/sls_det_storage/jungfrau_data1/jungfrau_ana_sophie/M%s_CalibAna", module_str.c_str()); //uncomment for SR
sprintf(savename,"%s", anadata_loc.c_str()); //uncomment for VH 210906
mkdir(savename, S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH); mkdir(savename, S_IRWXU | S_IRWXG | S_IROTH | S_IXOTH);
if ( save_to_afs ) { if ( save_to_afs ) {
@ -133,8 +127,6 @@ int main(int argc, char* argv[]) {
mapcanvas->SetTopMargin(0.08); mapcanvas->SetTopMargin(0.08);
mapcanvas->SetBottomMargin(0.15); mapcanvas->SetBottomMargin(0.15);
//TCanvas* c1 = new TCanvas("c1","");
/*****************************************************************************************************************/ /*****************************************************************************************************************/
/* NOTE: Since arrays, as they were used before, tend to cause memory leaks if they grow too large on the stack, */ /* NOTE: Since arrays, as they were used before, tend to cause memory leaks if they grow too large on the stack, */
/* it might be smarter to realize all array objects with std::vector instead. */ /* it might be smarter to realize all array objects with std::vector instead. */
@ -147,7 +139,7 @@ int main(int argc, char* argv[]) {
jungfrauPedestal* pedestalObject_SC = new jungfrauPedestal(); jungfrauPedestal* pedestalObject_SC = new jungfrauPedestal();
pedestalObject_SC->pedestalSetNFrames(100); pedestalObject_SC->pedestalSetNFrames(100);
vector<uint16_t> pedestals16_G0_start( NCH ); //I suppose, this is to track pedestal shifting over the course of data taking vector<uint16_t> pedestals16_G0_start( NCH ); //Track pedestal shifting over the course of data taking
vector<double> pedeRMS16_G0( NCH ); vector<double> pedeRMS16_G0( NCH );
//cout << "Last element of pedeRMS16_G0 = " << pedeRMS16_G0[NCH-1] << endl; //cout << "Last element of pedeRMS16_G0 = " << pedeRMS16_G0[NCH-1] << endl;
@ -169,7 +161,7 @@ int main(int argc, char* argv[]) {
for (int pedefilei = 0; pedefilei < pedefilen; ++pedefilei) { for (int pedefilei = 0; pedefilei < pedefilen; ++pedefilei) {
// open pede file // open pede file
sprintf( savename, "%s/%s_%%6.6d.dat", data_loc.c_str(), pede_file.c_str() ); //VH: note, this adds a double slash in the filepath sprintf( savename, "%s/%s_%%6.6d.dat", data_loc.c_str(), pede_file.c_str() );
thisfile->open( (char*)savename, pedefilei ); thisfile->open( (char*)savename, pedefilei );
//count events in file //count events in file
@ -254,8 +246,6 @@ int main(int argc, char* argv[]) {
TH1D* adcpc_spec; TH1D* adcpc_spec;
TH2F* pede_updates; TH2F* pede_updates;
TH2F* pede_diff; TH2F* pede_diff;
//cout << "Initialize vectors of histograms..." << endl;
//initialize //initialize
Char_t *_histoname = new Char_t[50]; Char_t *_histoname = new Char_t[50];
@ -266,12 +256,10 @@ int main(int argc, char* argv[]) {
snprintf( _histoname, 50, "pede_diff_sc%d", this_storagecell ); snprintf( _histoname, 50, "pede_diff_sc%d", this_storagecell );
pede_diff = new TH2F( _histoname, "", NC, -0.5, NC-0.5, NR, -0.5, NR-0.5 ); pede_diff = new TH2F( _histoname, "", NC, -0.5, NC-0.5, NR, -0.5, NR-0.5 );
//cout << "Done." << endl;
for (int filei = 0; filei < filen; filei++) { for (int filei = 0; filei < filen; filei++) {
// open data file // open data file
sprintf(savename,"%s/%s_%%6.6d.dat", data_loc.c_str(), data_file.c_str()); //VH: note, this adds a double slash in the filepath sprintf(savename,"%s/%s_%%6.6d.dat", data_loc.c_str(), data_file.c_str());
thisfile->open((char*)savename, filei); thisfile->open((char*)savename, filei);
while (thisfile->readNextFrame()) { while (thisfile->readNextFrame()) {
@ -355,8 +343,7 @@ int main(int argc, char* argv[]) {
} // end of file loop } // end of file loop
//sprintf(savename,"/mnt/sls_det_storage/jungfrau_data1/jungfrau_ana_sophie/M%s_CalibAna/CuFluo_%s_file0to%d.root", module_str.c_str(), gain_str.c_str(), filen-1); sprintf(savename,"%s/CuFluo_%s_sc%d_file0to%d.root", anadata_loc.c_str(), gain_str.c_str(), this_storagecell, filen-1);
sprintf(savename,"%s/CuFluo_%s_sc%d_file0to%d.root", anadata_loc.c_str(), gain_str.c_str(), this_storagecell, filen-1); //uncomment for VH 210906
TFile* saved_file = new TFile((const char *)(savename),"RECREATE"); TFile* saved_file = new TFile((const char *)(savename),"RECREATE");
for ( int adci = 0; adci < 8; ++adci ) { for ( int adci = 0; adci < 8; ++adci ) {
adc2d[adci]->Write(); adc2d[adci]->Write();
@ -370,7 +357,7 @@ int main(int argc, char* argv[]) {
saved_file->Close(); saved_file->Close();
} // end if i creatHistoFile } // end if creatHistoFile
} }