GrandCorrelator Class Reference

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

#include <GrandCorrelator.h>

Inheritance diagram for GrandCorrelator:

Collaboration diagram for GrandCorrelator:

Public Member Functions
	GrandCorrelator (const TString &name)
	Constructor with name.
	GrandCorrelator (const GrandCorrelator &name)
	~GrandCorrelator () 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
GrandCorrelator &	operator+= (const std::pair< TVectorD, TVectorD > &rhs)
GrandCorrelator &	operator+= (const GrandCorrelator &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, GrandCorrelator >
	~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, GrandCorrelator >
static VQwDataHandler *	Create (const std::string &name)
	Object creation.
static GrandCorrelator *	Cast (GrandCorrelator *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
std::vector< int >	fErrCounts_ALL
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
	GrandCorrelator ()

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	mNij
	new matrices and vectors
TMatrixD	mSij
TMatrixD	mMij
TMatrixD	mCij
TMatrixD	mVij
TMatrixD	mRij
TMatrixD	mVFULL
TMatrixD	mRFULL
TMatrixD	mSFULL
TMatrixD	sigma_ij
TMatrixD	sigma_ji
TMatrixD	mVFULL_clean
TMatrixD	mSFULL_clean
std::vector< const VQwHardwareChannel * >	fAllVar
std::vector< double >	fAllValues
std::vector< bool >	fAllGood
std::vector< std::string >	fAllName
TMatrixD	Axy
	slopes
TMatrixD	Ayx
TMatrixD	dAxy
TMatrixD	dAyx

Static Private Attributes
static bool	fPrintCorrelations = false

Friends
GrandCorrelator	operator+ (GrandCorrelator lhs, const GrandCorrelator &rhs)
std::ostream &	operator<< (std::ostream &stream, const GrandCorrelator &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 GrandCorrelator.h.

Constructor & Destructor Documentation

GrandCorrelator() [1/3]

GrandCorrelator::GrandCorrelator

(

const TString &

name

)

Constructor with name.

Definition at line 28 of file GrandCorrelator.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(), GrandCorrelator(), operator+, operator+=(), operator+=(), and operator<<.

Here is the call graph for this function:

Here is the caller graph for this function:

GrandCorrelator() [2/3]

GrandCorrelator::GrandCorrelator

(

const GrandCorrelator &

name

)

Definition at line 751 of file GrandCorrelator.cc.

: VQwDataHandler(source),
  fBlock(source.fBlock),
  fDisableHistos(source.fDisableHistos),
  fAlphaOutputFileBase(source.fAlphaOutputFileBase),
  fAlphaOutputFileSuff(source.fAlphaOutputFileSuff),
  fAlphaOutputPath(source.fAlphaOutputPath),
  fAlphaOutputFile(nullptr),
  fTree(nullptr),
  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 << "GrandCorrelator 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, GrandCorrelator(), nP, nY, QwMessage, QwWarning, and VQwDataHandler::VQwDataHandler().

Here is the call graph for this function:

~GrandCorrelator()

GrandCorrelator::~GrandCorrelator ( )

override

Definition at line 55 of file GrandCorrelator.cc.

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

References CloseAliasFile(), and CloseAlphaFile().

Here is the call graph for this function:

GrandCorrelator() [3/3]

GrandCorrelator::GrandCorrelator ( )

private

Definition at line 743 of file GrandCorrelator.cc.

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

References fErrorFlag, fGoodEventNumber, nP, and nY.

Member Function Documentation

AccumulateRunningSum()

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

overridevirtual

Reimplemented from VQwDataHandler.

Definition at line 169 of file GrandCorrelator.cc.

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

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

Here is the call graph for this function:

CalcCorrelations()

void GrandCorrelator::CalcCorrelations

(

)

Definition at line 181 of file GrandCorrelator.cc.

{
  // Check if any channels are active
  if (nP == 0 || nY == 0) {
    return;
  }
 
  QwMessage << "GrandCorrelator::CalcCorrelations(): name=" << GetName() << QwLog::endl;
 
  // Print entry summary
  QwVerbose << "GrandCorrelator: "
            << "total entries: " << fTotalCount << ", "
            << "good entries: " << fGoodEventNumber
            << QwLog::endl;
  // and warn if zero
  if (fTotalCount > 100 && fGoodEventNumber == 0) {
    QwWarning << "GrandCorrelator: "
              << "< 1% good events, "
              << fGoodEventNumber << " 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 (kTRUE || 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, fGoodEventNumber, 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 GrandCorrelator::clear

(

)

Definition at line 833 of file GrandCorrelator.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();
 
  mNij.Zero();
  mSij.Zero();
  mMij.Zero();
  mCij.Zero();
  mVij.Zero();
  mRij.Zero();
  sigma_ij.Zero();
  sigma_ji.Zero();
  mVFULL.Zero();
  mRFULL.Zero();
  mSFULL.Zero();
  mVFULL_clean.Zero();
  mSFULL_clean.Zero();
 
  fErrorFlag = -1;
  fGoodEventNumber = 0;
}

References Axy, Ayx, dAxy, dAyx, fErrorFlag, fGoodEventNumber, mCij, mMij, mMP, mMY, mMYp, mNij, mRFULL, mRij, mRPP, mRPY, mRYP, mRYY, mRYYp, mSFULL, mSFULL_clean, mSij, mSP, mSPP, mSPY, mSY, mSYP, mSYp, mSYY, mSYYp, mVFULL, mVFULL_clean, mVij, mVP, mVPP, mVPY, mVY, mVYP, mVYp, mVYY, mVYYp, sigma_ij, and sigma_ji.

Referenced by ClearEventData().

Here is the caller graph for this function:

ClearEventData()

void GrandCorrelator::ClearEventData ( )

overridevirtual

Reimplemented from VQwDataHandler.

Definition at line 153 of file GrandCorrelator.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;
  fGoodEventNumber = 0;
  fGoodEvent = -1;
 
  // Clear regression
  this->clear();
}

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

Referenced by GrandCorrelator(), and GrandCorrelator().

Here is the call graph for this function:

Here is the caller graph for this function:

CloseAliasFile()

void GrandCorrelator::CloseAliasFile ( )

protected

Definition at line 705 of file GrandCorrelator.cc.

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

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

Referenced by ~GrandCorrelator().

Here is the call graph for this function:

Here is the caller graph for this function:

CloseAlphaFile()

void GrandCorrelator::CloseAlphaFile ( )

protected

Definition at line 696 of file GrandCorrelator.cc.

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

References fAlphaOutputFile.

Referenced by ~GrandCorrelator().

Here is the caller graph for this function:

ConnectChannels()

Int_t GrandCorrelator::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 293 of file GrandCorrelator.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 << "GrandCorrelator::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);
 
  // Create vector list for all values
  fAllVar = fIndependentVar;
  fAllVar.insert(fAllVar.end(), fDependentVar.begin(), fDependentVar.end());
  fAllValues = fIndependentValues;
  fAllValues.insert(fAllValues.end(), fDependentValues.begin(), fDependentValues.end());
  fAllGood.resize(fAllValues.size(), true);
 
 
  return 0;
}

References QwLog::endl(), fAllGood, fAllValues, fAllVar, 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 GrandCorrelator::ConstructHistograms ( TDirectory * folder, TString & prefix )

overridevirtual

Construct the histograms in a folder with a prefix.

Reimplemented from VQwDataHandler.

Definition at line 488 of file GrandCorrelator.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 GrandCorrelator::ConstructTreeBranches ( QwRootFile * treerootfile, const std::string & treeprefix = "", const std::string & branchprefix = "" )

overridevirtual

Construct the tree branches.

Reimplemented from VQwDataHandler.

Definition at line 400 of file GrandCorrelator.cc.

{
  // Check if any channels are active
  if (nP == 0 || nY == 0) {
    return;
  }
 
  // Check if tree name is specified
  if (fTreeName == "") {
    QwWarning << "GrandCorrelator: 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"),  &fGoodEventNumber);
 
  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, 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 GrandCorrelator::DefineOptions ( QwOptions & options )

static

Definition at line 62 of file GrandCorrelator.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 GrandCorrelator::failed ( )

inline

Definition at line 182 of file GrandCorrelator.h.

{ return fGoodEventNumber < nP + 1; }

References fGoodEventNumber, and nP.

Referenced by CalcCorrelations().

Here is the caller graph for this function:

FillHistograms()

void GrandCorrelator::FillHistograms ( )

overridevirtual

Fill the histograms.

Reimplemented from VQwDataHandler.

Definition at line 572 of file GrandCorrelator.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 GrandCorrelator::FillTreeBranches ( QwRootFile * treerootfile )

inlineoverridevirtual

Fill the tree branches.

Reimplemented from VQwDataHandler.

Definition at line 62 of file GrandCorrelator.h.

{ };

FinishDataHandler()

void GrandCorrelator::FinishDataHandler ( )

inlineoverridevirtual

Reimplemented from VQwDataHandler.

Definition at line 51 of file GrandCorrelator.h.

                                   {
    CalcCorrelations();
  }

References CalcCorrelations().

Here is the call graph for this function:

getCovarianceP()

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

Get mean value of a variable, returns error code.

Definition at line 1049 of file GrandCorrelator.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 GrandCorrelator::getCovariancePY	(	int	ip,
		int	iy,
		Double_t &	covar ) const

Definition at line 1062 of file GrandCorrelator.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 GrandCorrelator::getCovarianceY	(	int	i,
		int	j,
		Double_t &	covar ) const

Definition at line 1075 of file GrandCorrelator.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 GrandCorrelator::getMeanP	(	const int	i,
		Double_t &	mean ) const

Get mean value of a variable, returns error code.

Definition at line 982 of file GrandCorrelator.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 GrandCorrelator::getMeanY	(	const int	i,
		Double_t &	mean ) const

Definition at line 993 of file GrandCorrelator.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 GrandCorrelator::getMeanYprime	(	const int	i,
		Double_t &	mean ) const

Definition at line 1004 of file GrandCorrelator.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 GrandCorrelator::getSigmaP	(	const int	i,
		Double_t &	sigma ) const

Get mean value of a variable, returns error code.

Definition at line 1015 of file GrandCorrelator.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 GrandCorrelator::getSigmaY	(	const int	i,
		Double_t &	sigma ) const

Definition at line 1027 of file GrandCorrelator.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 GrandCorrelator::getSigmaYprime	(	const int	i,
		Double_t &	sigma ) const

Definition at line 1038 of file GrandCorrelator.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 GrandCorrelator::getUsedEve ( ) const

inline

Definition at line 212 of file GrandCorrelator.h.

{ return fGoodEventNumber; };

References fGoodEventNumber.

Referenced by WriteAlphaFile().

Here is the caller graph for this function:

init()

void GrandCorrelator::init

(

)

Definition at line 778 of file GrandCorrelator.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);
 
  mNij.ResizeTo(nP+nY,nP+nY);
  mSij.ResizeTo(nP+nY,nP+nY);
  mMij.ResizeTo(nP+nY,nP+nY);
  mCij.ResizeTo(nP+nY,nP+nY);
  mVij.ResizeTo(nP+nY,nP+nY);
  mRij.ResizeTo(nP+nY,nP+nY);
  sigma_ij.ResizeTo(nP+nY,nP+nY);
  sigma_ji.ResizeTo(nP+nY,nP+nY);
  mVFULL.ResizeTo(nP+nY,nP+nY);
  mRFULL.ResizeTo(nP+nY,nP+nY);
  mSFULL.ResizeTo(nP+nY,nP+nY);
  mVFULL_clean.ResizeTo(nP+nY,nP+nY);
  mSFULL_clean.ResizeTo(nP+nY,nP+nY);
 
  fGoodEventNumber = 0;
}

