#include <bm12controller.h>

Inheritance diagram for BM12Controller:

Public Member Functions
	BM12Controller (MCAController ctrl, QwtPlot graph, QObject *parent=0)

virtual	~BM12Controller ()

EpicsMCA *	mca ()

EpicsA2DConverter *	a2d ()

EpicsPulseAnalyser *	pulseAnalyser (int chan)

EpicsPowerSupply *	powerSupply ()

double	calibratedValue (int bin)

void	selectSCA (int ch, int win)

void	save (int ch, int mem)

void	clear (int ch, int mem)

void	count (int ch, int win)

void	updateGraph ()

void	saveData (int ch)

void	saveEnergyWindows (int ch)

void	loadEnergyWindows (int ch)

Private Slots
void	mcaValueChanged ()

Private Attributes
int	m_NChannels

EpicsMCA	m_MCA

EpicsA2DConverter	m_ADC

EpicsPowerSupply	m_HVPS

QList< EpicsPulseAnalyser * >	m_TCA

QList< BM12Waveforms * >	m_Waveforms

MCAController *	m_Controller

QwtPlot *	m_Graph

int	m_CurrentChannel

int	m_CurrentWindow

QString	m_SaveDataPath

QString	m_SaveWindowsPath

Detailed Description

Definition at line 15 of file bm12controller.h.

Constructor & Destructor Documentation

BM12Controller::BM12Controller	(	MCAController *	ctrl,
		QwtPlot *	graph,
		QObject *	parent = `0`
	)

Definition at line 11 of file bm12controller.cpp.

References m_ADC, m_HVPS, m_MCA, m_NChannels, m_TCA, m_Waveforms, mcaValueChanged(), EpicsA2DConverter::setGain(), EpicsPowerSupply::setLimit(), EpicsA2DConverter::setLowerLevelDiscrim(), EpicsA2DConverter::setOffset(), EpicsA2DConverter::setRange(), EpicsA2DConverter::setUpperLevelDiscrim(), and EpicsPowerSupply::setVoltage().

   : QObject(parent),
     m_NChannels(13),
     m_MCA("bm12det:aim_adc1"),
     m_ADC("bm12det:adc1"),
     m_HVPS("bm12det:hvps1"),
     m_Controller(ctrl),
     m_Graph(graph),
     m_CurrentChannel(0),
     m_CurrentWindow(0),
     m_SaveDataPath(QDir::currentPath()),
     m_SaveWindowsPath(QDir::currentPath())
 {
 //   printf("BM12Controller::BM12Controller\n");
 
   for (int i=0; i<m_NChannels; i++) {
     m_TCA.append(new EpicsPulseAnalyser(QString("bm12det:tca%1").arg(i+1)));
     m_Waveforms.append(new BM12Waveforms());
   }
 
   connect(&m_MCA, SIGNAL(valueChanged()), this, SLOT(mcaValueChanged()));
 
   m_ADC.setGain(2048);
   m_ADC.setRange(2048);
   m_ADC.setOffset(0);
   m_ADC.setLowerLevelDiscrim(5);
   m_ADC.setUpperLevelDiscrim(100);
 
   m_HVPS.setVoltage(500);
   m_HVPS.setLimit(500);
 }

BM12Controller::~BM12Controller ( )

virtual

Definition at line 43 of file bm12controller.cpp.

References m_TCA, and m_Waveforms.

 {
   while (!m_TCA.isEmpty()) {
     delete m_TCA.takeFirst();
   }
 
   while (!m_Waveforms.isEmpty()) {
     delete m_Waveforms.takeFirst();
   }
 }

Member Function Documentation

EpicsA2DConverter * BM12Controller::a2d ( )

Definition at line 59 of file bm12controller.cpp.

References m_ADC.

 {
   return &m_ADC;
 }

double BM12Controller::calibratedValue ( int bin )

Definition at line 215 of file bm12controller.cpp.

Referenced by saveData(), and updateGraph().

 {
   return bin*0.0538287-0.060349;
 }

void BM12Controller::clear	(	int	ch,
		int	mem
	)

Definition at line 108 of file bm12controller.cpp.

References EpicsMCA::erase(), m_MCA, m_Waveforms, and updateGraph().

