RPCSim.cc

Go to the documentation of this file.

#include "RPCSim.h"
#include "SimMuon/RPCDigitizer/src/RPCSimSetUp.h"
#include "SimDataFormats/EncodedEventId/interface/EncodedEventId.h"

RPCSim::RPCSim(const edm::ParameterSet& config)
{

}

void
RPCSim::fillDigis(int rollDetId, RPCDigiCollection& digis)
{
  //  theRpcDigiSimLinks.clear();

//  std::vector<std::pair<int,int> > vdigi;
//  vdigi.clear();

  for (std::set< std::pair<int,int> >::iterator i=strips.begin();
       i!=strips.end(); i++){
    if(i->second != -999){
      RPCDigi rpcDigi(i->first,i->second);
 
      //NCA
      digis.insertDigi(RPCDetId(rollDetId),rpcDigi);
      this->addLinks(i->first,i->second);
    }
  }
  strips.clear();

  for (auto it: irpc_digis){
     if(it.bx() != -999){
        digis.insertDigi(RPCDetId(rollDetId),it);
        this->addLinks(it.strip(),it.bx());
     }
  }
  irpc_digis.clear();
}

void RPCSim::addLinks(unsigned int strip, int bx) {

  std::pair<unsigned int, int > digi(strip, bx);
  std::pair<DetectorHitMap::iterator, DetectorHitMap::iterator> channelHitItr 
     = theDetectorHitMap.equal_range(digi);

  for( DetectorHitMap::iterator hitItr = channelHitItr.first; 
                                hitItr != channelHitItr.second; ++hitItr){
    const PSimHit * hit = (hitItr->second);

    if(hit != nullptr) {
      theRpcDigiSimLinks.push_back( RPCDigiSimLink(digi, hit->entryPoint(),hit->momentumAtEntry(),
                           hit->timeOfFlight(),hit->energyLoss(),hit->particleType(), 
                           hit->detUnitId(), hit->trackId(), hit->eventId(), hit->processType() ) );

    }
  }
}