References Axy, Ayx, dAxy, dAyx, fGoodEventNumber, mCij, mMij, mMP, mMY, mMYp, mNij, mRFULL, mRij, mRPP, mRPY, mRYP, mRYY, mRYYp, mSFULL, mSFULL_clean, mSij, mSP, mSPP, mSPY, mSY, mSYP, mSYp, mSYY, mSYYp, mVFULL, mVFULL_clean, mVij, mVP, mVPP, mVPY, mVY, mVYP, mVYp, mVYY, mVYYp, nP, nY, sigma_ij, and sigma_ji.

Referenced by ConnectChannels().

Here is the caller graph for this function:

LoadChannelMap()

Int_t GrandCorrelator::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 252 of file GrandCorrelator.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 GrandCorrelator 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 GrandCorrelator::OpenAliasFile ( const std::string & prefix )

protected

Definition at line 678 of file GrandCorrelator.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 << "GrandCorrelator: 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 GrandCorrelator::OpenAlphaFile ( const std::string & prefix )

protected

Definition at line 664 of file GrandCorrelator.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 << "GrandCorrelator 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]

GrandCorrelator & GrandCorrelator::operator+=

(

const GrandCorrelator &

rhs

)

Definition at line 941 of file GrandCorrelator.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, GrandCorrelator(), mMP, mMY, mVPP, mVPY, and mVYY.

