Concrete hardware channel for HAPPEX 18-bit ADC modules. More...

#include <QwADC18_Channel.h>

Inheritance diagram for QwADC18_Channel:

Collaboration diagram for QwADC18_Channel:

Public Member Functions
	QwADC18_Channel ()

	QwADC18_Channel (TString name, TString datatosave="raw")

	QwADC18_Channel (const QwADC18_Channel &value)

	QwADC18_Channel (const QwADC18_Channel &value, VQwDataElement::EDataToSave datatosave)

	~QwADC18_Channel () override

VQwHardwareChannel *	Clone (VQwDataElement::EDataToSave datatosave) const override

void	InitializeChannel (TString name, TString datatosave) override
	Initialize the fields in this object.

void	InitializeChannel (TString subsystem, TString instrumenttype, TString name, TString datatosave) override
	Initialize the fields in this object.

void	LoadChannelParameters (QwParameterFile &paramfile) override

void	SetDefaultSampleSize (size_t num_samples_map)

void	ClearEventData () override
	Clear the event data in this element.

void	RandomizeEventData (int helicity=0.0, double time=0.0) override
	Internally generate random event data.

void	ForceMapfileSampleSize ()

void	SmearByResolution (double resolution) override

void	SetEventData (Double_t value)

void	SetRawEventData () override

void	EncodeEventData (std::vector< UInt_t > &buffer) override
	Encode the event data into a CODA buffer.

Bool_t	IsHeaderWord (UInt_t word) const
	Decode the event data from a CODA buffer.

Int_t	ProcessDataWord (UInt_t word)

Int_t	ProcessEvBuffer (UInt_t *buffer, UInt_t num_words_left, UInt_t index=0) override
	Process the CODA event buffer for this element.

void	ProcessEvent () override
	Process the event data according to pedestal and calibration factor.

QwADC18_Channel &	operator= (const QwADC18_Channel &value)

void	AssignScaledValue (const QwADC18_Channel &value, Double_t scale)

void	AssignValueFrom (const VQwDataElement *valueptr) override

void	AddValueFrom (const VQwHardwareChannel *valueptr) override

void	SubtractValueFrom (const VQwHardwareChannel *valueptr) override

void	MultiplyBy (const VQwHardwareChannel *valueptr) override

void	DivideBy (const VQwHardwareChannel *valueptr) override

void	ArcTan (const QwADC18_Channel &value)

QwADC18_Channel &	operator+= (const QwADC18_Channel &value)

QwADC18_Channel &	operator-= (const QwADC18_Channel &value)

QwADC18_Channel &	operator*= (const QwADC18_Channel &value)

VQwHardwareChannel &	operator+= (const VQwHardwareChannel &input) override

VQwHardwareChannel &	operator-= (const VQwHardwareChannel &input) override

VQwHardwareChannel &	operator*= (const VQwHardwareChannel &input) override

VQwHardwareChannel &	operator/= (const VQwHardwareChannel &input) override

const QwADC18_Channel	operator+ (const QwADC18_Channel &value) const

const QwADC18_Channel	operator- (const QwADC18_Channel &value) const

const QwADC18_Channel	operator* (const QwADC18_Channel &value) const

void	Sum (const QwADC18_Channel &value1, const QwADC18_Channel &value2)

void	Difference (const QwADC18_Channel &value1, const QwADC18_Channel &value2)

void	Ratio (const QwADC18_Channel &numer, const QwADC18_Channel &denom)

void	Product (const QwADC18_Channel &value1, const QwADC18_Channel &value2)

void	DivideBy (const QwADC18_Channel &denom)

void	AddChannelOffset (Double_t Offset)

void	Scale (Double_t Offset) override

void	AccumulateRunningSum (const QwADC18_Channel &value, Int_t count=0, Int_t ErrorMask=0xFFFFFFF)

void	AccumulateRunningSum (const VQwHardwareChannel *value, Int_t count=0, Int_t ErrorMask=0xFFFFFFF) override

void	DeaccumulateRunningSum (const QwADC18_Channel &value, Int_t ErrorMask=0xFFFFFFF)

void	CalculateRunningAverage () override

Bool_t	MatchSequenceNumber (size_t seqnum)

Bool_t	MatchNumberOfSamples (size_t numsamp)

Bool_t	ApplySingleEventCuts (Double_t LL, Double_t UL)

Bool_t	ApplySingleEventCuts () override

void	PrintErrorCounters () const override
	report number of events failed due to HW and event cut failure

void	SetADC18SaturationLimt (Double_t sat_volts=8.5)

Double_t	GetADC18SaturationLimt ()

Int_t	ApplyHWChecks () override

void	IncrementErrorCounters () override

Int_t	GetRawValue (size_t element) const override

Double_t	GetValue (size_t element) const override

Double_t	GetValueM2 (size_t element) const override

Double_t	GetValueError (size_t element) const override

void	ConstructHistograms (TDirectory *folder, TString &prefix) override
	Construct the histograms for this data element.

void	FillHistograms () override
	Fill the histograms for this data element.

void	ConstructBranchAndVector (TTree *tree, TString &prefix, QwRootTreeBranchVector &values) override

void	ConstructBranch (TTree *tree, TString &prefix) override

void	FillTreeVector (QwRootTreeBranchVector &values) const override

Double_t	GetAverageVolts () const

size_t	GetSequenceNumber () const

size_t	GetNumberOfSamples () const

void	SetCalibrationToVolts ()

void	PrintValue () const override
	Print single line of value and error of this data element.

void	PrintInfo () const override
	Print multiple lines of information about this data element.

void	Blind (const QwBlinder *blinder)
	Blind this channel as an asymmetry.

void	Blind (const QwBlinder *blinder, const QwADC18_Channel &yield)
	Blind this channel as a difference.

void	ScaledAdd (Double_t scale, const VQwHardwareChannel *value) override

virtual void	LoadMockDataParameters (QwParameterFile &paramfile)
	Load the mock data parameters from the current line in the param file.

Int_t	GetRawValue () const

Double_t	GetValue () const

Double_t	GetValueM2 () const

Double_t	GetValueError () const

Double_t	GetValueWidth () const

Double_t	GetValueWidth (size_t element) const

virtual void	DeaccumulateRunningSum (const VQwHardwareChannel *value, Int_t ErrorMask=0xFFFFFFF)

virtual VQwHardwareChannel *	Clone () const

Public Member Functions inherited from VQwHardwareChannel
	VQwHardwareChannel ()

	VQwHardwareChannel (const VQwHardwareChannel &value)

	VQwHardwareChannel (const VQwHardwareChannel &value, VQwDataElement::EDataToSave datatosave)

	~VQwHardwareChannel () override

void	CopyFrom (const VQwHardwareChannel &value)

void	ProcessOptions ()

size_t	GetNumberOfDataWords ()
	Get the number of data words in this data element.

size_t	GetNumberOfSubelements ()
	Get the number of subelements in this data element.

Int_t	GetRawValue () const

Double_t	GetValue () const

Double_t	GetValueM2 () const

Double_t	GetValueError () const

Double_t	GetValueWidth () const

Double_t	GetValueWidth (size_t element) const

void	ClearEventData () override
	Clear the event data in this element.

void	InitializeChannel (TString name)
	Initialize the fields in this object.

void	SetEventCutMode (Int_t bcuts)

virtual Bool_t	CheckForBurpFail (const VQwHardwareChannel *event)

void	SetSingleEventCuts (Double_t min, Double_t max)
	Set the upper and lower limits (fULimit and fLLimit) for this channel.

void	SetSingleEventCuts (UInt_t errorflag, Double_t min, Double_t max, Double_t stability=-1.0, Double_t BurpLevel=-1.0)
	Inherited from VQwDataElement to set the upper and lower limits (fULimit and fLLimit), stability % and the error flag on this channel.

Double_t	GetEventCutUpperLimit () const

Double_t	GetEventCutLowerLimit () const

Double_t	GetStabilityLimit () const

UInt_t	UpdateErrorFlag () override
	Update the error flag based on the error flags of internally contained objects Return parameter is the "Eventcut Error Flag".

void	UpdateErrorFlag (const VQwHardwareChannel &elem)

virtual UInt_t	GetErrorCode () const

VQwHardwareChannel &	operator= (const VQwHardwareChannel &value)
	Arithmetic assignment operator: Should only copy event-based data.

void	AssignScaledValue (const VQwHardwareChannel &value, Double_t scale)

virtual void	Ratio (const VQwHardwareChannel numer, const VQwHardwareChannel denom)

void	SetPedestal (Double_t ped)

Double_t	GetPedestal () const

void	SetCalibrationFactor (Double_t factor)

Double_t	GetCalibrationFactor () const

void	AddEntriesToList (std::vector< QwDBInterface > &row_list)

virtual void	AddErrEntriesToList (std::vector< QwErrDBInterface > &)

void	ConstructBranch (TTree *tree, TString &prefix, QwParameterFile &modulelist)

virtual void	CopyParameters (const VQwHardwareChannel *)

void	UpdateErrorFlag (const UInt_t &error)

Public Member Functions inherited from VQwDataElement
	VQwDataElement ()
	Default constructor.

	VQwDataElement (const VQwDataElement &value)
	Copy constructor.

	~VQwDataElement () override
	Virtual destructor.

void	CopyFrom (const VQwDataElement &value)

Bool_t	IsNameEmpty () const
	Is the name of this element empty?

void	SetElementName (const TString &name)
	Set the name of this element.

virtual const TString &	GetElementName () const
	Get the name of this element.

size_t	GetNumberOfDataWords ()
	Get the number of data words in this data element.

UInt_t	GetGoodEventCount () const

VQwDataElement &	operator+= (const VQwDataElement &)
	Addition-assignment operator.

VQwDataElement &	operator-= (const VQwDataElement &)
	Subtraction-assignment operator.

void	Sum (const VQwDataElement &, const VQwDataElement &)
	Sum operator (base class fallback throws runtime error)

void	Difference (const VQwDataElement &, const VQwDataElement &)
	Difference operator (base class fallback throws runtime error)

void	Ratio (const VQwDataElement &, const VQwDataElement &)
	Ratio operator (base class fallback throws runtime error)

void	SetSingleEventCuts (UInt_t, Double_t, Double_t, Double_t)
	set the upper and lower limits (fULimit and fLLimit), stability % and the error flag on this channel

Bool_t	CheckForBurpFail (const VQwDataElement *)

virtual UInt_t	GetEventcutErrorFlag ()
	return the error flag on this channel/device

virtual void	SetNeedsExternalClock (Bool_t)

virtual Bool_t	NeedsExternalClock ()

virtual std::string	GetExternalClockName ()

virtual void	SetExternalClockPtr (const VQwHardwareChannel *)

virtual void	SetExternalClockName (const std::string)

virtual Double_t	GetNormClockValue ()

TString	GetSubsystemName () const
	Return the name of the inheriting subsystem name.

void	SetSubsystemName (TString sysname)
	Set the name of the inheriting subsystem name.

TString	GetModuleType () const
	Return the type of the beam instrument.

void	SetModuleType (TString ModuleType)
	set the type of the beam instrument

Public Member Functions inherited from MQwHistograms
void	ShareHistograms (const MQwHistograms *source)
	Share histogram pointers between objects.

Public Member Functions inherited from MQwMockable
	MQwMockable ()

virtual	~MQwMockable ()

void	SetRandomEventDriftParameters (Double_t amplitude, Double_t phase, Double_t frequency)
	Set a single set of harmonic drift parameters.

void	AddRandomEventDriftParameters (Double_t amplitude, Double_t phase, Double_t frequency)
	Add drift parameters to the internal set.

void	SetRandomEventParameters (Double_t mean, Double_t sigma)
	Set the normal random event parameters.

void	SetRandomEventAsymmetry (Double_t asymmetry)
	Set the helicity asymmetry.

Double_t	GetRandomValue ()

void	UseExternalRandomVariable ()
	Set the flag to use an externally provided random variable.

void	SetExternalRandomVariable (Double_t random_variable)
	Set the externally provided random variable.

void	SetMockDataAsDiff ()

Static Public Member Functions
static Int_t	GetBufferOffset (Int_t moduleindex, Int_t channelindex)

static void	PrintErrorCounterHead ()

static void	PrintErrorCounterTail ()

static void	SetBurpHoldoff (Int_t holdoff)

Static Public Attributes
static const Double_t	kTimePerSample = 2.0 * Qw::us

Protected Member Functions
QwADC18_Channel &	operator/= (const QwADC18_Channel &value)

Protected Member Functions inherited from VQwHardwareChannel
void	SetNumberOfDataWords (const UInt_t &numwords)
	Set the number of data words in this data element.

void	SetNumberOfSubElements (const size_t elements)
	Set the number of data words in this data element.

void	SetDataToSave (TString datatosave)
	Set the flag indicating if raw or derived values are in this data element.

void	SetDataToSave (VQwDataElement::EDataToSave datatosave)
	Set the flag indicating if raw or derived values are in this data element.

void	SetDataToSaveByPrefix (const TString &prefix)
	Set the flag indicating if raw or derived values are in this data element based on prefix.

void	RangeCheck (size_t element) const
	Checks that the requested element is in range, to be used in accesses to subelements similar to std::vector::at().

Protected Member Functions inherited from VQwDataElement
void	SetNumberOfDataWords (const UInt_t &numwords)
	Set the number of data words in this data element.

VQwDataElement &	operator= (const VQwDataElement &value)
	Arithmetic assignment operator: Should only copy event-based data.

void	UpdateErrorFlag (const UInt_t &error)

Protected Member Functions inherited from MQwHistograms
	MQwHistograms ()
	Default constructor.

	MQwHistograms (const MQwHistograms &source)
	Copy constructor.

virtual	~MQwHistograms ()
	Virtual destructor.

