QwSubsystemArrayParity Class Reference

Virtual base class for the parity subsystems. More...

#include <QwSubsystemArrayParity.h>

Inheritance diagram for QwSubsystemArrayParity:

Collaboration diagram for QwSubsystemArrayParity:

Public Member Functions
	QwSubsystemArrayParity (QwOptions &options)
	Constructor with options. More...
	QwSubsystemArrayParity (const QwSubsystemArrayParity &source)
	Copy constructor by reference. More...
virtual	~QwSubsystemArrayParity ()
	Default destructor. More...
VQwSubsystemParity *	GetSubsystemByName (const TString &name)
	Get the subsystem with the specified name. More...
void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)
	Construct a branch and vector for this subsystem with a prefix. More...
void	FillTreeVector (std::vector< Double_t > &values) const
	Fill the vector for this subsystem. More...
void	FillHistograms ()
	Fill the histograms for this subsystem. More...
void	FillDB_MPS (QwParityDB *db, TString type)
	Fill the database with MPS-based variables Note that most subsystems don't need to do this. More...
void	FillDB (QwParityDB *db, TString type)
	Fill the database. More...
void	FillErrDB (QwParityDB *db, TString type)
QwSubsystemArrayParity &	operator= (const QwSubsystemArrayParity &value)
	Assignment operator. More...
QwSubsystemArrayParity &	operator+= (const QwSubsystemArrayParity &value)
	Addition-assignment operator. More...
QwSubsystemArrayParity &	operator-= (const QwSubsystemArrayParity &value)
	Subtraction-assignment operator. More...
void	Sum (const QwSubsystemArrayParity &value1, const QwSubsystemArrayParity &value2)
	Sum of two subsystem arrays. More...
void	Difference (const QwSubsystemArrayParity &value1, const QwSubsystemArrayParity &value2)
	Difference of two subsystem arrays. More...
void	Ratio (const QwSubsystemArrayParity &numer, const QwSubsystemArrayParity &denom)
	Ratio of two subsystem arrays. More...
void	Scale (Double_t factor)
	Scale this subsystem array. More...
void	AccumulateRunningSum (const QwSubsystemArrayParity &value)
	Update the running sums for devices accumulated for the global error non-zero events/patterns. More...
void	AccumulateAllRunningSum (const QwSubsystemArrayParity &value)
	Update the running sums for devices check only the error flags at the channel level. Only used for stability checks. More...
void	DeaccumulateRunningSum (const QwSubsystemArrayParity &value)
	Remove the entry value from the running sums for devices. More...
void	CalculateRunningAverage ()
	Calculate the average for all good events. More...
void	Blind (const QwBlinder *blinder)
	Blind the asymmetry of this subsystem. More...
void	UnBlind (const QwBlinder *blinder)
	Unblind the asymmetry of this subsystem. More...
void	Blind (const QwBlinder *blinder, const QwSubsystemArrayParity &yield)
	Blind the difference of this subsystem. More...
void	UnBlind (const QwBlinder *blinder, const QwSubsystemArrayParity &yield)
	Unblind the difference of this subsystem. More...
Bool_t	ApplySingleEventCuts ()
	Apply the single event cuts. More...
void	IncrementErrorCounters ()
	Update the data elements' error counters based on their internal error flags. More...
void	PrintErrorCounters () const
	Report the number of events failed due to HW and event cut failures. More...
UInt_t	GetEventcutErrorFlag () const
	Return the error flag to the main routine. More...
void	UpdateErrorFlag ()
	Update the error flag internally from all the subsystems. More...
void	UpdateErrorFlag (const QwSubsystemArrayParity &ev_error)
	update the error flag for each channel in the subsystem array with the corresponding value in the ev_error subsystem array More...
void	PrintValue () const
	Print value of all channels. More...
void	WritePromptSummary (QwPromptSummary *ps, TString type)
virtual Bool_t	CheckForEndOfBurst () const
Public Member Functions inherited from QwSubsystemArray
	QwSubsystemArray (QwOptions &options, CanContainFn myCanContain)
	Constructor with options. More...
	QwSubsystemArray (const QwSubsystemArray &source)
	Copy constructor by reference. More...
virtual	~QwSubsystemArray ()
	Virtual destructor. More...
void	SetCodaRunNumber (UInt_t runnum)
	Set the internal record of the CODA run number. More...
void	SetCodaSegmentNumber (UInt_t segnum)
	Set the internal record of the CODA segment number. More...
void	SetCodaEventNumber (UInt_t evtnum)
	Set the internal record of the CODA event number. More...
void	SetCodaEventType (UInt_t evttype)
	Set the internal record of the CODA event type. More...
UInt_t	GetCodaRunNumber () const
	Get the internal record of the CODA run number. More...
UInt_t	GetCodaSegmentNumber () const
	Get the internal record of the CODA segment number. More...
UInt_t	GetCodaEventNumber () const
	Get the internal record of the CODA event number. More...
