QwCombinedPMT Class Reference

Combines multiple integration PMTs into weighted sum/average. More...

#include <QwCombinedPMT.h>

Inheritance diagram for QwCombinedPMT:

Collaboration diagram for QwCombinedPMT:

Public Member Functions
	QwCombinedPMT ()
	QwCombinedPMT (TString name)
	QwCombinedPMT (TString subsystemname, TString name)
	QwCombinedPMT (const QwCombinedPMT &source)
	~QwCombinedPMT () override
void	InitializeChannel (TString name, TString datatosave)
	Initialize the combined PMT with a name and data-saving mode.
void	InitializeChannel (TString subsystem, TString name, TString datatosave)
	Initialize the combined PMT with subsystem and name.
void	LinkChannel (TString name)
	Link internal sum channel names to the given detector name.
const QwMollerADC_Channel *	GetChannel (const TString name) const
void	Add (QwIntegrationPMT *pmt, Double_t weight)
	Add a PMT channel to this combination with a weight.
Int_t	ProcessEvBuffer (UInt_t *buffer, UInt_t word_position_in_buffer, UInt_t subelement=0) override
	Process the CODA event buffer for this element.
void	ClearEventData () override
	Clear event-scoped data for the sum channel.
void	PrintErrorCounters ()
void	CalculateSumAndAverage ()
	Compute the weighted sum (and average) from member PMTs.
void	SetRandomEventParameters (Double_t mean, Double_t sigma)
void	SetRandomEventAsymmetry (Double_t asymmetry)
void	RandomizeEventData (int helicity)
void	SetHardwareSum (Double_t hwsum, UInt_t sequencenumber=0)
	Set the hardware-level sum for a sequence (unused for combo).
void	SetEventData (Double_t *block, UInt_t sequencenumber)
	Set the block data for the current event sequence.
void	EncodeEventData (std::vector< UInt_t > &buffer)
void	ProcessEvent ()
	Process event by computing the weighted average of members.
Bool_t	ApplyHWChecks ()
	Apply hardware checks (none needed at combiner level).
Bool_t	ApplySingleEventCuts ()
void	PrintErrorCounters () const override
	Print error counters aggregated by the sum ADC.
void	SetSingleEventCuts (UInt_t errorflag, Double_t LL, Double_t UL, Double_t stability, Double_t burplevel)
	Inherited from VQwDataElement to set the upper and lower limits (fULimit and fLLimit), stability % and the error flag on this channel.
void	SetDefaultSampleSize (Int_t sample_size)
	Set default sample size on the sum ADC.
void	SetEventCutMode (Int_t bcuts)
UInt_t	GetEventcutErrorFlag () override
	return the error flag on this channel/device
void	IncrementErrorCounters ()
Bool_t	CheckForBurpFail (const VQwDataElement *ev_error)
	Check for burp failures by delegating to the sum ADC channel.
UInt_t	UpdateErrorFlag () override
	Update the sum ADC error flag from member PMTs.
void	UpdateErrorFlag (const QwCombinedPMT *ev_error)
	Merge error flags from a reference combined PMT.
void	PrintInfo () const override
	Print multiple lines of information about this data element.
void	PrintValue () const override
	Print single line of value and error of this data element.
QwCombinedPMT &	operator= (const QwCombinedPMT &value)
QwCombinedPMT &	operator+= (const QwCombinedPMT &value)
QwCombinedPMT &	operator-= (const QwCombinedPMT &value)
void	Sum (const QwCombinedPMT &value1, const QwCombinedPMT &value2)
void	Difference (const QwCombinedPMT &value1, const QwCombinedPMT &value2)
void	Ratio (QwCombinedPMT &numer, QwCombinedPMT &denom)
void	Scale (Double_t factor)
void	Normalize (VQwDataElement *denom)
void	AccumulateRunningSum (const QwCombinedPMT &value, Int_t count=0, Int_t ErrorMask=0xFFFFFFF)
void	DeaccumulateRunningSum (QwCombinedPMT &value, Int_t ErrorMask=0xFFFFFFF)
void	CalculateRunningAverage ()
void	SetBlindability (Bool_t isblindable)
void	SetNormalizability (Bool_t isnormalizable)
void	Blind (const QwBlinder *blinder)
	Blind the asymmetry.
void	Blind (const QwBlinder *blinder, const QwCombinedPMT &yield)
	Blind the difference using the yield.
void	SetPedestal (Double_t ped)
void	SetCalibrationFactor (Double_t calib)
void	ConstructHistograms (TDirectory *folder, TString &prefix) override
	Construct the histograms for this data element.
void	FillHistograms () override
	Fill the histograms for this data element.
void	ConstructBranchAndVector (TTree *tree, TString &prefix, QwRootTreeBranchVector &values)
void	ConstructBranch (TTree *tree, TString &prefix)
void	ConstructBranch (TTree *tree, TString &prefix, QwParameterFile &modulelist)
void	FillTreeVector (QwRootTreeBranchVector &values) const
std::vector< QwDBInterface >	GetDBEntry ()
std::vector< QwErrDBInterface >	GetErrDBEntry ()
Public Member Functions inherited from VQwDataElement
	VQwDataElement ()
	Default constructor.
	VQwDataElement (const VQwDataElement &value)
	Copy constructor.
	~VQwDataElement () override
	Virtual destructor.