MQwHistograms &	operator= (const MQwHistograms &value)

void	Fill_Pointer (TH1_ptr hist_ptr, Double_t value)

void	AddHistogram (TH1 *h)
	Register a histogram.

Private Attributes
UInt_t	fDiff_Raw

UInt_t	fBase_Raw

UInt_t	fPeak_Raw

UInt_t	fValue_Raw

Double_t	fValue

Double_t	fValueM2

Double_t	fValueError

QwADC18_Channel *	fRunningSum
	Pointer to the running sum for this channel.

QwADC18_Channel *	fDAC
	Pointer to the DAC channel for this channel.

Static Private Attributes
static const Bool_t	kDEBUG = kFALSE

static const Int_t	kHeaderWordsPerBank = 1

static const Int_t	kFooterWordsPerBank = 1

static const Int_t	kHeaderWordsPerModule = 12

static const Int_t	kFooterWordsPerModule = 1

static const Int_t	kDataWordsPerChannel = 3

static const Int_t	kMaxChannels = 4

static const UInt_t	mask31x = 0x80000000

static const UInt_t	mask3029x = 0x60000000

static const UInt_t	mask2625x = 0x06000000

static const UInt_t	mask2422x = 0x01c00000

static const UInt_t	mask21x = 0x00200000

static const UInt_t	mask200x = 0x001fffff

static const UInt_t	mask2118x = 0x003c0000

static const UInt_t	mask170x = 0x0003ffff

static const UInt_t	mask150x = 0x0000ffff

ADC Calibration
static const Double_t	kADC18_VoltsPerBit = (20./(1<<18))

Friends
std::ostream &	operator<< (std::ostream &stream, const QwADC18_Channel &channel)

Channel configuration data members
UInt_t	fSequenceNumber
	Event sequence number for this channel.

UInt_t	fPreviousSequenceNumber
	Previous event sequence number for this channel.

UInt_t	fNumberOfSamples
	Number of samples read through the module.

UInt_t	fNumberOfSamples_map
	Number of samples in the expected to read through the module. This value is set in the QwBeamline map file.

Int_t	fErrorCount_HWSat
	check to see ADC channel is saturated

Int_t	fErrorCount_sample
	for sample size check

Int_t	fErrorCount_SW_HW
	HW_sum==SW_sum check.

Int_t	fErrorCount_Sequence
	sequence number check

Int_t	fErrorCount_SameHW
	check to see ADC returning same HW value

Int_t	fErrorCount_ZeroHW
	check to see ADC returning zero

Int_t	fNumEvtsWithEventCutsRejected
	Counts the Event cut rejected events.

Int_t	fADC_Same_NumEvt
	Keep track of how many events with same ADC value returned.

Int_t	fSequenceNo_Prev
	Keep the sequence number of the last event.

Int_t	fSequenceNo_Counter
	Internal counter to keep track of the sequence number.

Double_t	fPrev_HardwareBlockSum
	Previous Module-based sum of the four sub-blocks.

Double_t	fSaturationABSLimit
	absolute value of the ADC18 saturation volt

Bool_t	bHw_sum

Bool_t	bHw_sum_raw

Bool_t	bBlock

Bool_t	bBlock_raw

Bool_t	bNum_samples

Bool_t	bDevice_Error_Code

Bool_t	bSequence_number

static const Bool_t	bDEBUG =kFALSE
	debugging display purposes

Additional Inherited Members
Public Types inherited from VQwDataElement
enum	EDataToSave { kRaw = 0 , kDerived , kMoments }

Protected Attributes inherited from VQwHardwareChannel
UInt_t	fNumberOfDataWords
	Number of raw data words in this data element.

UInt_t	fNumberOfSubElements
	Number of subelements in this data element.

EDataToSave	fDataToSave

size_t	fTreeArrayIndex

size_t	fTreeArrayNumEntries

Double_t	fPedestal

Double_t	fCalibrationFactor

Bool_t	kFoundPedestal

Bool_t	kFoundGain

Int_t	bEVENTCUTMODE

Double_t	fULimit

Double_t	fLLimit

Double_t	fStability

Double_t	fBurpThreshold

Int_t	fBurpCountdown

Protected Attributes inherited from VQwDataElement
TString	fElementName
	Name of this data element.

UInt_t	fNumberOfDataWords
	Number of raw data words in this data element.

UInt_t	fGoodEventCount
	Number of good events accumulated in this element.

TString	fSubsystemName

TString	fModuleType

UInt_t	fErrorFlag
	This the standard error code generated for the channel that contains the global/local/stability flags and the Device error code (Unique error code for HW failures)

UInt_t	fErrorConfigFlag
	contains the global/local/stability flags

Protected Attributes inherited from MQwHistograms
std::vector< TH1_ptr >	fHistograms
	Histograms associated with this data element.

bool	fUseExternalRandomVariable
	Flag to use an externally provided normal random variable.

double	fExternalRandomVariable
	Externally provided normal random variable.

bool	fCalcMockDataAsDiff

Double_t	fMockAsymmetry
	Helicity asymmetry.

Double_t	fMockGaussianMean
	Mean of normal distribution.

Double_t	fMockGaussianSigma
	Sigma of normal distribution.

std::vector< Double_t >	fMockDriftAmplitude
	Harmonic drift amplitude.

std::vector< Double_t >	fMockDriftFrequency
	Harmonic drift frequency.

std::vector< Double_t >	fMockDriftPhase
	Harmonic drift phase.

static Int_t	fBurpHoldoff = 10

static boost::mt19937	fRandomnessGenerator
	Internal randomness generator.

static boost::normal_distribution< double >	fNormalDistribution
	Internal normal probability distribution.

static boost::variate_generator< boost::mt19937, boost::normal_distribution< double > >	fNormalRandomVariable
	Internal normal random variable.

Detailed Description

Concrete hardware channel for HAPPEX 18-bit ADC modules.

Decodes and processes data from HAPPEX 18-bit ADC channels, providing access to raw and calibrated values, statistical moments, single-event cuts, and running statistics. Implements the dual-operator pattern for both type-specific and polymorphic operations.

Definition at line 34 of file QwADC18_Channel.h.

Constructor & Destructor Documentation

◆ QwADC18_Channel() [1/4]

QwADC18_Channel::QwADC18_Channel ( )

inline

Definition at line 61 of file QwADC18_Channel.h.

                   : MQwMockable() {
    InitializeChannel("","");
    SetADC18SaturationLimt(8.5);//FIXME set the default saturation limit
  };

References InitializeChannel(), MQwMockable::MQwMockable(), and SetADC18SaturationLimt().

Referenced by AccumulateRunningSum(), AccumulateRunningSum(), AddValueFrom(), ArcTan(), AssignScaledValue(), AssignValueFrom(), Blind(), Clone(), DeaccumulateRunningSum(), Difference(), DivideBy(), DivideBy(), MultiplyBy(), operator*(), operator*=(), operator*=(), operator+(), operator+=(), operator+=(), operator-(), operator-=(), operator-=(), operator/=(), operator/=(), operator<<, operator=(), Product(), QwADC18_Channel(), QwADC18_Channel(), Ratio(), ScaledAdd(), SubtractValueFrom(), and Sum().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ QwADC18_Channel() [2/4]

QwADC18_Channel::QwADC18_Channel	(	TString	name,
		TString	datatosave = "raw" )

inline

Definition at line 65 of file QwADC18_Channel.h.

                                                           : MQwMockable() {
    InitializeChannel(name, datatosave);
    SetADC18SaturationLimt(8.5);//FIXME set the default saturation limit
  };

References InitializeChannel(), MQwMockable::MQwMockable(), and SetADC18SaturationLimt().

Here is the call graph for this function:

◆ QwADC18_Channel() [3/4]

QwADC18_Channel::QwADC18_Channel ( const QwADC18_Channel & value )

inline

Definition at line 69 of file QwADC18_Channel.h.

                                               :
    VQwHardwareChannel(value), MQwMockable(value),
    fNumberOfSamples_map(value.fNumberOfSamples_map),
    fSaturationABSLimit(value.fSaturationABSLimit)
  {
    *this = value;
  };

References fNumberOfSamples_map, fSaturationABSLimit, MQwMockable::MQwMockable(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ QwADC18_Channel() [4/4]

QwADC18_Channel::QwADC18_Channel	(	const QwADC18_Channel &	value,
		VQwDataElement::EDataToSave	datatosave )

inline

Definition at line 76 of file QwADC18_Channel.h.

                                                                                     :
    VQwHardwareChannel(value,datatosave), MQwMockable(value),
    fNumberOfSamples_map(value.fNumberOfSamples_map),
    fSaturationABSLimit(value.fSaturationABSLimit)
  {
    *this = value;
  };

References fNumberOfSamples_map, fSaturationABSLimit, MQwMockable::MQwMockable(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ ~QwADC18_Channel()

QwADC18_Channel::~QwADC18_Channel ( )

inlineoverride

Definition at line 83 of file QwADC18_Channel.h.

83{ };

Member Function Documentation

◆ AccumulateRunningSum() [1/2]

void QwADC18_Channel::AccumulateRunningSum	(	const QwADC18_Channel &	value,
		Int_t	count = 0,
		Int_t	ErrorMask = 0xFFFFFFF )

Accumulate event values into the running sum with optional scaling.

Parameters

value	Source channel to accumulate from.
count	Event count scaling (0 means use value.fGoodEventCount).
ErrorMask	Bit mask of error flags to exclude when accumulating.

Moments and uncertainty calculation on the running sums and averages The calculation of the first and second moments of the running sum is not completely straightforward due to numerical instabilities associated with small variances and large average values. The naive computation taking the difference of the square of the average and the average of the squares leads to the subtraction of two very large numbers to get a small number.

Alternative algorithms (including for higher order moments) are supplied in Pebay, Philippe (2008), "Formulas for Robust, One-Pass Parallel Computation of Covariances and Arbitrary-Order Statistical Moments", Technical Report SAND2008-6212, Sandia National Laboratories. http://infoserve.sandia.gov/sand_doc/2008/086212.pdf

In the following formulas the moments \( M^1 \) and \( M^2 \) are defined

\begin{eqnarray} M^1 & = & \frac{1}{n} \sum^n y \\ M^2 & = & \sum^n (y - \mu) \end{eqnarray}

Recurrence relations for the addition of a single event:

\begin{eqnarray} M^1_n & = & M^1_{n-1} + \frac{y - M^1_{n-1}}{n} \\ M^2_n & = & M^2_{n-1} + (y - M^1_{n-1})(y - M^1_n) \end{eqnarray}

For the addition of an already accumulated sum:

\begin{eqnarray} M^1 & = & M^1_1 + n_2 \frac{M^1_2 - M^1_1}{n} \\ M^2 & = & M^2_1 + M^2_2 + n_1 n_2 \frac{(M^1_2 - M^1_1)^2}{n} \end{eqnarray}

In these recursive formulas we start from \( M^1 = y \) and \( M^2 = 0 \).

To calculate the mean and standard deviation we use

\begin{eqnarray} \mu & = & M^1 \\ \sigma^2 & = & \frac{1}{n} M^2 \end{eqnarray}

The standard deviation is a biased estimator, but this is what ROOT uses. Better would be to divide by \( (n-1) \).

We use the formulas provided there for the calculation of the first and second moments (i.e. average and variance).

Definition at line 1108 of file QwADC18_Channel.cc.

{
 
  if(count==0){
    count = value.fGoodEventCount;
  }
  
  // Moment calculations
  Int_t n1 = fGoodEventCount;
  Int_t n2 = count;
 
  // If there are no good events, check whether device HW is good
  if (n2 == 0 && value.fErrorFlag == 0) {
    n2 = 1;
  }
 
  if (ErrorMask ==  kPreserveError){
    //n = 1;
    if (n2 == 0) {
      n2 = 1;
    }
    if (count == -1) {
      n2 = -1;
    }
  }
 
  Int_t n = n1 + n2;
 
  // Set up variables
  Double_t M11 = fValue;
  Double_t M12 = value.fValue;
  Double_t M22 = value.fValueM2;
  if (n2 == 0) {
    // no good events for addition
    return;
  } else if (n2 == 1) {
    // simple version for addition of single event
    fGoodEventCount++;
    fValue += (M12 - M11) / n;
    fValueM2 += (M12 - M11) * (M12 - fValue); // note: using updated mean
  } else if (n2 > 1) {
    // general version for addition of multi-event sets
    fGoodEventCount += n2;
    fValue += n2 * (M12 - M11) / n;
    fValueM2 += M22 + n1 * n2 * (M12 - M11) * (M12 - M11) / n;
  }
 
  // Nanny
  if (fValue != fValue)
    QwWarning << "Angry Nanny: NaN detected in " << GetElementName() << QwLog::endl;
}

References QwLog::endl(), VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fValue, fValueM2, VQwDataElement::GetElementName(), kPreserveError, QwADC18_Channel(), and QwWarning.

Referenced by AccumulateRunningSum(), and DeaccumulateRunningSum().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ AccumulateRunningSum() [2/2]

void QwADC18_Channel::AccumulateRunningSum	(	const VQwHardwareChannel *	value,
		Int_t	count = 0,
		Int_t	ErrorMask = 0xFFFFFFF )

inlineoverridevirtual

Reimplemented from VQwHardwareChannel.

Definition at line 168 of file QwADC18_Channel.h.

                                                                                                               {
    const QwADC18_Channel *tmp_ptr = dynamic_cast<const QwADC18_Channel*>(value);
    if (tmp_ptr != NULL) {
      AccumulateRunningSum(*tmp_ptr, count, ErrorMask);
    } else {
      throw std::invalid_argument("Standard exception from QwADC18_Channel::AccumulateRunningSum: incompatible hardware channel type");
    }
  };

References AccumulateRunningSum(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ AddChannelOffset()

void QwADC18_Channel::AddChannelOffset ( Double_t offset )

This function will add a offset to the hw_sum and add the same offset for blocks.

Definition at line 1039 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty())
    {
      fValue += offset;
    }
}

References fValue, and VQwDataElement::IsNameEmpty().

Here is the call graph for this function:

◆ AddValueFrom()

void QwADC18_Channel::AddValueFrom ( const VQwHardwareChannel * valueptr )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 777 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(valueptr);
  if (tmpptr!=NULL){
    *this += *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::AddValueFrom = "
      +valueptr->GetElementName()+" is an incompatible type.";
    throw std::invalid_argument(loc.Data());
  }
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ ApplyHWChecks()

Int_t QwADC18_Channel::ApplyHWChecks ( )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 93 of file QwADC18_Channel.cc.

{
  Bool_t bStatus;
  if (bEVENTCUTMODE>0){//Global switch to ON/OFF event cuts set at the event cut file
 
    if (bDEBUG)
      QwWarning<<" QwQWVK_Channel "<<GetElementName()<<"  "<<GetNumberOfSamples()<<QwLog::endl;
 
    // Sample size check
    bStatus = MatchNumberOfSamples(fNumberOfSamples_map);//compare the default sample size with no.of samples read by the module
    if (!bStatus) {
      fErrorFlag |= kErrorFlag_sample;
    }
 
    //check sequence number
    fSequenceNo_Prev++;
    if (fSequenceNo_Counter==0 || GetSequenceNumber()==0){//starting the data run
      fSequenceNo_Prev=GetSequenceNumber();
    }
 
    if (!MatchSequenceNumber(fSequenceNo_Prev)){//we have a sequence number error
      fErrorFlag|=kErrorFlag_Sequence;
      if (bDEBUG)       QwWarning<<" QwQWVK_Channel "<<GetElementName()<<" Sequence number previous value = "<<fSequenceNo_Prev<<" Current value= "<< GetSequenceNumber()<<QwLog::endl;
    }
  }
  else {
    fGoodEventCount = 1;
    fErrorFlag = 0;
  }
 
  return fErrorFlag;
}

References bDEBUG, VQwHardwareChannel::bEVENTCUTMODE, QwLog::endl(), VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fNumberOfSamples_map, fSequenceNo_Counter, fSequenceNo_Prev, VQwDataElement::GetElementName(), GetNumberOfSamples(), GetSequenceNumber(), kErrorFlag_sample, kErrorFlag_Sequence, MatchNumberOfSamples(), MatchSequenceNumber(), and QwWarning.

Here is the call graph for this function:

◆ ApplySingleEventCuts() [1/2]

Bool_t QwADC18_Channel::ApplySingleEventCuts ( )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 1271 of file QwADC18_Channel.cc.

{
  Bool_t status;
 
  if (bEVENTCUTMODE>=2){//Global switch to ON/OFF event cuts set at the event cut file
 
    if (fULimit < fLLimit){
      status=kTRUE;
    } else  if (GetValue()<=fULimit && GetValue()>=fLLimit){
      if ((fErrorFlag)==0)
    status=kTRUE;
      else
    status=kFALSE;//If the device HW is failed
    }
    else{
      if (GetValue()> fULimit)
    fErrorFlag|=kErrorFlag_EventCut_U;
      else
    fErrorFlag|=kErrorFlag_EventCut_L;
      status=kFALSE;
    }
 
    if (bEVENTCUTMODE==3){
      status=kTRUE; //Update the event cut fail flag but pass the event.
    }
 
 
  }
  else{
    status=kTRUE;
    //fErrorFlag=0;//we need to keep the device error codes 
  }
 
  return status;
}

References VQwHardwareChannel::bEVENTCUTMODE, VQwDataElement::fErrorFlag, VQwHardwareChannel::fLLimit, VQwHardwareChannel::fULimit, GetValue(), kErrorFlag_EventCut_L, and kErrorFlag_EventCut_U.

Here is the call graph for this function:

◆ ApplySingleEventCuts() [2/2]

Bool_t QwADC18_Channel::ApplySingleEventCuts	(	Double_t	LL,
		Double_t	UL )

Definition at line 1255 of file QwADC18_Channel.cc.

{
  Bool_t status = kFALSE;
 
  if (UL < LL){
    status=kTRUE;
  } else  if (GetValue()<=UL && GetValue()>=LL){
    if ((fErrorFlag & kPreserveError)!=0)
      status=kTRUE;
    else
      status=kFALSE;//If the device HW is failed
  }
  std::cout<<(this->fErrorFlag & kPreserveError)<<std::endl;
  return status;
}

References VQwDataElement::fErrorFlag, GetValue(), and kPreserveError.

Here is the call graph for this function:

◆ ArcTan()

void QwADC18_Channel::ArcTan ( const QwADC18_Channel & value )

Definition at line 1013 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()) {
    this->fValue = atan(value.fValue);
  }
  // FIXME stuff missing here
}