Referenced by BM12ControlPanel::clear1(), BM12ControlPanel::clear2(), BM12ControlPanel::clear3(), and BM12ControlPanel::clearTotal().

 {
   if (ch >= 0 && ch < 13) {
     BM12Waveforms *wf = m_Waveforms[ch];
 
     if (wf) {
       switch (win) {
       case 0:
         m_MCA.erase();
         wf -> YData().clear();
         break;
 
       case 1:
         wf -> YDataSaved1().clear();
         break;
 
       case 2:
         wf -> YDataSaved2().clear();
         break;
 
       case 3:
         wf -> YDataSaved3().clear();
         break;
       }
     }
 
     updateGraph();
   }
 }

void BM12Controller::count	(	int	ch,
		int	win
	)

Definition at line 168 of file bm12controller.cpp.

References EpicsMCA::eraseAndStart(), m_CurrentChannel, m_CurrentWindow, m_MCA, and mcaValueChanged().

Referenced by BM12ControlPanel::countMCA().

 {
 //   printf("BM12Controller::count(%d,%d)\n", ch,win);
 
   m_CurrentChannel = ch;
   m_CurrentWindow  = win;
 
   m_MCA.eraseAndStart();
 
   mcaValueChanged();
 }

void BM12Controller::loadEnergyWindows ( int ch )

Definition at line 379 of file bm12controller.cpp.

References m_SaveWindowsPath, pulseAnalyser(), EpicsPulseAnalyser::setLowerLevel(), and EpicsPulseAnalyser::setUpperLevel().

Referenced by BM12ControlPanel::loadEnergyWindows().

 {
   QString name = QFileDialog::getOpenFileName(NULL, "Load Energy Windows...", m_SaveWindowsPath);
 
   if (name.length()) {
     m_SaveWindowsPath = name;
     
 //     printf("load windows from file %s\n", qPrintable(name));
 
     FILE *ifile = fopen(qPrintable(name),"r");
     QTextStream istr(ifile);
     QString line;
     QRegExp match("epics_put\\(\\s*\"bm12det:tca(\\d+)SCA(\\d+)(\\S+)\",\\s*(\\S+)\\s*\\)\\;");
 
     do {
       line = istr.readLine();
 
       if (match.exactMatch(line)) {
         int chan = match.cap(1).toInt()-1;
         int win  = match.cap(2).toInt()-1;
         int islow = match.cap(3)=="_LOW";
         int ishigh = match.cap(3)=="_HI";
 
         double val = match.cap(4).toDouble();
                         
         if (islow) {
 //        printf("set chan %d, win %d _LOW = %g\n", chan, win, val);
           pulseAnalyser(chan)->setLowerLevel(win,val);
         } else if (ishigh) {
 //        printf("set chan %d, win %d _HI = %g\n", chan, win, val);
           pulseAnalyser(chan)->setUpperLevel(win,val);
         }
       }
 
     } while (!line.isNull());
   }
 }

EpicsMCA * BM12Controller::mca ( )

Definition at line 54 of file bm12controller.cpp.

References m_MCA.

Referenced by BM12ControlPanel::BM12ControlPanel().

 {
   return &m_MCA;
 }

void BM12Controller::mcaValueChanged ( )

privateslot

Definition at line 180 of file bm12controller.cpp.

References m_CurrentChannel, m_CurrentWindow, m_MCA, m_Waveforms, updateGraph(), and EpicsArray::value().

Referenced by BM12Controller(), and count().

 {
 //   printf("BM12Controller::mcaValueChanged() %d, %d\n", m_CurrentChannel, m_CurrentWindow);
 
   if (m_CurrentChannel >= 0 && m_CurrentChannel < 13) {
     BM12Waveforms *wf = m_Waveforms[m_CurrentChannel];
 
     if (wf) {
       switch (m_CurrentWindow) {
         
       case 0:
         wf -> YData() = m_MCA.value();
         break;
 
       case 1:
         wf -> YDataSCA1() = m_MCA.value();
         break;
 
       case 2:
         wf -> YDataSCA2() = m_MCA.value();
         break;
 
       case 3:
         wf -> YDataSCA3() = m_MCA.value();
         break;
 
       default:
         break;
       }
     }
   }
 
   updateGraph();
 }

EpicsPowerSupply * BM12Controller::powerSupply ( )

Definition at line 69 of file bm12controller.cpp.

References m_HVPS.

