finished UI for the coincident analyzer, need to check backward event build
This commit is contained in:
splitPoleDAQ
parent
a34bf75d5f
commit
ab54bbbc09
|
|
@ -739,6 +739,12 @@ void MainWindow::CloseDigitizers(){
|
|||
histThread = nullptr;
|
||||
}
|
||||
|
||||
if( onlineAnalyzer ){
|
||||
onlineAnalyzer->close();
|
||||
delete onlineAnalyzer;
|
||||
onlineAnalyzer = nullptr;
|
||||
}
|
||||
|
||||
|
||||
if( canvas ){
|
||||
canvas->close();
|
||||
|
|
|
|||
|
|
@ -255,7 +255,7 @@ void MultiBuilder::BuildEvents(bool isFinal, bool skipTrace, bool verbose){
|
|||
|
||||
void MultiBuilder::BuildEventsBackWard(int maxNumEvent, bool verbose){
|
||||
|
||||
//skip trace, and only build for 100 events max
|
||||
//skip trace, and only build for maxNumEvent events max
|
||||
|
||||
// remember the end of DataIndex, prevent over build
|
||||
for( int k = 0; k < nData; k++){
|
||||
|
|
|
|||
|
|
@ -24,8 +24,15 @@ public:
|
|||
influx = new InfluxDB("https://fsunuc.physics.fsu.edu/influx/");
|
||||
dataBaseName = "testing";
|
||||
|
||||
allowSignalSlot = false;
|
||||
SetUpCanvas();
|
||||
|
||||
LoadHistRange();
|
||||
|
||||
}
|
||||
|
||||
~CoincidentAnalyzer(){
|
||||
SaveHistRange();
|
||||
}
|
||||
|
||||
void SetUpCanvas();
|
||||
|
|
@ -40,6 +47,8 @@ private:
|
|||
|
||||
MultiBuilder *evtbder;
|
||||
|
||||
bool allowSignalSlot;
|
||||
|
||||
// declaie histograms
|
||||
Histogram2D * h2D;
|
||||
Histogram1D * h1;
|
||||
|
|
@ -64,6 +73,9 @@ private:
|
|||
RComboBox * aDigi;
|
||||
RComboBox * aCh;
|
||||
|
||||
void SaveHistRange();
|
||||
void LoadHistRange();
|
||||
|
||||
};
|
||||
|
||||
|
||||
|
|
@ -78,6 +90,7 @@ inline void CoincidentAnalyzer::SetUpCanvas(){
|
|||
boxLayout->setAlignment(Qt::AlignTop | Qt::AlignLeft);
|
||||
box->setLayout(boxLayout);
|
||||
|
||||
{
|
||||
chkRunAnalyzer = new QCheckBox("Run Analyzer", this);
|
||||
boxLayout->addWidget(chkRunAnalyzer, 0, 0);
|
||||
|
||||
|
|
@ -95,7 +108,7 @@ inline void CoincidentAnalyzer::SetUpCanvas(){
|
|||
chkBackWardBuilding = new QCheckBox("Use Backward builder", this);
|
||||
boxLayout->addWidget(chkBackWardBuilding, 1, 0);
|
||||
|
||||
QLabel * lbBKWindow = new QLabel("No. Backward Event", this);
|
||||
QLabel * lbBKWindow = new QLabel("Max No. Backward Event", this);
|
||||
lbBKWindow->setAlignment(Qt::AlignRight | Qt::AlignCenter);
|
||||
boxLayout->addWidget(lbBKWindow, 1, 1);
|
||||
sbBackwardCount = new RSpinBox(this, 0);
|
||||
|
|
@ -127,7 +140,9 @@ inline void CoincidentAnalyzer::SetUpCanvas(){
|
|||
connect(sbBuildWindow, &RSpinBox::valueChanged, this, [=](double value){
|
||||
evtbder->SetTimeWindow((int)value);
|
||||
});
|
||||
}
|
||||
|
||||
{
|
||||
QFrame *separator = new QFrame(box);
|
||||
separator->setFrameShape(QFrame::HLine);
|
||||
separator->setFrameShadow(QFrame::Sunken);
|
||||
|
|
@ -146,7 +161,7 @@ inline void CoincidentAnalyzer::SetUpCanvas(){
|
|||
lbXCh->setAlignment(Qt::AlignRight | Qt::AlignCenter);
|
||||
boxLayout->addWidget(lbXCh, 4, 2);
|
||||
xCh = new RComboBox(this);
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) xCh->addItem("Ch-" + QString::number(i));
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) xCh->addItem("Ch-" + QString::number(i), i);
|
||||
boxLayout->addWidget(xCh, 4, 3);
|
||||
|
||||
QLabel * lbYDigi = new QLabel("Y-Digi", this);
|
||||
|
|
@ -162,15 +177,59 @@ inline void CoincidentAnalyzer::SetUpCanvas(){
|
|||
lbYCh->setAlignment(Qt::AlignRight | Qt::AlignCenter);
|
||||
boxLayout->addWidget(lbYCh, 5, 2);
|
||||
yCh = new RComboBox(this);
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) yCh->addItem("Ch-" + QString::number(i));
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) yCh->addItem("Ch-" + QString::number(i), i);
|
||||
boxLayout->addWidget(yCh, 5, 3);
|
||||
|
||||
connect(xDigi, &RComboBox::currentIndexChanged, this, [=](){
|
||||
allowSignalSlot = false;
|
||||
xCh->clear();
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) xCh->addItem("Ch-" + QString::number(i), i);
|
||||
allowSignalSlot = true;
|
||||
|
||||
int bd = xDigi->currentData().toInt();
|
||||
int ch = xCh->currentData().toInt();
|
||||
h2D->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h2D->UpdatePlot();
|
||||
|
||||
});
|
||||
|
||||
connect(xCh, &RComboBox::currentIndexChanged, this, [=](){
|
||||
if( !allowSignalSlot) return;
|
||||
int bd = xDigi->currentData().toInt();
|
||||
int ch = xCh->currentData().toInt();
|
||||
h2D->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h2D->UpdatePlot();
|
||||
});
|
||||
|
||||
connect(yDigi, &RComboBox::currentIndexChanged, this, [=](){
|
||||
allowSignalSlot = false;
|
||||
yCh->clear();
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) yCh->addItem("Ch-" + QString::number(i), i);
|
||||
allowSignalSlot = true;
|
||||
|
||||
int bd = yDigi->currentData().toInt();
|
||||
int ch = yCh->currentData().toInt();
|
||||
h2D->SetYTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h2D->UpdatePlot();
|
||||
|
||||
});
|
||||
|
||||
connect(yCh, &RComboBox::currentIndexChanged, this, [=](){
|
||||
if( !allowSignalSlot) return;
|
||||
int bd = yDigi->currentData().toInt();
|
||||
int ch = yCh->currentData().toInt();
|
||||
h2D->SetYTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h2D->UpdatePlot();
|
||||
});
|
||||
|
||||
}
|
||||
|
||||
{
|
||||
QFrame *separator1 = new QFrame(box);
|
||||
separator1->setFrameShape(QFrame::HLine);
|
||||
separator1->setFrameShadow(QFrame::Sunken);
|
||||
boxLayout->addWidget(separator1, 6, 0, 1, 4);
|
||||
|
||||
|
||||
QLabel * lbaDigi = new QLabel("ID-Digi", this);
|
||||
lbaDigi->setAlignment(Qt::AlignRight | Qt::AlignCenter);
|
||||
boxLayout->addWidget(lbaDigi, 7, 0);
|
||||
|
|
@ -178,15 +237,60 @@ inline void CoincidentAnalyzer::SetUpCanvas(){
|
|||
for(unsigned int i = 0; i < nDigi; i ++ ){
|
||||
aDigi->addItem("Digi-" + QString::number(digi[i]->GetSerialNumber()), i);
|
||||
}
|
||||
boxLayout->addWidget(yDigi, 7, 1);
|
||||
boxLayout->addWidget(aDigi, 7, 1);
|
||||
|
||||
QLabel * lbaCh = new QLabel("1D-Ch", this);
|
||||
lbaCh->setAlignment(Qt::AlignRight | Qt::AlignCenter);
|
||||
boxLayout->addWidget(lbaCh, 7, 2);
|
||||
aCh = new RComboBox(this);
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) aCh->addItem("Ch-" + QString::number(i));
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) aCh->addItem("Ch-" + QString::number(i), i);
|
||||
boxLayout->addWidget(aCh, 7, 3);
|
||||
|
||||
connect(aDigi, &RComboBox::currentIndexChanged, this, [=](){
|
||||
allowSignalSlot = false;
|
||||
aCh->clear();
|
||||
for( int i = 0; i < digi[0]->GetNumInputCh(); i++) aCh->addItem("Ch-" + QString::number(i), i);
|
||||
allowSignalSlot = true;
|
||||
|
||||
int bd = aDigi->currentData().toInt();
|
||||
int ch = aCh->currentData().toInt();
|
||||
h1->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h1->UpdatePlot();
|
||||
h1g->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h1g->UpdatePlot();
|
||||
|
||||
});
|
||||
|
||||
connect(aCh, &RComboBox::currentIndexChanged, this, [=](){
|
||||
if( !allowSignalSlot) return;
|
||||
int bd = aDigi->currentData().toInt();
|
||||
int ch = aCh->currentData().toInt();
|
||||
h1->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h1->UpdatePlot();
|
||||
h1g->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h1g->UpdatePlot();
|
||||
});
|
||||
|
||||
}
|
||||
|
||||
{
|
||||
QFrame *separator1 = new QFrame(box);
|
||||
separator1->setFrameShape(QFrame::HLine);
|
||||
separator1->setFrameShadow(QFrame::Sunken);
|
||||
boxLayout->addWidget(separator1, 8, 0, 1, 4);
|
||||
|
||||
QPushButton * bnClearHist = new QPushButton("Clear All Hist.");
|
||||
boxLayout->addWidget(bnClearHist, 9, 1);
|
||||
|
||||
connect(bnClearHist, &QPushButton::clicked, this, [=](){
|
||||
h2D->Clear();
|
||||
h1->Clear();
|
||||
h1g->Clear();
|
||||
hMulti->Clear();
|
||||
});
|
||||
|
||||
}
|
||||
|
||||
|
||||
}
|
||||
|
||||
|
|
@ -199,12 +303,27 @@ inline void CoincidentAnalyzer::SetUpCanvas(){
|
|||
//layout is inheriatge from Analyzer
|
||||
layout->addWidget(h2D, 1, 0, 2, 1);
|
||||
|
||||
h1 = new Histogram1D("1D Plot", "XXX", 300, 30, 70, this);
|
||||
int bd = xDigi->currentData().toInt();
|
||||
int ch = xCh->currentData().toInt();
|
||||
h2D->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
bd = yDigi->currentData().toInt();
|
||||
ch = yCh->currentData().toInt();
|
||||
h2D->SetYTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h2D->UpdatePlot();
|
||||
|
||||
|
||||
h1 = new Histogram1D("1D Plot", "XXX", 300, 0, 5000, this);
|
||||
h1->SetColor(Qt::darkGreen);
|
||||
h1->AddDataList("Test", Qt::red); // add another histogram in h1, Max Data List is 10
|
||||
bd = aDigi->currentData().toInt();
|
||||
ch = aCh->currentData().toInt();
|
||||
h1->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h1->UpdatePlot();
|
||||
layout->addWidget(h1, 1, 1);
|
||||
|
||||
h1g = new Histogram1D("1D Plot (PID gated)", "XXX", 300, -2, 10, this);
|
||||
h1g = new Histogram1D("1D Plot (PID gated)", "XXX", 300, 0, 5000, this);
|
||||
h1g->SetXTitle("Digi-" + QString::number(digi[bd]->GetSerialNumber()) + ", Ch-" + QString::number(ch));
|
||||
h1g->UpdatePlot();
|
||||
layout->addWidget(h1g, 2, 1);
|
||||
|
||||
layout->setColumnStretch(0, 1);
|
||||
|
|
@ -229,6 +348,16 @@ inline void CoincidentAnalyzer::UpdateHistograms(){
|
|||
unsigned long long tMin[nCut] = {0xFFFFFFFFFFFFFFFF}, tMax[nCut] = {0};
|
||||
unsigned int count[nCut]={0};
|
||||
|
||||
//============ Get the channel to plot
|
||||
int a_bd = aDigi->currentData().toInt();
|
||||
int a_ch = aCh->currentData().toInt();
|
||||
|
||||
int x_bd = xDigi->currentData().toInt();
|
||||
int x_ch = xCh->currentData().toInt();
|
||||
|
||||
int y_bd = yDigi->currentData().toInt();
|
||||
int y_ch = yCh->currentData().toInt();
|
||||
|
||||
//============ Processing data and fill histograms
|
||||
long eventIndex = evtbder->eventIndex;
|
||||
long eventStart = eventIndex - eventBuilt + 1;
|
||||
|
|
@ -236,27 +365,39 @@ inline void CoincidentAnalyzer::UpdateHistograms(){
|
|||
|
||||
for( long i = eventStart ; i <= eventIndex; i ++ ){
|
||||
std::vector<Hit> event = evtbder->events[i];
|
||||
//printf("-------------- %ld\n", i);
|
||||
|
||||
hMulti->Fill((int) event.size());
|
||||
//if( event.size() < 9 ) return;
|
||||
if( event.size() == 0 ) return;
|
||||
|
||||
|
||||
int aE = -1;
|
||||
int xE = -1, yE = -1;
|
||||
unsigned long long xT = 0;
|
||||
for( int k = 0; k < (int) event.size(); k++ ){
|
||||
//event[k].Print();
|
||||
if( event[k].bd == a_bd && event[k].ch == a_ch) {
|
||||
h1->Fill(event[k].energy);
|
||||
aE = event[k].energy;
|
||||
}
|
||||
|
||||
//check events inside any Graphical cut and extract the rate, using tSR only
|
||||
if( event[k].bd == x_bd && event[k].ch == x_ch) {
|
||||
xE = event[k].energy;
|
||||
xT = event[k].timestamp;
|
||||
}
|
||||
if( event[k].bd == y_bd && event[k].ch == y_ch) yE = event[k].energy;
|
||||
}
|
||||
|
||||
if( xE >= 0 && yE >= 0 ) h2D->Fill(xE, yE);
|
||||
|
||||
//check events inside any Graphical cut and extract the rate
|
||||
for(int p = 0; p < cutList.count(); p++ ){
|
||||
if( cutList[p].isEmpty() ) continue;
|
||||
// if( cutList[p].containsPoint(QPointF(hit.eSL, hit.eSR), Qt::OddEvenFill) ){
|
||||
// if( hit.tSR < tMin[p] ) tMin[p] = hit.tSR;
|
||||
// if( hit.tSR > tMax[p] ) tMax[p] = hit.tSR;
|
||||
// count[p] ++;
|
||||
// //printf(".... %d \n", count[p]);
|
||||
// if( p == 0 ) h1g->Fill(hit.eSR);
|
||||
// }
|
||||
if( cutList[p].containsPoint(QPointF(xE, yE), Qt::OddEvenFill) && xE >= 0 && yE >= 0 ){
|
||||
if( xT < tMin[p] ) tMin[p] = xT;
|
||||
if( xT > tMax[p] ) tMax[p] = xT;
|
||||
count[p] ++;
|
||||
//printf(".... %d \n", count[p]);
|
||||
if( p == 0 && aE >= 0 ) h1g->Fill(aE); // only for the 1st gate
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
|
@ -265,9 +406,9 @@ inline void CoincidentAnalyzer::UpdateHistograms(){
|
|||
hMulti->UpdatePlot();
|
||||
h1g->UpdatePlot();
|
||||
|
||||
QList<QString> cutNameList = h2D->GetCutNameList();
|
||||
for( int p = 0; p < cutList.count(); p ++){
|
||||
if( cutList[p].isEmpty() ) continue;
|
||||
// QList<QString> cutNameList = h2D->GetCutNameList();
|
||||
// for( int p = 0; p < cutList.count(); p ++){
|
||||
// if( cutList[p].isEmpty() ) continue;
|
||||
// double dT = (tMax[p]-tMin[p]) * tick2ns / 1e9; // tick to sec
|
||||
// double rate = count[p]*1.0/(dT);
|
||||
//printf("%llu %llu, %f %d\n", tMin[p], tMax[p], dT, count[p]);
|
||||
|
|
@ -276,8 +417,15 @@ inline void CoincidentAnalyzer::UpdateHistograms(){
|
|||
// influx->AddDataPoint("Cut,name=" + cutNameList[p].toStdString()+ " value=" + std::to_string(rate));
|
||||
// influx->WriteData(dataBaseName);
|
||||
// influx->ClearDataPointsBuffer();
|
||||
// }
|
||||
|
||||
}
|
||||
|
||||
inline void CoincidentAnalyzer::SaveHistRange(){
|
||||
|
||||
}
|
||||
inline void CoincidentAnalyzer::LoadHistRange(){
|
||||
|
||||
}
|
||||
|
||||
|
||||
|
|
|
|||
Loading…
Reference in New Issue
Block a user