ApvTimingAnalysis Class Reference

Analysis for timing run using APV tick marks. More...

#include <ApvTimingAnalysis.h>

Inheritance diagram for ApvTimingAnalysis:

Public Member Functions
void	addErrorCode (const std::string &error) override
	ApvTimingAnalysis (const uint32_t &key)
	ApvTimingAnalysis ()
const float &	base () const
const float &	delay () const
const float &	error () const
bool	foundTickMark () const
uint16_t	frameFindingThreshold () const
const float &	height () const
bool	isValid () const override
const float &	optimumSamplingPoint () const
const float &	peak () const
void	print (std::stringstream &, uint32_t not_used=0) override
const float &	refTime () const
void	refTime (const float &time, const float &targetDelay=-1)
void	reset () override
const float &	time () const
	~ApvTimingAnalysis () override
Public Member Functions inherited from CommissioningAnalysis
	CommissioningAnalysis (const uint32_t &key, const std::string &my_name)
	CommissioningAnalysis (const std::string &my_name)
const uint32_t &	dcuId () const
void	dcuId (const uint32_t &)
const uint32_t &	detId () const
void	detId (const uint32_t &)
const uint32_t &	fecKey () const
void	fecKey (const uint32_t &)
const uint32_t &	fedKey () const
void	fedKey (const uint32_t &)
const VString &	getErrorCodes () const
virtual void	header (std::stringstream &) const
const std::string &	myName () const
void	setErrorCodes (const VString &errors)
virtual void	summary (std::stringstream &) const
virtual	~CommissioningAnalysis ()

Static Public Attributes
static const float	frameFindingThreshold_ = (2. / 3.)
static const float	optimumSamplingPoint_ = 15.
static const float	tickMarkHeightThreshold_ = 50.

Private Attributes
float	base_
float	delay_
float	error_
float	height_
float	peak_
bool	synchronized_
float	time_

Static Private Attributes
static float	refTime_ = 1. * sistrip::invalid_

Friends
class	ApvTimingAlgorithm

Additional Inherited Members
Public Types inherited from CommissioningAnalysis
typedef std::vector< bool >	VBool
typedef std::vector< float >	VFloat
typedef std::vector< uint16_t >	VInt
typedef std::vector< std::string >	VString
typedef std::vector< VBool >	VVBool
typedef std::vector< VFloat >	VVFloat
typedef std::vector< VInt >	VVInt

Detailed Description

Analysis for timing run using APV tick marks.

Author

M. Wingham, R.Bainbridge

Definition at line 15 of file ApvTimingAnalysis.h.

Constructor & Destructor Documentation

ApvTimingAnalysis() [1/2]

ApvTimingAnalysis::ApvTimingAnalysis

(

const uint32_t &

key

)

Definition at line 30 of file ApvTimingAnalysis.cc.

    : CommissioningAnalysis(key, sistrip::apvTimingAnalysis_),
      time_(1. * sistrip::invalid_),
      error_(1. * sistrip::invalid_),
      delay_(1. * sistrip::invalid_),
      height_(1. * sistrip::invalid_),
      base_(1. * sistrip::invalid_),
      peak_(1. * sistrip::invalid_),
      synchronized_(false) {
  ;
}

ApvTimingAnalysis() [2/2]

ApvTimingAnalysis::ApvTimingAnalysis

(

)

Definition at line 44 of file ApvTimingAnalysis.cc.

    : CommissioningAnalysis(sistrip::apvTimingAnalysis_),
      time_(1. * sistrip::invalid_),
      error_(1. * sistrip::invalid_),
      delay_(1. * sistrip::invalid_),
      height_(1. * sistrip::invalid_),
      base_(1. * sistrip::invalid_),
      peak_(1. * sistrip::invalid_),
      synchronized_(false) {
  ;
}

~ApvTimingAnalysis()

ApvTimingAnalysis::~ApvTimingAnalysis ( )

inlineoverride

Definition at line 23 of file ApvTimingAnalysis.h.

{ ; }

Member Function Documentation

addErrorCode()

void ApvTimingAnalysis::addErrorCode ( const std::string &  error )

inlineoverridevirtual

Adds error codes for analysis (overrides private base).

Reimplemented from CommissioningAnalysis.

Definition at line 125 of file ApvTimingAnalysis.h.

References CommissioningAnalysis::addErrorCode(), and error().

Referenced by refTime().

{ CommissioningAnalysis::addErrorCode(error); }

base()

const float & ApvTimingAnalysis::base ( ) const

inline

Baseline level of tick mark [ADC].

