#include <QwRegression.h>

Inheritance diagram for QwRegression:

Collaboration diagram for QwRegression:

Public Types
enum	EQwRegType { kRegTypeUnknown =0, kRegTypeMps, kRegTypeAsym, kRegTypeDiff }
	Type of regression variable. More...

typedef std::vector< std::pair < VQwHardwareChannel , VQwHardwareChannel > >::iterator	PairIterator

typedef std::vector< std::pair < VQwHardwareChannel , VQwHardwareChannel > >::const_iterator	PairConstIterator

Public Member Functions
	QwRegression (QwOptions &options, QwSubsystemArrayParity &event, QwHelicityPattern &helicitypattern)
	Constructor with single event and helicity pattern. More...

	QwRegression (QwOptions &options, QwSubsystemArrayParity &event)
	Constructor with single event. More...

	QwRegression (QwOptions &options, QwHelicityPattern &helicitypattern)
	Constructor with helicity pattern. More...

	QwRegression (QwOptions &options)
	Constructor with only options. More...

	QwRegression (const QwRegression &source)
	Copy constructor. More...

virtual	~QwRegression ()
	Virtual destructor. More...

void	ProcessOptions (QwOptions &options)
	Process the configuration options. More...

Int_t	LoadChannelMap (const std::string &mapfile)
	Load the channels and sensitivities. More...

Int_t	ConnectChannels (QwSubsystemArrayParity &event, QwSubsystemArrayParity &asym, QwSubsystemArrayParity &diff)
	Connect to channels. More...

Int_t	ConnectChannels (QwSubsystemArrayParity &event)
	Connect to Channels (event only) More...

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

void	LinearRegression (EQwRegType type)
	Linear regression. More...

void	AccumulateRunningSum (QwRegression value)

void	CalculateRunningAverage ()

void	PrintValue () const

void	FillDB (QwParityDB *db, TString datatype)

Tree and vector construction and maintenance
void	ConstructBranchAndVector (TTree *tree, std::vector< Double_t > &values)
	Construct the tree and vector for this subsystem. More...

void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)
	Construct a branch and vector for this subsystem with a prefix. More...

void	FillTreeVector (std::vector< Double_t > &values) const
	Fill the vector for this subsystem. More...

Static Public Member Functions
static void	DefineOptions (QwOptions &options)
	Define the configuration options. More...

Protected Member Functions
std::pair< EQwRegType, std::string >	ParseRegressionVariable (const std::string &variable)
	Parse regression variable in type and name. More...

Protected Attributes
bool	fEnableRegression
	Flag to enable regression. More...

std::string	fRegressionMapFile
	Regression dependent and independent variable map. More...

QwSubsystemArrayParity *	fSubsystemArray
	Single event pointer. More...

QwHelicityPattern *	fHelicityPattern
	Helicity pattern pointer. More...

UInt_t	fErrorFlag

std::vector< EQwRegType >	fDependentType
	List of channels to use in the regression. More...

std::vector< std::string >	fDependentName

std::vector< std::vector < EQwRegType > >	fIndependentType

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

std::vector< std::vector < Double_t > >	fSensitivity

std::vector< EQwRegType >	fDependentVarType

std::vector< std::pair < VQwHardwareChannel , VQwHardwareChannel > >	fDependentVar

std::vector< std::vector < std::pair< Double_t, const VQwHardwareChannel * > > >	fIndependentVar

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

Detailed Description

... info to follow ... really, I promise ...

Definition at line 19 of file QwRegression.h.

Member Typedef Documentation

typedef std::vector< std::pair< VQwHardwareChannel*,VQwHardwareChannel*> >::const_iterator QwRegression::PairConstIterator

Definition at line 29 of file QwRegression.h.

typedef std::vector< std::pair< VQwHardwareChannel*,VQwHardwareChannel*> >::iterator QwRegression::PairIterator

Definition at line 28 of file QwRegression.h.

Member Enumeration Documentation

enum QwRegression::EQwRegType

Type of regression variable.

Enumerator
kRegTypeUnknown
kRegTypeMps
kRegTypeAsym
kRegTypeDiff

Definition at line 24 of file QwRegression.h.

                     {
       kRegTypeUnknown=0, kRegTypeMps, kRegTypeAsym, kRegTypeDiff
     };

Constructor & Destructor Documentation

QwRegression::QwRegression	(	QwOptions &	options,
		QwSubsystemArrayParity &	event,
		QwHelicityPattern &	helicitypattern
	)

Constructor with single event and helicity pattern.

Constructor with single event and helicity pattern

Parameters

options	Options object
event	Single event
helicitypattern	Helicity pattern

Definition at line 25 of file QwRegression.cc.

