Data handler computing correlations and linear-regression coefficients. More...

#include <QwCorrelatorNew.h>

Inheritance diagram for QwCorrelatorNew:

Collaboration diagram for QwCorrelatorNew:

Public Member Functions
	QwCorrelatorNew (const TString &name)
	Constructor with name.

	QwCorrelatorNew (const QwCorrelatorNew &name)

	~QwCorrelatorNew () override

void	ProcessOptions (QwOptions &options)

void	ParseConfigFile (QwParameterFile &file) override

Int_t	LoadChannelMap (const std::string &mapfile) override

Int_t	ConnectChannels (QwSubsystemArrayParity &asym, QwSubsystemArrayParity &diff) override
	Connect to Channels (asymmetry/difference only)

void	ProcessData () override

void	FinishDataHandler () override

void	CalcCorrelations ()

void	ConstructTreeBranches (QwRootFile *treerootfile, const std::string &treeprefix="", const std::string &branchprefix="") override
	Construct the tree branches.

void	FillTreeBranches (QwRootFile *treerootfile) override
	Fill the tree branches.

void	ConstructHistograms (TDirectory *folder, TString &prefix) override
	Construct the histograms in a folder with a prefix.

void	FillHistograms () override
	Fill the histograms.

void	ClearEventData () override

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

void	solve ()

bool	failed ()

void	printSummaryP () const

void	printSummaryY () const

void	printSummaryYP () const

void	printSummaryAlphas () const

void	printSummaryMeansWithUnc () const

void	printSummaryMeansWithUncCorrected () const

void	print ()

void	init ()

void	clear ()

void	setDims (int a, int b)

Int_t	getMeanP (const int i, Double_t &mean) const
	Get mean value of a variable, returns error code.

Int_t	getMeanY (const int i, Double_t &mean) const

Int_t	getMeanYprime (const int i, Double_t &mean) const

Int_t	getSigmaP (const int i, Double_t &sigma) const
	Get mean value of a variable, returns error code.

Int_t	getSigmaY (const int i, Double_t &sigma) const

Int_t	getSigmaYprime (const int i, Double_t &sigma) const

Int_t	getCovarianceP (int i, int j, Double_t &covar) const
	Get mean value of a variable, returns error code.

Int_t	getCovariancePY (int ip, int iy, Double_t &covar) const

Int_t	getCovarianceY (int i, int j, Double_t &covar) const

double	getUsedEve () const

QwCorrelatorNew &	operator+= (const std::pair< TVectorD, TVectorD > &rhs)

QwCorrelatorNew &	operator+= (const QwCorrelatorNew &rhs)

Public Member Functions inherited from VQwDataHandler
	VQwDataHandler (const TString &name)

	VQwDataHandler (const VQwDataHandler &source)

void	SetPointer (QwHelicityPattern *ptr)

void	SetPointer (QwSubsystemArrayParity *ptr)

virtual Int_t	ConnectChannels (QwSubsystemArrayParity &, QwSubsystemArrayParity &asym, QwSubsystemArrayParity &diff)

virtual Int_t	ConnectChannels (QwSubsystemArrayParity &)

Int_t	ConnectChannels (QwHelicityPattern &helicitypattern)

virtual void	UpdateBurstCounter (Short_t burstcounter)

	~VQwDataHandler () override

TString	GetName ()

void	InitRunningSum ()

void	AccumulateRunningSum ()

void	CalculateRunningAverage ()

void	PrintValue () const

void	WritePromptSummary (QwPromptSummary *ps, TString type)

virtual void	ConstructNTupleFields (QwRootFile *treerootfile, const std::string &treeprefix="", const std::string &branchprefix="")
	RNTuple methods.

virtual void	FillNTupleFields (QwRootFile *treerootfile)

void	FillTreeVector (QwRootTreeBranchVector &values) const

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

void	SetRunLabel (TString x)

Int_t	LoadChannelMap ()

Bool_t	PublishInternalValues () const override
	Publish all variables of the subsystem.

Bool_t	PublishByRequest (TString device_name) override
	Try to publish an internal variable matching the submitted name.

Public Member Functions inherited from MQwPublishable_child< QwDataHandlerArray, VQwDataHandler >
	MQwPublishable_child ()
	Default constructor Initializes the child object and sets up self-reference for publishing.

	MQwPublishable_child (const MQwPublishable_child &source)
	Copy constructor.

virtual	~MQwPublishable_child ()
	Virtual destructor.

void	SetParent (QwDataHandlerArray *parent)
	Set the parent container for this child object.

QwDataHandlerArray *	GetParent () const
	Get the parent container for this child object.

Public Member Functions inherited from MQwCloneable< VQwDataHandler, QwCorrelatorNew >
	~MQwCloneable () override
	Virtual destructor.

VQwDataHandler *	Clone () const override
	Concrete clone method.

const VQwFactory< VQwDataHandler > *	Factory () const override
	Factory getter.

Public Member Functions inherited from VQwCloneable< VQwDataHandler >
virtual	~VQwCloneable ()
	Virtual destructor.

std::string	GetClassName () const
	Get demangled name of this class.

Static Public Member Functions
static void	DefineOptions (QwOptions &options)

Static Public Member Functions inherited from MQwCloneable< VQwDataHandler, QwCorrelatorNew >
static VQwDataHandler *	Create (const std::string &name)
	Object creation.

static QwCorrelatorNew *	Cast (QwCorrelatorNew *type)
	Object dynamic cast.

Protected Member Functions
void	OpenAlphaFile (const std::string &prefix)

void	WriteAlphaFile ()

void	CloseAlphaFile ()

void	OpenAliasFile (const std::string &prefix)

void	WriteAliasFile ()

void	CloseAliasFile ()

Protected Member Functions inherited from VQwDataHandler
	VQwDataHandler ()

void	SetEventcutErrorFlagPointer (const UInt_t *errorflagptr)

UInt_t	GetEventcutErrorFlag () const

std::pair< EQwHandleType, std::string >	ParseHandledVariable (const std::string &variable)

void	CalcOneOutput (const VQwHardwareChannel dv, VQwHardwareChannel output, std::vector< const VQwHardwareChannel * > &ivs, std::vector< Double_t > &sens)

Protected Member Functions inherited from MQwPublishable_child< QwDataHandlerArray, VQwDataHandler >
Bool_t	RequestExternalValue (const TString &name, VQwHardwareChannel *value) const
	Retrieve the variable name from other subsystem arrays Get the value corresponding to some variable name from a different data array.

const VQwHardwareChannel *	RequestExternalPointer (const TString &name) const
	Retrieve a pointer to an external variable by name Requests a direct pointer to a variable from sibling subsystems via the parent container.

Bool_t	PublishInternalValue (const TString name, const TString desc, const VQwHardwareChannel *element) const
	Publish a variable from this child into the parent container.

Protected Attributes
Int_t	fBlock

bool	fDisableHistos

std::vector< std::string >	fIndependentFull

std::vector< EQwHandleType >	fIndependentType

std::vector< std::string >	fIndependentName

std::vector< const VQwHardwareChannel * >	fIndependentVar

std::vector< Double_t >	fIndependentValues

std::string	fAlphaOutputFileBase

std::string	fAlphaOutputFileSuff

std::string	fAlphaOutputPath

TFile *	fAlphaOutputFile

TTree *	fTree

std::string	fAliasOutputFileBase

std::string	fAliasOutputFileSuff

std::string	fAliasOutputPath

std::ofstream	fAliasOutputFile

int	fTotalCount

int	fGoodCount

int	fErrCounts_EF

std::vector< int >	fErrCounts_IV

std::vector< int >	fErrCounts_DV

unsigned int	fGoodEvent

Protected Attributes inherited from VQwDataHandler
Int_t	fPriority

Short_t	fBurstCounter
	When a datahandler array is processed, handlers with lower priority will be processed before handlers with higher priority.

TString	fName

std::string	fMapFile

std::string	fTreeName

std::string	fTreeComment

std::string	fPrefix

TString	run_label

const UInt_t *	fErrorFlagPtr
	Error flag pointer.

QwSubsystemArrayParity *	fSubsystemArray
	Single event pointer.

QwHelicityPattern *	fHelicityPattern
	Helicity pattern pointer.

std::vector< std::string >	fDependentFull

std::vector< EQwHandleType >	fDependentType

std::vector< std::string >	fDependentName

std::vector< const VQwHardwareChannel * >	fDependentVar

std::vector< Double_t >	fDependentValues

std::vector< VQwHardwareChannel * >	fOutputVar

std::vector< Double_t >	fOutputValues

std::vector< std::vector< TString > >	fPublishList

std::string	ParseSeparator

Bool_t	fKeepRunningSum

Bool_t	fRunningsumFillsTree

VQwDataHandler *	fRunningsum

Private Member Functions
	QwCorrelatorNew ()

Private Attributes
TString	fNameNoSpaces

int	nP

int	nY

std::vector< TH1D >	fHnames

std::vector< TH1D >	fH1iv

std::vector< TH1D >	fH1dv

std::vector< std::vector< TH2D > >	fH2iv

std::vector< std::vector< TH2D > >	fH2dv

Int_t	fCycleCounter

Int_t	fErrorFlag
	is information valid

Long64_t	fGoodEventNumber
	accumulated so far

TMatrixD	mRPY
	correlations

TMatrixD	mRYP

TMatrixD	mRPP

TMatrixD	mRYY

TMatrixD	mRYYp

TMatrixD	mVPY
	unnormalized covariances

TMatrixD	mVYP

TMatrixD	mVPP

TMatrixD	mVYY

TMatrixD	mVYYp

TVectorD	mVP
	variances

TVectorD	mVY

TVectorD	mVYp

TMatrixD	mSPY
	normalized covariances

TMatrixD	mSYP

TMatrixD	mSPP

TMatrixD	mSYY

TMatrixD	mSYYp

TVectorD	mSP
	sigmas

TVectorD	mSY

TVectorD	mSYp

TVectorD	mMP
	mean values

TVectorD	mMY

TVectorD	mMYp

TMatrixD	Axy
	slopes

TMatrixD	Ayx

TMatrixD	dAxy

TMatrixD	dAyx

Static Private Attributes
static bool	fPrintCorrelations = false

Friends
QwCorrelatorNew	operator+ (QwCorrelatorNew lhs, const QwCorrelatorNew &rhs)

std::ostream &	operator<< (std::ostream &stream, const QwCorrelatorNew &h)
	Output stream operator.

Additional Inherited Members
Public Types inherited from VQwDataHandler
enum	EQwHandleType { kHandleTypeUnknown =0 , kHandleTypeMps , kHandleTypeAsym , kHandleTypeDiff , kHandleTypeYield }

typedef std::vector< VQwHardwareChannel * >::iterator	Iterator_HdwChan

typedef std::vector< VQwHardwareChannel * >::const_iterator	ConstIterator_HdwChan

Detailed Description

Data handler computing correlations and linear-regression coefficients.

Uses Bevington/Pebay algorithms to estimate correlations between independent and dependent variables selected from subsystem arrays. Produces summary histograms and optional output trees/files for further analysis.