Here is the call graph for this function:

operator+=() [2/2]

GrandCorrelator & GrandCorrelator::operator+=

(

const std::pair< TVectorD, TVectorD > &

rhs

)

Definition at line 903 of file GrandCorrelator.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, GrandCorrelator(), mMP, mMY, mVPP, mVPY, and mVYY.

Referenced by AccumulateRunningSum().

Here is the call graph for this function:

Here is the caller graph for this function:

ParseConfigFile()

void GrandCorrelator::ParseConfigFile ( QwParameterFile & file )

overridevirtual

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

Reimplemented from VQwDataHandler.

Definition at line 73 of file GrandCorrelator.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 << "GrandCorrelator: 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 GrandCorrelator::print

(

)

Definition at line 889 of file GrandCorrelator.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 GrandCorrelator::printSummaryAlphas

(

)

const

Definition at line 1144 of file GrandCorrelator.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 GrandCorrelator::printSummaryMeansWithUnc

(

)

const

Definition at line 1194 of file GrandCorrelator.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 GrandCorrelator::printSummaryMeansWithUncCorrected

(

)

const

Definition at line 1207 of file GrandCorrelator.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 GrandCorrelator::printSummaryP

(

)

const

Definition at line 1088 of file GrandCorrelator.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 GrandCorrelator::printSummaryY

