ME0DigiModel.h

Go to the documentation of this file.

#ifndef SimMuon_GEMDigitizer_ME0DigiModel_h
#define SimMuon_GEMDigitizer_ME0DigiModel_h

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "DataFormats/GEMDigi/interface/ME0DigiCollection.h"
#include "DataFormats/Common/interface/DetSet.h"
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "SimDataFormats/TrackerDigiSimLink/interface/StripDigiSimLink.h"
#include "SimDataFormats/GEMDigiSimLink/interface/ME0DigiSimLink.h"
#include "SimDataFormats/CrossingFrame/interface/MixCollection.h"

#include <map>
#include <set>

namespace CLHEP {
  class HepRandomEngine;
}

class ME0EtaPartition;
class ME0Geometry;
class PSimHit;

class ME0DigiModel {
public:
  typedef edm::DetSet<StripDigiSimLink> StripDigiSimLinks;
  typedef edm::DetSet<ME0DigiSimLink> ME0DigiSimLinks;

  virtual ~ME0DigiModel() {}

  void setGeometry(const ME0Geometry* geom) { geometry_ = geom; }

  const ME0Geometry* getGeometry() { return geometry_; }

  virtual void simulateSignal(const ME0EtaPartition*, const edm::PSimHitContainer&, CLHEP::HepRandomEngine* engine) = 0;

  virtual void simulateNoise(const ME0EtaPartition*, CLHEP::HepRandomEngine* engine) = 0;

  virtual std::vector<std::pair<int, int> > simulateClustering(const ME0EtaPartition*,
                                                               const PSimHit*,
                                                               const int,
                                                               CLHEP::HepRandomEngine* engine) = 0;

  void fillDigis(int rollDetId, ME0DigiCollection&);

  virtual void setup() = 0;

  const StripDigiSimLinks& stripDigiSimLinks() const { return stripDigiSimLinks_; }
  const ME0DigiSimLinks& me0DigiSimLinks() const { return theME0DigiSimLinks_; }

protected:
  ME0DigiModel(const edm::ParameterSet&) {}

  const ME0Geometry* geometry_;

  std::set<std::pair<int, int> > strips_;

  void addLinks(unsigned int strip, int bx);
  void addLinksWithPartId(unsigned int strip, int bx);

  // keeps track of which hits contribute to which channels
  typedef std::multimap<std::pair<unsigned int, int>, const PSimHit*, std::less<std::pair<unsigned int, int> > >
      DetectorHitMap;

  DetectorHitMap detectorHitMap_;
  StripDigiSimLinks stripDigiSimLinks_;
  ME0DigiSimLinks theME0DigiSimLinks_;
};
#endif