UInt_t	GetCodaEventType () const
	Get the internal record of the CODA event type. More...
void	SetCleanParameters (Double_t cleanparameter[3])
	Set the internal record of the CODA event number. More...
void	SetEventTypeMask (const UInt_t mask)
	Set event type mask. More...
UInt_t	GetEventTypeMask () const
	Get event type mask. More...
UInt_t	UpdateEventTypeMask ()
	Update the event type mask from the subsystems. More...
void	SetDataLoaded (const Bool_t flag)
	Set data loaded flag. More...
Bool_t	HasDataLoaded () const
	Get data loaded flag. More...
void	ProcessOptionsToplevel (QwOptions &options)
	Process configuration options for the subsystem array itself. More...
void	ProcessOptionsSubsystems (QwOptions &options)
	Process configuration options for all subsystems in the array. More...
void	ProcessOptions (QwOptions &options)
	Process configuration options (default behavior) More...
void	push_back (VQwSubsystem *subsys)
	Add the subsystem to this array. More...
virtual std::vector < VQwSubsystem * >	GetSubsystemByType (const std::string &type)
	Get the list of subsystems of the specified type. More...
void	ClearEventData ()
Int_t	ProcessConfigurationBuffer (const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
	Process the event buffer for configuration events. More...
Int_t	ProcessEvBuffer (const UInt_t event_type, const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
	Process the event buffer for events. More...
void	RandomizeEventData (int helicity=0, double time=0.0)
	Randomize the data in this event. More...
void	EncodeEventData (std::vector< UInt_t > &buffer)
	Encode the data in this event. More...
void	ProcessEvent ()
	Process the decoded data in this event. More...
void	AtEndOfEventLoop ()
	Perform actions at the end of the event loop. More...
Bool_t	RequestExternalValue (const TString &name, VQwHardwareChannel *value) const
	Retrieve the variable name from other subsystem arrays. More...
const VQwHardwareChannel *	ReturnInternalValue (const TString &name) const
	Retrieve the variable name from subsystems in this subsystem array. More...
Bool_t	ReturnInternalValue (const TString &name, VQwHardwareChannel *value) const
	Retrieve the variable name from subsystems in this subsystem array. More...
Bool_t	PublishInternalValue (const TString name, const TString desc, const VQwSubsystem *subsys, const VQwHardwareChannel *element)
	Publish the value name with description from a subsystem in this array. More...
void	ListPublishedValues () const
	List the published values and description in this subsystem array. More...
void	PrintParamFileList () const
	Print list of parameter files. More...
TList *	GetParamFileNameList (TString name) const
	Get list of parameter files. More...
void	PrintInfo () const
	Print some information about the subsystem. More...
void	push_back (boost::shared_ptr< VQwSubsystem > subsys)
void	ConstructHistograms ()
	Construct the histograms for this subsystem. More...
void	ConstructHistograms (TDirectory *folder)
	Construct the histograms for this subsystem in a folder. More...
void	ConstructHistograms (TDirectory *folder, TString &prefix)
	Construct the histograms for this subsystem in a folder with a prefix. More...
void	FillHistograms ()
	Fill the histograms for this subsystem. More...
void	ShareHistograms (const QwSubsystemArray &source)
	Share the histograms with another subsystem. More...
void	ConstructBranchAndVector (TTree *tree, std::vector< Double_t > &values)
	Construct the tree and vector for this subsystem. More...
void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)
	Construct a branch and vector for this subsystem with a prefix. More...
void	ConstructBranch (TTree *tree, TString &prefix)
	Construct a branch for this subsystem with a prefix. More...
void	ConstructBranch (TTree *tree, TString &prefix, QwParameterFile &trim_file)
	Construct a branch for this subsystem with a prefix after tree leave trimming. More...
void	FillTreeVector (std::vector< Double_t > &values) const
	Fill the vector for this subsystem. More...
void	ConstructTree ()
	Construct the tree for this subsystem. More...
void	ConstructTree (TDirectory *folder)
	Construct the tree for this subsystem in a folder. More...
void	ConstructTree (TDirectory *folder, TString &prefix)
	Construct the tree for this subsystem in a folder with a prefix. More...
void	FillTree ()
	Fill the tree for this subsystem. More...
void	DeleteTree ()
	Delete the tree for this subsystem. More...

Data Fields
const QwSubsystemArrayParity *	dummy_source
Data Fields inherited from std::vector< T >
T	elements
	STL member. More...

Static Protected Member Functions
static Bool_t	CanContain (VQwSubsystem *subsys)
	Test whether this subsystem array can contain a particular subsystem. More...
Static Protected Member Functions inherited from QwSubsystemArray
static Bool_t	CanContain (VQwSubsystem *subsys)
	Test whether this subsystem array can contain a particular subsystem. More...

Protected Attributes
UInt_t	fErrorFlag
Int_t	fErrorFlagTreeIndex
Protected Attributes inherited from QwSubsystemArray
size_t	fTreeArrayIndex
UInt_t	fCodaRunNumber
	Index of this data element in root tree. More...
UInt_t	fCodaSegmentNumber
	CODA segment number as provided by QwEventBuffer. More...
UInt_t	fCodaEventNumber
	CODA event number as provided by QwEventBuffer. More...
UInt_t	fCodaEventType
	CODA event type as provided by QwEventBuffer. More...
Double_t	fCleanParameter [3]
UInt_t	fEventTypeMask
	Mask of event types. More...
Bool_t	fHasDataLoaded
	Has this array gotten data to be processed? More...
CanContainFn	fnCanContain
	Function to determine which subsystems we can accept. More...

Private Member Functions
	QwSubsystemArrayParity ()
	Private default constructor. More...

Additional Inherited Members
Public Types inherited from QwSubsystemArray
typedef Bool_t(*	CanContainFn )(VQwSubsystem *)
Static Public Member Functions inherited from QwSubsystemArray
static void	DefineOptions (QwOptions &options)
	Define configuration options for global array. More...
Protected Member Functions inherited from QwSubsystemArray
void	LoadSubsystemsFromParameterFile (QwParameterFile &detectors)

Detailed Description

Virtual base class for the parity subsystems.

Virtual base class for the classes containing the event-based information from each parity subsystem. This will define the interfaces used in communicating with the CODA routines.

Definition at line 34 of file QwSubsystemArrayParity.h.

Constructor & Destructor Documentation

QwSubsystemArrayParity::QwSubsystemArrayParity ( )

private

Private default constructor.

QwSubsystemArrayParity::QwSubsystemArrayParity ( QwOptions &  options )

inline

Constructor with options.

Definition at line 44 of file QwSubsystemArrayParity.h.

: QwSubsystemArray(options, CanContain),fErrorFlag(0),fErrorFlagTreeIndex(-1) { };

QwSubsystemArrayParity::QwSubsystemArrayParity

(

const QwSubsystemArrayParity &

source

)

Copy constructor by reference.

Copy constructor.

Definition at line 19 of file QwSubsystemArrayParity.cc.

References dummy_source, fErrorFlag, and fErrorFlagTreeIndex.

: QwSubsystemArray(source)
{
  // Store pointer to source
  dummy_source = &source;

  // Copy error flags
  fErrorFlag = source.fErrorFlag;
  fErrorFlagTreeIndex = source.fErrorFlagTreeIndex;
}

QwSubsystemArrayParity::~QwSubsystemArrayParity ( )

virtual

Default destructor.

Destructor.

Definition at line 33 of file QwSubsystemArrayParity.cc.

{
  // nothing
}

Member Function Documentation

void QwSubsystemArrayParity::AccumulateAllRunningSum

(

const QwSubsystemArrayParity &

value

)

Update the running sums for devices check only the error flags at the channel level. Only used for stability checks.

Definition at line 316 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::AccumulateRunningSum(), QwLog::endl(), and QwError.

Referenced by QwEventRing::push().

{
  if (!value.empty()) {
    if (this->size() == value.size()) {
      //if (value.GetEventcutErrorFlag()==0){//do running sum only if error flag is zero. This way will prevent any Beam Trip(in ev mode 3) related events going into the running sum.
  for (size_t i = 0; i < value.size(); i++) {
    if (value.at(i)==NULL || this->at(i)==NULL) {
      //  Either the value or the destination subsystem
      //  are null
    } else {
      VQwSubsystemParity *ptr1 =
        dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
      if (typeid(*ptr1) == typeid(*(value.at(i).get()))) {
        ptr1->AccumulateRunningSum(value.at(i).get());
      } else {
        QwError << "QwSubsystemArrayParity::AccumulateRunningSum here where types don't match" << QwLog::endl;
        QwError << " typeid(ptr1)=" << typeid(ptr1).name()
          << " but typeid(value.at(i)))=" << typeid(value.at(i)).name()
          << QwLog::endl;
        //  Subsystems don't match
      }
    }
  }
  //}//else if ((value.fErrorFlag& 512)==512){
      //QwMessage << " AccumulateRunningSum "<<(value.fErrorFlag & 0x2FF)<<" - "<<value.GetCodaEventNumber()<< QwLog::endl;
      //}
    } else {
      //  Array sizes don't match

    }
  } else {
    //  The value is empty
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::AccumulateRunningSum

(

const QwSubsystemArrayParity &

value

)

Update the running sums for devices accumulated for the global error non-zero events/patterns.

Definition at line 282 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::AccumulateRunningSum(), QwLog::endl(), GetEventcutErrorFlag(), and QwError.

Referenced by QwHelicityPattern::AccumulateBurstSum(), QwHelicityPattern::AccumulateRunningBurstSum(), QwHelicityPattern::AccumulateRunningSum(), and main().

{
  if (!value.empty()) {
    if (this->size() == value.size()) {
      if (value.GetEventcutErrorFlag()==0){//do running sum only if error flag is zero. This way will prevent any Beam Trip(in ev mode 3) related events going into the running sum.
  for (size_t i = 0; i < value.size(); i++) {
    if (value.at(i)==NULL || this->at(i)==NULL) {
      //  Either the value or the destination subsystem
      //  are null
    } else {
      VQwSubsystemParity *ptr1 =
        dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
      if (typeid(*ptr1) == typeid(*(value.at(i).get()))) {
        ptr1->AccumulateRunningSum(value.at(i).get());
      } else {
        QwError << "QwSubsystemArrayParity::AccumulateRunningSum here where types don't match" << QwLog::endl;
        QwError << " typeid(ptr1)=" << typeid(ptr1).name()
          << " but typeid(value.at(i)))=" << typeid(value.at(i)).name()
          << QwLog::endl;
        //  Subsystems don't match
      }
    }
  }
      }

    } else {
      //  Array sizes don't match

    }
  } else {
    //  The value is empty
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

Bool_t QwSubsystemArrayParity::ApplySingleEventCuts

(

)

Apply the single event cuts.

Definition at line 478 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::ApplySingleEventCuts(), fErrorFlag, VQwSubsystemParity::GetEventcutErrorFlag(), kEventCutMode3, kGlobalCut, and UpdateErrorFlag().

Referenced by main().

                                                   {
  Int_t CountFalse;
  Bool_t status;
  UInt_t ErrorFlag;
  fErrorFlag=0;
  VQwSubsystemParity *subsys_parity;
  CountFalse=0;
  if (!empty()){
    for (iterator subsys = begin(); subsys != end(); ++subsys){
      subsys_parity=dynamic_cast<VQwSubsystemParity*>((subsys)->get());
      status=subsys_parity->ApplySingleEventCuts();
      ErrorFlag = subsys_parity->GetEventcutErrorFlag();
      if ((ErrorFlag & kEventCutMode3)==kEventCutMode3)//we only care about the event cut flag in event cut mode 3
  fErrorFlag |= ErrorFlag; 
      if (!status)
      {
  if ((ErrorFlag&kGlobalCut)==kGlobalCut){
    CountFalse++;
    fErrorFlag |= ErrorFlag; //we need the error code for failed events in event mode 2 for beam trips and etc.
  }
      }
      

    }
  }
  if (CountFalse > 0)
    status = kFALSE;
  else
    status = kTRUE;

  //  Propagate all error codes to derived objects in the subsystems.
  UpdateErrorFlag();

  return status;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::Blind

(

const QwBlinder *

blinder

)

Blind the asymmetry of this subsystem.

Definition at line 392 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::Blind(), QwLog::endl(), and QwError.

Referenced by QwBlinder::Blind().

{
  // Loop over subsystem array
  for (size_t i = 0; i < this->size(); i++) {
    // Cast into parity subsystems
    VQwSubsystemParity* subsys = dynamic_cast<VQwSubsystemParity*>(this->at(i).get());

    // Check for null pointers
    if (this->at(i) == 0) {
      QwError << "QwSubsystemArrayParity::Blind: "
              << "parity subsystem null pointer!" << QwLog::endl;
      return;
    }

    // Apply blinding
    subsys->Blind(blinder);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Blind the difference of this subsystem.

Definition at line 411 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::Blind(), QwLog::endl(), and QwError.

{
  // Check for array size
  if (this->size() != yield.size()) {
    QwError << "QwSubsystemArrayParity::Blind: "
            << "diff and yield array dimension mismatch!" << QwLog::endl;
    return;
  }

  // Loop over subsystem array
  for (size_t i = 0; i < this->size(); i++) {
    // Cast into parity subsystems
    VQwSubsystemParity* subsys_diff  = dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
    VQwSubsystemParity* subsys_yield = dynamic_cast<VQwSubsystemParity*>(yield.at(i).get());

    // Check for null pointers
    if (subsys_diff == 0 || subsys_yield == 0) {
      QwError << "QwSubsystemArrayParity::Blind: "
              << "diff or yield parity subsystem null pointer!" << QwLog::endl;
      return;
    }

    // Apply blinding
    subsys_diff->Blind(blinder, subsys_yield);
  }
}

Here is the call graph for this function:

void QwSubsystemArrayParity::CalculateRunningAverage

(

)

Calculate the average for all good events.

Definition at line 272 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::CalculateRunningAverage().

Referenced by QwHelicityPattern::CalculateBurstAverage(), QwHelicityPattern::CalculateRunningAverage(), QwHelicityPattern::CalculateRunningBurstAverage(), main(), and QwEventRing::push().

{
  for (iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    subsys_parity->CalculateRunningAverage();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

static Bool_t QwSubsystemArrayParity::CanContain ( VQwSubsystem *  subsys )

inlinestaticprotected

Test whether this subsystem array can contain a particular subsystem.

Definition at line 135 of file QwSubsystemArrayParity.h.

                                                   {
      return (dynamic_cast<VQwSubsystemParity*>(subsys) != 0);
    };

Bool_t QwSubsystemArrayParity::CheckForEndOfBurst ( ) const

virtual

Definition at line 260 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::CheckForEndOfBurst().

Referenced by main().

{
  Bool_t status = kFALSE;
  for (const_iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    status |= subsys_parity->CheckForEndOfBurst();
  }
  return status;
};

Here is the call graph for this function:

Here is the caller graph for this function:

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

Construct a branch and vector for this subsystem with a prefix.

Definition at line 590 of file QwSubsystemArrayParity.cc.

References QwSubsystemArray::ConstructBranchAndVector(), and fErrorFlagTreeIndex.

Referenced by QwHelicityPattern::ConstructBranchAndVector(), and main().

                                                                                                               {
  QwSubsystemArray::ConstructBranchAndVector(tree, prefix, values);
  if (prefix=="yield_" || prefix==""){
    values.push_back(0.0);
    fErrorFlagTreeIndex = values.size()-1;
    tree->Branch("ErrorFlag",&(values[fErrorFlagTreeIndex]),"ErrorFlag/D");
  } else {
    fErrorFlagTreeIndex = -1;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::DeaccumulateRunningSum

(

const QwSubsystemArrayParity &

value

)

Remove the entry value from the running sums for devices.

Definition at line 352 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::DeaccumulateRunningSum(), QwLog::endl(), and QwError.

Referenced by QwEventRing::pop().

{
  //Bool_t berror=kTRUE;//only needed for deaccumulation (stability check purposes)
  //if (value.fErrorFlag>0){//check the error is global
  //berror=((value.fErrorFlag & 0x2FF) == 0); //The operation value.fErrorFlag & 0x2FF clear everything else but the HW errors + event cut errors + blinder error    
  //}
  if (!value.empty()) {
    if (this->size() == value.size()) {
      //if (value.GetEventcutErrorFlag()==0){//do derunningsum only if error flag is zero. 
  for (size_t i = 0; i < value.size(); i++) {
    if (value.at(i)==NULL || this->at(i)==NULL) {
      //  Either the value or the destination subsystem
      //  are null
    } else {
      VQwSubsystemParity *ptr1 =
        dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
      if (typeid(*ptr1) == typeid(*(value.at(i).get()))) {
        ptr1->DeaccumulateRunningSum(value.at(i).get());
      } else {
        QwError << "QwSubsystemArrayParity::AccumulateRunningSum here where types don't match" << QwLog::endl;
        QwError << " typeid(ptr1)=" << typeid(ptr1).name()
          << " but typeid(value.at(i)))=" << typeid(value.at(i)).name()
          << QwLog::endl;
        //  Subsystems don't match
      }
    }
  }
  //}//else if ((value.fErrorFlag & 268435968)==268435968){
      //QwMessage << " DeaccumulateRunningSum "<<(value.fErrorFlag & 0x2FF)<<" - "<<value.GetCodaEventNumber()<< QwLog::endl;
      //}
    } else {
      //  Array sizes don't match

    }
  } else {
    //  The value is empty
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Difference of two subsystem arrays.

Difference of two subsystem arrays

Parameters

value1	First subsystem array
value2	Second subsystem array

Definition at line 221 of file QwSubsystemArrayParity.cc.

Referenced by QwHelicityPattern::CalculateAsymmetry().

{

  if (!value1.empty()&& !value2.empty()){
    *(this)  =  value1;
    *(this) -= value2;
  } else {
    //  The source is empty
  }
}

Here is the caller graph for this function:

void QwSubsystemArrayParity::FillDB	(	QwParityDB *	db,
		TString	type
	)

Fill the database.

Definition at line 55 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::FillDB().

Referenced by QwHelicityPattern::FillDB().

{
  for (iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    subsys_parity->FillDB(db, type);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::FillDB_MPS	(	QwParityDB *	db,
		TString	type
	)

Fill the database with MPS-based variables Note that most subsystems don't need to do this.

Definition at line 47 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::FillDB_MPS().

Referenced by main().

{
  for (iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    subsys_parity->FillDB_MPS(db, type);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::FillErrDB	(	QwParityDB *	db,
		TString	type
	)

Definition at line 63 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::FillErrDB().

Referenced by QwHelicityPattern::FillErrDB().

{
  //  for (const_iterator subsys = dummy_source->begin(); subsys != dummy_source->end(); ++subsys) {
  for (iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    subsys_parity->FillErrDB(db, type);
  }
  return;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::FillHistograms

(

)

Fill the histograms for this subsystem.

Definition at line 611 of file QwSubsystemArrayParity.cc.

References QwSubsystemArray::FillHistograms(), and GetEventcutErrorFlag().

Referenced by QwHelicityPattern::FillHistograms(), and main().

{
  if (GetEventcutErrorFlag()==0)
    QwSubsystemArray::FillHistograms();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::FillTreeVector

(

std::vector< Double_t > &

values

)

const

Fill the vector for this subsystem.

Definition at line 601 of file QwSubsystemArrayParity.cc.

References fErrorFlag, fErrorFlagTreeIndex, and QwSubsystemArray::FillTreeVector().

Referenced by QwHelicityPattern::FillTreeVector(), and main().

{
  QwSubsystemArray::FillTreeVector(values);
  if (fErrorFlagTreeIndex>=0 && fErrorFlagTreeIndex<static_cast<int>(values.size())){
    values.at(fErrorFlagTreeIndex) = fErrorFlag;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

UInt_t QwSubsystemArrayParity::GetEventcutErrorFlag ( ) const

inline

Return the error flag to the main routine.

Definition at line 114 of file QwSubsystemArrayParity.h.

References fErrorFlag.

Referenced by AccumulateRunningSum(), FillHistograms(), QwHelicityPattern::GetEventcutErrorFlag(), and QwRegression::LinearRegression().

                                       {
      return fErrorFlag;
    };

Here is the caller graph for this function:

VQwSubsystemParity * QwSubsystemArrayParity::GetSubsystemByName ( const TString &  name )

virtual

Get the subsystem with the specified name.

Reimplemented from QwSubsystemArray.

Definition at line 40 of file QwSubsystemArrayParity.cc.

References QwSubsystemArray::GetSubsystemByName().

Referenced by main().

{
  return dynamic_cast<VQwSubsystemParity*>(QwSubsystemArray::GetSubsystemByName(name));
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::IncrementErrorCounters

(

)

Update the data elements' error counters based on their internal error flags.

Definition at line 516 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::IncrementErrorCounters().

Referenced by QwHelicityPattern::CalculateAsymmetry(), and main().

{
  VQwSubsystemParity *subsys_parity;
  if (!empty()){
    for (iterator subsys = begin(); subsys != end(); ++subsys){
      subsys_parity=dynamic_cast<VQwSubsystemParity*>((subsys)->get());
      subsys_parity->IncrementErrorCounters();
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

QwSubsystemArrayParity & QwSubsystemArrayParity::operator+=

(

const QwSubsystemArrayParity &

value

)

Addition-assignment operator.

Addition-assignment operator

Parameters

value

Subsystem array to add to this array

Returns

This subsystem array after addition

Definition at line 130 of file QwSubsystemArrayParity.cc.

References QwLog::endl(), fErrorFlag, and QwError.

{
  if (!value.empty()){

    if (this->size() == value.size()){
      this->fErrorFlag|=value.fErrorFlag;
      for(size_t i=0;i<value.size();i++){
  if (value.at(i)==NULL || this->at(i)==NULL){
    //  Either the value or the destination subsystem
    //  are null
  } else {
    VQwSubsystemParity *ptr1 =
      dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
    if (typeid(*ptr1)==typeid(*(value.at(i).get()))){
      *(ptr1) += value.at(i).get();
      //std::cout<<"QwSubsystemArrayParity::operator+ here where types match \n";
    } else {
      QwError << "QwSubsystemArrayParity::operator+ here where types don't match" << QwLog::endl;
      QwError << " typeid(ptr1)=" << typeid(ptr1).name()
                    << " but typeid(value.at(i)))=" << typeid(value.at(i)).name()
                    << QwLog::endl;
      //  Subsystems don't match
    }
  }
      }
    } else {
      //  Array sizes don't match
    }
  } else {
    //  The vsource is empty
  }
  return *this;
}

Here is the call graph for this function:

QwSubsystemArrayParity & QwSubsystemArrayParity::operator-=

(

const QwSubsystemArrayParity &

value

)

Subtraction-assignment operator.

Subtraction-assignment operator

Parameters

value

Subsystem array to subtract from this array

Returns

This array after subtraction

Definition at line 169 of file QwSubsystemArrayParity.cc.

References fErrorFlag.

{
  if (!value.empty()){

    if (this->size() == value.size()){
      this->fErrorFlag|=value.fErrorFlag;
      for(size_t i=0;i<value.size();i++){
  if (value.at(i)==NULL || this->at(i)==NULL){
    //  Either the value or the destination subsystem
    //  are null
  } else {
    VQwSubsystemParity *ptr1 =
      dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
    if (typeid(*ptr1)==typeid(*(value.at(i).get()))){
      *(ptr1) -= value.at(i).get();
    } else {
      //  Subsystems don't match
    }
  }
      }
    } else {
      //  Array sizes don't match
    }
  } else {
    //  The vsource is empty
  }
  return *this;
}

QwSubsystemArrayParity & QwSubsystemArrayParity::operator=

(

const QwSubsystemArrayParity &

source

)

Assignment operator.

Assignment operator

Parameters

source

Subsystem array to assign to this array

Returns

This subsystem array after assignment

Definition at line 88 of file QwSubsystemArrayParity.cc.

References QwLog::endl(), QwSubsystemArray::fCodaEventNumber, fErrorFlag, and QwError.

{
  Bool_t localdebug=kFALSE;
  if(localdebug)  std::cout<<"QwSubsystemArrayParity::operator= \n";
  if (!source.empty()){
    if (this->size() == source.size()){
      this->fErrorFlag=source.fErrorFlag;
      this->fCodaEventNumber=source.fCodaEventNumber;
      for(size_t i=0;i<source.size();i++){
  if (source.at(i)==NULL || this->at(i)==NULL){
    //  Either the source or the destination subsystem
    //  are null
  } else {
    VQwSubsystemParity *ptr1 =
      dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
    if (typeid(*ptr1)==typeid(*(source.at(i).get()))){
      if(localdebug) std::cout<<" here in QwSubsystemArrayParity::operator= types mach \n";
      *(ptr1) = source.at(i).get();
    } else {
      //  Subsystems don't match
        QwError << " QwSubsystemArrayParity::operator= types do not mach" << QwLog::endl;
        QwError << " typeid(ptr1)=" << typeid(*ptr1).name()
                      << " but typeid(*(source.at(i).get()))=" << typeid(*(source.at(i).get())).name()
                      << QwLog::endl;
    }
  }
      }
    } else {
      //  Array sizes don't match
    }
  } else {
    //  The source is empty
  }
  return *this;
}

Here is the call graph for this function:

void QwSubsystemArrayParity::PrintErrorCounters

(

)

const

Report the number of events failed due to HW and event cut failures.

Definition at line 528 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::PrintErrorCounters().

Referenced by main().

                                                     {// report number of events failed due to HW and event cut faliure
  const VQwSubsystemParity *subsys_parity;
  if (!empty()){
    for (const_iterator subsys = begin(); subsys != end(); ++subsys){
      subsys_parity=dynamic_cast<const VQwSubsystemParity*>((subsys)->get());
      subsys_parity->PrintErrorCounters();
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::PrintValue

(

)

const

Print value of all channels.

Definition at line 248 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::PrintValue().

Referenced by main(), QwHelicityPattern::PrintBurstAverage(), QwHelicityPattern::PrintRunningAverage(), and QwHelicityPattern::PrintRunningBurstAverage().

{
  for (const_iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    subsys_parity->PrintValue();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Ratio of two subsystem arrays.

Definition at line 438 of file QwSubsystemArrayParity.cc.

References QwLog::endl(), fErrorFlag, QwError, and VQwSubsystemParity::Ratio().

Referenced by QwHelicityPattern::AccumulateRunningBurstSum(), and QwHelicityPattern::CalculateAsymmetry().

{
  Bool_t localdebug=kFALSE;

  if(localdebug) std::cout<<"QwSubsystemArrayParity::Ratio \n";
  *this=numer;
  if ( !denom.empty()){
    this->fErrorFlag=(numer.fErrorFlag|denom.fErrorFlag);
    if (this->size() == denom.size() ){
      for(size_t i=0;i<denom.size();i++){
        if (denom.at(i)==NULL || this->at(i)==NULL){
          //  Either the value or the destination subsystem  are null
    if(localdebug) std::cout<<"Either the value or the destination subsystem  are null\n";
        } else {
    VQwSubsystemParity *ptr1 =
      dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
    if (typeid(*ptr1)==typeid(*(denom.at(i).get())))
            {
              ptr1->Ratio(numer.at(i).get(),denom.at(i).get());
            } else {
              //  Subsystems don't match
              QwError << "subsystem #" << i
                      << " type do not match : ratio computation aborted" << QwLog::endl;
            }
        }
      }
    } else {
      QwError << "array size do not match : ratio computation aborted" << QwLog::endl;
      //  Array sizes don't match
    }
  } else {
    QwError << "source empty : ratio computation aborted" << QwLog::endl;
    //  The source is empty
  }
  if(localdebug) std::cout<<"I am out of it \n";

}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::Scale

(

Double_t

factor

)

Scale this subsystem array.

Scale this subsystem array

Parameters

factor

Scale factor

Definition at line 238 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::Scale().

Referenced by QwHelicityPattern::CalculateAsymmetry().

{
  for (iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    subsys_parity->Scale(factor);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Sum of two subsystem arrays.

Sum of two subsystem arrays

Parameters

value1	First subsystem array
value2	Second subsystem array

Definition at line 203 of file QwSubsystemArrayParity.cc.

Referenced by QwHelicityPattern::CalculateAsymmetry().

{

  if (!value1.empty()&& !value2.empty()){
    *(this)  = value1;
    *(this) += value2;
  } else {
    //  The source is empty
  }
}

Here is the caller graph for this function:

void QwSubsystemArrayParity::UnBlind ( const QwBlinder *  blinder )

inline

Unblind the asymmetry of this subsystem.

Definition at line 96 of file QwSubsystemArrayParity.h.

Referenced by QwBlinder::UnBlind().

      { /* Not yet implemented */ };

Here is the caller graph for this function:

void QwSubsystemArrayParity::UnBlind ( const QwBlinder *  blinder, const QwSubsystemArrayParity &  yield  )

inline

Unblind the difference of this subsystem.

Definition at line 101 of file QwSubsystemArrayParity.h.

      { /* Not yet implemented */ };

void QwSubsystemArrayParity::UpdateErrorFlag

(

)

Update the error flag internally from all the subsystems.

Definition at line 574 of file QwSubsystemArrayParity.cc.

References fErrorFlag, and VQwSubsystemParity::UpdateErrorFlag().

Referenced by ApplySingleEventCuts(), QwHelicityPattern::CalculateAsymmetry(), and QwEventRing::push().

                                             { 
  //this routine will refresh the global error flag after stability cut check
  //by default at the ApplySingleEventCuts routine fErrorFlag is updated properly and a const GetEventcutErrorFlag() routine 
  //returns the fErrorFlag value
  fErrorFlag=0;
  
  VQwSubsystemParity *subsys_parity;
  if (!empty()){
    for (iterator subsys = begin(); subsys != end(); ++subsys){
      subsys_parity=dynamic_cast<VQwSubsystemParity*>((subsys)->get());
      //Update the error flag of the parity subsystem
      fErrorFlag|=subsys_parity->UpdateErrorFlag();
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSubsystemArrayParity::UpdateErrorFlag

(

const QwSubsystemArrayParity &

ev_error

)

update the error flag for each channel in the subsystem array with the corresponding value in the ev_error subsystem array

Definition at line 538 of file QwSubsystemArrayParity.cc.

References QwLog::endl(), fErrorFlag, QwError, and VQwSubsystemParity::UpdateErrorFlag().

                                                                                  {
  Bool_t localdebug=kFALSE;//kTRUE;
  if(localdebug)  std::cout<<"QwSubsystemArrayParity::UpdateErrorFlag \n";
  if (!ev_error.empty()){
    if (this->size() == ev_error.size()){
      this->fErrorFlag |= ev_error.fErrorFlag;
      for(size_t i=0;i<ev_error.size();i++){
  if (ev_error.at(i)==NULL || this->at(i)==NULL){
    //  Either the source or the destination subsystem
    //  are null
  } else {
    VQwSubsystemParity *ptr1 =
      dynamic_cast<VQwSubsystemParity*>(this->at(i).get());
    if (typeid(*ptr1)==typeid(*(ev_error.at(i).get()))){
      if(localdebug) std::cout<<" here in QwSubsystemArrayParity::UpdateErrorFlag types mach \n";
      //*(ptr1) = ev_error.at(i).get();//when =operator is used
      //pass the correct subsystem to update the errorflags at subsystem to devices to channel levels
      ptr1->UpdateErrorFlag(ev_error.at(i).get());
    } else {
      //  Subsystems don't match
        QwError << " QwSubsystemArrayParity::UpdateErrorFlag types do not mach" << QwLog::endl;
        QwError << " typeid(ptr1)=" << typeid(*ptr1).name()
                      << " but typeid(*(ev_error.at(i).get()))=" << typeid(*(ev_error.at(i).get())).name()
                      << QwLog::endl;
    }
  }
      }
    } else {
      //  Array sizes don't match
    }
  } else {
    //  The source is empty
  }  
};

Here is the call graph for this function:

void QwSubsystemArrayParity::WritePromptSummary	(	QwPromptSummary *	ps,
		TString	type
	)

Definition at line 73 of file QwSubsystemArrayParity.cc.

References VQwSubsystemParity::WritePromptSummary().

Referenced by QwHelicityPattern::WritePromptSummary().

{
  for (const_iterator subsys = begin(); subsys != end(); ++subsys) {
    VQwSubsystemParity* subsys_parity = dynamic_cast<VQwSubsystemParity*>(subsys->get());
    subsys_parity->WritePromptSummary(ps, type);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

Field Documentation

const QwSubsystemArrayParity* QwSubsystemArrayParity::dummy_source

Definition at line 65 of file QwSubsystemArrayParity.h.

Referenced by QwSubsystemArrayParity().

UInt_t QwSubsystemArrayParity::fErrorFlag

protected

Definition at line 137 of file QwSubsystemArrayParity.h.

Referenced by ApplySingleEventCuts(), FillTreeVector(), GetEventcutErrorFlag(), operator+=(), operator-=(), operator=(), QwSubsystemArrayParity(), Ratio(), and UpdateErrorFlag().

Int_t QwSubsystemArrayParity::fErrorFlagTreeIndex

protected

Definition at line 140 of file QwSubsystemArrayParity.h.

Referenced by ConstructBranchAndVector(), FillTreeVector(), and QwSubsystemArrayParity().

---

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

• QwSubsystemArrayParity.h
• QwSubsystemArrayParity.cc