DTDigitizer Class Reference

#include <DTDigitizer.h>

Inheritance diagram for DTDigitizer:

Classes
struct	hitLessT

Public Member Functions
	DTDigitizer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

Private Types
typedef std::map< DTWireId, std::vector< const PSimHit * > >	DTWireIdMap
typedef DTWireIdMap::const_iterator	DTWireIdMapConstIter
typedef DTWireIdMap::iterator	DTWireIdMapIter
typedef std::pair< const PSimHit *, float >	hitAndT
typedef std::vector< hitAndT >	TDContainer

Private Member Functions
float	asymGausSmear (double mean, double sigmaLeft, double sigmaRight, CLHEP::HepRandomEngine *) const
std::pair< float, bool >	computeTime (const DTLayer *layer, const DTWireId &wireId, const PSimHit *hit, const LocalVector &BLoc, CLHEP::HepRandomEngine *)
std::pair< float, bool >	driftTimeFromParametrization (float x, float alpha, float By, float Bz, CLHEP::HepRandomEngine *) const
std::pair< float, bool >	driftTimeFromTimeMap () const
void	dumpHit (const PSimHit *hit, float xEntry, float xExit, const DTTopology &topo)
float	externalDelays (const DTLayer *layer, const DTWireId &wireId, const PSimHit *hit) const
void	storeDigis (DTWireId &wireId, TDContainer &hits, DTDigiCollection &output, DTDigiSimLinkCollection &outputLinks)

Private Attributes
edm::EDGetTokenT< CrossingFrame< PSimHit > >	cf_token
std::string	collection_for_XF
float	deadTime
bool	debug
std::string	geometryType
bool	IdealModel
bool	interpolate
float	LinksTimeWindow
std::string	mix_
bool	MultipleLinks
bool	onlyMuHits
float	smearing
std::string	syncName
float	theConstVDrift
std::unique_ptr< DTDigiSyncBase >	theSync
double	vPropWire

Friends
class	DTDigitizerAnalysis

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Digitize the muon drift tubes. The parametrisation function in DTDriftTimeParametrization from P.G.Abia, J.Puerta is used in all cases where it is applicable.

Authors

: G. Bevilacqua, N. Amapane, G. Cerminara, R. Bellan

Definition at line 45 of file DTDigitizer.h.

Member Typedef Documentation

typedef std::map<DTWireId, std::vector<const PSimHit*> > DTDigitizer::DTWireIdMap

private

Definition at line 57 of file DTDigitizer.h.

typedef DTWireIdMap::const_iterator DTDigitizer::DTWireIdMapConstIter

private

Definition at line 59 of file DTDigitizer.h.

typedef DTWireIdMap::iterator DTDigitizer::DTWireIdMapIter

private

Definition at line 58 of file DTDigitizer.h.

typedef std::pair<const PSimHit*,float> DTDigitizer::hitAndT

private

Definition at line 54 of file DTDigitizer.h.

typedef std::vector<hitAndT> DTDigitizer::TDContainer

private

Definition at line 55 of file DTDigitizer.h.

Constructor & Destructor Documentation

DTDigitizer::DTDigitizer ( const edm::ParameterSet &  conf_ )

explicit

Definition at line 57 of file DTDigitizer.cc.

