DTDigitizer.h

Go to the documentation of this file.

#ifndef SimMuon_DTDigitizer_h
#define SimMuon_DTDigitizer_h

#include "FWCore/Framework/interface/EDProducer.h"

#include "DataFormats/DTDigi/interface/DTDigiCollection.h"
#include "SimDataFormats/DigiSimLinks/interface/DTDigiSimLinkCollection.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"

#include "DataFormats/GeometryVector/interface/LocalVector.h"

#include <vector>

namespace CLHEP {
  class RandGaussQ;
  class RandFlat;
}

class DTLayer;
class PSimHit;
class DTWireType;
class DTBaseDigiSync;
class DTTopology;
class DTDigiSyncBase;


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

class DTDigitizer : public edm::EDProducer {
  
 public:

  explicit DTDigitizer(const edm::ParameterSet&);
  ~DTDigitizer();
  virtual void produce(edm::Event&, const edm::EventSetup&);
  
 private:
  typedef std::pair<const PSimHit*,float> hitAndT; // hit & corresponding time
  typedef std::vector<hitAndT> TDContainer; // hits & times for one wire

  typedef std::map<DTWireId, std::vector<const PSimHit*> > DTWireIdMap; 
  typedef DTWireIdMap::iterator DTWireIdMapIter;  
  typedef DTWireIdMap::const_iterator DTWireIdMapConstIter;  

  // Sort hits container by time.
  struct hitLessT {
    bool operator()(const  hitAndT & h1, const hitAndT & h2) {
      if (h1.second < h2.second) return true;
      return false;
    }
  };

  // Calculate the drift time for one hit. 
  // if status flag == false, hit has to be discarded.
  std::pair<float,bool> computeTime(const DTLayer* layer,const DTWireId &wireId, 
                    const PSimHit *hit, 
                    const LocalVector &BLoc); //FIXME?? 
  
  // Calculate the drift time using the GARFIELD cell parametrization,
  // taking care of all conversions from CMSSW local coordinates
  // to the conventions used for the parametrization.
  std::pair<float,bool> driftTimeFromParametrization(float x, float alpha, float By,
                             float Bz) const;
  
  // Calculate the drift time for the cases where it is not possible
  // to use the GARFIELD cell parametrization.
  std::pair<float,bool> driftTimeFromTimeMap() const;
  
  // Add all delays other than drift times (signal propagation along the wire, 
  // TOF etc.; subtract calibration time.
  float externalDelays(const DTLayer* layer,
               const DTWireId &wireId, 
               const PSimHit *hit) const;

  // Store digis for one wire, taking into account the dead time.
  //FiXME put alias for the map.
  void storeDigis(DTWireId &wireId, 
          TDContainer &hits,
          DTDigiCollection &output, DTDigiSimLinkCollection &outputLinks);
  
  // Debug output
  void dumpHit(const PSimHit * hit, float xEntry, float xExit, const DTTopology &topo);
  
  // Double half-gaussian smearing.
  float asymGausSmear(double mean, double sigmaLeft, double sigmaRight) const;
  
  // Allow debugging and testing.
  friend class DTDigitizerAnalysis;

  //Its Atributes:
  double vPropWire;
  float deadTime;
  float smearing;
  bool debug;
  bool interpolate;
  bool onlyMuHits;

  std::string syncName;
  DTDigiSyncBase *theSync;

  std::string geometryType;

  // Ideal model. Used for debug
  bool IdealModel;
  float theConstVDrift;  

  // the random generator
  CLHEP::RandGaussQ* theGaussianDistribution;
  CLHEP::RandFlat* theFlatDistribution;

  // to configure the creation of Digi-Sim links
  bool MultipleLinks;
  float LinksTimeWindow;
  
  //Name of Collection use for create the XF 
  std::string collection_for_XF;

};
#endif