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;
}

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_ = new PileupMixingContent(numInteractionList, bunchCrossingList, TrueInteractionList, eventInfoList, bunchSpacing);
  } 

  PileupMixingContent* getEventPileupInfo() override { return PileupInfo_; } 

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

  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_; 


  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;
  std::vector<CLHEP::HepRandomEngine*> randomEngines_;
  std::vector<std::pair<int,std::bitset<6>>> theAffectedAPVvector;

  PileupMixingContent* PileupInfo_;

};

#endif