Referenced by BM12ControlPanel::BM12ControlPanel().

 {
   return &m_HVPS;
 }

EpicsPulseAnalyser * BM12Controller::pulseAnalyser ( int chan )

Definition at line 64 of file bm12controller.cpp.

References m_TCA.

Referenced by BM12ControlPanel::BM12ControlPanel(), loadEnergyWindows(), and selectSCA().

 {
   return m_TCA.value(chan);
 }

void BM12Controller::save	(	int	ch,
		int	mem
	)

Definition at line 138 of file bm12controller.cpp.

References m_Waveforms, and updateGraph().

Referenced by BM12ControlPanel::save1(), BM12ControlPanel::save2(), and BM12ControlPanel::save3().

 {
  //  printf("BM12Controller::save(%d,%d)\n", ch,win);
 
 //   m_CurrentChannel = ch;
 //   m_CurrentWindow  = win;
 
   if (ch >= 0 && ch < 13) {
     BM12Waveforms *wf = m_Waveforms[ch];
 
     if (wf) {
       switch (win) {
       case 1:
         wf -> YDataSaved1() = wf -> YData();
         break;
 
       case 2:
         wf -> YDataSaved2() = wf -> YData();
         break;
 
       case 3:
         wf -> YDataSaved3() = wf -> YData();
         break;
       }
     }
 
     updateGraph();
   }
 }

void BM12Controller::saveData ( int ch )

Definition at line 316 of file bm12controller.cpp.

References calibratedValue(), m_SaveDataPath, m_Waveforms, max(), BM12Waveforms::YData(), BM12Waveforms::YDataSaved1(), BM12Waveforms::YDataSaved2(), BM12Waveforms::YDataSaved3(), BM12Waveforms::YDataSCA1(), BM12Waveforms::YDataSCA2(), and BM12Waveforms::YDataSCA3().

Referenced by BM12ControlPanel::saveData().

 {
   int maxsize = 0;
   int ncolumns = 0;
 
   BM12Waveforms *wf;
 
   foreach(wf, m_Waveforms) {
     max(ncolumns, maxsize, wf->YData().size());
     max(ncolumns, maxsize, wf->YDataSCA1().size());
     max(ncolumns, maxsize, wf->YDataSCA2().size());
     max(ncolumns, maxsize, wf->YDataSCA3().size());
     max(ncolumns, maxsize, wf->YDataSaved1().size());
     max(ncolumns, maxsize, wf->YDataSaved2().size());
     max(ncolumns, maxsize, wf->YDataSaved3().size());
   }
 
   QString name = QFileDialog::getSaveFileName(NULL, "Save Data...", m_SaveDataPath);
 
   if (name.length()) {
     m_SaveDataPath = name;
     
 //     printf("save mca data in file %s\n", qPrintable(name));
 //     printf("dirName %s\n", qPrintable(m_SaveDataPath));
 
     FILE *ofile = fopen(qPrintable(name),"a+");
 
     fprintf(ofile, "#S 1\n");
     fprintf(ofile, "#N %d\n", ncolumns+2);
     fprintf(ofile, "#L  point  percent");
 
     for (int i=0; i<m_Waveforms.size(); i++) {
       wf = m_Waveforms[i];
       if (wf->YData().size())      fprintf(ofile,"  ch%d", i+1);
       if (wf->YDataSCA1().size())  fprintf(ofile,"  ch%dsca1", i+1);
       if (wf->YDataSCA2().size())  fprintf(ofile,"  ch%dsca2", i+1);
       if (wf->YDataSCA3().size())  fprintf(ofile,"  ch%dsca3", i+1);
       if (wf->YDataSaved1().size()) fprintf(ofile, "  ch%dsav1", i+1);
       if (wf->YDataSaved2().size()) fprintf(ofile, "  ch%dsav2", i+1);
       if (wf->YDataSaved3().size()) fprintf(ofile, "  ch%dsav3", i+1);
     }
     fprintf(ofile, "\n");
 
     for (int row=0; row<maxsize; row++) {
       fprintf(ofile,"%d\t%g", row, calibratedValue(row));
 
       for (int i=0; i<m_Waveforms.size(); i++) {
         wf = m_Waveforms[i];
         if (wf->YData().size()>row)      fprintf(ofile,"\t%g", wf->YData()[row]);
         if (wf->YDataSCA1().size()>row)  fprintf(ofile,"\t%g", wf->YDataSCA1()[row]);
         if (wf->YDataSCA2().size()>row)  fprintf(ofile,"\t%g", wf->YDataSCA2()[row]);
         if (wf->YDataSCA3().size()>row)  fprintf(ofile,"\t%g", wf->YDataSCA3()[row]);
         if (wf->YDataSaved1().size()>row) fprintf(ofile, "\t%g", wf->YDataSaved1()[row]);
         if (wf->YDataSaved2().size()>row) fprintf(ofile, "\t%g", wf->YDataSaved2()[row]);
         if (wf->YDataSaved3().size()>row) fprintf(ofile, "\t%g", wf->YDataSaved3()[row]);
       }
       fprintf(ofile, "\n");
     }
 
     fclose(ofile);
   }
 }

