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

SC analysis compiles, intermediate switching to std::vector

Browse Source
This commit is contained in:
hinger_v
2022-05-09 17:14:48 +02:00
parent 76183e9a6a
commit 25fc057033
1 changed files with 14 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
CuFluo_analysis_sc.cpp
View File

@ -119,21 +119,23 @@ int main(int argc, char* argv[]) {
jungfrauFile *thisfile = new jungfrauFile();
jungfrauPedestal *pedestalObject_SC[16];
TH2F *pedestalsG0[16];
TH2F *pedeRMSG0[16];
vector<jungfrauPedestal*> pedestalObject_SC(16); //should call default constructor 16 times
vector<TH2F*> pedestalsG0;
pedestalsG0.reserve(16);
vector<TH2F*> pedeRMSG0;
pedeRMSG0.reserve(16);
for (int sci = 0; sci < 16; ++sci) {
pedestalObject_SC[sci] = new jungfrauPedestal();
//pedestalObject_SC[sci] = new jungfrauPedestal();
pedestalObject_SC[sci]->pedestalSetNFrames(100);
Char_t *pedehistoname = new Char_t[50];
snprintf( pedehistoname, 50, "pedestalsG0_sc%d", sci );
pedestalsG0[sci] = new TH2F( pedehistoname, "", NC, -0.5, NC-0.5, NR, -0.5, NR-0.5 );
pedestalsG0.push_back( new TH2F( pedehistoname, "", NC, -0.5, NC-0.5, NR, -0.5, NR-0.5 ) );
snprintf( pedehistoname, 50, "pedeRMSG0_sc%d", sci );
pedeRMSG0[sci] = new TH2F( pedehistoname, "", NC, -0.5, NC-0.5, NR, -0.5, NR-0.5 );
pedeRMSG0.push_back( new TH2F( pedehistoname, "", NC, -0.5, NC-0.5, NR, -0.5, NR-0.5 ) );
}
static uint16_t pedestals16_G0_start[16][NCH]; //I suppose, this is to track pedestal shifting over the course of data taking
static double pedeRMS16_G0[16][NCH];
static vector< vector<uint16_t> > pedestals16_G0_start( 16, vector<uint16_t>(NCH) ); //I suppose, this is to track pedestal shifting over the course of data taking
static vector< vector<double> > pedeRMS16_G0( 16, vector<double>(NCH) );
for (int pedefilei = 0; pedefilei < pedefilen; ++pedefilei) {
@ -168,7 +170,7 @@ int main(int argc, char* argv[]) {
}
}
pedestalObject_SC[sci]->pedestalData((uint16_t*)(&pedestals16_G0_start[sci])); //this loads the pedestal data into the array pedestals16_G0_start
pedestalObject_SC[sci]->pedestalRMSData(pedeRMS16_G0[sci]); //same here
pedestalObject_SC[sci]->pedestalRMSData(pedeRMS16_G0[sci].data()); //same here
pedestalObject_SC[sci]->pedestalResetUpdates();
//pedestalObject[sci]->pedestalClear(); //I don't need to clear if I only calculate pedestal once.
}
@ -201,6 +203,8 @@ int main(int argc, char* argv[]) {
int adc2d_nbin= 1200;
if (isJF11) adc2d_nbin= 1600;
//HERE!!!
TH2I *adc2d[8][16];
/*
@ -242,7 +246,7 @@ int main(int argc, char* argv[]) {
uint16_t* imagedptr = thisfile->getFrameDataHandle();
uint64_t scnumber = thisfile->currentSCnumber();
pedestalObject_SC[scnumber]->addG0FrameToPedestalCalculationWThreshold( imagedptr, pedestalObject_SC[scnumber], pedeRMS16_G0[scnumber] );
pedestalObject_SC[scnumber]->addG0FrameToPedestalCalculationWThreshold( imagedptr, pedestalObject_SC[scnumber], pedeRMS16_G0[scnumber].data() );
for (int i = 0; i < NCH; i++) {