/data/refman/pasoursint/CMSSW_5_3_10_patch2/src/SimMuon/Neutron/interface/SubsystemNeutronWriter.h

Go to the documentation of this file.

#ifndef SubsystemNeutronWriter_h
#define SubsystemNeutronWriter_h

#include <vector>
#include <map>
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "CLHEP/Random/RandFlat.h"


class NeutronWriter;

class SubsystemNeutronWriter : public edm::EDProducer
{
public:

  explicit SubsystemNeutronWriter(edm::ParameterSet const& pset);

  virtual ~SubsystemNeutronWriter();

  void printStats();

  virtual void produce(edm::Event & e, edm::EventSetup const& c);

  virtual int localDetId(int globalDetId) const = 0;

  virtual int chamberType(int globalDetId) const = 0;

  virtual int chamberId(int globalDetId) const = 0;

  virtual bool accept(const edm::PSimHitContainer & cluster) const = 0;

  void initialize(int chamberType);

protected:


  virtual void writeHits(int chamberType, edm::PSimHitContainer & chamberHits);

  void writeCluster(int chamberType, const edm::PSimHitContainer & cluster);

  void adjust(PSimHit & h, float timeOffset, float smearing);

  void updateCount(int chamberType);

private:
  NeutronWriter * theHitWriter;
  CLHEP::RandFlat * theRandFlat;
  edm::InputTag theInputTag;
  double theNeutronTimeCut;
  double theTimeWindow;
  double theT0;
  int theNEvents;
  bool initialized;
  // true means to translate DetId into just layer number, e.g., 1-6 in CSC
  bool useLocalDetId_;
  std::map<int, int> theCountPerChamberType;
};

#endif