References ConnectChannels(), QwHelicityPattern::fAsymmetry, QwHelicityPattern::fDifference, fEnableRegression, fHelicityPattern, fRegressionMapFile, fSubsystemArray, LoadChannelMap(), and ProcessOptions().

 {
   fEnableRegression = false;
   ProcessOptions(options);
   fSubsystemArray = &event;
   fHelicityPattern = &helicitypattern;
   LoadChannelMap(fRegressionMapFile);
   QwSubsystemArrayParity& asym = helicitypattern.fAsymmetry;
   QwSubsystemArrayParity& diff = helicitypattern.fDifference;
   ConnectChannels(event,asym,diff);
 }

Here is the call graph for this function:

QwRegression::QwRegression	(	QwOptions &	options,
		QwSubsystemArrayParity &	event
	)

Constructor with single event.

Constructor with single event

Parameters

options	Options object
event	Single event

Definition at line 46 of file QwRegression.cc.

References ConnectChannels(), fEnableRegression, fRegressionMapFile, fSubsystemArray, LoadChannelMap(), and ProcessOptions().

                                                                            :fHelicityPattern(NULL)
 {
     fEnableRegression = false;
     ProcessOptions(options);
     fSubsystemArray = &event;
     LoadChannelMap(fRegressionMapFile);
     ConnectChannels(event);
       
 }

Here is the call graph for this function:

QwRegression::QwRegression	(	QwOptions &	options,
		QwHelicityPattern &	helicitypattern
	)

Constructor with helicity pattern.

Constructor with only helicity pattern

Parameters

options	Obtions object
helicitypattern	Helicity Pattern

Definition at line 62 of file QwRegression.cc.

References ConnectChannels(), QwHelicityPattern::fAsymmetry, QwHelicityPattern::fDifference, fEnableRegression, fHelicityPattern, fRegressionMapFile, LoadChannelMap(), and ProcessOptions().

                                                                                 :fSubsystemArray(NULL)
 {
   fEnableRegression = false;
   ProcessOptions(options);
   LoadChannelMap(fRegressionMapFile);
   fHelicityPattern = &helicitypattern;
   QwSubsystemArrayParity& asym = helicitypattern.fAsymmetry;
   QwSubsystemArrayParity& diff = helicitypattern.fDifference;
   ConnectChannels(asym,diff);
 }

Here is the call graph for this function:

QwRegression::QwRegression ( QwOptions & options )

Constructor with only options.

Constructor with options

Parameters

options Obtions object

Definition at line 78 of file QwRegression.cc.

References fEnableRegression, fRegressionMapFile, LoadChannelMap(), and ProcessOptions().

 {
   fEnableRegression = false;
   ProcessOptions(options);
   LoadChannelMap(fRegressionMapFile);
 } 

Here is the call graph for this function:

QwRegression::QwRegression ( const QwRegression & source )

Copy constructor.

Definition at line 86 of file QwRegression.cc.

References fDependentVar, fDependentVarType, and VQwDataElement::kDerived.

                                                     :
   fEnableRegression(source.fEnableRegression),
   fRegressionMapFile(source.fRegressionMapFile),
   fErrorFlag(source.fErrorFlag)
 {
     this->fDependentVar.resize(source.fDependentVar.size());
     fDependentVarType.resize(source.fDependentVar.size());
     for (size_t i = 0; i < this->fDependentVar.size(); i++) {
       const QwVQWK_Channel* vqwk = dynamic_cast<const QwVQWK_Channel*>(source.fDependentVar[i].second);
       this->fDependentVar[i].first = NULL;
       this->fDependentVar[i].second = new QwVQWK_Channel(*vqwk, VQwDataElement::kDerived);
       fDependentVarType[i] = source.fDependentVarType[i];
     }
 }

QwRegression::~QwRegression ( )

virtual

Virtual destructor.

Destructor.

Definition at line 102 of file QwRegression.cc.

References fDependentVar.

 {
   PairIterator element;
   for (element = fDependentVar.begin();
       element != fDependentVar.end(); element++) {
       if (element->second != NULL){
          delete element->second;
       }
   }
   fDependentVar.clear();
 }

QwRegression::QwRegression ( )

inlineprivate

Private default constructor.

Definition at line 94 of file QwRegression.h.

94 { };

Member Function Documentation

void QwRegression::AccumulateRunningSum ( QwRegression value )

Definition at line 612 of file QwRegression.cc.

References fDependentVar, and fErrorFlag.

Referenced by QwRegressionSubsystem::AccumulateRunningSum(), and main().

 {
   if (value.fErrorFlag==0){
     for (size_t i = 0; i < value.fDependentVar.size(); i++){
       this->fDependentVar[i].second->AccumulateRunningSum(value.fDependentVar[i].second);
     }
   }
 }

Here is the caller graph for this function:

void QwRegression::CalculateRunningAverage ( )

Definition at line 622 of file QwRegression.cc.

References fDependentVar.

Referenced by QwRegressionSubsystem::CalculateRunningAverage(), and main().

 {
   for(size_t i = 0; i < fDependentVar.size(); i++) {
     fDependentVar[i].second->CalculateRunningAverage();
   }
   
   return;
 }