References fValue, VQwDataElement::IsNameEmpty(), and QwADC18_Channel().

Here is the call graph for this function:

◆ AssignScaledValue()

void QwADC18_Channel::AssignScaledValue	(	const QwADC18_Channel &	value,
		Double_t	scale )

Definition at line 751 of file QwADC18_Channel.cc.

{
  if(this == &value) return;
 
  if (!IsNameEmpty()) {
    this->fValue      = value.fValue * scale;
    this->fValueError = value.fValueError;
    this->fValueM2    = value.fValueM2 * scale * scale;
    this->fErrorFlag  = value.fErrorFlag;//error code is updated.
    this->fGoodEventCount = value.fGoodEventCount;
  }
}

References VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fValue, fValueError, fValueM2, VQwDataElement::IsNameEmpty(), and QwADC18_Channel().

Here is the call graph for this function:

◆ AssignValueFrom()

void QwADC18_Channel::AssignValueFrom ( const VQwDataElement * valueptr )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 765 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(valueptr);
  if (tmpptr!=NULL){
    *this = *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::AssignValueFrom = "
      +valueptr->GetElementName()+" is an incompatible type.";
    throw std::invalid_argument(loc.Data());
  }
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwDataElement::VQwDataElement().

Here is the call graph for this function:

◆ Blind() [1/2]

void QwADC18_Channel::Blind ( const QwBlinder * blinder )

Blind this channel as an asymmetry.

Blind this channel as an asymmetry

Parameters

blinder Blinder

Definition at line 1199 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()) {
    if (blinder->IsBlinderOkay() && ((fErrorFlag)==0) ){
      blinder->BlindValue(fValue);
    } else {
      blinder->ModifyThisErrorCode(fErrorFlag);
      fValue = QwBlinder::kValue_BlinderFail;
    }
  }
}

References QwBlinder::BlindValue(), VQwDataElement::fErrorFlag, fValue, QwBlinder::IsBlinderOkay(), VQwDataElement::IsNameEmpty(), QwBlinder::kValue_BlinderFail, and QwBlinder::ModifyThisErrorCode().

Here is the call graph for this function:

◆ Blind() [2/2]

void QwADC18_Channel::Blind	(	const QwBlinder *	blinder,
		const QwADC18_Channel &	yield )

Blind this channel as a difference.

Blind this channel as a difference with specified yield

Parameters

blinder	Blinder
yield	Corresponding yield

Definition at line 1216 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()) {
    if (blinder->IsBlinderOkay() && ((fErrorFlag) ==0) ){
      blinder->BlindValue(fValue, yield.fValue);
    } else {
      blinder->ModifyThisErrorCode(fErrorFlag);//update the HW error code
      fValue = QwBlinder::kValue_BlinderFail * yield.fValue;
    }
  }
}

References QwBlinder::BlindValue(), VQwDataElement::fErrorFlag, fValue, QwBlinder::IsBlinderOkay(), VQwDataElement::IsNameEmpty(), QwBlinder::kValue_BlinderFail, QwBlinder::ModifyThisErrorCode(), and QwADC18_Channel().

Here is the call graph for this function:

◆ CalculateRunningAverage()

void QwADC18_Channel::CalculateRunningAverage ( )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 1161 of file QwADC18_Channel.cc.

{
  //  See notes in QwVQWK_Channel;  we are using:
  //         error = sqrt(M2)/n,
  //  or alternately we could use the unbiased estimator for both
  //  the sigma and error on the mean:
  //          error = sqrt(M2)/(n-1)
  if(fGoodEventCount > 0) {
    fValueError = sqrt(fValueM2) / fGoodEventCount;
  } else {
    fValueError = 0.0;
  }
}

References VQwDataElement::fGoodEventCount, fValueError, and fValueM2.

◆ ClearEventData()

void QwADC18_Channel::ClearEventData ( )

overridevirtual

Clear the event data in this element.

Reimplemented from VQwDataElement.

Definition at line 234 of file QwADC18_Channel.cc.

{
  fDiff_Raw    = 0;
  fBase_Raw    = 0;
  fPeak_Raw    = 0;
  fValue       = 0.0;
  fValueM2     = 0.0;
  fValueError  = 0.0;
  fSequenceNumber   = 0;
  fNumberOfSamples  = 0;
  fGoodEventCount   = 0;
  fErrorFlag = 0;
}

References fBase_Raw, fDiff_Raw, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fNumberOfSamples, fPeak_Raw, fSequenceNumber, fValue, fValueError, and fValueM2.

Referenced by InitializeChannel().

Here is the caller graph for this function:

◆ Clone() [1/2]

virtual VQwHardwareChannel * VQwHardwareChannel::Clone ( ) const

inlinevirtual

Reimplemented from VQwHardwareChannel.

Definition at line 112 of file VQwHardwareChannel.h.

                                           {
    return Clone(this->fDataToSave);
  };

◆ Clone() [2/2]

VQwHardwareChannel * QwADC18_Channel::Clone ( VQwDataElement::EDataToSave datatosave ) const

inlineoverridevirtual

Implements VQwHardwareChannel.

Definition at line 87 of file QwADC18_Channel.h.

                                                                                {
    return new QwADC18_Channel(*this,datatosave);
  };

References QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ ConstructBranch()

void QwADC18_Channel::ConstructBranch	(	TTree *	tree,
		TString &	prefix )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 603 of file QwADC18_Channel.cc.

{
  if (IsNameEmpty()){
    //  This channel is not used, so skip setting up the tree.
  } else {
    TString basename = prefix + GetElementName();
    tree->Branch(basename, &fValue, basename+"/D");
  }
}

References fValue, VQwDataElement::GetElementName(), and VQwDataElement::IsNameEmpty().

Here is the call graph for this function:

◆ ConstructBranchAndVector()

void QwADC18_Channel::ConstructBranchAndVector	(	TTree *	tree,
		TString &	prefix,
		QwRootTreeBranchVector &	values )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 572 of file QwADC18_Channel.cc.

{
  if (IsNameEmpty()){
    //  This channel is not used, so skip setting up the tree.
  } else {
    //  Decide what to store based on prefix
    SetDataToSaveByPrefix(prefix);
 
    TString basename = prefix(0, (prefix.First("|") >= 0)? prefix.First("|"): prefix.Length()) + GetElementName();
    fTreeArrayIndex  = values.size();
 
    values.push_back("value", 'D');
    if (fDataToSave == kMoments) {
      values.push_back("value_m2", 'D');
      values.push_back("value_err", 'D');
    }
 
    values.push_back("Device_Error_Code", 'i');
    if (fDataToSave == kRaw){
      values.push_back("raw", 'I');
      values.push_back("diff", 'I');
      values.push_back("peak", 'I');
      values.push_back("base", 'I');
    }
 
    fTreeArrayNumEntries = values.size() - fTreeArrayIndex;
    if (gQwHists.MatchDeviceParamsFromList(basename.Data()))
      tree->Branch(basename, &(values[fTreeArrayIndex]), values.LeafList(fTreeArrayIndex).c_str());
  }
}

References VQwHardwareChannel::fDataToSave, VQwHardwareChannel::fTreeArrayIndex, VQwHardwareChannel::fTreeArrayNumEntries, VQwDataElement::GetElementName(), gQwHists, VQwDataElement::IsNameEmpty(), VQwDataElement::kMoments, VQwDataElement::kRaw, QwRootTreeBranchVector::LeafList(), QwRootTreeBranchVector::push_back(), VQwHardwareChannel::SetDataToSaveByPrefix(), and QwRootTreeBranchVector::size().

Here is the call graph for this function:

◆ ConstructHistograms()

void QwADC18_Channel::ConstructHistograms	(	TDirectory *	folder,
		TString &	prefix )

overridevirtual

Construct the histograms for this data element.

Implements VQwDataElement.

Definition at line 511 of file QwADC18_Channel.cc.

{
  //  If we have defined a subdirectory in the ROOT file, then change into it.
  if (folder != NULL) folder->cd();
 
  if (IsNameEmpty()){
    //  This channel is not used, so skip filling the histograms.
  } else {
    //  Now create the histograms.
    if (prefix == TString("asym_")
       || prefix == TString("diff_")
       || prefix == TString("yield_"))
      fDataToSave=kDerived;
 
    TString basename, fullname;
    basename = prefix + GetElementName();
 
    if(fDataToSave==kRaw)
      {
    fHistograms.resize(2, NULL);
    size_t index=0;
    fHistograms[index++] = gQwHists.Construct1DHist(basename);
    fHistograms[index++] = gQwHists.Construct1DHist(basename+Form("_raw"));
      }
    else if(fDataToSave==kDerived)
      {
    fHistograms.resize(1, NULL);
    Int_t index=0;
    fHistograms[index++] = gQwHists.Construct1DHist(basename);
      }
    else
      {
    // this is not recognized
      }
  }
}