Definition at line 26 of file QwCorrelatorNew.h.

Constructor & Destructor Documentation

◆ QwCorrelatorNew() [1/3]

QwCorrelatorNew::QwCorrelatorNew ( const TString & name )

Constructor with name.

Definition at line 21 of file QwCorrelatorNew.cc.

: VQwDataHandler(name),
  fBlock(-1),
  fDisableHistos(true),
  fAlphaOutputFileBase("blueR"),
  fAlphaOutputFileSuff("new.slope.root"),
  fAlphaOutputPath("."),
  fAlphaOutputFile(0),
  fTree(0),
  fAliasOutputFileBase("regalias_"),
  fAliasOutputFileSuff(""),
  fAliasOutputPath("."),
  fNameNoSpaces(name),
  nP(0),nY(0),
  fCycleCounter(0)
{
  fNameNoSpaces.ReplaceAll(" ","_");
  // Set default tree name and descriptions (in VQwDataHandler)
  fTreeName = "lrb";
  fTreeComment = "Correlations";
  // Parsing separator
  ParseSeparator = "_";
 
  // Clear all data
  ClearEventData();
}

References ClearEventData(), fAliasOutputFileBase, fAliasOutputFileSuff, fAliasOutputPath, fAlphaOutputFile, fAlphaOutputFileBase, fAlphaOutputFileSuff, fAlphaOutputPath, fBlock, fCycleCounter, fDisableHistos, fNameNoSpaces, fTree, VQwDataHandler::fTreeComment, VQwDataHandler::fTreeName, nP, nY, VQwDataHandler::ParseSeparator, and VQwDataHandler::VQwDataHandler().

Referenced by AccumulateRunningSum(), operator+, operator+=(), operator+=(), operator<<, and QwCorrelatorNew().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ QwCorrelatorNew() [2/3]

QwCorrelatorNew::QwCorrelatorNew ( const QwCorrelatorNew & name )

Definition at line 709 of file QwCorrelatorNew.cc.

: VQwDataHandler(source),
  fBlock(source.fBlock),
  fDisableHistos(source.fDisableHistos),
  fAlphaOutputFileBase(source.fAlphaOutputFileBase),
  fAlphaOutputFileSuff(source.fAlphaOutputFileSuff),
  fAlphaOutputPath(source.fAlphaOutputPath),
  fAlphaOutputFile(0),
  fTree(0),
  fAliasOutputFileBase(source.fAliasOutputFileBase),
  fAliasOutputFileSuff(source.fAliasOutputFileSuff),
  fAliasOutputPath(source.fAliasOutputPath),
  nP(source.nP),nY(source.nY),
  fCycleCounter(source.fCycleCounter),
  fErrorFlag(-1),
  fGoodEventNumber(0)
{
  QwMessage << fGoodEventNumber << QwLog::endl;
 
  QwWarning << "QwCorrelatorNew copy constructor required but untested" << QwLog::endl;
 
  // Clear all data
  ClearEventData();
}

References ClearEventData(), QwLog::endl(), fAliasOutputFileBase, fAliasOutputFileSuff, fAliasOutputPath, fAlphaOutputFile, fAlphaOutputFileBase, fAlphaOutputFileSuff, fAlphaOutputPath, fBlock, fCycleCounter, fDisableHistos, fErrorFlag, fGoodEventNumber, fTree, nP, nY, QwCorrelatorNew(), QwMessage, QwWarning, and VQwDataHandler::VQwDataHandler().

Here is the call graph for this function:

◆ ~QwCorrelatorNew()

QwCorrelatorNew::~QwCorrelatorNew ( )

override

Definition at line 48 of file QwCorrelatorNew.cc.

{
  // Close alpha and alias file
  CloseAlphaFile();
  CloseAliasFile();
}

References CloseAliasFile(), and CloseAlphaFile().

Here is the call graph for this function:

◆ QwCorrelatorNew() [3/3]

QwCorrelatorNew::QwCorrelatorNew ( )

private

Definition at line 701 of file QwCorrelatorNew.cc.

: nP(0),nY(0),
  fErrorFlag(-1),
  fGoodEventNumber(0)
{ }

References fErrorFlag, fGoodEventNumber, nP, and nY.

Member Function Documentation

◆ AccumulateRunningSum()

void QwCorrelatorNew::AccumulateRunningSum	(	VQwDataHandler &	value,
		Int_t	count = 0,
		Int_t	ErrorMask = 0xFFFFFFF )

overridevirtual

Reimplemented from VQwDataHandler.

Definition at line 132 of file QwCorrelatorNew.cc.

{
  QwCorrelatorNew* correlator = dynamic_cast<QwCorrelatorNew*>(&value);
  if (correlator) {
    operator+=(*correlator);
  } else {
    QwWarning << "QwCorrelatorNew::AccumulateRunningSum "
              << "can only accept other QwCorrelatorNew objects."
              << QwLog::endl;
  }
}

References QwLog::endl(), operator+=(), QwCorrelatorNew(), QwWarning, and VQwDataHandler::VQwDataHandler().

Here is the call graph for this function:

◆ CalcCorrelations()

void QwCorrelatorNew::CalcCorrelations ( )

Definition at line 144 of file QwCorrelatorNew.cc.

{
  // Check if any channels are active
  if (nP == 0 || nY == 0) {
    return;
  }
 
  QwMessage << "QwCorrelatorNew::CalcCorrelations(): name=" << GetName() << QwLog::endl;
 
  // Print entry summary
  QwVerbose << "QwCorrelatorNew: "
            << "total entries: " << fTotalCount << ", "
            << "good entries: " << fGoodCount
            << QwLog::endl;
  // and warn if zero
  if (fTotalCount > 100 && fGoodCount == 0) {
    QwWarning << "QwCorrelatorNew: "
              << "< 1% good events, "
              << fGoodCount << " of " << fTotalCount
              << QwLog::endl;
  }
 
  // Event error flag
  if (fErrCounts_EF > 0) {
    QwVerbose << "   Entries failed due to error flag: "
              << fErrCounts_EF << QwLog::endl;
  }
  // Dependent variable error codes
  for (size_t i = 0; i < fDependentVar.size(); ++i) {
    if (fErrCounts_DV.at(i) > 0) {
      QwVerbose << "   Entries failed due to " << fDependentVar.at(i)->GetElementName()
                << ": " <<  fErrCounts_DV.at(i) << QwLog::endl;
    }
  }
  // Independent variable error codes
  for (size_t i = 0; i < fIndependentVar.size(); ++i) {
    if (fErrCounts_IV.at(i) > 0) {
      QwVerbose << "   Entries failed due to " << fIndependentVar.at(i)->GetElementName()
                << ": " <<  fErrCounts_IV.at(i) << QwLog::endl;
    }
  }
 
  if (! this->failed()) {
 
    if (fPrintCorrelations) {
      this->printSummaryP();
      this->printSummaryY();
    }
 
    this->solve();
 
    if (fPrintCorrelations) {
      this->printSummaryAlphas();
      this->printSummaryMeansWithUnc();
      this->printSummaryMeansWithUncCorrected();
    }
  }
 
  // Fill tree
  if (fTree) fTree->Fill();
  else QwWarning << "No tree" << QwLog::endl;
 
  // Write alpha and alias file
  WriteAlphaFile();
  WriteAliasFile();
}

References QwLog::endl(), failed(), VQwDataHandler::fDependentVar, fErrCounts_DV, fErrCounts_EF, fErrCounts_IV, fGoodCount, fIndependentVar, fPrintCorrelations, fTotalCount, fTree, VQwDataHandler::GetName(), nP, nY, printSummaryAlphas(), printSummaryMeansWithUnc(), printSummaryMeansWithUncCorrected(), printSummaryP(), printSummaryY(), QwMessage, QwVerbose, QwWarning, solve(), WriteAliasFile(), and WriteAlphaFile().

Referenced by FinishDataHandler().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ clear()

void QwCorrelatorNew::clear ( )

Definition at line 775 of file QwCorrelatorNew.cc.

{
  mMP.Zero();
  mMY.Zero();
  mMYp.Zero();
 
  mVPP.Zero();
  mVPY.Zero();
  mVYP.Zero();
  mVYY.Zero();
  mVYYp.Zero();
  mVP.Zero();
  mVY.Zero();
  mVYp.Zero();
 
  mSPP.Zero();
  mSPY.Zero();
  mSYP.Zero();
  mSYY.Zero();
  mSYYp.Zero();
 
  Axy.Zero();
  Ayx.Zero();
  dAxy.Zero();
  dAyx.Zero();
 
  mSP.Zero();
  mSY.Zero();
  mSYp.Zero();
 
  mRPP.Zero();
  mRPY.Zero();
  mRYP.Zero();
  mRYY.Zero();
  mRYYp.Zero();
 
  fErrorFlag = -1;
  fGoodEventNumber = 0;
}

References Axy, Ayx, dAxy, dAyx, fErrorFlag, fGoodEventNumber, mMP, mMY, mMYp, mRPP, mRPY, mRYP, mRYY, mRYYp, mSP, mSPP, mSPY, mSY, mSYP, mSYp, mSYY, mSYYp, mVP, mVPP, mVPY, mVY, mVYP, mVYp, mVYY, and mVYYp.

Referenced by ClearEventData().

Here is the caller graph for this function:

◆ ClearEventData()

void QwCorrelatorNew::ClearEventData ( )

overridevirtual

Reimplemented from VQwDataHandler.

Definition at line 116 of file QwCorrelatorNew.cc.

{
  // Clear error counters
  fErrCounts_EF = 0;
  std::fill(fErrCounts_DV.begin(), fErrCounts_DV.end(), 0);
  std::fill(fErrCounts_IV.begin(), fErrCounts_IV.end(), 0);
 
  // Clear event counts
  fTotalCount = 0;
  fGoodCount = 0;
  fGoodEvent = -1;
 
  // Clear regression
  this->clear();
}

References clear(), fErrCounts_DV, fErrCounts_EF, fErrCounts_IV, fGoodCount, fGoodEvent, and fTotalCount.

Referenced by QwCorrelatorNew(), and QwCorrelatorNew().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ CloseAliasFile()

void QwCorrelatorNew::CloseAliasFile ( )

protected

Definition at line 657 of file QwCorrelatorNew.cc.

{
  // Close alias output file
  if (fAliasOutputFile.good()) {
    fAliasOutputFile << "}" << std::endl << std::endl;
    fAliasOutputFile.close();
  } else {
    QwWarning << "QwCorrelatorNew: Unable to close alias output file." << QwLog::endl;
  }
}

References QwLog::endl(), fAliasOutputFile, and QwWarning.

Referenced by ~QwCorrelatorNew().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ CloseAlphaFile()

void QwCorrelatorNew::CloseAlphaFile ( )

protected

Definition at line 648 of file QwCorrelatorNew.cc.

{
  // Close slopes output file
  if (fAlphaOutputFile) {
    fAlphaOutputFile->Write();
    fAlphaOutputFile->Close();
  }
}

References fAlphaOutputFile.