void	CopyFrom (const VQwDataElement &value)
Bool_t	IsNameEmpty () const
	Is the name of this element empty?
void	SetElementName (const TString &name)
	Set the name of this element.
virtual const TString &	GetElementName () const
	Get the name of this element.
virtual void	LoadChannelParameters (QwParameterFile &)
virtual void	LoadMockDataParameters (QwParameterFile &)
size_t	GetNumberOfDataWords ()
	Get the number of data words in this data element.
UInt_t	GetGoodEventCount () const
virtual void	AssignValueFrom (const VQwDataElement *)
VQwDataElement &	operator+= (const VQwDataElement &)
	Addition-assignment operator.
VQwDataElement &	operator-= (const VQwDataElement &)
	Subtraction-assignment operator.
void	Sum (const VQwDataElement &, const VQwDataElement &)
	Sum operator (base class fallback throws runtime error)
void	Difference (const VQwDataElement &, const VQwDataElement &)
	Difference operator (base class fallback throws runtime error)
void	Ratio (const VQwDataElement &, const VQwDataElement &)
	Ratio operator (base class fallback throws runtime error)
void	SetSingleEventCuts (UInt_t, Double_t, Double_t, Double_t)
	set the upper and lower limits (fULimit and fLLimit), stability % and the error flag on this channel
Bool_t	CheckForBurpFail (const VQwDataElement *)
virtual void	SetNeedsExternalClock (Bool_t)
virtual Bool_t	NeedsExternalClock ()
virtual std::string	GetExternalClockName ()
virtual void	SetExternalClockPtr (const VQwHardwareChannel *)
virtual void	SetExternalClockName (const std::string)
virtual Double_t	GetNormClockValue ()
TString	GetSubsystemName () const
	Return the name of the inheriting subsystem name.
void	SetSubsystemName (TString sysname)
	Set the name of the inheriting subsystem name.
TString	GetModuleType () const
	Return the type of the beam instrument.
void	SetModuleType (TString ModuleType)
	set the type of the beam instrument
Public Member Functions inherited from MQwHistograms
void	ShareHistograms (const MQwHistograms *source)
	Share histogram pointers between objects.

Private Attributes
Int_t	fDataToSave
Double_t	fCalibration
Double_t	fULimit
Double_t	fLLimit
Double_t	fSequenceNo_Prev
Bool_t	fGoodEvent
std::vector< QwIntegrationPMT * >	fElement
	used to validate sequence number in the IsGoodEvent()
std::vector< Double_t >	fWeights
QwIntegrationPMT	fSumADC
Int_t	fDevice_flag
Int_t	fDeviceErrorCode
Bool_t	bEVENTCUTMODE

Static Private Attributes
static const Bool_t	bDEBUG =kFALSE

Additional Inherited Members
Public Types inherited from VQwDataElement
enum	EDataToSave { kRaw = 0 , kDerived , kMoments }
Protected Member Functions inherited from VQwDataElement
void	SetNumberOfDataWords (const UInt_t &numwords)
	Set the number of data words in this data element.
VQwDataElement &	operator= (const VQwDataElement &value)
	Arithmetic assignment operator: Should only copy event-based data.
void	UpdateErrorFlag (const UInt_t &error)
Protected Member Functions inherited from MQwHistograms
	MQwHistograms ()
	Default constructor.
	MQwHistograms (const MQwHistograms &source)
	Copy constructor.
virtual	~MQwHistograms ()
	Virtual destructor.
MQwHistograms &	operator= (const MQwHistograms &value)
void	Fill_Pointer (TH1_ptr hist_ptr, Double_t value)
void	AddHistogram (TH1 *h)
	Register a histogram.
Protected Attributes inherited from VQwDataElement
TString	fElementName
	Name of this data element.
UInt_t	fNumberOfDataWords
	Number of raw data words in this data element.
UInt_t	fGoodEventCount
	Number of good events accumulated in this element.
TString	fSubsystemName
TString	fModuleType
UInt_t	fErrorFlag
	This the standard error code generated for the channel that contains the global/local/stability flags and the Device error code (Unique error code for HW failures)
UInt_t	fErrorConfigFlag
	contains the global/local/stability flags
Protected Attributes inherited from MQwHistograms
std::vector< TH1_ptr >	fHistograms
	Histograms associated with this data element.

Detailed Description

Combines multiple integration PMTs into weighted sum/average.

Definition at line 31 of file QwCombinedPMT.h.

Constructor & Destructor Documentation

QwCombinedPMT() [1/4]

QwCombinedPMT::QwCombinedPMT ( )

inline

Definition at line 34 of file QwCombinedPMT.h.

                 {
    InitializeChannel("", "derived");
  };

References InitializeChannel().

Referenced by AccumulateRunningSum(), Blind(), CheckForBurpFail(), DeaccumulateRunningSum(), Difference(), operator+=(), operator-=(), operator=(), QwCombinedPMT(), Ratio(), Sum(), and UpdateErrorFlag().