Definition at line 123 of file ApvTimingAnalysis.h.

References base_.

{ return base_; }

delay()

const float & ApvTimingAnalysis::delay ( ) const

inline

Delay required to sync w.r.t. reference tick mark [ns].

Definition at line 121 of file ApvTimingAnalysis.h.

References delay_.

{ return delay_; }

error()

const float & ApvTimingAnalysis::error ( ) const

inline

Error on time of tick mark rising edge [ns].

Definition at line 118 of file ApvTimingAnalysis.h.

References error_.

Referenced by addErrorCode(), and Page1Parser.Page1Parser::check_for_whole_start_tag().

{ return error_; }

foundTickMark()

bool ApvTimingAnalysis::foundTickMark

(

)

const

Identifies if tick mark is found or not.

Definition at line 133 of file ApvTimingAnalysis.cc.

References base_, relativeConstraints::empty, frameFindingThreshold(), CommissioningAnalysis::getErrorCodes(), height_, peak_, time_, and sistrip::valid_.

Referenced by print().

                                            {
  return ((getErrorCodes().empty() || getErrorCodes()[0] == "TickMarkRecovered") && time_ < sistrip::valid_ &&
          base_ < sistrip::valid_ && peak_ < sistrip::valid_ && height_ < sistrip::valid_ &&
          frameFindingThreshold() < sistrip::valid_);
}

frameFindingThreshold()

uint16_t ApvTimingAnalysis::frameFindingThreshold

(

)

const

FED frame-finding threshold [ADC] (returns 65535 if invalid).

Definition at line 121 of file ApvTimingAnalysis.cc.

References base_, relativeConstraints::empty, frameFindingThreshold_, CommissioningAnalysis::getErrorCodes(), height_, sistrip::invalid_, peak_, tickMarkHeightThreshold_, time_, and sistrip::valid_.

Referenced by foundTickMark(), isValid(), and print().

                                                        {
  if ((getErrorCodes().empty() || getErrorCodes()[0] == "TickMarkRecovered") && time_ < sistrip::valid_ &&
      base_ < sistrip::valid_ && peak_ < sistrip::valid_ && height_ < sistrip::valid_ &&
      height_ > tickMarkHeightThreshold_) {
    return ((static_cast<uint16_t>(base_ + height_ * ApvTimingAnalysis::frameFindingThreshold_) / 32) * 32);
  } else {
    return sistrip::invalid_;
  }
}

height()

const float & ApvTimingAnalysis::height ( ) const

inline

Height of tick mark [ADC].

Definition at line 122 of file ApvTimingAnalysis.h.

References height_.

{ return height_; }

isValid()

bool ApvTimingAnalysis::isValid ( void  ) const

overridevirtual

Identifies if analysis is valid or not.

Reimplemented from CommissioningAnalysis.

Definition at line 141 of file ApvTimingAnalysis.cc.

References base_, delay_, relativeConstraints::empty, frameFindingThreshold(), CommissioningAnalysis::getErrorCodes(), height_, peak_, refTime_, synchronized_, time_, and sistrip::valid_.

Referenced by ntupleDataFormat._Object::_checkIsValid(), print(), and core.AutoHandle.AutoHandle::ReallyLoad().

                                      {
  return ((getErrorCodes().empty() || getErrorCodes()[0] == "TickMarkRecovered") && time_ < sistrip::valid_ &&
          base_ < sistrip::valid_ && peak_ < sistrip::valid_ && height_ < sistrip::valid_ &&
          frameFindingThreshold() < sistrip::valid_ && synchronized_ && refTime_ < sistrip::valid_ &&
          delay_ < sistrip::valid_);
}

optimumSamplingPoint()

const float & ApvTimingAnalysis::optimumSamplingPoint ( ) const

inline

Optimum sampling point, defined w.r.t. rising edge [ns].

Definition at line 119 of file ApvTimingAnalysis.h.

References optimumSamplingPoint_.

{ return optimumSamplingPoint_; }

peak()

const float & ApvTimingAnalysis::peak ( ) const

inline

Level of tick mark top [ADC].

Definition at line 124 of file ApvTimingAnalysis.h.

References peak_.

{ return peak_; }

print()

void ApvTimingAnalysis::print ( std::stringstream &  ss, uint32_t  not_used = 0  )

overridevirtual

Prints analysis results.

Implements CommissioningAnalysis.

Definition at line 150 of file ApvTimingAnalysis.cc.