Referenced by ~QwCorrelatorNew().

Here is the caller graph for this function:

◆ ConnectChannels()

Int_t QwCorrelatorNew::ConnectChannels	(	QwSubsystemArrayParity &	asym,
		QwSubsystemArrayParity &	diff )

overridevirtual

Connect to Channels (asymmetry/difference only)

Parameters

asym	Subsystem array with asymmetries
diff	Subsystem array with differences

Returns: 0 on success, non-zero on failure

Fill vector of pointers to the relevant data elements

Reimplemented from VQwDataHandler.

Definition at line 259 of file QwCorrelatorNew.cc.

{
  SetEventcutErrorFlagPointer(asym.GetEventcutErrorFlagPointer());
 
  // Return if correlator is not enabled
 
  /// Fill vector of pointers to the relevant data elements
  for (size_t dv = 0; dv < fDependentName.size(); dv++) {
    // Get the dependent variables
 
    const VQwHardwareChannel* dv_ptr = 0;
 
    if (fDependentType.at(dv)==kHandleTypeMps){
      //  Quietly ignore the MPS type when we're connecting the asym & diff
      continue;
    }else{
      dv_ptr = this->RequestExternalPointer(fDependentFull.at(dv));
      if (dv_ptr==NULL){
    switch (fDependentType.at(dv)) {
        case kHandleTypeAsym:
          dv_ptr = asym.RequestExternalPointer(fDependentName.at(dv));
          break;
        case kHandleTypeDiff:
          dv_ptr = diff.RequestExternalPointer(fDependentName.at(dv));
          break;
        default:
          QwWarning << "QwCorrelatorNew::ConnectChannels(QwSubsystemArrayParity& asym, QwSubsystemArrayParity& diff): "
                    << "Dependent variable, " << fDependentName.at(dv)
                    << ", for asym/diff correlator does not have proper type, type=="
                    << fDependentType.at(dv) << "." << QwLog::endl;
          break;
        }
      }
      if (dv_ptr == NULL){
    QwWarning << "QwCombiner::ConnectChannels(QwSubsystemArrayParity& asym, QwSubsystemArrayParity& diff):  Dependent variable, "
          << fDependentName.at(dv)
          << ", was not found (fullname=="
          << fDependentFull.at(dv)<< ")." << QwLog::endl;
     continue;
      }
    }
 
    // pair creation
    if(dv_ptr != NULL){
      // fDependentVarType.push_back(fDependentType.at(dv));
      fDependentVar.push_back(dv_ptr);
    }
 
  }
 
  // Add independent variables
  for (size_t iv = 0; iv < fIndependentName.size(); iv++) {
    // Get the independent variables
    const VQwHardwareChannel* iv_ptr = 0;
    iv_ptr = this->RequestExternalPointer(fIndependentFull.at(iv));
    if (iv_ptr==NULL){
      switch (fIndependentType.at(iv)) {
      case kHandleTypeAsym:
        iv_ptr = asym.RequestExternalPointer(fIndependentName.at(iv));
        break;
      case kHandleTypeDiff:
        iv_ptr = diff.RequestExternalPointer(fIndependentName.at(iv));
        break;
      default:
        QwWarning << "Independent variable for correlator has unknown type."
                  << QwLog::endl;
        break;
      }
    }
    if (iv_ptr) {
      fIndependentVar.push_back(iv_ptr);
    } else {
      QwWarning << "Independent variable " << fIndependentName.at(iv) << " for correlator could not be found."
                << QwLog::endl;
    }
 
  }
  fIndependentValues.resize(fIndependentVar.size());
  fDependentValues.resize(fDependentVar.size());
 
  nP = fIndependentName.size();
  nY = fDependentName.size();
 
  this->setDims(nP, nY);
  this->init();
 
  fErrCounts_IV.resize(fIndependentVar.size(),0);
  fErrCounts_DV.resize(fDependentVar.size(),0);
 
  return 0;
}

References QwLog::endl(), VQwDataHandler::fDependentFull, VQwDataHandler::fDependentName, VQwDataHandler::fDependentType, VQwDataHandler::fDependentValues, VQwDataHandler::fDependentVar, fErrCounts_DV, fErrCounts_IV, fIndependentFull, fIndependentName, fIndependentType, fIndependentValues, fIndependentVar, QwSubsystemArrayParity::GetEventcutErrorFlagPointer(), init(), VQwDataHandler::kHandleTypeAsym, VQwDataHandler::kHandleTypeDiff, VQwDataHandler::kHandleTypeMps, nP, nY, QwWarning, MQwPublishable< U, T >::RequestExternalPointer(), MQwPublishable_child< QwDataHandlerArray, VQwDataHandler >::RequestExternalPointer(), setDims(), and VQwDataHandler::SetEventcutErrorFlagPointer().

Here is the call graph for this function:

◆ ConstructHistograms()

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

overridevirtual

Construct the histograms in a folder with a prefix.

Reimplemented from VQwDataHandler.

Definition at line 440 of file QwCorrelatorNew.cc.

{
  // Skip if disabled
  if (fDisableHistos) return;
 
  // Check if any channels are active
  if (nP == 0 || nY == 0) {
    return;
  }
 
  // Go to directory
  TString name(fName);
  name.ReplaceAll(" ","_");
  folder->mkdir(name)->cd();
 
  //..... 1D,  iv
  fH1iv.resize(nP);
  for (int i = 0; i < nP; i++) {
    fH1iv[i] = TH1D(
        Form("P%d",i),
        Form("iv P%d=%s, pass=%s ;iv=%s (ppm)",i,fIndependentName[i].c_str(),fName.Data(),fIndependentName[i].c_str()),
        128,0.,0.);
    fH1iv[i].GetXaxis()->SetNdivisions(4);
  }
 
  //..... 2D,  iv correlations
  Double_t x1 = 0;
  fH2iv.resize(nP);
  for (int i = 0; i < nP; i++) {
    fH2iv[i].resize(nP);
    for (int j = i+1; j < nP; j++) { // not all are used
      fH2iv[i][j] = TH2D(
          Form("P%d_P%d",i,j),
          Form("iv correlation  P%d_P%d, pass=%s ;P%d=%s (ppm);P%d=%s   (ppm)  ",
              i,j,fName.Data(),i,fIndependentName[i].c_str(),j,fIndependentName[j].c_str()),
          64,-x1,x1,
          64,-x1,x1);
      fH2iv[i][j].GetXaxis()->SetTitleColor(kBlue);
      fH2iv[i][j].GetYaxis()->SetTitleColor(kBlue);
      fH2iv[i][j].GetXaxis()->SetNdivisions(4);
      fH2iv[i][j].GetYaxis()->SetNdivisions(4);
    }
  }
 
  //..... 1D,  dv
  fH1dv.resize(nY);
  for (int i = 0; i < nY; i++) {
    fH1dv[i] = TH1D(
        Form("Y%d",i),
        Form("dv Y%d=%s, pass=%s ;dv=%s (ppm)",i,fDependentName[i].c_str(),fName.Data(),fDependentName[i].c_str()),
        128,0.,0.);
    fH1dv[i].GetXaxis()->SetNdivisions(4);
  }
 
  //..... 2D,  dv-iv correlations
  Double_t y1 = 0;
  fH2dv.resize(nP);
  for (int i = 0; i < nP; i++) {
    fH2dv[i].resize(nY);
    for (int j = 0; j < nY; j++) {
      fH2dv[i][j] = TH2D(
          Form("P%d_Y%d",i,j),
          Form("iv-dv correlation  P%d_Y%d, pass=%s ;P%d=%s (ppm);Y%d=%s   (ppm)  ",
              i,j,fName.Data(),i,fIndependentName[i].c_str(),j,fDependentName[j].c_str()),
          64,-x1,x1,
          64,-y1,y1);
      fH2dv[i][j].GetXaxis()->SetTitleColor(kBlue);
      fH2dv[i][j].GetYaxis()->SetTitleColor(kBlue);
      fH2dv[i][j].GetXaxis()->SetNdivisions(4);
      fH2dv[i][j].GetYaxis()->SetNdivisions(4);
    }
  }
 
  // store list of names to be archived
  fHnames.resize(2);
  fHnames[0] = TH1D("NamesIV",Form("IV name list nIV=%d",nP),nP,0,1);
  for (int i = 0; i < nP; i++)
    fHnames[0].Fill(fIndependentName[i].c_str(),1.*i);
  fHnames[1] = TH1D("NamesDV",Form("DV name list nIV=%d",nY),nY,0,1);
  for (int i = 0; i < nY; i++)
    fHnames[1].Fill(fDependentName[i].c_str(),i*1.);
}

References VQwDataHandler::fDependentName, fDisableHistos, fH1dv, fH1iv, fH2dv, fH2iv, fHnames, fIndependentName, VQwDataHandler::fName, nP, and nY.

◆ ConstructTreeBranches()

void QwCorrelatorNew::ConstructTreeBranches	(	QwRootFile *	treerootfile,
		const std::string &	treeprefix = "",
		const std::string &	branchprefix = "" )

overridevirtual

Construct the tree branches.

Reimplemented from VQwDataHandler.

Definition at line 352 of file QwCorrelatorNew.cc.

{
  // Check if any channels are active
  if (nP == 0 || nY == 0) {
    return;
  }
 
  // Check if tree name is specified
  if (fTreeName == "") {
    QwWarning << "QwCorrelatorNew: no tree name specified, use 'tree-name = value'" << QwLog::endl;
    return;
  }
 
  // Create alpha and alias files before trying to create the tree
  OpenAlphaFile(treeprefix);
  OpenAliasFile(treeprefix);
 
  // Construct tree name and create new tree
  const std::string name = treeprefix + fTreeName;
  treerootfile->NewTree(name, fTreeComment.c_str());
  fTree = treerootfile->GetTree(name);
  // Check to make sure the tree was created successfully
  if (fTree == NULL) return;
 
  // Set up branches
  fTree->Branch(TString(branchprefix + "total_count"), &fTotalCount);
  fTree->Branch(TString(branchprefix + "good_count"),  &fGoodCount);
 
  fTree->Branch(TString(branchprefix + "n"), &(this->fGoodEventNumber));
  fTree->Branch(TString(branchprefix + "ErrorFlag"), &(this->fErrorFlag));
 
  auto bn = [&](const TString& n) {
    return TString(branchprefix + n);
  };
  auto pm = [](TMatrixD& m) {
    return m.GetMatrixArray();
  };
  auto lm = [](TMatrixD& m, const TString& n) {
    return Form("%s[%d][%d]/D", n.Data(), m.GetNrows(), m.GetNcols());
  };
  auto branchm = [&](TTree* tree, TMatrixD& m, const TString& n) {
    tree->Branch(bn(n),pm(m),lm(m,n));
  };
  auto pv = [](TVectorD& v) {
    return v.GetMatrixArray();
  };
  auto lv = [](TVectorD& v, const TString& n) {
    return Form("%s[%d]/D", n.Data(), v.GetNrows());
  };
  auto branchv = [&](TTree* tree, TVectorD& v, const TString& n) {
    tree->Branch(bn(n),pv(v),lv(v,n));
  };
 
  branchm(fTree,this->Axy,  "A");
  branchm(fTree,this->dAxy, "dA");
 
  branchm(fTree,this->mVPP,  "VPP");
  branchm(fTree,this->mVPY,  "VPY");
  branchm(fTree,this->mVYP,  "VYP");
  branchm(fTree,this->mVYY,  "VYY");
  branchm(fTree,this->mVYYp, "VYYp");
 
  branchm(fTree,this->mSPP,  "SPP");
  branchm(fTree,this->mSPY,  "SPY");
  branchm(fTree,this->mSYP,  "SYP");
  branchm(fTree,this->mSYY,  "SYY");
  branchm(fTree,this->mSYYp, "SYYp");
 
  branchm(fTree,this->mRPP,  "RPP");
  branchm(fTree,this->mRPY,  "RPY");
  branchm(fTree,this->mRYP,  "RYP");
  branchm(fTree,this->mRYY,  "RYY");
  branchm(fTree,this->mRYYp, "RYYp");
 
  branchv(fTree,this->mMP,  "MP");   // Parameter mean
  branchv(fTree,this->mMY,  "MY");   // Uncorrected mean
  branchv(fTree,this->mMYp, "MYp");  // Corrected mean
 
  branchv(fTree,this->mSP,  "dMP");  // Parameter mean error
  branchv(fTree,this->mSY,  "dMY");  // Uncorrected mean error
  branchv(fTree,this->mSYp, "dMYp"); // Corrected mean error
 
}

