#include <QwVQWK_Channel.h>

Inheritance diagram for QwVQWK_Channel:

Collaboration diagram for QwVQWK_Channel:

Public Member Functions
	QwVQWK_Channel ()

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

	QwVQWK_Channel (const QwVQWK_Channel &value)

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

virtual	~QwVQWK_Channel ()

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

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

void	LoadChannelParameters (QwParameterFile &paramfile)

void	SetDefaultSampleSize (size_t num_samples_map)

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

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

void	SetHardwareSum (Double_t hwsum, UInt_t sequencenumber=0)

void	SetEventData (Double_t *block, UInt_t sequencenumber=0)

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

Int_t	ProcessEvBuffer (UInt_t *buffer, UInt_t num_words_left, UInt_t index=0)
	Decode the event data from a CODA buffer. More...

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

QwVQWK_Channel &	operator= (const QwVQWK_Channel &value)

void	AssignScaledValue (const QwVQWK_Channel &value, Double_t scale)

void	AssignValueFrom (const VQwDataElement *valueptr)

void	AddValueFrom (const VQwHardwareChannel *valueptr)

void	SubtractValueFrom (const VQwHardwareChannel *valueptr)

void	MultiplyBy (const VQwHardwareChannel *valueptr)

void	DivideBy (const VQwHardwareChannel *valueptr)

QwVQWK_Channel &	operator+= (const QwVQWK_Channel &value)

QwVQWK_Channel &	operator-= (const QwVQWK_Channel &value)

QwVQWK_Channel &	operator*= (const QwVQWK_Channel &value)

VQwHardwareChannel &	operator+= (const VQwHardwareChannel *input)

VQwHardwareChannel &	operator-= (const VQwHardwareChannel *input)

VQwHardwareChannel &	operator= (const VQwHardwareChannel input)

VQwHardwareChannel &	operator/= (const VQwHardwareChannel *input)

const QwVQWK_Channel	operator+ (const QwVQWK_Channel &value) const

const QwVQWK_Channel	operator- (const QwVQWK_Channel &value) const

const QwVQWK_Channel	operator* (const QwVQWK_Channel &value) const

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

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

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

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

void	DivideBy (const QwVQWK_Channel &denom)

void	AddChannelOffset (Double_t Offset)

void	Scale (Double_t Offset)

void	AccumulateRunningSum (const QwVQWK_Channel &value)

void	AccumulateRunningSum (const QwVQWK_Channel &value, Int_t count)

void	AccumulateRunningSum (const VQwHardwareChannel *value, Int_t count)

void	DeaccumulateRunningSum (const QwVQWK_Channel &value)

void	CalculateRunningAverage ()

Bool_t	MatchSequenceNumber (size_t seqnum)

Bool_t	MatchNumberOfSamples (size_t numsamp)

Bool_t	ApplySingleEventCuts (Double_t LL, Double_t UL)

Bool_t	ApplySingleEventCuts ()

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

void	SetVQWKSaturationLimt (Double_t sat_volts=8.5)

Double_t	GetVQWKSaturationLimt ()

Int_t	ApplyHWChecks ()

void	IncrementErrorCounters ()

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

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

void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)

void	ConstructBranch (TTree *tree, TString &prefix)

void	FillTreeVector (std::vector< Double_t > &values) const

Int_t	GetRawValue (size_t element) const

Double_t	GetValue (size_t element) const

Double_t	GetValueM2 (size_t element) const

Double_t	GetValueError (size_t element) const

Double_t	GetAverageVolts () const

size_t	GetSequenceNumber () const

size_t	GetNumberOfSamples () const

void	SetCalibrationToVolts ()

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

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

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

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

void	ScaledAdd (Double_t scale, const VQwHardwareChannel *value)

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

Public Member Functions inherited from VQwHardwareChannel
	VQwHardwareChannel ()

	VQwHardwareChannel (const VQwHardwareChannel &value)

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

virtual	~VQwHardwareChannel ()

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

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

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	InitializeChannel (TString name)
	Initialize the fields in this object. More...

void	SetEventCutMode (Int_t bcuts)

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

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

Double_t	GetEventCutUpperLimit () const

Double_t	GetEventCutLowerLimit () const

Double_t	GetStabilityLimit () const

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

void	UpdateErrorFlag (const VQwHardwareChannel &elem)

virtual UInt_t	GetErrorCode () const

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

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	AccumulateRunningSum (const VQwHardwareChannel *value)

virtual void	DeaccumulateRunningSum (const VQwHardwareChannel *value)

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

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

	VQwDataElement (const VQwDataElement &value)
	Copy constructor. More...

virtual	~VQwDataElement ()
	Virtual destructor. More...

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

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

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

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

UInt_t	GetGoodEventCount () const

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

virtual VQwDataElement &	operator-= (const VQwDataElement &value)
	Subtraction-assignment operator. More...

virtual void	Sum (const VQwDataElement &value1, const VQwDataElement &value2)
	Sum operator. More...

virtual void	Difference (const VQwDataElement &value1, const VQwDataElement &value2)
	Difference operator. More...

virtual void	Ratio (const VQwDataElement &numer, const VQwDataElement &denom)
	Ratio operator. More...

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

virtual void	SetNeedsExternalClock (Bool_t needed)

virtual Bool_t	NeedsExternalClock ()

virtual std::string	GetExternalClockName ()

virtual void	SetExternalClockPtr (const VQwHardwareChannel *clock)

virtual void	SetExternalClockName (const std::string name)

virtual Double_t	GetNormClockValue ()

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

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

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

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

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

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. More...

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

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

void	SetRandomEventAsymmetry (Double_t asymmetry)
	Set the helicity asymmetry. More...

Double_t	GetRandomValue ()

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

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

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

static void	PrintErrorCounterHead ()

static void	PrintErrorCounterTail ()

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

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

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

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

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

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(). More...

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

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

virtual void	UpdateErrorFlag (const UInt_t &error)

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

	MQwHistograms (const MQwHistograms &source)
	Copy constructor. More...

virtual	~MQwHistograms ()
	Virtual destructor. More...

virtual MQwHistograms &	operator= (const MQwHistograms &value)

void	Fill_Pointer (TH1_ptr hist_ptr, Double_t value)

void	AddHistogram (TH1 *h)
	Register a histogram. More...

Private Member Functions
Double_t	GetBlockValue (size_t blocknum) const

Double_t	GetBlockErrorValue (size_t blocknum) const

Double_t	GetHardwareSum () const

Double_t	GetHardwareSumM2 () const

Double_t	GetHardwareSumWidth () const

Double_t	GetHardwareSumError () const

Int_t	GetRawBlockValue (size_t blocknum) const

Int_t	GetRawHardwareSum () const

Int_t	GetRawSoftwareSum () const

Private Attributes
QwVQWK_Channel *	fRunningSum
	Pointer to the running sum for this channel. More...

size_t	fSequenceNumber
	Event sequence number for this channel. More...

size_t	fPreviousSequenceNumber
	Previous event sequence number for this channel. More...

size_t	fNumberOfSamples
	Number of samples read through the module. More...

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

Int_t	fErrorCount_HWSat
	check to see ADC channel is saturated More...

Int_t	fErrorCount_sample
	for sample size check More...

Int_t	fErrorCount_SW_HW
	HW_sum==SW_sum check. More...

Int_t	fErrorCount_Sequence
	sequence number check More...

Int_t	fErrorCount_SameHW
	check to see ADC returning same HW value More...

Int_t	fErrorCount_ZeroHW
	check to see ADC returning zero More...

Int_t	fNumEvtsWithEventCutsRejected
	Counts the Event cut rejected events. More...

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

Int_t	fSequenceNo_Prev
	Keep the sequence number of the last event. More...

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

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

Double_t	fSaturationABSLimit
	absolute value of the VQWK saturation volt More...

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

Channel configuration data members
Short_t	fBlocksPerEvent

Event data members---Raw values
Int_t	fBlock_raw [4]
	Array of the sub-block data as read from the module. More...

Int_t	fHardwareBlockSum_raw
	Module-based sum of the four sub-blocks as read from the module. More...

Int_t	fSoftwareBlockSum_raw
	Sum of the data in the four sub-blocks raw. More...

Event data members---Potentially calibrated values
Double_t	fBlock [4]
	Array of the sub-block data. More...

Double_t	fHardwareBlockSum
	Module-based sum of the four sub-blocks. More...

Calculation of the statistical moments
Double_t	fBlockM2 [4]
	Second moment of the sub-block. More...

Double_t	fBlockError [4]
	Uncertainty on the sub-block. More...

Double_t	fHardwareBlockSumM2
	Second moment of the hardware sum. More...

Double_t	fHardwareBlockSumError
	Uncertainty on the hardware sum. More...

Static Private Attributes
static const Bool_t	kDEBUG = kFALSE

static const Int_t	kWordsPerChannel = 6

static const Int_t	kMaxChannels = 8

static const Bool_t	bDEBUG =kFALSE
	debugging display purposes More...

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

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

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

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

UInt_t	fNumberOfSubElements
	Number of subelements in this data element. More...

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

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

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

Int_t	fGoodEventCount
	Number of good events accumulated in this element. More...

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) More...

UInt_t	fErrorConfigFlag
	contains the global/local/stability flags More...

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

Protected Attributes inherited from MQwMockable
bool	fUseExternalRandomVariable
	Flag to use an externally provided normal random variable. More...

double	fExternalRandomVariable
	Externally provided normal random variable. More...

Double_t	fMockAsymmetry
	Helicity asymmetry. More...

Double_t	fMockGaussianMean
	Mean of normal distribution. More...

Double_t	fMockGaussianSigma
	Sigma of normal distribution. More...

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

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

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

Static Protected Attributes inherited from MQwMockable
static boost::mt19937	fRandomnessGenerator
	Internal randomness generator. More...

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

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

Detailed Description

Definition at line 30 of file QwVQWK_Channel.h.

Constructor & Destructor Documentation

QwVQWK_Channel::QwVQWK_Channel ( )

inline

Definition at line 54 of file QwVQWK_Channel.h.

References InitializeChannel(), and SetVQWKSaturationLimt().

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

Here is the call graph for this function:

QwVQWK_Channel::QwVQWK_Channel	(	TString	name,
		TString	datatosave = `"raw"`
	)

inline

Definition at line 58 of file QwVQWK_Channel.h.

References InitializeChannel(), and SetVQWKSaturationLimt().

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

Here is the call graph for this function:

QwVQWK_Channel::QwVQWK_Channel ( const QwVQWK_Channel & value )

inline

Definition at line 62 of file QwVQWK_Channel.h.

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

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

inline

Definition at line 70 of file QwVQWK_Channel.h.

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

virtual QwVQWK_Channel::~QwVQWK_Channel ( )

inlinevirtual

Definition at line 78 of file QwVQWK_Channel.h.

78 { };

Member Function Documentation

void QwVQWK_Channel::AccumulateRunningSum ( const QwVQWK_Channel & value )

inline

Definition at line 143 of file QwVQWK_Channel.h.

References VQwDataElement::fGoodEventCount.

Referenced by QwEnergyCalculator::AccumulateRunningSum(), QwQPD::AccumulateRunningSum(), QwBPMCavity::AccumulateRunningSum(), QwLinearDiodeArray::AccumulateRunningSum(), QwIntegrationPMT::AccumulateRunningSum(), AccumulateRunningSum(), and DeaccumulateRunningSum().

                                                                {
     AccumulateRunningSum(value, value.fGoodEventCount);
   }

Here is the caller graph for this function:

void QwVQWK_Channel::AccumulateRunningSum	(	const QwVQWK_Channel &	value,
		Int_t	count
	)

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). Accumulate the running moments M1 and M2

Parameters

value	Object (single event or accumulated) to add to running moments
count	Number of good events in value

Definition at line 1232 of file QwVQWK_Channel.cc.

