HGCalRecHitSimpleAlgo.h

Go to the documentation of this file.

#ifndef RecoLocalCalo_HGCalRecAlgos_HGCalRecHitSimpleAlgo_HH
#define RecoLocalCalo_HGCalRecAlgos_HGCalRecHitSimpleAlgo_HH

#include "FWCore/Utilities/interface/Exception.h"
#include "RecoLocalCalo/HGCalRecAlgos/interface/HGCalRecHitAbsAlgo.h"
#include "DataFormats/HGCDigi/interface/HGCDataFrame.h"
#include "DataFormats/ForwardDetId/interface/HGCalDetId.h"
#include "DataFormats/ForwardDetId/interface/HGCScintillatorDetId.h"
#include "DataFormats/ForwardDetId/interface/HGCSiliconDetId.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include <iostream>

class HGCalRecHitSimpleAlgo : public HGCalRecHitAbsAlgo {
 public:
  // default ctor
  HGCalRecHitSimpleAlgo() {
    adcToGeVConstant_ = -1;
    adcToGeVConstantIsSet_ = false;
  }

  void setLayerWeights(const std::vector<float>& weights) override {
    weights_ = weights;
  }
  
  void setADCToGeVConstant(const float value) override {
    adcToGeVConstant_ = value;
    adcToGeVConstantIsSet_ = true;
  }
  
  
  // destructor
  ~HGCalRecHitSimpleAlgo() override { };
  
  HGCRecHit makeRecHit(const HGCUncalibratedRecHit& uncalibRH,
                               const uint32_t& flags = 0) const override {
    
    if(!adcToGeVConstantIsSet_) {
      throw cms::Exception("HGCalRecHitSimpleAlgoBadConfig") 
        << "makeRecHit: adcToGeVConstant_ not set before calling this method!";
    }
    
    DetId baseid = uncalibRH.id();
    unsigned layer = 0;
    if( DetId::Hcal == baseid.det() && HcalEndcap == baseid.subdetId() ) {
      layer =  HcalDetId(baseid).depth() + 40;
    } else if ( DetId::Forward == baseid.det() && HGCEE == baseid.subdetId() ) {
      layer = HGCalDetId(baseid).layer();
    } else if ( DetId::Forward == baseid.det() && HGCHEF == baseid.subdetId() ) {
      layer = HGCalDetId(baseid).layer() + 28;
    } else if ( DetId::HGCalEE == baseid.det() ) {
      layer = HGCSiliconDetId(baseid).layer();
    } else if ( DetId::HGCalHSi == baseid.det() ) {
      layer = HGCSiliconDetId(baseid).layer() + 28;
    } else if ( DetId::HGCalHSc == baseid.det() ) {
      layer = HGCScintillatorDetId(baseid).layer() + 28;
    } else {
      throw cms::Exception("InvalidRecHit")
    << "HGCalRecHitSimpleAlgo encountered a non-HGCal det id: " << baseid.det() << ' ' << baseid.subdetId() << ' ' << baseid.rawId();
    }

    //    float clockToNsConstant = 25;    
    float energy = uncalibRH.amplitude() * weights_[layer] * 0.001f;
    float time   = uncalibRH.jitter();

    //if(time<0) time   = 0; // fast-track digi conversion
    
    HGCRecHit rh( uncalibRH.id(), energy, time );
    
    // Now fill flags
    // all rechits from the digitizer are "good" at present
    rh.setFlag(HGCRecHit::kGood);
    
    return rh;
  }

private:
  float adcToGeVConstant_;
  bool  adcToGeVConstantIsSet_;
  std::vector<float> weights_;
};
#endif