References Axy, dAxy, QwLog::endl(), fErrorFlag, fGoodCount, fGoodEventNumber, fTotalCount, fTree, VQwDataHandler::fTreeComment, VQwDataHandler::fTreeName, QwRootFile::GetTree(), mMP, mMY, mMYp, mRPP, mRPY, mRYP, mRYY, mRYYp, mSP, mSPP, mSPY, mSY, mSYP, mSYp, mSYY, mSYYp, mVPP, mVPY, mVYP, mVYY, mVYYp, QwRootFile::NewTree(), nP, nY, OpenAliasFile(), OpenAlphaFile(), and QwWarning.

Here is the call graph for this function:

◆ DefineOptions()

void QwCorrelatorNew::DefineOptions ( QwOptions & options )

static

Definition at line 55 of file QwCorrelatorNew.cc.

{
  options.AddOptions()("print-correlations",
      po::value<bool>(&fPrintCorrelations)->default_bool_value(false),
      "print correlations after determining them");
}

References QwOptions::AddOptions(), and fPrintCorrelations.

Here is the call graph for this function:

◆ failed()

bool QwCorrelatorNew::failed ( )

inline

Definition at line 167 of file QwCorrelatorNew.h.

167{ return fGoodEventNumber < nP + 1; }

References fGoodEventNumber, and nP.

Referenced by CalcCorrelations().

Here is the caller graph for this function:

◆ FillHistograms()

void QwCorrelatorNew::FillHistograms ( )

overridevirtual

Fill the histograms.

Reimplemented from VQwDataHandler.

Definition at line 524 of file QwCorrelatorNew.cc.

{
  // Skip if disabled
  if (fDisableHistos) return;
 
  // Check if any channels are active
  if (nP == 0 || nY == 0) {
    return;
  }
 
  // Skip if bad event
  if (fGoodEvent != 0) return;
 
  // Fill histograms
  for (size_t i = 0; i < fIndependentValues.size(); i++) {
    fH1iv[i].Fill(fIndependentValues[i]);
    for (size_t j = i+1; j < fIndependentValues.size(); j++)
      fH2iv[i][j].Fill(fIndependentValues[i], fIndependentValues[j]);
  }
  for (size_t j = 0; j < fDependentValues.size(); j++) {
    fH1dv[j].Fill(fDependentValues[j]);
    for (size_t i = 0; i < fIndependentValues.size(); i++)
      fH2dv[i][j].Fill(fIndependentValues[i], fDependentValues[j]);
  }
}

References VQwDataHandler::fDependentValues, fDisableHistos, fGoodEvent, fH1dv, fH1iv, fH2dv, fH2iv, fIndependentValues, nP, and nY.

◆ FillTreeBranches()

void QwCorrelatorNew::FillTreeBranches ( QwRootFile * treerootfile )

inlineoverridevirtual

Fill the tree branches.

Reimplemented from VQwDataHandler.

Definition at line 62 of file QwCorrelatorNew.h.

62{ };

◆ FinishDataHandler()

void QwCorrelatorNew::FinishDataHandler ( )

inlineoverridevirtual

Reimplemented from VQwDataHandler.

Definition at line 51 of file QwCorrelatorNew.h.

                                   {
    CalcCorrelations();
  }

References CalcCorrelations().

Here is the call graph for this function:

◆ getCovarianceP()

Int_t QwCorrelatorNew::getCovarianceP	(	int	i,
		int	j,
		Double_t &	covar ) const

Get mean value of a variable, returns error code.

Definition at line 979 of file QwCorrelatorNew.cc.

{
    covar=-1e50;
    if( i>j) { int k=i; i=j; j=k; }//swap i & j
    //... now we need only upper right triangle
    if(i<0 || i >= nP ) return -11;
    if( fGoodEventNumber<2) return -14;
    covar=mVPP(i,j)/(fGoodEventNumber-1.);
    return 0;
}

References fGoodEventNumber, mVPP, and nP.

Referenced by printSummaryP().

Here is the caller graph for this function:

◆ getCovariancePY()

Int_t QwCorrelatorNew::getCovariancePY	(	int	ip,
		int	iy,
		Double_t &	covar ) const

Definition at line 992 of file QwCorrelatorNew.cc.

{
    covar=-1e50;
    //... now we need only upper right triangle
    if(ip<0 || ip >= nP ) return -11;
    if(iy<0 || iy >= nY ) return -12;
    if( fGoodEventNumber<2) return -14;
    covar=mVPY(ip,iy)/(fGoodEventNumber-1.);
    return 0;
}

References fGoodEventNumber, mVPY, nP, and nY.

Referenced by printSummaryYP().

Here is the caller graph for this function:

◆ getCovarianceY()

Int_t QwCorrelatorNew::getCovarianceY	(	int	i,
		int	j,
		Double_t &	covar ) const

Definition at line 1005 of file QwCorrelatorNew.cc.

{
    covar=-1e50;
    if( i>j) { int k=i; i=j; j=k; }//swap i & j
    //... now we need only upper right triangle
    if(i<0 || i >= nY ) return -11;
    if( fGoodEventNumber<2) return -14;
    covar=mVYY(i,j)/(fGoodEventNumber-1.);
    return 0;
}

References fGoodEventNumber, mVYY, and nY.

◆ getMeanP()

Int_t QwCorrelatorNew::getMeanP	(	const int	i,
		Double_t &	mean ) const

Get mean value of a variable, returns error code.

Definition at line 912 of file QwCorrelatorNew.cc.

{
   mean=-1e50;
   if(i<0 || i >= nP ) return -1;
   if( fGoodEventNumber<1) return -3;
   mean = mMP(i);    return 0;
}

References fGoodEventNumber, mMP, and nP.

Referenced by printSummaryP().

Here is the caller graph for this function:

◆ getMeanY()

Int_t QwCorrelatorNew::getMeanY	(	const int	i,
		Double_t &	mean ) const

Definition at line 923 of file QwCorrelatorNew.cc.

{
  mean=-1e50;
  if(i<0 || i >= nY ) return -1;
  if( fGoodEventNumber<1) return -3;
  mean = mMY(i);    return 0;
}

References fGoodEventNumber, mMY, and nY.

Referenced by printSummaryY(), and printSummaryYP().

Here is the caller graph for this function:

◆ getMeanYprime()

Int_t QwCorrelatorNew::getMeanYprime	(	const int	i,
		Double_t &	mean ) const

Definition at line 934 of file QwCorrelatorNew.cc.

{
  mean=-1e50;
  if(i<0 || i >= nY ) return -1;
  if( fGoodEventNumber<1) return -3;
  mean = mMYp(i);    return 0;
}

References fGoodEventNumber, mMYp, and nY.

◆ getSigmaP()

Int_t QwCorrelatorNew::getSigmaP	(	const int	i,
		Double_t &	sigma ) const

Get mean value of a variable, returns error code.

Definition at line 945 of file QwCorrelatorNew.cc.

{
  sigma=-1e50;
  if(i<0 || i >= nP ) return -1;
  if( fGoodEventNumber<2) return -3;
  sigma=sqrt(mVPP(i,i)/(fGoodEventNumber-1.));
  return 0;
}

References fGoodEventNumber, mVPP, and nP.

Referenced by printSummaryP(), and printSummaryYP().

Here is the caller graph for this function:

◆ getSigmaY()

Int_t QwCorrelatorNew::getSigmaY	(	const int	i,
		Double_t &	sigma ) const

Definition at line 957 of file QwCorrelatorNew.cc.

{
  sigma=-1e50;
  if(i<0 || i >= nY ) return -1;
  if( fGoodEventNumber<2) return -3;
  sigma=sqrt(mVYY(i,i)/(fGoodEventNumber-1.));
  return 0;
}

References fGoodEventNumber, mVYY, and nY.

Referenced by printSummaryY(), and printSummaryYP().

Here is the caller graph for this function:

◆ getSigmaYprime()

Int_t QwCorrelatorNew::getSigmaYprime	(	const int	i,
		Double_t &	sigma ) const

Definition at line 968 of file QwCorrelatorNew.cc.

{
  sigma=-1e50;
  if(i<0 || i >= nY ) return -1;
  if( fGoodEventNumber<2) return -3;
  sigma=sqrt(mVYYp(i,i)/(fGoodEventNumber-1.));
  return 0;
}

References fGoodEventNumber, mVYYp, and nY.

◆ getUsedEve()

double QwCorrelatorNew::getUsedEve ( ) const

inline

Definition at line 197 of file QwCorrelatorNew.h.

197{ return fGoodEventNumber; };

References fGoodEventNumber.

Referenced by WriteAlphaFile().

Here is the caller graph for this function:

◆ init()

void QwCorrelatorNew::init ( )

Definition at line 736 of file QwCorrelatorNew.cc.

{
  mMP.ResizeTo(nP);
  mMY.ResizeTo(nY);
  mMYp.ResizeTo(nY);
 
  mVPP.ResizeTo(nP,nP);
  mVPY.ResizeTo(nP,nY);
  mVYP.ResizeTo(nY,nP);
  mVYY.ResizeTo(nY,nY);
  mVYYp.ResizeTo(nY,nY);
  mVP.ResizeTo(nP);
  mVY.ResizeTo(nY);
  mVYp.ResizeTo(nY);
 
  mSPP.ResizeTo(mVPP);
  mSPY.ResizeTo(mVPY);
  mSYP.ResizeTo(mVYP);
  mSYY.ResizeTo(mVYY);
  mSYYp.ResizeTo(mVYYp);
 
  Axy.ResizeTo(nP,nY);
  Ayx.ResizeTo(nY,nP);
  dAxy.ResizeTo(Axy);
  dAyx.ResizeTo(Ayx);
 
  mSP.ResizeTo(nP);
  mSY.ResizeTo(nY);
  mSYp.ResizeTo(nY);
 
  mRPP.ResizeTo(mVPP);
  mRPY.ResizeTo(mVPY);
  mRYP.ResizeTo(mVYP);
  mRYY.ResizeTo(mVYY);
  mRYYp.ResizeTo(mVYYp);
 
  fGoodEventNumber = 0;
}