References beamerCreator::create(), muonDTDigis_cfi::deadTime, debug, Exception, hgcalDigitizer_cfi::geometryType, reco::get(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), muonDTDigis_cfi::IdealModel, HiCaloJetParameters_cff::interpolate, edm::Service< T >::isAvailable(), muonDTDigis_cfi::LinksTimeWindow, muonDTDigis_cfi::MultipleLinks, muonDTDigis_cfi::onlyMuHits, AlCaHLTBitMon_QueryRunRegistry::string, and DTLinearDriftAlgo_cfi::vPropWire.

                                                  {
  
  // Set verbose output
  debug=conf_.getUntrackedParameter<bool>("debug"); 
    
  if (debug) LogPrint("DTDigitizer")<<"Creating a DTDigitizer"<<endl;
  
  //register the Producer with a label
  //produces<DTDigiCollection>("MuonDTDigis"); // FIXME: Do I pass it by ParameterSet?
  produces<DTDigiCollection>(); // FIXME: Do I pass it by ParameterSet?
  //  produces<DTDigiSimLinkCollection>("MuonDTDigiSimLinks");
  produces<DTDigiSimLinkCollection>();
  
  //Parameters:

  // build digis only for mu hits (for debug purposes) 
  onlyMuHits=conf_.getParameter<bool>("onlyMuHits");
  
  // interpolate parametrization function
  interpolate=conf_.getParameter<bool>("interpolate");
  
  // Velocity of signal propagation along the wire (cm/ns)
  // For the default value
  // cfr. CMS-IN 2000-021:   (2.56+-0.17)x1e8 m/s
  //      CMS NOTE 2003-17:  (0.244)  m/ns
  vPropWire=conf_.getParameter<double>("vPropWire"); //24.4

  // Dead time for signals on the same wire (number from M. Pegoraro)  
  deadTime=conf_.getParameter<double>("deadTime"); //150
  
  // further configurable smearing
  smearing=conf_.getParameter<double>("Smearing"); // 3.

  // Sync Algo
  syncName = conf_.getParameter<string>("SyncName");
  theSync.reset( DTDigiSyncFactory::get()->create(syncName,conf_.getParameter<ParameterSet>("pset")) );

  // Debug flag to switch to the Ideal model
  // it uses a constant drift velocity and doesn't set any external delay
  IdealModel = conf_.getParameter<bool>("IdealModel");

  // Constant drift velocity needed by the above flag
  if(IdealModel)
    theConstVDrift = conf_.getParameter<double>("IdealModelConstantDriftVelocity"); // 55 um/ns
  else theConstVDrift = 55.;

  // get random engine
  edm::Service<edm::RandomNumberGenerator> rng;
  if ( ! rng.isAvailable()) {
    throw cms::Exception("Configuration")
      << "RandomNumberGeneratorService for DTDigitizer missing in cfg file";
  }

  // MultipleLinks=false ==> one-to-one correspondence between digis and SimHits 
  MultipleLinks = conf_.getParameter<bool>("MultipleLinks");
  // MultipleLinks=true ==> association of SimHits within a time window LinksTimeWindow 
  // (of the order of the resolution)
  LinksTimeWindow = conf_.getParameter<double>("LinksTimeWindow"); // (10 ns)

  //Name of Collection used for create the XF 
  mix_ = conf_.getParameter<std::string>("mixLabel");
  collection_for_XF = conf_.getParameter<std::string>("InputCollection");
  cf_token = consumes<CrossingFrame<PSimHit> >( edm::InputTag(mix_, collection_for_XF) );

  //String to choice between ideal (the deafult) and (mis)aligned geometry for the digitization step 
  geometryType = conf_.getParameter<std::string>("GeometryType");

}

Member Function Documentation

float DTDigitizer::asymGausSmear ( double  mean, double  sigmaLeft, double  sigmaRight, CLHEP::HepRandomEngine *  engine  ) const

private

Definition at line 442 of file DTDigitizer.cc.

References f, and lumiQTWidget::t.

                                                                     {

  double f = sigmaLeft/(sigmaLeft+sigmaRight);
  double t;

  if (CLHEP::RandFlat::shoot(engine) <= f) {
    t = CLHEP::RandGaussQ::shoot(engine, mean, sigmaLeft);
    t = mean - fabs(t - mean);
  } else {
    t = CLHEP::RandGaussQ::shoot(engine, mean, sigmaRight);
    t = mean + fabs(t - mean);
  }
  return static_cast<float>(t);
}

pair< float, bool > DTDigitizer::computeTime ( const DTLayer *  layer, const DTWireId &  wireId, const PSimHit *  hit, const LocalVector &  BLoc, CLHEP::HepRandomEngine *  engine  )

private

Definition at line 229 of file DTDigitizer.cc.

