DataFormats/HcalRecHit/interface/HFPreRecHit.h - master - githubmirror/cmssw

DataFormats/HcalRecHit/interface/HFPreRecHit.h
90 строк3 KB
Matti Kortelainen
Move stored types in io_v1 namespace, reset class versions to 3
22 май 2026, 18:18
22 май 2026, 18:18909f324
#ifndef DATAFORMATS_HCALRECHIT_HFPRERECHIT_H
#define DATAFORMATS_HCALRECHIT_HFPRERECHIT_H

#include <utility>

#include "DataFormats/HcalRecHit/interface/HFQIE10Info.h"

/** \class HFPreRecHit
*
* Class to contain the info needed to perform HF reconstruction
* using QIE10 chips and dual-anode readout. Groups the information
* provided by a single PMT.
*/
namespace io_v1 {
  class HFPreRecHit {
  public:
    typedef HcalDetId key_type;

    constexpr HFPreRecHit() : hasInfo_{false, false} {}

    // In the constructor below, either "first" or "second"
    // pointer can be nullptr
    constexpr HFPreRecHit(const HcalDetId& id, const HFQIE10Info* first, const HFQIE10Info* second)
        : id_(id), hasInfo_{false, false} {
      if (first) {
        hfQIE10Info_[0] = *first;
        hasInfo_[0] = true;
      }
      if (second) {
        hfQIE10Info_[1] = *second;
        hasInfo_[1] = true;
      }
    }

    constexpr HcalDetId id() const { return id_; }

    // Get a pointer to the QIE10 info. nullptr will be returned
    // if the info with the given index does not exist or if the
    // index is out of range.
    constexpr HFQIE10Info const* getHFQIE10Info(unsigned index) const {
      if (index < 2 && hasInfo_[index])
        return &hfQIE10Info_[index];
      else
        return nullptr;
    }

    // Quantities simply added from both anodes
    constexpr float charge() const {
      float q = 0.f;
      for (unsigned i = 0; i < 2; ++i)
        if (hasInfo_[i])
          q += hfQIE10Info_[i].charge();
      return q;
    }
    constexpr float energy() const {
      float e = 0.f;
      for (unsigned i = 0; i < 2; ++i)
        if (hasInfo_[i])
          e += hfQIE10Info_[i].energy();
      return e;
    }

    // The following function returns a pair.
    // The first element of the pair is the charge asymmetry,
    // if calculated. The second element of the pair indicates
    // whether the asymmetry is actually calculated. "true"
    // means yes, it is. The asymmetry is calculated if all
    // of the following conditions are satisfied:
    //
    //   1) The data is available for both PMT anodes.
    //   2) The sum of the charges is positive.
    //   3) The sum of the charges is at least "chargeThreshold".
    //
    // If the asymmetry is not calculated, the first element of
    // the pair is set to 0 and the second to "false".
    //
    std::pair<float, bool> chargeAsymmetry(float chargeThreshold) const;

    // Similar function for the energy asymmetry
    std::pair<float, bool> energyAsymmetry(float energyThreshold) const;

  private:
    HcalDetId id_;

    HFQIE10Info hfQIE10Info_[2];
    bool hasInfo_[2];
  };
}  // namespace io_v1
using HFPreRecHit = io_v1::HFPreRecHit;
#endif  // DATAFORMATS_HCALRECHIT_HFPRERECHIT_H