HcalSD.h

#ifndef SIMCORE_HCALSD_H
#define SIMCORE_HCALSD_H
 
/*~~~~~~~~~~~~~~*/
/*   DetDescr   */
/*~~~~~~~~~~~~~~*/
#include <algorithm>
#include <string>
#include <vector>
 
#include "DetDescr/HcalID.h"
#include "DetDescr/PackedIndex.h"
#include "SimCore/Event/SimCalorimeterHit.h"
#include "SimCore/G4User/TrackingAction.h"
#include "SimCore/SensitiveDetector.h"
#include "SimCore/TrackMap.h"
 
namespace simcore {
 
class HcalSD : public SensitiveDetector {
 public:
  static const std::string COLLECTION_NAME;
 
  HcalSD(const std::string& name, simcore::ConditionsInterface& ci,
         const framework::config::Parameters& params);
 
  virtual ~HcalSD() = default;
 
  bool isSensDet(G4LogicalVolume* volume) const override {
    auto region = volume->GetRegion();
    if (region and region->GetName().contains("CalorimeterRegion")) {
      const auto name{volume->GetName()};
      return std::find_if(std::begin(gdmlIdentifiers_),
                          std::end(gdmlIdentifiers_),
                          [&name](const auto& identifier) {
                            return name.contains(identifier);
                          }) != std::end(gdmlIdentifiers_);
    }
    return false;
  }
 
  ldmx::HcalID decodeCopyNumber(const std::uint32_t copyNumber,
                                const G4ThreeVector& localPosition,
                                const G4Box* scint);
 
  virtual G4bool ProcessHits(G4Step* aStep,
                             G4TouchableHistory* ROhist) override;
 
  virtual void saveHits(framework::Event& event) override {
    event.add(COLLECTION_NAME, hits_);
  }
 
  virtual void OnFinishedEvent() override { hits_.clear(); }
 
 private:
  // A list of identifiers used to find out whether or not a given logical
  // volume is one of the Hcal sensitive detector volumes. Any volume that is
  // part of the CalorimeterRegion region and has a name which contains at least
  // one of the identifiers in here will be considered a sensitive detector in
  // the Hcal.
  std::vector<std::string> gdmlIdentifiers_;
  // TODO: document!
  double birksc1_;
 
  // TODO: document!
  double birksc2_;
 
  // collection of hits to write to event bus
  std::vector<ldmx::SimCalorimeterHit> hits_;
 
};  // HcalSD
 
}  // namespace simcore
 
#endif  // SIMCORE_HCALSD_H