References Vector3DBase< T, FrameTag >::cross(), edmIntegrityCheck::d, debug, delta, dir, Vector3DBase< T, FrameTag >::dot(), PSimHit::entryPoint(), PSimHit::exitPoint(), muonDTDigis_cfi::IdealModel, PSimHit::localDirection(), PSimHit::localPosition(), M_PI, PV3DBase< T, PVType, FrameType >::mag(), PSimHit::momentumAtEntry(), DTTopology::none, DTTopology::onWhichBorder(), PSimHit::pabs(), PSimHit::particleType(), EnergyCorrector::pt, PVValHelper::pT, DTLayer::specificTopology(), mathSSE::sqrt(), theta(), Vector3DBase< T, FrameTag >::unit(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), DTTopology::xMax, DTTopology::xMin, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                       {

  LocalPoint entryP = hit->entryPoint();
  LocalPoint exitP = hit->exitPoint();
  int partType = hit->particleType();

  const DTTopology &topo = layer->specificTopology();

  // Pay attention: in CMSSW the rf of the SimHit is in the layer's rf
  
  if(debug)  LogPrint("DTDigitizer")<<"Hit local entry point: "<<entryP<<endl
         <<"Hit local exit point: "<<exitP<<endl;

  float xwire = topo.wirePosition(wireId.wire()); 
  float xEntry = entryP.x() - xwire;
  float xExit  = exitP.x() - xwire;

  if(debug) LogPrint("DTDigitizer")<<"wire position: "<<xwire
        <<" x entry in cell rf: "<<xEntry
        <<" x exit in cell rf: "<<xExit<<endl;
  
  DTTopology::Side entrySide = topo.onWhichBorder(xEntry,entryP.y(),entryP.z());
  DTTopology::Side exitSide  = topo.onWhichBorder(xExit,exitP.y(),exitP.z());

   if (debug) dumpHit(hit, xEntry, xExit,topo);

  // The bolean is used to flag the drift time computation
  pair<float,bool> driftTime(0.,false);  

  // if delta in gas->ignore, since it is included in the parametrisation.
  // FIXME: should check that it is actually a delta ray produced by a nearby
  // muon hit. 

  if (partType == 11 && entrySide == DTTopology::none) {
     if (debug) LogPrint("DTDigitizer") << "    e- hit in gas; discarding " << endl;
    return driftTime;
  }

  float By = BLoc.y();
  float Bz = BLoc.z();

  // Radius and sagitta according to direction of momentum
  // (just for printing)
  // NOTE: in cmsim, d is always taken // pHat!
  LocalVector d = (exitP-entryP);
  LocalVector pHat = hit->localDirection().unit();
  LocalVector hHat = (d.cross(pHat.cross(d))).unit();
  float cosAlpha = hHat.dot(pHat);
  float sinAlpha = sqrt(1.-cosAlpha*cosAlpha);
  float radius_P = (d.mag())/(2.*cosAlpha);
  float sagitta_P = radius_P*(1.-sinAlpha);

  // Radius, sagitta according to field bending
  // (just for printing)
  float halfd = d.mag()/2.;
  float BMag = BLoc.mag();
  LocalVector pT = (pHat - (BLoc.unit()*pHat.dot(BLoc.unit())))*(hit->pabs());
  float radius_B = (pT.mag()/(0.3*BMag))*100.;
  float sagitta_B;
  if (radius_B > halfd) {
    sagitta_B = radius_B - sqrt(radius_B*radius_B - halfd*halfd);
  } else {
    sagitta_B = radius_B;
  }

  // cos(delta), delta= angle between direction at entry and hit segment
  // (just for printing)
  float delta = pHat.dot(d.unit());
  if (debug) LogPrint("DTDigitizer") << "   delta                 = " << delta  << endl
          << "   cosAlpha              = " << cosAlpha << endl
          << "   sinAlpha              = " << sinAlpha << endl
          << "   pMag                  = " << pT.mag() << endl
          << "   bMag                  = " << BMag << endl
          << "   pT                    = " << pT << endl
          << "   halfd                 = " << halfd << endl
          << "   radius_P  (cm)        = " << radius_P << endl
          << "   sagitta_P (um)        = " << sagitta_P*10000. << endl
          << "   radius_B  (cm)        = " << radius_B << endl
          << "   sagitta_B (um)        = " << sagitta_B*10000. << endl;

  // Select cases where parametrization can not be used.
  bool noParametrisation = 
    ( ( entrySide == DTTopology::none || exitSide == DTTopology::none ) // case # 2,3,8,9 or 11
      || (entrySide == exitSide)                   // case # 4 or 10
      || ((entrySide == DTTopology::xMin && exitSide == DTTopology::xMax) || 
      (entrySide == DTTopology::xMax && exitSide == DTTopology::xMin)) // Hit is case # 7
      );

  // FIXME: now, debug warning only; consider treating those 
  // with TM algo. 
  if ( delta < 0.99996 // Track is not straight. FIXME: use sagitta?
       && (noParametrisation == false)) {
    if (debug) LogPrint("DTDigitizer") << "*** WARNING: hit is not straight, type = " << partType << endl;
  }

  //************ 5A ***************

  if (!noParametrisation) {
   
    LocalVector dir = hit->momentumAtEntry(); // ex Measurement3DVector dir = hit->measurementDirection(); //FIXME?
    float theta = atan(dir.x()/-dir.z())*180/M_PI;

    // FIXME: use dir if M.S. is included as GARFIELD option...
    //        otherwise use hit segment dirction.
    //    LocalVector dir0 = (exitP-entryP).unit();
    //    float theta = atan(dir0.x()/-dir0.z())*180/M_PI;
    float x;

    Local3DPoint pt = hit->localPosition();  //ex Measurement3DPoint pt = hit->measurementPosition(); // FIXME?
     
    if(fabs(pt.z()) < 0.002) { 
      // hit center within 20 um from z=0, no need to extrapolate.
      x = pt.x() - xwire;
    } else {
      x = xEntry - (entryP.z()*(xExit-xEntry))/(exitP.z()-entryP.z());
    }

    if(IdealModel) return make_pair(fabs(x)/theConstVDrift,true);
    else driftTime = driftTimeFromParametrization(x, theta, By, Bz, engine);

  }

 
  if ((driftTime.second)==false) {
    // Parametrisation not applicable, or failed. Use time map.
    driftTime = driftTimeFromTimeMap();
  }
  
  //************ 5B ***************

  // Signal propagation, TOF etc.
  if (driftTime.second) {
    driftTime.first += externalDelays(layer,wireId,hit);
  }
  return driftTime;
} 

