HGCalTBRecHitProducer.cc

Go to the documentation of this file.

#include "HGCal/Reco/plugins/HGCalTBRecHitProducer.h"

HGCalTBRecHitProducer::HGCalTBRecHitProducer(const edm::ParameterSet& cfg)
    : outputCollectionName(cfg.getParameter<std::string>("OutputCollectionName")),
      _digisToken(consumes<HGCalTBDigiCollection>(cfg.getParameter<edm::InputTag>("digiCollection"))),
      _pedestalLow_filename(cfg.getParameter<std::string>("pedestalLow")),
      _pedestalHigh_filename(cfg.getParameter<std::string>("pedestalHigh")),
      _gainsLow_filename(cfg.getParameter<std::string>("gainLow")),
      _gainsHigh_filename(cfg.getParameter<std::string>("gainHigh")),
      _adcSaturation(cfg.getParameter<unsigned int>("adcSaturation"))
{
    produces <HGCalTBRecHitCollection>(outputCollectionName);

}

void HGCalTBRecHitProducer::produce(edm::Event& event, const edm::EventSetup& iSetup)
{

    std::auto_ptr<HGCalTBRecHitCollection> rechits(new HGCalTBRecHitCollection); //auto_ptr are automatically deleted when out of scope

    edm::Handle<HGCalTBDigiCollection> digisHandle;
    event.getByToken(_digisToken, digisHandle);

#ifdef ESPRODUCER
// the conditions should be defined by ESProduers and available in the iSetup
    edm::ESHandle<HGCalCondPedestals> pedestalsHandle;
    //iSetup.get<HGCalCondPedestals>().get(pedestalsHandle); ///\todo should we support such a syntax?
    HGCalCondPedestals pedestals = iSetup.get<HGCalCondPedestals>();

    edm::ESHandle<HGCalCondGains>     adcToGeVHandle;
    //iSetup.get<HGCalCondGains>().get(adcToGeVHandle);
    HGCalCondGains adcToGeV = iSetup.get<HGCalCondGains>();
#else
    HGCalCondGains adcToGeV_low, adcToGeV_high;
    HGCalCondPedestals pedestals_low, pedestals_high;
    HGCalCondObjectTextIO condIO(HGCalTBNumberingScheme::scheme());
    //HGCalElectronicsMap emap;
//    assert(io.load("mapfile.txt",emap)); ///\todo to be trasformed into exception
    if(_pedestalLow_filename != "") assert(condIO.load(_pedestalLow_filename, pedestals_low));
    if(_pedestalHigh_filename != "")    assert(condIO.load(_pedestalHigh_filename, pedestals_high));
    if(_gainsLow_filename != "")    assert(condIO.load(_gainsLow_filename, adcToGeV_low));
    if(_gainsHigh_filename != "")   assert(condIO.load(_gainsHigh_filename, adcToGeV_high));
    ///\todo check if reading the conditions from file some channels are not in the file!
#endif

    for(auto digi_itr = digisHandle->begin(); digi_itr != digisHandle->end(); ++digi_itr) {
#ifdef DEBUG
        std::cout << "[RECHIT PRODUCER: digi]" << *digi_itr << std::endl;
#endif

//------------------------------ this part should go into HGCalTBRecoAlgo class

        SKIROC2DataFrame digi(*digi_itr);
        unsigned int nSamples = digi.samples();

        // if there are more than 1 sample, we need to define a reconstruction algorithm
        // now taking the first sample
        for(unsigned int iSample = 0; iSample < nSamples; ++iSample) {

            float pedestal_low_value = (pedestals_low.size() > 0) ? pedestals_low.get(digi.detid())->value : 0;
            float pedestal_high_value = (pedestals_high.size() > 0) ? pedestals_high.get(digi.detid())->value : 0;
            float adcToGeV_low_value = (adcToGeV_low.size() > 0) ? adcToGeV_low.get(digi.detid())->value : 1;
            float adcToGeV_high_value = (adcToGeV_high.size() > 0) ? adcToGeV_high.get(digi.detid())->value : 1;

            float energyLow = digi[iSample].adcLow() - pedestal_low_value * adcToGeV_low_value;
            float energyHigh = digi[iSample].adcHigh() - pedestal_high_value * adcToGeV_high_value;

            HGCalTBRecHit recHit(digi.detid(), energyLow, energyHigh, digi[iSample].tdc()); ///\todo use time calibration!

            if(digi[iSample].adcHigh() > _adcSaturation) recHit.setFlag(kHighGainSaturated);
            if(digi[iSample].adcLow() > _adcSaturation) recHit.setFlag(kLowGainSaturated);
            

#ifdef DEBUG
            std::cout << recHit << std::endl;
#endif
            if(iSample == 0) rechits->push_back(recHit); ///\todo define an algorithm for the energy if more than 1 sample, code inefficient
        }
    }
    event.put(rechits, outputCollectionName);
}
// Should there be a destructor ??
DEFINE_FWK_MODULE(HGCalTBRecHitProducer);