Here is the call graph for this function:

Here is the caller graph for this function:

QwCombinedPMT() [2/4]

QwCombinedPMT::QwCombinedPMT ( TString name )

inline

Definition at line 38 of file QwCombinedPMT.h.

                             {
    InitializeChannel(name, "derived");
  };

References InitializeChannel().

Here is the call graph for this function:

QwCombinedPMT() [3/4]

QwCombinedPMT::QwCombinedPMT ( TString subsystemname, TString name )

inline

Definition at line 42 of file QwCombinedPMT.h.

                                                    {
    SetSubsystemName(subsystemname);
    InitializeChannel(subsystemname, name, "derived");
  };

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

Here is the call graph for this function:

QwCombinedPMT() [4/4]

QwCombinedPMT::QwCombinedPMT ( const QwCombinedPMT & source )

inline

Definition at line 46 of file QwCombinedPMT.h.

  : VQwDataElement(source),
    fCalibration(source.fCalibration),
    fElement(source.fElement),
    fWeights(source.fWeights),
    fSumADC(source.fSumADC)
  { }

References fCalibration, fElement, fSumADC, fWeights, QwCombinedPMT(), and VQwDataElement::VQwDataElement().

Here is the call graph for this function:

~QwCombinedPMT()

QwCombinedPMT::~QwCombinedPMT ( )

inlineoverride

Definition at line 53 of file QwCombinedPMT.h.

{ };

Member Function Documentation

AccumulateRunningSum()

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

Definition at line 350 of file QwCombinedPMT.cc.

{
  fSumADC.AccumulateRunningSum(value.fSumADC, count, ErrorMask);
}

References fSumADC, and QwCombinedPMT().

Here is the call graph for this function:

Add()

void QwCombinedPMT::Add	(	QwIntegrationPMT *	pmt,
		Double_t	weight )

Add a PMT channel to this combination with a weight.

Parameters

pmt	Pointer to the PMT to include.
weight	Weight applied to the PMT in the sum/average.

Definition at line 26 of file QwCombinedPMT.cc.

{
  //std::cout<<"QwCombinedPMT: Got "<<pmt->GetElementName()<<"  and weight ="<<weight<<"\n";
  fElement.push_back(pmt);
  fWeights.push_back(weight);
}

References fElement, and fWeights.

ApplyHWChecks()

Bool_t QwCombinedPMT::ApplyHWChecks

(

)

Apply hardware checks (none needed at combiner level).

Definition at line 159 of file QwCombinedPMT.cc.

{
  Bool_t eventokay=kTRUE;
 
 
  return eventokay;
}

ApplySingleEventCuts()

Bool_t QwCombinedPMT::ApplySingleEventCuts

(

)

Definition at line 261 of file QwCombinedPMT.cc.

                                          {
  return fSumADC.ApplySingleEventCuts();
}

References fSumADC.

Blind() [1/2]

void QwCombinedPMT::Blind

(

const QwBlinder *

blinder

)

Blind the asymmetry.

Definition at line 388 of file QwCombinedPMT.cc.

{
  fSumADC.Blind(blinder);
//  fAvgADC.Blind(blinder);
}

References fSumADC.

Blind() [2/2]

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

Blind the difference using the yield.

Definition at line 394 of file QwCombinedPMT.cc.

{
  fSumADC.Blind(blinder, yield.fSumADC);
//  fAvgADC.Blind(blinder, yield.fAvgADC);
}

References fSumADC, and QwCombinedPMT().

Here is the call graph for this function:

CalculateRunningAverage()

void QwCombinedPMT::CalculateRunningAverage

(

)

Definition at line 382 of file QwCombinedPMT.cc.

{
  fSumADC.CalculateRunningAverage();
//  fAvgADC.CalculateRunningAverage();
}

References fSumADC.

CalculateSumAndAverage()

void QwCombinedPMT::CalculateSumAndAverage

(

)

Compute the weighted sum (and average) from member PMTs.

Definition at line 109 of file QwCombinedPMT.cc.

{
 
  Bool_t ldebug = kFALSE;
  Double_t  total_weights=0.0;
 
  fSumADC.ClearEventData();
  static QwIntegrationPMT tmpADC("tmpADC");
 
  for (size_t i=0;i<fElement.size();i++)
    {
      //std::cout<<"=========fElement["<<i<<"]=========="<<std::endl;
      //fElement[i]->Print();
      tmpADC = *(fElement[i]);
      //std::cout<<"=========tmpADC========="<<std::endl;
      //tmpADC->Print();
      tmpADC.Scale(fWeights[i]);
      fSumADC += tmpADC;
      total_weights += fWeights[i];
    }
 
//   fAvgADC = fSumADC;
//   if (total_weights!=0.0)
//     fAvgADC.Scale(1/total_weights);
 
   if (total_weights!=0.0)
     fSumADC.Scale(1/total_weights);
 
 
  if (ldebug)
    {
      std::cout<<"QwCombinedPMT::CalculateAverage()"<<std::endl;
//      fAvgADC.PrintInfo();
      fSumADC.PrintInfo();
 
//      std::cout<<"QwCombinedPMT: "<<GetElementName()
//      <<"\nweighted average of hardware sums = "<<fAvgADC.GetValue()<<"\n";
      std::cout<<"QwCombinedPMT: "<<GetElementName()
      <<"\nweighted average of hardware sums = "<<fSumADC.GetValue()<<"\n";
 
      for (size_t i=0;i<4;i++)
        {
          std::cout<<"weighted average of block["<<i<<"] = "<<fSumADC.GetValue(i+1)<<"\n";
        }
    }
 
}