References base_, delay_, relativeConstraints::empty, alignBH_cfg::fixed, foundTickMark(), frameFindingThreshold(), CommissioningAnalysis::getErrorCodes(), CommissioningAnalysis::header(), height_, sistrip::invalid_, isValid(), optimumSamplingPoint_, peak_, refTime_, contentValuesCheck::ss, time_, and sistrip::valid_.

                                                                    {
  header(ss);

  float sampling1 = sistrip::invalid_;
  if (time_ <= sistrip::valid_) {
    sampling1 = time_ + optimumSamplingPoint_;
  }

  float sampling2 = sistrip::invalid_;
  if (refTime_ <= sistrip::valid_) {
    sampling2 = refTime_ + optimumSamplingPoint_;
  }

  float adjust = sistrip::invalid_;
  if (sampling1 <= sistrip::valid_ && delay_ <= sistrip::valid_) {
    adjust = sampling1 + delay_;
  }

  ss << std::fixed << std::setprecision(2) << " Tick mark: time of rising edge     [ns] : " << time_
     << std::endl
     //<< " Error on time of rising edge     [ns] : " << error_ << std::endl
     << " Tick mark: time of sampling point  [ns] : " << sampling1 << std::endl
     << " Ref tick: time of rising edge      [ns] : " << refTime_ << std::endl
     << " Ref tick: time of sampling point   [ns] : " << sampling2 << std::endl
     << " Ref tick: adjusted sampling point  [ns] : " << adjust << std::endl
     << " Delay required to synchronise      [ns] : " << delay_ << std::endl
     << " Tick mark bottom (baseline)       [ADC] : " << base_ << std::endl
     << " Tick mark top                     [ADC] : " << peak_ << std::endl
     << " Tick mark height                  [ADC] : " << height_ << std::endl
     << " Frame finding threshold           [ADC] : " << frameFindingThreshold() << std::endl
     << std::boolalpha << " Tick mark found                         : " << foundTickMark() << std::endl
     << " isValid                                 : " << isValid() << std::endl
     << std::noboolalpha << " Error codes (found " << std::setw(3) << std::setfill(' ') << getErrorCodes().size()
     << ")                 : ";
  if (getErrorCodes().empty()) {
    ss << "(none)";
  } else {
    VString::const_iterator istr = getErrorCodes().begin();
    VString::const_iterator jstr = getErrorCodes().end();
    for (; istr != jstr; ++istr) {
      ss << *istr << " ";
    }
  }
  ss << std::endl;
}

refTime() [1/2]

const float & ApvTimingAnalysis::refTime ( ) const

inline

Sampling point of "reference" tick mark [ns].

Definition at line 120 of file ApvTimingAnalysis.h.

References refTime_.

{ return refTime_; }

refTime() [2/2]

void ApvTimingAnalysis::refTime	(	const float &	time,
		const float &	targetDelay = -1
	)

Static method to set global reference time [ns].

Definition at line 70 of file ApvTimingAnalysis.cc.

References addErrorCode(), delay_, sistrip::invalid_, sistrip::invalidDelayTime_, sistrip::invalidRefTime_, sistrip::mlCommissioning_, CommissioningAnalysis::myName(), optimumSamplingPoint_, position, refTime_, synchronized_, time(), time_, and sistrip::valid_.

                                                                           {
  // Checks synchronization to reference time is done only once
  if (synchronized_) {
    edm::LogWarning(mlCommissioning_) << "[" << myName() << "::" << __func__ << "]"
                                      << " Attempting to re-synchronize with reference time!"
                                      << " Not allowed!";
    return;
  }
  synchronized_ = true;

  // Set reference time and check if tick mark time is valid
  refTime_ = time;
  if (time_ > sistrip::valid_) {
    return;
  }

  // Calculate position of "sampling point" of last tick;
  int32_t position;
  if (targetDelay == -1) {  // by default use latest tick
    position = static_cast<int32_t>(rint(refTime_ + optimumSamplingPoint_));
  } else {
    position = static_cast<int32_t>(rint(targetDelay + optimumSamplingPoint_));
  }

  // Calculate adjustment so that sampling point is multiple of 25 (ie, synched with FED sampling)
  float adjustment = 25 - position % 25;

  // Calculate delay required to synchronise with this adjusted sampling position
  if (targetDelay == -1) {  // by default align forward to the latest tick
    delay_ = (refTime_ + adjustment) - time_;
  } else {  // otherwise use the supplied target delay
    if (adjustment > 25 / 2)
      adjustment -= 25;  // go as close as possible to desired target
    delay_ = (targetDelay + adjustment) - time_;
  }

  // Check reference time
  if (refTime_ < 0. || refTime_ > sistrip::valid_) {
    refTime_ = sistrip::invalid_;
    addErrorCode(sistrip::invalidRefTime_);
  }

  // Check delay is valid
  if (delay_ < -sistrip::valid_ || delay_ > sistrip::valid_) {
    delay_ = sistrip::invalid_;
    addErrorCode(sistrip::invalidDelayTime_);
  }
}