References Axy, Ayx, dAxy, dAyx, fGoodEventNumber, mMP, mMY, mMYp, mRPP, mRPY, mRYP, mRYY, mRYYp, mSP, mSPP, mSPY, mSY, mSYP, mSYp, mSYY, mSYYp, mVP, mVPP, mVPY, mVY, mVYP, mVYp, mVYY, mVYYp, nP, and nY.

Referenced by ConnectChannels().

Here is the caller graph for this function:

◆ LoadChannelMap()

Int_t QwCorrelatorNew::LoadChannelMap ( const std::string & mapfile )

overridevirtual

Load the channel map

Parameters

mapfile Filename of map file

Returns: Zero when success

Reimplemented from VQwDataHandler.

Definition at line 217 of file QwCorrelatorNew.cc.

{
  // Open the file
  QwParameterFile map(mapfile);
 
  // Read the sections of dependent variables
  std::pair<EQwHandleType,std::string> type_name;
 
  // Add independent variables and sensitivities
  while (map.ReadNextLine()) {
    // Throw away comments, whitespace, empty lines
    map.TrimComment();
    map.TrimWhitespace();
    if (map.LineIsEmpty()) continue;
    // Get first token: label (dv or iv), second token is the name like "asym_blah"
    string primary_token = map.GetNextToken(" ");
    string current_token = map.GetNextToken(" ");
    // Parse current token into independent variable type and name
    type_name = ParseHandledVariable(current_token);
 
    if (primary_token == "iv") {
      fIndependentType.push_back(type_name.first);
      fIndependentName.push_back(type_name.second);
      fIndependentFull.push_back(current_token);
    }
    else if (primary_token == "dv") {
      fDependentType.push_back(type_name.first);
      fDependentName.push_back(type_name.second);
      fDependentFull.push_back(current_token);
    }
    else if (primary_token == "treetype") {
      QwMessage << "Tree Type read, ignoring." << QwLog::endl;
    }
    else {
      QwError << "LoadChannelMap in QwCorrelatorNew read invalid primary_token " << primary_token << QwLog::endl;
    }
  }
 
  return 0;
}

References QwLog::endl(), VQwDataHandler::fDependentFull, VQwDataHandler::fDependentName, VQwDataHandler::fDependentType, fIndependentFull, fIndependentName, fIndependentType, QwParameterFile::GetNextToken(), QwParameterFile::LineIsEmpty(), VQwDataHandler::ParseHandledVariable(), QwError, QwMessage, QwParameterFile::ReadNextLine(), QwParameterFile::TrimComment(), and QwParameterFile::TrimWhitespace().

Here is the call graph for this function:

◆ OpenAliasFile()

void QwCorrelatorNew::OpenAliasFile ( const std::string & prefix )

protected

Definition at line 630 of file QwCorrelatorNew.cc.

{
  // Turn "." into "_" in run_label (no "." allowed in function name, and must
  // agree with the filename)
  std::string label(run_label);
  std::replace(label.begin(), label.end(), '.', '_');
  // Create old-style regalias script
  std::string name = prefix + fAliasOutputFileBase + label + fAliasOutputFileSuff;
  std::string path = fAliasOutputPath + "/";
  std::string file = path + name + ".C"; // add extension outside of file suffix
  fAliasOutputFile.open(file, std::ofstream::out);
  if (fAliasOutputFile.good()) {
    fAliasOutputFile << Form("void %s(int i = 0) {", name.c_str()) << std::endl;
  } else {
    QwWarning << "QwCorrelatorNew: Could not write to alias output file " << QwLog::endl;
  }
}

References QwLog::endl(), fAliasOutputFile, fAliasOutputFileBase, fAliasOutputFileSuff, fAliasOutputPath, QwWarning, and VQwDataHandler::run_label.

Referenced by ConstructTreeBranches().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ OpenAlphaFile()

void QwCorrelatorNew::OpenAlphaFile ( const std::string & prefix )

protected

Definition at line 616 of file QwCorrelatorNew.cc.

{
  // Create old-style blueR ROOT file
  std::string name = prefix + fAlphaOutputFileBase + run_label.Data() + fAlphaOutputFileSuff;
  std::string path = fAlphaOutputPath + "/";
  std::string file = path + name;
  fAlphaOutputFile = new TFile(TString(file), "RECREATE", "correlation coefficients");
  if (! fAlphaOutputFile->IsWritable()) {
    QwError << "QwCorrelatorNew could not create output file " << file << QwLog::endl;
    delete fAlphaOutputFile;
    fAlphaOutputFile = 0;
  }
}

References QwLog::endl(), fAlphaOutputFile, fAlphaOutputFileBase, fAlphaOutputFileSuff, fAlphaOutputPath, QwError, and VQwDataHandler::run_label.

Referenced by ConstructTreeBranches().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ operator+=() [1/2]

QwCorrelatorNew & QwCorrelatorNew::operator+= ( const QwCorrelatorNew & rhs )

Definition at line 871 of file QwCorrelatorNew.cc.

{
  // If set X = A + B, then
  //   Cov[X] = Cov[A] + Cov[B]
  //          + (E[x_A] - E[x_B]) * (E[y_A] - E[y_B]) * n_A * n_B / n_X
  // Ref: E. Schubert, M. Gertz (9 July 2018).
  // "Numerically stable parallel computation of (co-)variance".
  // SSDBM '18 Proceedings of the 30th International Conference
  // on Scientific and Statistical Database Management.
  // https://doi.org/10.1145/3221269.3223036
 
  if(fGoodEventNumber + rhs.fGoodEventNumber == 0)
    return *this;
 
  // Deviations from mean
  TVectorD delta_y(mMY - rhs.mMY);
  TVectorD delta_p(mMP - rhs.mMP);
 
  // Update covariances
  Double_t alpha = fGoodEventNumber * rhs.fGoodEventNumber
                / (fGoodEventNumber + rhs.fGoodEventNumber);
  mVYY += rhs.mVYY;
  mVYY.Rank1Update(delta_y, alpha);
  mVPY += rhs.mVPY;
  mVPY.Rank1Update(delta_p, delta_y, alpha);
  mVPP += rhs.mVPP;
  mVPP.Rank1Update(delta_p, alpha);
 
  // Update means
  Double_t beta = rhs.fGoodEventNumber / (fGoodEventNumber + rhs.fGoodEventNumber);
  mMY += delta_y * beta;
  mMP += delta_p * beta;
 
  fGoodEventNumber += rhs.fGoodEventNumber;
 
  return *this;
}

References fGoodEventNumber, mMP, mMY, mVPP, mVPY, mVYY, and QwCorrelatorNew().

Here is the call graph for this function:

◆ operator+=() [2/2]

QwCorrelatorNew & QwCorrelatorNew::operator+= ( const std::pair< TVectorD, TVectorD > & rhs )

Definition at line 832 of file QwCorrelatorNew.cc.

{
  // Get independent and dependent components
  const TVectorD& P = rhs.first;
  const TVectorD& Y = rhs.second;
 
  // Update number of events
  fGoodEventNumber++;
 
  if (fGoodEventNumber <= 1) {
    // First event, set covariances to zero and means to first value
    mVPP.Zero();
    mVPY.Zero();
    mVYY.Zero();
    mMP = P;
    mMY = Y;
  } else {
    // Deviations from mean
    TVectorD delta_y(Y - mMY);
    TVectorD delta_p(P - mMP);
 
    // Update covariances
    Double_t alpha = (fGoodEventNumber - 1.0) / fGoodEventNumber;
    mVPP.Rank1Update(delta_p, alpha);
    mVPY.Rank1Update(delta_p, delta_y, alpha);
    mVYY.Rank1Update(delta_y, alpha);
 
    // Update means
    Double_t beta = 1.0 / fGoodEventNumber;
    mMP += delta_p * beta;
    mMY += delta_y * beta;
  }
 
  return *this;
}

References fGoodEventNumber, mMP, mMY, mVPP, mVPY, mVYY, and QwCorrelatorNew().

Referenced by AccumulateRunningSum(), and ProcessData().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ParseConfigFile()

void QwCorrelatorNew::ParseConfigFile ( QwParameterFile & file )

overridevirtual

Parse configuration file for map file, priority, tree settings.

Reimplemented from VQwDataHandler.

Definition at line 66 of file QwCorrelatorNew.cc.

{
  VQwDataHandler::ParseConfigFile(file);
  file.PopValue("slope-file-base", fAlphaOutputFileBase);
  file.PopValue("slope-file-suff", fAlphaOutputFileSuff);
  file.PopValue("slope-path", fAlphaOutputPath);
  file.PopValue("alias-file-base", fAliasOutputFileBase);
  file.PopValue("alias-file-suff", fAliasOutputFileSuff);
  file.PopValue("alias-path", fAliasOutputPath);
  file.PopValue("disable-histos", fDisableHistos);
  file.PopValue("block", fBlock);
  if (fBlock >= 4)
    QwWarning << "QwCorrelatorNew: expect 0 <= block <= 3 but block = "
              << fBlock << QwLog::endl;
}

References QwLog::endl(), fAliasOutputFileBase, fAliasOutputFileSuff, fAliasOutputPath, fAlphaOutputFileBase, fAlphaOutputFileSuff, fAlphaOutputPath, fBlock, fDisableHistos, VQwDataHandler::ParseConfigFile(), QwParameterFile::PopValue(), and QwWarning.

Here is the call graph for this function:

◆ print()

void QwCorrelatorNew::print ( )

Definition at line 817 of file QwCorrelatorNew.cc.

{
  QwMessage << "LinReg dims:  nP=" << nP << " nY=" << nY << QwLog::endl;
 
  QwMessage << "MP:"; mMP.Print();
  QwMessage << "MY:"; mMY.Print();
  QwMessage << "VPP:"; mVPP.Print();
  QwMessage << "VPY:"; mVPY.Print();
  QwMessage << "VYY:"; mVYY.Print();
  QwMessage << "VYYprime:"; mVYYp.Print();
}

References QwLog::endl(), mMP, mMY, mVPP, mVPY, mVYY, mVYYp, nP, nY, and QwMessage.

Here is the call graph for this function:

◆ printSummaryAlphas()

void QwCorrelatorNew::printSummaryAlphas ( ) const

Definition at line 1075 of file QwCorrelatorNew.cc.