(

)

const

Definition at line 1124 of file GrandCorrelator.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 GrandCorrelator::printSummaryYP

(

)

const

Definition at line 1164 of file GrandCorrelator.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 GrandCorrelator::ProcessData ( )

overridevirtual

Copy dependent variables to output variables (default processing).

Reimplemented from VQwDataHandler.

Definition at line 90 of file GrandCorrelator.cc.

{
  fTotalCount++;
  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))++;
  }
 
  for(size_t i = 0; i < fAllVar.size(); ++i){
    fAllGood[i] = (fAllVar[i]->GetErrorCode() == 0);
    fAllValues[i] = fAllVar[i]->GetValue(fBlock+1);
    //if(!fAllGood[i]) fErrCounts_IV[i]++;
  }
  if(fGoodEvent != 0){
    return;
  }
  fGoodEventNumber++;
  //TVectorD P(fIndependentValues.size(), fIndependentValues.data());
  //TVectorD Y(fDependentValues.size(),   fDependentValues.data());
  //operator+= (std::make_pair(P, Y));
 
  for(size_t i = 0; i<fAllVar.size(); ++i){
    for(size_t j = i; j<fAllVar.size(); ++j){
        if(!fAllGood[i] || !fAllGood[j]) continue;
 
        double xi = fAllValues[i];
        double xj = fAllValues[j];
        double delta_i = xi - mMij(i,j);
        double delta_j = xj - mMij(j,i);
 
        mNij(i,j) += 1;
        mMij(i,j) += delta_i / mNij(i,j);
        mSij(i,j) += delta_i * (xi-mMij(i,j));
        if (mNij(i,j)<=1){
          mCij(i,j) = 0.0;
        } else {
          mCij(i,j) += ((mNij(i,j) - 1) / mNij(i,j))  * (delta_i) * (delta_j);
        }
 
        if (j>i){
          mNij(j,i) = mNij(i,j);
          mMij(j,i) += delta_j / mNij(j,i);
          mSij(j,i) += delta_j * (xj - mMij(j,i));
          mCij(j,i) = mCij(i,j);
        }
        
 
    }
  }
}

References fAllGood, fAllValues, fAllVar, fBlock, VQwDataHandler::fDependentValues, VQwDataHandler::fDependentVar, fErrCounts_DV, fErrCounts_EF, fErrCounts_IV, fGoodEvent, fGoodEventNumber, fIndependentValues, fIndependentVar, fTotalCount, VQwDataHandler::GetEventcutErrorFlag(), mCij, mMij, mNij, and mSij.

Here is the call graph for this function:

ProcessOptions()

void GrandCorrelator::ProcessOptions

(

QwOptions &

options

)

Definition at line 69 of file GrandCorrelator.cc.

{
}

setDims()

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

inline

Definition at line 195 of file GrandCorrelator.h.

{ nP = a; nY = b;}

References nP, and nY.

Referenced by ConnectChannels().

Here is the caller graph for this function:

solve()

void GrandCorrelator::solve

(

)

Definition at line 1219 of file GrandCorrelator.cc.