References QwLog::endl(), fBlock, fBlockM2, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fHardwareBlockSum, fHardwareBlockSumM2, fRunningSum, VQwDataElement::GetElementName(), PrintValue(), and QwWarning.

 {
   // Store pointer to running sum
   if (value.fRunningSum == 0) {
     QwVQWK_Channel* vqwk = const_cast<QwVQWK_Channel*>(&value);
     vqwk->fRunningSum = this;
   }
 
   // Moment calculations
   Bool_t berror=kTRUE;//only needed for deaccumulation (stability check purposes)
 
   /*
     note:
     The AccumulateRunningSum is called on a dedicated subsystem array object and for the standard running avg computations
     we only need value.fErrorFlag==0 events to be included in the running avg. So the "berror" conditions is only used for the stability check purposes.
 
     The need for this check below came due to fact that when routine DeaccumulateRunningSum is called the errorflag is updated with 
     the kBeamStabilityError flag (+ configuration flags for global errors) and need to make sure we remove this flag and any configuration flags before 
     checking the (fErrorFlag != 0) condition
     
     See how the stability check is implemented in the QwEventRing class
 
     Rakitha
   */
 
   if (count==-1 && value.fErrorFlag>0){
     berror=(((value.fErrorFlag) & 0xFFFFFFF) == 0); //The operation value.fErrorFlag & 0xFFFFFFF set the stability failed error bit to zero //-value.fErrorConfigFlag
     if (GetElementName()=="qwk_enegy"){//qwk_target_EffectiveCharge
       PrintValue();
     }
   }
 
   
   
   Int_t n1 = fGoodEventCount;
   Int_t n2 = count;
   // If there are no good events, check the error flag
   if (n2 == 0 && (value.fErrorFlag) == 0) {
     n2 = 1;
     //one event is removed from the sum (Deaccumulation)
   }else if (n2 == -1 && berror) { //check only single event cut errors except stability fail flag since by the time the value is deaccumulated this could be flagged as stability failed error. 
     n2 = -1;
   }else
     n2 = -100;//ignore it
   Int_t n = n1 + n2;
   // Set up variables
   Double_t M11 = fHardwareBlockSum;
   Double_t M12 = value.fHardwareBlockSum;
   Double_t M22 = value.fHardwareBlockSumM2;
   if (n2 == 0) {
     // no good events for addition
 
     return;
   } else if (n2 == -1) {
     // simple version for removal of single event from the sum
     fGoodEventCount--;
     //QwMessage<<"Deaccumulate before "<<QwLog::endl;
     if (n>0){
       fHardwareBlockSum -= (M12 - M11) / n;
       fHardwareBlockSumM2 -= (M12 - M11)
   * (M12 - fHardwareBlockSum); // note: using updated mean
       // and for individual blocks
       for (Int_t i = 0; i < 4; i++) {
   M11 = fBlock[i];
   M12 = value.fBlock[i];
   M22 = value.fBlockM2[i];
   fBlock[i] -= (M12 - M11) / n;
   fBlockM2[i] -= (M12 - M11) * (M12 - fBlock[i]); // note: using updated mean
       }
     }else if (n==0){
       //QwMessage<<"Deaccumulate at zero "<<QwLog::endl;
       /* 
       //Need any fail safe check for Deaccumulation????
       */
     }
 
   } else if (n2 == 1) {
     // simple version for addition of single event
     fGoodEventCount++;
     fHardwareBlockSum += (M12 - M11) / n;
     fHardwareBlockSumM2 += (M12 - M11)
          * (M12 - fHardwareBlockSum); // note: using updated mean
     // and for individual blocks
     for (Int_t i = 0; i < 4; i++) {
       M11 = fBlock[i];
       M12 = value.fBlock[i];
       M22 = value.fBlockM2[i];
       fBlock[i] += (M12 - M11) / n;
       fBlockM2[i] += (M12 - M11) * (M12 - fBlock[i]); // note: using updated mean
     }
   } else if (n2 > 1) {
     // general version for addition of multi-event sets
     fGoodEventCount += n2;
     fHardwareBlockSum += n2 * (M12 - M11) / n;
     fHardwareBlockSumM2 += M22 + n1 * n2 * (M12 - M11) * (M12 - M11) / n;
     // and for individual blocks
     for (Int_t i = 0; i < 4; i++) {
       M11 = fBlock[i];
       M12 = value.fBlock[i];
       M22 = value.fBlockM2[i];
       fBlock[i] += n2 * (M12 - M11) / n;
       fBlockM2[i] += M22 + n1 * n2 * (M12 - M11) * (M12 - M11) / n;
     }
   }
 
   // Nanny
   if (fHardwareBlockSum != fHardwareBlockSum)
     QwWarning << "Angry Nanny: NaN detected in " << GetElementName() << QwLog::endl;
 
 
    
 
 }

Here is the call graph for this function:

void QwVQWK_Channel::AccumulateRunningSum	(	const VQwHardwareChannel *	value,
		Int_t	count
	)

inlinevirtual

Implements VQwHardwareChannel.

Definition at line 147 of file QwVQWK_Channel.h.

References AccumulateRunningSum().

                                                                          {
     const QwVQWK_Channel *tmp_ptr = dynamic_cast<const QwVQWK_Channel*>(value);
     if (tmp_ptr != NULL) AccumulateRunningSum(*tmp_ptr, count);
   };

Here is the call graph for this function:

void QwVQWK_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 1150 of file QwVQWK_Channel.cc.

References fBlock, fBlocksPerEvent, fHardwareBlockSum, and VQwDataElement::IsNameEmpty().