References fElement, fSumADC, fWeights, VQwDataElement::GetElementName(), and QwIntegrationPMT::Scale().

Referenced by ProcessEvent().

Here is the call graph for this function:

Here is the caller graph for this function:

CheckForBurpFail()

Bool_t QwCombinedPMT::CheckForBurpFail

(

const VQwDataElement *

ev_error

)

Check for burp failures by delegating to the sum ADC channel.

Parameters

ev_error

Reference combined PMT to compare against.

Returns

kTRUE if a burp failure was detected; otherwise kFALSE.

Definition at line 207 of file QwCombinedPMT.cc.

                                                                    {
  Bool_t burpstatus = kFALSE;
  try {
    if(typeid(*ev_error)==typeid(*this)) {
      //std::cout<<" Here in QwCombinedPMT::CheckForBurpFail \n";
      if (this->GetElementName()!="") {
        const QwCombinedPMT* value_pmt = dynamic_cast<const QwCombinedPMT* >(ev_error);
        burpstatus |= fSumADC.CheckForBurpFail(&(value_pmt->fSumADC)); 
      }
    } else {
      TString loc="Standard exception from QwCombinedPMT::CheckForBurpFail :"+
        ev_error->GetElementName()+" "+this->GetElementName()+" are not of the "
        +"same type";
      throw std::invalid_argument(loc.Data());
    }
  } catch (std::exception& e) {
    std::cerr<< e.what()<<std::endl;
  }
  return burpstatus;
};

References fSumADC, VQwDataElement::GetElementName(), QwCombinedPMT(), and VQwDataElement::VQwDataElement().

Here is the call graph for this function:

ClearEventData()

void QwCombinedPMT::ClearEventData ( )

overridevirtual

Clear event-scoped data for the sum channel.

Reimplemented from VQwDataElement.

Definition at line 88 of file QwCombinedPMT.cc.

{
  fSumADC.ClearEventData();
//  fAvgADC.ClearEventData();
}

References fSumADC.

ConstructBranch() [1/2]

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

Definition at line 461 of file QwCombinedPMT.cc.

{
  if (GetElementName()=="")
    {
      //  This channel is not used, so skip
    }
  else
    {
      TString sumprefix =  prefix+"";
      fSumADC.ConstructBranch(tree, sumprefix);
    }
  return;
}

References fSumADC, and VQwDataElement::GetElementName().

Here is the call graph for this function:

ConstructBranch() [2/2]

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

Definition at line 475 of file QwCombinedPMT.cc.

{
  TString devicename;
  devicename=GetElementName();
  devicename.ToLower();
  if (GetElementName()=="")
    {
      //  This channel is not used, so skip
    }
  else
    {
      if (modulelist.HasValue(devicename)){
    TString sumprefix =  prefix+"";
    fSumADC.ConstructBranch(tree, sumprefix);
    QwMessage <<" Tree leave added to "<<devicename<<QwLog::endl;
      }
    }
  return;
}

References QwLog::endl(), fSumADC, VQwDataElement::GetElementName(), QwParameterFile::HasValue(), and QwMessage.

Here is the call graph for this function:

ConstructBranchAndVector()

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

Definition at line 446 of file QwCombinedPMT.cc.

{
  if (GetElementName()=="")
    {
      //  This channel is not used, so skip
    }
  else
    {
      TString sumprefix =  prefix+"";
      fSumADC.ConstructBranchAndVector(tree, sumprefix,values);
 
    }
  return;
}

References fSumADC, and VQwDataElement::GetElementName().

Here is the call graph for this function:

ConstructHistograms()

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

overridevirtual

Construct the histograms for this data element.

Implements VQwDataElement.

Definition at line 413 of file QwCombinedPMT.cc.

{
  if (GetElementName()=="")
    {
      //  This channel is not used, so skip filling the histograms.
    }
  else
    {
      TString sumprefix = prefix+TString("");
      fSumADC.ConstructHistograms(folder, sumprefix);
//      TString avgprefix = prefix+TString("");
//      fAvgADC.ConstructHistograms(folder, avgprefix);
    }
  return;
 
}

References fSumADC, and VQwDataElement::GetElementName().

Here is the call graph for this function:

DeaccumulateRunningSum()

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

Definition at line 355 of file QwCombinedPMT.cc.

{
  fSumADC.DeaccumulateRunningSum(value.fSumADC, ErrorMask);
}

References fSumADC, and QwCombinedPMT().

Here is the call graph for this function:

