AlgAdapt.h

Go to the documentation of this file.

#ifndef ANP_ALGADAPT_H
#define ANP_ALGADAPT_H

// 
// $Id: AlgAdapt.h,v 1.12 2009/10/15 18:30:13 jyuko Exp $
//
// ABC for algorithm that changes true or reco quantities in a fit

// Inheriting algorithms have access to 5 base vectors:
//    fCurErr:  current error vector<double>
//    fCurPar:  current parameter vector<double>
//    fDefErr:  default error vector<double>
//    fDefPar:  default parameter vector<double>
//    fParName: optional names of parameters
//

// C++
#include <iostream>
#include <sstream>
#include <string>
#include <vector>

// MINOS
#include "Util/UtilString.h"

// Local
#include "PhysicsNtuple/Handle.h"
#include "PhysicsNtuple/RecoNu.h"
#include "PhysicsNtuple/TrueNu.h"

class Registry;
class TDirectory;

namespace Anp
{
   class AlgAdapt: public virtual Base
   {      
   public:      
      
      AlgAdapt();
      virtual ~AlgAdapt();

      virtual bool Init() = 0;
      
      virtual long double GetPT() const;

      virtual unsigned int GetNPar() const = 0;

      virtual void InitReco(RecoNu &) const;
      virtual void InitTrue(TrueNu &) const;

      virtual double GetWeight(const TrueNu &) const;

      virtual void Vary(RecoNu &, const TrueNu &) const;

      virtual const std::vector<double>& Get(const std::string &key) const;

      virtual void Set(const std::vector<double> &dvec, const std::string &key);

      virtual const std::vector<std::string>& GetParName() const;

      virtual void Config(const Registry &) = 0;
     
      virtual void Print(std::ostream &os = std::cout) const;

      virtual void Read(const std::string &fname);
      virtual void Save(const std::string &fname);

   protected:

      void ReadVec(const std::string &prefix, const Registry &reg);

      const std::vector<double> GetVec(const std::string &key, const Registry &reg) const;

      void PrintVec(std::ostream &o, int w=7, int p=5) const;

      void PrintVec(const std::vector<double> &v, std::ostream &o, int w=7, int p=5) const;

   protected:

      //
      // !!!Derived classes must not redefine variables defined below!!!
      //
      std::vector<double> fCurErr;
      std::vector<double> fCurPar;
      std::vector<double> fDefErr;
      std::vector<double> fDefPar;
      
      std::vector<std::string> fParName;
   };

   //
   // Global helper functions
   //
   const std::string ConvertVector2String(const std::vector<double> &dvec);
   const std::vector<double> ConvertString2Vector(const std::string &dstr);

   //
   // Inlined member functions
   //
   inline const std::vector<std::string>& AlgAdapt::GetParName() const
   {
      return fParName;
   }

   inline void AlgAdapt::InitReco(RecoNu &) const {}
   inline void AlgAdapt::InitTrue(TrueNu &) const {}

   inline double AlgAdapt::GetWeight(const TrueNu &) const { return 1.0; }
   inline void AlgAdapt::Vary(RecoNu &, const TrueNu &) const {}

   //
   // Template global functions
   //

   template<class T>
      const std::string ConvertVectorToString(const std::vector<T> &dvec)
   {
      std::stringstream dstr;
      for(typename std::vector<T>::const_iterator dit = dvec.begin(); dit != dvec.end(); ++dit)
      {
         dstr << *dit;
         if(dit + 1 != dvec.end()) dstr << ", ";
      }
      
      return dstr.str();
   }

  template<class T>
     const std::vector<T> ConvertStringToVector(const std::string &dstr)
  {
     std::vector<T> dvec;

     //
     // Parse into individual strings
     //
     std::vector<std::string> svec;
     UtilString::StringTok(svec, dstr, ", ");
     
     if(svec.empty())
     {
        return dvec;
     }

     for(unsigned int i = 0; i < svec.size(); ++i)
     {
        std::stringstream temp;
        temp << svec[i];
        
        T value;
        temp >> value;
        
        if(!temp.fail())
        {
           dvec.push_back(value);
        }
        else
        {
           std::cerr << "Anp::ConvertString2Vector - stringstream >> failed for " << svec[i] << std::endl;
        }
     }
     
     return dvec;
  }
}

#endif

---