References VQwHardwareChannel::fDataToSave, MQwHistograms::fHistograms, VQwDataElement::GetElementName(), gQwHists, VQwDataElement::IsNameEmpty(), VQwDataElement::kDerived, and VQwDataElement::kRaw.

Here is the call graph for this function:

◆ DeaccumulateRunningSum() [1/2]

void QwADC18_Channel::DeaccumulateRunningSum	(	const QwADC18_Channel &	value,
		Int_t	ErrorMask = 0xFFFFFFF )

inline

Definition at line 177 of file QwADC18_Channel.h.

                                                                                             {
    AccumulateRunningSum(value, -1, ErrorMask);
  };

References AccumulateRunningSum(), and QwADC18_Channel().

Here is the call graph for this function:

◆ DeaccumulateRunningSum() [2/2]

virtual void VQwHardwareChannel::DeaccumulateRunningSum	(	const VQwHardwareChannel *	value,
		Int_t	ErrorMask = 0xFFFFFFF )

inlinevirtual

Reimplemented from VQwHardwareChannel.

Definition at line 255 of file VQwHardwareChannel.h.

                                                                                                 {
    AccumulateRunningSum(value, -1, ErrorMask);
  };

◆ Difference()

void QwADC18_Channel::Difference	(	const QwADC18_Channel &	value1,
		const QwADC18_Channel &	value2 )

Definition at line 948 of file QwADC18_Channel.cc.

{
  *this  = value1;
  *this -= value2;
}

References QwADC18_Channel().

Here is the call graph for this function:

◆ DivideBy() [1/2]

void QwADC18_Channel::DivideBy ( const QwADC18_Channel & denom )

Definition at line 1056 of file QwADC18_Channel.cc.

{
  *this /= denom;
}

References QwADC18_Channel().

Here is the call graph for this function:

◆ DivideBy() [2/2]

void QwADC18_Channel::DivideBy ( const VQwHardwareChannel * valueptr )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 813 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(valueptr);
  if (tmpptr!=NULL){
    *this /= *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::DivideBy = "
      +valueptr->GetElementName()+" is an incompatible type.";
    throw std::invalid_argument(loc.Data());
  }
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ EncodeEventData()

void QwADC18_Channel::EncodeEventData ( std::vector< UInt_t > & buffer )

overridevirtual

Encode the event data into a CODA buffer.

Implements MQwMockable.

Definition at line 295 of file QwADC18_Channel.cc.

{
  UInt_t localbuf[kDataWordsPerChannel] = {0};
 
  if (IsNameEmpty()) {
    //  This channel is not used, but is present in the data stream.
    //  Fill in with zero.
    localbuf[0] = 0;
  } else {
    localbuf[0] |=   mask31x & (0 << 31);
    localbuf[0] |= mask3029x & (0 << 29);
    localbuf[0] |= mask2625x & (0 << 25);
    UInt_t exp_dtype = 1;
    switch (exp_dtype) {
      case 0:
        break;
      case 1:
      case 2:
        localbuf[0] |= mask2422x & (2 << 22);
        localbuf[0] |=  mask170x & fPeak_Raw;
        break;
      case 4:
        break;
      default:
        QwError << "QwADC18_Channel::EncodeEventData: Unknown data type" << QwLog::endl;
    }
  }
 
  for (Int_t i = 0; i < kDataWordsPerChannel; i++) {
    buffer.push_back(localbuf[i]);
  }
}

References QwLog::endl(), fPeak_Raw, VQwDataElement::IsNameEmpty(), kDataWordsPerChannel, mask170x, mask2422x, mask2625x, mask3029x, mask31x, and QwError.

Here is the call graph for this function:

◆ FillHistograms()

void QwADC18_Channel::FillHistograms ( )

overridevirtual

Fill the histograms for this data element.

Implements VQwDataElement.

Definition at line 548 of file QwADC18_Channel.cc.

{
  Int_t index=0;
 
  if (IsNameEmpty())
    {
      //  This channel is not used, so skip creating the histograms.
    } else
      {
    if(fDataToSave==kRaw)
      {
        if (fHistograms[index] != NULL && (fErrorFlag)==0)
          fHistograms[index++]->Fill(GetValue());
        if (fHistograms[index] != NULL && (fErrorFlag)==0)
          fHistograms[index++]->Fill(GetRawValue());
      }
    else if(fDataToSave==kDerived)
      {
        if (fHistograms[index] != NULL && (fErrorFlag)==0)
          fHistograms[index++]->Fill(GetValue());
      }
    }
}

References VQwHardwareChannel::fDataToSave, VQwDataElement::fErrorFlag, MQwHistograms::fHistograms, GetRawValue(), GetValue(), VQwDataElement::IsNameEmpty(), VQwDataElement::kDerived, and VQwDataElement::kRaw.

Here is the call graph for this function:

◆ FillTreeVector()

void QwADC18_Channel::FillTreeVector ( QwRootTreeBranchVector & values ) const

overridevirtual

Implements VQwHardwareChannel.

Definition at line 613 of file QwADC18_Channel.cc.

{
  if (IsNameEmpty()) {
    //  This channel is not used, so skip setting up the tree.
  } else if (fTreeArrayNumEntries == 0) {
    static bool warned = false;
    if (!warned) {
      QwError << "QwADC18_Channel::FillTreeVector:  fTreeArrayNumEntries=="
              << fTreeArrayNumEntries << " (no branch constructed?)" << QwLog::endl;
      QwError << "Offending element is " << GetElementName() << QwLog::endl;
      warned = true;
    }
  } else if (values.size() < fTreeArrayIndex+fTreeArrayNumEntries) {
    QwError << "QwADC18_Channel::FillTreeVector:  values.size()=="
            << values.size() << " name: " << fElementName
            << "; fTreeArrayIndex+fTreeArrayNumEntries=="
            << fTreeArrayIndex << '+' << fTreeArrayNumEntries << '='
            << fTreeArrayIndex+fTreeArrayNumEntries
            << QwLog::endl;
  } else {
    size_t index = fTreeArrayIndex;
    values.SetValue(index++, this->fValue);
    if (fDataToSave == kMoments) {
      values.SetValue(index++, fValueM2);
      values.SetValue(index++, fValueError);
    }
    values.SetValue(index++, this->fErrorFlag);
    if(fDataToSave==kRaw){
      values.SetValue(index++, this->fValue_Raw);
      values.SetValue(index++, this->fDiff_Raw);
      values.SetValue(index++, this->fPeak_Raw);
      values.SetValue(index++, this->fBase_Raw);
    }
  }
}

References QwLog::endl(), fBase_Raw, VQwHardwareChannel::fDataToSave, fDiff_Raw, VQwDataElement::fElementName, VQwDataElement::fErrorFlag, fPeak_Raw, VQwHardwareChannel::fTreeArrayIndex, VQwHardwareChannel::fTreeArrayNumEntries, fValue, fValue_Raw, fValueError, fValueM2, VQwDataElement::GetElementName(), VQwDataElement::IsNameEmpty(), VQwDataElement::kMoments, VQwDataElement::kRaw, QwError, QwRootTreeBranchVector::SetValue(), and QwRootTreeBranchVector::size().

Here is the call graph for this function:

◆ ForceMapfileSampleSize()

void QwADC18_Channel::ForceMapfileSampleSize ( )

inline

Forces the event "number of samples" variable to be what was expected from the mapfile. NOTE: this should only be used in mock data generation!

Definition at line 112 of file QwADC18_Channel.h.

112{fNumberOfSamples = fNumberOfSamples_map;};

References fNumberOfSamples, and fNumberOfSamples_map.

◆ GetADC18SaturationLimt()

Double_t QwADC18_Channel::GetADC18SaturationLimt ( )

inline

Definition at line 210 of file QwADC18_Channel.h.

                                   {
    return fSaturationABSLimit;
  }

References fSaturationABSLimit.

◆ GetAverageVolts()

Double_t QwADC18_Channel::GetAverageVolts ( ) const

Definition at line 490 of file QwADC18_Channel.cc.

{
  return fValue * kADC18_VoltsPerBit / fNumberOfSamples;
}

References fNumberOfSamples, fValue, and kADC18_VoltsPerBit.

◆ GetBufferOffset()

Int_t QwADC18_Channel::GetBufferOffset	(	Int_t	moduleindex,
		Int_t	channelindex )

static

Class: QwADC18_Channel Base class containing decoding functions for the HAPPEX 18-bit ADC The functions in this class will decode a single channel worth of ADC18_Channel data and provide the components through member functions.

Static member function to return the word offset within a data buffer given the module number index and the channel number index.

Parameters

moduleindex	Module index within this buffer; counts from zero
channelindex	Channel index within this module; counts from -1 (header)

Returns: The number of words offset to the beginning of this channel's data from the beginning of the ADC18 buffer.

Definition at line 62 of file QwADC18_Channel.cc.

{
  Int_t offset = -1;
  if (moduleindex<0 ){
    QwError << "QwADC18_Channel::GetBufferOffset:  Invalid module index,"
            << moduleindex
            << ".  Must be zero or greater."
            << QwLog::endl;
  } else if (channelindex<-1 || channelindex>kMaxChannels){
    QwError << "QwADC18_Channel::GetBufferOffset:  Invalid channel index,"
            << channelindex
            << ".  Must be in range [-1," << kMaxChannels << "]."
            << QwLog::endl;
  } else {
    if (channelindex == -1) {
      // Header
      offset = kHeaderWordsPerBank
             + moduleindex * (kHeaderWordsPerModule + kMaxChannels * kDataWordsPerChannel + kFooterWordsPerModule);
    } else {
      // Data word
      offset = kHeaderWordsPerBank
             + moduleindex * (kHeaderWordsPerModule + kMaxChannels * kDataWordsPerChannel + kFooterWordsPerModule)
             + kHeaderWordsPerModule
             + channelindex * kDataWordsPerChannel;
    }
  }
  return offset;
}

References QwLog::endl(), kDataWordsPerChannel, kFooterWordsPerModule, kHeaderWordsPerBank, kHeaderWordsPerModule, kMaxChannels, and QwError.

Referenced by QwBeamDetectorID::QwBeamDetectorID().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetNumberOfSamples()

size_t QwADC18_Channel::GetNumberOfSamples ( ) const

inline

Definition at line 242 of file QwADC18_Channel.h.

242{return (fNumberOfSamples);};

References fNumberOfSamples.

Referenced by ApplyHWChecks().

Here is the caller graph for this function:

◆ GetRawValue() [1/2]

Int_t VQwHardwareChannel::GetRawValue ( ) const

inline

Definition at line 125 of file VQwHardwareChannel.h.

125{return this->GetRawValue(0);};

Referenced by FillHistograms().

Here is the caller graph for this function:

◆ GetRawValue() [2/2]

Int_t QwADC18_Channel::GetRawValue ( size_t element ) const

inlineoverridevirtual

Implements VQwHardwareChannel.

Definition at line 223 of file QwADC18_Channel.h.

223{ return fValue_Raw; };

References fValue_Raw.

◆ GetSequenceNumber()

size_t QwADC18_Channel::GetSequenceNumber ( ) const

inline

Definition at line 241 of file QwADC18_Channel.h.

241{return (fSequenceNumber);};

References fSequenceNumber.

Referenced by ApplyHWChecks().

Here is the caller graph for this function:

◆ GetValue() [1/2]

Double_t VQwHardwareChannel::GetValue ( ) const

inline

Definition at line 126 of file VQwHardwareChannel.h.

126{return this->GetValue(0);};

Referenced by ApplySingleEventCuts(), ApplySingleEventCuts(), FillHistograms(), and PrintValue().

Here is the caller graph for this function:

◆ GetValue() [2/2]

Double_t QwADC18_Channel::GetValue ( size_t element ) const

inlineoverridevirtual

Implements VQwHardwareChannel.

Definition at line 224 of file QwADC18_Channel.h.

224{ return fValue; };

References fValue.

Referenced by operator<<.

Here is the caller graph for this function:

◆ GetValueError() [1/2]

Double_t VQwHardwareChannel::GetValueError ( ) const

inline

Definition at line 128 of file VQwHardwareChannel.h.

128{return this->GetValueError(0);};

QwADC18_Channel::GetValueError

Double_t GetValueError() const

Definition VQwHardwareChannel.h:128

Referenced by PrintValue().

Here is the caller graph for this function:

◆ GetValueError() [2/2]

Double_t QwADC18_Channel::GetValueError ( size_t element ) const

inlineoverridevirtual

Implements VQwHardwareChannel.

Definition at line 226 of file QwADC18_Channel.h.

226{ return fValueError; };

References fValueError.

◆ GetValueM2() [1/2]

Double_t VQwHardwareChannel::GetValueM2 ( ) const

inline

Definition at line 127 of file VQwHardwareChannel.h.

127{return this->GetValueM2(0);};

QwADC18_Channel::GetValueM2

Double_t GetValueM2() const

Definition VQwHardwareChannel.h:127

◆ GetValueM2() [2/2]

Double_t QwADC18_Channel::GetValueM2 ( size_t element ) const

inlineoverridevirtual

Implements VQwHardwareChannel.

Definition at line 225 of file QwADC18_Channel.h.

225{ return fValueM2; };

References fValueM2.

◆ GetValueWidth() [1/2]

Double_t VQwHardwareChannel::GetValueWidth ( ) const

inline

Definition at line 129 of file VQwHardwareChannel.h.

129{return this->GetValueWidth(0);};

QwADC18_Channel::GetValueWidth

Double_t GetValueWidth() const

Definition VQwHardwareChannel.h:129

Referenced by PrintValue().

Here is the caller graph for this function:

◆ GetValueWidth() [2/2]

Double_t VQwHardwareChannel::GetValueWidth ( size_t element ) const

inline

