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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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(), HLT_2018_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();
 }