SiStripDigitizer.h

Go to the documentation of this file.

#ifndef SiStripDigitizer_h
#define SiStripDigitizer_h

#include <map>
#include <memory>
#include <string>
#include <vector>
#include <bitset>
#include "SimGeneral/MixingModule/interface/DigiAccumulatorMixMod.h"
#include "SimDataFormats/PileupSummaryInfo/interface/PileupMixingContent.h"
#include "FWCore/Framework/interface/ESHandle.h"

class TrackerTopology;

namespace CLHEP {
  class HepRandomEngine;
}

namespace edm {
  class ConsumesCollector;
  class ProducerBase;
  class Event;
  class EventSetup;
  class ParameterSet;
  template <typename T>
  class Handle;
  class StreamID;
}  // namespace edm

class MagneticField;
class PileUpEventPrincipal;
class PSimHit;
class SiStripDigitizerAlgorithm;
class StripGeomDetUnit;
class TrackerGeometry;

class SiStripDigitizer : public DigiAccumulatorMixMod {
public:
  explicit SiStripDigitizer(const edm::ParameterSet& conf, edm::ProducerBase& mixMod, edm::ConsumesCollector& iC);

  ~SiStripDigitizer() override;

  void initializeEvent(edm::Event const& e, edm::EventSetup const& c) override;
  void accumulate(edm::Event const& e, edm::EventSetup const& c) override;
  void accumulate(PileUpEventPrincipal const& e, edm::EventSetup const& c, edm::StreamID const&) override;
  void finalizeEvent(edm::Event& e, edm::EventSetup const& c) override;

  void StorePileupInformation(std::vector<int>& numInteractionList,
                              std::vector<int>& bunchCrossingList,
                              std::vector<float>& TrueInteractionList,
                              std::vector<edm::EventID>& eventInfoList,
                              int bunchSpacing) override {
    PileupInfo_ = std::make_unique<PileupMixingContent>(
        numInteractionList, bunchCrossingList, TrueInteractionList, eventInfoList, bunchSpacing);
  }

  PileupMixingContent* getEventPileupInfo() override { return PileupInfo_.get(); }

private:
  void accumulateStripHits(edm::Handle<std::vector<PSimHit>>,
                           const TrackerTopology* tTopo,
                           size_t globalSimHitIndex,
                           const unsigned int tofBin);

  typedef std::vector<std::string> vstring;
  typedef std::map<unsigned int, std::vector<std::pair<const PSimHit*, int>>, std::less<unsigned int>> simhit_map;
  typedef simhit_map::iterator simhit_map_iterator;

  const std::string gainLabel;
  const std::string hitsProducer;
  const vstring trackerContainers;
  const std::string ZSDigi;
  const std::string SCDigi;
  const std::string VRDigi;
  const std::string PRDigi;
  const std::string geometryType;
  const bool useConfFromDB;
  const bool zeroSuppression;
  const bool makeDigiSimLinks_;
  const bool includeAPVSimulation_;
  const double fracOfEventsToSimAPV_;


  std::map<std::string, size_t> crossingSimHitIndexOffset_;

  std::unique_ptr<SiStripDigitizerAlgorithm> theDigiAlgo;
  std::map<uint32_t, std::vector<int>> theDetIdList;
  edm::ESHandle<TrackerGeometry> pDD;
  edm::ESHandle<MagneticField> pSetup;
  std::map<unsigned int, StripGeomDetUnit const*> detectorUnits;
  CLHEP::HepRandomEngine* randomEngine_ = nullptr;
  std::vector<std::pair<int, std::bitset<6>>> theAffectedAPVvector;

  std::unique_ptr<PileupMixingContent> PileupInfo_;
};

#endif