Definition at line 134 of file VQwHardwareChannel.h.

                                               {
    RangeCheck(element);
    Double_t width;
    if (fGoodEventCount>0){
      width = (GetValueError(element)*std::sqrt(Double_t(fGoodEventCount))); 
    } else {
      width = 0.0;
    } 
    return width;
  };

◆ IncrementErrorCounters()

void QwADC18_Channel::IncrementErrorCounters ( )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 128 of file QwADC18_Channel.cc.

                                            {
  if ( (kErrorFlag_sample &  fErrorFlag)==kErrorFlag_sample)
    fErrorCount_sample++; //increment the hw error counter
  if ( (kErrorFlag_SW_HW &  fErrorFlag)==kErrorFlag_SW_HW)
    fErrorCount_SW_HW++; //increment the hw error counter
  if ( (kErrorFlag_Sequence &  fErrorFlag)==kErrorFlag_Sequence)
    fErrorCount_Sequence++; //increment the hw error counter
  if ( (kErrorFlag_SameHW &  fErrorFlag)==kErrorFlag_SameHW)
    fErrorCount_SameHW++; //increment the hw error counter
  if ( (kErrorFlag_ZeroHW &  fErrorFlag)==kErrorFlag_ZeroHW)
    fErrorCount_ZeroHW++; //increment the hw error counter
  if ( ((kErrorFlag_EventCut_L &  fErrorFlag)==kErrorFlag_EventCut_L) 
       || ((kErrorFlag_EventCut_U &  fErrorFlag)==kErrorFlag_EventCut_U)){
    fNumEvtsWithEventCutsRejected++; //increment the event cut error counter
  }
}

References fErrorCount_SameHW, fErrorCount_sample, fErrorCount_Sequence, fErrorCount_SW_HW, fErrorCount_ZeroHW, VQwDataElement::fErrorFlag, fNumEvtsWithEventCutsRejected, kErrorFlag_EventCut_L, kErrorFlag_EventCut_U, kErrorFlag_SameHW, kErrorFlag_sample, kErrorFlag_Sequence, kErrorFlag_SW_HW, and kErrorFlag_ZeroHW.

◆ InitializeChannel() [1/2]

void QwADC18_Channel::InitializeChannel	(	TString	name,
		TString	datatosave )

overridevirtual

Initialize the fields in this object.

Implements VQwHardwareChannel.

Definition at line 147 of file QwADC18_Channel.cc.

{
  SetElementName(name);
  SetDataToSave(datatosave);
  SetNumberOfDataWords(1);
  SetNumberOfSubElements(1);
 
  kFoundPedestal = 0;
  kFoundGain = 0;
 
  fPedestal            = 0.0;
  fCalibrationFactor   = 1.0;
 
  fTreeArrayIndex      = 0;
  fTreeArrayNumEntries = 0;
 
  ClearEventData();
 
  fPreviousSequenceNumber = 0;
  fNumberOfSamples_map    = 0;
  fNumberOfSamples        = 0;
 
  // Use internal random variable by default
  fUseExternalRandomVariable = false;
 
  // Mock drifts
  fMockDriftAmplitude.clear();
  fMockDriftFrequency.clear();
  fMockDriftPhase.clear();
 
  // Mock asymmetries
  fMockAsymmetry     = 0.0;
  fMockGaussianMean  = 0.0;
  fMockGaussianSigma = 0.0;
 
  // Event cuts
  fULimit=-1;
  fLLimit=1;
  fNumEvtsWithEventCutsRejected = 0;
 
  fErrorFlag=0;               //Initialize the error flag
  fErrorConfigFlag=0;         //Initialize the error config. flag
 
  //init error counters//
  fErrorCount_sample     = 0;
  fErrorCount_SW_HW      = 0;
  fErrorCount_Sequence   = 0;
  fErrorCount_SameHW     = 0;
  fErrorCount_ZeroHW     = 0;
  fErrorCount_HWSat      = 0;
 
  fRunningSum            = 0;
 
  fADC_Same_NumEvt       = 0;
  fSequenceNo_Prev       = 0;
  fSequenceNo_Counter    = 0;
 
  fGoodEventCount        = 0;
 
  bEVENTCUTMODE          = 0;
 
  //std::cout<< "name = "<<name<<" error count same _HW = "<<fErrorCount_SameHW <<std::endl;
  return;
}

Referenced by InitializeChannel(), QwADC18_Channel(), and QwADC18_Channel().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ InitializeChannel() [2/2]

void QwADC18_Channel::InitializeChannel	(	TString	subsystem,
		TString	instrumenttype,
		TString	name,
		TString	datatosave )

overridevirtual

Initialize the fields in this object.

Implements VQwHardwareChannel.

Definition at line 214 of file QwADC18_Channel.cc.

                                                                                                                  {
  InitializeChannel(name,datatosave);
  SetSubsystemName(subsystem);
  SetModuleType(instrumenttype);
  //PrintInfo();
}

References InitializeChannel(), VQwDataElement::SetModuleType(), and VQwDataElement::SetSubsystemName().

Here is the call graph for this function:

◆ IsHeaderWord()

Bool_t QwADC18_Channel::IsHeaderWord ( UInt_t word ) const

Decode the event data from a CODA buffer.

Definition at line 328 of file QwADC18_Channel.cc.

{
  return ((rawd & mask31x) != 0);
}

References mask31x.

Referenced by ProcessEvBuffer().

Here is the caller graph for this function:

◆ LoadChannelParameters()

void QwADC18_Channel::LoadChannelParameters ( QwParameterFile & paramfile )

overridevirtual

Reimplemented from VQwDataElement.

Definition at line 221 of file QwADC18_Channel.cc.

                                                                     {
  UInt_t value = 0;
  if (paramfile.ReturnValue("sample_size",value)){
    SetDefaultSampleSize(value);
  } else {
    QwWarning << "ADC18 Channel "
          << GetElementName()
          << " cannot set the default sample size."
          << QwLog::endl;
  }
};

References QwLog::endl(), VQwDataElement::GetElementName(), QwWarning, QwParameterFile::ReturnValue(), and SetDefaultSampleSize().

Here is the call graph for this function:

◆ LoadMockDataParameters()

void MQwMockable::LoadMockDataParameters ( QwParameterFile & paramfile )

virtual

Load the mock data parameters from the current line in the param file.

Reimplemented from VQwDataElement.

Definition at line 60 of file MQwMockable.cc.

                                                                  {
  Bool_t   ldebug=kFALSE;
  Double_t asym=0.0, mean=0.0, sigma=0.0;
  Double_t amplitude=0.0, phase=0.0, frequency=0.0;
 
  //Check to see if this line contains "drift"
  if (paramfile.GetLine().find("drift")!=std::string::npos){
    //  "drift" appears somewhere in the line.  Assume it is the next token and move on...
    paramfile.GetNextToken();  //Throw away this token.  Now the next three should be the drift parameters.
    //read 3 parameters
    amplitude = paramfile.GetTypedNextToken<Double_t>();
    phase     = paramfile.GetTypedNextToken<Double_t>();
    frequency = paramfile.GetTypedNextToken<Double_t>();  //  The frequency we read in should be in Hz.
    this->AddRandomEventDriftParameters(amplitude, phase, frequency*Qw::Hz);
    // std::cout << "In MQwMockable::LoadMockDataParameters: amp = " << amplitude << "\t phase = " << phase << "\t freq = " << frequency << std::endl;
  } 
  else {
    asym  = paramfile.GetTypedNextToken<Double_t>();
    mean  = paramfile.GetTypedNextToken<Double_t>(); 
    sigma = paramfile.GetTypedNextToken<Double_t>();      
    if (ldebug==1) {
      std::cout << "#################### \n";
      std::cout << "asym, mean, sigma \n" << std::endl;
      std::cout << asym                   << " / "
            << mean                   << " / "
            << sigma                  << " / "
                << std::endl;  
    }
    this->SetRandomEventParameters(mean, sigma);
    this->SetRandomEventAsymmetry(asym);
  }
}

References MQwMockable::fMockDriftAmplitude, MQwMockable::fMockDriftFrequency, and MQwMockable::fMockDriftPhase.

◆ MatchNumberOfSamples()

Bool_t QwADC18_Channel::MatchNumberOfSamples ( size_t numsamp )

Definition at line 1237 of file QwADC18_Channel.cc.

{
  Bool_t status = kTRUE;
  if (!IsNameEmpty()){
    status = (fNumberOfSamples==numsamp);
    if (! status){
      if (bDEBUG)
    QwError << "QwADC18_Channel::MatchNumberOfSamples:  Channel "
        << GetElementName()
        << " had fNumberOfSamples==" << fNumberOfSamples
        << " and was supposed to have " << numsamp
        << std::endl;
    }
  }
  return status;
}

References bDEBUG, fNumberOfSamples, VQwDataElement::GetElementName(), VQwDataElement::IsNameEmpty(), and QwError.

Referenced by ApplyHWChecks().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ MatchSequenceNumber()

Bool_t QwADC18_Channel::MatchSequenceNumber ( size_t seqnum )

Definition at line 1228 of file QwADC18_Channel.cc.

{
  Bool_t status = kTRUE;
  if (!IsNameEmpty()){
    status = (fSequenceNumber == seqnum);
    }
  return status;
}

References fSequenceNumber, and VQwDataElement::IsNameEmpty().

Referenced by ApplyHWChecks().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ MultiplyBy()

void QwADC18_Channel::MultiplyBy ( const VQwHardwareChannel * valueptr )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 801 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(valueptr);
  if (tmpptr!=NULL){
    *this *= *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::MultiplyBy = "
      +valueptr->GetElementName()+" is an incompatible type.";
    throw std::invalid_argument(loc.Data());
  }
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ operator*()

const QwADC18_Channel QwADC18_Channel::operator* ( const QwADC18_Channel & value ) const

Definition at line 861 of file QwADC18_Channel.cc.

{
  QwADC18_Channel result = *this;
  result *= value;
  return result;
}

References QwADC18_Channel().

Here is the call graph for this function:

◆ operator*=() [1/2]

QwADC18_Channel & QwADC18_Channel::operator*= ( const QwADC18_Channel & value )

Definition at line 868 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()){
    this->fValue     *= value.fValue;
    this->fDiff_Raw   = 0;
    this->fPeak_Raw   = 0;
    this->fBase_Raw   = 0;
    this->fValueM2    = 0.0;
    this->fErrorFlag |= (value.fErrorFlag);//error code is ORed.
  }
  return *this;
}

References fBase_Raw, fDiff_Raw, VQwDataElement::fErrorFlag, fPeak_Raw, fValue, fValueM2, VQwDataElement::IsNameEmpty(), and QwADC18_Channel().

Here is the call graph for this function:

◆ operator*=() [2/2]

VQwHardwareChannel & QwADC18_Channel::operator*= ( const VQwHardwareChannel & input )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 911 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(&source);
  if (tmpptr!=NULL){
    *this *= *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::operator*= "
        +source.GetElementName()+" "
        +this->GetElementName()+" are not of the same type";
    throw(std::invalid_argument(loc.Data()));
  }
  return *this;
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ operator+()

const QwADC18_Channel QwADC18_Channel::operator+ ( const QwADC18_Channel & value ) const

Definition at line 827 of file QwADC18_Channel.cc.

{
  QwADC18_Channel result = *this;
  result += value;
  return result;
}

References QwADC18_Channel().

Here is the call graph for this function:

◆ operator+=() [1/2]

QwADC18_Channel & QwADC18_Channel::operator+= ( const QwADC18_Channel & value )

Definition at line 834 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()) {
    this->fValue           += value.fValue;
    this->fValueM2          = 0.0;
    this->fErrorFlag       |= value.fErrorFlag;//error code is ORed.
  }
  return *this;
}

References VQwDataElement::fErrorFlag, fValue, fValueM2, VQwDataElement::IsNameEmpty(), and QwADC18_Channel().

Here is the call graph for this function:

◆ operator+=() [2/2]

VQwHardwareChannel & QwADC18_Channel::operator+= ( const VQwHardwareChannel & input )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 881 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(&source);
  if (tmpptr!=NULL){
    *this += *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::operator+= "
        +source.GetElementName()+" "
        +this->GetElementName()+" are not of the same type";
    throw(std::invalid_argument(loc.Data()));
  }
  return *this;
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ operator-()

const QwADC18_Channel QwADC18_Channel::operator- ( const QwADC18_Channel & value ) const

Definition at line 844 of file QwADC18_Channel.cc.

{
  QwADC18_Channel result = *this;
  result -= value;
  return result;
}

References QwADC18_Channel().

Here is the call graph for this function:

◆ operator-=() [1/2]

QwADC18_Channel & QwADC18_Channel::operator-= ( const QwADC18_Channel & value )

Definition at line 851 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()){
    this->fValue           -= value.fValue;
    this->fValueM2          = 0.0;
    this->fErrorFlag       |= (value.fErrorFlag);//error code is ORed.
  }
  return *this;
}

References VQwDataElement::fErrorFlag, fValue, fValueM2, VQwDataElement::IsNameEmpty(), and QwADC18_Channel().

Here is the call graph for this function:

◆ operator-=() [2/2]

VQwHardwareChannel & QwADC18_Channel::operator-= ( const VQwHardwareChannel & input )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 896 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(&source);
  if (tmpptr!=NULL){
    *this -= *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::operator-= "
        +source.GetElementName()+" "
        +this->GetElementName()+" are not of the same type";
    throw(std::invalid_argument(loc.Data()));
  }
  return *this;
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ operator/=() [1/2]

QwADC18_Channel & QwADC18_Channel::operator/= ( const QwADC18_Channel & value )

protected

Definition at line 971 of file QwADC18_Channel.cc.

