Inheritance diagram for TotemTimingLocalTrackFitter:

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

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 ( const edm::ParameterSet & iConfig )

explicit

Definition at line 45 of file TotemTimingLocalTrackFitter.cc.

References edm::ParameterSet::getParameter(), triggerObjects_cff::id, and trk_algo_map_.

                                                                                          :
   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

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

static

Definition at line 109 of file TotemTimingLocalTrackFitter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   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 );
 }

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

overrideprivate

Definition at line 60 of file TotemTimingLocalTrackFitter.cc.

References edm::Event::getByToken(), maxPlaneActiveChannels_, objects.autophobj::motherId, eostools::move(), MillePedeFileConverter_cfg::out, edm::Event::put(), recHitsToken_, TotemTimingDetId::setChannel(), CTPPSDetId::setStation(), and trk_algo_map_.

 {
   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();
 }