void BM12Controller::saveEnergyWindows ( int ch )

Definition at line 417 of file bm12controller.cpp.

References EpicsPulseAnalyser::channelCount(), EpicsPulseAnalyser::lowerLevel(), m_SaveWindowsPath, m_TCA, EpicsObject::pvName(), EpicsPulseAnalyser::upperLevel(), and EpicsDouble::value().

Referenced by BM12ControlPanel::saveEnergyWindows().

 {
   QString cmd;
   EpicsPulseAnalyser *tca;
 
   foreach(tca, m_TCA) {
     for (int i=0; i<tca->channelCount(); i++) {
       EpicsDouble *lll = tca->lowerLevel(i);
       EpicsDouble *hll = tca->upperLevel(i);
 
       cmd += QString("epics_put( \"%1\", %2 );\n").arg(lll->pvName()).arg(lll->value());
       cmd += QString("epics_put( \"%1\", %2 );\n").arg(hll->pvName()).arg(hll->value());
     }
   }
   
 //   printf("Energy Windows\n%s\n", qPrintable(cmd));
 
   QString name = QFileDialog::getSaveFileName(NULL, "Save Energy Windows...", m_SaveWindowsPath);
 
   if (name.length()) {
     m_SaveWindowsPath = name;
 
 //     printf("save windows in file %s\n", qPrintable(name));
 
     FILE *ofile = fopen(qPrintable(name), "w");
     fputs(qPrintable(cmd), ofile);
     fclose(ofile);
   }
 }

void BM12Controller::selectSCA	(	int	ch,
		int	win
	)

Definition at line 74 of file bm12controller.cpp.

References EpicsPulseAnalyser::channelCount(), m_ADC, m_CurrentChannel, m_CurrentWindow, m_NChannels, pulseAnalyser(), EpicsA2DConverter::setCoincidence(), EpicsPulseAnalyser::setGated(), EpicsPulseAnalyser::setSelected(), EpicsObject::settle(), and updateGraph().

Referenced by BM12ControlPanel::selectSCA(), and BM12ControlPanel::setChannel().

 {
 //   printf("BM12Controller::selectSCA(%d,%d)\n", ch,win);
 
 //   pulseAnalyser(m_CurrentChannel)->setSelected(false);
 //   pulseAnalyser(m_CurrentChannel)->setGated(m_CurrentWindow, false);
 
   m_CurrentChannel = ch;
   m_CurrentWindow  = win;
 
 //   pulseAnalyser(m_CurrentChannel)->setSelected(true);
 //   pulseAnalyser(m_CurrentChannel)->setGated(m_CurrentWindow, true);
 
   for (int i=0; i<m_NChannels; i++) {
     EpicsPulseAnalyser *p = pulseAnalyser(i);
 
     p->setSelected(i==ch);
 
     for (int j=0; j<p->channelCount(); j++) {
       p->setGated(j, j==(win-1));
     }
   }
 
   if (win == 0) {
     m_ADC.setCoincidence(1);
   } else {
     m_ADC.setCoincidence(0);
   }
 
   pulseAnalyser(0)->settle(0.25);
 
   updateGraph();
 }

void BM12Controller::updateGraph ( )

Definition at line 220 of file bm12controller.cpp.

References calibratedValue(), m_Controller, m_CurrentChannel, m_Graph, m_Waveforms, BM12Waveforms::XData(), BM12Waveforms::YData(), BM12Waveforms::YDataSaved1(), BM12Waveforms::YDataSaved2(), BM12Waveforms::YDataSaved3(), BM12Waveforms::YDataSCA1(), BM12Waveforms::YDataSCA2(), and BM12Waveforms::YDataSCA3().