Difference()

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

Definition at line 345 of file QwCombinedPMT.cc.

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

References QwCombinedPMT().

Here is the call graph for this function:

EncodeEventData()

void QwCombinedPMT::EncodeEventData

(

std::vector< UInt_t > &

buffer

)

FillHistograms()

void QwCombinedPMT::FillHistograms ( )

overridevirtual

Fill the histograms for this data element.

Implements VQwDataElement.

Definition at line 430 of file QwCombinedPMT.cc.

{
  if (GetElementName()=="")
    {
      //  This channel is not used, so skip filling the histograms.
    }
  else
    {
      fSumADC.FillHistograms();
//      fAvgADC.FillHistograms();
    }
 
 
  return;
}

References fSumADC, and VQwDataElement::GetElementName().

Here is the call graph for this function:

FillTreeVector()

void QwCombinedPMT::FillTreeVector

(

QwRootTreeBranchVector &

values

)

const

Definition at line 496 of file QwCombinedPMT.cc.

{
  if (GetElementName()=="") {
    //  This channel is not used, so skip filling the histograms.
  } else {
    fSumADC.FillTreeVector(values);
//  fAvgADC.FillTreeVector(values);
  }
}

References fSumADC, and VQwDataElement::GetElementName().

Here is the call graph for this function:

GetChannel()

const QwMollerADC_Channel * QwCombinedPMT::GetChannel ( const TString name ) const

inline

Definition at line 60 of file QwCombinedPMT.h.

                                                                  {
    if (fSumADC.GetElementName() == name) return fSumADC.GetChannel(name);
    else return 0;
  };

References fSumADC.

Referenced by VQwDetectorArray::PublishInternalValues().

Here is the caller graph for this function:

GetDBEntry()

std::vector< QwDBInterface > QwCombinedPMT::GetDBEntry

(

)

GetErrDBEntry()

std::vector< QwErrDBInterface > QwCombinedPMT::GetErrDBEntry

(

)

GetEventcutErrorFlag()

UInt_t QwCombinedPMT::GetEventcutErrorFlag ( )

inlineoverridevirtual

return the error flag on this channel/device

Reimplemented from VQwDataElement.

Definition at line 93 of file QwCombinedPMT.h.

                                        {//return the error flag
    return fSumADC.GetEventcutErrorFlag();
  }

References fSumADC.

Referenced by UpdateErrorFlag().

Here is the caller graph for this function:

IncrementErrorCounters()

void QwCombinedPMT::IncrementErrorCounters ( )

inline

Definition at line 97 of file QwCombinedPMT.h.

                               {
    fSumADC.IncrementErrorCounters();
  }

References fSumADC.

InitializeChannel() [1/2]

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

Initialize the combined PMT with a name and data-saving mode.

Parameters

name	Detector name used for branches and histograms.
datatosave	Storage mode (e.g., "raw" or "derived").

Definition at line 39 of file QwCombinedPMT.cc.

{
  SetElementName(name);
  //SetPedestal(0.);
  //SetCalibrationFactor(1.);
  if (datatosave=="raw") fDataToSave=kRaw;
  else
    if (datatosave=="derived") fDataToSave=kDerived;
 
  fSumADC.InitializeChannel(name, datatosave);
  SetBlindability(kTRUE);
  return;
}

References fDataToSave, fSumADC, VQwDataElement::kDerived, VQwDataElement::kRaw, SetBlindability(), and VQwDataElement::SetElementName().

Referenced by QwCombinedPMT(), QwCombinedPMT(), and QwCombinedPMT().

Here is the call graph for this function:

Here is the caller graph for this function:

InitializeChannel() [2/2]

void QwCombinedPMT::InitializeChannel	(	TString	subsystemname,
		TString	name,
		TString	datatosave )

Initialize the combined PMT with subsystem and name.

Parameters

subsystemname	Subsystem identifier.
name	Detector name used for branches and histograms.
datatosave	Storage mode (e.g., "raw" or "derived").

Definition at line 59 of file QwCombinedPMT.cc.

{
  SetElementName(name);
  //SetPedestal(0.);
  //SetCalibrationFactor(1.);
  if (datatosave=="raw") fDataToSave=kRaw;
  else
    if (datatosave=="derived") fDataToSave=kDerived;
 
  fSumADC.InitializeChannel(subsystemname, "QwCombinedPMT", name, datatosave);
  SetBlindability(kTRUE);
 
  return;
}

References fDataToSave, fSumADC, VQwDataElement::kDerived, VQwDataElement::kRaw, SetBlindability(), and VQwDataElement::SetElementName().

Here is the call graph for this function:

LinkChannel()

void QwCombinedPMT::LinkChannel

(

TString

name

)

Link internal sum channel names to the given detector name.

Definition at line 75 of file QwCombinedPMT.cc.

{
  Bool_t local_debug = false;
  SetElementName(name);
//  TString sumstr = name+TString("_sum");
  TString sumstr = name+TString("");
  fSumADC.SetElementName(sumstr);
//   TString avgstr = name+TString("_avg");
//   fAvgADC.SetElementName(avgstr);
  if(local_debug) std::cout<<"linked combined PMT channel "<< GetElementName()<<std::endl;
}