reset()

void ApvTimingAnalysis::reset ( void  )

overridevirtual

Resets analysis member data.

Implements CommissioningAnalysis.

Definition at line 58 of file ApvTimingAnalysis.cc.

References base_, delay_, error_, height_, sistrip::invalid_, peak_, synchronized_, and time_.

Referenced by MatrixReader.MatrixReader::__init__(), and MatrixReader.MatrixReader::showRaw().

                              {
  time_ = 1. * sistrip::invalid_;
  error_ = 1. * sistrip::invalid_;
  delay_ = 1. * sistrip::invalid_;
  height_ = 1. * sistrip::invalid_;
  base_ = 1. * sistrip::invalid_;
  peak_ = 1. * sistrip::invalid_;
  synchronized_ = false;
}

time()

const float & ApvTimingAnalysis::time ( ) const

inline

Time of tick mark rising edge [ns].

Definition at line 117 of file ApvTimingAnalysis.h.

References time_.

Referenced by refTime().

{ return time_; }

Friends And Related Function Documentation

ApvTimingAlgorithm

friend class ApvTimingAlgorithm

friend

Definition at line 25 of file ApvTimingAnalysis.h.

Member Data Documentation

base_

float ApvTimingAnalysis::base_

private

Baseline level of tick mark [ADC].

Definition at line 106 of file ApvTimingAnalysis.h.

Referenced by base(), foundTickMark(), frameFindingThreshold(), isValid(), print(), and reset().

delay_

float ApvTimingAnalysis::delay_

private

Delay required to sync w.r.t. reference tick mark [ns].

Definition at line 100 of file ApvTimingAnalysis.h.

Referenced by delay(), isValid(), print(), refTime(), and reset().

error_

float ApvTimingAnalysis::error_

private

Error on time of tick mark rising edge [ns].

Definition at line 94 of file ApvTimingAnalysis.h.

Referenced by error(), and reset().

frameFindingThreshold_

const float ApvTimingAnalysis::frameFindingThreshold_ = (2. / 3.)

static

Threshold for FED frame finding (fraction of tick height).

Definition at line 85 of file ApvTimingAnalysis.h.

Referenced by frameFindingThreshold().

height_

float ApvTimingAnalysis::height_

private

Height of tick mark [ADC].

Definition at line 103 of file ApvTimingAnalysis.h.

Referenced by foundTickMark(), frameFindingThreshold(), height(), isValid(), print(), and reset().

optimumSamplingPoint_

const float ApvTimingAnalysis::optimumSamplingPoint_ = 15.

static

Optimum sampling point, defined w.r.t. rising edge [ns].

Definition at line 79 of file ApvTimingAnalysis.h.

Referenced by optimumSamplingPoint(), print(), and refTime().

peak_

float ApvTimingAnalysis::peak_

private

Level of tick mark top [ADC].

Definition at line 109 of file ApvTimingAnalysis.h.

Referenced by foundTickMark(), frameFindingThreshold(), isValid(), peak(), print(), and reset().

refTime_

float ApvTimingAnalysis::refTime_ = 1. * sistrip::invalid_

staticprivate

Sampling point of "reference" tick mark [ns].

Definition at line 97 of file ApvTimingAnalysis.h.

Referenced by isValid(), print(), and refTime().

synchronized_

bool ApvTimingAnalysis::synchronized_

private

Checks synchronization to ref time is done only once.

Definition at line 112 of file ApvTimingAnalysis.h.

Referenced by isValid(), refTime(), and reset().

tickMarkHeightThreshold_

const float ApvTimingAnalysis::tickMarkHeightThreshold_ = 50.

static

Threshold defining minimum tick mark height [ADC].

Definition at line 82 of file ApvTimingAnalysis.h.

Referenced by ApvTimingAlgorithm::analyse(), ApvTimingHistosUsingDb::create(), and frameFindingThreshold().

time_

float ApvTimingAnalysis::time_

private

Time of tick mark rising edge [ns].

Definition at line 91 of file ApvTimingAnalysis.h.

Referenced by foundTickMark(), frameFindingThreshold(), isValid(), print(), refTime(), reset(), and time().