{
//==========================================================
//Solve step 1
for(size_t i = 0; i < fAllVar.size(); ++i){
    for(size_t j = i; j < fAllVar.size(); ++j){
      if(mNij(i,j) >= 2){
        mVij(i,j) = mCij(i,j) / (mNij(i,j) - 1.);
        sigma_ij(i,j) = sqrt((mSij(i,j)) / (mNij(i,j) - 1.));
        sigma_ji(j,i) = sqrt((mSij(j,i)) / (mNij(j,i) - 1.));
 
        if(sigma_ij(i,j) > 0.0 && sigma_ji(j,i) > 0.0){
        mRij(i,j) = mVij(i,j) / (sigma_ij(i,j) * sigma_ji(j,i));
        } else {
            mRij(i,j) = 0.0;
        }
 
        if(j > i){
        mVij(j,i) = mVij(i,j);
        sigma_ij(j,i) = sqrt((mSij(j,i)) / (mNij(j,i) - 1.));
        sigma_ji(i,j) = sqrt((mSij(i,j)) / (mNij(i,j) - 1.));
 
        if(sigma_ij(j,i) > 0.0 && sigma_ji(i,j) > 0.0){
        mRij(j,i) = mVij(j,i) / (sigma_ij(j,i) * sigma_ji(i,j));
        } else {
            mRij(i,j) = 0.0;
        }
        }
      }
      else {
        mVij(i,j) = 0;
        sigma_ij(i,j) = 0;
        sigma_ji(j,i) = 0;
        mRij(i,j) = 0;
      }
 
      //Bottom half matrix
 
      mVij(j,i) = mVij(i,j);
 
        // Correlation
        if (sigma_ij(i,j) > 0.0 && sigma_ji(j,i) > 0.0) {
            mRij(i,j) = mVij(i,j) / (sigma_ij(i,j) * sigma_ji(j,i));
            mRij(j,i) = mRij(i,j); // mirror
        } else {
            mRij(i,j) = 0.0;
            mRij(j,i) = 0.0;
        }
    }
}
 
 
//==========================================================
//Solve step 2
  mVFULL = mCij;
  mVPY = mVFULL.GetSub(0,nP-1,nP,nP+nY-1);
  mVPP = mVFULL.GetSub(0,nP-1,0,nP-1);
  mVYY = mVFULL.GetSub(nP,nP+nY-1,nP,nP+nY-1);
 
  mRFULL = mRij;
  mRPY = mRFULL.GetSub(0,nP-1,nP,nP+nY-1);
  mRPP = mRFULL.GetSub(0,nP-1,0,nP-1);
  mRYY = mRFULL.GetSub(nP,nP+nY-1,nP,nP+nY-1);
 
  mSFULL = mVij;
  mSPY = mSFULL.GetSub(0,nP-1,nP,nP+nY-1);
  mSPP = mSFULL.GetSub(0,nP-1,0,nP-1);
  mSYY = mSFULL.GetSub(nP,nP+nY-1,nP,nP+nY-1);
 
 
  // off-diagonal raw covariance
  mVYP.Transpose(mVPY);
 
  // diagonal variances
  TMatrixD sigmaP = sigma_ij.GetSub(0,nP-1,0,nP-1);
  TMatrixD sigmaY = sigma_ij.GetSub(nP,nP+nY-1,nP,nP+nY-1);
 
  mVP = TMatrixDDiag(sigmaP);
  mVY = TMatrixDDiag(sigmaY);
 
  // "Clean" matrices
  for(size_t i = 0; i < fAllVar.size(); ++i){
    for(size_t j = i; j < fAllVar.size(); ++j){
      if(i < 5 && j < 5){
      }
        mVFULL_clean(i,j) = mRij(i,j) * sigma_ij(i,j) * sigma_ji(j,i) * (fGoodEventNumber - 1);
        mVFULL_clean(j,i) = mVFULL_clean(i,j);
        mSFULL_clean(i,j) = mRij(i,j) * sigma_ij(i,j) * sigma_ji(j,i);
        mSFULL_clean(j,i) = mSFULL_clean(i,j);
    }
}
  TMatrixD mVPY_clean = mVFULL_clean.GetSub(0,nP-1,nP,nP+nY-1);
  TMatrixD mVPP_clean = mVFULL_clean.GetSub(0,nP-1,0,nP-1);
  TMatrixD mVYY_clean = mVFULL_clean.GetSub(nP,nP+nY-1,nP,nP+nY-1);
  TMatrixD mVYP_clean(TMatrixD::kTransposed, mVPY_clean);
 
 
  TMatrixD mSPY_clean = mSFULL_clean.GetSub(0,nP-1,nP,nP+nY-1);
  TMatrixD mSPP_clean = mSFULL_clean.GetSub(0,nP-1,0,nP-1);
  TMatrixD mSYY_clean = mSFULL_clean.GetSub(nP,nP+nY-1,nP,nP+nY-1);
 
 
  TMatrixD mSYP_clean(TMatrixD::kTransposed, mSPY_clean);
 
  TVectorD mSP_clean = TMatrixDDiag(mSPP_clean);
  mSP_clean.Sqrt();
 
  TVectorD mSY_clean = TMatrixDDiag(mSYY_clean);
  mSY_clean.Sqrt();
 
  // Check for goodness and then get rid of bad columns
 
  // 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_clean.Print();
      QwMessage << "Correlation matrix: " << QwLog::endl; mRPP.Print();
    }
    QwWarning << "LRB: solving failed (this happens when only few events)."
              << QwLog::endl;
    return;
  }
 
  //==========================================================
  //Solve Step 3
  // slopes
  TMatrixD invRPP(TMatrixD::kInverted, mRPP);
  Axy = TMatrixD(invRPP, TMatrixD::kMult, mRPY);
 
  Axy.NormByColumn(mSP_clean); // divide
 
  Axy.NormByRow(mSY_clean, ""); // mult
 
  Ayx.Transpose(Axy);
 
  // new means
  mMYp = mMY - Ayx * mMP;
 
  // new raw covariance
  mVYYp = mVYY_clean + Ayx * mVPP_clean * Axy - (Ayx * mVPY_clean + mVYP_clean * Axy);
 
  // new variances
  mVYp = TMatrixDDiag(mVYYp); 
  for (int i = 0; i < mVYp.GetNrows(); i++) {
    if (mVYp(i) < 0 && fabs(mVYp(i)) < 1e-12) {
        mVYp(i) = 0;
    }
  }
 
  mVYp.Sqrt();
 
  // new normalized covariance
  mSYYp = mSYY_clean + Ayx * mSPP_clean * Axy - (Ayx * mSPY_clean + mSYP_clean * 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_clean); // divide
 
  dAxy.NormByRow(mSYp, ""); // mult
 
  dAyx.Transpose(dAxy);
 
 
  fErrorFlag = 0;
}