{
  //  In this function, leave the "raw" variables untouched.
  Double_t ratio;
  Double_t variance;
  if (!IsNameEmpty()) {
    // The variances are calculated using the following formula:
    //   Var[ratio] = ratio^2 (Var[numer] / numer^2 + Var[denom] / denom^2)
    //
    // This requires that both the numerator and denominator are non-zero!
    //
    if (this->fValue != 0.0 && denom.fValue != 0.0){
      ratio = (this->fValue) / (denom.fValue);
      variance =  ratio * ratio *
        (this->fValueM2 / this->fValue / this->fValue
         + denom.fValueM2 / denom.fValue / denom.fValue);
      fValue   = ratio;
      fValueM2 = variance;
    } else if (this->fValue == 0.0) {
      fValue   = 0.0;
      fValueM2 = 0.0;
    } else {
      QwVerbose << "Attempting to divide by zero in " 
                << GetElementName() << QwLog::endl;
      fValue   = 0.0;
      fValueM2 = 0.0;
    }
 
    // Remaining variables
    //  Don't change fGoodEventCount.
    //  'OR' the device error codes together.
    fErrorFlag |= denom.fErrorFlag;
  }
 
  // Nanny
  if (fValue != fValue)
    QwWarning << "Angry Nanny: NaN detected in " << GetElementName() << QwLog::endl;
  return *this;
}

References QwLog::endl(), VQwDataElement::fErrorFlag, fValue, fValueM2, VQwDataElement::GetElementName(), VQwDataElement::IsNameEmpty(), QwADC18_Channel(), QwVerbose, and QwWarning.

Here is the call graph for this function:

◆ operator/=() [2/2]

VQwHardwareChannel & QwADC18_Channel::operator/= ( const VQwHardwareChannel & input )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 926 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(&source);
  if (tmpptr!=NULL){
    *this /= *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::operator/= "
        +source.GetElementName()+" "
        +this->GetElementName()+" are not of the same type";
    throw(std::invalid_argument(loc.Data()));
  }
  return *this;
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ operator=()

QwADC18_Channel & QwADC18_Channel::operator= ( const QwADC18_Channel & value )

Definition at line 734 of file QwADC18_Channel.cc.

{
  if (this == &value) return *this;
 
  if (!IsNameEmpty()) {
    VQwHardwareChannel::operator=(value);
    this->fDiff_Raw   = value.fDiff_Raw;
    this->fPeak_Raw   = value.fPeak_Raw;
    this->fBase_Raw   = value.fBase_Raw;
    this->fValue_Raw  = value.fValue_Raw;
    this->fValue      = value.fValue;
    this->fValueError = value.fValueError;
    this->fValueM2    = value.fValueM2;
  }
  return *this;
}

References fBase_Raw, fDiff_Raw, fPeak_Raw, fValue, fValue_Raw, fValueError, fValueM2, VQwDataElement::IsNameEmpty(), VQwHardwareChannel::operator=(), and QwADC18_Channel().

Here is the call graph for this function:

◆ PrintErrorCounterHead()

void QwADC18_Channel::PrintErrorCounterHead ( )

static

Definition at line 1307 of file QwADC18_Channel.cc.

{
  TString message;
  message  = Form("%30s","Device name");
  message += Form("%9s", "HW Sat");
  message += Form("%9s", "Sample");
  message += Form("%9s", "SW_HW");
  message += Form("%9s", "Sequence");
  message += Form("%9s", "SameHW");
  message += Form("%9s", "ZeroHW");
  message += Form("%9s", "EventCut");
  QwMessage << "---------------------------------------------------------------------------------------------" << QwLog::endl;
  QwMessage << message << QwLog::endl; 
  QwMessage << "---------------------------------------------------------------------------------------------" << QwLog::endl;
}

References QwLog::endl(), and QwMessage.

Here is the call graph for this function:

◆ PrintErrorCounters()

void QwADC18_Channel::PrintErrorCounters ( ) const

overridevirtual

report number of events failed due to HW and event cut failure

Reimplemented from VQwDataElement.

Definition at line 1328 of file QwADC18_Channel.cc.

{
  TString message;
  if (fErrorCount_sample || fErrorCount_SW_HW || fErrorCount_Sequence || fErrorCount_SameHW || fErrorCount_ZeroHW || fErrorCount_HWSat || fNumEvtsWithEventCutsRejected) {
    message  = Form("%30s", GetElementName().Data());
    message += Form("%9d", fErrorCount_HWSat);
    message += Form("%9d", fErrorCount_sample);
    message += Form("%9d", fErrorCount_SW_HW);
    message += Form("%9d", fErrorCount_Sequence);
    message += Form("%9d", fErrorCount_SameHW);
    message += Form("%9d", fErrorCount_ZeroHW);
    message += Form("%9d", fNumEvtsWithEventCutsRejected);
 
    if((fDataToSave == kRaw) && (!kFoundPedestal||!kFoundGain)){
      message += " >>>>> No Pedestal or Gain in map file";
    }
 
    QwMessage << message << QwLog::endl;
  }
}

References QwLog::endl(), VQwHardwareChannel::fDataToSave, fErrorCount_HWSat, fErrorCount_SameHW, fErrorCount_sample, fErrorCount_Sequence, fErrorCount_SW_HW, fErrorCount_ZeroHW, fNumEvtsWithEventCutsRejected, VQwDataElement::GetElementName(), VQwHardwareChannel::kFoundGain, VQwHardwareChannel::kFoundPedestal, VQwDataElement::kRaw, and QwMessage.

Here is the call graph for this function:

◆ PrintErrorCounterTail()

void QwADC18_Channel::PrintErrorCounterTail ( )

static

Definition at line 1323 of file QwADC18_Channel.cc.

{
  QwMessage << "---------------------------------------------------------------------------------------------" << QwLog::endl;
}

References QwLog::endl(), and QwMessage.

Here is the call graph for this function:

◆ PrintInfo()

void QwADC18_Channel::PrintInfo ( ) const

overridevirtual

Print multiple lines of information about this data element.

Reimplemented from VQwDataElement.

Definition at line 495 of file QwADC18_Channel.cc.

{
  QwMessage<<"***************************************"<<QwLog::endl;
  QwMessage<<"Subsystem "<<GetSubsystemName()<<QwLog::endl;
  QwMessage<<"Beam Instrument Type: "<<GetModuleType()<<QwLog::endl;
  QwMessage<<"QwADC18 channel: "<<GetElementName()<<QwLog::endl;
  QwMessage<<"fPedestal= "<< fPedestal<<QwLog::endl;
  QwMessage<<"fCalibrationFactor= "<<fCalibrationFactor<<QwLog::endl;
  QwMessage<<"fSequenceNumber= "<<fSequenceNumber<<QwLog::endl;
  QwMessage<<"fNumberOfSamples= "<<fNumberOfSamples<<QwLog::endl;
  QwMessage<<"fDiff_Raw= "<<fDiff_Raw<<QwLog::endl;
  QwMessage<<"fPeak_Raw= "<<fPeak_Raw<<QwLog::endl;
  QwMessage<<"fBase_Raw= "<<fBase_Raw<<QwLog::endl;
  QwMessage<<"fValue = "<<std::setprecision(8) <<fValue << QwLog::endl;
}

References QwLog::endl(), fBase_Raw, VQwHardwareChannel::fCalibrationFactor, fDiff_Raw, fNumberOfSamples, fPeak_Raw, VQwHardwareChannel::fPedestal, fSequenceNumber, fValue, VQwDataElement::GetElementName(), VQwDataElement::GetModuleType(), VQwDataElement::GetSubsystemName(), and QwMessage.

Here is the call graph for this function:

◆ PrintValue()

void QwADC18_Channel::PrintValue ( ) const

overridevirtual

Print single line of value and error of this data element.

Reimplemented from VQwDataElement.

Definition at line 1176 of file QwADC18_Channel.cc.

{
  QwMessage << std::setprecision(8)
            << std::setw(18) << std::left << GetSubsystemName()    << ""
            << std::setw(18) << std::left << GetModuleType()       << ""
            << std::setw(18) << std::left << GetElementName()      << ""
            << std::setw(12) << std::left << GetValue()            << " +/- "
            << std::setw(12) << std::left << GetValueError()       << " sig "
        << std::setw(12) << std::left << GetValueWidth()       << ""
            << std::setw(12) << std::left << fGoodEventCount       << ""
            << QwLog::endl;
}

References QwLog::endl(), VQwDataElement::fGoodEventCount, VQwDataElement::GetElementName(), VQwDataElement::GetModuleType(), VQwDataElement::GetSubsystemName(), GetValue(), GetValueError(), GetValueWidth(), and QwMessage.

Here is the call graph for this function:

◆ ProcessDataWord()

Int_t QwADC18_Channel::ProcessDataWord ( UInt_t word )

Definition at line 333 of file QwADC18_Channel.cc.

{
  // "Actual" values from data word
  UInt_t act_dtype  = (rawd & mask2422x) >> 22;
  [[maybe_unused]]
  UInt_t act_chan   = (act_dtype != 4) ?
                      ((rawd & mask3029x) >> 29) : 0;
  UInt_t act_dvalue = (rawd & mask2625x) >> 25;
  UInt_t act_snum   = (act_dtype == 1 || act_dtype == 2) ?
                      ((rawd & mask2118x) >> 18) : 0;
 
  // Interpret by data type
  UInt_t value_raw = 0;
  switch (act_dtype) {
    case 0: // Diff word
      static UInt_t prev_dvalue = act_dvalue;
      if (act_dvalue != prev_dvalue) {
        QwError << "QwADC18_Channel::ProcessEvBuffer: Number of samples changed " << act_dvalue << " " << prev_dvalue << QwLog::endl;
        return 0;
      }
      value_raw = rawd & mask200x;
      if (rawd & mask21x) value_raw = -((~value_raw & 0x1fffffff) + 1);
      fNumberOfSamples = (1 << act_dvalue);
      return value_raw;
      break;
    case 1: // Peak word
    case 2: // Base word
      if (act_snum != fNumberOfSamples) {
        QwError << "QwADC18_Channel::ProcessEvBuffer: Number of samples changed " << act_snum << " " << fNumberOfSamples << QwLog::endl;
        return 0;
      }
      if (act_dvalue != 0) {
        QwError << "QwADC18_Channel::ProcessEvBuffer: Divider value non-zero 0x" << std::hex << rawd << std::dec << QwLog::endl;
        return 0;
      }
      return rawd & mask170x;
      break;
    case 4: // DAC word
      if (act_dvalue != 0) {
        QwError << "QwADC18_Channel::ProcessEvBuffer: Divider value non-zero 0x" << std::hex << rawd << std::dec << QwLog::endl;
        return 0;
      }
      return rawd & mask150x;
      break;
    default:
      QwError << "QwADC18_Channel::ProcessEvBuffer: Unknown data type 0x" << std::hex << rawd << std::dec << QwLog::endl;
      return 0;
  }
}

References QwLog::endl(), fNumberOfSamples, mask150x, mask170x, mask200x, mask2118x, mask21x, mask2422x, mask2625x, mask3029x, and QwError.

Referenced by ProcessEvBuffer().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ProcessEvBuffer()

Int_t QwADC18_Channel::ProcessEvBuffer	(	UInt_t *	buffer,
		UInt_t	num_words_left,
		UInt_t	subelement = 0 )

overridevirtual

Process the CODA event buffer for this element.

Implements VQwDataElement.

Definition at line 384 of file QwADC18_Channel.cc.

{
  Bool_t debug = false;
 
  // Print buffer
  if (debug) {
    QwOut << GetElementName() << " : " << QwLog::endl << std::hex;
    UInt_t n = 25;
    for (size_t i = 0; i < num_words_left; i++) {
      QwOut << "0x" << std::setfill('0') << std::setw(8) << buffer[i] << " ";
      if (i % n == n - 1) QwOut << QwLog::endl;
    }
    QwOut << std::dec << std::setfill(' ') << std::setw(0) << QwLog::endl;
  }
 
  UInt_t words_read = 0;
  if (IsNameEmpty()) {
    //  This channel is not used, but is present in the data stream.
    //  Skip over this data.
    words_read = kDataWordsPerChannel;
  } else if (num_words_left >= fNumberOfDataWords) {
 
    // Is this a header word?
    if (IsHeaderWord(buffer[0])) {
 
      // Debug output
      if (debug) {
        QwOut << " : header " << std::hex;
        UInt_t n = kHeaderWordsPerModule;
        for (size_t i = 0; i < n && i < num_words_left; i++) {
          QwOut << "0x" << std::setfill('0') << std::setw(8) << buffer[i] << " ";
        }
        QwOut << std::dec << std::setfill(' ') << std::setw(0) << QwLog::endl;
      }
 
      // Check if enough words left
      if (num_words_left < kHeaderWordsPerModule) {
        QwError << "QwADC18_Channel::ProcessEvBuffer: Not enough words left!" << QwLog::endl;
        return num_words_left;
      }
 
      // FIXME Catch 0xfa180bad
 
      // Header word: read DAC value
      fValue_Raw = ProcessDataWord(buffer[kHeaderWordsPerModule-1]);
 
      words_read = kHeaderWordsPerModule;
 
    } else {
 
      // Debug output
      if (debug) {
        QwOut << " : channel " << std::hex;
        UInt_t n = kDataWordsPerChannel;
        for (size_t i = 0; i < n && i < num_words_left; i++) {
          QwOut << "0x" << std::setfill('0') << std::setw(8) << buffer[i] << " ";
        }
        QwOut << std::dec << std::setfill(' ') << std::setw(0) << QwLog::endl;
      }
 
      // Check if enough words left
      if (num_words_left < kDataWordsPerChannel) {
        QwError << "QwADC18_Channel::ProcessEvBuffer: Not enough words left!" << QwLog::endl;
        return num_words_left;
      }
 
      // Data channel words: read diff, peak, base
      fDiff_Raw = ProcessDataWord(buffer[0]);
      fPeak_Raw = ProcessDataWord(buffer[1]);
      fBase_Raw = ProcessDataWord(buffer[2]);
 
      words_read = kDataWordsPerChannel;
    }
 
  } else {
    QwError << "QwADC18_Channel::ProcessEvBuffer: Not enough words!" << QwLog::endl;
  }
 
  return words_read;
}