References fSumADC, VQwDataElement::GetElementName(), and VQwDataElement::SetElementName().

Here is the call graph for this function:

Normalize()

void QwCombinedPMT::Normalize

(

VQwDataElement *

denom

)

Definition at line 374 of file QwCombinedPMT.cc.

{
  fSumADC.Normalize(denom);
 
  return;
}

References fSumADC, and VQwDataElement::VQwDataElement().

Here is the call graph for this function:

operator+=()

QwCombinedPMT & QwCombinedPMT::operator+=

(

const QwCombinedPMT &

value

)

Definition at line 302 of file QwCombinedPMT.cc.

{
  //std::cout<<"Calling QwCombinedPMT::operator+="<<std::endl;
  if ( (GetElementName()!="") && (this->fElement.size()==value.fElement.size()) )
    {
      for (size_t i=0; i<value.fElement.size(); i++)
        {
          *(this->fElement[i]) += *(value.fElement[i]);
          this->fWeights[i] += value.fWeights[i];
        }
 
      this->fSumADC+=value.fSumADC;
//      this->fAvgADC+=value.fAvgADC;
    }
 
  return *this;
}

References fElement, fSumADC, fWeights, VQwDataElement::GetElementName(), and QwCombinedPMT().

Here is the call graph for this function:

operator-=()

QwCombinedPMT & QwCombinedPMT::operator-=

(

const QwCombinedPMT &

value

)

Definition at line 320 of file QwCombinedPMT.cc.

{
  //std::cout<<"Calling QwCombinedPMT::operator-="<<std::endl;
  if ( (GetElementName()!="") && (this->fElement.size()==value.fElement.size()) )
    {
      for (size_t i=0; i<value.fElement.size(); i++)
        {
          *(this->fElement[i]) -= *(value.fElement[i]);
          this->fWeights[i] -= value.fWeights[i];
        }
 
      this->fSumADC-=value.fSumADC;
//      this->fAvgADC-=value.fAvgADC;
      //value.fSumADC.Print();
      //fSumADC.Print();
    }
 
  return *this;
}

References fElement, fSumADC, fWeights, VQwDataElement::GetElementName(), and QwCombinedPMT().

Here is the call graph for this function:

operator=()

QwCombinedPMT & QwCombinedPMT::operator=

(

const QwCombinedPMT &

value

)

Definition at line 273 of file QwCombinedPMT.cc.

{
  //std::cout<<"Calling QwCombinedPMT::operator="<<std::endl;
  if ( (GetElementName()!="") && (this->fElement.size()==value.fElement.size()) )
    {
      for (size_t i=0; i<value.fElement.size(); i++)
        {
          //std::cout<<"value.fWeight["<<i<<"]="<<value.fWeights[i]<<std::endl;
          //std::cout<<"=========value.fElement["<<i<<"]=========="<<std::endl;
          //value.fElement[i]->Print();
          (this->fElement[i]) = (value.fElement[i]);
          this->fWeights[i] = value.fWeights[i];
          //std::cout<<"this->fWeight["<<i<<"]="<<fWeights[i]<<std::endl;
          //std::cout<<"=========this->fElement["<<i<<"]=========="<<std::endl;
          //fElement[i]->Print();
        }
 
      this->fSumADC=value.fSumADC;
//      this->fAvgADC=value.fAvgADC;
 
      //std::cout<<"value.fSumADC"<<std::endl;
      //value.fSumADC.Print();
      //std::cout<<"this->fSumADC"<<std::endl;
      //fSumADC.Print();
    }
 
  return *this;
}

References fElement, fSumADC, fWeights, VQwDataElement::GetElementName(), and QwCombinedPMT().

Here is the call graph for this function:

PrintErrorCounters() [1/2]

void QwCombinedPMT::PrintErrorCounters

(

)

PrintErrorCounters() [2/2]

void QwCombinedPMT::PrintErrorCounters ( ) const

overridevirtual

Print error counters aggregated by the sum ADC.

Reimplemented from VQwDataElement.

Definition at line 197 of file QwCombinedPMT.cc.

{
  fSumADC.PrintErrorCounters();
}

References fSumADC.

PrintInfo()

void QwCombinedPMT::PrintInfo ( ) const

overridevirtual

Print multiple lines of information about this data element.

Reimplemented from VQwDataElement.

Definition at line 406 of file QwCombinedPMT.cc.

{
  fSumADC.PrintInfo();
//  fAvgADC.PrintInfo();
}

References fSumADC.

PrintValue()

void QwCombinedPMT::PrintValue ( ) const

overridevirtual

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

Reimplemented from VQwDataElement.

Definition at line 400 of file QwCombinedPMT.cc.

{
  fSumADC.PrintValue();
//  fAvgADC.PrintValue();
}

References fSumADC.

ProcessEvBuffer()

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

overridevirtual

Process the CODA event buffer for this element.

Implements VQwDataElement.

Definition at line 267 of file QwCombinedPMT.cc.

{
  return 0;
}