References Axy, Ayx, dAxy, dAyx, QwLog::endl(), fAllVar, fErrorFlag, fGoodEventNumber, mCij, mMP, mMY, mMYp, mNij, mRFULL, mRij, mRPP, mRPY, mRYY, mRYYp, mSFULL, mSFULL_clean, mSij, mSPP, mSPY, mSYp, mSYY, mSYYp, mVFULL, mVFULL_clean, mVij, mVP, mVPP, mVPY, mVY, mVYP, mVYp, mVYY, mVYYp, nP, nY, QwMessage, QwWarning, sigma_ij, and sigma_ji.

Referenced by CalcCorrelations().

Here is the call graph for this function:

Here is the caller graph for this function:

WriteAliasFile()

void GrandCorrelator::WriteAliasFile ( )

protected

Definition at line 716 of file GrandCorrelator.cc.

{
  // Ensure output file is open
  if (fAliasOutputFile.bad()) {
    QwWarning << "GrandCorrelator: 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 GrandCorrelator::WriteAlphaFile ( )

protected

Definition at line 598 of file GrandCorrelator.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+

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

friend

Definition at line 220 of file GrandCorrelator.h.

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

References GrandCorrelator().

operator<<

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

friend

Output stream operator.

Definition at line 232 of file GrandCorrelator.h.

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

References fGoodEventNumber, and GrandCorrelator().

Field Documentation

Axy

TMatrixD GrandCorrelator::Axy

private

slopes

Definition at line 177 of file GrandCorrelator.h.

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

Ayx

TMatrixD GrandCorrelator::Ayx

private

Definition at line 177 of file GrandCorrelator.h.

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

dAxy

TMatrixD GrandCorrelator::dAxy

private

Definition at line 177 of file GrandCorrelator.h.

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

dAyx

TMatrixD GrandCorrelator::dAyx

private

Definition at line 177 of file GrandCorrelator.h.

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

fAliasOutputFile

std::ofstream GrandCorrelator::fAliasOutputFile

protected

Definition at line 100 of file GrandCorrelator.h.

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

fAliasOutputFileBase

std::string GrandCorrelator::fAliasOutputFileBase

protected

Definition at line 97 of file GrandCorrelator.h.

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

fAliasOutputFileSuff

std::string GrandCorrelator::fAliasOutputFileSuff

protected

Definition at line 98 of file GrandCorrelator.h.

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

fAliasOutputPath

std::string GrandCorrelator::fAliasOutputPath

protected

Definition at line 99 of file GrandCorrelator.h.

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

fAllGood

std::vector<bool> GrandCorrelator::fAllGood

private

Definition at line 171 of file GrandCorrelator.h.

Referenced by ConnectChannels(), and ProcessData().

fAllName

std::vector<std::string> GrandCorrelator::fAllName

private

Definition at line 172 of file GrandCorrelator.h.

fAllValues

std::vector<double> GrandCorrelator::fAllValues

private

Definition at line 170 of file GrandCorrelator.h.

Referenced by ConnectChannels(), and ProcessData().

fAllVar

std::vector<const VQwHardwareChannel*> GrandCorrelator::fAllVar

private

Definition at line 169 of file GrandCorrelator.h.

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

fAlphaOutputFile

TFile* GrandCorrelator::fAlphaOutputFile

protected

Definition at line 90 of file GrandCorrelator.h.

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

fAlphaOutputFileBase

std::string GrandCorrelator::fAlphaOutputFileBase

protected

Definition at line 87 of file GrandCorrelator.h.

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

fAlphaOutputFileSuff

std::string GrandCorrelator::fAlphaOutputFileSuff

protected

Definition at line 88 of file GrandCorrelator.h.

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

fAlphaOutputPath

std::string GrandCorrelator::fAlphaOutputPath

protected

Definition at line 89 of file GrandCorrelator.h.

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

fBlock

Int_t GrandCorrelator::fBlock

protected

Definition at line 74 of file GrandCorrelator.h.

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

fCycleCounter

Int_t GrandCorrelator::fCycleCounter

private

Definition at line 127 of file GrandCorrelator.h.

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

fDisableHistos

bool GrandCorrelator::fDisableHistos

protected

Definition at line 76 of file GrandCorrelator.h.

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

fErrCounts_ALL

std::vector<int> GrandCorrelator::fErrCounts_ALL

protected

Definition at line 111 of file GrandCorrelator.h.

fErrCounts_DV

std::vector<int> GrandCorrelator::fErrCounts_DV

protected

Definition at line 110 of file GrandCorrelator.h.

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

fErrCounts_EF

int GrandCorrelator::fErrCounts_EF

protected

Definition at line 108 of file GrandCorrelator.h.

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

fErrCounts_IV

std::vector<int> GrandCorrelator::fErrCounts_IV

protected

Definition at line 109 of file GrandCorrelator.h.

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

fErrorFlag

Int_t GrandCorrelator::fErrorFlag

private

is information valid

Definition at line 134 of file GrandCorrelator.h.

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

fGoodCount

int GrandCorrelator::fGoodCount

protected

Definition at line 106 of file GrandCorrelator.h.

fGoodEvent

unsigned int GrandCorrelator::fGoodEvent

protected

Definition at line 113 of file GrandCorrelator.h.

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

fGoodEventNumber

Long64_t GrandCorrelator::fGoodEventNumber

private

accumulated so far

Definition at line 135 of file GrandCorrelator.h.

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

fH1dv

std::vector<TH1D> GrandCorrelator::fH1dv

private

Definition at line 123 of file GrandCorrelator.h.

Referenced by ConstructHistograms(), and FillHistograms().

fH1iv

std::vector<TH1D> GrandCorrelator::fH1iv

private

Definition at line 122 of file GrandCorrelator.h.

Referenced by ConstructHistograms(), and FillHistograms().

fH2dv

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

private

Definition at line 125 of file GrandCorrelator.h.

Referenced by ConstructHistograms(), and FillHistograms().

fH2iv

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

private

Definition at line 124 of file GrandCorrelator.h.

Referenced by ConstructHistograms(), and FillHistograms().

fHnames

std::vector<TH1D> GrandCorrelator::fHnames

private

Definition at line 121 of file GrandCorrelator.h.

Referenced by ConstructHistograms().

fIndependentFull

std::vector< std::string > GrandCorrelator::fIndependentFull

protected

Definition at line 78 of file GrandCorrelator.h.

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

fIndependentName

std::vector< std::string > GrandCorrelator::fIndependentName

protected

Definition at line 82 of file GrandCorrelator.h.

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

fIndependentType

std::vector< EQwHandleType > GrandCorrelator::fIndependentType

protected

Definition at line 81 of file GrandCorrelator.h.

Referenced by ConnectChannels(), and LoadChannelMap().

fIndependentValues

std::vector< Double_t > GrandCorrelator::fIndependentValues

protected

Definition at line 85 of file GrandCorrelator.h.

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

fIndependentVar

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

protected

Definition at line 84 of file GrandCorrelator.h.

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

fNameNoSpaces

TString GrandCorrelator::fNameNoSpaces

private

Definition at line 117 of file GrandCorrelator.h.

Referenced by GrandCorrelator().

fPrintCorrelations

bool GrandCorrelator::fPrintCorrelations = false

staticprivate

Definition at line 35 of file GrandCorrelator.h.

Referenced by CalcCorrelations(), and DefineOptions().

fTotalCount

int GrandCorrelator::fTotalCount

protected

Definition at line 105 of file GrandCorrelator.h.

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

fTree

TTree* GrandCorrelator::fTree

protected

Definition at line 95 of file GrandCorrelator.h.

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

mCij

TMatrixD GrandCorrelator::mCij

private

Definition at line 164 of file GrandCorrelator.h.

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

mMij

TMatrixD GrandCorrelator::mMij

private

Definition at line 163 of file GrandCorrelator.h.

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

mMP

TVectorD GrandCorrelator::mMP

private

mean values

Definition at line 159 of file GrandCorrelator.h.

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

mMY

TVectorD GrandCorrelator::mMY

private

Definition at line 159 of file GrandCorrelator.h.

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

mMYp

TVectorD GrandCorrelator::mMYp

private

Definition at line 159 of file GrandCorrelator.h.

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

mNij

TMatrixD GrandCorrelator::mNij

private

new matrices and vectors

Definition at line 163 of file GrandCorrelator.h.

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

mRFULL

TMatrixD GrandCorrelator::mRFULL

private

Definition at line 165 of file GrandCorrelator.h.

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

mRij

TMatrixD GrandCorrelator::mRij

private

Definition at line 164 of file GrandCorrelator.h.

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

mRPP

TMatrixD GrandCorrelator::mRPP

private

Definition at line 139 of file GrandCorrelator.h.

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

mRPY

TMatrixD GrandCorrelator::mRPY

private

correlations

Definition at line 138 of file GrandCorrelator.h.

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

mRYP

TMatrixD GrandCorrelator::mRYP

private

Definition at line 138 of file GrandCorrelator.h.

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

mRYY

TMatrixD GrandCorrelator::mRYY

private

Definition at line 139 of file GrandCorrelator.h.

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

mRYYp

TMatrixD GrandCorrelator::mRYYp

private

Definition at line 140 of file GrandCorrelator.h.

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

mSFULL

TMatrixD GrandCorrelator::mSFULL

private

Definition at line 165 of file GrandCorrelator.h.

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

mSFULL_clean

TMatrixD GrandCorrelator::mSFULL_clean

private

Definition at line 168 of file GrandCorrelator.h.

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

mSij

TMatrixD GrandCorrelator::mSij

private

Definition at line 163 of file GrandCorrelator.h.

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

mSP

TVectorD GrandCorrelator::mSP

private

sigmas

Definition at line 155 of file GrandCorrelator.h.

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

mSPP

TMatrixD GrandCorrelator::mSPP

private

Definition at line 152 of file GrandCorrelator.h.

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

mSPY

TMatrixD GrandCorrelator::mSPY

private

normalized covariances

Definition at line 151 of file GrandCorrelator.h.

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

mSY

TVectorD GrandCorrelator::mSY

private

Definition at line 155 of file GrandCorrelator.h.

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

mSYP

TMatrixD GrandCorrelator::mSYP

private

Definition at line 151 of file GrandCorrelator.h.

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

mSYp

TVectorD GrandCorrelator::mSYp

private

Definition at line 156 of file GrandCorrelator.h.

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

mSYY

TMatrixD GrandCorrelator::mSYY

private

Definition at line 152 of file GrandCorrelator.h.

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

mSYYp

TMatrixD GrandCorrelator::mSYYp

private

Definition at line 153 of file GrandCorrelator.h.

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

mVFULL

TMatrixD GrandCorrelator::mVFULL

private

Definition at line 165 of file GrandCorrelator.h.

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

mVFULL_clean

TMatrixD GrandCorrelator::mVFULL_clean

private

Definition at line 167 of file GrandCorrelator.h.

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

mVij

TMatrixD GrandCorrelator::mVij

private

Definition at line 164 of file GrandCorrelator.h.

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

mVP

TVectorD GrandCorrelator::mVP

private

variances

Definition at line 147 of file GrandCorrelator.h.

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

mVPP

TMatrixD GrandCorrelator::mVPP

private

Definition at line 144 of file GrandCorrelator.h.

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

mVPY

TMatrixD GrandCorrelator::mVPY

private

unnormalized covariances

Definition at line 143 of file GrandCorrelator.h.

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

mVY

TVectorD GrandCorrelator::mVY

private

Definition at line 147 of file GrandCorrelator.h.

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

mVYP

TMatrixD GrandCorrelator::mVYP

private

Definition at line 143 of file GrandCorrelator.h.

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

mVYp

TVectorD GrandCorrelator::mVYp

private

Definition at line 148 of file GrandCorrelator.h.

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

mVYY

TMatrixD GrandCorrelator::mVYY

private

Definition at line 144 of file GrandCorrelator.h.

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

mVYYp

TMatrixD GrandCorrelator::mVYYp

private

Definition at line 145 of file GrandCorrelator.h.

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

nP

int GrandCorrelator::nP

private

Definition at line 118 of file GrandCorrelator.h.

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

nY

int GrandCorrelator::nY

private

Definition at line 118 of file GrandCorrelator.h.

Referenced by CalcCorrelations(), ConnectChannels(), ConstructHistograms(), ConstructTreeBranches(), FillHistograms(), getCovariancePY(), getCovarianceY(), getMeanY(), getMeanYprime(), getSigmaY(), getSigmaYprime(), GrandCorrelator(), GrandCorrelator(), GrandCorrelator(), init(), print(), printSummaryAlphas(), printSummaryMeansWithUnc(), printSummaryMeansWithUncCorrected(), printSummaryY(), printSummaryYP(), setDims(), solve(), WriteAliasFile(), and WriteAlphaFile().

sigma_ij

TMatrixD GrandCorrelator::sigma_ij

private

Definition at line 166 of file GrandCorrelator.h.

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

sigma_ji

TMatrixD GrandCorrelator::sigma_ji

private

Definition at line 166 of file GrandCorrelator.h.

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

---

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

• GrandCorrelator.h
• GrandCorrelator.cc