{
  QwMessage << Form("\nLinRegBevPeb::printSummaryAlphas seen good eve=%lld",fGoodEventNumber)<<QwLog::endl;
  QwMessage << Form("\n  j                slope         sigma     mean/sigma\n");
  for (int iy = 0; iy <nY; iy++) {
    QwMessage << Form("dv=Y%d: ",iy)<<QwLog::endl;
    for (int j = 0; j < nP; j++) {
      double val=Axy(j,iy);
      double err=dAxy(j,iy);
      double nSig=val/err;
      char x=' ';
      if(fabs(nSig)>3.) x='*';
      QwMessage << Form("  slope_%d = %11.3g  +/-%11.3g  (nSig=%.2f) %c\n",j,val, err,nSig,x);
    }
  }
}

References Axy, dAxy, QwLog::endl(), fGoodEventNumber, nP, nY, and QwMessage.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ printSummaryMeansWithUnc()

void QwCorrelatorNew::printSummaryMeansWithUnc ( ) const

Definition at line 1125 of file QwCorrelatorNew.cc.

{
  QwMessage << "Uncorrected Y values:" << QwLog::endl;
  QwMessage << "     mean          sig" << QwLog::endl;
  for (int i = 0; i < nY; i++){
    QwMessage << "Y" << i << ":  " << mMY(i) << " +- " << mSY(i) << QwLog::endl;
  }
  QwMessage << QwLog::endl;
}

References QwLog::endl(), mMY, mSY, nY, and QwMessage.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ printSummaryMeansWithUncCorrected()

void QwCorrelatorNew::printSummaryMeansWithUncCorrected ( ) const

Definition at line 1138 of file QwCorrelatorNew.cc.

{
  QwMessage << "Corrected Y values:" << QwLog::endl;
  QwMessage << "     mean          sig" << QwLog::endl;
  for (int i = 0; i < nY; i++){
    QwMessage << "Y" << i << ":  " << mMYp(i) << " +- " << mSYp(i) << QwLog::endl;
  }
  QwMessage << QwLog::endl;
}

References QwLog::endl(), mMYp, mSYp, nY, and QwMessage.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ printSummaryP()

void QwCorrelatorNew::printSummaryP ( ) const

Definition at line 1018 of file QwCorrelatorNew.cc.

{
  QwMessage << Form("\nLinRegBevPeb::printSummaryP seen good eve=%lld",fGoodEventNumber)<<QwLog::endl;
 
  size_t dim=nP;
  if(fGoodEventNumber>2) { // print full matrix
    QwMessage << Form("\nname:                                               ");
    for (size_t i = 1; i <dim; i++) {
      QwMessage << Form("P%d%11s",(int)i," ");
    }
    QwMessage << Form("\n           mean     sig(distrib)   nSig(mean)       correlation-matrix ....\n");
    for (size_t i = 0; i <dim; i++) {
      double meanI,sigI;
      if (getMeanP(i,meanI) < 0) QwWarning << "LRB::getMeanP failed" << QwLog::endl;
      if (getSigmaP(i,sigI) < 0) QwWarning << "LRB::getSigmaP failed" << QwLog::endl;
      double nSig=-1;
      double err=sigI/sqrt(fGoodEventNumber);
      if(sigI>0.) nSig=meanI/err;
 
      QwMessage << Form("P%d:  %+12.4g  %12.3g      %.1f    ",(int)i,meanI,sigI,nSig);
      for (size_t j = 1; j <dim; j++) {
        if( j<=i) { QwMessage << Form("  %12s","._._._."); continue;}
        double sigJ,cov;
        if (getSigmaP(j,sigJ) < 0) QwWarning << "LRB::getSigmaP failed" << QwLog::endl;
        if (getCovarianceP(i,j,cov) < 0) QwWarning << "LRB::getCovarianceP failed" << QwLog::endl;
        double corel=cov / sigI / sigJ;
 
        QwMessage << Form("  %12.3g",corel);
      }
      QwMessage << Form("\n");
    }
  }
}

References QwLog::endl(), fGoodEventNumber, getCovarianceP(), getMeanP(), getSigmaP(), nP, QwMessage, and QwWarning.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ printSummaryY()

void QwCorrelatorNew::printSummaryY ( ) const

Definition at line 1055 of file QwCorrelatorNew.cc.

{
  QwMessage << Form("\nLinRegBevPeb::printSummaryY seen good eve=%lld  (CSV-format)",fGoodEventNumber)<<QwLog::endl;
  QwMessage << Form("  j,       mean,     sig(mean),   nSig(mean),  sig(distribution)    \n");
 
  for (int i = 0; i <nY; i++) {
    double meanI,sigI;
    if (getMeanY(i,meanI) < 0) QwWarning << "LRB::getMeanY failed" << QwLog::endl;
    if (getSigmaY(i,sigI) < 0) QwWarning << "LRB::getSigmaY failed" << QwLog::endl;
    double err = sigI / sqrt(fGoodEventNumber);
    double nSigErr = meanI / err;
    QwMessage << Form("Y%02d, %+11.4g, %12.4g, %8.1f, %12.4g "" ",i,meanI,err,nSigErr,sigI)<<QwLog::endl;
 
  }
}

References QwLog::endl(), fGoodEventNumber, getMeanY(), getSigmaY(), nY, QwMessage, and QwWarning.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ printSummaryYP()

void QwCorrelatorNew::printSummaryYP ( ) const

Definition at line 1095 of file QwCorrelatorNew.cc.

{
  QwMessage << Form("\nLinRegBevPeb::printSummaryYP seen good eve=%lld",fGoodEventNumber)<<QwLog::endl;
 
  if(fGoodEventNumber<2) { QwMessage<<"  too few events, skip"<<QwLog::endl; return;}
  QwMessage << Form("\n         name:             ");
  for (int i = 0; i <nP; i++) {
    QwMessage << Form(" %10sP%d "," ",i);
  }
  QwMessage << Form("\n  j                   meanY         sigY      correlation with Ps ....\n");
  for (int iy = 0; iy <nY; iy++) {
    double meanI,sigI;
    if (getMeanY(iy,meanI) < 0) QwWarning << "LRB::getMeanY failed" << QwLog::endl;
    if (getSigmaY(iy,sigI) < 0) QwWarning << "LRB::getSigmaY failed" << QwLog::endl;
 
    QwMessage << Form(" %3d %6sY%d:  %+12.4g  %12.4g ",iy," ",iy,meanI,sigI);
    for (int ip = 0; ip <nP; ip++) {
      double sigJ,cov;
      if (getSigmaP(ip,sigJ) < 0) QwWarning << "LRB::getSigmaP failed" << QwLog::endl;
      if (getCovariancePY(ip,iy,cov) < 0) QwWarning << "LRB::getCovariancePY failed" << QwLog::endl;
      double corel = cov / sigI / sigJ;
      QwMessage << Form("  %12.3g",corel);
    }
    QwMessage << Form("\n");
  }
}

References QwLog::endl(), fGoodEventNumber, getCovariancePY(), getMeanY(), getSigmaP(), getSigmaY(), nP, nY, QwMessage, and QwWarning.

Here is the call graph for this function:

◆ ProcessData()

void QwCorrelatorNew::ProcessData ( )

overridevirtual

Copy dependent variables to output variables (default processing).

Reimplemented from VQwDataHandler.

Definition at line 82 of file QwCorrelatorNew.cc.

{
  // Add to total count
  fTotalCount++;
 
  // Start as good event
  fGoodEvent = 0;
 
  // Event error flag
  fGoodEvent |= GetEventcutErrorFlag();
  if ( GetEventcutErrorFlag() != 0) fErrCounts_EF++;
  // Dependent variable error codes
  for (size_t i = 0; i < fDependentVar.size(); ++i) {
    fGoodEvent |= fDependentVar.at(i)->GetErrorCode();
    fDependentValues.at(i) = (fDependentVar[i]->GetValue(fBlock+1));
    if (fDependentVar.at(i)->GetErrorCode() !=0)  (fErrCounts_DV.at(i))++;
  }
  // Independent variable error codes
  for (size_t i = 0; i < fIndependentVar.size(); ++i) {
    fGoodEvent |= fIndependentVar.at(i)->GetErrorCode();
    fIndependentValues.at(i) = (fIndependentVar[i]->GetValue(fBlock+1));
    if (fIndependentVar.at(i)->GetErrorCode() !=0)  (fErrCounts_IV.at(i))++;
  }
 
  // If good, process event
  if (fGoodEvent == 0) {
    fGoodCount++;
 
    TVectorD P(fIndependentValues.size(), fIndependentValues.data());
    TVectorD Y(fDependentValues.size(),   fDependentValues.data());
    operator+= (std::make_pair(P, Y));
  }
}

References fBlock, VQwDataHandler::fDependentValues, VQwDataHandler::fDependentVar, fErrCounts_DV, fErrCounts_EF, fErrCounts_IV, fGoodCount, fGoodEvent, fIndependentValues, fIndependentVar, fTotalCount, VQwDataHandler::GetEventcutErrorFlag(), and operator+=().

Here is the call graph for this function:

◆ ProcessOptions()

void QwCorrelatorNew::ProcessOptions ( QwOptions & options )

Definition at line 62 of file QwCorrelatorNew.cc.

63{

64}

◆ setDims()

void QwCorrelatorNew::setDims	(	int	a,
		int	b )

inline

Definition at line 180 of file QwCorrelatorNew.h.

180{ nP = a; nY = b;}

References nP, and nY.

Referenced by ConnectChannels().

Here is the caller graph for this function:

◆ solve()

void QwCorrelatorNew::solve ( )

normalized covariances

Definition at line 1151 of file QwCorrelatorNew.cc.