pair< float, bool > DTDigitizer::driftTimeFromParametrization ( float  x, float  alpha, float  By, float  Bz, CLHEP::HepRandomEngine *  engine  ) const

private

Definition at line 370 of file DTDigitizer.cc.

References debug, GeomDetEnumerators::DT, f, RemoveAddSevLevel::flag, HiCaloJetParameters_cff::interpolate, edm::max(), DTDriftTimeParametrization::MB_DT_drift_time(), DTDriftTimeParametrization::drift_time::t_drift, DTDriftTimeParametrization::drift_time::t_width_m, DTDriftTimeParametrization::drift_time::t_width_p, theta(), and ntuplemaker::time.

                                                                                               {

  // Convert from CMSSW frame/units r.f. to parametrization ones.
  x *= 10.;  //cm -> mm 

  // FIXME: Current parametrisation can extrapolate above 21 mm,
  // however a detailed study is needed before using this.
   if (fabs(x) > 21.) {
     if (debug) LogPrint("DTDigitizer") << "*** WARNING: parametrisation: x out of range = "
             << x << ", skipping" << endl;
     return pair<float,bool>(0.f,false);
  }

  // Different r.f. of the parametrization:
  // X_par = X_ORCA; Y_par=Z_ORCA; Z_par = -Y_ORCA  

  float By_par = Bz;  // Bnorm
  float Bz_par = -By; // Bwire
  float theta_par = theta;

  // Parametrisation uses interpolation up to |theta|=45 deg,
  // |Bwire|=0.4, |Bnorm|=0.75; extrapolation above.
  if (fabs(theta_par)>45.) {
    if (debug) LogPrint("DTDigitizer") << "*** WARNING: extrapolating theta > 45: "
            << theta << endl;
    // theta_par = min(fabs(theta_par),45.f)*((theta_par<0.)?-1.:1.);
  }
  if (fabs(By_par)>0.75) {
    if (debug) LogPrint("DTDigitizer") << "*** WARNING: extrapolating Bnorm > 0.75: "
            << By_par << endl;
    // By_par = min(fabs(By_par),0.75f)*((By_par<0.)?-1.:1.);
  }
  if (fabs(Bz_par)>0.4) {
    if (debug) LogPrint("DTDigitizer") << "*** WARNING: extrapolating Bwire >0.4: "
            << Bz_par << endl;
    // Bz_par = min(fabs(Bz_par),0.4)*((Bz_par<0.)?-1.:1.);
  }

  DTDriftTimeParametrization::drift_time DT;
  static const DTDriftTimeParametrization par;
  unsigned short flag = par.MB_DT_drift_time (x, theta_par, By_par, Bz_par, 0, &DT, interpolate);

  if (debug) {
    LogPrint("DTDigitizer") << "    Parametrisation: x, theta, Bnorm, Bwire = "
     << x << " " <<  theta_par << " " << By_par << " " << Bz_par << endl
     << "  time=" << DT.t_drift
     << "  sigma_m=" <<  DT.t_width_m
     << "  sigma_p=" <<  DT.t_width_p << endl;
    if (flag!=1) {
      LogPrint("DTDigitizer") << "*** WARNING: call to parametrisation failed" << endl;
      return pair<float,bool>(0.f,false); 
    }
  }

  // Double half-gaussian smearing
  float time = asymGausSmear(DT.t_drift, DT.t_width_m, DT.t_width_p, engine);

  // Do not allow the smearing to lead to negative values
  time = max(time,0.f);

  // Apply a Gaussian smearing to account for electronic effects (cf. 2004 TB analysis)
  // The width of the Gaussian can be configured with the "Smearing" parameter

  double u = CLHEP::RandGaussQ::shoot(engine, 0., smearing);
  time += u;

  if (debug) LogPrint("DTDigitizer") << "  drift time = " << time << endl;
  
  return pair<float,bool>(time,true); 
}

