HcalElectronicsSim.cc

Go to the documentation of this file.

#include "SimCalorimetry/HcalSimAlgos/interface/HcalElectronicsSim.h"
#include "SimCalorimetry/HcalSimAlgos/interface/HcalAmplifier.h"
#include "SimCalorimetry/HcalSimAlgos/interface/HcalCoderFactory.h"
#include "DataFormats/HcalDigi/interface/HBHEDataFrame.h"
#include "DataFormats/HcalDigi/interface/HODataFrame.h"
#include "DataFormats/HcalDigi/interface/HFDataFrame.h"
#include "DataFormats/HcalDigi/interface/ZDCDataFrame.h"
#include "DataFormats/HcalDigi/interface/HcalUpgradeDataFrame.h"
#include "DataFormats/HcalDigi/interface/QIE10DataFrame.h"
#include "CLHEP/Random/RandFlat.h"
#include <math.h>

HcalElectronicsSim::HcalElectronicsSim(HcalAmplifier * amplifier, const HcalCoderFactory * coderFactory, bool PreMixing)
  : theAmplifier(amplifier),
    theCoderFactory(coderFactory),
    theStartingCapId(0),
    theStartingCapIdIsRandom(true),
    PreMixDigis(PreMixing)
{
}


HcalElectronicsSim::~HcalElectronicsSim() {
}


void HcalElectronicsSim::setDbService(const HcalDbService * service) {
  //  theAmplifier->setDbService(service);
  theTDC.setDbService(service);
}

template<class Digi> 
void HcalElectronicsSim::convert(CaloSamples & frame, Digi & result, CLHEP::HepRandomEngine* engine) {
  result.setSize(frame.size());
  theAmplifier->amplify(frame, engine);
  theCoderFactory->coder(frame.id())->fC2adc(frame, result, theStartingCapId);
}

template<> 
void HcalElectronicsSim::convert<QIE10DataFrame>(CaloSamples & frame, QIE10DataFrame & result, CLHEP::HepRandomEngine* engine) {
  theAmplifier->amplify(frame, engine);
  theCoderFactory->coder(frame.id())->fC2adc(frame, result, theStartingCapId);
}

void HcalElectronicsSim::analogToDigital(CLHEP::HepRandomEngine* engine, CaloSamples & lf, HBHEDataFrame & result) {
  convert<HBHEDataFrame>(lf, result, engine);
  if(PreMixDigis) {
    for(int isample = 0; isample !=lf.size(); ++isample) {
      uint16_t theADC = round(10.0*lf[isample]);
      unsigned capId = result[isample].capid();

      if(theADC > 126) {
    uint16_t keepADC = result[isample].adc();
    HcalQIESample mysamp(keepADC, capId, 0, 0, true, true); // set error bit as a flag
    result.setSample(isample, HcalQIESample(keepADC, capId, 0, 0, true, true) );
      }
      else {
    result.setSample(isample, HcalQIESample(theADC, capId, 0, 0) ); // preserve fC, no noise
    HcalQIESample mysamp(theADC, capId, 0, 0);
      }
    }
  }
}


void HcalElectronicsSim::analogToDigital(CLHEP::HepRandomEngine* engine, CaloSamples & lf, HODataFrame & result) {
  convert<HODataFrame>(lf, result, engine);
  if(PreMixDigis) {
    for(int isample = 0; isample !=lf.size(); ++isample) {
      uint16_t theADC = round(10.0*lf[isample]);
      unsigned capId = result[isample].capid();

      if(theADC > 126) {
    uint16_t keepADC = result[isample].adc();
    HcalQIESample mysamp(keepADC, capId, 0, 0, true, true);// set error bit as a flag
    result.setSample(isample, HcalQIESample(keepADC, capId, 0, 0, true, true) );
      }
      else {
    result.setSample(isample, HcalQIESample(theADC, capId, 0, 0) ); // preserve fC, no noise
    HcalQIESample mysamp(theADC, capId, 0, 0);
      }
    }
  }
}


void HcalElectronicsSim::analogToDigital(CLHEP::HepRandomEngine* engine, CaloSamples & lf, HFDataFrame & result) {
  convert<HFDataFrame>(lf, result, engine);
  if(PreMixDigis) {
    for(int isample = 0; isample !=lf.size(); ++isample) {
      uint16_t theADC = round(10.0*lf[isample]);
      unsigned capId = result[isample].capid();

      if(theADC > 126) {
    uint16_t keepADC = result[isample].adc();
    HcalQIESample mysamp(keepADC, capId, 0, 0, true, true);// set error bit as a flag
    result.setSample(isample, HcalQIESample(keepADC, capId, 0, 0, true, true) );
      }
      else {
    result.setSample(isample, HcalQIESample(theADC, capId, 0, 0)  );  // preserve fC, no noise
    HcalQIESample mysamp(theADC, capId, 0, 0);
      }
    }
  }
}

void HcalElectronicsSim::analogToDigital(CLHEP::HepRandomEngine* engine, CaloSamples & lf, ZDCDataFrame & result) {
  convert<ZDCDataFrame>(lf, result, engine);
  if(PreMixDigis) {
    for(int isample = 0; isample !=lf.size(); ++isample) {
      uint16_t theADC = round(10.0*lf[isample]);
      unsigned capId = result[isample].capid();

      if(theADC > 126) {
    uint16_t keepADC = result[isample].adc();
    HcalQIESample mysamp(keepADC, capId, 0, 0, true, true);//set error bit as a flag
    result.setSample(isample, HcalQIESample(keepADC, capId, 0, 0, true, true) );
      }
      else {
    result.setSample(isample, HcalQIESample(theADC, capId, 0, 0)  );  // preserve fC, no noise
    HcalQIESample mysamp(theADC, capId, 0, 0);
      }
    }
  }
}


void HcalElectronicsSim::analogToDigital(CLHEP::HepRandomEngine* engine, CaloSamples & lf,
                     HcalUpgradeDataFrame & result) {
  convert<HcalUpgradeDataFrame>(lf, result, engine);
//   std::cout << HcalDetId(lf.id()) << ' ' << lf;
  theTDC.timing(lf, result, engine);
}

void HcalElectronicsSim::analogToDigital(CLHEP::HepRandomEngine* engine, CaloSamples & lf,
                     QIE10DataFrame & result) {
  convert<QIE10DataFrame>(lf, result, engine);
  //TODO:
  //HcalTDC extension for QIE10?
  //PreMixDigis for QIE10?
}

void HcalElectronicsSim::newEvent(CLHEP::HepRandomEngine* engine) {
  // pick a new starting Capacitor ID
  if(theStartingCapIdIsRandom)
  {
    theStartingCapId = CLHEP::RandFlat::shootInt(engine, 4);
    theAmplifier->setStartingCapId(theStartingCapId);
  }
}


void HcalElectronicsSim::setStartingCapId(int startingCapId)
{
  theStartingCapId = startingCapId;
  theAmplifier->setStartingCapId(theStartingCapId);
  // turns off random capIDs forever for this instance
  theStartingCapIdIsRandom = false;
}