Here is the caller graph for this function:

Int_t QwRegression::ConnectChannels	(	QwSubsystemArrayParity &	event,
		QwSubsystemArrayParity &	asym,
		QwSubsystemArrayParity &	diff
	)

Connect to channels.

Connect to the dependent and independent channels

Parameters

event	Helicity event structure
asym	Asymmetry event structure
diff	Difference event structure

Returns: Zero on success

Fill vector of pointers to the relevant data elements

Definition at line 222 of file QwRegression.cc.

References QwLog::endl(), fDependentName, fDependentType, fDependentVar, fDependentVarType, fEnableRegression, fIndependentName, fIndependentType, fIndependentVar, fSensitivity, VQwDataElement::GetElementName(), VQwDataElement::kDerived, kRegTypeAsym, kRegTypeDiff, kRegTypeMps, QwMessage, QwWarning, QwSubsystemArray::ReturnInternalValue(), and VQwDataElement::SetElementName().

Referenced by QwRegression().

 {
   // Return if regression is not enabled
   if (! fEnableRegression) return 0;
 
   /// Fill vector of pointers to the relevant data elements
   for (size_t dv = 0; dv < fDependentName.size(); dv++) {
     // Get the dependent variables
 
     VQwHardwareChannel* dv_ptr = 0;
     QwVQWK_Channel* new_vqwk = NULL;
     QwVQWK_Channel* vqwk = NULL;
     string name = "";
     string reg = "reg_";
     
     if(fDependentName.at(dv).at(0) == '@' ){
         name = fDependentName.at(dv).substr(1,fDependentName.at(dv).length());
     }else{
         switch (fDependentType.at(dv)) {
           case kRegTypeMps:
             dv_ptr = event.ReturnInternalValueForFriends(fDependentName.at(dv));
             break;
           case kRegTypeAsym:
             dv_ptr = asym.ReturnInternalValueForFriends(fDependentName.at(dv));
             break;
           case kRegTypeDiff:
             dv_ptr = diff.ReturnInternalValueForFriends(fDependentName.at(dv));
             break;
           default:
             QwWarning << "Dependent variable for regression has unknown type."
                       << QwLog::endl;
             break;
         }
 
         vqwk = dynamic_cast<QwVQWK_Channel*>(dv_ptr);
         name = vqwk->GetElementName().Data();
         name.insert(0, reg);
         new_vqwk = new QwVQWK_Channel(*vqwk, VQwDataElement::kDerived);
         new_vqwk->SetElementName(name);
     }
 
     // alias
     if(fDependentName.at(dv).at(0) == '@'){
         QwMessage << "dv: " << name << QwLog::endl;
         new_vqwk = new QwVQWK_Channel(name, VQwDataElement::kDerived);
     }
     // defined type
     else if(dv_ptr!=NULL){
         QwMessage << "dv: " << fDependentName.at(dv) << QwLog::endl;
     }else {
         QwWarning << "Dependent variable " << fDependentName.at(dv) << " could not be found, "
                 << "or is not a VQWK channel." << QwLog::endl;
                 continue; 
     }
 
     // pair creation
     if(new_vqwk != NULL){
       fDependentVarType.push_back(fDependentType.at(dv));
       fDependentVar.push_back(std::make_pair(vqwk, new_vqwk));
     }
 
       // Add independent variables
       fIndependentVar.resize(fDependentVar.size());
       for (size_t iv = 0; iv < fIndependentName.at(dv).size(); iv++) {
         // Get the independent variables
         const VQwHardwareChannel* iv_ptr = 0;
         switch (fIndependentType.at(dv).at(iv)) {
           case kRegTypeMps:
             iv_ptr = event.ReturnInternalValue(fIndependentName.at(dv).at(iv));
             break;
           case kRegTypeAsym:
             iv_ptr = asym.ReturnInternalValue(fIndependentName.at(dv).at(iv));
             break;
           case kRegTypeDiff:
             iv_ptr = diff.ReturnInternalValue(fIndependentName.at(dv).at(iv));
             break;
           default:
             QwWarning << "Independent variable for regression has unknown type."
                       << QwLog::endl;
             break;
         }
         if (iv_ptr) {
           QwMessage << " iv: " << fIndependentName.at(dv).at(iv) << " (sens = "
                     << fSensitivity.at(dv).at(iv) << ")" << QwLog::endl;
           fIndependentVar.back().push_back(std::make_pair(fSensitivity.at(dv).at(iv),
                 iv_ptr));
         } else {
           QwWarning << "Independent variable " << fIndependentName.at(dv).at(iv) << " for regression of "
                     << "dependent variable " << fDependentName.at(dv) << " could not be found."
                     << QwLog::endl;
         }
       }
   }
   
   return 0;
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwRegression::ConnectChannels ( QwSubsystemArrayParity & event )

Connect to Channels (event only)

Connect to the dependent and independent channels

Parameters

event Helicity event structure

Returns: Zero on success

Fill vector of pointers to the relevant data elements

Definition at line 432 of file QwRegression.cc.

References QwLog::endl(), fDependentName, fDependentType, fDependentVar, fDependentVarType, fEnableRegression, fIndependentName, fIndependentType, fIndependentVar, fSensitivity, VQwDataElement::GetElementName(), VQwDataElement::kDerived, kRegTypeAsym, kRegTypeDiff, kRegTypeMps, QwMessage, QwWarning, and VQwDataElement::SetElementName().

 {
   // Return if regression is not enabled
   if (! fEnableRegression) return 0;
 
   /// Fill vector of pointers to the relevant data elements
   for (size_t dv = 0; dv < fDependentName.size(); dv++) {
     // Get the dependent variables
 
     VQwHardwareChannel* dv_ptr = 0;
     QwVQWK_Channel* new_vqwk = NULL;
     QwVQWK_Channel* vqwk = NULL;
     string name = " s";
     string reg = "reg_";
 
     if (fDependentType.at(dv)==kRegTypeAsym || fDependentType.at(dv)==kRegTypeDiff){
       //  Quietly skip the asymmetry or difference types.
       continue;
     } else if(fDependentType.at(dv) != kRegTypeMps){
       QwWarning << "QwRegression::ConnectChannels(QwSubsystemArrayParity& event):  Dependent variable, "
         << fDependentName.at(dv)
         << ", for MPS regression does not have MPS type, type=="
         << fDependentType.at(dv) << "."<< QwLog::endl;
       continue;
     } else {
       if(fDependentName.at(dv).at(0) == '@' ){
         name = fDependentName.at(dv).substr(1,fDependentName.at(dv).length());
         new_vqwk = new QwVQWK_Channel(name, VQwDataElement::kDerived);
       } else {
         dv_ptr = event.ReturnInternalValueForFriends(fDependentName.at(dv));
 
   vqwk = dynamic_cast<QwVQWK_Channel*>(dv_ptr);
         name = vqwk->GetElementName().Data();
         name.insert(0,reg);
         new_vqwk = new QwVQWK_Channel(*vqwk, VQwDataElement::kDerived);
         new_vqwk->SetElementName(name);
       }
     }
 
     // alias
     if(new_vqwk==NULL){
       QwWarning << "Dependent variable " << fDependentName.at(dv) << " could not be found, "
                 << "or is not a VQWK channel." << QwLog::endl;
       continue; 
     } else {
       QwMessage << "dv: " << new_vqwk->GetElementName() << QwLog::endl;
       // pair creation
       fDependentVarType.push_back(fDependentType.at(dv));
       fDependentVar.push_back(std::make_pair(vqwk, new_vqwk));
     }
 
       // Add independent variables
       fIndependentVar.resize(fDependentVar.size());
       for (size_t iv = 0; iv < fIndependentName.at(dv).size(); iv++) {
         // Get the independent variables
         const VQwHardwareChannel* iv_ptr = 0;
         if(fIndependentType.at(dv).at(iv) == kRegTypeMps){
             iv_ptr = event.ReturnInternalValue(fIndependentName.at(dv).at(iv));
   }else{
             QwWarning << "Independent variable for MPS regression has unknown type."
                       << QwLog::endl;
         }
         if (iv_ptr) {
           QwMessage << " iv: " << fIndependentName.at(dv).at(iv) << " (sens = "
                     << fSensitivity.at(dv).at(iv) << ")" << QwLog::endl;
           fIndependentVar.back().push_back(std::make_pair(fSensitivity.at(dv).at(iv),
                 iv_ptr));
         } else {
           QwWarning << "Independent variable " << fIndependentName.at(dv).at(iv) << " for regression of "
                     << "dependent variable " << fDependentName.at(dv) << " could not be found."
                     << QwLog::endl;
         }
       }
   }
   
   return 0;
 }

Here is the call graph for this function:

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

Connect to Channels (asymmetry/difference only)

Connect to the dependent and independent channels

Parameters

asym	Asymmetry event structure
diff	Difference event structure

Returns: Zero on success

Fill vector of pointers to the relevant data elements

Definition at line 328 of file QwRegression.cc.

References QwLog::endl(), fDependentName, fDependentType, fDependentVar, fDependentVarType, fEnableRegression, fIndependentName, fIndependentType, fIndependentVar, fSensitivity, VQwDataElement::GetElementName(), VQwDataElement::kDerived, kRegTypeAsym, kRegTypeDiff, kRegTypeMps, QwMessage, QwWarning, QwSubsystemArray::ReturnInternalValue(), and VQwDataElement::SetElementName().

 {
   // Return if regression is not enabled
   if (! fEnableRegression) return 0;
 
   /// Fill vector of pointers to the relevant data elements
   for (size_t dv = 0; dv < fDependentName.size(); dv++) {
     // Get the dependent variables
 
     VQwHardwareChannel* dv_ptr = 0;
     QwVQWK_Channel* new_vqwk = NULL;
     QwVQWK_Channel* vqwk = NULL;
     string name = "";
     string reg = "reg_";
     
     if (fDependentType.at(dv)==kRegTypeMps){
       //  Quietly ignore the MPS type when we're connecting the asym & diff
       continue;
     } else if(fDependentName.at(dv).at(0) == '@' ){
         name = fDependentName.at(dv).substr(1,fDependentName.at(dv).length());
     }else{
         switch (fDependentType.at(dv)) {
           case kRegTypeAsym:
             dv_ptr = asym.ReturnInternalValueForFriends(fDependentName.at(dv));
             break;
           case kRegTypeDiff:
             dv_ptr = diff.ReturnInternalValueForFriends(fDependentName.at(dv));
             break;
           default:
             QwWarning << "QwRegression::ConnectChannels(QwSubsystemArrayParity& asym, QwSubsystemArrayParity& diff):  Dependent variable, "
           << fDependentName.at(dv)
           << ", for asym/diff regression does not have proper type, type=="
           << fDependentType.at(dv) << "."<< QwLog::endl;
             break;
         }
 
         vqwk = dynamic_cast<QwVQWK_Channel*>(dv_ptr);
         name = vqwk->GetElementName().Data();
         name.insert(0, reg);
         new_vqwk = new QwVQWK_Channel(*vqwk, VQwDataElement::kDerived);
         new_vqwk->SetElementName(name);
     }
 
     // alias
     if(fDependentName.at(dv).at(0) == '@'){
         QwMessage << "dv: " << name << QwLog::endl;
         new_vqwk = new QwVQWK_Channel(name, VQwDataElement::kDerived);
     }
     // defined type
     else if(dv_ptr!=NULL){
         QwMessage << "dv: " << fDependentName.at(dv) << QwLog::endl;
     }else {
         QwWarning << "Dependent variable " << fDependentName.at(dv) << " could not be found, "
                 << "or is not a VQWK channel." << QwLog::endl;
                 continue; 
     }
 
     // pair creation
     if(new_vqwk != NULL){
       fDependentVarType.push_back(fDependentType.at(dv));
       fDependentVar.push_back(std::make_pair(vqwk, new_vqwk));
     }
 
       // Add independent variables
       fIndependentVar.resize(fDependentVar.size());
       for (size_t iv = 0; iv < fIndependentName.at(dv).size(); iv++) {
         // Get the independent variables
         const VQwHardwareChannel* iv_ptr = 0;
         switch (fIndependentType.at(dv).at(iv)) {
           case kRegTypeAsym:
             iv_ptr = asym.ReturnInternalValue(fIndependentName.at(dv).at(iv));
             break;
           case kRegTypeDiff:
             iv_ptr = diff.ReturnInternalValue(fIndependentName.at(dv).at(iv));
             break;
           default:
             QwWarning << "Independent variable for regression has unknown type."
                       << QwLog::endl;
             break;
         }
         if (iv_ptr) {
           QwMessage << " iv: " << fIndependentName.at(dv).at(iv) << " (sens = "
                     << fSensitivity.at(dv).at(iv) << ")" << QwLog::endl;
           fIndependentVar.back().push_back(std::make_pair(fSensitivity.at(dv).at(iv),
                 iv_ptr));
         } else {
           QwWarning << "Independent variable " << fIndependentName.at(dv).at(iv) << " for regression of "
                     << "dependent variable " << fDependentName.at(dv) << " could not be found."
                     << QwLog::endl;
         }
       }
   }
   
   return 0;
 
 }

Here is the call graph for this function:

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

inline

Construct the tree and vector for this subsystem.

Definition at line 74 of file QwRegression.h.

Referenced by QwRegressionSubsystem::ConstructBranchAndVector().

                                                                              {
       TString tmpstr("");
       ConstructBranchAndVector(tree,tmpstr,values);
     };

Here is the caller graph for this function:

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

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

Construct the branch and tree vector

Parameters

tree	Tree
prefix	Prefix
values	Vector of values

Definition at line 584 of file QwRegression.cc.

References fDependentVar.

 {
   PairConstIterator element;
   for (element = fDependentVar.begin();
        element != fDependentVar.end(); ++element) {
     element->second->ConstructBranchAndVector(tree, prefix, values);
   }
 };

void QwRegression::DefineOptions ( QwOptions & options )

static

Define the configuration options.

Defines configuration options using QwOptions functionality.

Parameters

options Options object

Definition at line 514 of file QwRegression.cc.

References QwOptions::AddOptions().

Referenced by DefineOptionsParity().

 {
   options.AddOptions("Linear regression")
     ("enable-regression", po::value<bool>()->zero_tokens()->default_value(false),
      "enable linear regression");
   options.AddOptions("Linear regression")
     ("regression-map", po::value<std::string>()->default_value("regression.map"),
      "variables and sensitivities for regression");
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwRegression::FillDB	(	QwParityDB *	db,
		TString	datatype
	)

Definition at line 640 of file QwRegression.cc.

References QwDatabase::Connect(), QwDBInterface::DetermineMeasurementTypeID(), QwDatabase::Disconnect(), QwLog::endl(), fDependentVar, fEnableRegression, QwParityDB::GetAnalysisID(), QwDBInterface::kQwDBI_BeamTable, QwDBInterface::kQwDBI_LumiTable, QwDBInterface::kQwDBI_MDTable, QwDBInterface::kQwDBI_OtherTable, QwDatabase::Query(), QwDebug, QwError, and QwMessage.

Referenced by main().

 {
   if (! fEnableRegression){
     QwDebug << " QwRegression::FillDB: Regression is not enabled!" 
       << QwLog::endl;
     return;
   }
   Bool_t local_print_flag = kTRUE;
 
   UInt_t analysis_id = db->GetAnalysisID();
 
   TString measurement_type;
   measurement_type = QwDBInterface::DetermineMeasurementTypeID(datatype);
 
   std::vector<QwDBInterface> interface;
 
   std::vector<QwParitySSQLS::beam>      beamlist;
   std::vector<QwParitySSQLS::md_data>   mdlist;
   std::vector<QwParitySSQLS::lumi_data> lumilist;
 
   QwDBInterface::EQwDBIDataTableType tabletype;
 
 
   for(size_t i = 0; i < fDependentVar.size(); i++) {
     interface.clear();
     fDependentVar[i].second->AddEntriesToList(interface);
     for(size_t j=0; j<interface.size(); j++) {
       interface.at(j).SetAnalysisID( analysis_id ) ;
       interface.at(j).SetMeasurementTypeID( measurement_type );
       tabletype = interface.at(j).SetDetectorID( db );
       if (tabletype==QwDBInterface::kQwDBI_OtherTable){
   TString tmp_name = interface.at(j).GetDeviceName();
   tmp_name.Remove(0,5);
   interface.at(j).SetDetectorName(tmp_name);
   tabletype = interface.at(j).SetDetectorID( db );
       }
       if (tabletype==QwDBInterface::kQwDBI_BeamTable){
   interface.at(j).AddThisEntryToList( beamlist );
       } else if (tabletype==QwDBInterface::kQwDBI_MDTable){
   interface.at(j).AddThisEntryToList( mdlist );
       } else if (tabletype==QwDBInterface::kQwDBI_LumiTable){
   interface.at(j).AddThisEntryToList( lumilist );
       } else {
   QwError << "QwRegression::FillDB:  Unrecognized detector name:  " 
     << interface.at(j).GetDeviceName() << QwLog::endl;
       }
       interface.at(j).PrintStatus( local_print_flag);
     }
   }
 
   db->Connect();
   // Check the entrylist size, if it isn't zero, start to query..
   if( beamlist.size() ) {
     mysqlpp::Query query= db->Query();
     query.insert(beamlist.begin(), beamlist.end());
     query.execute();
   } else {
     QwMessage << "QwRegression::FillDB :: This is the case when the beamlist contains nothing for type="<< measurement_type.Data() 
         << QwLog::endl;
   }
   if( mdlist.size() ) {
     mysqlpp::Query query= db->Query();
     query.insert(mdlist.begin(), mdlist.end());
     query.execute();
   } else {
     QwMessage << "QwRegression::FillDB :: This is the case when the mdlist contains nothing for type="<< measurement_type.Data() 
         << QwLog::endl;
   }
   if( lumilist.size() ) {
     mysqlpp::Query query= db->Query();
     query.insert(lumilist.begin(), lumilist.end());
     query.execute();
   } else {
     QwMessage << "QwRegression::FillDB :: This is the case when the lumilist contains nothing for type="<< measurement_type.Data() 
         << QwLog::endl;
   }
   db->Disconnect();
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwRegression::FillTreeVector ( std::vector< Double_t > & values ) const

Fill the vector for this subsystem.

Fill the tree vector

Parameters

values Vector of values

Definition at line 600 of file QwRegression.cc.

References fDependentVar.

Referenced by QwRegressionSubsystem::FillTreeVector().

 {
   // Fill the data element
   PairConstIterator element;
   for (element = fDependentVar.begin();
        element != fDependentVar.end(); ++element) {
     if (element->second == NULL) continue;
     element->second->FillTreeVector(values);
   }
 };

Here is the caller graph for this function:

void QwRegression::LinearRegression ( EQwRegType type )

Linear regression.

Do the linear regression.

Definition at line 536 of file QwRegression.cc.

References QwLog::endl(), fDependentVar, fDependentVarType, fEnableRegression, fErrorFlag, fHelicityPattern, fIndependentVar, fSubsystemArray, QwSubsystemArrayParity::GetEventcutErrorFlag(), QwHelicityPattern::GetEventcutErrorFlag(), QwDebug, and QwError.

Referenced by QwRegressionSubsystem::LinearRegression(), and main().

 {
   // Return if regression is not enabled
   if (! fEnableRegression){
     QwDebug << "Regression is not enabled!" << QwLog::endl;
     return;
   }
   // Get error flag from QwHelicityPattern
   if (fHelicityPattern != NULL){
     fErrorFlag = fHelicityPattern->GetEventcutErrorFlag();
   } else if (fSubsystemArray != NULL){
     fErrorFlag = fSubsystemArray->GetEventcutErrorFlag();
   } else {
     QwError << "QwRegression::LinearRegression: Can't set fErrorFlag" << QwLog::endl;
     fErrorFlag = 0;
   }
   // Linear regression for each dependent variable
   for (size_t dv = 0; dv < fDependentVar.size(); dv++) {
     // if second is NULL, can't do regression
     if (fDependentVar.at(dv).second == NULL){
         QwError<<"Second is value is NULL, unable to calculate regression."<<QwLog::endl;
         continue;
     }
     // For correct type (asym, diff, mps)
     if (fDependentVarType.at(dv) != type) continue;
 
     // Clear data in second, if first is NULL
     if (fDependentVar.at(dv).first == NULL){
         fDependentVar.at(dv).second->ClearEventData();
     }else{
         // Update second value
         fDependentVar.at(dv).second->AssignValueFrom(fDependentVar.at(dv).first);
     }
     // Add corrections
     for (size_t iv = 0; iv < fIndependentVar.at(dv).size(); iv++) {
       fDependentVar.at(dv).second->ScaledAdd(fIndependentVar.at(dv).at(iv).first,
                fIndependentVar.at(dv).at(iv).second);
     }
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Load the channels and sensitivities.

Load the channel map

Parameters

mapfile Filename of map file

Returns: Zero when success

Definition at line 152 of file QwRegression.cc.

References QwLog::endl(), fDependentName, fDependentType, fEnableRegression, fIndependentName, fIndependentType, fSensitivity, QwParameterFile::GetNextToken(), QwParameterFile::LineHasSectionHeader(), QwParameterFile::LineIsEmpty(), ParseRegressionVariable(), QwError, QwParameterFile::ReadNextLine(), QwParameterFile::ReadNextSection(), QwParameterFile::TrimComment(), QwParameterFile::TrimSectionHeader(), and QwParameterFile::TrimWhitespace().

Referenced by QwRegression().

 {
   // Return if regression is not enabled
   if (! fEnableRegression) return 0;
 
   // Open the file
   QwParameterFile map(mapfile);
 
   // Read the sections of dependent variables
   bool keep_header = true;
   std::string section_name;
   QwParameterFile* section = 0;
   std::pair<EQwRegType,std::string> type_name;
   while ((section = map.ReadNextSection(section_name,keep_header))) {
 
     // Store index to the current position in the dv vector
     size_t current_dv_start = fDependentName.size();
 
     // Add dependent variables from the section header
     section->ReadNextLine();
     if (section->LineHasSectionHeader()) {
       section->TrimSectionHeader();
       section->TrimWhitespace();
       // Parse section header into tokens separated by a comma
       std::string current_token;
       std::string previous_token;
       do {
         previous_token = current_token;
         current_token = section->GetNextToken(",");
         if (current_token.size() == 0) continue;
         // Parse current token into dependent variable type and name
         type_name = ParseRegressionVariable(current_token);
         fDependentType.push_back(type_name.first);
         fDependentName.push_back(type_name.second);
         // Resize the vectors of sensitivities and independent variables
         fSensitivity.resize(fDependentName.size());
         fIndependentType.resize(fDependentName.size());
         fIndependentName.resize(fDependentName.size());
       } while (current_token.size() != 0);
     } else QwError << "Section does not start with header." << QwLog::endl;
 
     // Add independent variables and sensitivities
     while (section->ReadNextLine()) {
       // Throw away comments, whitespace, empty lines
       section->TrimComment();
       section->TrimWhitespace();
       if (section->LineIsEmpty()) continue;
       // Get first token: independent variable
       std::string current_token = section->GetNextToken(",");
       // Parse current token into independent variable type and name
       type_name = ParseRegressionVariable(current_token);
       // Loop over dependent variables to set sensitivities
       for (size_t dv = current_dv_start; dv < fDependentName.size(); dv++) {
         Double_t sensitivity = atof(section->GetNextToken(",").c_str());
         fSensitivity.at(dv).push_back(sensitivity);
         fIndependentType.at(dv).push_back(type_name.first);
         fIndependentName.at(dv).push_back(type_name.second);
       }
     }
   }
   return 0;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

std::pair< QwRegression::EQwRegType, std::string > QwRegression::ParseRegressionVariable ( const std::string & variable )

protected

Parse regression variable in type and name.

Parse the variable type and name from a section header

Parameters

variable String with variable type and name

Returns: Pair of variable type and name

Definition at line 120 of file QwRegression.cc.

References kRegTypeAsym, kRegTypeDiff, kRegTypeMps, and kRegTypeUnknown.

Referenced by LoadChannelMap().

 {
   std::pair<EQwRegType,std::string> type_name;
   size_t len = variable.length();
   size_t pos1 = variable.find_first_of(":");
   size_t pos2 = variable.find_first_not_of(":",pos1);
   if (pos1 == std::string::npos) {
     type_name.first  = kRegTypeUnknown;
     type_name.second = variable;
   } else {
     std::string type = variable.substr(0,pos1);
     std::string name = variable.substr(pos2,len-pos2);
     if (type == "asym")
       type_name.first = kRegTypeAsym;
     else if (type == "diff")
       type_name.first = kRegTypeDiff;
     else if (type == "mps")
       type_name.first = kRegTypeMps;
     else
       type_name.first = kRegTypeUnknown;
     type_name.second = name;
   }
   return type_name;
 }

Here is the caller graph for this function:

void QwRegression::PrintValue ( ) const

Definition at line 632 of file QwRegression.cc.

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

Referenced by main(), and QwRegressionSubsystem::PrintValue().

 {
   QwMessage<<"=== QwRegression ==="<<QwLog::endl<<QwLog::endl;
   for(size_t i = 0; i < fDependentVar.size(); i++) {
     fDependentVar[i].second->PrintValue();
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwRegression::ProcessOptions ( QwOptions & options )

Process the configuration options.

Process configuration options using QwOptions functionality.

Parameters

options Options object

Definition at line 528 of file QwRegression.cc.

References fEnableRegression, fRegressionMapFile, and QwOptions::GetValue().

Referenced by QwRegressionSubsystem::ProcessOptions(), and QwRegression().

 {
   fEnableRegression = options.GetValue<bool>("enable-regression");
   fRegressionMapFile = options.GetValue<std::string>("regression-map");
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Field Documentation

std::vector< std::string > QwRegression::fDependentName

protected

Definition at line 115 of file QwRegression.h.

Referenced by ConnectChannels(), and LoadChannelMap().

std::vector< EQwRegType > QwRegression::fDependentType

protected

List of channels to use in the regression.

Definition at line 114 of file QwRegression.h.

Referenced by ConnectChannels(), and LoadChannelMap().

std::vector< std::pair< VQwHardwareChannel*, VQwHardwareChannel*> > QwRegression::fDependentVar

protected

std::vector< EQwRegType > QwRegression::fDependentVarType

protected

Definition at line 120 of file QwRegression.h.

Referenced by ConnectChannels(), LinearRegression(), and QwRegression().

bool QwRegression::fEnableRegression

protected

Flag to enable regression.

Definition at line 94 of file QwRegression.h.

Referenced by ConnectChannels(), FillDB(), LinearRegression(), LoadChannelMap(), ProcessOptions(), and QwRegression().

UInt_t QwRegression::fErrorFlag

protected

Definition at line 108 of file QwRegression.h.

Referenced by AccumulateRunningSum(), and LinearRegression().

QwHelicityPattern* QwRegression::fHelicityPattern

protected

Helicity pattern pointer.

Definition at line 106 of file QwRegression.h.

Referenced by LinearRegression(), and QwRegression().

std::vector< std::vector< std::string > > QwRegression::fIndependentName

protected

Definition at line 117 of file QwRegression.h.

Referenced by ConnectChannels(), and LoadChannelMap().

std::vector< std::vector< EQwRegType > > QwRegression::fIndependentType

protected

Definition at line 116 of file QwRegression.h.

Referenced by ConnectChannels(), and LoadChannelMap().

std::vector< std::vector< std::pair< Double_t, const VQwHardwareChannel* > > > QwRegression::fIndependentVar

protected

Definition at line 122 of file QwRegression.h.

Referenced by ConnectChannels(), and LinearRegression().

std::string QwRegression::fRegressionMapFile

protected

Regression dependent and independent variable map.

Definition at line 101 of file QwRegression.h.

Referenced by ProcessOptions(), and QwRegression().

std::vector< std::vector< Double_t> > QwRegression::fSensitivity

protected

Definition at line 118 of file QwRegression.h.

Referenced by ConnectChannels(), and LoadChannelMap().

QwSubsystemArrayParity* QwRegression::fSubsystemArray

protected

Single event pointer.

Definition at line 104 of file QwRegression.h.

Referenced by LinearRegression(), and QwRegression().

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

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Field Documentation