References QwLog::endl(), fBase_Raw, fDiff_Raw, VQwHardwareChannel::fNumberOfDataWords, fPeak_Raw, fValue_Raw, VQwDataElement::GetElementName(), IsHeaderWord(), VQwDataElement::IsNameEmpty(), kDataWordsPerChannel, kHeaderWordsPerModule, ProcessDataWord(), QwError, and QwOut.

Here is the call graph for this function:

◆ ProcessEvent()

void QwADC18_Channel::ProcessEvent ( )

overridevirtual

Process the event data according to pedestal and calibration factor.

Implements VQwHardwareChannel.

Definition at line 467 of file QwADC18_Channel.cc.

{
  if (fNumberOfSamples == 0 && fDiff_Raw == 0) {
    //  There isn't valid data for this channel.  Just flag it and move on.
    fValue = 0.0;
    fValueM2 = 0.0;
    fErrorFlag |= kErrorFlag_sample;
  } else if (fNumberOfSamples == 0) {
    //  This is probably a more serious problem.
    QwWarning << "QwADC18_Channel::ProcessEvent:  Channel "
              << GetElementName()
              << " has fNumberOfSamples == 0 but has valid data in the value filed."
              << "Flag this as an error."
              << QwLog::endl;
    fValue = 0.0;
    fValueM2 = 0.0;
    fErrorFlag |= kErrorFlag_sample;
  } else {
    fValue = fCalibrationFactor * ( (Double_t(fDiff_Raw) / fNumberOfSamples) - fPedestal );
    fValueM2 = 0.0; // second moment is zero for single events
  }
}

References QwLog::endl(), VQwHardwareChannel::fCalibrationFactor, fDiff_Raw, VQwDataElement::fErrorFlag, fNumberOfSamples, VQwHardwareChannel::fPedestal, fValue, fValueM2, VQwDataElement::GetElementName(), kErrorFlag_sample, and QwWarning.

Here is the call graph for this function:

◆ Product()

void QwADC18_Channel::Product	(	const QwADC18_Channel &	value1,
		const QwADC18_Channel &	value2 )

Definition at line 1022 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()){
    fValue = value1.fValue * value2.fValue;
    fDiff_Raw = 0;
    fPeak_Raw = 0;
    fBase_Raw = 0;
 
    // Remaining variables
    fGoodEventCount  = value2.fGoodEventCount;
    fErrorFlag = (value1.fErrorFlag | value2.fErrorFlag);//error code is ORed.
  }
}

References fBase_Raw, fDiff_Raw, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fPeak_Raw, fValue, VQwDataElement::IsNameEmpty(), and QwADC18_Channel().

Here is the call graph for this function:

◆ RandomizeEventData()

void QwADC18_Channel::RandomizeEventData	(	int	helicity = 0.0,
		double	time = 0.0 )

overridevirtual

Internally generate random event data.

Implements MQwMockable.

Definition at line 248 of file QwADC18_Channel.cc.

{
  // Calculate drift (if time is not specified, it stays constant at zero)
  Double_t drift = 0.0;
  for (UInt_t i = 0; i < fMockDriftFrequency.size(); i++) {
    drift += fMockDriftAmplitude[i] * sin(2.0 * Qw::pi * fMockDriftFrequency[i] * time + fMockDriftPhase[i]);
  }
 
  Double_t value = fMockGaussianMean * (1 + helicity * fMockAsymmetry)
    + fMockGaussianSigma * GetRandomValue()
    + drift;
 
  fValue     = value;
  fSequenceNumber = 0;
  fNumberOfSamples = fNumberOfSamples_map;
}

References MQwMockable::fMockAsymmetry, MQwMockable::fMockDriftAmplitude, MQwMockable::fMockDriftFrequency, MQwMockable::fMockDriftPhase, MQwMockable::fMockGaussianMean, MQwMockable::fMockGaussianSigma, fNumberOfSamples, fNumberOfSamples_map, fSequenceNumber, fValue, MQwMockable::GetRandomValue(), and Qw::pi.

Here is the call graph for this function:

◆ Ratio()

void QwADC18_Channel::Ratio	(	const QwADC18_Channel &	numer,
		const QwADC18_Channel &	denom )

Definition at line 954 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty()) {
    *this  = numer;
    *this /= denom;
 
    //  Set the raw values to zero.
    fDiff_Raw = 0;
    fPeak_Raw = 0;
    fBase_Raw = 0;
 
    // Remaining variables
    fGoodEventCount  = denom.fGoodEventCount;
    fErrorFlag = (numer.fErrorFlag|denom.fErrorFlag);//error code is ORed.
  }
}

References fBase_Raw, fDiff_Raw, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fPeak_Raw, VQwDataElement::IsNameEmpty(), and QwADC18_Channel().

Here is the call graph for this function:

◆ Scale()

void QwADC18_Channel::Scale ( Double_t Offset )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 1047 of file QwADC18_Channel.cc.

{
  if (!IsNameEmpty())
    {
      fValue   *= scale;
      fValueM2 *= scale * scale;
    }
}

References fValue, fValueM2, and VQwDataElement::IsNameEmpty().

Here is the call graph for this function:

◆ ScaledAdd()

void QwADC18_Channel::ScaledAdd	(	Double_t	scale,
		const VQwHardwareChannel *	value )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 1349 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* input = dynamic_cast<const QwADC18_Channel*>(value);
 
  // follows same steps as += but w/ scaling factor
  if(input!=NULL && !IsNameEmpty()){
    //     QwWarning << "Adding " << input->GetElementName()
    //        << " to " << GetElementName()
    //        << " with scale factor " << scale
    //        << QwLog::endl;
    //     PrintValue();
    //     input->PrintValue();
    this->fDiff_Raw = 0;
    this->fPeak_Raw = 0;
    this->fBase_Raw = 0;
    this->fValue += scale * input->fValue;
    this->fValueM2 = 0.0;
    this->fNumberOfSamples += input->fNumberOfSamples;
    this->fSequenceNumber  =  0;
    this->fErrorFlag       |= (input->fErrorFlag);   
  } else if (input == NULL && value != NULL) {
    TString loc="Standard exception from QwADC18_Channel::ScaledAdd "
        +value->GetElementName()+" "
        +this->GetElementName()+" are not of the same type";
    throw(std::invalid_argument(loc.Data()));
  }
}

