CTPPSDiamondRecHitProducerAlgorithm Class Reference

#include <CTPPSDiamondRecHitProducerAlgorithm.h>

Public Member Functions
void	build (const CTPPSGeometry *, const edm::DetSetVector< CTPPSDiamondDigi > &, edm::DetSetVector< CTPPSDiamondRecHit > &)
	CTPPSDiamondRecHitProducerAlgorithm (const edm::ParameterSet &conf)

Private Attributes
double	ts_to_ns_
	Conversion constant between HPTDC time slice and absolute time (in ns) More...

Detailed Description

Definition at line 23 of file CTPPSDiamondRecHitProducerAlgorithm.h.

Constructor & Destructor Documentation

CTPPSDiamondRecHitProducerAlgorithm::CTPPSDiamondRecHitProducerAlgorithm

(

const edm::ParameterSet &

conf

)

Definition at line 13 of file CTPPSDiamondRecHitProducerAlgorithm.cc.

                                                                                                         :
  ts_to_ns_( iConfig.getParameter<double>( "timeSliceNs" ) )
{}

Member Function Documentation

void CTPPSDiamondRecHitProducerAlgorithm::build	(	const CTPPSGeometry *	geom,
		const edm::DetSetVector< CTPPSDiamondDigi > &	input,
		edm::DetSetVector< CTPPSDiamondRecHit > &	output
	)

Definition at line 18 of file CTPPSDiamondRecHitProducerAlgorithm.cc.

References edm::DetSetVector< T >::find_or_insert(), CTPPSGeometry::getSensor(), DetGeomDesc::params(), DetGeomDesc::parents(), edm::DetSet< T >::push_back(), lumiQTWidget::t, cscNeutronWriter_cfi::t0, CTPPSDiamondRecHit::TIMESLICE_WITHOUT_LEADING, DetGeomDesc::translation(), and ts_to_ns_.

Referenced by CTPPSDiamondRecHitProducer::produce().

{
  for ( const auto& vec : input ) {
    const CTPPSDiamondDetId detid( vec.detId() );

    if ( detid.channel() > 20 ) continue; // VFAT-like information, to be ignored

    // retrieve the geometry element associated to this DetID
    const DetGeomDesc* det = geom->getSensor( detid );

    const float x_pos = det->translation().x(),
                y_pos = det->translation().y();
    float z_pos = 0.;
    if ( det->parents().empty() )
      edm::LogWarning("CTPPSDiamondRecHitProducerAlgorithm") << "The geometry element for " << detid << " has no parents. Check the geometry hierarchy!";
    else
      z_pos = det->parents()[det->parents().size()-1].absTranslation().z(); // retrieve the plane position;

    const float x_width = 2.0 * det->params().at( 0 ), // parameters stand for half the size
                y_width = 2.0 * det->params().at( 1 ),
                z_width = 2.0 * det->params().at( 2 );

    edm::DetSet<CTPPSDiamondRecHit>& rec_hits = output.find_or_insert( detid );

    for ( const auto& digi : vec ) {
      if ( digi.getLeadingEdge() == 0 && digi.getTrailingEdge() == 0 ) continue;

      const int t = digi.getLeadingEdge();
      const int t0 = t % 1024;
      const int time_slice = ( t != 0 ) ? t / 1024 : CTPPSDiamondRecHit::TIMESLICE_WITHOUT_LEADING;

      int tot = 0;
      if ( t != 0 && digi.getTrailingEdge() != 0 ) tot = ( (int)digi.getTrailingEdge() ) - t;

      rec_hits.push_back( CTPPSDiamondRecHit( x_pos, x_width, y_pos, y_width, z_pos, z_width, // spatial information
                                              ( t0 * ts_to_ns_ ),
                                              ( tot * ts_to_ns_),
                                              0., // time precision
                                              time_slice,
                                              digi.getHPTDCErrorFlags(),
                                              digi.getMultipleHit() ) );
    }
  }
}

Member Data Documentation

double CTPPSDiamondRecHitProducerAlgorithm::ts_to_ns_

private

Conversion constant between HPTDC time slice and absolute time (in ns)

Definition at line 32 of file CTPPSDiamondRecHitProducerAlgorithm.h.

Referenced by build().