References QwBPMStripline< T >::subelement.

ProcessEvent()

void QwCombinedPMT::ProcessEvent

(

)

Process event by computing the weighted average of members.

Definition at line 179 of file QwCombinedPMT.cc.

{
//Calculate the weigted averages of the hardware sum and each of the four blocks.
  CalculateSumAndAverage();
  //fSumADC.ProcessEvent(); //This is not necessary for combined devices-Rakitha 11-15-2010
 
  return;
}

References CalculateSumAndAverage().

Here is the call graph for this function:

RandomizeEventData()

void QwCombinedPMT::RandomizeEventData

(

int

helicity

)

Ratio()

void QwCombinedPMT::Ratio	(	QwCombinedPMT &	numer,
		QwCombinedPMT &	denom )

Definition at line 360 of file QwCombinedPMT.cc.

{
  //std::cout<<"Calling QwCombinedPMT::Ratio"<<std::endl;
  fSumADC.Ratio(numer.fSumADC,denom.fSumADC);
//  fAvgADC.Ratio(numer.fAvgADC,denom.fAvgADC);
  return;
}

References fSumADC, and QwCombinedPMT().

Here is the call graph for this function:

Scale()

void QwCombinedPMT::Scale

(

Double_t

factor

)

Definition at line 368 of file QwCombinedPMT.cc.

{
  fSumADC.Scale(factor);
//  fAvgADC.Scale(factor);
  return;
}

References fSumADC.

SetBlindability()

void QwCombinedPMT::SetBlindability ( Bool_t isblindable )

inline

Definition at line 121 of file QwCombinedPMT.h.

{fSumADC.SetBlindability(isblindable);};

References fSumADC.

Referenced by InitializeChannel(), InitializeChannel(), and VQwDetectorArray::LoadChannelMap().

Here is the caller graph for this function:

SetCalibrationFactor()

void QwCombinedPMT::SetCalibrationFactor

(

Double_t

calib

)

SetDefaultSampleSize()

void QwCombinedPMT::SetDefaultSampleSize

(

Int_t

sample_size

)

Set default sample size on the sum ADC.

Definition at line 190 of file QwCombinedPMT.cc.

{
  fSumADC.SetDefaultSampleSize((size_t)sample_size);
}

References fSumADC.

SetEventCutMode()

void QwCombinedPMT::SetEventCutMode ( Int_t bcuts )

inline

Definition at line 89 of file QwCombinedPMT.h.

                                   {
    bEVENTCUTMODE=bcuts;
    fSumADC.SetEventCutMode(bcuts);
  }

References bEVENTCUTMODE, and fSumADC.

SetEventData()

void QwCombinedPMT::SetEventData	(	Double_t *	block,
		UInt_t	sequencenumber )

Set the block data for the current event sequence.

Definition at line 102 of file QwCombinedPMT.cc.

{
  fSumADC.SetEventData(block, sequencenumber);
//  fAvgADC.SetEventData(block, sequencenumber);
}

References fSumADC.

SetHardwareSum()

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

Set the hardware-level sum for a sequence (unused for combo).

Definition at line 96 of file QwCombinedPMT.cc.

{
 
}

SetNormalizability()

void QwCombinedPMT::SetNormalizability ( Bool_t isnormalizable )

inline

Definition at line 122 of file QwCombinedPMT.h.

{fSumADC.SetNormalizability(isnormalizable);}; 

References fSumADC.

Referenced by VQwDetectorArray::LoadChannelMap().

Here is the caller graph for this function:

SetPedestal()

void QwCombinedPMT::SetPedestal

(

Double_t

ped

)

SetRandomEventAsymmetry()

void QwCombinedPMT::SetRandomEventAsymmetry

(

Double_t

asymmetry

)

SetRandomEventParameters()

void QwCombinedPMT::SetRandomEventParameters	(	Double_t	mean,
		Double_t	sigma )

SetSingleEventCuts()

void QwCombinedPMT::SetSingleEventCuts	(	UInt_t	errorflag,
		Double_t	LL = 0,
		Double_t	UL = 0,
		Double_t	stability = 0,
		Double_t	burplevel = 0 )

Inherited from VQwDataElement to set the upper and lower limits (fULimit and fLLimit), stability % and the error flag on this channel.

Configure detailed single-event cuts forwarded to the sum ADC.

Definition at line 169 of file QwCombinedPMT.cc.

                                                                                                                                {
  //set the unique tag to identify device type (Int.PMT & Comb. PMT)
  //errorflag|=kPMTErrorFlag;
  QwMessage<<"QwCombinedPMT Error Code passing to QwIntegrationPMT "<<errorflag<<QwLog::endl;
  fSumADC.SetSingleEventCuts(errorflag,LL,UL,stability,burplevel);
  
}

References QwLog::endl(), fSumADC, and QwMessage.

Here is the call graph for this function:

Sum()

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

Definition at line 340 of file QwCombinedPMT.cc.

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

References QwCombinedPMT().

Here is the call graph for this function:

UpdateErrorFlag() [1/2]

UInt_t QwCombinedPMT::UpdateErrorFlag ( )