References fBase_Raw, fDiff_Raw, VQwDataElement::fErrorFlag, fNumberOfSamples, fPeak_Raw, fSequenceNumber, fValue, fValueM2, VQwDataElement::GetElementName(), VQwDataElement::IsNameEmpty(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ SetADC18SaturationLimt()

void QwADC18_Channel::SetADC18SaturationLimt ( Double_t sat_volts = 8.5 )

inline

Definition at line 205 of file QwADC18_Channel.h.

                                                       {
    fSaturationABSLimit = sat_volts;
  }

References fSaturationABSLimit.

Referenced by QwADC18_Channel(), and QwADC18_Channel().

Here is the caller graph for this function:

◆ SetCalibrationToVolts()

void QwADC18_Channel::SetCalibrationToVolts ( )

inline

Definition at line 244 of file QwADC18_Channel.h.

244{SetCalibrationFactor(kADC18_VoltsPerBit);};

VQwHardwareChannel::SetCalibrationFactor

void SetCalibrationFactor(Double_t factor)

Definition VQwHardwareChannel.h:241

References kADC18_VoltsPerBit, and VQwHardwareChannel::SetCalibrationFactor().

Here is the call graph for this function:

◆ SetDefaultSampleSize()

void QwADC18_Channel::SetDefaultSampleSize ( size_t num_samples_map )

inline

Definition at line 100 of file QwADC18_Channel.h.

                                                    {
    // This will be checked against the no.of samples read by the module
    fNumberOfSamples_map = num_samples_map;
  };

References fNumberOfSamples_map.

Referenced by LoadChannelParameters().

Here is the caller graph for this function:

◆ SetEventData()

void QwADC18_Channel::SetEventData ( Double_t value )

Definition at line 277 of file QwADC18_Channel.cc.

{
  fValue     = value;
  fDiff_Raw = (UInt_t) value;
  fPeak_Raw = (UInt_t) value;
  fBase_Raw = (UInt_t) 0;
  fSequenceNumber = 0;
  fNumberOfSamples = fNumberOfSamples_map;
}

References fBase_Raw, fDiff_Raw, fNumberOfSamples, fNumberOfSamples_map, fPeak_Raw, fSequenceNumber, and fValue.

◆ SetRawEventData()

void QwADC18_Channel::SetRawEventData ( )

overridevirtual

Implements MQwMockable.

Definition at line 287 of file QwADC18_Channel.cc.

{
  fNumberOfSamples = fNumberOfSamples_map;
  fDiff_Raw = Int_t(fValue / fCalibrationFactor + fPedestal) * fNumberOfSamples;
  fPeak_Raw = Int_t(fValue / fCalibrationFactor + fPedestal) * fNumberOfSamples;
};

References VQwHardwareChannel::fCalibrationFactor, fDiff_Raw, fNumberOfSamples, fNumberOfSamples_map, fPeak_Raw, VQwHardwareChannel::fPedestal, and fValue.

Referenced by SmearByResolution().

Here is the caller graph for this function:

◆ SmearByResolution()

void QwADC18_Channel::SmearByResolution ( double resolution )

overridevirtual

Implements MQwMockable.

Definition at line 265 of file QwADC18_Channel.cc.

{
  fValue = resolution*GetRandomValue();
  fValueM2 = 0.0;
  fNumberOfSamples = fNumberOfSamples_map;
  SetRawEventData();
}

References fNumberOfSamples, fNumberOfSamples_map, fValue, fValueM2, MQwMockable::GetRandomValue(), and SetRawEventData().

Here is the call graph for this function:

◆ SubtractValueFrom()

void QwADC18_Channel::SubtractValueFrom ( const VQwHardwareChannel * valueptr )

overridevirtual

Implements VQwHardwareChannel.

Definition at line 789 of file QwADC18_Channel.cc.

{
  const QwADC18_Channel* tmpptr;
  tmpptr = dynamic_cast<const QwADC18_Channel*>(valueptr);
  if (tmpptr!=NULL){
    *this -= *tmpptr;
  } else {
    TString loc="Standard exception from QwADC18_Channel::SubtractValueFrom = "
      +valueptr->GetElementName()+" is an incompatible type.";
    throw std::invalid_argument(loc.Data());
  }
}

References VQwDataElement::GetElementName(), QwADC18_Channel(), and VQwHardwareChannel::VQwHardwareChannel().

Here is the call graph for this function:

◆ Sum()

void QwADC18_Channel::Sum	(	const QwADC18_Channel &	value1,
		const QwADC18_Channel &	value2 )

Definition at line 942 of file QwADC18_Channel.cc.

{
  *this  = value1;
  *this += value2;
}

References QwADC18_Channel().

Here is the call graph for this function:

Friends And Related Symbol Documentation

◆ operator<<

std::ostream & operator<<	(	std::ostream &	stream,
		const QwADC18_Channel &	channel )

friend

Definition at line 1189 of file QwADC18_Channel.cc.

{
  stream << channel.GetValue();
  return stream;
}

References GetValue(), and QwADC18_Channel().

Field Documentation

◆ bBlock

Bool_t QwADC18_Channel::bBlock

private

Definition at line 346 of file QwADC18_Channel.h.

◆ bBlock_raw

Bool_t QwADC18_Channel::bBlock_raw

private

Definition at line 347 of file QwADC18_Channel.h.

◆ bDEBUG

const Bool_t QwADC18_Channel::bDEBUG =kFALSE

staticprivate

debugging display purposes

For ADC18 data element trimming uses

Definition at line 341 of file QwADC18_Channel.h.

Referenced by ApplyHWChecks(), and MatchNumberOfSamples().

◆ bDevice_Error_Code

Bool_t QwADC18_Channel::bDevice_Error_Code

private

Definition at line 349 of file QwADC18_Channel.h.

◆ bHw_sum

Bool_t QwADC18_Channel::bHw_sum

private

Definition at line 344 of file QwADC18_Channel.h.

◆ bHw_sum_raw

Bool_t QwADC18_Channel::bHw_sum_raw

private

Definition at line 345 of file QwADC18_Channel.h.

◆ bNum_samples

Bool_t QwADC18_Channel::bNum_samples

private

Definition at line 348 of file QwADC18_Channel.h.

◆ bSequence_number

Bool_t QwADC18_Channel::bSequence_number

private

Definition at line 350 of file QwADC18_Channel.h.

◆ fADC_Same_NumEvt

Int_t QwADC18_Channel::fADC_Same_NumEvt

private

Keep track of how many events with same ADC value returned.

Definition at line 331 of file QwADC18_Channel.h.

Referenced by InitializeChannel().

◆ fBase_Raw

UInt_t QwADC18_Channel::fBase_Raw

private

Definition at line 268 of file QwADC18_Channel.h.

Referenced by ClearEventData(), FillTreeVector(), operator*=(), operator=(), PrintInfo(), ProcessEvBuffer(), Product(), Ratio(), ScaledAdd(), and SetEventData().

◆ fDAC

QwADC18_Channel* QwADC18_Channel::fDAC

private

Pointer to the DAC channel for this channel.

Definition at line 301 of file QwADC18_Channel.h.

◆ fDiff_Raw

UInt_t QwADC18_Channel::fDiff_Raw

private

Definition at line 267 of file QwADC18_Channel.h.

Referenced by ClearEventData(), FillTreeVector(), operator*=(), operator=(), PrintInfo(), ProcessEvBuffer(), ProcessEvent(), Product(), Ratio(), ScaledAdd(), SetEventData(), and SetRawEventData().

◆ fErrorCount_HWSat

Int_t QwADC18_Channel::fErrorCount_HWSat

private

check to see ADC channel is saturated

Definition at line 320 of file QwADC18_Channel.h.

Referenced by InitializeChannel(), and PrintErrorCounters().

◆ fErrorCount_SameHW

Int_t QwADC18_Channel::fErrorCount_SameHW

private

check to see ADC returning same HW value

Definition at line 324 of file QwADC18_Channel.h.

Referenced by IncrementErrorCounters(), InitializeChannel(), and PrintErrorCounters().

◆ fErrorCount_sample

Int_t QwADC18_Channel::fErrorCount_sample

private

for sample size check

Definition at line 321 of file QwADC18_Channel.h.

Referenced by IncrementErrorCounters(), InitializeChannel(), and PrintErrorCounters().

◆ fErrorCount_Sequence

Int_t QwADC18_Channel::fErrorCount_Sequence

private

sequence number check

Definition at line 323 of file QwADC18_Channel.h.

Referenced by IncrementErrorCounters(), InitializeChannel(), and PrintErrorCounters().

◆ fErrorCount_SW_HW

Int_t QwADC18_Channel::fErrorCount_SW_HW

private

HW_sum==SW_sum check.

Definition at line 322 of file QwADC18_Channel.h.

Referenced by IncrementErrorCounters(), InitializeChannel(), and PrintErrorCounters().

◆ fErrorCount_ZeroHW

Int_t QwADC18_Channel::fErrorCount_ZeroHW

private

check to see ADC returning zero

Definition at line 325 of file QwADC18_Channel.h.

Referenced by IncrementErrorCounters(), InitializeChannel(), and PrintErrorCounters().

◆ fNumberOfSamples

UInt_t QwADC18_Channel::fNumberOfSamples

private

Number of samples read through the module.

Definition at line 316 of file QwADC18_Channel.h.

Referenced by ClearEventData(), ForceMapfileSampleSize(), GetAverageVolts(), GetNumberOfSamples(), InitializeChannel(), MatchNumberOfSamples(), PrintInfo(), ProcessDataWord(), ProcessEvent(), RandomizeEventData(), ScaledAdd(), SetEventData(), SetRawEventData(), and SmearByResolution().

◆ fNumberOfSamples_map

UInt_t QwADC18_Channel::fNumberOfSamples_map

private

Number of samples in the expected to read through the module. This value is set in the QwBeamline map file.

Definition at line 317 of file QwADC18_Channel.h.

Referenced by ApplyHWChecks(), ForceMapfileSampleSize(), InitializeChannel(), QwADC18_Channel(), QwADC18_Channel(), RandomizeEventData(), SetDefaultSampleSize(), SetEventData(), SetRawEventData(), and SmearByResolution().

◆ fNumEvtsWithEventCutsRejected

Int_t QwADC18_Channel::fNumEvtsWithEventCutsRejected

private

Counts the Event cut rejected events.

Definition at line 327 of file QwADC18_Channel.h.

Referenced by IncrementErrorCounters(), InitializeChannel(), and PrintErrorCounters().

◆ fPeak_Raw

UInt_t QwADC18_Channel::fPeak_Raw

private

Definition at line 269 of file QwADC18_Channel.h.

Referenced by ClearEventData(), EncodeEventData(), FillTreeVector(), operator*=(), operator=(), PrintInfo(), ProcessEvBuffer(), Product(), Ratio(), ScaledAdd(), SetEventData(), and SetRawEventData().

◆ fPrev_HardwareBlockSum

Double_t QwADC18_Channel::fPrev_HardwareBlockSum

private

Previous Module-based sum of the four sub-blocks.

Definition at line 334 of file QwADC18_Channel.h.

◆ fPreviousSequenceNumber

UInt_t QwADC18_Channel::fPreviousSequenceNumber

private

Previous event sequence number for this channel.

Definition at line 315 of file QwADC18_Channel.h.

Referenced by InitializeChannel().

◆ fRunningSum

QwADC18_Channel* QwADC18_Channel::fRunningSum

private

Pointer to the running sum for this channel.

Definition at line 298 of file QwADC18_Channel.h.

Referenced by InitializeChannel().

◆ fSaturationABSLimit

Double_t QwADC18_Channel::fSaturationABSLimit

private

absolute value of the ADC18 saturation volt

Definition at line 338 of file QwADC18_Channel.h.

Referenced by GetADC18SaturationLimt(), QwADC18_Channel(), QwADC18_Channel(), and SetADC18SaturationLimt().

◆ fSequenceNo_Counter

Int_t QwADC18_Channel::fSequenceNo_Counter

private

Internal counter to keep track of the sequence number.

Definition at line 333 of file QwADC18_Channel.h.

Referenced by ApplyHWChecks(), and InitializeChannel().

◆ fSequenceNo_Prev

Int_t QwADC18_Channel::fSequenceNo_Prev

private

Keep the sequence number of the last event.

Definition at line 332 of file QwADC18_Channel.h.

Referenced by ApplyHWChecks(), and InitializeChannel().

◆ fSequenceNumber

UInt_t QwADC18_Channel::fSequenceNumber

private

Event sequence number for this channel.

Definition at line 314 of file QwADC18_Channel.h.

Referenced by ClearEventData(), GetSequenceNumber(), MatchSequenceNumber(), PrintInfo(), RandomizeEventData(), ScaledAdd(), and SetEventData().

◆ fValue

Double_t QwADC18_Channel::fValue

private

Definition at line 272 of file QwADC18_Channel.h.

Referenced by AccumulateRunningSum(), AddChannelOffset(), ArcTan(), AssignScaledValue(), Blind(), Blind(), ClearEventData(), ConstructBranch(), FillTreeVector(), GetAverageVolts(), GetValue(), operator*=(), operator+=(), operator-=(), operator/=(), operator=(), PrintInfo(), ProcessEvent(), Product(), RandomizeEventData(), Scale(), ScaledAdd(), SetEventData(), SetRawEventData(), and SmearByResolution().

◆ fValue_Raw

UInt_t QwADC18_Channel::fValue_Raw

private

Definition at line 270 of file QwADC18_Channel.h.

Referenced by FillTreeVector(), GetRawValue(), operator=(), and ProcessEvBuffer().

◆ fValueError

Double_t QwADC18_Channel::fValueError

private

Definition at line 274 of file QwADC18_Channel.h.

Referenced by AssignScaledValue(), CalculateRunningAverage(), ClearEventData(), FillTreeVector(), GetValueError(), and operator=().

◆ fValueM2

Double_t QwADC18_Channel::fValueM2

private

Definition at line 273 of file QwADC18_Channel.h.

Referenced by AccumulateRunningSum(), AssignScaledValue(), CalculateRunningAverage(), ClearEventData(), FillTreeVector(), GetValueM2(), operator*=(), operator+=(), operator-=(), operator/=(), operator=(), ProcessEvent(), Scale(), ScaledAdd(), and SmearByResolution().

◆ kADC18_VoltsPerBit

const Double_t QwADC18_Channel::kADC18_VoltsPerBit = (20./(1<<18))

staticprivate

Conversion factor to translate the average bit count in an ADC channel into average voltage. The base factor is roughly 76 uV per count, and zero counts corresponds to zero voltage. Store as the exact value for 20 V range, 18 bit ADC.

Definition at line 305 of file QwADC18_Channel.h.

Referenced by GetAverageVolts(), and SetCalibrationToVolts().

◆ kDataWordsPerChannel

const Int_t QwADC18_Channel::kDataWordsPerChannel = 3

staticprivate

Definition at line 282 of file QwADC18_Channel.h.

Referenced by EncodeEventData(), GetBufferOffset(), and ProcessEvBuffer().

◆ kDEBUG

const Bool_t QwADC18_Channel::kDEBUG = kFALSE

staticprivate

Definition at line 277 of file QwADC18_Channel.h.

◆ kFooterWordsPerBank

const Int_t QwADC18_Channel::kFooterWordsPerBank = 1

staticprivate

Definition at line 279 of file QwADC18_Channel.h.

◆ kFooterWordsPerModule

const Int_t QwADC18_Channel::kFooterWordsPerModule = 1

staticprivate

Definition at line 281 of file QwADC18_Channel.h.

Referenced by GetBufferOffset().

◆ kHeaderWordsPerBank

const Int_t QwADC18_Channel::kHeaderWordsPerBank = 1

staticprivate

Definition at line 278 of file QwADC18_Channel.h.

Referenced by GetBufferOffset().

◆ kHeaderWordsPerModule

const Int_t QwADC18_Channel::kHeaderWordsPerModule = 12

staticprivate

Definition at line 280 of file QwADC18_Channel.h.

Referenced by GetBufferOffset(), and ProcessEvBuffer().

◆ kMaxChannels

const Int_t QwADC18_Channel::kMaxChannels = 4

staticprivate

Definition at line 283 of file QwADC18_Channel.h.

Referenced by GetBufferOffset().

◆ kTimePerSample

const Double_t QwADC18_Channel::kTimePerSample = 2.0 * Qw::us

static

Definition at line 47 of file QwADC18_Channel.h.

◆ mask150x

const UInt_t QwADC18_Channel::mask150x = 0x0000ffff

staticprivate

Definition at line 294 of file QwADC18_Channel.h.

Referenced by ProcessDataWord().

◆ mask170x

const UInt_t QwADC18_Channel::mask170x = 0x0003ffff

staticprivate

Definition at line 293 of file QwADC18_Channel.h.

Referenced by EncodeEventData(), and ProcessDataWord().

◆ mask200x

const UInt_t QwADC18_Channel::mask200x = 0x001fffff

staticprivate

Definition at line 291 of file QwADC18_Channel.h.

Referenced by ProcessDataWord().

◆ mask2118x

const UInt_t QwADC18_Channel::mask2118x = 0x003c0000

staticprivate

Definition at line 292 of file QwADC18_Channel.h.

Referenced by ProcessDataWord().

◆ mask21x

const UInt_t QwADC18_Channel::mask21x = 0x00200000

staticprivate

Definition at line 290 of file QwADC18_Channel.h.

Referenced by ProcessDataWord().

◆ mask2422x

const UInt_t QwADC18_Channel::mask2422x = 0x01c00000

staticprivate

Definition at line 289 of file QwADC18_Channel.h.

Referenced by EncodeEventData(), and ProcessDataWord().

◆ mask2625x

const UInt_t QwADC18_Channel::mask2625x = 0x06000000

staticprivate

Definition at line 288 of file QwADC18_Channel.h.

Referenced by EncodeEventData(), and ProcessDataWord().

◆ mask3029x

const UInt_t QwADC18_Channel::mask3029x = 0x60000000

staticprivate

Definition at line 287 of file QwADC18_Channel.h.

Referenced by EncodeEventData(), and ProcessDataWord().

◆ mask31x

const UInt_t QwADC18_Channel::mask31x = 0x80000000

staticprivate

Definition at line 286 of file QwADC18_Channel.h.

Referenced by EncodeEventData(), and IsHeaderWord().

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Private Attributes

Static Private Attributes

Friends

Channel configuration data members

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ QwADC18_Channel() [1/4]

◆ QwADC18_Channel() [2/4]

◆ QwADC18_Channel() [3/4]

◆ QwADC18_Channel() [4/4]

◆ ~QwADC18_Channel()

Member Function Documentation

◆ AccumulateRunningSum() [1/2]

◆ AccumulateRunningSum() [2/2]

◆ AddChannelOffset()

◆ AddValueFrom()

◆ ApplyHWChecks()

◆ ApplySingleEventCuts() [1/2]

◆ ApplySingleEventCuts() [2/2]

◆ ArcTan()

◆ AssignScaledValue()

◆ AssignValueFrom()

◆ Blind() [1/2]

◆ Blind() [2/2]

◆ CalculateRunningAverage()

◆ ClearEventData()

◆ Clone() [1/2]

◆ Clone() [2/2]

◆ ConstructBranch()

◆ ConstructBranchAndVector()

◆ ConstructHistograms()

◆ DeaccumulateRunningSum() [1/2]

◆ DeaccumulateRunningSum() [2/2]

◆ Difference()

◆ DivideBy() [1/2]

◆ DivideBy() [2/2]

◆ EncodeEventData()

◆ FillHistograms()

◆ FillTreeVector()

◆ ForceMapfileSampleSize()

◆ GetADC18SaturationLimt()

◆ GetAverageVolts()

◆ GetBufferOffset()

◆ GetNumberOfSamples()

◆ GetRawValue() [1/2]

◆ GetRawValue() [2/2]

◆ GetSequenceNumber()

◆ GetValue() [1/2]

◆ GetValue() [2/2]

◆ GetValueError() [1/2]

◆ GetValueError() [2/2]

◆ GetValueM2() [1/2]

◆ GetValueM2() [2/2]

◆ GetValueWidth() [1/2]

◆ GetValueWidth() [2/2]

◆ IncrementErrorCounters()

◆ InitializeChannel() [1/2]

◆ InitializeChannel() [2/2]

◆ IsHeaderWord()

◆ LoadChannelParameters()

◆ LoadMockDataParameters()

◆ MatchNumberOfSamples()

◆ MatchSequenceNumber()

◆ MultiplyBy()

◆ operator*()

◆ operator*=() [1/2]

◆ operator*=() [2/2]

◆ operator+()

◆ operator+=() [1/2]

◆ operator+=() [2/2]

◆ operator-()

◆ operator-=() [1/2]

◆ operator-=() [2/2]

◆ operator/=() [1/2]

◆ operator/=() [2/2]