Referenced by clear(), mcaValueChanged(), save(), and selectSCA().

 {
   if (m_CurrentChannel >= 0 && m_CurrentChannel < 13) {
     BM12Waveforms *wf = m_Waveforms[m_CurrentChannel];
 
     if (wf) {
       int xlen = wf->XData().size();
       int ylen = wf->YData().size();
 
 //       if (ylen != xlen) {
       wf->XData().resize(ylen);
       for (int i=0; i<ylen; i++) {
         wf->XData()[i] = calibratedValue(i);
       }
 //       }
 
       m_Graph -> detachItems(QwtPlotItem::Rtti_PlotCurve);
       m_Graph -> detachItems(QwtPlotItem::Rtti_PlotMarker);
 
       QwtPlotCurve *pc = new QwtPlotCurve("Total");
       pc -> setSamples(wf->XData().data(), wf->YData().data(), ylen);
       QPen pen(Qt::black);
       pc -> setPen(pen);
       pc -> setSymbol(new QwtSymbol(QwtSymbol::Ellipse,QBrush(),pen,QSize(3,3)));
       pc -> attach(m_Graph);
       
       if (ylen == wf->YDataSCA1().size()) {
         QwtPlotCurve *pc1 = new QwtPlotCurve("SCA1");
         QPen pen(Qt::red);
         pc1 -> setSamples(wf->XData().data(), wf->YDataSCA1().data(), ylen);
         pc1 -> setPen(pen);
         pc1 -> setSymbol(new QwtSymbol(QwtSymbol::Ellipse,QBrush(),pen,QSize(3,3)));
         pc1 -> setBrush(QBrush(QColor(255,0,0,20)));
         pc1 -> attach(m_Graph);
       }
 
       if (ylen == wf->YDataSCA2().size()) {
         QwtPlotCurve *pc1 = new QwtPlotCurve("SCA2");
         QPen pen(Qt::green);
         pc1 -> setSamples(wf->XData().data(), wf->YDataSCA2().data(), ylen);
         pc1 -> setPen(pen);
         pc1 -> setSymbol(new QwtSymbol(QwtSymbol::Ellipse,QBrush(),pen,QSize(3,3)));
         pc1 -> setBrush(QBrush(QColor(0,255,0,20)));
         pc1 -> attach(m_Graph);
       }
 
       if (ylen == wf->YDataSCA3().size()) {
         QwtPlotCurve *pc1 = new QwtPlotCurve("SCA3");
         QPen pen(Qt::blue);
         pc1 -> setSamples(wf->XData().data(), wf->YDataSCA3().data(), ylen);
         pc1 -> setPen(pen);
         pc1 -> setSymbol(new QwtSymbol(QwtSymbol::Ellipse,QBrush(),pen,QSize(3,3)));
         pc1 -> setBrush(QBrush(QColor(0,0,255,20)));
         pc1 -> attach(m_Graph);
       }
 
       if (ylen == wf->YDataSaved1().size()) {
         QwtPlotCurve *pc1 = new QwtPlotCurve("Saved A");
         QPen pen(Qt::darkRed);
         pc1 -> setSamples(wf->XData().data(), wf->YDataSaved1().data(), ylen);
         pc1 -> setPen(pen);
         pc1 -> attach(m_Graph);
       }
 
       if (ylen == wf->YDataSaved2().size()) {
         QwtPlotCurve *pc1 = new QwtPlotCurve("Saved B");
         QPen pen(Qt::darkGreen);
         pc1 -> setSamples(wf->XData().data(), wf->YDataSaved2().data(), ylen);
         pc1 -> setPen(pen);
         pc1 -> attach(m_Graph);
       }
 
       if (ylen == wf->YDataSaved3().size()) {
         QwtPlotCurve *pc1 = new QwtPlotCurve("Saved C");
         QPen pen(Qt::darkBlue);
         pc1 -> setSamples(wf->XData().data(), wf->YDataSaved3().data(), ylen);
         pc1 -> setPen(pen);
         pc1 -> attach(m_Graph);
       }
 
       m_Controller -> replot();
     }
   }
 }