pair< float, bool > DTDigitizer::driftTimeFromTimeMap ( ) const

private

Definition at line 458 of file DTDigitizer.cc.

References debug.

                                                         {
  // FIXME: not yet implemented.
  if (debug) LogPrint("DTDigitizer") << "   TimeMap " << endl;
  return pair<float,bool>(0.,false);
}

void DTDigitizer::dumpHit ( const PSimHit *  hit, float  xEntry, float  xExit, const DTTopology &  topo  )

private

Definition at line 573 of file DTDigitizer.cc.

References DTTopology::cellHeight(), DTTopology::cellLenght(), DTTopology::cellWidth(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), createfilelist::int, mag(), PSimHit::momentumAtEntry(), DTTopology::onWhichBorder(), PSimHit::pabs(), PSimHit::particleType(), PSimHit::processType(), PSimHit::trackId(), Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                  {
  
  LocalPoint entryP = hit->entryPoint();
  LocalPoint exitP = hit->exitPoint();
  
  DTTopology::Side entrySide = topo.onWhichBorder(xEntry,entryP.y(),entryP.z());
  DTTopology::Side exitSide  = topo.onWhichBorder(xExit,exitP.y(),exitP.z());
  //  ProcessTypeEnumerator pTypes;
  
  LogPrint("DTDigitizer") << endl
       << "------- SimHit: " << endl
       << "   Particle type         = " << hit->particleType() << endl
       << "   process type          = " << hit->processType() << endl
       << "   process type          = " << hit->processType() << endl
    // << "   packedTrackId         = " << hit->packedTrackId() << endl
       << "   trackId               = " << hit->trackId() << endl // new,is the same as the
                                                                  // previous?? FIXME-Check
       << "   |p|                   = " << hit->pabs() << endl
       << "   Energy loss           = " << hit->energyLoss() << endl
    // << "   timeOffset            = " << hit->timeOffset() << endl
    // << "   measurementPosition   = " << hit->measurementPosition() << endl
    // << "   measurementDirection  = " << hit->measurementDirection() << endl      //FIXME
       << "   localDirection        = " << hit->momentumAtEntry().unit() << endl    //FIXME is it a versor?
       << "   Entry point           = " << entryP << " cell x = " << xEntry << endl
       << "   Exit point            = " << exitP << " cell x = " << xExit << endl
       << "   DR =                  = " << (exitP-entryP).mag() << endl
       << "   Dx =                  = " << (exitP-entryP).x() << endl
       << "   Cell w,h,l            = (" << topo.cellWidth()
       << " , " << topo.cellHeight() 
       << " , " << topo.cellLenght() << ") cm" << endl
       << "   DY entry from edge    = " << topo.cellLenght()/2.-fabs(entryP.y())
       << "   DY exit  from edge    = " << topo.cellLenght()/2.-fabs(exitP.y())
       << "   entrySide = "  << (int)entrySide
       << " ; exitSide = " << (int)exitSide << endl;

}

