Inheritance diagram for TotemTimingLocalTrackFitter:

Public Member Functions
	TotemTimingLocalTrackFitter (const edm::ParameterSet &)

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &)

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
const int	maxPlaneActiveChannels_

const edm::EDGetTokenT< edm::DetSetVector< TotemTimingRecHit > >	recHitsToken_

std::map< TotemTimingDetId, TotemTimingTrackRecognition >	trk_algo_map_

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

Definition at line 30 of file TotemTimingLocalTrackFitter.cc.

Constructor & Destructor Documentation

◆ TotemTimingLocalTrackFitter()

TotemTimingLocalTrackFitter::TotemTimingLocalTrackFitter ( const edm::ParameterSet & iConfig )

explicit

Definition at line 44 of file TotemTimingLocalTrackFitter.cc.

     : recHitsToken_(consumes<edm::DetSetVector<TotemTimingRecHit> >(iConfig.getParameter<edm::InputTag>("recHitsTag"))),
       maxPlaneActiveChannels_(iConfig.getParameter<int>("maxPlaneActiveChannels")) {
   produces<edm::DetSetVector<TotemTimingLocalTrack> >();
  
   for (unsigned short armNo = 0; armNo < 2; armNo++)
     for (unsigned short rpNo = 0; rpNo < 2; rpNo++) {
       TotemTimingDetId id(armNo, 1, rpNo, 0, 0);
       TotemTimingTrackRecognition trk_algo(iConfig.getParameter<edm::ParameterSet>("trackingAlgorithmParams"));
       trk_algo_map_.insert(std::make_pair(id, trk_algo));
     }
 }

Member Function Documentation

◆ fillDescriptions()

void TotemTimingLocalTrackFitter::fillDescriptions ( edm::ConfigurationDescriptions & descr )

static

Definition at line 104 of file TotemTimingLocalTrackFitter.cc.

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("recHitsTag", edm::InputTag("totemTimingRecHits"))
       ->setComment("input rechits collection to retrieve");
   desc.add<int>("maxPlaneActiveChannels", 2)->setComment("threshold for discriminating noisy planes");
  
   edm::ParameterSetDescription trackingAlgoParams;
   trackingAlgoParams.add<double>("threshold", 1.5)
       ->setComment("minimal number of rechits to be observed before launching the track recognition algorithm");
   trackingAlgoParams.add<double>("thresholdFromMaximum", 0.5)
       ->setComment("threshold relative to hit profile function local maximum for determining the width of the track");
   trackingAlgoParams.add<double>("resolution", 0.01 /* mm */)
       ->setComment("spatial resolution on the horizontal coordinate (in mm)");
   trackingAlgoParams.add<double>("sigma", 0.)
       ->setComment("pixel efficiency function parameter determining the smoothness of the step");
   trackingAlgoParams.add<double>("tolerance", 0.1 /* mm */)
       ->setComment("tolerance used for checking if the track contains certain hit");
  
   trackingAlgoParams.add<std::string>("pixelEfficiencyFunction", "(x>[0]-0.5*[1]-0.05)*(x<[0]+0.5*[1]-0.05)+0*[2]")
       ->setComment(
           "efficiency function for single pixel\n"
           "can be defined as:\n"
           " * Precise: "
           "(TMath::Erf((x-[0]+0.5*([1]-0.05))/([2]/4)+2)+1)*TMath::Erfc((x-[0]-0.5*([1]-0.05))/([2]/4)-2)/4\n"
           " * Fast: "
           "(x>[0]-0.5*([1]-0.05))*(x<[0]+0.5*([1]-0.05))+((x-[0]+0.5*([1]-0.05)+[2])/"
           "[2])*(x>[0]-0.5*([1]-0.05)-[2])*(x<[0]-0.5*([1]-0.05))+(2-(x-[0]-0.5*([1]-0.05)+[2])/"
           "[2])*(x>[0]+0.5*([1]-0.05))*(x<[0]+0.5*([1]-0.05)+[2])\n"
           " * Legacy: (1/(1+exp(-(x-[0]+0.5*([1]-0.05))/[2])))*(1/(1+exp((x-[0]-0.5*([1]-0.05))/[2])))\n"
           " * Default (sigma ignored): (x>[0]-0.5*[1]-0.05)*(x<[0]+0.5*[1]-0.05)+0*[2]\n"
           "with:\n"
           "  [0]: centre of pad\n"
           "  [1]: width of pad\n"
           "  [2]: sigma: distance between efficiency ~100 -> 0 outside width");
  
   trackingAlgoParams.add<double>("yPosition", 0.0)->setComment("vertical offset of the outcoming track centre");
   trackingAlgoParams.add<double>("yWidth", 0.0)->setComment("vertical track width");
   desc.add<edm::ParameterSetDescription>("trackingAlgorithmParams", trackingAlgoParams)
       ->setComment("list of parameters associated to the track recognition algorithm");
  
   descr.add("totemTimingLocalTracks", desc);
 }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, edm::ParameterDescriptionNode::setComment(), and AlCaHLTBitMon_QueryRunRegistry::string.

◆ produce()

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

overrideprivate

Definition at line 57 of file TotemTimingLocalTrackFitter.cc.

                                                                                        {
   auto pOut = std::make_unique<edm::DetSetVector<TotemTimingLocalTrack> >();
  
   edm::Handle<edm::DetSetVector<TotemTimingRecHit> > recHits;
   iEvent.getByToken(recHitsToken_, recHits);
  
   for (const auto& trk_algo_entry : trk_algo_map_)
     pOut->find_or_insert(trk_algo_entry.first);
  
   std::map<TotemTimingDetId, int> planeActivityMap;
  
   auto motherId = [](const edm::det_id_type& detid) {
     TotemTimingDetId out(detid);
     out.setStation(1);
     out.setChannel(0);
     return out;
   };
  
   for (const auto& vec : *recHits)
     planeActivityMap[motherId(vec.detId())] += vec.size();
  
   // feed hits to the track producers
   for (const auto& vec : *recHits) {
     auto detId = motherId(vec.detId());
     if (planeActivityMap[detId] > maxPlaneActiveChannels_)
       continue;
  
     detId.setPlane(0);
     for (const auto& hit : vec) {
       if (trk_algo_map_.find(detId) == trk_algo_map_.end())
         throw cms::Exception("TotemTimingLocalTrackFitter")
             << "Invalid detId for rechit: arm=" << detId.arm() << ", rp=" << detId.rp();
       trk_algo_map_.find(detId)->second.addHit(hit);
     }
   }
  
   // retrieves tracks for all hit sets
   for (auto& trk_algo_entry : trk_algo_map_)
     trk_algo_entry.second.produceTracks(pOut->operator[](trk_algo_entry.first));
  
   iEvent.put(std::move(pOut));
  
   // remove all hits from the track producers to prepare for the next event
   for (auto& trk_algo_entry : trk_algo_map_)
     trk_algo_entry.second.clear();
 }