MuonSimHitProducer.h

Go to the documentation of this file.

#ifndef FastSimulation_MuonSimHitProducer_MuonSimHitProducer_h
#define FastSimulation_MuonSimHitProducer_MuonSimHitProducer_h

//
// Package:    MuonSimHitProducer
// Class:      MuonSimHitProducer
// 
//
//  Author:  Martijn Mulders
// Created:  Wed July 11 12:37:24 CET 2007
// $Id: MuonSimHitProducer.h,v 1.10 2010/05/13 15:23:21 aperrott Exp $
//


// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"

// FastSimulation headers
class MagneticField;
class TrackerGeometry;
class DTGeometry;
class CSCGeometry;
class RPCGeometry;
class MuonServiceProxy;
class MaterialEffects;
class TrajectoryStateOnSurface;
class Propagator;
class RandomEngineAndDistribution;
class SimTrack;
class SimVertex;
/*
namespace reco { 
  class Muon;
}
*/

namespace edm { 
  class ParameterSet;
  class Event;
  class EventSetup;
}

//
// class declaration
//

class MuonSimHitProducer : public edm::stream::EDProducer <> {
   public:

      explicit MuonSimHitProducer(const edm::ParameterSet&);
      ~MuonSimHitProducer() override;

   private:

      MuonServiceProxy *theService;
      Chi2MeasurementEstimator theEstimator;

      const MagneticField*  magfield;
      const DTGeometry*     dtGeom;
      const CSCGeometry*    cscGeom;
      const RPCGeometry*    rpcGeom;
      const Propagator*     propagatorWithMaterial;
            Propagator* propagatorWithoutMaterial;

      MaterialEffects* theMaterialEffects;
  
      void beginRun(edm::Run const& run, const edm::EventSetup & es) override;
      void produce(edm::Event&, const edm::EventSetup&) override;
      void readParameters(const edm::ParameterSet&, 
              const edm::ParameterSet&,
              const edm::ParameterSet& );

      // Parameters to emulate the muonSimHit association inefficiency due to delta's
      double kDT;
      double fDT;
      double kCSC;
      double fCSC;

      void applyMaterialEffects(TrajectoryStateOnSurface& tsosWithdEdx,
                TrajectoryStateOnSurface& tsos,
                double radPath,
                                RandomEngineAndDistribution const*,
                                HepPDT::ParticleDataTable const&);

          
      // ----------- parameters ---------------------------- 
      bool fullPattern_;
      bool doL1_ , doL3_ , doGL_;

      // tags
      edm::InputTag simMuonLabel;
      edm::InputTag simVertexLabel;

      // tokens
      edm::EDGetTokenT<std::vector<SimTrack> > simMuonToken;
      edm::EDGetTokenT<std::vector<SimVertex> > simVertexToken;
      
      
};

#endif