overridevirtual

Update the sum ADC error flag from member PMTs.

Reimplemented from VQwDataElement.

Definition at line 231 of file QwCombinedPMT.cc.

{
  for (size_t i=0;i<fElement.size();i++) {
    fSumADC.UpdateErrorFlag(fElement[i]);
  }
  return GetEventcutErrorFlag();
}

References fElement, fSumADC, and GetEventcutErrorFlag().

Here is the call graph for this function:

UpdateErrorFlag() [2/2]

void QwCombinedPMT::UpdateErrorFlag

(

const QwCombinedPMT *

ev_error

)

Merge error flags from a reference combined PMT.

Definition at line 241 of file QwCombinedPMT.cc.

                                                                {
  try {
    if(typeid(*ev_error)==typeid(*this)) {
      // std::cout<<" Here in QwCombinedPMT::UpdateErrorFlag \n";
      if (this->GetElementName()!="") {
    fSumADC.UpdateErrorFlag(&(ev_error->fSumADC));
      }
    } else {
      TString loc="Standard exception from QwCombinedPMT::UpdateErrorFlag :"+
        ev_error->GetElementName()+" "+this->GetElementName()+" are not of the "
        +"same type";
      throw std::invalid_argument(loc.Data());
    }
  } catch (std::exception& e) {
    std::cerr<< e.what()<<std::endl;
  }  
};

References fSumADC, VQwDataElement::GetElementName(), and QwCombinedPMT().

Here is the call graph for this function:

Field Documentation

bDEBUG

const Bool_t QwCombinedPMT::bDEBUG =kFALSE

staticprivate

If this set to kFALSE then Event cuts do not depend on HW checks. This is set externally through the qweak_beamline_eventcuts.map

Definition at line 177 of file QwCombinedPMT.h.

bEVENTCUTMODE

Bool_t QwCombinedPMT::bEVENTCUTMODE

private

keep the device HW status using a unique code from the QwMollerADC_Channel::fDeviceErrorCode

Definition at line 173 of file QwCombinedPMT.h.

Referenced by SetEventCutMode().

fCalibration

Double_t QwCombinedPMT::fCalibration

private

Definition at line 154 of file QwCombinedPMT.h.

Referenced by QwCombinedPMT().

fDataToSave

Int_t QwCombinedPMT::fDataToSave

private

Definition at line 153 of file QwCombinedPMT.h.

Referenced by InitializeChannel(), and InitializeChannel().

fDevice_flag

Int_t QwCombinedPMT::fDevice_flag

private

Definition at line 166 of file QwCombinedPMT.h.

fDeviceErrorCode

Int_t QwCombinedPMT::fDeviceErrorCode

private

sets the event cut level for the device fDevice_flag=1 Event cuts & HW check, fDevice_flag=0 HW check, fDevice_flag=-1 no check

Definition at line 170 of file QwCombinedPMT.h.

fElement

std::vector<QwIntegrationPMT*> QwCombinedPMT::fElement

private

used to validate sequence number in the IsGoodEvent()

Definition at line 160 of file QwCombinedPMT.h.

Referenced by Add(), CalculateSumAndAverage(), operator+=(), operator-=(), operator=(), QwCombinedPMT(), and UpdateErrorFlag().

fGoodEvent

Bool_t QwCombinedPMT::fGoodEvent

private

Definition at line 158 of file QwCombinedPMT.h.

fLLimit

Double_t QwCombinedPMT::fLLimit

private

Definition at line 155 of file QwCombinedPMT.h.

fSequenceNo_Prev

Double_t QwCombinedPMT::fSequenceNo_Prev

private

Definition at line 156 of file QwCombinedPMT.h.

fSumADC

QwIntegrationPMT QwCombinedPMT::fSumADC

private

Definition at line 163 of file QwCombinedPMT.h.

Referenced by AccumulateRunningSum(), ApplySingleEventCuts(), Blind(), Blind(), CalculateRunningAverage(), CalculateSumAndAverage(), CheckForBurpFail(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), DeaccumulateRunningSum(), FillHistograms(), FillTreeVector(), GetChannel(), GetEventcutErrorFlag(), IncrementErrorCounters(), InitializeChannel(), InitializeChannel(), LinkChannel(), Normalize(), operator+=(), operator-=(), operator=(), PrintErrorCounters(), PrintInfo(), PrintValue(), QwCombinedPMT(), Ratio(), Scale(), SetBlindability(), SetDefaultSampleSize(), SetEventCutMode(), SetEventData(), SetNormalizability(), SetSingleEventCuts(), UpdateErrorFlag(), and UpdateErrorFlag().

fULimit

Double_t QwCombinedPMT::fULimit

private

Definition at line 155 of file QwCombinedPMT.h.

fWeights

std::vector<Double_t> QwCombinedPMT::fWeights

private

Definition at line 161 of file QwCombinedPMT.h.

Referenced by Add(), CalculateSumAndAverage(), operator+=(), operator-=(), operator=(), and QwCombinedPMT().

---

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

• QwCombinedPMT.h
• QwCombinedPMT.cc