source: tools/ddd/Functions.cpp@ 60

Last change on this file since 60 was 60, checked in by qweitzel, 15 years ago
comments added for M0 display code
File size: 13.7 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 CloseDatafile(); // Close previous file if open
22
23 // Write run header to temporary file
24 FILE *Tmpfile = tmpfile();
25 if(Tmpfile==NULL) {
26 QMessageBox::warning(this, "ddd Message","Could not open temporary file.",QMessageBox::Ok);
27 CloseDatafile();
28 return;
29 }
30
31 switch (RD->OpenDataFile(FilenameBox->text().toAscii().data(), Tmpfile)) {
32 case CTX_FOPEN: QMessageBox::warning(this, "ddd Message","Could not open file.",QMessageBox::Ok);
33 return;
34 case CTX_RHEADER: QMessageBox::warning(this, "ddd Message","Could not read run header.",QMessageBox::Ok);
35 return;
36 case CTX_BSTRUCT: QMessageBox::warning(this, "ddd Message","Could not read board structures.",QMessageBox::Ok);
37 return;
38 default: break;
39 }
40 if (RD->RHeader->MagicNum == MAGICNUM_OPEN) {
41 QMessageBox::warning(this, "ddd Message","Magic number in run header indicates that the file has not been closed properly.",QMessageBox::Ok);
42 }
43 if (RD->RHeader->MagicNum == MAGICNUM_ERROR) {
44 QMessageBox::warning(this, "ddd Message","Magic number in run header indicates that an error occurred while writing the file.",QMessageBox::Ok);
45 }
46
47 rewind(Tmpfile);
48 QTextStream in(Tmpfile);
49 QString text = in.readAll();
50 RunHeaderDisplay->setPlainText(text);
51 fclose(Tmpfile);
52
53 // Enable spin boxes, set ranges and display first event
54 EventNo->setEnabled(true);
55 ChannelNo->setEnabled(true);
56 BoardNo->setEnabled(true);
57 PixelID->setEnabled(true);
58 EventNo->setValue(1);
59 EventNo->setRange(1, RD->RHeader->Events);
60 ChannelNo->setValue(0);
61 ChannelNo->setRange(0, RD->RHeader->NChannels*RD->RHeader->NChips-1);
62 BoardNo->setValue(0);
63 BoardNo->setRange(0, RD->RHeader->NBoards-1);
64 DisplayEvent();
65}
66
67// +++ Close data file file if open, delete event header and disable spin boxes and displays +++
68void ddd::CloseDatafile() {
69 if(RD->CloseDataFile()==CTX_OK) {
70 EventNo->setEnabled(false);
71 ChannelNo->setEnabled(false);
72 BoardNo->setEnabled(false);
73 PixelID->setEnabled(false);
74 RunHeaderDisplay->clear();
75 EventHeaderDisplay->clear();
76 Signal->hide();
77 }
78}
79
80// +++ Read event header and display event (only called if Datafile is open) +++
81void ddd::DisplayEvent(int) {
82
83 PixelID->setText(PixMap->DRS_to_Pixel(BoardNo->value(),ChannelNo->value()/10,ChannelNo->value()%10).c_str());
84 if(Socket->state() == QAbstractSocket::ConnectedState) return; // do not execute if socket is open
85
86 // Read event
87 FILE *Tmpfile = tmpfile();
88
89 switch(RD->ReadEvent(EventNo->value(), Tmpfile)) {
90 case !CTX_OK:
91 QMessageBox::warning(this, "ddd Warning","Could not read event.",QMessageBox::Ok);
92 EventHeaderDisplay->clear();
93 break;
94 default:
95 rewind(Tmpfile);
96 QTextStream in(Tmpfile);
97 QString text = in.readAll();
98 EventHeaderDisplay->setPlainText(text);
99 // Case data in double format required by qwt library
100 double* x = new double [RD->RHeader->Samples];
101 double* y = new double [RD->RHeader->Samples];
102
103 for (unsigned int i=0; i<RD->RHeader->Samples; i++) {
104 x[i] = (double) (i/RD->BStruct[BoardNo->value()].NomFreq);
105 y[i] = (double) RD->Data[BoardNo->value()*RD->RHeader->NChips*RD->RHeader->NChannels *
106 RD->RHeader->Samples+ChannelNo->value()*RD->RHeader->Samples+i]*RD->BStruct[BoardNo->value()].ScaleFactor;
107 }
108 Signal->setData(x, y, RD->RHeader->Samples);
109 Signal->show();
110 Zoomer->setZoomBase(Signal->boundingRect());
111
112 //Get data for M0 display (event based)
113
114 double z[6][6];//36 pixels
115
116 for(unsigned int i=0; i<RD->RHeader->NBoards; i++) {//board loop
117 for(unsigned int j=0; j<RD->RHeader->NChips; j++) {//chip loop
118 for(unsigned int k=0; k<RD->RHeader->NChannels; k++) {//channel loop
119
120 //only interested in M0 data
121 if( ( (i==0 || i==1) && (j<=1) && (k<=7) ) || ( (i==2) && (j==0) && (k<=3) ) ) {
122
123 //get module, superpixel and pixel number from pixel name
124
125 std::string pixelname = PixMap->DRS_to_Pixel(i,j,k);
126 char pixelname_copy[256];
127 memset(pixelname_copy,'\0',256);
128 pixelname.copy(pixelname_copy, 256);
129
130 char delim[] = "-";
131 char *buffer = NULL;
132 int module = -1;
133 int superpixel = -1;
134 int pixel = -1;
135
136 buffer = strtok(pixelname_copy, delim);
137 module = atoi(buffer);
138 buffer = strtok(NULL, delim);
139 superpixel = atoi(buffer);
140 buffer = strtok(NULL, delim);
141 pixel = atoi(buffer);
142
143 int binx = 5-(int((superpixel-1)/3)*2)-(int((pixel%4)/2));
144 int biny = 5-(((superpixel-1)%3)*2)-(int((pixel-1)/2));
145
146 //search maximum sample amplitude within following window -> should not be hardcoded!
147 int StartBin = 100;
148 int StopBin = 120;
149
150 for(int l=StartBin; l<(StopBin+1); l++){
151
152 float sample = RD->Data[i*RD->RHeader->NChips*RD->RHeader->NChannels*RD->RHeader->Samples+
153 j*RD->RHeader->NChannels*RD->RHeader->Samples+
154 k*RD->RHeader->Samples+
155 l]*RD->BStruct[i].ScaleFactor;
156
157 if (sample > z[binx][biny]) {
158 z[binx][biny]=sample;
159 }
160
161 }//sample loop
162
163 }//only M0 data
164
165 }//channel loop
166 }//chip loop
167 }//board loop
168
169 //fill data to M0 display (event based)
170 Signal2D->setData(SpectrogramDataM0(z));
171 Graph2D->axisWidget(QwtPlot::yRight)->setColorMap(Signal2D->data().range(),Signal2D->colorMap());
172 Graph2D->setAxisScale(QwtPlot::yRight,Signal2D->data().range().minValue(),Signal2D->data().range().maxValue() );
173 Graph2D->replot();
174 //Signal2D->show();
175
176 delete[] x; delete[] y;
177
178 }
179 if(Tmpfile!=NULL) fclose(Tmpfile);
180}
181
182// +++ Open sub window handling the socket interface +++
183void ddd::OpenSocketWindow() {
184
185 if(SocketWindow->isVisible()) SocketWindow->hide();
186 else SocketWindow->show();
187}
188
189// +++ Open sub window for M0 Display +++
190void ddd::OpenM0Window() {
191
192 if(M0Window->isVisible()) M0Window->hide();
193 else M0Window->show();
194}
195
196
197// +++ Acquire data through socket (acquire botton only available if socket exists) +++
198void ddd::GetSignalFromSocket() {
199 char Command[MAX_COM_SIZE];
200
201 GetButton->setEnabled(false);
202 WaitForData = true;
203 sprintf(Command, "read %d %d %d restart", BoardNo->value(), ChannelNo->value()/10, ChannelNo->value()%10);
204 Socket->write(Command);
205}
206
207// Quit application when clicking close button on window
208void ddd::closeEvent(QCloseEvent *) {
209 qApp->quit();
210}
211
212// +++ Connecting or disconnecting from client +++
213void ddd::MakeConnection() {
214
215 if(Socket->state() == QAbstractSocket::ConnectedState) {
216 ManualDisconnect = true;
217 Socket->disconnectFromHost();
218 }
219 else {
220 if (RD->IsFileOpen() && QMessageBox::question(this, "ddd Request","Connecting will close current data file. Proceed?",
221 QMessageBox::No, QMessageBox::Yes) != QMessageBox::Yes) return;
222 Socket->connectToHost(IPAddress->text(),Port->value());
223 Connect->setEnabled(false); // While waiting for connection, button not available
224 Socket->waitForConnected(SOCKET_TIMEOUT);
225 Connect->setEnabled(true);
226 if(Socket->state() != QAbstractSocket::ConnectedState)
227 QMessageBox::warning(this, "ddd Message","Could not connect to host.",QMessageBox::Ok);
228 else {
229 Connect->setText("Disconnect");
230 ConnectAction->setText("Disconnect");
231 Port->setEnabled(false);
232 IPAddress->setEnabled(false);
233 Command->setEnabled(true);
234 GetButton->setEnabled(true);
235 ManualDisconnect = false;
236
237 OpenAction->setEnabled(false);
238 FilenameBox->setEnabled(false);
239 LoadButton->setEnabled(false);
240 FilenameBox->clear();
241 CloseDatafile();
242
243 ChannelNo->setEnabled(true);
244 BoardNo->setEnabled(true);
245 PixelID->setEnabled(true);
246 ChannelNo->setRange(0, 65535);
247 BoardNo->setRange(0, 65535);
248
249 TabWidget->setTabEnabled(1,false);
250 TabWidget->setTabEnabled(2,false);
251
252 RunHeaderDisplay->clear();
253 EventHeaderDisplay->clear();
254 Signal->hide();
255 }
256 }
257}
258
259// +++ Send command to socket (command button available only if socket existing) +++
260void ddd::SendToSocket() {
261 Socket->write(Command->text().toAscii());
262 Command->clear();
263}
264
265// +++ Read data from socket and display +++
266void ddd::ReadFromSocket() {
267 // Check if channel data expected and error message arrived
268 QByteArray Data = Socket->readAll();
269 if (WaitForData && Data.contains("Error")) {
270 WaitForData = false;
271 GetButton->setEnabled(true);
272 QMessageBox::warning(this, "ddd Message","Could not read waveform data from socket.",QMessageBox::Ok);
273 return;
274 }
275
276 // Check if channel data were transmitted, if yes and waiting for data, extract and plot them
277 SocketOutput->insertPlainText(Data);
278 QString Text = SocketOutput->document()->toPlainText().trimmed();
279 if (WaitForData && Text.endsWith(QLatin1String("==END=="))) {
280 // Extract text between ==START== and ==END==
281 QByteArray Data=Text.mid(Text.lastIndexOf("==START==")+9,Text.length()-Text.lastIndexOf("==START==")-16).toAscii();
282
283 char *NextNumber = strtok(Data.data()," "); // Number of entries that follow
284 int Count=0, NumberOfEntries = atoi(NextNumber);
285 double *x = new double [NumberOfEntries];
286 double *y = new double [NumberOfEntries];
287
288 // Convert all entries (separated by a whitespace) to numbers
289 while((NextNumber=strtok(NULL, " "))!=NULL && Count<NumberOfEntries)
290 *(y+Count++) = atof(NextNumber);
291 if (Count==NumberOfEntries && NextNumber!=0)
292 QMessageBox::warning(this, "ddd Message","Found too many numbers in data block, truncated.",QMessageBox::Ok);
293 // Apply sampling frequency and scaling factor
294 for(int i=2; i<Count; i++) {
295 x[i] = (i-2) / y[0];
296 y[i] = y[i] * y[1];
297 }
298 if(NumberOfEntries>2) {
299 Signal->setData(x+2, y+2, NumberOfEntries-2); // Copies data, arrays can be deleted afterwards
300 Signal->show();
301 Zoomer->setZoomBase(Signal->boundingRect());
302 }
303 delete[] x; delete[] y;
304
305 if(ContinuousBox->isChecked()) {
306 usleep(100000); // Wait to limit maximum update rate
307 GetSignalFromSocket();
308 }
309 else {
310 WaitForData = false;
311 GetButton->setEnabled(true);
312 }
313 }
314}
315
316// +++ Disconnect from socket +++
317void ddd::GotDisconnected() {
318 Connect->setText("Connect");
319 ConnectAction->setText("Connect");
320 Port->setEnabled(true);
321 IPAddress->setEnabled(true);
322 Command->setEnabled(false);
323
324 GetButton->setEnabled(false);
325 FilenameBox->setEnabled(true);
326 OpenAction->setEnabled(true);
327 LoadButton->setEnabled(true);
328 ChannelNo->setEnabled(false);
329 BoardNo->setEnabled(false);
330 PixelID->setEnabled(false);
331 Signal->hide();
332 TabWidget->setTabEnabled(1, true);
333 TabWidget->setTabEnabled(2, true);
334
335 SocketOutput->clear();
336 if(!ManualDisconnect) QMessageBox::warning(this, "ddd Message","Socket disconnected, maybe host gone.",QMessageBox::Ok);
337}
338
339// +++ Translate pixel ID +++
340void ddd::TranslatePixelID() {
341
342 int Board = PixMap->Pixel_to_DRSboard(PixelID->text().toStdString());
343 int Chip = PixMap->Pixel_to_DRSchip(PixelID->text().toStdString());
344 int Channel = PixMap->Pixel_to_DRSchannel(PixelID->text().toStdString());
345
346 if(Board>=BoardNo->minimum() && Board<=BoardNo->maximum() &&
347 (Chip*10+Channel)>=ChannelNo->minimum() && (Chip*10+Channel)<=ChannelNo->maximum()) {
348 BoardNo->setValue(Board);
349 ChannelNo->setValue(Chip*10+Channel);
350 }
351 else if(Board==999999999) QMessageBox::warning(this, "ddd Message","Pixel ID unknown.",QMessageBox::Ok);
352 else QMessageBox::warning(this, "ddd Message","Pixel ID out of current range.",QMessageBox::Ok);
353}
354
355
356//------------------------------------------------------------------
357//**************************** All menus ***************************
358//------------------------------------------------------------------
359
360void ddd::MenuSave() {
361 QString Filename = QFileDialog::getSaveFileName(this,
362 "Filename of image", "/home/ogrimm/ddd", "Image files (*.bmp *.jpg *.png *.ppm *.tiff *.xbm *.xpm);;All files (*)");
363 if (Filename.length()>0) {
364 QPixmap Pixmap = QPixmap::grabWidget(Graph);
365 if(!Pixmap.save(Filename)) {
366 QMessageBox::warning(this, "ddd Message","Could not write image file.",QMessageBox::Ok);
367 remove(Filename.toAscii().data());
368 }
369 }
370}
371
372void ddd::MenuPrint() {
373 QPrinter *Printer = new QPrinter;
374 QPrintDialog *PrintDialog = new QPrintDialog(Printer, this);
375 if (PrintDialog->exec() == QDialog::Accepted) {
376 QPainter Painter(Printer);
377 QPixmap Pixmap = QPixmap::grabWidget(Graph);
378 Painter.drawPixmap(0, 0, Pixmap);
379 }
380 delete Printer; delete PrintDialog;
381}
382
383void ddd::MenuHelp() {
384 QMessageBox Message;
385 Message.setText("The DRS Data Display program can be used for two purposes:\n\n"
386 "1. Reading and displaying the content of a raw data file written by the drsdaq program\n"
387 "2. Acquiring and displaying online data from a running drsdaq program via the socket interface\n\n"
388 "With an established socket connection, displaying of raw data files is disabled.");
389 Message.setWindowTitle("ddd Help");
390 Message.exec();
391}
392
393void ddd::MenuAbout() {
394 QMessageBox::about(this, "ddd About","DRS Data Display\n\n"
395 "Written by Oliver Grimm, IPP, ETH Zurich\n"
396 "This version compiled "__DATE__".\n\n"
397 "Graphical user interface implemented with Qt.\n"
398 "Bug reports to oliver.grimm@phys.ethz.ch.");
399}
Note: See TracBrowser for help on using the repository browser.