source: tools/ddd/Functions.cpp@ 197

Last change on this file since 197 was 197, checked in by qweitzel, 14 years ago
M0 display adapted to new DRS4 data
File size: 17.1 KB
Line 
1
2#include "GUI.h"
3
4//---------------------------------------------------------------------
5//************************ All functions ****************************
6//-------------------------------------------------------------------
7
8// +++ Open file dialog +++
9void ddd::FileDialog(void) {
10QString Filename = QFileDialog::getOpenFileName(this,
11 "Open raw file", INITIAL_DIRECTORY, "Raw data files (*.raw);; All files (*)");
12 if (Filename != NULL) {
13 FilenameBox->setText(Filename);
14 OpenDatafile();
15 }
16}
17
18// +++ Open selected file and read run header +++
19void ddd::OpenDatafile() {
20
21 if(Socket->state() == QAbstractSocket::ConnectedState) { // do not execute if socket is open
22 MakeConnection();
23 }
24
25 CloseDatafile(); // Close previous file if open
26
27 // Write run header to temporary file
28 ftruncate(fileno(Tmpfile),0);
29 rewind(Tmpfile);
30 switch (RD->OpenDataFile(FilenameBox->text().toAscii().data(), Tmpfile)) {
31 case CTX_FOPEN: QMessageBox::warning(this, "ddd Message","Could not open file.",QMessageBox::Ok);
32 return;
33 case CTX_RHEADER: QMessageBox::warning(this, "ddd Message","Could not read run header.",QMessageBox::Ok);
34 return;
35 case CTX_BSTRUCT: QMessageBox::warning(this, "ddd Message","Could not read board structures.",QMessageBox::Ok);
36 return;
37 default: break;
38 }
39 if (RD->RHeader->MagicNum == MAGICNUM_OPEN) {
40 QMessageBox::warning(this, "ddd Message","Magic number in run header indicates that the file has not been closed properly.",QMessageBox::Ok);
41 }
42 if (RD->RHeader->MagicNum == MAGICNUM_ERROR) {
43 QMessageBox::warning(this, "ddd Message","Magic number in run header indicates that an error occurred while writing the file.",QMessageBox::Ok);
44 }
45
46 rewind(Tmpfile);
47 QTextStream in(Tmpfile);
48 QString text = in.readAll();
49 RunHeaderDisplay->setPlainText(text);
50
51 // Enable spin boxes, set ranges and display first event
52 EventNo->setEnabled(true);
53 ChannelNo->setEnabled(true);
54 BoardNo->setEnabled(true);
55 PixelID->setEnabled(true);
56 M0Start->setEnabled(true);
57 M0Stop->setEnabled(true);
58 PhysPipeAction->setEnabled(true);
59 EventNo->setRange(1, RD->RHeader->Events);
60 M0Display->setEnabled(true);
61 M0Start->setRange(0,(RD->RHeader->Samples)-1);
62 M0Stop->setRange(0,(RD->RHeader->Samples)-1);
63 ChannelNo->setRange(0, RD->RHeader->NChannels*RD->RHeader->NChips-1);
64 BoardNo->setRange(0, RD->RHeader->NBoards-1);
65 DisplayEvent();
66}
67
68// +++ Close data file file if open, delete event header and disable spin boxes and displays +++
69void ddd::CloseDatafile() {
70 if(RD->CloseDataFile()==CTX_OK) {
71 EventNo->setEnabled(false);
72 ChannelNo->setEnabled(false);
73 BoardNo->setEnabled(false);
74 PixelID->setEnabled(false);
75 M0Start->setEnabled(false);
76 M0Stop->setEnabled(false);
77 PhysPipeAction->setEnabled(false);
78 RunHeaderDisplay->clear();
79 EventHeaderDisplay->clear();
80 Signal->hide();
81 }
82}
83
84// +++ Read event header and display event (only called if Datafile is open) +++
85void ddd::DisplayEvent(int) {
86
87 PixelID->setText(PixMap->DRS_to_Pixel(BoardNo->value(), ChannelNo->value()/RD->RHeader->NChannels,
88 ChannelNo->value()%RD->RHeader->NChannels).c_str()); // Translate to pixel ID
89 if(Socket->state() == QAbstractSocket::ConnectedState) return; // do not execute if socket is open
90
91 if (PhysPipeAction->isChecked()) {
92 Graph->setAxisTitle(QwtPlot::xBottom, "Time from start of physical pipeline (ns)");
93 }
94 else {
95 Graph->setAxisTitle(QwtPlot::xBottom, "Time from trigger (ns)");
96 }
97
98 // Read event
99 ftruncate(fileno(Tmpfile),0);
100 rewind(Tmpfile);
101 if (RD->ReadEvent(EventNo->value(), Tmpfile) != CTX_OK) {
102 QMessageBox::warning(this, "ddd Warning","Could not read event.",QMessageBox::Ok);
103 EventHeaderDisplay->clear();
104 return;
105 }
106
107 // Print event header and trigger cell information from event data
108 rewind(Tmpfile);
109 QTextStream in(Tmpfile);
110 QString text = in.readAll();
111
112 text.append("\nTrigger cells: ");
113 for (unsigned int i=0; i<RD->RHeader->NBoards*RD->RHeader->NChips; i++) {
114 QString a;
115 text.append(a.sprintf("%d ", *((int *)RD->Data + i)));
116 }
117 EventHeaderDisplay->setPlainText(text);
118
119 // Case data in double format required by qwt library
120 double *x = new double [RD->RHeader->Samples];
121 double *y = new double [RD->RHeader->Samples];
122
123 for (unsigned int i=0; i<RD->RHeader->Samples; i++) {
124 x[i] = (double) (i/RD->BStruct[BoardNo->value()].NomFreq);
125 if (PhysPipeAction->isChecked()) {
126 y[(i + *((int *) RD->Data + BoardNo->value()*RD->RHeader->NChips+ChannelNo->value()/RD->RHeader->NChannels))%1024] = (double) *((short *) (RD->Data + RD->RHeader->NBoards*RD->RHeader->NChips*sizeof(int)) + BoardNo->value()*RD->RHeader->NChips*RD->RHeader->NChannels * RD->RHeader->Samples + ChannelNo->value()*RD->RHeader->Samples + i)*RD->BStruct[BoardNo->value()].ScaleFactor;
127 }
128 else {
129 y[i] = (double) *((short *) (RD->Data + RD->RHeader->NBoards*RD->RHeader->NChips*sizeof(int)) + BoardNo->value()*RD->RHeader->NChips*RD->RHeader->NChannels * RD->RHeader->Samples+ChannelNo->value()*RD->RHeader->Samples+(i))*RD->BStruct[BoardNo->value()].ScaleFactor;
130 }
131 }
132
133 Signal->setData(x, y, (int) RD->RHeader->Samples);
134 Signal->show();
135 Zoomer->setZoomBase(Signal->boundingRect());
136
137 delete[] x; delete[] y;
138
139 // ************************************
140 // Get data for M0 display (event based)
141 // ************************************
142
143 if (!(M0Display->isEnabled())) return;
144
145 double z[6][6] = {{0}};//36 pixels
146 bool IDerror = false;
147
148 //only interested in M0 data using DRS2
149 //if (RD->RHeader->NChannels == 9) IDerror = true;
150
151 //only interested in M0 data using DRS4
152 if (RD->RHeader->NChannels == 10) IDerror = true;
153
154 for(unsigned int i=0; i<RD->RHeader->NBoards; i++) {//board loop
155 for(unsigned int j=0; j<RD->RHeader->NChips; j++) {//chip loop
156 for(unsigned int k=0; k<RD->RHeader->NChannels; k++) {//channel loop
157
158 //for DRS2 data (requires the correct pixelmap!!!
159 //if( ( (i==0 || i==1) && (j<=1) && (k<=7) ) || ( (i==2) && (j==0) && (k<=3) ) ) {
160
161 //for DRS4 data (requires the correct pixelmap!!!
162 if( ( (i==0 || i==1) && (j<=3) && (k==0 || k==2 || k==4 || k==6) ) || ( (i==2) && (j==0) && (k==0 || k==2 || k==4 || k==6) ) ) {
163
164 //get module, superpixel and pixel number from pixel name
165 std::string pixelname = PixMap->DRS_to_Pixel(i,j,k);
166 if (pixelname == "") {
167 IDerror = true;
168 continue;
169 }
170 char pixelname_copy[256];
171 memset(pixelname_copy,'\0',256);
172 pixelname.copy(pixelname_copy, 256);
173
174 char delim[] = "-";
175 char *buffer = NULL;
176 int module = -1;
177 int superpixel = -1;
178 int pixel = -1;
179
180 buffer = strtok(pixelname_copy, delim);
181 module = atoi(buffer);
182 buffer = strtok(NULL, delim);
183 superpixel = atoi(buffer);
184 buffer = strtok(NULL, delim);
185 pixel = atoi(buffer);
186
187 //usual M0 mapping
188 //int binx = 5-(int((superpixel-1)/3)*2)-(int((pixel%4)/2));
189 //int biny = 5-(((superpixel-1)%3)*2)-(int((pixel-1)/2));
190
191 //M0 upside down
192 int binx = 5-(5-(int((superpixel-1)/3)*2)-(int((pixel%4)/2)));
193 int biny = 5-(5-(((superpixel-1)%3)*2)-(int((pixel-1)/2)));
194
195 //search maximum sample amplitude within user specified window
196 //start bin is always smaller than stop bin (taken care of by updated ranges)
197 int StartBin = (int)(M0Start->value());
198 int StopBin = (int)(M0Stop->value());
199
200 for(int l=StartBin; l<=StopBin; l++){
201
202 float sample = *((short *) (RD->Data + RD->RHeader->NBoards*RD->RHeader->NChips*sizeof(int)) +
203 i*RD->RHeader->NChips*RD->RHeader->NChannels*RD->RHeader->Samples+
204 j*RD->RHeader->NChannels*RD->RHeader->Samples+
205 k*RD->RHeader->Samples+
206 l)*RD->BStruct[i].ScaleFactor;
207
208 if (fabs(sample) > z[binx][biny]) {
209 z[binx][biny]=fabs(sample);
210 }
211
212 }//sample loop
213 }//only M0 data
214 }//channel loop
215 }//chip loop
216 }//board loop
217
218 if(IDerror){
219 QMessageBox::warning(this, "ddd Message","Sorry! The M0 display is not available for this data file because of a pixel ID mismatch.",QMessageBox::Ok);
220 if(M0Window->isVisible()) M0Window->hide();
221 M0Display->setEnabled(false);
222 return;
223 }
224
225 //fill data to M0 display (event based)
226 Signal2D->setData(SpectrogramDataM0(z));
227 Graph2D->axisWidget(QwtPlot::yRight)->setColorMap(Signal2D->data().range(),Signal2D->colorMap());
228 Graph2D->setAxisScale(QwtPlot::yRight,Signal2D->data().range().minValue(),Signal2D->data().range().maxValue() );
229 Graph2D->replot();
230
231 //update ranges for start and stop bin to avoid startbin > stopbin
232 M0Start->setRange(0, M0Stop->value());
233 M0Stop->setRange(M0Start->value(),(RD->RHeader->Samples)-1);
234}
235
236// +++ Open sub window handling the socket interface +++
237void ddd::OpenSocketWindow() {
238
239 if(SocketWindow->isVisible()) SocketWindow->hide();
240 else SocketWindow->show();
241}
242
243// +++ Open sub window for M0 Display +++
244void ddd::OpenM0Window() {
245
246 if(M0Window->isVisible()) M0Window->hide();
247 else M0Window->show();
248}
249
250
251// +++ Acquire data through socket (acquire botton only available if socket exists) +++
252void ddd::GetSignalFromSocket() {
253 char Command[MAX_COM_SIZE];
254
255 GetButton->setEnabled(false);
256 WaitForData = true;
257 sprintf(Command, "read %d %d %d restart", BoardNo->value(), ChannelNo->value()/10, ChannelNo->value()%10);
258 Socket->write(Command);
259}
260
261// Quit application when clicking close button on window
262void ddd::closeEvent(QCloseEvent *) {
263 qApp->quit();
264}
265
266// +++ Connecting or disconnecting from client +++
267void ddd::MakeConnection() {
268
269 if(Socket->state() == QAbstractSocket::ConnectedState) {
270 ManualDisconnect = true;
271 Socket->disconnectFromHost();
272 }
273 else {
274 if (RD->IsFileOpen() && QMessageBox::question(this, "ddd Request","Connecting will close current data file. Proceed?",
275 QMessageBox::No, QMessageBox::Yes) != QMessageBox::Yes) return;
276
277 PhysPipeAction->setChecked(false);
278
279 Socket->connectToHost(IPAddress->text(),Port->value());
280 Connect->setEnabled(false); // While waiting for connection, button not available
281 Socket->waitForConnected(SOCKET_TIMEOUT);
282 Connect->setEnabled(true);
283 if(Socket->state() != QAbstractSocket::ConnectedState)
284 QMessageBox::warning(this, "ddd Message","Could not connect to host.",QMessageBox::Ok);
285 else {
286 Connect->setText("Disconnect");
287 ConnectAction->setText("Disconnect");
288 Port->setEnabled(false);
289 IPAddress->setEnabled(false);
290 Command->setEnabled(true);
291 GetButton->setEnabled(true);
292 ManualDisconnect = false;
293
294 FilenameBox->clear();
295 CloseDatafile();
296
297 M0Window->hide();
298 M0Display->setEnabled(false);
299
300 ChannelNo->setEnabled(true);
301 BoardNo->setEnabled(true);
302 PixelID->setEnabled(true);
303 ChannelNo->setRange(0, 65535);
304 BoardNo->setRange(0, 65535);
305
306 TabWidget->setTabEnabled(1,false);
307 TabWidget->setTabEnabled(2,false);
308
309 RunHeaderDisplay->clear();
310 EventHeaderDisplay->clear();
311 Signal->hide();
312 }
313 }
314}
315
316// +++ Send command to socket (command button available only if socket existing) +++
317void ddd::SendToSocket() {
318 Socket->write(Command->text().toAscii());
319 Command->clear();
320}
321
322// +++ Read data from socket and display +++
323void ddd::ReadFromSocket() {
324 // Check if channel data expected and error message arrived
325 QByteArray Data = Socket->readAll();
326 if (WaitForData && Data.contains("Error")) {
327 WaitForData = false;
328 GetButton->setEnabled(true);
329 QMessageBox::warning(this, "ddd Message","Could not read waveform data from socket.",QMessageBox::Ok);
330 return;
331 }
332
333 // Check if channel data were transmitted, if yes and waiting for data, extract and plot them
334 SocketOutput->insertPlainText(Data);
335 QString Text = SocketOutput->document()->toPlainText().trimmed();
336 if (WaitForData && Text.endsWith(QLatin1String("==END=="))) {
337 // Extract text between ==START== and ==END==
338 QByteArray Data = Text.mid(Text.lastIndexOf("==START==")+9, Text.length() - Text.lastIndexOf("==START==")-16).toAscii();
339
340 char *NextNumber = strtok(Data.data()," "); // Number of entries that follow
341 int Count=0, NumberOfEntries = atoi(NextNumber);
342 double *x = new double [NumberOfEntries];
343 double *y = new double [NumberOfEntries];
344
345 // Convert all entries (separated by a whitespace) to numbers
346 while((NextNumber=strtok(NULL, " "))!=NULL && Count<NumberOfEntries)
347 *(y+Count++) = atof(NextNumber);
348 if (Count==NumberOfEntries && NextNumber!=0)
349 QMessageBox::warning(this, "ddd Message","Found too many numbers in data block, truncated.",QMessageBox::Ok);
350 // Apply sampling frequency and scaling factor
351 for(int i=2; i<Count; i++) {
352 x[i] = (i-2) / y[0];
353 y[i] = y[i] * y[1];
354 }
355 if(NumberOfEntries>2) {
356 Signal->setData(x+2, y+2, NumberOfEntries-2); // Copies data, arrays can be deleted afterwards
357 Signal->show();
358 Zoomer->setZoomBase(Signal->boundingRect());
359 Graph->setAxisTitle(QwtPlot::xBottom, "Time from trigger (ns)");
360 }
361 delete[] x; delete[] y;
362
363 if(ContinuousBox->isChecked()) {
364 usleep(100000); // Wait to limit maximum update rate
365 GetSignalFromSocket();
366 }
367 else {
368 WaitForData = false;
369 GetButton->setEnabled(true);
370 }
371 }
372}
373
374// +++ Reset graph axes to autoscale when fully unzoomed +++
375void ddd::HandleZoom(const QwtDoubleRect &) {
376 if(Zoomer->zoomRectIndex() == 0) {
377 Graph->setAxisAutoScale(QwtPlot::xBottom);
378 Graph->setAxisAutoScale(QwtPlot::yLeft);
379 }
380}
381
382// +++ Disconnect from socket +++
383void ddd::GotDisconnected() {
384 Connect->setText("Connect");
385 ConnectAction->setText("Connect");
386 Port->setEnabled(true);
387 IPAddress->setEnabled(true);
388 Command->setEnabled(false);
389
390 GetButton->setEnabled(false);
391 ChannelNo->setEnabled(false);
392 BoardNo->setEnabled(false);
393 PixelID->setEnabled(false);
394 Signal->hide();
395 TabWidget->setTabEnabled(1, true);
396 TabWidget->setTabEnabled(2, true);
397
398 SocketOutput->clear();
399 if(!ManualDisconnect) QMessageBox::warning(this, "ddd Message","Socket disconnected, maybe host gone.",QMessageBox::Ok);
400}
401
402// +++ Translate pixel ID +++
403void ddd::TranslatePixelID() {
404
405 int Board = PixMap->Pixel_to_DRSboard(PixelID->text().toStdString());
406 int Chip = PixMap->Pixel_to_DRSchip(PixelID->text().toStdString());
407 int Channel = PixMap->Pixel_to_DRSchannel(PixelID->text().toStdString());
408
409 if(Board>=BoardNo->minimum() && Board<=BoardNo->maximum() &&
410 (Chip*10+Channel)>=ChannelNo->minimum() && (Chip*10+Channel)<=ChannelNo->maximum()) {
411 BoardNo->setValue(Board);
412 ChannelNo->setValue(Chip*10+Channel);
413 }
414 else if(Board==999999999) QMessageBox::warning(this, "ddd Message","Pixel ID unknown.",QMessageBox::Ok);
415 else QMessageBox::warning(this, "ddd Message","Pixel ID out of current range.",QMessageBox::Ok);
416}
417
418
419//------------------------------------------------------------------
420//**************************** All menus ***************************
421//------------------------------------------------------------------
422
423void ddd::MenuSave() {
424 QString Filename = QFileDialog::getSaveFileName(this,
425 "Filename of image", "/home/ogrimm/ddd", "Image files (*.bmp *.jpg *.png *.ppm *.tiff *.xbm *.xpm);;All files (*)");
426 if (Filename.length()>0) {
427 QPixmap Pixmap = QPixmap::grabWidget(Graph);
428 if(!Pixmap.save(Filename)) {
429 QMessageBox::warning(this, "ddd Message","Could not write image file.",QMessageBox::Ok);
430 remove(Filename.toAscii().data());
431 }
432 }
433}
434
435void ddd::MenuPrint() {
436 QPrinter *Printer = new QPrinter;
437 QPrintDialog *PrintDialog = new QPrintDialog(Printer, this);
438 if (PrintDialog->exec() == QDialog::Accepted) {
439 QPainter Painter(Printer);
440 QPixmap Pixmap = QPixmap::grabWidget(Graph);
441 Painter.drawPixmap(0, 0, Pixmap);
442 }
443 delete Printer; delete PrintDialog;
444}
445
446void ddd::MenuSaveASCII() {
447 QString Filename = QFileDialog::getSaveFileName(this,
448 "Filename", ".", "Text files (*.txt *.ascii *.asc);;All files (*)");
449 if (Filename.length()>0) {
450 QFile File(Filename);
451 if (File.open(QFile::WriteOnly | QIODevice::Text | QFile::Truncate)) {
452 QTextStream Stream(&File);
453 for (int i=0; i<Signal->dataSize(); i++) {
454 Stream << Signal->y(i) << endl;
455 }
456 }
457 else {
458 QMessageBox::warning(this, "ddd Message","Could not write data to file.",QMessageBox::Ok);
459 }
460 }
461}
462
463void ddd::MenuHelp() {
464 QMessageBox Message;
465 Message.setText("The DRS Data Display program can be used for two purposes:\n\n"
466 "1. Reading and displaying the content of a raw data file written by the drsdaq program\n"
467 "2. Acquiring and displaying online data from a running drsdaq program via the socket interface\n\n"
468 "With an established socket connection, displaying of raw data files is disabled.\n\n"
469 "Navigation in signal display: Left mouse button to zoom, right to unzoom fully, middle to unzoom one level. Left mouse button plus shift key to pan.\n"
470 "When unzoomed, the axes are rescaled automatically.");
471
472 Message.setWindowTitle("ddd Help");
473 Message.exec();
474}
475
476void ddd::MenuAbout() {
477 QMessageBox::about(this, "ddd About","DRS Data Display\n\n"
478 "Written by Oliver Grimm, IPP, ETH Zurich\n"
479 "Event display by Quirin Weitzel.\n\n"
480 "This version compiled "__DATE__".\n\n"
481 "Graphical user interface implemented with Qt.\n"
482 "Bug reports to oliver.grimm@phys.ethz.ch.");
483}
Note: See TracBrowser for help on using the repository browser.