{
  // off-diagonal raw covariance
  mVYP.Transpose(mVPY);
 
  // diagonal variances
  mVP = TMatrixDDiag(mVPP); mVP.Sqrt();
  mVY = TMatrixDDiag(mVYY); mVY.Sqrt();
 
  // correlation matrices
  mRPP = mVPP; mRPP.NormByColumn(mVP); mRPP.NormByRow(mVP);
  mRYY = mVYY; mRYY.NormByColumn(mVY); mRYY.NormByRow(mVY);
  mRPY = mVPY; mRPY.NormByColumn(mVP); mRPY.NormByRow(mVY);
 
  /// normalized covariances
  mSYY = mVYY * (1.0 / (fGoodEventNumber - 1.));
  mSPP = mVPP * (1.0 / (fGoodEventNumber - 1.));
  mSPY = mVPY * (1.0 / (fGoodEventNumber - 1.));
  mSYP.Transpose(mSPY);
 
  // uncertainties on the means
  mSP = TMatrixDDiag(mSPP); mSP.Sqrt();
  mSY = TMatrixDDiag(mSYY); mSY.Sqrt();
 
  // Warn if correlation matrix determinant close to zero (heuristic)
  if (mRPP.Determinant() < std::pow(10,-(2*nP))) {
    QwWarning << "LRB: correlation matrix nearly singular, "
              << "determinant = " << mRPP.Determinant()
              << " (set includes highly correlated variable pairs)"
              << QwLog::endl;
    if (fGoodEventNumber > 10*nP) {
      QwMessage << fGoodEventNumber << " events" << QwLog::endl;
      QwMessage << "Covariance matrix: " << QwLog::endl; mVPP.Print();
      QwMessage << "Correlation matrix: " << QwLog::endl; mRPP.Print();
    }
    QwWarning << "LRB: solving failed (this happens when only few events)."
              << QwLog::endl;
    return;
  }
  // slopes
  TMatrixD invRPP(TMatrixD::kInverted, mRPP);
  Axy = TMatrixD(invRPP, TMatrixD::kMult, mRPY);
  Axy.NormByColumn(mSP); // divide
  Axy.NormByRow(mSY, ""); // mult
  Ayx.Transpose(Axy);
 
  // new means
  mMYp = mMY - Ayx * mMP;
 
  // new raw covariance
  mVYYp = mVYY + Ayx * mVPP * Axy - (Ayx * mVPY + mVYP * Axy);
  // new variances
  mVYp = TMatrixDDiag(mVYYp); mVYp.Sqrt();
 
  // new normalized covariance
  mSYYp = mSYY + Ayx * mSPP * Axy - (Ayx * mSPY + mSYP * Axy);
  // uncertainties on the new means
  mSYp = TMatrixDDiag(mSYYp); mSYp.Sqrt();
 
  // new correlation matrix
  mRYYp = mVYYp; mRYYp.NormByColumn(mVYp); mRYYp.NormByRow(mVYp);
 
  // slope uncertainties
  double norm = 1. / (fGoodEventNumber - nP - 1);
  dAxy.Zero();
  dAxy.Rank1Update(TMatrixDDiag(invRPP), TMatrixDDiag(mRYYp), norm); // diag mRYYp = row of ones
  dAxy.Sqrt();
  dAxy.NormByColumn(mSP); // divide
  dAxy.NormByRow(mSYp, ""); // mult
  dAyx.Transpose(dAxy);
 
  fErrorFlag = 0;
}

References Axy, Ayx, dAxy, dAyx, QwLog::endl(), fErrorFlag, fGoodEventNumber, mMP, mMY, mMYp, mRPP, mRPY, mRYY, mRYYp, mSP, mSPP, mSPY, mSY, mSYP, mSYp, mSYY, mSYYp, mVP, mVPP, mVPY, mVY, mVYP, mVYp, mVYY, mVYYp, nP, QwMessage, and QwWarning.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ WriteAliasFile()

void QwCorrelatorNew::WriteAliasFile ( )

protected

Definition at line 668 of file QwCorrelatorNew.cc.

{
  // Ensure output file is open
  if (fAliasOutputFile.bad()) {
    QwWarning << "QwCorrelatorNew: Could not write to alias output file " << QwLog::endl;
    return;
  }
 
  fAliasOutputFile << " if (i == " << fCycleCounter << ") {" << std::endl;
  fAliasOutputFile << Form("  TTree* tree = (TTree*) gDirectory->Get(\"mul\");") << std::endl;
  for (int i = 0; i < nY; i++) {
    fAliasOutputFile << Form("  tree->SetAlias(\"reg_%s\",",fDependentFull[i].c_str()) << std::endl;
    fAliasOutputFile << Form("         \"%s",fDependentFull[i].c_str());
    for (int j = 0; j < nP; j++) {
      fAliasOutputFile << Form("%+.4e*%s", -this->Axy(j,i), fIndependentFull[j].c_str());
    }
    fAliasOutputFile << "\");" << std::endl;
  }
  fAliasOutputFile << " }" << std::endl;
 
  // Increment call counter
  fCycleCounter++;
}

References Axy, QwLog::endl(), fAliasOutputFile, fCycleCounter, VQwDataHandler::fDependentFull, fIndependentFull, nP, nY, and QwWarning.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ WriteAlphaFile()

void QwCorrelatorNew::WriteAlphaFile ( )

protected

Definition at line 550 of file QwCorrelatorNew.cc.

{
  // Ensure in output file
  if (fAlphaOutputFile) fAlphaOutputFile->cd();
 
  // Write objects
  this->Axy.Write("slopes");
  this->dAxy.Write("sigSlopes");
 
  this->mRPP.Write("IV_IV_correlation");
  this->mRPY.Write("IV_DV_correlation");
  this->mRYY.Write("DV_DV_correlation");
  this->mRYYp.Write("DV_DV_correlation_prime");
 
  this->mMP.Write("IV_mean");
  this->mMY.Write("DV_mean");
  this->mMYp.Write("DV_mean_prime");
 
  // number of events
  TMatrixD Mstat(1,1);
  Mstat(0,0)=this->getUsedEve();
  Mstat.Write("MyStat");
 
  //... IVs
  TH1D hiv("IVname","names of IVs",nP,-0.5,nP-0.5);
  for (int i=0;i<nP;i++) hiv.Fill(fIndependentFull[i].c_str(),i);
  hiv.Write();
 
  //... DVs
  TH1D hdv("DVname","names of IVs",nY,-0.5,nY-0.5);
  for (int i=0;i<nY;i++) hdv.Fill(fDependentFull[i].c_str(),i);
  hdv.Write();
 
  // sigmas
  this->mSP.Write("IV_sigma");
  this->mSY.Write("DV_sigma");
  this->mSYp.Write("DV_sigma_prime");
 
  // raw covariances
  this->mVPP.Write("IV_IV_rawVariance");
  this->mVPY.Write("IV_DV_rawVariance");
  this->mVYY.Write("DV_DV_rawVariance");
  this->mVYYp.Write("DV_DV_rawVariance_prime");
  TVectorD mVY2(TMatrixDDiag(this->mVYY));
  mVY2.Write("DV_rawVariance");
  TVectorD mVP2(TMatrixDDiag(this->mVPP));
  mVP2.Write("IV_rawVariance");
  TVectorD mVY2prime(TMatrixDDiag(this->mVYYp));
  mVY2prime.Write("DV_rawVariance_prime");
 
  // normalized covariances
  this->mSPP.Write("IV_IV_normVariance");
  this->mSPY.Write("IV_DV_normVariance");
  this->mSYY.Write("DV_DV_normVariance");
  this->mSYYp.Write("DV_DV_normVariance_prime");
  TVectorD sigY2(TMatrixDDiag(this->mSYY));
  sigY2.Write("DV_normVariance");
  TVectorD sigX2(TMatrixDDiag(this->mSPP));
  sigX2.Write("IV_normVariance");
  TVectorD sigY2prime(TMatrixDDiag(this->mSYYp));
  sigY2prime.Write("DV_normVariance_prime");
 
  this->Axy.Write("A_xy");
  this->Ayx.Write("A_yx");
}

References Axy, Ayx, dAxy, fAlphaOutputFile, VQwDataHandler::fDependentFull, fIndependentFull, getUsedEve(), mMP, mMY, mMYp, mRPP, mRPY, mRYY, mRYYp, mSP, mSPP, mSPY, mSY, mSYp, mSYY, mSYYp, mVPP, mVPY, mVYY, mVYYp, nP, and nY.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

Friends And Related Symbol Documentation

◆ operator+

QwCorrelatorNew operator+	(	QwCorrelatorNew	lhs,
		const QwCorrelatorNew &	rhs )

friend

Definition at line 205 of file QwCorrelatorNew.h.

  {
    lhs += rhs; // reuse compound assignment
    return lhs; // return the result by value (uses move constructor)
  }

References QwCorrelatorNew().

◆ operator<<

std::ostream & operator<<	(	std::ostream &	stream,
		const QwCorrelatorNew &	h )

friend

Output stream operator.

Definition at line 217 of file QwCorrelatorNew.h.

{
  stream << "LRB: " << h.fGoodEventNumber << " events";
  return stream;
}

References fGoodEventNumber, and QwCorrelatorNew().

Field Documentation

◆ Axy

TMatrixD QwCorrelatorNew::Axy

private

slopes

Definition at line 162 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), printSummaryAlphas(), solve(), WriteAliasFile(), and WriteAlphaFile().

◆ Ayx

TMatrixD QwCorrelatorNew::Ayx

private

Definition at line 162 of file QwCorrelatorNew.h.

Referenced by clear(), init(), solve(), and WriteAlphaFile().

◆ dAxy

TMatrixD QwCorrelatorNew::dAxy

private

Definition at line 162 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), printSummaryAlphas(), solve(), and WriteAlphaFile().

◆ dAyx

TMatrixD QwCorrelatorNew::dAyx

private

Definition at line 162 of file QwCorrelatorNew.h.

Referenced by clear(), init(), and solve().

◆ fAliasOutputFile

std::ofstream QwCorrelatorNew::fAliasOutputFile

protected

Definition at line 100 of file QwCorrelatorNew.h.

Referenced by CloseAliasFile(), OpenAliasFile(), and WriteAliasFile().

◆ fAliasOutputFileBase

std::string QwCorrelatorNew::fAliasOutputFileBase

protected

Definition at line 97 of file QwCorrelatorNew.h.

