MTDRecHitAlgo.cc

Go to the documentation of this file.

#include "RecoLocalFastTime/FTLCommonAlgos/interface/MTDRecHitAlgoBase.h"

#include "RecoLocalFastTime/Records/interface/MTDTimeCalibRecord.h"
#include "RecoLocalFastTime/FTLCommonAlgos/interface/MTDTimeCalib.h"

class MTDRecHitAlgo : public MTDRecHitAlgoBase {
public:
  MTDRecHitAlgo(const edm::ParameterSet& conf, edm::ConsumesCollector& sumes);

  ~MTDRecHitAlgo() override {}

  void getEvent(const edm::Event&) final {}
  void getEventSetup(const edm::EventSetup&) final;

  FTLRecHit makeRecHit(const FTLUncalibratedRecHit& uRecHit, uint32_t& flags) const final;

private:
  double thresholdToKeep_, calibration_;
  const MTDTimeCalib* time_calib_;
  edm::ESGetToken<MTDTimeCalib, MTDTimeCalibRecord> tcToken_;
};

MTDRecHitAlgo::MTDRecHitAlgo(const edm::ParameterSet& conf, edm::ConsumesCollector& sumes)
    : MTDRecHitAlgoBase(conf, sumes),
      thresholdToKeep_(conf.getParameter<double>("thresholdToKeep")),
      calibration_(conf.getParameter<double>("calibrationConstant")) {
  tcToken_ = sumes.esConsumes<MTDTimeCalib, MTDTimeCalibRecord>(edm::ESInputTag("", "MTDTimeCalib"));
}

void MTDRecHitAlgo::getEventSetup(const edm::EventSetup& es) {
  auto pTC = es.getHandle(tcToken_);
  time_calib_ = pTC.product();
}

FTLRecHit MTDRecHitAlgo::makeRecHit(const FTLUncalibratedRecHit& uRecHit, uint32_t& flags) const {
  unsigned char flagsWord = uRecHit.flags();
  float timeError = uRecHit.timeError();

  float energy = 0.;
  float time = 0.;

  float position = -1.f;
  float positionError = -1.f;

  switch (flagsWord) {
    // BTL bar geometry with only the right SiPM information available
    case 0x2: {
      energy = uRecHit.amplitude().second;
      time = uRecHit.time().second;

      break;
    }
    // BTL bar geometry with left and right SiPMs information available
    case 0x3: {
      energy = 0.5 * (uRecHit.amplitude().first + uRecHit.amplitude().second);
      time = 0.5 * (uRecHit.time().first + uRecHit.time().second);

      position = uRecHit.position();
      positionError = uRecHit.positionError();

      break;
    }
    // ETL, BTL tile geometry, BTL bar geometry with only the left SiPM information available
    default: {
      energy = uRecHit.amplitude().first;  //for ETL, it is the time_over_threshold
      time = uRecHit.time().first;

      break;
    }
  }

  // --- Energy calibration: for the time being this is just a conversion pC --> MeV
  energy *= calibration_;

  // --- Time calibration: for the time being just removes a time offset in BTL
  time += time_calib_->getTimeCalib(uRecHit.id());

  FTLRecHit rh(uRecHit.id(), uRecHit.row(), uRecHit.column(), energy, time, timeError, position, positionError);

  // Now fill flags
  // all rechits from the digitizer are "good" at present
  if (energy > thresholdToKeep_) {
    flags = FTLRecHit::kGood;
    rh.setFlag(flags);
  } else {
    flags = FTLRecHit::kKilled;
    rh.setFlag(flags);
  }

  return rh;
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_EDM_PLUGIN(MTDRecHitAlgoFactory, MTDRecHitAlgo, "MTDRecHitAlgo");