float DTDigitizer::externalDelays ( const DTLayer *  layer, const DTWireId &  wireId, const PSimHit *  hit  ) const

private

Definition at line 466 of file DTDigitizer.cc.

References DTTopology::cellLenght(), debug, PSimHit::localPosition(), DTLayer::specificTopology(), PSimHit::tof(), DTLinearDriftAlgo_cfi::vPropWire, and PV3DBase< T, PVType, FrameType >::y().

                                            {
  
  // Time of signal propagation along wire.
  
  float wireCoord = hit->localPosition().y();
  float halfL     = (layer->specificTopology().cellLenght())/2.;
  float propgL = halfL - wireCoord; // the FE is always located at the pos coord.

  float propDelay = propgL/vPropWire;

  // Real TOF.
  float tof = hit->tof();  

  // Delays and t0 according to theSync

  double sync= theSync->digitizerOffset(&wireId,layer);


  if (debug) {
    LogPrint("DTDigitizer") << "    propDelay =" << propDelay
     << "; TOF=" << tof
     << "; sync= " << sync
     << endl;
  }
  
  return propDelay + tof + sync;
}

void DTDigitizer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 127 of file DTDigitizer.cc.

References debug, edm::EventID::event(), hgcalDigitizer_cfi::geometryType, edm::EventSetup::get(), edm::Event::getByToken(), edm::RandomNumberGenerator::getEngine(), edm::EventBase::id(), MagneticField::inTesla(), DTGeometry::layer(), DTWireId::layerId(), eostools::move(), convertSQLitetoXML_cfg::output, edm::Handle< T >::product(), edm::Event::put(), edm::EventID::run(), rpcPointValidation_cfi::simHit, trackerHits::simHits, edm::Event::streamID(), GeomDet::surface(), ntuplemaker::time, Surface::toGlobal(), and GloballyPositioned< T >::toLocal().

                                                                {

  edm::Service<edm::RandomNumberGenerator> rng;
  CLHEP::HepRandomEngine* engine = &rng->getEngine(iEvent.streamID());

  if(debug)
    LogPrint("DTDigitizer") << "--- Run: " << iEvent.id().run()
     << " Event: " << iEvent.id().event() << endl;
  
  //************ 1 ***************
   // create the container for the SimHits
  //  Handle<PSimHitContainer> simHits; 
  //  iEvent.getByLabel("g4SimHits","MuonDTHits",simHits);
    
  // use MixCollection instead of the previous
  Handle<CrossingFrame<PSimHit> > xFrame;
  iEvent.getByToken(cf_token, xFrame);
  
  unique_ptr<MixCollection<PSimHit> >
    simHits( new MixCollection<PSimHit>(xFrame.product()) );

   // create the pointer to the Digi container
  unique_ptr<DTDigiCollection> output(new DTDigiCollection());
   // pointer to the DigiSimLink container
  unique_ptr<DTDigiSimLinkCollection> outputLinks(new DTDigiSimLinkCollection());
  
  // Muon Geometry
  ESHandle<DTGeometry> muonGeom;
  iSetup.get<MuonGeometryRecord>().get(geometryType,muonGeom);

  // Magnetic Field  
  ESHandle<MagneticField> magnField;
  iSetup.get<IdealMagneticFieldRecord>().get(magnField);

  //************ 2 ***************

  // These are sorted by DetId, i.e. by layer and then by wire #
  //  map<DTDetId, vector<const PSimHit*> > wireMap;     
  DTWireIdMap wireMap;     
  
  for(MixCollection<PSimHit>::MixItr simHit = simHits->begin();
       simHit != simHits->end(); simHit++){
    
    // Create the id of the wire, the simHits in the DT known also the wireId
 
    DTWireId wireId( (*simHit).detUnitId() );
    // Fill the map
    wireMap[wireId].push_back(&(*simHit));
  }
  
  pair<float,bool> time(0.,false);

  //************ 3 ***************
  // Loop over the wires
  for(DTWireIdMapConstIter wire = wireMap.begin(); wire!=wireMap.end(); wire++){
    // SimHit Container associated to the wire
    const vector<const PSimHit*> & vhit = (*wire).second; 
    if(!vhit.empty()) {
      TDContainer tdCont; // It is a vector<pair<const PSimHit*,float> >;
      
      //************ 4 ***************
      DTWireId wireId = (*wire).first;

      //const DTLayer* layer = dynamic_cast< const DTLayer* > (muonGeom->idToDet(wireId.layerId()));
      const DTLayer* layer = muonGeom->layer(wireId.layerId()); 

      // Loop on the hits of this wire    
      for (vector<const PSimHit*>::const_iterator hit=vhit.begin();
       hit != vhit.end(); hit++){
    //************ 5 ***************
    LocalPoint locPos = (*hit)->localPosition();

    const LocalVector BLoc=layer->surface().toLocal(magnField->inTesla(layer->surface().toGlobal(locPos)));

    time = computeTime(layer, wireId, *hit, BLoc, engine);

    //************ 6 ***************
    if (time.second) {
      tdCont.push_back(make_pair((*hit),time.first));
    } else {
      if (debug) LogPrint("DTDigitizer") << "hit discarded" << endl;
    }
      }

      //************ 7 ***************

      // the loading must be done by layer but
      // the digitization must be done by wire (in order to take into account the dead time)

      storeDigis(wireId,tdCont,*output,*outputLinks);
    }
    
  }

  //************ 8 ***************  
  // Load the Digi Container in the Event
  //iEvent.put(std::move(output),"MuonDTDigis");
  iEvent.put(std::move(output));
  iEvent.put(std::move(outputLinks));

}