Referenced by OpenAliasFile(), ParseConfigFile(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fAliasOutputFileSuff

std::string QwCorrelatorNew::fAliasOutputFileSuff

protected

Definition at line 98 of file QwCorrelatorNew.h.

Referenced by OpenAliasFile(), ParseConfigFile(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fAliasOutputPath

std::string QwCorrelatorNew::fAliasOutputPath

protected

Definition at line 99 of file QwCorrelatorNew.h.

Referenced by OpenAliasFile(), ParseConfigFile(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fAlphaOutputFile

TFile* QwCorrelatorNew::fAlphaOutputFile

protected

Definition at line 90 of file QwCorrelatorNew.h.

Referenced by CloseAlphaFile(), OpenAlphaFile(), QwCorrelatorNew(), QwCorrelatorNew(), and WriteAlphaFile().

◆ fAlphaOutputFileBase

std::string QwCorrelatorNew::fAlphaOutputFileBase

protected

Definition at line 87 of file QwCorrelatorNew.h.

Referenced by OpenAlphaFile(), ParseConfigFile(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fAlphaOutputFileSuff

std::string QwCorrelatorNew::fAlphaOutputFileSuff

protected

Definition at line 88 of file QwCorrelatorNew.h.

Referenced by OpenAlphaFile(), ParseConfigFile(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fAlphaOutputPath

std::string QwCorrelatorNew::fAlphaOutputPath

protected

Definition at line 89 of file QwCorrelatorNew.h.

Referenced by OpenAlphaFile(), ParseConfigFile(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fBlock

Int_t QwCorrelatorNew::fBlock

protected

Definition at line 74 of file QwCorrelatorNew.h.

Referenced by ParseConfigFile(), ProcessData(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fCycleCounter

Int_t QwCorrelatorNew::fCycleCounter

private

Definition at line 126 of file QwCorrelatorNew.h.

Referenced by QwCorrelatorNew(), QwCorrelatorNew(), and WriteAliasFile().

◆ fDisableHistos

bool QwCorrelatorNew::fDisableHistos

protected

Definition at line 76 of file QwCorrelatorNew.h.

Referenced by ConstructHistograms(), FillHistograms(), ParseConfigFile(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ fErrCounts_DV

std::vector<int> QwCorrelatorNew::fErrCounts_DV

protected

Definition at line 110 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ClearEventData(), ConnectChannels(), and ProcessData().

◆ fErrCounts_EF

int QwCorrelatorNew::fErrCounts_EF

protected

Definition at line 108 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ClearEventData(), and ProcessData().

◆ fErrCounts_IV

std::vector<int> QwCorrelatorNew::fErrCounts_IV

protected

Definition at line 109 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ClearEventData(), ConnectChannels(), and ProcessData().

◆ fErrorFlag

Int_t QwCorrelatorNew::fErrorFlag

private

is information valid

Definition at line 133 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), QwCorrelatorNew(), QwCorrelatorNew(), and solve().

◆ fGoodCount

int QwCorrelatorNew::fGoodCount

protected

Definition at line 106 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ClearEventData(), ConstructTreeBranches(), and ProcessData().

◆ fGoodEvent

unsigned int QwCorrelatorNew::fGoodEvent

protected

Definition at line 112 of file QwCorrelatorNew.h.

Referenced by ClearEventData(), FillHistograms(), and ProcessData().

◆ fGoodEventNumber

Long64_t QwCorrelatorNew::fGoodEventNumber

private

accumulated so far

Definition at line 134 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), failed(), getCovarianceP(), getCovariancePY(), getCovarianceY(), getMeanP(), getMeanY(), getMeanYprime(), getSigmaP(), getSigmaY(), getSigmaYprime(), getUsedEve(), init(), operator+=(), operator+=(), operator<<, printSummaryAlphas(), printSummaryP(), printSummaryY(), printSummaryYP(), QwCorrelatorNew(), QwCorrelatorNew(), and solve().

◆ fH1dv

std::vector<TH1D> QwCorrelatorNew::fH1dv

private

Definition at line 122 of file QwCorrelatorNew.h.

Referenced by ConstructHistograms(), and FillHistograms().

◆ fH1iv

std::vector<TH1D> QwCorrelatorNew::fH1iv

private

Definition at line 121 of file QwCorrelatorNew.h.

Referenced by ConstructHistograms(), and FillHistograms().

◆ fH2dv

std::vector<std::vector<TH2D> > QwCorrelatorNew::fH2dv

private

Definition at line 124 of file QwCorrelatorNew.h.

Referenced by ConstructHistograms(), and FillHistograms().

◆ fH2iv

std::vector<std::vector<TH2D> > QwCorrelatorNew::fH2iv

private

Definition at line 123 of file QwCorrelatorNew.h.

Referenced by ConstructHistograms(), and FillHistograms().

◆ fHnames

std::vector<TH1D> QwCorrelatorNew::fHnames

private

Definition at line 120 of file QwCorrelatorNew.h.

Referenced by ConstructHistograms().

◆ fIndependentFull

std::vector< std::string > QwCorrelatorNew::fIndependentFull

protected

Definition at line 78 of file QwCorrelatorNew.h.

Referenced by ConnectChannels(), LoadChannelMap(), WriteAliasFile(), and WriteAlphaFile().

◆ fIndependentName

std::vector< std::string > QwCorrelatorNew::fIndependentName

protected

Definition at line 82 of file QwCorrelatorNew.h.

Referenced by ConnectChannels(), ConstructHistograms(), and LoadChannelMap().

◆ fIndependentType

std::vector< EQwHandleType > QwCorrelatorNew::fIndependentType

protected

Definition at line 81 of file QwCorrelatorNew.h.

Referenced by ConnectChannels(), and LoadChannelMap().

◆ fIndependentValues

std::vector< Double_t > QwCorrelatorNew::fIndependentValues

protected

Definition at line 85 of file QwCorrelatorNew.h.

Referenced by ConnectChannels(), FillHistograms(), and ProcessData().

◆ fIndependentVar

std::vector< const VQwHardwareChannel* > QwCorrelatorNew::fIndependentVar

protected

Definition at line 84 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ConnectChannels(), and ProcessData().

◆ fNameNoSpaces

TString QwCorrelatorNew::fNameNoSpaces

private

Definition at line 116 of file QwCorrelatorNew.h.

Referenced by QwCorrelatorNew().

◆ fPrintCorrelations

bool QwCorrelatorNew::fPrintCorrelations = false

staticprivate

Definition at line 35 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), and DefineOptions().

◆ fTotalCount

int QwCorrelatorNew::fTotalCount

protected

Definition at line 105 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ClearEventData(), ConstructTreeBranches(), and ProcessData().

◆ fTree

TTree* QwCorrelatorNew::fTree

protected

Definition at line 95 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ConstructTreeBranches(), QwCorrelatorNew(), and QwCorrelatorNew().

◆ mMP

TVectorD QwCorrelatorNew::mMP

private

mean values

Definition at line 158 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), getMeanP(), init(), operator+=(), operator+=(), print(), solve(), and WriteAlphaFile().

◆ mMY

TVectorD QwCorrelatorNew::mMY

private

Definition at line 158 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), getMeanY(), init(), operator+=(), operator+=(), print(), printSummaryMeansWithUnc(), solve(), and WriteAlphaFile().

◆ mMYp

TVectorD QwCorrelatorNew::mMYp

private

Definition at line 158 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), getMeanYprime(), init(), printSummaryMeansWithUncCorrected(), solve(), and WriteAlphaFile().

◆ mRPP

TMatrixD QwCorrelatorNew::mRPP

private

Definition at line 138 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mRPY

TMatrixD QwCorrelatorNew::mRPY

private

correlations

Definition at line 137 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mRYP

TMatrixD QwCorrelatorNew::mRYP

private

Definition at line 137 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), and init().

◆ mRYY

TMatrixD QwCorrelatorNew::mRYY

private

Definition at line 138 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mRYYp

TMatrixD QwCorrelatorNew::mRYYp

private

Definition at line 139 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mSP

TVectorD QwCorrelatorNew::mSP

private

sigmas

Definition at line 154 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mSPP

TMatrixD QwCorrelatorNew::mSPP

private

Definition at line 151 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mSPY

TMatrixD QwCorrelatorNew::mSPY

private

normalized covariances

Definition at line 150 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mSY

TVectorD QwCorrelatorNew::mSY

private

Definition at line 154 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), printSummaryMeansWithUnc(), solve(), and WriteAlphaFile().

◆ mSYP

TMatrixD QwCorrelatorNew::mSYP

private

Definition at line 150 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), and solve().

◆ mSYp

TVectorD QwCorrelatorNew::mSYp

private

Definition at line 155 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), printSummaryMeansWithUncCorrected(), solve(), and WriteAlphaFile().

◆ mSYY

TMatrixD QwCorrelatorNew::mSYY

private

Definition at line 151 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mSYYp

TMatrixD QwCorrelatorNew::mSYYp

private

Definition at line 152 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), solve(), and WriteAlphaFile().

◆ mVP

TVectorD QwCorrelatorNew::mVP

private

variances

Definition at line 146 of file QwCorrelatorNew.h.

Referenced by clear(), init(), and solve().

◆ mVPP

TMatrixD QwCorrelatorNew::mVPP

private

Definition at line 143 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), getCovarianceP(), getSigmaP(), init(), operator+=(), operator+=(), print(), solve(), and WriteAlphaFile().

◆ mVPY

TMatrixD QwCorrelatorNew::mVPY

private

unnormalized covariances

Definition at line 142 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), getCovariancePY(), init(), operator+=(), operator+=(), print(), solve(), and WriteAlphaFile().

◆ mVY

TVectorD QwCorrelatorNew::mVY

private

Definition at line 146 of file QwCorrelatorNew.h.

Referenced by clear(), init(), and solve().

◆ mVYP

TMatrixD QwCorrelatorNew::mVYP

private

Definition at line 142 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), init(), and solve().

◆ mVYp

TVectorD QwCorrelatorNew::mVYp

private

Definition at line 147 of file QwCorrelatorNew.h.

Referenced by clear(), init(), and solve().

◆ mVYY

TMatrixD QwCorrelatorNew::mVYY

private

Definition at line 143 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), getCovarianceY(), getSigmaY(), init(), operator+=(), operator+=(), print(), solve(), and WriteAlphaFile().

◆ mVYYp

TMatrixD QwCorrelatorNew::mVYYp

private

Definition at line 144 of file QwCorrelatorNew.h.

Referenced by clear(), ConstructTreeBranches(), getSigmaYprime(), init(), print(), solve(), and WriteAlphaFile().

◆ nP

int QwCorrelatorNew::nP

private

Definition at line 117 of file QwCorrelatorNew.h.

Referenced by CalcCorrelations(), ConnectChannels(), ConstructHistograms(), ConstructTreeBranches(), failed(), FillHistograms(), getCovarianceP(), getCovariancePY(), getMeanP(), getSigmaP(), init(), print(), printSummaryAlphas(), printSummaryP(), printSummaryYP(), QwCorrelatorNew(), QwCorrelatorNew(), QwCorrelatorNew(), setDims(), solve(), WriteAliasFile(), and WriteAlphaFile().

◆ nY

int QwCorrelatorNew::nY

private

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ QwCorrelatorNew() [1/3]

◆ QwCorrelatorNew() [2/3]

◆ ~QwCorrelatorNew()

◆ QwCorrelatorNew() [3/3]

Member Function Documentation

◆ AccumulateRunningSum()

◆ CalcCorrelations()

◆ clear()

◆ ClearEventData()

◆ CloseAliasFile()

◆ CloseAlphaFile()

◆ ConnectChannels()

◆ ConstructHistograms()

◆ ConstructTreeBranches()

◆ DefineOptions()

◆ failed()

◆ FillHistograms()

◆ FillTreeBranches()

◆ FinishDataHandler()

◆ getCovarianceP()

◆ getCovariancePY()

◆ getCovarianceY()

◆ getMeanP()

◆ getMeanY()

◆ getMeanYprime()

◆ getSigmaP()

◆ getSigmaY()

◆ getSigmaYprime()

◆ getUsedEve()

◆ init()

◆ LoadChannelMap()

◆ OpenAliasFile()

◆ OpenAlphaFile()

◆ operator+=() [1/2]

◆ operator+=() [2/2]

◆ ParseConfigFile()

◆ print()

◆ printSummaryAlphas()

◆ printSummaryMeansWithUnc()

◆ printSummaryMeansWithUncCorrected()

◆ printSummaryP()

◆ printSummaryY()

◆ printSummaryYP()

◆ ProcessData()

◆ ProcessOptions()

◆ setDims()

◆ solve()

◆ WriteAliasFile()

◆ WriteAlphaFile()

Friends And Related Symbol Documentation

◆ operator+

◆ operator<<

Field Documentation

◆ Axy

◆ Ayx

◆ dAxy

◆ dAyx

◆ fAliasOutputFile

◆ fAliasOutputFileBase

◆ fAliasOutputFileSuff

◆ fAliasOutputPath

◆ fAlphaOutputFile

◆ fAlphaOutputFileBase

◆ fAlphaOutputFileSuff

◆ fAlphaOutputPath

◆ fBlock

◆ fCycleCounter

◆ fDisableHistos

◆ fErrCounts_DV