Quadrant photodiode BPM computing X/Y and effective charge. More...

#include <QwQPD.h>

Inheritance diagram for QwQPD:

Collaboration diagram for QwQPD:

Public Member Functions
	QwQPD ()

	QwQPD (TString name)

	QwQPD (TString subsystemname, TString name)

	QwQPD (const QwQPD &source)

	~QwQPD () override

void	InitializeChannel (TString name)

void	InitializeChannel (TString subsystem, TString name)

void	LoadChannelParameters (QwParameterFile &paramfile) override

void	GetCalibrationFactors (Double_t AlphaX, Double_t AlphaY)

void	ClearEventData () override

Int_t	ProcessEvBuffer (UInt_t *buffer, UInt_t word_position_in_buffer, UInt_t indexnumber) override

void	ProcessEvent () override

const VQwHardwareChannel *	GetPosition (EBeamPositionMonitorAxis axis) const override

const VQwHardwareChannel *	GetEffectiveCharge () const override

TString	GetSubElementName (Int_t subindex) override

void	GetAbsolutePosition () override

Bool_t	ApplyHWChecks ()

Bool_t	ApplySingleEventCuts () override

void	SetSingleEventCuts (TString ch_name, UInt_t errorflag, Double_t minX, Double_t maxX, 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	SetEventCutMode (Int_t bcuts) override

void	IncrementErrorCounters () override

void	PrintErrorCounters () const override

UInt_t	GetEventcutErrorFlag () override

UInt_t	UpdateErrorFlag () override

void	UpdateErrorFlag (const VQwBPM *ev_error) override

Bool_t	CheckForBurpFail (const VQwDataElement *ev_error) override

void	SetDefaultSampleSize (Int_t sample_size) override

void	SetRandomEventParameters (Double_t meanX, Double_t sigmaX, Double_t meanY, Double_t sigmaY) override

void	RandomizeEventData (int helicity=0, double time=0.0) override

void	SetEventData (Double_t *block, UInt_t sequencenumber)

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

void	SetSubElementPedestal (Int_t j, Double_t value) override

void	SetSubElementCalibrationFactor (Int_t j, Double_t value) override

void	Ratio (QwQPD &numer, QwQPD &denom)

void	Scale (Double_t factor) override

VQwBPM &	operator= (const VQwBPM &value) override

VQwBPM &	operator+= (const VQwBPM &value) override

VQwBPM &	operator-= (const VQwBPM &value) override

virtual QwQPD &	operator= (const QwQPD &value)

virtual QwQPD &	operator+= (const QwQPD &value)

virtual QwQPD &	operator-= (const QwQPD &value)

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

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

void	DeaccumulateRunningSum (VQwBPM &value, Int_t ErrorMask=0xFFFFFFF) override

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

void	CalculateRunningAverage () override

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

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

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

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

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

void	FillTreeVector (QwRootTreeBranchVector &values) const override

void	MakeQPDList ()

Public Member Functions inherited from VQwBPM
	VQwBPM ()

	VQwBPM (TString &)

	VQwBPM (const VQwBPM &source)

	~VQwBPM () override

void	InitializeChannel (TString name)
	Initialize common BPM state and set the element name.

virtual void	GetProjectedPosition (VQwBPM *)

virtual size_t	GetNumberOfElements ()

virtual void	FillRawEventData ()

void	GetSurveyOffsets (Double_t Xoffset, Double_t Yoffset, Double_t Zoffset)
	Store geometry/survey offsets for absolute position calibration.

void	GetElectronicFactors (Double_t BSENfactor, Double_t AlphaX, Double_t AlphaY)
	Apply per-detector electronic calibration and relative gains.

void	SetRotation (Double_t)
	Set detector rotation angle and update cached trigonometric values.

void	SetRotationOff ()

void	SetSingleEventCuts (TString, Double_t, Double_t)

void	SetSingleEventCuts (TString, UInt_t, Double_t, Double_t, Double_t, Double_t)

void	SetGains (TString pos, Double_t value)

void	SetRootSaveStatus (TString &prefix)

Double_t	GetPositionInZ () const

virtual void	Ratio (VQwBPM &, VQwBPM &)

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

virtual const VQwHardwareChannel *	GetAngleX () const

virtual const VQwHardwareChannel *	GetAngleY () const

virtual void	SetBPMForCombo (const VQwBPM *, Double_t, Double_t, Double_t, Double_t)

virtual void	SetResolution (Double_t resolutionX, Double_t resolutionY)

virtual void	SetRandomEventParameters (Double_t, Double_t, Double_t, Double_t, Double_t, Double_t, Double_t, Double_t)

virtual void	SetRandomEventAsymmetry (Double_t)

virtual void	ApplyResolutionSmearing ()

virtual void	ApplyResolutionSmearing (EBeamPositionMonitorAxis)

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

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

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

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.

Static Public Member Functions
static UInt_t	GetSubElementIndex (TString subname)

Static Public Member Functions inherited from VQwBPM
static VQwBPM *	CreateStripline (TString subsystemname, TString type, TString name)
	A fast way of creating a BPM stripline of specified type.

static VQwBPM *	CreateStripline (const VQwBPM &source)

static VQwBPM *	CreateCombo (TString subsystemname, TString type, TString name)
	A fast way of creating a BPM stripline of specified type.

static VQwBPM *	CreateCombo (const VQwBPM &source)

Protected Member Functions
VQwHardwareChannel *	GetSubelementByName (TString ch_name) override

Protected Member Functions inherited from VQwBPM
virtual VQwHardwareChannel *	GetPosition (EBeamPositionMonitorAxis axis)

VQwHardwareChannel *	GetSubelementByIndex (size_t index)

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
std::array< QwVQWK_Channel, 4 >	fPhotodiode

std::array< QwVQWK_Channel, 2 >	fRelPos

std::array< QwVQWK_Channel, 2 >	fAbsPos

QwVQWK_Channel	fEffectiveCharge

std::vector< QwVQWK_Channel >	fQPDElementList

Protected Attributes inherited from VQwBPM
std::vector< TString >	fSubelementNames

Double_t	fPositionCenter [3]

Double_t	fQwStriplineCalibration

Double_t	fQwStriplineCorrection

Double_t	fRelativeGains [2]

Double_t	fGains [2]

Bool_t	bRotated

Double_t	fRotationAngle

Double_t	fCosRotation

Double_t	fSinRotation

Double_t	fResolution [2]

Bool_t	fGoodEvent

Bool_t	bFullSave

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.

Private Attributes
Double_t	fQwQPDCalibration [2]

Static Private Attributes
static const TString	subelement [4] ={"BR","TR","BL","TL"}

Additional Inherited Members
Public Types inherited from VQwBPM
enum	EBeamPositionMonitorAxis { kXAxis =0 , kYAxis , kNumAxes }

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

Static Protected Attributes inherited from VQwBPM
static const TString	kAxisLabel [2] ={"X","Y"}

static const TString	axis [3]

static const Bool_t	bDEBUG =kFALSE

Detailed Description

Quadrant photodiode BPM computing X/Y and effective charge.

Definition at line 38 of file QwQPD.h.

Constructor & Destructor Documentation

◆ QwQPD() [1/4]

QwQPD::QwQPD ( )

inline

Definition at line 43 of file QwQPD.h.

43 {

44 };

Referenced by AccumulateRunningSum(), AccumulateRunningSum(), CheckForBurpFail(), DeaccumulateRunningSum(), DeaccumulateRunningSum(), operator+=(), operator+=(), operator-=(), operator-=(), operator=(), operator=(), QwQPD(), Ratio(), and UpdateErrorFlag().

Here is the caller graph for this function:

◆ QwQPD() [2/4]

QwQPD::QwQPD ( TString name )

inline

Definition at line 45 of file QwQPD.h.

                     :VQwBPM(name){
    InitializeChannel(name);
  };

References InitializeChannel(), and VQwBPM::VQwBPM().

Here is the call graph for this function:

◆ QwQPD() [3/4]

QwQPD::QwQPD	(	TString	subsystemname,
		TString	name )

inline

Definition at line 48 of file QwQPD.h.

                                            :VQwBPM(name){
    SetSubsystemName(subsystemname);
    InitializeChannel(subsystemname, name);
    fQwQPDCalibration[0] = 1.0;
    fQwQPDCalibration[1] = 1.0;
  };

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

Here is the call graph for this function:

◆ QwQPD() [4/4]

QwQPD::QwQPD ( const QwQPD & source )

inline

Definition at line 54 of file QwQPD.h.

  : VQwBPM(source),
    fEffectiveCharge(source.fEffectiveCharge)
  {
    QwCopyArray(source.fPhotodiode, fPhotodiode);
    QwCopyArray(source.fRelPos, fRelPos);
    QwCopyArray(source.fAbsPos, fAbsPos);
  }

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, QwCopyArray(), QwQPD(), and VQwBPM::VQwBPM().

Here is the call graph for this function:

◆ ~QwQPD()

QwQPD::~QwQPD ( )

inlineoverride

Definition at line 62 of file QwQPD.h.

62{ };

Member Function Documentation

◆ AccumulateRunningSum() [1/2]

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

Definition at line 702 of file QwQPD.cc.

{
  Short_t i = 0;
  for(i=0;i<4;i++) fPhotodiode[i].AccumulateRunningSum(value.fPhotodiode[i], count, ErrorMask);
  for (i = 0; i < 2; i++){    
    fRelPos[i].AccumulateRunningSum(value.fRelPos[i], count, ErrorMask);
    fAbsPos[i].AccumulateRunningSum(value.fAbsPos[i], count, ErrorMask);
  }
  fEffectiveCharge.AccumulateRunningSum(value.fEffectiveCharge, count, ErrorMask);
  return;
 
 
}

References AccumulateRunningSum(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, and QwQPD().

Referenced by AccumulateRunningSum(), and AccumulateRunningSum().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ AccumulateRunningSum() [2/2]

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

overridevirtual

Reimplemented from VQwBPM.

Definition at line 697 of file QwQPD.cc.

{
  AccumulateRunningSum(*dynamic_cast<const QwQPD* >(&value), count, ErrorMask);
}

References AccumulateRunningSum(), QwQPD(), and VQwBPM::VQwBPM().

Here is the call graph for this function:

◆ ApplyHWChecks()

Bool_t QwQPD::ApplyHWChecks ( )

Apply hardware-level checks to all four photodiode channels and aggregate their status.

Returns: true if all photodiodes pass hardware checks.

Definition at line 129 of file QwQPD.cc.

{
  Bool_t eventokay=kTRUE;
 
  UInt_t deviceerror=0;
  for(Short_t i=0;i<4;i++)
    {
      deviceerror|= fPhotodiode[i].ApplyHWChecks();  //OR the error code from each wire
      eventokay &= (deviceerror & 0x0);//AND with 0 since zero means HW is good.
 
      if (bDEBUG) std::cout<<" Inconsistent within QPD terminals photodiode[ "<<i<<" ] "<<std::endl;
      if (bDEBUG) std::cout<<" photodiode[ "<<i<<" ] sequence num "<<fPhotodiode[i].GetSequenceNumber()<<" sample size "<<fPhotodiode[i].GetNumberOfSamples()<<std::endl;
    }
 
  return eventokay;
}

References VQwBPM::bDEBUG, and fPhotodiode.

Referenced by ProcessEvent().

Here is the caller graph for this function:

◆ ApplySingleEventCuts()

Bool_t QwQPD::ApplySingleEventCuts ( )

overridevirtual

Apply single-event cuts to photodiodes and derived channels, propagating error codes to position and charge outputs.

Returns: true if all channels pass their cuts.

Implements VQwBPM.

Definition at line 219 of file QwQPD.cc.

{
  Bool_t status=kTRUE;
  Int_t i=0;
  UInt_t error_code = 0;
  //Event cuts for four pads
  for(i=0;i<4;i++){
    if (fPhotodiode[i].ApplySingleEventCuts()){ //for RelX
      status&=kTRUE;
    }
    else{
      status&=kFALSE;
      if (bDEBUG) std::cout<<" single pad "<<fPhotodiode[i].GetElementName()<<" event cut failed ";
    }
    //Get the Event cut error flag for the pads
    error_code|=fPhotodiode[i].GetErrorCode();
  }
 
  for(i=0;i<2;i++){
    fRelPos[i].UpdateErrorFlag(error_code);//To update the rel/abs error code from the channels/wires event cut error codes
    status &= fRelPos[i].ApplySingleEventCuts();
 
    fAbsPos[i].UpdateErrorFlag(error_code);//To update the rel/abs error code from the channels/wires event cut error codes
    status &= fAbsPos[i].ApplySingleEventCuts();
  }
  fEffectiveCharge.UpdateErrorFlag(error_code);// To update the eff-charge error code from the channels/wires event cut error codes
  status &= fEffectiveCharge.ApplySingleEventCuts();
  
  return status;
}

References ApplySingleEventCuts(), VQwBPM::bDEBUG, fAbsPos, fEffectiveCharge, fPhotodiode, and fRelPos.

Referenced by ApplySingleEventCuts().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ CalculateRunningAverage()

void QwQPD::CalculateRunningAverage ( )

overridevirtual

Implements VQwBPM.

Definition at line 685 of file QwQPD.cc.

{
  Short_t i = 0;
  for(i=0;i<4;i++) fPhotodiode[i].CalculateRunningAverage();
  for (i = 0; i < 2; i++){
    fRelPos[i].CalculateRunningAverage();
    fAbsPos[i].CalculateRunningAverage();
  }
  fEffectiveCharge.CalculateRunningAverage();
  return;
}

References CalculateRunningAverage(), fAbsPos, fEffectiveCharge, fPhotodiode, and fRelPos.

Referenced by CalculateRunningAverage().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ CheckForBurpFail()

Bool_t QwQPD::CheckForBurpFail ( const VQwDataElement * ev_error )

overridevirtual

Polymorphic burp check against a reference QPD, delegating to all photodiodes, position, and charge channels.

Returns: true if any subchannel reports a burp failure.

Exceptions

std::invalid_argument if ev_error is not a QwQPD.

Reimplemented from VQwBPM.

Definition at line 383 of file QwQPD.cc.

                                                            {
  Short_t i=0;
  Bool_t burpstatus = kFALSE;
  try {
    if(typeid(*ev_error)==typeid(*this)) {
      //std::cout<<" Here in QwQPD::CheckForBurpFail \n";
      if (this->GetElementName()!="") {
        const QwQPD* value_qpd = dynamic_cast<const QwQPD* >(ev_error);
        for(i=0;i<4;i++){
          burpstatus |= fPhotodiode[i].CheckForBurpFail(&(value_qpd->fPhotodiode[i]));
        }
        for(i=kXAxis;i<kNumAxes;i++) {
          burpstatus |= fRelPos[i].CheckForBurpFail(&(value_qpd->fRelPos[i]));
          burpstatus |= fAbsPos[i].CheckForBurpFail(&(value_qpd->fAbsPos[i])); 
        }
        burpstatus |= fEffectiveCharge.CheckForBurpFail(&(value_qpd->fEffectiveCharge)); 
      }
    } else {
      TString loc="Standard exception from QwQPD::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 fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, QwQPD(), and VQwDataElement::VQwDataElement().

Here is the call graph for this function:

◆ ClearEventData()

void QwQPD::ClearEventData ( )

overridevirtual

Clear event data for all photodiodes and derived position/charge channels.

Reimplemented from VQwDataElement.

Definition at line 106 of file QwQPD.cc.

{
  Short_t i=0;
 
  for(i=0;i<4;i++) fPhotodiode[i].ClearEventData();
 
  for(i=kXAxis;i<kNumAxes;i++) {
    fRelPos[i].ClearEventData();
    fAbsPos[i].ClearEventData();
  }
  
  fEffectiveCharge.ClearEventData();
 
 return;
}

References ClearEventData(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwBPM::kNumAxes, and VQwBPM::kXAxis.

Referenced by ClearEventData().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ConstructBranch() [1/2]

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

overridevirtual

Implements VQwBPM.

Definition at line 807 of file QwQPD.cc.

{
  if (GetElementName()==""){
    //  This channel is not used, so skip constructing trees.
  }
  else {
    TString thisprefix=prefix;
    if(prefix.Contains("asym_"))
      thisprefix.ReplaceAll("asym_","diff_");
 
    SetRootSaveStatus(prefix);
    
    fEffectiveCharge.ConstructBranch(tree,prefix);
    Short_t i = 0;
    if(bFullSave) {
      for(i=0;i<4;i++) fPhotodiode[i].ConstructBranch(tree,thisprefix);
    }
    for(i=kXAxis;i<kNumAxes;i++) {
      fAbsPos[i].ConstructBranch(tree,thisprefix);
      fRelPos[i].ConstructBranch(tree,thisprefix);
    }
 
  }
  return;
}

References VQwBPM::bFullSave, ConstructBranch(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, and VQwBPM::SetRootSaveStatus().

Referenced by ConstructBranch(), and ConstructBranch().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ConstructBranch() [2/2]

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

overridevirtual

Implements VQwBPM.

Definition at line 833 of file QwQPD.cc.

{
  TString devicename;
 
  devicename=GetElementName();
  devicename.ToLower();
  if (GetElementName()==""){
    //  This channel is not used, so skip filling the histograms.
  } else
    {
      if (modulelist.HasValue(devicename)){
    TString thisprefix=prefix;
    if(prefix.Contains("asym_"))
      thisprefix.ReplaceAll("asym_","diff_");
 
    SetRootSaveStatus(prefix);
 
    fEffectiveCharge.ConstructBranch(tree,prefix);
    Short_t i = 0;
    if(bFullSave) {
      for(i=0;i<4;i++) fPhotodiode[i].ConstructBranch(tree,thisprefix);
    }
    for(i=kXAxis;i<kNumAxes;i++) {
      fAbsPos[i].ConstructBranch(tree,thisprefix);
      fRelPos[i].ConstructBranch(tree,thisprefix);
    }
 
    QwMessage <<" Tree leaves added to "<<devicename<<" Corresponding channels"<<QwLog::endl;
      }
      // this functions doesn't do anything yet
    }
 
 
 
 
 
  return;
}

References VQwBPM::bFullSave, ConstructBranch(), QwLog::endl(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), QwParameterFile::HasValue(), VQwBPM::kNumAxes, VQwBPM::kXAxis, QwMessage, and VQwBPM::SetRootSaveStatus().

Here is the call graph for this function:

◆ ConstructBranchAndVector()

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

overridevirtual

Implements VQwBPM.

Definition at line 781 of file QwQPD.cc.

{
  if (GetElementName()==""){
    //  This channel is not used, so skip constructing trees.
  }
  else {
    TString thisprefix=prefix;
    if(prefix.Contains("asym_"))
      thisprefix.ReplaceAll("asym_","diff_");
 
    SetRootSaveStatus(prefix);
 
    fEffectiveCharge.ConstructBranchAndVector(tree,prefix,values);
    Short_t i = 0;
    if(bFullSave) {
      for(i=0;i<4;i++) fPhotodiode[i].ConstructBranchAndVector(tree,thisprefix,values);
    }
    for(i=kXAxis;i<kNumAxes;i++) {
      fAbsPos[i].ConstructBranchAndVector(tree,thisprefix,values);
      fRelPos[i].ConstructBranchAndVector(tree,thisprefix,values);
    }
 
  }
  return;
}

References VQwBPM::bFullSave, ConstructBranchAndVector(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, and VQwBPM::SetRootSaveStatus().

Referenced by ConstructBranchAndVector().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ConstructHistograms()

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

overridevirtual

Construct the histograms for this data element.

Implements VQwBPM.

Definition at line 738 of file QwQPD.cc.

{
 
  if (GetElementName()=="") {
    //  This channel is not used, so skip filling the histograms.
  }  else {
    fEffectiveCharge.ConstructHistograms(folder, prefix);
    TString thisprefix=prefix;
 
    if(prefix.Contains("asym_"))
      thisprefix.ReplaceAll("asym_","diff_");
    SetRootSaveStatus(prefix);
    Short_t i = 0;
    if(bFullSave) {
      for(i=0;i<4;i++) fPhotodiode[i].ConstructHistograms(folder, thisprefix);
    }
    for(i=kXAxis;i<kNumAxes;i++) {
      fAbsPos[i].ConstructHistograms(folder, thisprefix);
      fRelPos[i].ConstructHistograms(folder, thisprefix);
    }
  }
  return;
}

References VQwBPM::bFullSave, ConstructHistograms(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, and VQwBPM::SetRootSaveStatus().

Referenced by ConstructHistograms().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ DeaccumulateRunningSum() [1/2]

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

Definition at line 721 of file QwQPD.cc.

{
 
  Short_t i = 0;
  for(i=0;i<4;i++) fPhotodiode[i].DeaccumulateRunningSum(value.fPhotodiode[i]);
  for (i = 0; i < 2; i++){    
    fRelPos[i].DeaccumulateRunningSum(value.fRelPos[i], ErrorMask);
    fAbsPos[i].DeaccumulateRunningSum(value.fAbsPos[i], ErrorMask);
  }
  fEffectiveCharge.DeaccumulateRunningSum(value.fEffectiveCharge, ErrorMask);
  return;
 
 
}

References DeaccumulateRunningSum(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, and QwQPD().

Here is the call graph for this function:

◆ DeaccumulateRunningSum() [2/2]

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

overridevirtual

Implements VQwBPM.

Definition at line 716 of file QwQPD.cc.

{
  DeaccumulateRunningSum(*dynamic_cast<QwQPD* >(&value), ErrorMask);
}

References DeaccumulateRunningSum(), QwQPD(), and VQwBPM::VQwBPM().

Referenced by DeaccumulateRunningSum(), and DeaccumulateRunningSum().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ EncodeEventData()

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

overridevirtual

Reimplemented from VQwBPM.

Definition at line 1043 of file QwQPD.cc.

{
  for (Short_t i=0; i<4; i++) fPhotodiode[i].EncodeEventData(buffer);
}

References EncodeEventData(), and fPhotodiode.

Referenced by EncodeEventData().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ FillHistograms()

void QwQPD::FillHistograms ( )

overridevirtual

Fill the histograms for this data element.

Implements VQwBPM.

Definition at line 762 of file QwQPD.cc.

{
  if (GetElementName()=="") {
    //  This channel is not used, so skip filling the histograms.
  }
  else {
    fEffectiveCharge.FillHistograms();
    Short_t i = 0;
    if(bFullSave) {
      for(i=0;i<4;i++) fPhotodiode[i].FillHistograms();
    }
    for(i=kXAxis;i<kNumAxes;i++){
      fAbsPos[i].FillHistograms();
      fRelPos[i].FillHistograms();
    }
  }
  return;
}

References VQwBPM::bFullSave, fAbsPos, fEffectiveCharge, FillHistograms(), fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, and VQwBPM::kXAxis.

Referenced by FillHistograms().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ FillTreeVector()

void QwQPD::FillTreeVector ( QwRootTreeBranchVector & values ) const

overridevirtual

Implements VQwBPM.

Definition at line 872 of file QwQPD.cc.

{
  if (GetElementName()=="") {
    //  This channel is not used, so skip filling the tree.
  }
  else {
    fEffectiveCharge.FillTreeVector(values);
    Short_t i = 0;
    if(bFullSave) {
      for(i=0;i<4;i++) fPhotodiode[i].FillTreeVector(values);
    }
 
    for(i=kXAxis;i<kNumAxes;i++){
      fAbsPos[i].FillTreeVector(values);
      fRelPos[i].FillTreeVector(values);
    }
  }
  return;
}

References VQwBPM::bFullSave, fAbsPos, fEffectiveCharge, FillTreeVector(), fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, and VQwBPM::kXAxis.

Referenced by FillTreeVector().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetAbsolutePosition()

void QwQPD::GetAbsolutePosition ( )

inlineoverridevirtual

Reimplemented from VQwBPM.

Definition at line 91 of file QwQPD.h.

91{};

◆ GetCalibrationFactors()

void QwQPD::GetCalibrationFactors	(	Double_t	AlphaX,
		Double_t	AlphaY )

Apply QPD-specific calibration factors converting ADC counts to mm. These come from geometry/calibration files for each axis.

Definition at line 84 of file QwQPD.cc.

{
  // Read in the calibration factors from the injector_beamline_geometry.map
  // for the QPD, AlphaX and AlphaY gives the conversion from adc counts to mm.
 
  Bool_t ldebug = kFALSE;
 
  fQwQPDCalibration[0]=1.0/AlphaX;
  fQwQPDCalibration[1]=1.0/AlphaY;
 
  if(ldebug){
    std::cout<<"\n%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n";
    std::cout<<this->GetElementName();
    std::cout<<"\nfQwQPDCalibration[0]  = "<<fQwQPDCalibration[0]<<std::endl;
    std::cout<<"\nfQwQPDCalibration[1]  = "<<fQwQPDCalibration[1]<<std::endl;
    std::cout<<"AlphaX = "<<fRelativeGains[0]<<std::endl;
    std::cout<<"AlphaY = "<<fRelativeGains[1]<<std::endl;    
  }
  return;
}

References fQwQPDCalibration, VQwBPM::fRelativeGains, and VQwDataElement::GetElementName().

Here is the call graph for this function:

◆ GetEffectiveCharge()

const VQwHardwareChannel * QwQPD::GetEffectiveCharge ( ) const

inlineoverridevirtual

Implements VQwBPM.

Definition at line 88 of file QwQPD.h.

88{return &fEffectiveCharge;}

References fEffectiveCharge.

◆ GetEventcutErrorFlag()

UInt_t QwQPD::GetEventcutErrorFlag ( )

overridevirtual

Return OR of event-cut error flags across all photodiodes and derived channels.

Reimplemented from VQwDataElement.

Definition at line 173 of file QwQPD.cc.

{
  Short_t i=0;
  UInt_t error=0;
  for(i=0;i<4;i++) 
    error|=fPhotodiode[i].GetEventcutErrorFlag();
 
  for(i=kXAxis;i<kNumAxes;i++) {
    error|=fRelPos[i].GetEventcutErrorFlag();
    error|=fAbsPos[i].GetEventcutErrorFlag();
  }
 
  error|=fEffectiveCharge.GetEventcutErrorFlag();
  return error;
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, GetEventcutErrorFlag(), VQwBPM::kNumAxes, and VQwBPM::kXAxis.

Referenced by GetEventcutErrorFlag().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetPosition()

const VQwHardwareChannel * QwQPD::GetPosition ( EBeamPositionMonitorAxis axis ) const

inlineoverridevirtual

Reimplemented from VQwBPM.

Definition at line 80 of file QwQPD.h.

                                                                                      {
    if (axis<0 || axis>2){
      TString loc="QwQPD::GetPosition for "
        +this->GetElementName()+" failed for axis value "+Form("%d",axis);
      throw std::out_of_range(loc.Data());
    }
    return &fAbsPos[axis];
  }

References VQwBPM::axis, fAbsPos, and VQwDataElement::GetElementName().

Here is the call graph for this function:

◆ GetSubelementByName()

VQwHardwareChannel * QwQPD::GetSubelementByName ( TString ch_name )

overrideprotectedvirtual

Map a human-readable name to the corresponding channel pointer. Supports: tl/tr/br/bl (photodiodes), relx/rely, absx/x, absy/y, effectivecharge/charge.

Exceptions

std::invalid_argument on unrecognized names.

Implements VQwBPM.

Definition at line 257 of file QwQPD.cc.

{
  VQwHardwareChannel* tmpptr = NULL;
  ch_name.ToLower();
  if (ch_name=="tl"){
    tmpptr = &fPhotodiode[0];
  }else if (ch_name=="tr"){
    tmpptr = &fPhotodiode[1];
  }else if (ch_name=="br"){
    tmpptr = &fPhotodiode[2];
  }else if (ch_name=="bl"){
    tmpptr = &fPhotodiode[3];
  }else if (ch_name=="relx"){
    tmpptr = &fRelPos[0];    
  }else if (ch_name=="rely"){
    tmpptr = &fRelPos[1];
  }else  if (ch_name=="absx" || ch_name=="x" ){
    tmpptr = &fAbsPos[0];
  }else if (ch_name=="absy" || ch_name=="y"){
    tmpptr = &fAbsPos[1];
  }else if (ch_name=="effectivecharge" || ch_name=="charge"){
    tmpptr = &fEffectiveCharge;
  } else {
    TString loc="QwQPD::GetSubelementByName for"
      + this->GetElementName() + " was passed "
      + ch_name + ", which is an unrecognized subelement name.";
    throw std::invalid_argument(loc.Data());
  }
  return tmpptr;
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, and VQwDataElement::GetElementName().

Here is the call graph for this function:

◆ GetSubElementIndex()

UInt_t QwQPD::GetSubElementIndex ( TString subname )

static

Definition at line 584 of file QwQPD.cc.

{
  subname.ToUpper();
  UInt_t localindex=kInvalidSubelementIndex;
  for(Short_t i=0;i<4;i++) if(subname==subelement[i])localindex=i;
 
  if(localindex>3)
    std::cerr << "QwQPD::GetSubElementIndex is unable to associate the string -"
          <<subname<<"- to any index"<<std::endl;
 
  return localindex;
}

References kInvalidSubelementIndex, and subelement.

Referenced by QwBeamDetectorID::QwBeamDetectorID().

Here is the caller graph for this function:

◆ GetSubElementName()

TString QwQPD::GetSubElementName ( Int_t subindex )

overridevirtual

Reimplemented from VQwBPM.

Definition at line 572 of file QwQPD.cc.

{
  TString thisname;
  if(subindex<4&&subindex>-1)
    thisname=fPhotodiode[subindex].GetElementName();
  else
    std::cerr<<"QwQPD::GetSubElementName for "<<
      GetElementName()<<" this subindex doesn't exists \n";
 
  return thisname;
}

References fPhotodiode, and VQwDataElement::GetElementName().

Here is the call graph for this function:

◆ IncrementErrorCounters()

void QwQPD::IncrementErrorCounters ( )

overridevirtual

Increment persistent error counters for all channels.

Implements VQwBPM.

Definition at line 147 of file QwQPD.cc.

{
  Short_t i=0;
  for(i=0;i<4;i++) 
    fPhotodiode[i].IncrementErrorCounters();
  for(i=kXAxis;i<kNumAxes;i++) {
    fRelPos[i].IncrementErrorCounters();
    fAbsPos[i].IncrementErrorCounters();
  }
  fEffectiveCharge.IncrementErrorCounters();
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, IncrementErrorCounters(), VQwBPM::kNumAxes, and VQwBPM::kXAxis.

Referenced by IncrementErrorCounters().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ InitializeChannel() [1/2]

void QwQPD::InitializeChannel ( TString name )

Initialize QPD with 4 photodiode raw channels and derived position/charge.

Definition at line 31 of file QwQPD.cc.

{
  Short_t i=0;
  Bool_t localdebug = kFALSE;
 
  VQwBPM::InitializeChannel(name);
 
  fEffectiveCharge.InitializeChannel(name+"WS","derived");
 
  for(i=0;i<4;i++) {
    fPhotodiode[i].InitializeChannel(name+subelement[i],"raw");
    
    if(localdebug)
      std::cout<<" photodiode ["<<i<<"]="<<fPhotodiode[i].GetElementName()<<"\n";
  }
  
  for(i=kXAxis;i<kNumAxes;i++) {
    fRelPos[i].InitializeChannel(name+"Rel"+kAxisLabel[i],"derived");
    fAbsPos[i].InitializeChannel(name+kAxisLabel[i],"derived");
  }
  
  bFullSave=kTRUE;
 
  return;
}

References VQwBPM::bFullSave, fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwBPM::InitializeChannel(), VQwBPM::kAxisLabel, VQwBPM::kNumAxes, VQwBPM::kXAxis, and subelement.

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

Here is the call graph for this function:

Here is the caller graph for this function:

◆ InitializeChannel() [2/2]

void QwQPD::InitializeChannel	(	TString	subsystem,
		TString	name )

Initialize QPD with subsystem scoping for proper tree/histogram naming.

Definition at line 58 of file QwQPD.cc.

{
  Short_t i=0;
  Bool_t localdebug = kFALSE;
 
  VQwBPM::InitializeChannel(name);
 
  fEffectiveCharge.InitializeChannel(subsystem, "QwQPD", name+"WS","derived");
 
  for(i=0;i<4;i++) {
    fPhotodiode[i].InitializeChannel(subsystem, "QwQPD", name+subelement[i],"raw");
    if(localdebug)
      std::cout<<" photodiode ["<<i<<"]="<<fPhotodiode[i].GetElementName()<<"\n";
  }
  for(i=kXAxis;i<kNumAxes;i++) fRelPos[i].InitializeChannel(subsystem, "QwQPD", name+"Rel"+kAxisLabel[i],"derived");
  for(i=kXAxis;i<kNumAxes;i++) fAbsPos[i].InitializeChannel(subsystem, "QwQPD", name+kAxisLabel[i],"derived");
 
  bFullSave=kTRUE;
 
  return;
}

References VQwBPM::bFullSave, fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, InitializeChannel(), VQwBPM::InitializeChannel(), VQwBPM::kAxisLabel, VQwBPM::kNumAxes, VQwBPM::kXAxis, and subelement.

Here is the call graph for this function:

◆ LoadChannelParameters()

void QwQPD::LoadChannelParameters ( QwParameterFile & paramfile )

inlineoverridevirtual

Implements VQwBPM.

Definition at line 68 of file QwQPD.h.

                                                                 {
    for(Short_t i=0;i<4;i++)
      fPhotodiode[i].LoadChannelParameters(paramfile);
  }

References fPhotodiode, and LoadChannelParameters().

Referenced by LoadChannelParameters().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ MakeQPDList()

void QwQPD::MakeQPDList ( )

Definition at line 955 of file QwQPD.cc.

{
  for (size_t i = kXAxis; i < kNumAxes; i++) {
    QwVQWK_Channel abspos(fAbsPos[i]);
    abspos = fAbsPos[i];
    fQPDElementList.push_back(abspos);
  }
  QwVQWK_Channel effectivecharge(fEffectiveCharge);
  effectivecharge = fEffectiveCharge;
  fQPDElementList.push_back(effectivecharge);
}

References fAbsPos, fEffectiveCharge, fQPDElementList, VQwBPM::kNumAxes, and VQwBPM::kXAxis.

◆ operator+=() [1/2]

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

virtual

Definition at line 623 of file QwQPD.cc.

{
  if (GetElementName()!=""){
    Short_t i = 0;
    this->fEffectiveCharge+=value.fEffectiveCharge;
    for(i=0;i<4;i++) this->fPhotodiode[i]+=value.fPhotodiode[i];
 
    for(i=kXAxis;i<kNumAxes;i++){
      this->fRelPos[i]+=value.fRelPos[i];
      this->fAbsPos[i]+=value.fAbsPos[i];
    }
  }
  return *this;
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, and QwQPD().

Here is the call graph for this function:

◆ operator+=() [2/2]

VQwBPM & QwQPD::operator+= ( const VQwBPM & value )

overridevirtual

Implements VQwBPM.

Definition at line 617 of file QwQPD.cc.

{
  *(dynamic_cast<QwQPD*>(this)) += *(dynamic_cast<const QwQPD*>(&value));
  return *this;
}

References QwQPD(), and VQwBPM::VQwBPM().

Here is the call graph for this function:

◆ operator-=() [1/2]

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

virtual

Definition at line 644 of file QwQPD.cc.

{
  if (GetElementName()!=""){
    Short_t i = 0;
    this->fEffectiveCharge-=value.fEffectiveCharge;
    for(i=0;i<4;i++) this->fPhotodiode[i]-=value.fPhotodiode[i];
 
    for(i=kXAxis;i<kNumAxes;i++){
      this->fRelPos[i]-=value.fRelPos[i];
      this->fAbsPos[i]-=value.fAbsPos[i];
    }
  }
  return *this;
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, and QwQPD().

Here is the call graph for this function:

◆ operator-=() [2/2]

VQwBPM & QwQPD::operator-= ( const VQwBPM & value )

overridevirtual

Implements VQwBPM.

Definition at line 638 of file QwQPD.cc.

{
  *(dynamic_cast<QwQPD*>(this)) -= *(dynamic_cast<const QwQPD*>(&value));
  return *this;
}

References QwQPD(), and VQwBPM::VQwBPM().

Here is the call graph for this function:

◆ operator=() [1/2]

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

virtual

Definition at line 603 of file QwQPD.cc.

{
  if (GetElementName()!=""){
    Short_t i = 0;
    this->fEffectiveCharge=value.fEffectiveCharge;
    for(i=0;i<4;i++) this->fPhotodiode[i]=value.fPhotodiode[i];
    for(i=kXAxis;i<kNumAxes;i++){
      this->fRelPos[i]=value.fRelPos[i];
      this->fAbsPos[i]=value.fAbsPos[i];
    }
  }
  return *this;
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, and QwQPD().

Here is the call graph for this function:

◆ operator=() [2/2]

VQwBPM & QwQPD::operator= ( const VQwBPM & value )

overridevirtual

Implements VQwBPM.

Definition at line 597 of file QwQPD.cc.

{
  *(dynamic_cast<QwQPD*>(this)) = *(dynamic_cast<const QwQPD*>(&value));
  return *this;
}

References QwQPD(), and VQwBPM::VQwBPM().

Here is the call graph for this function:

◆ PrintErrorCounters()

void QwQPD::PrintErrorCounters ( ) const

overridevirtual

Print error counter summaries for all channels.

Reimplemented from VQwDataElement.

Definition at line 160 of file QwQPD.cc.

{
  Short_t i=0;
  for(i=0;i<4;i++) 
    fPhotodiode[i].PrintErrorCounters();
  for(i=kXAxis;i<kNumAxes;i++) {
    fRelPos[i].PrintErrorCounters();
    fAbsPos[i].PrintErrorCounters();
  }
  fEffectiveCharge.PrintErrorCounters();
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwBPM::kNumAxes, VQwBPM::kXAxis, and PrintErrorCounters().

Referenced by PrintErrorCounters().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ProcessEvBuffer()

Int_t QwQPD::ProcessEvBuffer	(	UInt_t *	buffer,
		UInt_t	word_position_in_buffer,
		UInt_t	index )

overridevirtual

Route raw buffer data to the specified photodiode channel.

Parameters

buffer	Raw event buffer
word_position_in_buffer	Position in buffer to read from
index	Target photodiode index (0-3)

Returns: Original word_position_in_buffer

Implements VQwDataElement.

Definition at line 557 of file QwQPD.cc.

{
  if(index<4)
    {
      fPhotodiode[index].ProcessEvBuffer(buffer,word_position_in_buffer);
    }
  else
    {
    std::cerr <<
      "QwQPD::ProcessEvBuffer(): attempt to fill in raw date for a wire that doesn't exist \n";
    }
  return word_position_in_buffer;
}

References fPhotodiode.

◆ ProcessEvent()

void QwQPD::ProcessEvent ( )

overridevirtual

Process the current event: apply HW checks, sum photodiodes for effective charge, and calculate X/Y positions using the standard QPD formula: X = ((TL-TR) + (BL-BR)) / (TL+TR+BL+BR) Y = ((TL-BL) + (TR-BR)) / (TL+TR+BL+BR)

The positions X and Y from a QPD are calculated using following equations,

TL ------— TR

| 4 | 2 |

BL | 3 | 1 | BR

     (4-2) + (3-1)                (4-3) + (2-1)
 X =  ------------           Y =  -------------
        1+2+3+4                      1+2+3+4

Implements VQwBPM.

Definition at line 453 of file QwQPD.cc.

{
  Bool_t localdebug = kFALSE;
  static QwVQWK_Channel numer[2];
  static QwVQWK_Channel tmp("tmp");
  static QwVQWK_Channel tmp1("tmp1");
  static QwVQWK_Channel tmp2("tmp2");
 
  numer[0].InitializeChannel("Xnumerator","raw");
  numer[1].InitializeChannel("Ynumerator","raw");
 
  Short_t i = 0;
 
  ApplyHWChecks();
  /* Frst apply HW checks and update HW  error flags. 
     Calling this routine here and not in ApplySingleEventCuts  
     makes a difference for a QPDs because they are derived devices.
  */
 
  fEffectiveCharge.ClearEventData();
 
  for(i=0;i<4;i++){
    fPhotodiode[i].ProcessEvent();
    fEffectiveCharge+=fPhotodiode[i];
  }
  
  if (localdebug) fEffectiveCharge.PrintInfo();
   
 
  /** The positions X and Y from a QPD are calculated using following equations,
 
 TL   ---------  TR
     | 4  |  2 |
      ---------
 BL  | 3  |  1 | BR
      ---------
 
 
          (4-2) + (3-1)                (4-3) + (2-1)
      X =  ------------           Y =  -------------
             1+2+3+4                      1+2+3+4
  **/
 
  if(localdebug){
    std::cout<<"#############################\n";
    std::cout<<" QPD name = "<<fElementName<<std::endl;
    std::cout<<" event number = "<<fPhotodiode[0].GetSequenceNumber()<<"\n";
    std::cout<<" hw  BR ="<<fPhotodiode[0].GetValue()<<"\n";
    std::cout<<" hw  TR ="<<fPhotodiode[1].GetValue()<<"\n";
    std::cout<<" hw  BL ="<<fPhotodiode[2].GetValue()<<"\n";
    std::cout<<" hw  TL ="<<fPhotodiode[3].GetValue()<<"\n\n";
  }
 
  // X numerator
  tmp1.ClearEventData();
  tmp1.Difference(fPhotodiode[3],fPhotodiode[1]);  // 4-2
  tmp2.ClearEventData();
  tmp2.Difference(fPhotodiode[2],fPhotodiode[0]);  // 3-1
  tmp.ClearEventData();
  numer[0].Sum(tmp1,tmp2);
 
  // Y numerator
  tmp1.ClearEventData();
  tmp1.Difference(fPhotodiode[3],fPhotodiode[2]);  // 4-3
  tmp2.ClearEventData();
  tmp2.Difference(fPhotodiode[1],fPhotodiode[0]);  // 2-1
  tmp.ClearEventData();
  numer[1].Sum(tmp1,tmp2);
 
  for(i=kXAxis;i<kNumAxes;i++){
    tmp.ClearEventData();
    tmp.Sum(fPhotodiode[0],fPhotodiode[1]);
    tmp1.ClearEventData();
    tmp1.Sum(fPhotodiode[2],fPhotodiode[3]);
    fEffectiveCharge.Sum(tmp,tmp1);
 
    // X/Y reading in ADC counts
    fRelPos[i].Ratio(numer[i],fEffectiveCharge);
 
    // X/Y reading in mm.
    fAbsPos[i]=fRelPos[i];
    fAbsPos[i].Scale(fQwQPDCalibration[i]);
 
    if(localdebug){
      std::cout<<" hw  numerator= "<<numer[i].GetValue()<<"  ";
      std::cout<<" hw  denominator (== Effective_Charge)= "<<fEffectiveCharge.GetValue()<<"\n";
      std::cout<<" hw  clibration factors= "<<fQwQPDCalibration[i]<<"\n";
      std::cout<<" hw  fRelPos["<<kAxisLabel[i]<<"]="<<fRelPos[i].GetValue()<<"\n \n";
      std::cout<<" hw  fAbsPos["<<kAxisLabel[i]<<"]="<<fAbsPos[i].GetValue()<<"\n \n";
    }
  }
  
  return;
}

References ApplyHWChecks(), QwVQWK_Channel::ClearEventData(), QwVQWK_Channel::Difference(), fAbsPos, fEffectiveCharge, VQwDataElement::fElementName, fPhotodiode, fQwQPDCalibration, fRelPos, QwVQWK_Channel::GetValue(), QwVQWK_Channel::InitializeChannel(), VQwBPM::kAxisLabel, VQwBPM::kNumAxes, VQwBPM::kXAxis, and QwVQWK_Channel::Sum().

Here is the call graph for this function:

◆ RandomizeEventData()

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

overridevirtual

Reimplemented from VQwBPM.

Definition at line 1024 of file QwQPD.cc.

{
  for (Short_t i=0; i<4; i++) fPhotodiode[i].RandomizeEventData(helicity, time);
 
  return;
}

References fPhotodiode, and RandomizeEventData().

Referenced by RandomizeEventData().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ Ratio()

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

Definition at line 660 of file QwQPD.cc.

{
  // this function is called when forming asymmetries. In this case what we actually want for the
  // QPD is the difference only not the asymmetries
 
  *this=numer;
  this->fEffectiveCharge.Ratio(numer.fEffectiveCharge,denom.fEffectiveCharge);
  return;
}

References fEffectiveCharge, and QwQPD().

Here is the call graph for this function:

◆ Scale()

void QwQPD::Scale ( Double_t factor )

overridevirtual

Reimplemented from VQwBPM.

Definition at line 672 of file QwQPD.cc.

{
  Short_t i = 0;
 
  for(i=0;i<4;i++) fPhotodiode[i].Scale(factor);
  for(i=kXAxis;i<kNumAxes;i++){
    fRelPos[i].Scale(factor);
    fAbsPos[i].Scale(factor);
  }
  fEffectiveCharge.Scale(factor);
  return;
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwBPM::kNumAxes, VQwBPM::kXAxis, and Scale().

Referenced by Scale().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ SetDefaultSampleSize()

void QwQPD::SetDefaultSampleSize ( Int_t sample_size )

overridevirtual

Reimplemented from VQwBPM.

Definition at line 1049 of file QwQPD.cc.

{
  for(Short_t i=0;i<4;i++) fPhotodiode[i].SetDefaultSampleSize((size_t)sample_size);
  return;
}

References fPhotodiode, and SetDefaultSampleSize().

Referenced by SetDefaultSampleSize().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ SetEventCutMode()

void QwQPD::SetEventCutMode ( Int_t bcuts )

overridevirtual

Implements VQwBPM.

Definition at line 942 of file QwQPD.cc.

{
  Short_t i = 0;
  // bEVENTCUTMODE=bcuts;
  for (i=0;i<4;i++) fPhotodiode[i].SetEventCutMode(bcuts);
  for (i=kXAxis;i<kNumAxes;i++) {
    fAbsPos[i].SetEventCutMode(bcuts);
    fRelPos[i].SetEventCutMode(bcuts);
  }
  fEffectiveCharge.SetEventCutMode(bcuts);
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwBPM::kNumAxes, VQwBPM::kXAxis, and SetEventCutMode().

Referenced by SetEventCutMode().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ SetEventData()

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

Definition at line 1032 of file QwQPD.cc.

{
  for (Short_t i=0; i<2; i++)
    {
      fAbsPos[i].SetHardwareSum(relpos[i], sequencenumber);
    }
 
  return;
}

References fAbsPos.

◆ SetRandomEventParameters()

void QwQPD::SetRandomEventParameters	(	Double_t	meanX,
		Double_t	sigmaX,
		Double_t	meanY,
		Double_t	sigmaY )

overridevirtual

Reimplemented from VQwBPM.

Definition at line 996 of file QwQPD.cc.

{
  // Average values of the signals in the stripline ADCs
  Double_t sumX = 1.1e8; // These are just guesses, but I made X and Y different
  Double_t sumY = 0.9e8; // to make it more interesting for the analyzer...
 
  
  // Determine the asymmetry from the position
  Double_t meanXP = (1.0 + meanX / fQwQPDCalibration[0]) * sumX / 2.0;
  Double_t meanXM = (1.0 - meanX / fQwQPDCalibration[0]) * sumX / 2.0; // = sumX - meanXP;
  Double_t meanYP = (1.0 + meanY / fQwQPDCalibration[1]) * sumY / 2.0;
  Double_t meanYM = (1.0 - meanY / fQwQPDCalibration[1]) * sumY / 2.0; // = sumY - meanYP;
 
  // Determine the spread of the asymmetry (this is not tested yet)
  // (negative sigma should work in the QwVQWK_Channel, but still using fabs)
  Double_t sigmaXP = fabs(sumX * sigmaX / meanX);
  Double_t sigmaXM = sigmaXP;
  Double_t sigmaYP = fabs(sumY * sigmaY / meanY);
  Double_t sigmaYM = sigmaYP;
 
  // Propagate these parameters to the ADCs
  fPhotodiode[0].SetRandomEventParameters(meanXP, sigmaXP);
  fPhotodiode[1].SetRandomEventParameters(meanXM, sigmaXM);
  fPhotodiode[2].SetRandomEventParameters(meanYP, sigmaYP);
  fPhotodiode[3].SetRandomEventParameters(meanYM, sigmaYM);
}

References fPhotodiode, and fQwQPDCalibration.

◆ SetSingleEventCuts()

void QwQPD::SetSingleEventCuts	(	TString	ch_name,
		UInt_t	errorflag,
		Double_t	minX,
		Double_t	maxX,
		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.

Configure single-event cuts for a specific subelement by name, including error flags, limits, stability, and burp detection parameters.

Parameters

ch_name	Subelement name (tl,tr,bl,br,relx,rely,absx,absy,effectivecharge)
errorflag	Error flag mask to set when the cut fails
minX	Lower limit
maxX	Upper limit
stability	Stability cut width
burplevel	Burp detection threshold

Definition at line 334 of file QwQPD.cc.

                                                                                                                                    {
  errorflag|=kBPMErrorFlag;//update the device flag
  if (ch_name=="tl"){
    QwMessage<<"TL LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fPhotodiode[0].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }else if (ch_name=="tr"){
    QwMessage<<"TR LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fPhotodiode[1].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }else if (ch_name=="br"){
    QwMessage<<"BR LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fPhotodiode[2].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }else if (ch_name=="bl"){
    QwMessage<<"BL LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fPhotodiode[3].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }else if (ch_name=="relx"){
    QwMessage<<"RelX LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fRelPos[0].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }else if (ch_name=="rely"){
    QwMessage<<"RelY LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fRelPos[1].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }else if (ch_name=="absx"){
    QwMessage<<"AbsX LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fAbsPos[0].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }else if (ch_name=="absy"){
    QwMessage<<"AbsY LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fAbsPos[1].SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
    
  }else if (ch_name=="effectivecharge"){
    QwMessage<<"EffectveQ LL " <<  minX <<" UL " << maxX <<QwLog::endl;
    fEffectiveCharge.SetSingleEventCuts(errorflag,minX,maxX,stability,burplevel); 
 
  }
 
}

References QwLog::endl(), fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, kBPMErrorFlag, and QwMessage.

Here is the call graph for this function:

◆ SetSubElementCalibrationFactor()

void QwQPD::SetSubElementCalibrationFactor	(	Int_t	j,
		Double_t	value )

overridevirtual

Reimplemented from VQwBPM.

Definition at line 1062 of file QwQPD.cc.

{
  fPhotodiode[j].SetCalibrationFactor(value);
  return;
}

References fPhotodiode.

◆ SetSubElementPedestal()

void QwQPD::SetSubElementPedestal	(	Int_t	j,
		Double_t	value )

overridevirtual

Reimplemented from VQwBPM.

Definition at line 1056 of file QwQPD.cc.

{
  fPhotodiode[j].SetPedestal(value);
  return;
}

References fPhotodiode.

◆ UpdateErrorFlag() [1/2]

UInt_t QwQPD::UpdateErrorFlag ( )

overridevirtual

Update derived channel error flags based on photodiode error codes and return aggregated event-cut flags.

Implements VQwBPM.

Definition at line 193 of file QwQPD.cc.

{
  Short_t i=0;
  UInt_t error1=0;
  UInt_t error2=0;
  for(i=0;i<4;i++){
    error1|=fPhotodiode[i].GetErrorCode();
    error2|=fPhotodiode[i].GetEventcutErrorFlag();
  }
  for(i=kXAxis;i<kNumAxes;i++) {
    fRelPos[i].UpdateErrorFlag(error1);
    fAbsPos[i].UpdateErrorFlag(error1);
    error2|=fRelPos[i].GetEventcutErrorFlag();
    error2|=fAbsPos[i].GetEventcutErrorFlag();
  }
  fEffectiveCharge.UpdateErrorFlag(error1);
  error2|=fEffectiveCharge.GetEventcutErrorFlag();
  return error2;
}

References fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwBPM::kNumAxes, and VQwBPM::kXAxis.

◆ UpdateErrorFlag() [2/2]

void QwQPD::UpdateErrorFlag ( const VQwBPM * ev_error )

overridevirtual

Copy error flags from a reference QPD to this instance.

Exceptions

std::invalid_argument if ev_error is not a QwQPD.

Implements VQwBPM.

Definition at line 417 of file QwQPD.cc.

                                                 {
  Short_t i=0;
  try {
    if(typeid(*ev_error)==typeid(*this)) {
      // std::cout<<" Here in QwQPD::UpdateErrorFlag \n";
      if (this->GetElementName()!="") {
        const QwQPD* value_bpm = dynamic_cast<const QwQPD* >(ev_error);
    for(i=0;i<4;i++){
      fPhotodiode[i].UpdateErrorFlag(value_bpm->fPhotodiode[i]);
    }
    for(i=kXAxis;i<kNumAxes;i++) {
      fRelPos[i].UpdateErrorFlag(value_bpm->fRelPos[i]);
      fAbsPos[i].UpdateErrorFlag(value_bpm->fAbsPos[i]);
    }
    fEffectiveCharge.UpdateErrorFlag(value_bpm->fEffectiveCharge);
      }
    } else {
      TString loc="Standard exception from QwQPD::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 fAbsPos, fEffectiveCharge, fPhotodiode, fRelPos, VQwDataElement::GetElementName(), VQwBPM::kNumAxes, VQwBPM::kXAxis, QwQPD(), and VQwBPM::VQwBPM().

Here is the call graph for this function:

Field Documentation

◆ fAbsPos

std::array<QwVQWK_Channel,2> QwQPD::fAbsPos

protected

◆ fEffectiveCharge

QwVQWK_Channel QwQPD::fEffectiveCharge

protected

◆ fPhotodiode

std::array<QwVQWK_Channel,4> QwQPD::fPhotodiode

protected

Definition at line 167 of file QwQPD.h.

Referenced by AccumulateRunningSum(), ApplyHWChecks(), ApplySingleEventCuts(), CalculateRunningAverage(), CheckForBurpFail(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), DeaccumulateRunningSum(), EncodeEventData(), FillHistograms(), FillTreeVector(), GetEventcutErrorFlag(), GetSubelementByName(), GetSubElementName(), IncrementErrorCounters(), InitializeChannel(), InitializeChannel(), LoadChannelParameters(), operator+=(), operator-=(), operator=(), PrintErrorCounters(), ProcessEvBuffer(), ProcessEvent(), QwQPD(), RandomizeEventData(), Scale(), SetDefaultSampleSize(), SetEventCutMode(), SetRandomEventParameters(), SetSingleEventCuts(), SetSubElementCalibrationFactor(), SetSubElementPedestal(), UpdateErrorFlag(), and UpdateErrorFlag().

◆ fQPDElementList

std::vector<QwVQWK_Channel> QwQPD::fQPDElementList

protected

Definition at line 175 of file QwQPD.h.

Referenced by MakeQPDList().

◆ fQwQPDCalibration

Double_t QwQPD::fQwQPDCalibration[2]

private

Definition at line 163 of file QwQPD.h.

Referenced by GetCalibrationFactors(), ProcessEvent(), QwQPD(), and SetRandomEventParameters().

◆ fRelPos

std::array<QwVQWK_Channel,2> QwQPD::fRelPos

protected

◆ subelement

const TString QwQPD::subelement ={"BR","TR","BL","TL"}

staticprivate

Definition at line 27 of file QwQPD.h.

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

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ QwQPD() [1/4]

◆ QwQPD() [2/4]

◆ QwQPD() [3/4]

◆ QwQPD() [4/4]

◆ ~QwQPD()

Member Function Documentation

◆ AccumulateRunningSum() [1/2]

◆ AccumulateRunningSum() [2/2]

◆ ApplyHWChecks()

◆ ApplySingleEventCuts()

◆ CalculateRunningAverage()

◆ CheckForBurpFail()

◆ ClearEventData()

◆ ConstructBranch() [1/2]

◆ ConstructBranch() [2/2]

◆ ConstructBranchAndVector()

◆ ConstructHistograms()

◆ DeaccumulateRunningSum() [1/2]

◆ DeaccumulateRunningSum() [2/2]

◆ EncodeEventData()

◆ FillHistograms()

◆ FillTreeVector()

◆ GetAbsolutePosition()

◆ GetCalibrationFactors()

◆ GetEffectiveCharge()

◆ GetEventcutErrorFlag()

◆ GetPosition()

◆ GetSubelementByName()

◆ GetSubElementIndex()

◆ GetSubElementName()

◆ IncrementErrorCounters()

◆ InitializeChannel() [1/2]

◆ InitializeChannel() [2/2]

◆ LoadChannelParameters()

◆ MakeQPDList()

◆ operator+=() [1/2]

◆ operator+=() [2/2]

◆ operator-=() [1/2]

◆ operator-=() [2/2]

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ PrintErrorCounters()

◆ ProcessEvBuffer()

◆ ProcessEvent()

| 4 | 2 |

BL | 3 | 1 | BR

◆ RandomizeEventData()

◆ Ratio()

◆ Scale()

◆ SetDefaultSampleSize()

◆ SetEventCutMode()

◆ SetEventData()

◆ SetRandomEventParameters()

◆ SetSingleEventCuts()

◆ SetSubElementCalibrationFactor()

◆ SetSubElementPedestal()

◆ UpdateErrorFlag() [1/2]

◆ UpdateErrorFlag() [2/2]

Field Documentation

◆ fAbsPos

◆ fEffectiveCharge

◆ fPhotodiode

◆ fQPDElementList

◆ fQwQPDCalibration

◆ fRelPos

◆ subelement