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;
}