/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Emc/EmcRawEvent/EmcRawEvent-00-02-06/src/EmcDigi.cxx

Go to the documentation of this file.

#include "EmcRawEvent/EmcDigi.h"
#include "EmcWaveform/EmcWaveform.h"
#include "RawEvent/RawDataUtil.h"

EmcDigi::EmcDigi(const Identifier& id, const unsigned int timeChannel, const unsigned int chargeChannel) :
  RawData(id, timeChannel, chargeChannel),
  m_measure(0) {
}

EmcDigi::EmcDigi(const Identifier& id, const unsigned int timeChannel, const unsigned int chargeChannel,
    const unsigned int measure) :
  RawData(id, timeChannel, chargeChannel),
  m_measure(measure) {
}

EmcDigi::EmcDigi(const unsigned int id) :
  RawData(id), 
  m_measure(0) {
}

//mixing two digi with EmcWaveform
EmcDigi& EmcDigi::operator+=(const EmcDigi &digi)
{
  if(m_id != digi.identify()) // do nothing
    return *this;
  
  double e1=RawDataUtil::EmcCharge(m_measure,m_chargeChannel);
  double t1=RawDataUtil::EmcTime(m_timeChannel);
  double e2=RawDataUtil::EmcCharge(digi.getMeasure(),digi.getChargeChannel());
  double t2=RawDataUtil::EmcTime(digi.getTimeChannel());

  EmcWaveform wave1,wave2;
  wave1.makeWaveform(e1,t1);
  wave2.makeWaveform(e2,t2);

  //mixing
  wave1 += wave2;

  double energy;
  int time;
  energy = wave1.max(time);

  m_chargeChannel = RawDataUtil::EmcChargeChannel(energy);
  m_timeChannel = RawDataUtil::EmcTimeChannel(time);
  m_measure = RawDataUtil::EmcChargeMeasure(energy);

  return *this;
}

---