void DTDigitizer::storeDigis ( DTWireId &  wireId, TDContainer &  hits, DTDigiCollection &  output, DTDigiSimLinkCollection &  outputLinks  )

private

Definition at line 499 of file DTDigitizer.cc.

References funct::abs(), DTDigi::countsTDC(), muonDTDigis_cfi::deadTime, debug, MuonDigiCollection< IndexType, DigiType >::insertDigi(), DTWireId::layerId(), muonDTDigis_cfi::LinksTimeWindow, muonDTDigis_cfi::MultipleLinks, DTDigi::number(), muonDTDigis_cfi::onlyMuHits, funct::pow(), DTDigi::time(), ntuplemaker::time, and DTWireId::wire().

                                                                                {

  //************ 7A ***************

  // sort signal times
  sort(hits.begin(),hits.end(),hitLessT());

  //************ 7B ***************

  float wakeTime = -999999.0;
  float resolTime = -999999.0;
  int digiN = -1; // Digi identifier within the cell (for multiple digis)
  DTDigi digi;

  // loop over signal times and drop signals inside dead time
  for ( TDContainer::const_iterator hit = hits.begin() ; hit != hits.end() ; 
    hit++ ) {

    if (onlyMuHits && abs((*hit).first->particleType())!=13) continue;

    //************ 7C ***************
    
    float time = (*hit).second;
    if ( time > wakeTime ) {
      // Note that digi is constructed with a float value (in ns)
      int wireN = wireId.wire();
      digiN++;
      digi = DTDigi(wireN, time, digiN);

      // Add association between THIS digi and the corresponding SimTrack
      unsigned int SimTrackId = (*hit).first->trackId();
      EncodedEventId evId = (*hit).first->eventId();
      DTDigiSimLink digisimLink(wireN, digiN, time, SimTrackId, evId);

      if(debug) {
    LogPrint("DTDigitizer")<<endl<<"---- DTDigitizer ----"<<endl;
    LogPrint("DTDigitizer")<<"wireId: "<<wireId<<endl;
    LogPrint("DTDigitizer")<<"sim. time = "<<time<<endl;
    LogPrint("DTDigitizer")<<"digi number = "<< digi.number()<<", digi time = "<<digi.time()
        <<", linked to SimTrack Id = "<<SimTrackId<<endl;
      }

      //************ 7D ***************
      if(digi.countsTDC() < pow(2.,16)){
    DTLayerId layerID = wireId.layerId();  //taking the layer in which reside the wire
        output.insertDigi(layerID, digi); // ordering Digis by layer
        outputLinks.insertDigi(layerID, digisimLink);
    wakeTime = time + deadTime;
    resolTime = time + LinksTimeWindow;
      } 
      else {
    digiN--;
      }
    }
    else if (MultipleLinks && time < resolTime){
      int wireN = wireId.wire();
      unsigned int SimTrackId = (*hit).first->trackId();
      EncodedEventId evId = (*hit).first->eventId();
      DTDigiSimLink digisimLink(wireN, digiN, time, SimTrackId, evId);
      DTLayerId layerID = wireId.layerId(); 
      outputLinks.insertDigi(layerID, digisimLink);

      if(debug) { 
    LogPrint("DTDigitizer")<<"\nAdded multiple link: \n"
        <<"digi number = "<<digi.number()<<", digi time = "<<digi.time()<<" (sim. time = "<<time<<")"
        <<", linked to SimTrack Id = "<<SimTrackId<<endl;
      }
    }
  }
  
}