Referenced by QwBPMCavity::GetAbsolutePosition(), QwLinearDiodeArray::GetAbsolutePosition(), QwEnergyCalculator::ProcessEvent(), and QwBPMCavity::ProcessEvent().

 {
   if (!IsNameEmpty()){
     fHardwareBlockSum += offset;
     for (Int_t i=0; i<fBlocksPerEvent; i++) fBlock[i] += offset;
   }
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::AddErrEntriesToList ( std::vector< QwErrDBInterface > & row_list )

virtual

Reimplemented from VQwHardwareChannel.

Definition at line 1618 of file QwVQWK_Channel.cc.

References fErrorCount_HWSat, fErrorCount_SameHW, fErrorCount_sample, fErrorCount_Sequence, fErrorCount_SW_HW, fErrorCount_ZeroHW, fNumEvtsWithEventCutsRejected, VQwDataElement::GetElementName(), QwErrDBInterface::Reset(), QwErrDBInterface::SetDeviceName(), QwErrDBInterface::SetErrorCodeId(), and QwErrDBInterface::SetN().

Referenced by QwEnergyCalculator::GetErrDBEntry(), QwQPD::GetErrDBEntry(), QwBPMCavity::GetErrDBEntry(), QwLinearDiodeArray::GetErrDBEntry(), and QwIntegrationPMT::GetErrDBEntry().

 {
 
   // TString message;
   // 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);
   // QwMessage << message << QwLog::endl;
 
   // kErrorFlag_VQWK_Sat   =0x1;    //VQWK Saturation Cut. Currently saturation limit is set to +/-8.5V
   // kErrorFlag_sample     =0x2;    //If sample size mis-matches with the default value in the map file.
   // kErrorFlag_SW_HW      =0x4;    //If software sum and hardware sum are not equal.
   // kErrorFlag_Sequence   =0x8;    //If the ADC sequence number is not incrementing properly
   // kErrorFlag_SameHW     =0x10;   //If ADC value keep returning the same value
   // kErrorFlag_ZeroHW     =0x20;   //Check to see ADC is returning zero
   
 
   
   // kErrorFlag_EventCut_L =0x40;   //Flagged if lower limit of the event cut has failed
   // kErrorFlag_EventCut_U =0x80;   //Flagged if upper limit of the event cut has failed
   // >>>>>>  fNumEvtsWithEventCutsRejected
   
   
   // outside QwVQWK_Channel
   // kErrorFlag_BlinderFail = 0x0200;// in Decimal  512 to identify the blinder fail flag
   // kStabilityCutError     = 0x10000000;// in Decimal 2^28 to identify the stability cut failure
   
   // This is my modified mysql DB, Thursday, December  8 16:40:36 EST 2011, jhlee
   // Error code must be matched to MySQL DB
   // 
   // mysql> select * from error_code;
   // +---------------+------------------------------+
   // | error_code_id | quantity                     |
   // +---------------+------------------------------+
   // |             1 | kErrorFlag_VQWK_Sat          | 
   // |             2 | kErrorFlag_sample            | 
   // |             3 | kErrorFlag_SW_HW             | 
   // |             4 | kErrorFlag_Sequence          | 
   // |             5 | kErrorFlag_SameHW            | 
   // |             6 | kErrorFlag_ZeroHW            | 
   // |             7 | kErrorFlag_EventCut_Rejected | 
   // |             8 | kErrorFlag_EventCut_L        | 
   // |             9 | kErrorFlag_EventCut_U        | 
   // |            10 | kErrorFlag_BlinderFail       | 
   // |            11 | kStabilityCutError           | 
   // +---------------+------------------------------+
   // 11 rows in set (0.00 sec)
 
 
   QwErrDBInterface row;
   TString name    = GetElementName();
   
   row.Reset();
   row.SetDeviceName(name);
   row.SetErrorCodeId(1); 
   row.SetN(fErrorCount_HWSat);
   row_list.push_back(row);
   
   row.Reset();
   row.SetDeviceName(name);
   row.SetErrorCodeId(2);
   row.SetN(fErrorCount_sample);
   row_list.push_back(row);
   
   row.Reset();
   row.SetDeviceName(name);
   row.SetErrorCodeId(3);
   row.SetN(fErrorCount_SW_HW);
   row_list.push_back(row);
   
   
   row.Reset();
   row.SetDeviceName(name);
   row.SetErrorCodeId(4);
   row.SetN(fErrorCount_Sequence);
   row_list.push_back(row);
   
   
   row.Reset();
   row.SetDeviceName(name);
   row.SetErrorCodeId(5); 
   row.SetN(fErrorCount_SameHW);
   row_list.push_back(row);
   
   row.Reset();
   row.SetDeviceName(name);
   row.SetErrorCodeId(6); 
   row.SetN(fErrorCount_ZeroHW);
   row_list.push_back(row);
 
 
   row.Reset();
   row.SetDeviceName(name);
   row.SetErrorCodeId(7); 
   row.SetN(fNumEvtsWithEventCutsRejected);
   row_list.push_back(row);
   return;
   
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::AddValueFrom ( const VQwHardwareChannel * valueptr )

virtual

Implements VQwHardwareChannel.

Definition at line 846 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

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

Here is the call graph for this function:

Int_t QwVQWK_Channel::ApplyHWChecks ( )

virtual

Implements VQwHardwareChannel.

Definition at line 54 of file QwVQWK_Channel.cc.

References bDEBUG, VQwHardwareChannel::bEVENTCUTMODE, QwLog::endl(), fADC_Same_NumEvt, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fNumberOfSamples, fNumberOfSamples_map, fPrev_HardwareBlockSum, fSequenceNo_Counter, fSequenceNo_Prev, VQwDataElement::GetElementName(), GetNumberOfSamples(), GetRawHardwareSum(), GetRawSoftwareSum(), GetSequenceNumber(), GetVQWKSaturationLimt(), kErrorFlag_SameHW, kErrorFlag_sample, kErrorFlag_Sequence, kErrorFlag_SW_HW, kErrorFlag_VQWK_Sat, kErrorFlag_ZeroHW, kVQWK_VoltsPerBit, MatchNumberOfSamples(), MatchSequenceNumber(), and QwWarning.

Referenced by QwQPD::ApplyHWChecks(), QwBPMCavity::ApplyHWChecks(), QwLinearDiodeArray::ApplyHWChecks(), and QwIntegrationPMT::ApplyHWChecks().

 {
   Bool_t fEventIsGood=kTRUE;
   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 SW and HW return the same sum
     bStatus = (GetRawHardwareSum() == GetRawSoftwareSum());
     //fEventIsGood = bStatus;
     if (!bStatus) {
       fErrorFlag |= kErrorFlag_SW_HW;
     }
 
 
 
     //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
       fEventIsGood=kFALSE;
       fErrorFlag|=kErrorFlag_Sequence;
       if (bDEBUG)       QwWarning<<" QwQWVK_Channel "<<GetElementName()<<" Sequence number  previous value = "<<fSequenceNo_Prev<<" Current value= "<< GetSequenceNumber()<<QwLog::endl;
     }
 
     fSequenceNo_Counter++;
 
     //Checking for HW_sum is returning same value.
     if (fPrev_HardwareBlockSum != GetRawHardwareSum()){
       //std::cout<<" BCM hardware sum is different  "<<std::endl;
       fPrev_HardwareBlockSum = GetRawHardwareSum();
       fADC_Same_NumEvt=0;
     }else
       fADC_Same_NumEvt++;//hw_sum is same increment the counter
 
     //check for the hw_sum is giving the same value
     if (fADC_Same_NumEvt>0){//we have ADC stuck with same value
       if (bDEBUG) QwWarning<<" BCM hardware sum is same for more than  "<<fADC_Same_NumEvt<<" time consecutively  "<<QwLog::endl;
       fErrorFlag|=kErrorFlag_SameHW;
     }
 
     //check for the hw_sum is zero
     if (GetRawHardwareSum()==0){
       fErrorFlag|=kErrorFlag_ZeroHW;
     }
     if (!fEventIsGood)    
       fSequenceNo_Counter=0;//resetting the counter after ApplyHWChecks() a failure
 
     if ((TMath::Abs(GetRawHardwareSum())*kVQWK_VoltsPerBit/fNumberOfSamples) > GetVQWKSaturationLimt()){
       if (bDEBUG) 
   QwWarning << this->GetElementName()<<" "<<GetRawHardwareSum() << "Saturating VQWK invoked! " <<TMath::Abs(GetRawHardwareSum())*kVQWK_VoltsPerBit/fNumberOfSamples<<" Limit "<<GetVQWKSaturationLimt() << QwLog::endl;
       fErrorFlag|=kErrorFlag_VQWK_Sat; 
     }
 
   }
   else {
     fGoodEventCount = 1;
     fErrorFlag = 0;
   }
 
   return fErrorFlag;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 1492 of file QwVQWK_Channel.cc.

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

Referenced by QwEnergyCalculator::ApplySingleEventCuts(), QwQPD::ApplySingleEventCuts(), QwBPMCavity::ApplySingleEventCuts(), QwLinearDiodeArray::ApplySingleEventCuts(), and QwIntegrationPMT::ApplySingleEventCuts().

 {
   Bool_t status = kFALSE;
 
   if (LL==0 && UL==0){
     status=kTRUE;
   } else  if (GetHardwareSum()<=UL && GetHardwareSum()>=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;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Bool_t QwVQWK_Channel::ApplySingleEventCuts ( )

virtual

Implements VQwHardwareChannel.

Definition at line 1508 of file QwVQWK_Channel.cc.

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

 {
   Bool_t status;
 
   if (bEVENTCUTMODE>=2){//Global switch to ON/OFF event cuts set at the event cut file
 
     if (fULimit==0 && fLLimit==0){
       status=kTRUE;
     } else  if (GetHardwareSum()<=fULimit && GetHardwareSum()>=fLLimit){
       if ((fErrorFlag)==0)
   status=kTRUE;
       else
   status=kFALSE;//If the device HW is failed
     }
     else{
       if (GetHardwareSum()> 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;
 }

Here is the call graph for this function:

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

Definition at line 811 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumError, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, and VQwDataElement::IsNameEmpty().

 {
   if(this == &value) return;
 
   if (!IsNameEmpty()) {
     for (Int_t i=0; i<fBlocksPerEvent; i++){
       this->fBlock_raw[i] = value.fBlock_raw[i];   // Keep this?
       this->fBlock[i]   = value.fBlock[i]   * scale;
       this->fBlockM2[i] = value.fBlockM2[i] * scale * scale;
     }
     this->fHardwareBlockSum_raw = value.fHardwareBlockSum_raw;  // Keep this?
     this->fSoftwareBlockSum_raw = value.fSoftwareBlockSum_raw;  // Keep this?
     this->fHardwareBlockSum   = value.fHardwareBlockSum * scale;
     this->fHardwareBlockSumM2 = value.fHardwareBlockSumM2 * scale * scale;
     this->fHardwareBlockSumError = value.fHardwareBlockSumError;   // Keep this?
     this->fGoodEventCount  = value.fGoodEventCount;
     this->fNumberOfSamples = value.fNumberOfSamples;
     this->fSequenceNumber  = value.fSequenceNumber;
     this->fErrorFlag       = value.fErrorFlag;
   }
 }

Here is the call graph for this function:

void QwVQWK_Channel::AssignValueFrom ( const VQwDataElement * valueptr )

virtual

Implements VQwHardwareChannel.

Definition at line 834 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

Referenced by QwEnergyCalculator::ProcessEvent().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::Blind ( const QwBlinder * blinder )

Blind this channel as an asymmetry.

Blind this channel as an asymmetry

Parameters

blinder Blinder

Definition at line 1427 of file QwVQWK_Channel.cc.

References QwBlinder::BlindValue(), fBlock, fBlocksPerEvent, VQwDataElement::fErrorFlag, fHardwareBlockSum, QwBlinder::IsBlinderOkay(), VQwDataElement::IsNameEmpty(), and QwBlinder::ModifyThisErrorCode().

Referenced by QwIntegrationPMT::Blind().

 {
   if (!IsNameEmpty()) {
     if (blinder->IsBlinderOkay() && ((fErrorFlag)==0) ){
       for (Int_t i = 0; i < fBlocksPerEvent; i++)
   blinder->BlindValue(fBlock[i]);
       blinder->BlindValue(fHardwareBlockSum);
     } else {
       blinder->ModifyThisErrorCode(fErrorFlag);
       for (Int_t i = 0; i < fBlocksPerEvent; i++)  fBlock[i] = 0.0;
       fHardwareBlockSum = 0.0;
     }
   }
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::Blind	(	const QwBlinder *	blinder,
		const QwVQWK_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 1448 of file QwVQWK_Channel.cc.

References QwBlinder::BlindValue(), fBlock, fBlocksPerEvent, VQwDataElement::fErrorFlag, fHardwareBlockSum, QwBlinder::IsBlinderOkay(), VQwDataElement::IsNameEmpty(), and QwBlinder::ModifyThisErrorCode().

 {
   if (!IsNameEmpty()) {
     if (blinder->IsBlinderOkay() && ((fErrorFlag) ==0) ){
       for (Int_t i = 0; i < fBlocksPerEvent; i++)
   blinder->BlindValue(fBlock[i], yield.fBlock[i]);
       blinder->BlindValue(fHardwareBlockSum, yield.fHardwareBlockSum);
     } else {
       blinder->ModifyThisErrorCode(fErrorFlag);//update the HW error code
       for (Int_t i = 0; i < fBlocksPerEvent; i++)  fBlock[i] = 0.0;
       fHardwareBlockSum = 0.0;
     }
   }
   return;
 }

Here is the call graph for this function:

void QwVQWK_Channel::CalculateRunningAverage ( )

virtual

Implements VQwHardwareChannel.

Definition at line 1347 of file QwVQWK_Channel.cc.

References fBlockError, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorConfigFlag, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fHardwareBlockSumError, fHardwareBlockSumM2, VQwHardwareChannel::fStability, VQwHardwareChannel::GetValueWidth(), kBeamStabilityError, and kStabilityCut.

Referenced by QwEnergyCalculator::CalculateRunningAverage(), QwQPD::CalculateRunningAverage(), QwBPMCavity::CalculateRunningAverage(), QwLinearDiodeArray::CalculateRunningAverage(), and QwIntegrationPMT::CalculateRunningAverage().

 {
   /*
   if (GetElementName()=="qwk_engy"){//qwk_target_EffectiveCharge
       PrintValue();
   }
   */
 
   if (fGoodEventCount <= 0)
     {
       for (Int_t i = 0; i < fBlocksPerEvent; i++) {
         fBlockError[i] = 0.0;
       }
       fHardwareBlockSumError = 0.0;
     }
   else
     {
       // We use a biased estimator by dividing by n.  Use (n - 1) to get the
       // unbiased estimator for the standard deviation.
       //
       // Note we want to calculate the error here, not sigma:
       //    sigma = sqrt(M2 / n);
       //    error = sigma / sqrt (n) = sqrt(M2) / n;
       for (Int_t i = 0; i < fBlocksPerEvent; i++)
         fBlockError[i] = sqrt(fBlockM2[i]) / fGoodEventCount;
       fHardwareBlockSumError = sqrt(fHardwareBlockSumM2) / fGoodEventCount;
       //Stability check 83951872 
       if ((fStability>0) &&( (fErrorConfigFlag & kStabilityCut)==kStabilityCut)){//check to see the channel has stability cut activated in the event cut file
   /*
     //for Debugging
     PrintValue();
   */
   if (GetValueWidth()>fStability){//if the width is greater than the stability required flag the event
     fErrorFlag=kBeamStabilityError;
   }else
     fErrorFlag=0;
       }
     
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::ClearEventData ( )

virtual

Clear the event data in this element.

Reimplemented from VQwHardwareChannel.

Definition at line 241 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockError, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumError, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, and fSoftwareBlockSum_raw.

Referenced by QwEnergyCalculator::ClearEventData(), QwBPMCavity::ClearEventData(), QwLinearDiodeArray::ClearEventData(), QwQPD::ClearEventData(), QwIntegrationPMT::ClearEventData(), InitializeChannel(), QwEnergyCalculator::ProcessEvent(), QwQPD::ProcessEvent(), and QwLinearDiodeArray::ProcessEvent().

 {
   for (Int_t i = 0; i < fBlocksPerEvent; i++) {
     fBlock_raw[i] = 0;
     fBlock[i] = 0.0;
     fBlockM2[i] = 0.0;
     fBlockError[i] = 0.0;
   }
   fHardwareBlockSum_raw = 0;
   fSoftwareBlockSum_raw = 0;
   fHardwareBlockSum   = 0.0;
   fHardwareBlockSumM2 = 0.0;
   fHardwareBlockSumError = 0.0;
   fSequenceNumber   = 0;
   fNumberOfSamples  = 0;
   fGoodEventCount   = 0;
   fErrorFlag=0;
   return;
 }

Here is the caller graph for this function:

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

virtual

Implements VQwHardwareChannel.

Definition at line 687 of file QwVQWK_Channel.cc.

References fHardwareBlockSum, VQwHardwareChannel::fTreeArrayIndex, VQwHardwareChannel::fTreeArrayNumEntries, VQwDataElement::GetElementName(), VQwDataElement::IsNameEmpty(), and kDEBUG.

Referenced by QwEnergyCalculator::ConstructBranch(), QwQPD::ConstructBranch(), QwBPMCavity::ConstructBranch(), QwLinearDiodeArray::ConstructBranch(), and QwIntegrationPMT::ConstructBranch().

 {
   if (IsNameEmpty()){
     //  This channel is not used, so skip setting up the tree.
   } else {
     TString basename = prefix + GetElementName();
     tree->Branch(basename,&fHardwareBlockSum,basename+"/D");
     if (kDEBUG){
       std::cerr << "QwVQWK_Channel::ConstructBranchAndVector: fTreeArrayIndex==" << fTreeArrayIndex
     << "; fTreeArrayNumEntries==" << fTreeArrayNumEntries
     << std::endl;
     }
   }
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::ConstructBranchAndVector	(	TTree *	tree,
		TString &	prefix,
		std::vector< Double_t > &	values
	)

virtual

Implements VQwHardwareChannel.

Definition at line 591 of file QwVQWK_Channel.cc.

References bBlock, bBlock_raw, bDevice_Error_Code, bHw_sum, bHw_sum_raw, bNum_samples, bSequence_number, QwLog::endl(), VQwHardwareChannel::fDataToSave, VQwHardwareChannel::fTreeArrayIndex, VQwHardwareChannel::fTreeArrayNumEntries, VQwDataElement::GetElementName(), VQwDataElement::GetModuleType(), VQwDataElement::GetSubsystemName(), gQwHists, VQwDataElement::IsNameEmpty(), kDEBUG, VQwDataElement::kRaw, QwHistogramHelper::MatchDeviceParamsFromList(), QwHistogramHelper::MatchVQWKElementFromList(), and QwMessage.

Referenced by QwEnergyCalculator::ConstructBranchAndVector(), QwQPD::ConstructBranchAndVector(), QwBPMCavity::ConstructBranchAndVector(), QwLinearDiodeArray::ConstructBranchAndVector(), and QwIntegrationPMT::ConstructBranchAndVector().

 {
 
   if (IsNameEmpty()){
     //  This channel is not used, so skip setting up the tree.
   } else {
     TString basename = prefix + GetElementName();
     fTreeArrayIndex  = values.size();
 
     TString list="";
 
     bHw_sum=gQwHists.MatchVQWKElementFromList(GetSubsystemName().Data(), GetModuleType().Data(), "hw_sum");
     bHw_sum_raw=gQwHists.MatchVQWKElementFromList(GetSubsystemName().Data(), GetModuleType().Data(), "hw_sum_raw");
     bBlock=gQwHists.MatchVQWKElementFromList(GetSubsystemName().Data(), GetModuleType().Data(), "block");
     bBlock_raw=gQwHists.MatchVQWKElementFromList(GetSubsystemName().Data(), GetModuleType().Data(), "block_raw");
     bNum_samples=gQwHists.MatchVQWKElementFromList(GetSubsystemName().Data(), GetModuleType().Data(), "num_samples");
     bDevice_Error_Code=gQwHists.MatchVQWKElementFromList(GetSubsystemName().Data(), GetModuleType().Data(), "Device_Error_Code");
     bSequence_number=gQwHists.MatchVQWKElementFromList(GetSubsystemName().Data(), GetModuleType().Data(), "sequence_number");
 
     if (bHw_sum){
       values.push_back(0.0);
       list += "hw_sum/D";
     }
     if (bBlock){
       values.push_back(0.0);
       list += ":block0/D";
 
       values.push_back(0.0);
       list += ":block1/D";
 
       values.push_back(0.0);
       list += ":block2/D";
 
       values.push_back(0.0);
       list += ":block3/D";
     }
 
     if (bNum_samples){
       values.push_back(0.0);
       list += ":num_samples/D";
     }
 
     if (bDevice_Error_Code){
       values.push_back(0.0);
       list += ":Device_Error_Code/D";
     }
 
     if(fDataToSave==kRaw)
       {
   if (bHw_sum_raw){
     values.push_back(0.0);
     list += ":hw_sum_raw/D";
   }
   if (bBlock_raw){
     values.push_back(0.0);
     list += ":block0_raw/D";
 
     values.push_back(0.0);
     list += ":block1_raw/D";
 
     values.push_back(0.0);
     list += ":block2_raw/D";
 
     values.push_back(0.0);
     list += ":block3_raw/D";
   }
   if (bSequence_number){
     values.push_back(0.0);
     list += ":sequence_number/D";
   }
       }
 
     fTreeArrayNumEntries = values.size() - fTreeArrayIndex;
  
     if (gQwHists.MatchDeviceParamsFromList(basename.Data()) && (bHw_sum || bBlock || bNum_samples || bDevice_Error_Code || bHw_sum_raw || bBlock_raw || bSequence_number)){
 
       if (list=="hw_sum/D")//this is for the RT mode
   list=basename+"/D"; 
       
       if (kDEBUG)
   QwMessage <<"base name "<<basename<<" List "<<list<<  QwLog::endl;
       tree->Branch(basename, &(values[fTreeArrayIndex]), list);
     }
 
     if (kDEBUG){
       std::cerr << "QwVQWK_Channel::ConstructBranchAndVector: fTreeArrayIndex==" << fTreeArrayIndex
     << "; fTreeArrayNumEntries==" << fTreeArrayNumEntries
     << "; values.size()==" << values.size()
     << "; list==" << list
     << std::endl;
     }
   }
   //exit(1);
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Construct the histograms for this data element.

Implements VQwDataElement.

Definition at line 481 of file QwVQWK_Channel.cc.

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

Referenced by QwEnergyCalculator::ConstructHistograms(), QwQPD::ConstructHistograms(), QwBPMCavity::ConstructHistograms(), QwLinearDiodeArray::ConstructHistograms(), and QwIntegrationPMT::ConstructHistograms().

 {
   //  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(8+2+1, NULL);
   size_t index=0;
   for (Int_t i=0; i<fBlocksPerEvent; i++){
     fHistograms[index]   = gQwHists.Construct1DHist(basename+Form("_block%d_raw",i));
     fHistograms[index+1] = gQwHists.Construct1DHist(basename+Form("_block%d",i));
     index += 2;
   }
   fHistograms[index]   = gQwHists.Construct1DHist(basename+Form("_hw_raw"));
   fHistograms[index+1] = gQwHists.Construct1DHist(basename+Form("_hw"));
   index += 2;
   fHistograms[index]   = gQwHists.Construct1DHist(basename+Form("_sw-hw_raw"));
       }
     else if(fDataToSave==kDerived)
       {
   fHistograms.resize(4+1+1, NULL);
   Int_t index=0;
   for (Int_t i=0; i<fBlocksPerEvent; i++){
     fHistograms[index] = gQwHists.Construct1DHist(basename+Form("_block%d",i));
     index += 1;
   }
   fHistograms[index] = gQwHists.Construct1DHist(basename+Form("_hw"));
   index += 1;
   fHistograms[index] = gQwHists.Construct1DHist(basename+Form("_dev_err"));
   index += 1;
       }
     else
       {
   // this is not recognized
       }
 
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::DeaccumulateRunningSum ( const QwVQWK_Channel & value )

inline

Definition at line 152 of file QwVQWK_Channel.h.

References AccumulateRunningSum().

Referenced by QwEnergyCalculator::DeaccumulateRunningSum(), QwQPD::DeaccumulateRunningSum(), QwBPMCavity::DeaccumulateRunningSum(), QwLinearDiodeArray::DeaccumulateRunningSum(), and QwIntegrationPMT::DeaccumulateRunningSum().

                                                                  {
     AccumulateRunningSum(value, -1);
   };

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 1042 of file QwVQWK_Channel.cc.

Referenced by QwLinearDiodeArray::ProcessEvent(), and QwQPD::ProcessEvent().

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

Here is the caller graph for this function:

void QwVQWK_Channel::DivideBy ( const VQwHardwareChannel * valueptr )

virtual

Implements VQwHardwareChannel.

Definition at line 882 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

Referenced by QwIntegrationPMT::Normalize().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::DivideBy ( const QwVQWK_Channel & denom )

Definition at line 1173 of file QwVQWK_Channel.cc.

 {
   *this /= denom;
 }

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

virtual

Encode the event data into a CODA buffer.

Implements MQwMockable.

Definition at line 331 of file QwVQWK_Channel.cc.

References fBlock_raw, fHardwareBlockSum_raw, fNumberOfSamples, fSequenceNumber, and VQwDataElement::IsNameEmpty().

Referenced by QwIntegrationPMT::EncodeEventData().

 {
   Long_t localbuf[6] = {0};
 
   if (IsNameEmpty()) {
     //  This channel is not used, but is present in the data stream.
     //  Skip over this data.
   } else {
     //    localbuf[4] = 0;
     for (Int_t i = 0; i < 4; i++) {
       localbuf[i] = fBlock_raw[i];
       //        localbuf[4] += localbuf[i]; // fHardwareBlockSum_raw
     }
     // The following causes many rounding errors and skips due to the check
     // that fHardwareBlockSum_raw == fSoftwareBlockSum_raw in IsGoodEvent().
     localbuf[4] = fHardwareBlockSum_raw;
     localbuf[5] = (fNumberOfSamples << 16 & 0xFFFF0000)
                 | (fSequenceNumber  << 8  & 0x0000FF00);
 
     for (Int_t i = 0; i < 6; i++){
         buffer.push_back(localbuf[i]);
     }
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::FillHistograms ( )

virtual

Fill the histograms for this data element.

Implements VQwDataElement.

Definition at line 533 of file QwVQWK_Channel.cc.

References fBlocksPerEvent, VQwHardwareChannel::fDataToSave, VQwDataElement::fErrorFlag, MQwHistograms::fHistograms, GetBlockValue(), GetHardwareSum(), GetRawBlockValue(), GetRawHardwareSum(), GetRawSoftwareSum(), VQwDataElement::IsNameEmpty(), VQwDataElement::kDerived, kErrorFlag_SameHW, kErrorFlag_sample, kErrorFlag_Sequence, kErrorFlag_SW_HW, kErrorFlag_VQWK_Sat, kErrorFlag_ZeroHW, and VQwDataElement::kRaw.

Referenced by QwEnergyCalculator::FillHistograms(), QwQPD::FillHistograms(), QwBPMCavity::FillHistograms(), QwLinearDiodeArray::FillHistograms(), and QwIntegrationPMT::FillHistograms().

 {
   Int_t index=0;
   
   if (IsNameEmpty())
     {
       //  This channel is not used, so skip creating the histograms.
     } else
       {
   if(fDataToSave==kRaw)
     {
       for (Int_t i=0; i<fBlocksPerEvent; i++)
         {
     if (fHistograms[index] != NULL && (fErrorFlag)==0)
       fHistograms[index]->Fill(this->GetRawBlockValue(i));
     if (fHistograms[index+1] != NULL && (fErrorFlag)==0)
       fHistograms[index+1]->Fill(this->GetBlockValue(i));
     index+=2;
         }
       if (fHistograms[index] != NULL && (fErrorFlag)==0)
         fHistograms[index]->Fill(this->GetRawHardwareSum());
       if (fHistograms[index+1] != NULL && (fErrorFlag)==0)
         fHistograms[index+1]->Fill(this->GetHardwareSum());
       index+=2;
       if (fHistograms[index] != NULL && (fErrorFlag)==0)
         fHistograms[index]->Fill(this->GetRawSoftwareSum()-this->GetRawHardwareSum());
     }
   else if(fDataToSave==kDerived)
     {
       for (Int_t i=0; i<fBlocksPerEvent; i++)
         {
     if (fHistograms[index] != NULL && (fErrorFlag)==0)
       fHistograms[index]->Fill(this->GetBlockValue(i));
     index+=1;
         }
       if (fHistograms[index] != NULL && (fErrorFlag)==0)
         fHistograms[index]->Fill(this->GetHardwareSum());
       index+=1; 
       if (fHistograms[index] != NULL){
         if ( (kErrorFlag_sample &  fErrorFlag)==kErrorFlag_sample)
     fHistograms[index]->Fill(kErrorFlag_sample);
         if ( (kErrorFlag_SW_HW &  fErrorFlag)==kErrorFlag_SW_HW)
     fHistograms[index]->Fill(kErrorFlag_SW_HW);
         if ( (kErrorFlag_Sequence &  fErrorFlag)==kErrorFlag_Sequence)
     fHistograms[index]->Fill(kErrorFlag_Sequence);
         if ( (kErrorFlag_ZeroHW &  fErrorFlag)==kErrorFlag_ZeroHW)
     fHistograms[index]->Fill(kErrorFlag_ZeroHW);
         if ( (kErrorFlag_VQWK_Sat &  fErrorFlag)==kErrorFlag_VQWK_Sat)
     fHistograms[index]->Fill(kErrorFlag_VQWK_Sat);
         if ( (kErrorFlag_SameHW &  fErrorFlag)==kErrorFlag_SameHW)
     fHistograms[index]->Fill(kErrorFlag_SameHW);
       }
       
     }
  
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::FillTreeVector ( std::vector< Double_t > & values ) const

virtual

Implements VQwHardwareChannel.

Definition at line 704 of file QwVQWK_Channel.cc.

References bBlock, bBlock_raw, bDEBUG, bDevice_Error_Code, bHw_sum, bHw_sum_raw, bNum_samples, bSequence_number, fBlocksPerEvent, VQwHardwareChannel::fDataToSave, VQwDataElement::fErrorFlag, fNumberOfSamples, fSequenceNumber, VQwHardwareChannel::fTreeArrayIndex, VQwHardwareChannel::fTreeArrayNumEntries, GetBlockValue(), GetHardwareSum(), GetRawBlockValue(), GetRawHardwareSum(), VQwDataElement::IsNameEmpty(), and VQwDataElement::kRaw.

Referenced by QwEnergyCalculator::FillTreeVector(), QwBPMCavity::FillTreeVector(), QwQPD::FillTreeVector(), QwLinearDiodeArray::FillTreeVector(), and QwIntegrationPMT::FillTreeVector().

 {
   if (IsNameEmpty()) {
     //  This channel is not used, so skip filling the tree vector.
   } else if (fTreeArrayNumEntries <= 0) {
     if (bDEBUG) std::cerr << "QwVQWK_Channel::FillTreeVector:  fTreeArrayNumEntries=="
         << fTreeArrayNumEntries << std::endl;
   } else if (values.size() < fTreeArrayIndex+fTreeArrayNumEntries){
     if (bDEBUG) std::cerr << "QwVQWK_Channel::FillTreeVector:  values.size()=="
         << values.size()
         << "; fTreeArrayIndex+fTreeArrayNumEntries=="
         << fTreeArrayIndex+fTreeArrayNumEntries
         << std::endl;
   } else {
     UInt_t index = fTreeArrayIndex;
     //hw_sum
     if (bHw_sum)
       values[index++] = this->GetHardwareSum();
 
     if (bBlock){
       for (Int_t i=0; i<fBlocksPerEvent; i++){
   //blocki
   values[index++] = this->GetBlockValue(i);
       }
     }
 
     //num_samples
     if (bNum_samples)
       values[index++] = this->fNumberOfSamples;
 
     //Device_Error_Code
     if (bDevice_Error_Code)
       values[index++] = (this->fErrorFlag);
 
     if(fDataToSave==kRaw)
       {
   //hw_sum_raw
   if (bHw_sum_raw)
     values[index++] = this->GetRawHardwareSum();
 
   if (bBlock_raw){
     for (Int_t i=0; i<fBlocksPerEvent; i++){
       //blocki_raw
       values[index++] = this->GetRawBlockValue(i);
     }
   }
 
   //sequence_number
   if (bSequence_number)
     values[index++]=this->fSequenceNumber;
       }
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Double_t QwVQWK_Channel::GetAverageVolts ( ) const

Definition at line 442 of file QwVQWK_Channel.cc.

References fHardwareBlockSum, fNumberOfSamples, and kVQWK_VoltsPerBit.

 {
   //Double_t avgVolts = (fBlock[0]+fBlock[1]+fBlock[2]+fBlock[3])*kVQWK_VoltsPerBit/fNumberOfSamples;
   Double_t avgVolts = fHardwareBlockSum * kVQWK_VoltsPerBit / fNumberOfSamples;
   //std::cout<<"QwVQWK_Channel::GetAverageVolts() = "<<avgVolts<<std::endl;
   return avgVolts;
 
 }

Double_t QwVQWK_Channel::GetBlockErrorValue ( size_t blocknum ) const

inlineprivate

Definition at line 245 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetValueError().

Referenced by PrintValue().

245 { return GetValueError(blocknum+1);};

VQwHardwareChannel::GetValueError

Double_t GetValueError() const

Definition: VQwHardwareChannel.h:53

Here is the call graph for this function:

Here is the caller graph for this function:

Double_t QwVQWK_Channel::GetBlockValue ( size_t blocknum ) const

inlineprivate

Definition at line 244 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetValue().

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

244 { return GetValue(blocknum+1);};

VQwHardwareChannel::GetValue

Double_t GetValue() const

Definition: VQwHardwareChannel.h:51

Here is the call graph for this function:

Here is the caller graph for this function:

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

static

Class: QwVQWK_Channel Base class containing decoding functions for the VQWK_Channel 6 32-bit datawords. The functions in this class will decode a single channel worth of VQWK_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 zero

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

Definition at line 33 of file QwVQWK_Channel.cc.

References QwLog::endl(), kMaxChannels, kWordsPerChannel, and QwError.

Referenced by QwMainCerenkovDetector::LoadChannelMap(), QwIntegratedRaster::LoadChannelMap(), QwLumi::LoadChannelMap(), QwScanner::LoadChannelMap(), and QwBeamDetectorID::QwBeamDetectorID().

                                                                           {
     Int_t offset = -1;
     if (moduleindex<0 ){
       QwError << "QwVQWK_Channel::GetBufferOffset:  Invalid module index,"
         << moduleindex
         << ".  Must be zero or greater."
         << QwLog::endl;
     } else if (channelindex<0 || channelindex>kMaxChannels){
       QwError << "QwVQWK_Channel::GetBufferOffset:  Invalid channel index,"
         << channelindex
         << ".  Must be in range [0," << kMaxChannels << "]."
         << QwLog::endl;
     } else {
       offset = ( (moduleindex * kMaxChannels) + channelindex )
   * kWordsPerChannel;
     }
     return offset;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

Double_t QwVQWK_Channel::GetHardwareSum ( ) const

inlineprivate

Definition at line 247 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetValue().

Referenced by ApplySingleEventCuts(), FillHistograms(), FillTreeVector(), operator<<(), and PrintValue().

247 { return GetValue(0);};

VQwHardwareChannel::GetValue

Double_t GetValue() const

Definition: VQwHardwareChannel.h:51

Here is the call graph for this function:

Here is the caller graph for this function:

Double_t QwVQWK_Channel::GetHardwareSumError ( ) const

inlineprivate

Definition at line 250 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetValueError().

Referenced by PrintValue().

250 { return GetValueError(0); };

VQwHardwareChannel::GetValueError

Double_t GetValueError() const

Definition: VQwHardwareChannel.h:53

Here is the call graph for this function:

Here is the caller graph for this function:

Double_t QwVQWK_Channel::GetHardwareSumM2 ( ) const

inlineprivate

Definition at line 248 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetValueM2().

248 { return GetValueM2(0); };

VQwHardwareChannel::GetValueM2

Double_t GetValueM2() const

Definition: VQwHardwareChannel.h:52

Here is the call graph for this function:

Double_t QwVQWK_Channel::GetHardwareSumWidth ( ) const

inlineprivate

Definition at line 249 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetValueWidth().

Referenced by PrintValue().

249 { return GetValueWidth(0); };

VQwHardwareChannel::GetValueWidth

Double_t GetValueWidth() const

Definition: VQwHardwareChannel.h:54

Here is the call graph for this function:

Here is the caller graph for this function:

size_t QwVQWK_Channel::GetNumberOfSamples ( ) const

inline

Definition at line 219 of file QwVQWK_Channel.h.

References fNumberOfSamples.

Referenced by QwQPD::ApplyHWChecks(), QwBPMCavity::ApplyHWChecks(), QwLinearDiodeArray::ApplyHWChecks(), and ApplyHWChecks().

219 {return (fNumberOfSamples);};

QwVQWK_Channel::fNumberOfSamples

size_t fNumberOfSamples

Number of samples read through the module.

Definition: QwVQWK_Channel.h:307

Here is the caller graph for this function:

Int_t QwVQWK_Channel::GetRawBlockValue ( size_t blocknum ) const

inlineprivate

Definition at line 253 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetRawValue().

Referenced by FillHistograms(), and FillTreeVector().

253 {return GetRawValue(blocknum+1);};

VQwHardwareChannel::GetRawValue

Int_t GetRawValue() const

Definition: VQwHardwareChannel.h:50

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwVQWK_Channel::GetRawHardwareSum ( ) const

inlineprivate

Definition at line 254 of file QwVQWK_Channel.h.

References VQwHardwareChannel::GetRawValue().

Referenced by ApplyHWChecks(), FillHistograms(), and FillTreeVector().

254 { return GetRawValue(0);};

VQwHardwareChannel::GetRawValue

Int_t GetRawValue() const

Definition: VQwHardwareChannel.h:50

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwVQWK_Channel::GetRawSoftwareSum ( ) const

inlineprivate

Definition at line 255 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), and FillHistograms().

255 {return fSoftwareBlockSum_raw;};

QwVQWK_Channel::fSoftwareBlockSum_raw

Int_t fSoftwareBlockSum_raw

Sum of the data in the four sub-blocks raw.

Definition: QwVQWK_Channel.h:285

Here is the caller graph for this function:

Int_t QwVQWK_Channel::GetRawValue ( size_t element ) const

inlinevirtual

Implements VQwHardwareChannel.

Definition at line 194 of file QwVQWK_Channel.h.

References fBlock_raw, fHardwareBlockSum_raw, and VQwHardwareChannel::RangeCheck().

                                           {
     RangeCheck(element);
     if (element==0) return fHardwareBlockSum_raw;
     return fBlock_raw[element-1];
   }

Here is the call graph for this function:

size_t QwVQWK_Channel::GetSequenceNumber ( ) const

inline

Definition at line 218 of file QwVQWK_Channel.h.

References fSequenceNumber.

Referenced by QwQPD::ApplyHWChecks(), QwBPMCavity::ApplyHWChecks(), QwLinearDiodeArray::ApplyHWChecks(), ApplyHWChecks(), QwLinearDiodeArray::ProcessEvent(), and QwQPD::ProcessEvent().

218 {return (fSequenceNumber);};

QwVQWK_Channel::fSequenceNumber

size_t fSequenceNumber

Event sequence number for this channel.

Definition: QwVQWK_Channel.h:305

Here is the caller graph for this function:

Double_t QwVQWK_Channel::GetValue ( size_t element ) const

inlinevirtual

Implements VQwHardwareChannel.

Definition at line 199 of file QwVQWK_Channel.h.

References fBlock, fHardwareBlockSum, and VQwHardwareChannel::RangeCheck().

Referenced by QwIntegrationPMT::GetValue(), QwEnergyCalculator::ProcessEvent(), QwLinearDiodeArray::ProcessEvent(), QwQPD::ProcessEvent(), and QwBlinder::Update().

                                           {
     RangeCheck(element);
     if (element==0) return fHardwareBlockSum;
     return fBlock[element-1];
   }

Here is the call graph for this function:

Here is the caller graph for this function:

Double_t QwVQWK_Channel::GetValueError ( size_t element ) const

inlinevirtual

Implements VQwHardwareChannel.

Definition at line 209 of file QwVQWK_Channel.h.

References fBlockError, fHardwareBlockSumError, and VQwHardwareChannel::RangeCheck().

                                                {
     RangeCheck(element);
     if (element==0) return fHardwareBlockSumError;
     return fBlockError[element-1];
   }

Here is the call graph for this function:

Double_t QwVQWK_Channel::GetValueM2 ( size_t element ) const

inlinevirtual

Implements VQwHardwareChannel.

Definition at line 204 of file QwVQWK_Channel.h.

References fBlockM2, fHardwareBlockSumM2, and VQwHardwareChannel::RangeCheck().

                                             {
     RangeCheck(element);
     if (element==0) return fHardwareBlockSumM2;
     return fBlockM2[element-1];
   }

Here is the call graph for this function:

Double_t QwVQWK_Channel::GetVQWKSaturationLimt ( )

inline

Definition at line 176 of file QwVQWK_Channel.h.

References fSaturationABSLimit.

Referenced by ApplyHWChecks().

                                   {//Get the absolute staturation limit in volts.
     return fSaturationABSLimit;
   }

Here is the caller graph for this function:

void QwVQWK_Channel::IncrementErrorCounters ( )

Definition at line 130 of file QwVQWK_Channel.cc.

References fErrorCount_HWSat, 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, kErrorFlag_VQWK_Sat, and kErrorFlag_ZeroHW.

Referenced by QwEnergyCalculator::IncrementErrorCounters(), QwQPD::IncrementErrorCounters(), QwBPMCavity::IncrementErrorCounters(), QwIntegrationPMT::IncrementErrorCounters(), and QwLinearDiodeArray::IncrementErrorCounters().

                                            {
   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_VQWK_Sat &  fErrorFlag)==kErrorFlag_VQWK_Sat)
     fErrorCount_HWSat++; //increment the hw saturation error counter
   if ( ((kErrorFlag_EventCut_L &  fErrorFlag)==kErrorFlag_EventCut_L) 
        || ((kErrorFlag_EventCut_U &  fErrorFlag)==kErrorFlag_EventCut_U)){
     fNumEvtsWithEventCutsRejected++; //increment the event cut error counter
   }
 }

Here is the caller graph for this function:

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

virtual

Initialize the fields in this object.

Implements VQwHardwareChannel.

Definition at line 151 of file QwVQWK_Channel.cc.

Referenced by QwEnergyCalculator::InitializeChannel(), QwQPD::InitializeChannel(), QwIntegrationPMT::InitializeChannel(), QwBPMCavity::InitializeChannel(), QwLinearDiodeArray::InitializeChannel(), InitializeChannel(), QwBeamMod::LoadChannelMap(), QwEnergyCalculator::ProcessEvent(), QwLinearDiodeArray::ProcessEvent(), QwQPD::ProcessEvent(), and QwVQWK_Channel().

 {
   SetElementName(name);
   SetDataToSave(datatosave);
   SetNumberOfDataWords(6);
   SetNumberOfSubElements(5);
 
   kFoundPedestal = 0;
   kFoundGain = 0;
 
   fPedestal            = 0.0;
   fCalibrationFactor   = 1.0;
 
   fBlocksPerEvent      = 4;
 
   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=0;
   fLLimit=0;
   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;
   fPrev_HardwareBlockSum = 0.0;
 
   fGoodEventCount        = 0;
 
   bEVENTCUTMODE          = 0;
 
   //std::cout<< "name = "<<name<<" error count same _HW = "<<fErrorCount_SameHW <<std::endl;
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Initialize the fields in this object.

Implements VQwHardwareChannel.

Definition at line 221 of file QwVQWK_Channel.cc.

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

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

Here is the call graph for this function:

void QwVQWK_Channel::LoadChannelParameters ( QwParameterFile & paramfile )

virtual

Reimplemented from VQwDataElement.

Definition at line 228 of file QwVQWK_Channel.cc.

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

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

Here is the call graph for this function:

Bool_t QwVQWK_Channel::MatchNumberOfSamples ( size_t numsamp )

Definition at line 1474 of file QwVQWK_Channel.cc.

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

Referenced by ApplyHWChecks().

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

Here is the call graph for this function:

Here is the caller graph for this function:

Bool_t QwVQWK_Channel::MatchSequenceNumber ( size_t seqnum )

Definition at line 1464 of file QwVQWK_Channel.cc.

References fSequenceNumber, and VQwDataElement::IsNameEmpty().

Referenced by ApplyHWChecks().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::MultiplyBy ( const VQwHardwareChannel * valueptr )

virtual

Implements VQwHardwareChannel.

Definition at line 870 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

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

Here is the call graph for this function:

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

Definition at line 951 of file QwVQWK_Channel.cc.

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

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

Definition at line 958 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorFlag, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()){
     for (Int_t i=0; i<fBlocksPerEvent; i++){
       this->fBlock[i] *= value.fBlock[i];
       this->fBlock_raw[i] *= value.fBlock_raw[i];
       this->fBlockM2[i] = 0.0;
     }
     this->fHardwareBlockSum_raw *= value.fHardwareBlockSum_raw;
     this->fSoftwareBlockSum_raw *= value.fSoftwareBlockSum_raw;
     this->fHardwareBlockSum     *= value.fHardwareBlockSum;
     this->fHardwareBlockSumM2    = 0.0;
     this->fNumberOfSamples      *= value.fNumberOfSamples;
     this->fSequenceNumber        = 0;
     this->fErrorFlag            |= (value.fErrorFlag);
   }
 
   return *this;
 }

Here is the call graph for this function:

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

virtual

Implements VQwHardwareChannel.

Definition at line 1006 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

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

Here is the call graph for this function:

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

Definition at line 896 of file QwVQWK_Channel.cc.

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

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

Definition at line 903 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorFlag, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (Int_t i = 0; i < fBlocksPerEvent; i++) {
       this->fBlock[i] += value.fBlock[i];
       this->fBlock_raw[i] += value.fBlock_raw[i];
       this->fBlockM2[i] = 0.0;
     }
     this->fHardwareBlockSum_raw = value.fHardwareBlockSum_raw;
     this->fSoftwareBlockSum_raw = value.fSoftwareBlockSum_raw;
     this->fHardwareBlockSum    += value.fHardwareBlockSum;
     this->fHardwareBlockSumM2   = 0.0;
     this->fNumberOfSamples     += value.fNumberOfSamples;
     this->fSequenceNumber       = 0;
     this->fErrorFlag            |= (value.fErrorFlag);
 
   }
 
   return *this;
 }

Here is the call graph for this function:

VQwHardwareChannel & QwVQWK_Channel::operator+= ( const VQwHardwareChannel * input )

virtual

Implements VQwHardwareChannel.

Definition at line 978 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

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

Here is the call graph for this function:

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

Definition at line 924 of file QwVQWK_Channel.cc.

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

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

Definition at line 931 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorFlag, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()){
     for (Int_t i=0; i<fBlocksPerEvent; i++){
       this->fBlock[i] -= value.fBlock[i];
       this->fBlock_raw[i] = 0;
       this->fBlockM2[i] = 0.0;
     }
     this->fHardwareBlockSum_raw = 0;
     this->fSoftwareBlockSum_raw = 0;
     this->fHardwareBlockSum    -= value.fHardwareBlockSum;
     this->fHardwareBlockSumM2   = 0.0;
     this->fNumberOfSamples     += value.fNumberOfSamples;
     this->fSequenceNumber       = 0;
     this->fErrorFlag           |= (value.fErrorFlag);
   }
 
   return *this;
 }

Here is the call graph for this function:

VQwHardwareChannel & QwVQWK_Channel::operator-= ( const VQwHardwareChannel * input )

virtual

Implements VQwHardwareChannel.

Definition at line 992 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

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

Here is the call graph for this function:

VQwHardwareChannel & QwVQWK_Channel::operator/= ( const VQwHardwareChannel * input )

virtual

Implements VQwHardwareChannel.

Definition at line 1020 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

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

Here is the call graph for this function:

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

protected

Definition at line 1066 of file QwVQWK_Channel.cc.

References QwLog::endl(), fBlock, fBlockM2, VQwDataElement::fErrorFlag, fHardwareBlockSum, fHardwareBlockSumM2, VQwDataElement::GetElementName(), VQwDataElement::IsNameEmpty(), QwVerbose, and QwWarning.

 {
   //  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!
     //
     for (Int_t i = 0; i < 4; i++) {
       if (this->fBlock[i] != 0.0 && denom.fBlock[i] != 0.0){
   ratio    = (this->fBlock[i]) / (denom.fBlock[i]);
         variance =  ratio * ratio *
            (this->fBlockM2[i] / this->fBlock[i] / this->fBlock[i]
           + denom.fBlockM2[i] / denom.fBlock[i] / denom.fBlock[i]);
   fBlock[i]   = ratio;
   fBlockM2[i] = variance;
       } else if (this->fBlock[i] == 0.0) {
   this->fBlock[i]   = 0.0;
         this->fBlockM2[i] = 0.0;
       } else {
         QwVerbose << "Attempting to divide by zero block in " 
       << GetElementName() << QwLog::endl;
   fBlock[i]   = 0.0;
         fBlockM2[i] = 0.0;
       }
     }
     if (this->fHardwareBlockSum != 0.0 && denom.fHardwareBlockSum != 0.0){
       ratio    =  (this->fHardwareBlockSum) / (denom.fHardwareBlockSum);
       variance =  ratio * ratio *
   (this->fHardwareBlockSumM2 / this->fHardwareBlockSum / this->fHardwareBlockSum
    + denom.fHardwareBlockSumM2 / denom.fHardwareBlockSum / denom.fHardwareBlockSum);
       fHardwareBlockSum   = ratio;
       fHardwareBlockSumM2 = variance;
     } else if (this->fHardwareBlockSum == 0.0) {
       fHardwareBlockSum   = 0.0;
       fHardwareBlockSumM2 = 0.0;
     } else {
       QwVerbose << "Attempting to divide by zero sum in " 
     << GetElementName() << QwLog::endl;
       fHardwareBlockSumM2 = 0.0;
     }
     // Remaining variables
     //  Don't change fNumberOfSamples, fSequenceNumber, fGoodEventCount,
     //  'OR' the HW error codes in the fErrorFlag values together.
     fErrorFlag |= (denom.fErrorFlag);//mix only the hardware error codes
   }
 
   // Nanny
   if (fHardwareBlockSum != fHardwareBlockSum)
     QwWarning << "Angry Nanny: NaN detected in " << GetElementName() << QwLog::endl;
 
   return *this;
 }

Here is the call graph for this function:

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

Definition at line 787 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumError, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, VQwDataElement::IsNameEmpty(), and VQwHardwareChannel::operator=().

 {
   if(this ==&value) return *this;
 
   if (!IsNameEmpty()) {
     VQwHardwareChannel::operator=(value);
     for (Int_t i=0; i<fBlocksPerEvent; i++){
       this->fBlock_raw[i] = value.fBlock_raw[i];
       this->fBlock[i]     = value.fBlock[i];
       this->fBlockM2[i]   = value.fBlockM2[i];
     }
     this->fHardwareBlockSum_raw = value.fHardwareBlockSum_raw;
     this->fSoftwareBlockSum_raw = value.fSoftwareBlockSum_raw;
     this->fHardwareBlockSum = value.fHardwareBlockSum;
     this->fHardwareBlockSumM2 = value.fHardwareBlockSumM2;
     this->fHardwareBlockSumError = value.fHardwareBlockSumError;
     this->fNumberOfSamples = value.fNumberOfSamples;
     this->fSequenceNumber  = value.fSequenceNumber;
    
 
   }
   return *this;
 }

Here is the call graph for this function:

void QwVQWK_Channel::PrintErrorCounterHead ( )

static

Definition at line 1544 of file QwVQWK_Channel.cc.

References QwLog::endl(), and QwMessage.

Referenced by QwBeamLine::PrintErrorCounters(), QwMainCerenkovDetector::PrintErrorCounters(), QwIntegratedRaster::PrintErrorCounters(), QwBeamMod::PrintErrorCounters(), and QwLumi::PrintErrorCounters().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::PrintErrorCounters ( ) const

virtual

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

Reimplemented from VQwDataElement.

Definition at line 1567 of file QwVQWK_Channel.cc.

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.

Referenced by QwIntegrationPMT::PrintErrorCounters(), QwEnergyCalculator::PrintErrorCounters(), QwQPD::PrintErrorCounters(), QwBPMCavity::PrintErrorCounters(), and QwLinearDiodeArray::PrintErrorCounters().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::PrintErrorCounterTail ( )

static

Definition at line 1561 of file QwVQWK_Channel.cc.

References QwLog::endl(), and QwMessage.

Referenced by QwBeamLine::PrintErrorCounters(), QwMainCerenkovDetector::PrintErrorCounters(), QwIntegratedRaster::PrintErrorCounters(), QwBeamMod::PrintErrorCounters(), and QwLumi::PrintErrorCounters().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::PrintInfo ( ) const

virtual

Print multiple lines of information about this data element.

Reimplemented from VQwDataElement.

Definition at line 451 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlocksPerEvent, VQwHardwareChannel::fCalibrationFactor, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumError, fHardwareBlockSumM2, fNumberOfSamples, VQwHardwareChannel::fPedestal, fSequenceNumber, fSoftwareBlockSum_raw, VQwDataElement::GetElementName(), VQwDataElement::GetModuleType(), and VQwDataElement::GetSubsystemName().

Referenced by QwEnergyCalculator::PrintInfo(), QwBPMCavity::PrintInfo(), QwLinearDiodeArray::PrintInfo(), QwIntegrationPMT::PrintInfo(), QwBPMCavity::ProcessEvent(), and QwQPD::ProcessEvent().

 {
   std::cout<<"***************************************"<<"\n";
   std::cout<<"Subsystem "<<GetSubsystemName()<<"\n"<<"\n";
   std::cout<<"Beam Instrument Type: "<<GetModuleType()<<"\n"<<"\n";
   std::cout<<"QwVQWK channel: "<<GetElementName()<<"\n"<<"\n";
   std::cout<<"fPedestal= "<< fPedestal<<"\n";
   std::cout<<"fCalibrationFactor= "<<fCalibrationFactor<<"\n";
   std::cout<<"fBlocksPerEvent= "<<fBlocksPerEvent<<"\n"<<"\n";
   std::cout<<"fSequenceNumber= "<<fSequenceNumber<<"\n";
   std::cout<<"fNumberOfSamples= "<<fNumberOfSamples<<"\n";
   std::cout<<"fBlock_raw ";
 
   for (Int_t i = 0; i < fBlocksPerEvent; i++)
     std::cout << " : " << fBlock_raw[i];
   std::cout<<"\n";
   std::cout<<"fHardwareBlockSum_raw= "<<fHardwareBlockSum_raw<<"\n";
   std::cout<<"fSoftwareBlockSum_raw= "<<fSoftwareBlockSum_raw<<"\n";
   std::cout<<"fBlock ";
   for (Int_t i = 0; i < fBlocksPerEvent; i++)
     std::cout << " : " << fBlock[i];
   std::cout << std::endl;
 
   std::cout << "fHardwareBlockSum = "<<fHardwareBlockSum << std::endl;
   std::cout << "fHardwareBlockSumM2 = "<<fHardwareBlockSumM2 << std::endl;
   std::cout << "fHardwareBlockSumError = "<<fHardwareBlockSumError << std::endl;
 
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::PrintValue ( ) const

virtual

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

Reimplemented from VQwDataElement.

Definition at line 1389 of file QwVQWK_Channel.cc.

References QwLog::endl(), VQwDataElement::fGoodEventCount, GetBlockErrorValue(), GetBlockValue(), VQwDataElement::GetElementName(), VQwDataElement::GetGoodEventCount(), GetHardwareSum(), GetHardwareSumError(), GetHardwareSumWidth(), VQwDataElement::GetModuleType(), VQwDataElement::GetSubsystemName(), and QwMessage.

Referenced by AccumulateRunningSum(), QwEnergyCalculator::PrintValue(), QwBPMCavity::PrintValue(), QwLinearDiodeArray::PrintValue(), and QwIntegrationPMT::PrintValue().

 {
   QwMessage << std::setprecision(4)
             << std::setw(18) << std::left << GetSubsystemName()      << " "
             << std::setw(18) << std::left << GetModuleType()      << " "
             << std::setw(18) << std::left << GetElementName()      << " "
             << std::setw(12) << std::left << GetHardwareSum()      << "+/- "
             << std::setw(12) << std::left << GetHardwareSumError() << " sig "
       << std::setw(12) << std::left << GetHardwareSumWidth() << " "
             << std::setw(10) << std::left << GetGoodEventCount()   << " "
             << std::setw(12) << std::left << GetBlockValue(0)      << "+/- "
             << std::setw(12) << std::left << GetBlockErrorValue(0) << " "
             << std::setw(12) << std::left << GetBlockValue(1)      << "+/- "
             << std::setw(12) << std::left << GetBlockErrorValue(1) << " "
             << std::setw(12) << std::left << GetBlockValue(2)      << "+/- "
             << std::setw(12) << std::left << GetBlockErrorValue(2) << " "
             << std::setw(12) << std::left << GetBlockValue(3)      << "+/- "
             << std::setw(12) << std::left << GetBlockErrorValue(3) << " "
             << std::setw(12) << std::left << fGoodEventCount << " "
             << QwLog::endl;
   /*
     //for Debudding
             << std::setw(12) << std::left << fErrorFlag << " err "
             << std::setw(12) << std::left << fErrorConfigFlag << " c-err "
 
   */
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwVQWK_Channel::ProcessEvBuffer	(	UInt_t *	buffer,
		UInt_t	num_words_left,
		UInt_t	index = `0`
	)

virtual

Decode the event data from a CODA buffer.

Implements VQwDataElement.

Definition at line 358 of file QwVQWK_Channel.cc.

References fBlock_raw, fBlocksPerEvent, fHardwareBlockSum_raw, VQwHardwareChannel::fNumberOfDataWords, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, VQwDataElement::IsNameEmpty(), and kWordsPerChannel.

Referenced by QwIntegrationPMT::ProcessEvBuffer(), QwBPMCavity::ProcessEvBuffer(), QwLinearDiodeArray::ProcessEvBuffer(), and QwQPD::ProcessEvBuffer().

 {
   UInt_t words_read = 0;
   UInt_t localbuf[kWordsPerChannel] = {0};
   // The conversion from UInt_t to Double_t discards the sign, so we need an intermediate
   // static_cast from UInt_t to Int_t.
   Int_t localbuf_signed[kWordsPerChannel] = {0};
 
   if (IsNameEmpty()){
     //  This channel is not used, but is present in the data stream.
     //  Skip over this data.
     words_read = fNumberOfDataWords;
   } else if (num_words_left >= fNumberOfDataWords)
     {
       for (Int_t i=0; i<kWordsPerChannel; i++){
   localbuf[i] = buffer[i];
         localbuf_signed[i] = static_cast<Int_t>(localbuf[i]);
       }
 
       fSoftwareBlockSum_raw = 0;
       for (Int_t i=0; i<fBlocksPerEvent; i++){
   fBlock_raw[i] = localbuf_signed[i];
   fSoftwareBlockSum_raw += fBlock_raw[i];
       }
       fHardwareBlockSum_raw = localbuf_signed[4];
 
       /*  Permanent change in the structure of the 6th word of the ADC readout.
        *  The upper 16 bits are the number of samples, and the upper 8 of the
        *  lower 16 are the sequence number.  This matches the structure of
        *  the ADC readout in block read mode, and now also in register read mode.
        *  P.King, 2007sep04.
        */
       fSequenceNumber   = (localbuf[5]>>8)  & 0xFF;
       fNumberOfSamples  = (localbuf[5]>>16) & 0xFFFF;
 
       words_read = fNumberOfDataWords;
 
     } else
       {
   std::cerr << "QwVQWK_Channel::ProcessEvBuffer: Not enough words!"
       << std::endl;
       }
   return words_read;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::ProcessEvent ( )

Process the event data according to pedestal and calibration factor.

Definition at line 405 of file QwVQWK_Channel.cc.

References QwLog::endl(), fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwHardwareChannel::fCalibrationFactor, VQwDataElement::fErrorFlag, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumM2, fNumberOfSamples, VQwHardwareChannel::fPedestal, VQwDataElement::GetElementName(), kErrorFlag_sample, and QwWarning.

Referenced by QwBPMCavity::ProcessEvent(), QwQPD::ProcessEvent(), QwLinearDiodeArray::ProcessEvent(), and QwIntegrationPMT::ProcessEvent().

 {
   if (fNumberOfSamples == 0 && fHardwareBlockSum_raw == 0) {
     //  There isn't valid data for this channel.  Just flag it and
     //  move on.
     for (Int_t i = 0; i < fBlocksPerEvent; i++) {
       fBlock[i] = 0.0;
       fBlockM2[i] = 0.0;
     }
     fHardwareBlockSum = 0.0;
     fHardwareBlockSumM2 = 0.0;
     fErrorFlag |= kErrorFlag_sample;
   } else if (fNumberOfSamples == 0) {
     //  This is probably a more serious problem.
     QwWarning << "QwVQWK_Channel::ProcessEvent:  Channel "
         << this->GetElementName().Data()
         << " has fNumberOfSamples==0 but has valid data in the hardware sum.  "
         << "Flag this as an error."
         << QwLog::endl;
     for (Int_t i = 0; i < fBlocksPerEvent; i++) {
       fBlock[i] = 0.0;
       fBlockM2[i] = 0.0;
     }
     fHardwareBlockSum = 0.0;
     fHardwareBlockSumM2 = 0.0;
     fErrorFlag|=kErrorFlag_sample;
   } else {
     for (Int_t i = 0; i < fBlocksPerEvent; i++) {
       fBlock[i] = fCalibrationFactor * ( (1.0 * fBlock_raw[i] * fBlocksPerEvent / fNumberOfSamples) - fPedestal );
       fBlockM2[i] = 0.0; // second moment is zero for single events
     }
     fHardwareBlockSum = fCalibrationFactor * ( (1.0 * fHardwareBlockSum_raw / fNumberOfSamples) - fPedestal );
     fHardwareBlockSumM2 = 0.0; // second moment is zero for single events
   }
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 1124 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorFlag, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, and VQwDataElement::IsNameEmpty().

Referenced by QwLinearDiodeArray::ProcessEvent().

 {
   if (!IsNameEmpty()){
     for (Int_t i = 0; i < fBlocksPerEvent; i++) {
       this->fBlock[i] = (value1.fBlock[i]) * (value2.fBlock[i]);
       this->fBlock_raw[i] = 0;
       // For a single event the second moment is still zero
       this->fBlockM2[i] = 0.0;
     }
 
     // For a single event the second moment is still zero
     this->fHardwareBlockSumM2 = 0.0;
 
     this->fSoftwareBlockSum_raw = 0;
     this->fHardwareBlockSum_raw = value1.fHardwareBlockSum_raw * value2.fHardwareBlockSum_raw;
     this->fHardwareBlockSum = value1.fHardwareBlockSum * value2.fHardwareBlockSum;
     this->fNumberOfSamples = value1.fNumberOfSamples;
     this->fSequenceNumber  = 0;
     this->fErrorFlag       = (value1.fErrorFlag|value2.fErrorFlag);
   }
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Internally generate random event data.

Implements MQwMockable.

Definition at line 261 of file QwVQWK_Channel.cc.

References fBlocksPerEvent, MQwMockable::fMockAsymmetry, MQwMockable::fMockDriftAmplitude, MQwMockable::fMockDriftFrequency, MQwMockable::fMockDriftPhase, MQwMockable::fMockGaussianMean, MQwMockable::fMockGaussianSigma, MQwMockable::GetRandomValue(), Qw::pi, and SetEventData().

Referenced by QwIntegrationPMT::RandomizeEventData().

 {
   // The blocks are assumed to be independent measurements
   Double_t* block = new Double_t[fBlocksPerEvent];
   Double_t sqrt_fBlocksPerEvent = 0.0;
   sqrt_fBlocksPerEvent = sqrt(fBlocksPerEvent);
 
   // 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]);
   }
 
   // Calculate signal
   for (Int_t i = 0; i < fBlocksPerEvent; i++) {
     block[i] = 
       fMockGaussianMean * (1 + helicity * fMockAsymmetry) / fBlocksPerEvent
       + fMockGaussianSigma / sqrt_fBlocksPerEvent * GetRandomValue()
       + drift / fBlocksPerEvent;
   }
 
   SetEventData(block);
 
   delete block;
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 1048 of file QwVQWK_Channel.cc.

References fBlock_raw, fBlocksPerEvent, VQwDataElement::fErrorFlag, VQwDataElement::fGoodEventCount, fHardwareBlockSum_raw, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, and VQwDataElement::IsNameEmpty().

Referenced by QwLinearDiodeArray::ProcessEvent(), QwQPD::ProcessEvent(), QwQPD::Ratio(), QwBPMCavity::Ratio(), QwLinearDiodeArray::Ratio(), and QwIntegrationPMT::Ratio().

 {
   if (!IsNameEmpty()) {
     *this  = numer;
     *this /= denom;
 
     //  Set the raw values to zero.
     for (Int_t i = 0; i < fBlocksPerEvent; i++) fBlock_raw[i] = 0;
     fHardwareBlockSum_raw = 0;
     fSoftwareBlockSum_raw = 0;
     // Remaining variables
     fNumberOfSamples      = denom.fNumberOfSamples;
     fSequenceNumber       = 0;
     fGoodEventCount       = denom.fGoodEventCount;
     fErrorFlag            = (numer.fErrorFlag|denom.fErrorFlag);
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::Scale ( Double_t Offset )

virtual

Implements VQwHardwareChannel.

Definition at line 1159 of file QwVQWK_Channel.cc.

References fBlock, fBlockM2, fBlocksPerEvent, fHardwareBlockSum, fHardwareBlockSumM2, and VQwDataElement::IsNameEmpty().

Referenced by QwEnergyCalculator::ProcessEvent(), QwBPMCavity::ProcessEvent(), QwQPD::ProcessEvent(), QwLinearDiodeArray::ProcessEvent(), QwEnergyCalculator::Scale(), QwQPD::Scale(), QwBPMCavity::Scale(), QwLinearDiodeArray::Scale(), and QwIntegrationPMT::Scale().

 {
   if (!IsNameEmpty())
     {
       for (Int_t i = 0; i < fBlocksPerEvent; i++) {
         fBlock[i] *= scale;
         fBlockM2[i] *= scale * scale;
       }
       fHardwareBlockSum *= scale;
       fHardwareBlockSumM2 *= scale * scale;
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Implements VQwHardwareChannel.

Definition at line 1589 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwDataElement::fErrorFlag, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumM2, fNumberOfSamples, fSequenceNumber, fSoftwareBlockSum_raw, and VQwDataElement::IsNameEmpty().

 {
   const QwVQWK_Channel* input = dynamic_cast<const QwVQWK_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();
     for(Int_t i = 0; i < fBlocksPerEvent; i++){
       this -> fBlock[i] += scale * input->fBlock[i];
       this->fBlock_raw[i] = 0;
       this -> fBlockM2[i] = 0.0;
     }
     this->fHardwareBlockSum_raw = 0;
     this->fSoftwareBlockSum_raw = 0;
     this -> fHardwareBlockSum += scale * input->fHardwareBlockSum;
     this -> fHardwareBlockSumM2 = 0.0;
     this -> fNumberOfSamples += input->fNumberOfSamples;
     this -> fSequenceNumber  =  0;
     this -> fErrorFlag       |= (input->fErrorFlag);   
   }
   //   QwWarning << "Finsihed with addition"  << QwLog::endl;
   //   PrintValue();
 }

Here is the call graph for this function:

void QwVQWK_Channel::SetCalibrationToVolts ( )

inline

Definition at line 221 of file QwVQWK_Channel.h.

References kVQWK_VoltsPerBit, and VQwHardwareChannel::SetCalibrationFactor().

221 {SetCalibrationFactor(kVQWK_VoltsPerBit);};

VQwHardwareChannel::SetCalibrationFactor

void SetCalibrationFactor(Double_t factor)

Definition: VQwHardwareChannel.h:125

QwVQWK_Channel::kVQWK_VoltsPerBit

static const Double_t kVQWK_VoltsPerBit

Definition: QwVQWK_Channel.h:267

Here is the call graph for this function:

void QwVQWK_Channel::SetDefaultSampleSize ( size_t num_samples_map )

inline

Definition at line 89 of file QwVQWK_Channel.h.

References fNumberOfSamples_map.

Referenced by LoadChannelParameters(), and QwIntegrationPMT::SetDefaultSampleSize().

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

Here is the caller graph for this function:

void QwVQWK_Channel::SetEventData	(	Double_t *	block,
		UInt_t	sequencenumber = `0`
	)

Definition at line 303 of file QwVQWK_Channel.cc.

References fBlock, fBlock_raw, fBlockM2, fBlocksPerEvent, VQwHardwareChannel::fCalibrationFactor, fHardwareBlockSum, fHardwareBlockSum_raw, fHardwareBlockSumM2, fNumberOfSamples, VQwHardwareChannel::fPedestal, fSequenceNumber, and fSoftwareBlockSum_raw.

Referenced by RandomizeEventData(), QwIntegrationPMT::SetEventData(), and SetHardwareSum().

 {
   fHardwareBlockSum = 0.0;
   fHardwareBlockSumM2 = 0.0; // second moment is zero for single events
   for (Int_t i = 0; i < fBlocksPerEvent; i++) {
     fBlock[i] = block[i];
     fBlockM2[i] = 0.0; // second moment is zero for single events
     fHardwareBlockSum += block[i];
   }
 
   fSequenceNumber = sequencenumber;
   fNumberOfSamples = 16680;
 
   Double_t thispedestal = 0.0;
   thispedestal = fPedestal * fNumberOfSamples;
 
   fHardwareBlockSum_raw = 0;
   for (Int_t i = 0; i < fBlocksPerEvent; i++)
     {
       fBlock_raw[i] = Int_t(fBlock[i] / fCalibrationFactor +  thispedestal / (fBlocksPerEvent * 1.));
       fHardwareBlockSum_raw += fBlock_raw[i];
     }
   //  fHardwareBlockSum_raw = fHardwareBlockSum / fCalibrationFactor + thispedestal;
   fSoftwareBlockSum_raw = fHardwareBlockSum_raw;
 
   return;
 }

Here is the caller graph for this function:

void QwVQWK_Channel::SetHardwareSum	(	Double_t	hwsum,
		UInt_t	sequencenumber = `0`
	)

TODO: SetHardwareSum should be removed, and SetEventData should be made protected.

Definition at line 288 of file QwVQWK_Channel.cc.

References fBlocksPerEvent, and SetEventData().

Referenced by QwQPD::SetEventData(), QwBPMCavity::SetEventData(), QwLinearDiodeArray::SetEventData(), and QwIntegrationPMT::SetHardwareSum().

 {
   Double_t* block = new Double_t[fBlocksPerEvent];
   for (Int_t i = 0; i < fBlocksPerEvent; i++)
     block[i] = hwsum / fBlocksPerEvent;
   SetEventData(block);
 
   delete block;
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwVQWK_Channel::SetVQWKSaturationLimt ( Double_t sat_volts = 8.5 )

inline

Definition at line 172 of file QwVQWK_Channel.h.

References fSaturationABSLimit.

Referenced by QwVQWK_Channel().

                                                     {//Set the absolute staturation limit in volts.
     fSaturationABSLimit=sat_volts;
   }

Here is the caller graph for this function:

void QwVQWK_Channel::SubtractValueFrom ( const VQwHardwareChannel * valueptr )

virtual

Implements VQwHardwareChannel.

Definition at line 858 of file QwVQWK_Channel.cc.

References VQwDataElement::GetElementName().

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

Here is the call graph for this function:

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

Definition at line 1036 of file QwVQWK_Channel.cc.

Referenced by QwQPD::ProcessEvent().

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

Here is the caller graph for this function:

Friends And Related Function Documentation

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

friend

Definition at line 1417 of file QwVQWK_Channel.cc.

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

Field Documentation

Bool_t QwVQWK_Channel::bBlock

private

Definition at line 342 of file QwVQWK_Channel.h.

Referenced by ConstructBranchAndVector(), and FillTreeVector().

Bool_t QwVQWK_Channel::bBlock_raw

private

Definition at line 343 of file QwVQWK_Channel.h.

Referenced by ConstructBranchAndVector(), and FillTreeVector().

const Bool_t QwVQWK_Channel::bDEBUG =kFALSE

staticprivate

debugging display purposes

For VQWK data element trimming uses

Definition at line 337 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), FillTreeVector(), and MatchNumberOfSamples().

Bool_t QwVQWK_Channel::bDevice_Error_Code

private

Definition at line 345 of file QwVQWK_Channel.h.

Referenced by ConstructBranchAndVector(), and FillTreeVector().

Bool_t QwVQWK_Channel::bHw_sum

private

Definition at line 340 of file QwVQWK_Channel.h.

Referenced by ConstructBranchAndVector(), and FillTreeVector().

Bool_t QwVQWK_Channel::bHw_sum_raw

private

Definition at line 341 of file QwVQWK_Channel.h.

Referenced by ConstructBranchAndVector(), and FillTreeVector().

Bool_t QwVQWK_Channel::bNum_samples

private

Definition at line 344 of file QwVQWK_Channel.h.

Referenced by ConstructBranchAndVector(), and FillTreeVector().

Bool_t QwVQWK_Channel::bSequence_number

private

Definition at line 346 of file QwVQWK_Channel.h.

Referenced by ConstructBranchAndVector(), and FillTreeVector().

Int_t QwVQWK_Channel::fADC_Same_NumEvt

private

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

Definition at line 327 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), and InitializeChannel().

Double_t QwVQWK_Channel::fBlock[4]

private

Array of the sub-block data.

Definition at line 289 of file QwVQWK_Channel.h.

Referenced by AccumulateRunningSum(), AddChannelOffset(), AssignScaledValue(), Blind(), ClearEventData(), GetValue(), operator*=(), operator+=(), operator-=(), operator/=(), operator=(), PrintInfo(), ProcessEvent(), Product(), Scale(), ScaledAdd(), and SetEventData().

Int_t QwVQWK_Channel::fBlock_raw[4]

private

Array of the sub-block data as read from the module.

Definition at line 283 of file QwVQWK_Channel.h.

Referenced by AssignScaledValue(), ClearEventData(), EncodeEventData(), GetRawValue(), operator*=(), operator+=(), operator-=(), operator=(), PrintInfo(), ProcessEvBuffer(), ProcessEvent(), Product(), Ratio(), ScaledAdd(), and SetEventData().

Double_t QwVQWK_Channel::fBlockError[4]

private

Uncertainty on the sub-block.

Definition at line 298 of file QwVQWK_Channel.h.

Referenced by CalculateRunningAverage(), ClearEventData(), and GetValueError().

Double_t QwVQWK_Channel::fBlockM2[4]

private

Second moment of the sub-block.

Definition at line 297 of file QwVQWK_Channel.h.

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

Short_t QwVQWK_Channel::fBlocksPerEvent

private

Definition at line 277 of file QwVQWK_Channel.h.

Referenced by AddChannelOffset(), AssignScaledValue(), Blind(), CalculateRunningAverage(), ClearEventData(), ConstructHistograms(), FillHistograms(), FillTreeVector(), InitializeChannel(), operator*=(), operator+=(), operator-=(), operator=(), PrintInfo(), ProcessEvBuffer(), ProcessEvent(), Product(), RandomizeEventData(), Ratio(), Scale(), ScaledAdd(), SetEventData(), and SetHardwareSum().

Int_t QwVQWK_Channel::fErrorCount_HWSat

private

check to see ADC channel is saturated

Definition at line 313 of file QwVQWK_Channel.h.

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

Int_t QwVQWK_Channel::fErrorCount_SameHW

private

check to see ADC returning same HW value

Definition at line 317 of file QwVQWK_Channel.h.

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

Int_t QwVQWK_Channel::fErrorCount_sample

private

for sample size check

Definition at line 314 of file QwVQWK_Channel.h.

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

Int_t QwVQWK_Channel::fErrorCount_Sequence

private

sequence number check

Definition at line 316 of file QwVQWK_Channel.h.

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

Int_t QwVQWK_Channel::fErrorCount_SW_HW

private

HW_sum==SW_sum check.

Definition at line 315 of file QwVQWK_Channel.h.

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

Int_t QwVQWK_Channel::fErrorCount_ZeroHW

private

check to see ADC returning zero

Definition at line 318 of file QwVQWK_Channel.h.

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

Double_t QwVQWK_Channel::fHardwareBlockSum

private

Module-based sum of the four sub-blocks.

Definition at line 290 of file QwVQWK_Channel.h.

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

Int_t QwVQWK_Channel::fHardwareBlockSum_raw

private

Module-based sum of the four sub-blocks as read from the module.

Definition at line 284 of file QwVQWK_Channel.h.

Referenced by AssignScaledValue(), ClearEventData(), EncodeEventData(), GetRawValue(), operator*=(), operator+=(), operator-=(), operator=(), PrintInfo(), ProcessEvBuffer(), ProcessEvent(), Product(), Ratio(), ScaledAdd(), and SetEventData().

Double_t QwVQWK_Channel::fHardwareBlockSumError

private

Uncertainty on the hardware sum.

Definition at line 301 of file QwVQWK_Channel.h.

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

Double_t QwVQWK_Channel::fHardwareBlockSumM2

private

Second moment of the hardware sum.

Definition at line 300 of file QwVQWK_Channel.h.

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

size_t QwVQWK_Channel::fNumberOfSamples

private

Number of samples read through the module.

Definition at line 307 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), AssignScaledValue(), ClearEventData(), EncodeEventData(), FillTreeVector(), GetAverageVolts(), GetNumberOfSamples(), InitializeChannel(), MatchNumberOfSamples(), operator*=(), operator+=(), operator-=(), operator=(), PrintInfo(), ProcessEvBuffer(), ProcessEvent(), Product(), Ratio(), ScaledAdd(), and SetEventData().

size_t QwVQWK_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 308 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), InitializeChannel(), and SetDefaultSampleSize().

Int_t QwVQWK_Channel::fNumEvtsWithEventCutsRejected

private

Counts the Event cut rejected events.

Definition at line 320 of file QwVQWK_Channel.h.

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

Double_t QwVQWK_Channel::fPrev_HardwareBlockSum

private

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

Definition at line 330 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), and InitializeChannel().

size_t QwVQWK_Channel::fPreviousSequenceNumber

private

Previous event sequence number for this channel.

Definition at line 306 of file QwVQWK_Channel.h.

Referenced by InitializeChannel().

QwVQWK_Channel* QwVQWK_Channel::fRunningSum

private

Pointer to the running sum for this channel.

Definition at line 263 of file QwVQWK_Channel.h.

Referenced by AccumulateRunningSum(), and InitializeChannel().

Double_t QwVQWK_Channel::fSaturationABSLimit

private

absolute value of the VQWK saturation volt

Definition at line 334 of file QwVQWK_Channel.h.

Referenced by GetVQWKSaturationLimt(), and SetVQWKSaturationLimt().

Int_t QwVQWK_Channel::fSequenceNo_Counter

private

Internal counter to keep track of the sequence number.

Definition at line 329 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), and InitializeChannel().

Int_t QwVQWK_Channel::fSequenceNo_Prev

private

Keep the sequence number of the last event.

Definition at line 328 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), and InitializeChannel().

size_t QwVQWK_Channel::fSequenceNumber

private

Event sequence number for this channel.

Definition at line 305 of file QwVQWK_Channel.h.

Referenced by AssignScaledValue(), ClearEventData(), EncodeEventData(), FillTreeVector(), GetSequenceNumber(), MatchSequenceNumber(), operator*=(), operator+=(), operator-=(), operator=(), PrintInfo(), ProcessEvBuffer(), Product(), Ratio(), ScaledAdd(), and SetEventData().

Int_t QwVQWK_Channel::fSoftwareBlockSum_raw

private

Sum of the data in the four sub-blocks raw.

Definition at line 285 of file QwVQWK_Channel.h.

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

const Bool_t QwVQWK_Channel::kDEBUG = kFALSE

staticprivate

Definition at line 255 of file QwVQWK_Channel.h.

Referenced by ConstructBranch(), and ConstructBranchAndVector().

const Int_t QwVQWK_Channel::kMaxChannels = 8

staticprivate

Definition at line 260 of file QwVQWK_Channel.h.

Referenced by GetBufferOffset().

const Double_t QwVQWK_Channel::kVQWK_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 267 of file QwVQWK_Channel.h.

Referenced by ApplyHWChecks(), GetAverageVolts(), and SetCalibrationToVolts().

const Int_t QwVQWK_Channel::kWordsPerChannel = 6

staticprivate

Definition at line 259 of file QwVQWK_Channel.h.

Referenced by GetBufferOffset(), and ProcessEvBuffer().

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Field Documentation