Friends And Related Function Documentation

friend class DTDigitizerAnalysis

friend

Definition at line 105 of file DTDigitizer.h.

Member Data Documentation

edm::EDGetTokenT<CrossingFrame<PSimHit> > DTDigitizer::cf_token

private

Definition at line 132 of file DTDigitizer.h.

std::string DTDigitizer::collection_for_XF

private

Definition at line 130 of file DTDigitizer.h.

float DTDigitizer::deadTime

private

Definition at line 109 of file DTDigitizer.h.

bool DTDigitizer::debug

private

Definition at line 111 of file DTDigitizer.h.

std::string DTDigitizer::geometryType

private

Definition at line 118 of file DTDigitizer.h.

bool DTDigitizer::IdealModel

private

Definition at line 121 of file DTDigitizer.h.

bool DTDigitizer::interpolate

private

Definition at line 112 of file DTDigitizer.h.

float DTDigitizer::LinksTimeWindow

private

Definition at line 126 of file DTDigitizer.h.

std::string DTDigitizer::mix_

private

Definition at line 129 of file DTDigitizer.h.

bool DTDigitizer::MultipleLinks

private

Definition at line 125 of file DTDigitizer.h.

bool DTDigitizer::onlyMuHits

private

Definition at line 113 of file DTDigitizer.h.

float DTDigitizer::smearing

private

Definition at line 110 of file DTDigitizer.h.

std::string DTDigitizer::syncName

private

Definition at line 115 of file DTDigitizer.h.

float DTDigitizer::theConstVDrift

private

Definition at line 122 of file DTDigitizer.h.

std::unique_ptr<DTDigiSyncBase> DTDigitizer::theSync

private

Definition at line 116 of file DTDigitizer.h.

double DTDigitizer::vPropWire

private

Definition at line 108 of file DTDigitizer.h.