RooFactoryWSTool.cxx

Go to the documentation of this file.

/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitCore                                                       *
 * @(#)root/roofitcore:$Id: RooFactoryWSTool.cxx 37172 2010-12-02 08:44:43Z wouter $
 * Authors:                                                                  *
 *   WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu       *
 *   DK, David Kirkby,    UC Irvine,         dkirkby@uci.edu                 *
 *                                                                           *
 * Copyright (c) 2000-2005, Regents of the University of California          *
 *                          and Stanford University. All rights reserved.    *
 *                                                                           *
 * Redistribution and use in source and binary forms,                        *
 * with or without modification, are permitted according to the terms        *
 * listed in LICENSE (http://roofit.sourceforge.net/license.txt)             *
 *****************************************************************************/

//////////////////////////////////////////////////////////////////////////////
// 
// BEGIN_HTML
// RooFactoryWSTool is a clase like TTree::MakeClass() that generates
// skeleton code for RooAbsPdf and RooAbsReal functions given
// a list of input parameter names. The factory can also compile
// the generated code on the fly, and on request also immediate
// instantiate objects.
// END_HTML
//

#include "RooFit.h"

#include "RooFactoryWSTool.h"
#include "RooAbsReal.h"
#include "RooAbsCategory.h"
#include "RooArgList.h"
#include "RooRealVar.h"
#include "RooCategory.h"
#include "RooMsgService.h"
#include "RooWorkspace.h"
#include "TInterpreter.h"
#include "TClass.h"
#include "TClassTable.h"
#include "RooAbsPdf.h"
#include "RooGaussian.h"
#include <fstream>
#include <vector>
#include <string>
#include "RooGlobalFunc.h"
#include "RooDataSet.h"
#include "RooDataHist.h"
#include "RooCintUtils.h"
#include "RooAddPdf.h"
#include "RooProdPdf.h"
#include "RooSimultaneous.h"
#include "RooFFTConvPdf.h"
#include "RooNumConvPdf.h"
#include "RooResolutionModel.h"
#include "RooProduct.h"
#include "RooAddition.h"
#include "RooChi2Var.h"
#include "RooNLLVar.h"
#include "RooRealSumPdf.h"
#include "RooConstVar.h"
#include "RooDerivative.h"
#include "TROOT.h"

using namespace RooFit ;
using namespace std ;

ClassImp(RooFactoryWSTool) 
;

RooFactoryWSTool* RooFactoryWSTool::_of = 0 ;
map<string,RooFactoryWSTool::IFace*>* RooFactoryWSTool::_hooks=0 ;

static Int_t init()
{
  RooFactoryWSTool::IFace* iface = new RooFactoryWSTool::SpecialsIFace ;

  // Operator p.d.f.s
  RooFactoryWSTool::registerSpecial("SUM",iface) ;
  RooFactoryWSTool::registerSpecial("RSUM",iface) ;
  RooFactoryWSTool::registerSpecial("ASUM",iface) ;
  RooFactoryWSTool::registerSpecial("PROD",iface) ;
  RooFactoryWSTool::registerSpecial("SIMUL",iface) ;
  RooFactoryWSTool::registerSpecial("EXPR",iface) ;
  RooFactoryWSTool::registerSpecial("FCONV",iface) ;
  RooFactoryWSTool::registerSpecial("NCONV",iface) ;

  // Operator functions
  RooFactoryWSTool::registerSpecial("sum",iface) ;
  RooFactoryWSTool::registerSpecial("prod",iface) ;
  RooFactoryWSTool::registerSpecial("expr",iface) ;
  RooFactoryWSTool::registerSpecial("nconv",iface) ;

  // Test statistics
  RooFactoryWSTool::registerSpecial("nll",iface) ;
  RooFactoryWSTool::registerSpecial("chi2",iface) ;
  RooFactoryWSTool::registerSpecial("profile",iface) ;

  // Integration and derivation
  RooFactoryWSTool::registerSpecial("int",iface) ;
  RooFactoryWSTool::registerSpecial("deriv",iface) ;
  RooFactoryWSTool::registerSpecial("cdf",iface) ;
  RooFactoryWSTool::registerSpecial("PROJ",iface) ;

  // Miscellaneous
  RooFactoryWSTool::registerSpecial("dataobs",iface) ;
  RooFactoryWSTool::registerSpecial("set",iface) ;

  return 0 ;
}
static Int_t dummy = init() ;


#ifndef _WIN32
#include <strings.h>
#else

static char *strtok_r(char *s1, const char *s2, char **lasts)
{
  char *ret;
  
  if (s1 == NULL)
    s1 = *lasts;
  while(*s1 && strchr(s2, *s1))
    ++s1;
  if(*s1 == '\0')
    return NULL;
  ret = s1;
  while(*s1 && !strchr(s2, *s1))
    ++s1;
  if(*s1)
    *s1++ = '\0';
  *lasts = s1;
  return ret;
}

#endif



//_____________________________________________________________________________
RooFactoryWSTool::RooFactoryWSTool(RooWorkspace& inws) : _ws(&inws), _errorCount(0)
{
  // Default constructor  
}



//_____________________________________________________________________________
RooFactoryWSTool::~RooFactoryWSTool() 
{
  // Destructor
}




//_____________________________________________________________________________
RooRealVar* RooFactoryWSTool::createVariable(const char* name, Double_t xmin, Double_t xmax)
{
  // Low-level factory interface for creating a RooRealVar with a given range and initial value

  // First check if variable already exists
  if (_ws->var(name)) {
    coutE(ObjectHandling) << "RooFactoryWSTool::createFactory() ERROR: variable with name '" << name << "' already exists" << endl ;
    logError() ;
    return 0 ;
  }

  // Create variable
  RooRealVar var(name,name,xmin,xmax) ;

  // Put in workspace
  if (_ws->import(var,Silence())) logError() ;

  return _ws->var(name) ;
}



//_____________________________________________________________________________
RooCategory* RooFactoryWSTool::createCategory(const char* name, const char* stateNameList) 
{
  // Low-level factory interface for creating a RooCategory with a given list of state names. The State name list
  // can be of the form 'name1,name2,name3' or of the form 'name1=id1,name2=id2,name3=id3'



  // Create variable
  RooCategory cat(name,name) ;

  // Add listed state names
  if (stateNameList) {
    char *tmp = new char[strlen(stateNameList)+1] ;
    strlcpy(tmp,stateNameList,strlen(stateNameList)+1) ;
    char* save ;
    char* tok = strtok_r(tmp,",",&save) ;
    while(tok) {
      char* sep = strchr(tok,'=') ;
      if (sep) {
        *sep = 0 ;
        Int_t id = atoi(sep+1) ;
        cat.defineType(tok,id) ;
        *sep = '=' ;
      } else {
        cat.defineType(tok) ;
      }
      tok = strtok_r(0,",",&save) ;
    }
    delete[] tmp ;
  }

  cat.setStringAttribute("factory_tag",Form("%s[%s]",name,stateNameList)) ;

  // Put in workspace
  if (_ws->import(cat,Silence())) logError() ;  

  return _ws->cat(name) ;
}



//_____________________________________________________________________________
RooAbsArg* RooFactoryWSTool::createArg(const char* className, const char* objName, const char* varList) 
{
  // Low-level factory interface for creating a RooAbsPdf of a given class with a given list of input variables
  // The variable list varList should be of the form "a,b,c" where the interpretation of the argument is
  // dependent on the p.d.f. Set and List arguments can be passed by substituting a single argument with
  // the form (a,b,c), i.e. one can set varList to "x,(a0,a1,a2)" to pass a RooAbsReal and a RooArgSet as arguments.
  
  // Find class in ROOT class table
  TClass* tc = resolveClassName(className) ;
  if (!tc) {
    coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR class " << className << " not found in factory alias table, nor in ROOT class table" << endl ;
    logError() ;
    return 0 ;
  }

  className = tc->GetName() ;

  // Check that class inherits from RooAbsPdf
  if (!tc->InheritsFrom(RooAbsArg::Class())) {
    coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR class " << className << " does not inherit from RooAbsArg" << endl ;
    logError() ;
    return 0 ;
  }


  _args.clear() ;
  char tmp[1024] ;
  strlcpy(tmp,varList,1024) ;
  char* p=tmp ;
  char* tok=tmp ;
  Int_t blevel(0) ;
  Bool_t litmode(kFALSE) ;
  while(*p) {

    // Keep track of opening and closing brackets
    if (*p=='{' || *p=='(' || *p=='[') blevel++ ;
    if (*p=='}' || *p==')' || *p==']') blevel-- ;

    // Keep track of string literals
    if (*p=='"' || *p=='\'') litmode = !litmode ;
    
    // If we encounter a comma at zero bracket level
    // finalize the current token as a completed argument
    // and start the next token
    if (!litmode && blevel==0 && ((*p)==',')) {
      *p = 0 ;
      _args.push_back(tok) ;
      tok = p+1 ;
    }

    p++ ;
  }
  _args.push_back(tok) ;


  // Try CINT interface
  pair<list<string>,unsigned int> ca = RooCintUtils::ctorArgs(className,_args.size()+2) ;
  if (ca.first.size()==0) {
    coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR no suitable constructor found for class " << className << endl ;
    logError() ;
    return 0 ;
  }  


  // Check if number of provided args is in valid range (add two to accomodate name and title strings)
  if (_args.size()+2<ca.second || _args.size()+2>ca.first.size()) {
    if (ca.second==ca.first.size()) {
      coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR number of arguments provided (" << _args.size() << ") for class is invalid, " << className 
                            << " expects " << ca.first.size()-2 << endl ;
      logError() ;
    } else {
      coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR number of arguments provided (" << _args.size() << ") for class is invalid " << className 
                            << " expect number between " << ca.second-2 << " and " << ca.first.size()-2 << endl ;
      logError() ;
    }
    return 0 ;
  }

  // Now construct CINT constructor spec, start with mandatory name and title args
  string cintExpr(Form("new %s(\"%s\",\"%s\"",className,objName,objName)) ;

  // Install argument in static data member to be accessed below through static CINT interface functions  
  _of = this ;

  
  try {
    Int_t i(0) ;
    list<string>::iterator ti = ca.first.begin() ; ti++ ; ti++ ;
    for (vector<string>::iterator ai = _args.begin() ; ai != _args.end() ; ai++,ti++,i++) {
      if ((*ti)=="RooAbsReal&" || (*ti)=="const RooAbsReal&") {
        RooFactoryWSTool::as_FUNC(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_FUNC(%d)",i) ;
      } else if ((*ti)=="RooAbsArg&" || (*ti)=="const RooAbsArg&") {
        RooFactoryWSTool::as_ARG(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_ARG(%d)",i) ;
      } else if ((*ti)=="RooRealVar&" || (*ti)=="const RooRealVar&") {
        RooFactoryWSTool::as_VAR(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_VAR(%d)",i) ;
      } else if ((*ti)=="RooAbsRealLValue&" || (*ti)=="const RooAbsRealLValue&") {
        RooFactoryWSTool::as_VARLV(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_VARLV(%d)",i) ;
      } else if ((*ti)=="RooCategory&" || (*ti)=="const RooCategory&") {
        RooFactoryWSTool::as_CAT(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_CAT(%d)",i) ;
      } else if ((*ti)=="RooAbsCategory&" || (*ti)=="const RooAbsCategory&") {
        RooFactoryWSTool::as_CATFUNC(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_CATFUNC(%d)",i) ;
      } else if ((*ti)=="RooAbsCategoryLValue&" || (*ti)=="const RooAbsCategoryLValue&") {
        RooFactoryWSTool::as_CATLV(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_CATLV(%d)",i) ;
      } else if ((*ti)=="RooAbsPdf&" || (*ti)=="const RooAbsPdf&") {
        RooFactoryWSTool::as_PDF(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_PDF(%d)",i) ;
      } else if ((*ti)=="RooResolutionModel&" || (*ti)=="const RooResolutionModel&") {
        RooFactoryWSTool::as_RMODEL(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_RMODEL(%d)",i) ;
      } else if ((*ti)=="RooAbsData&" || (*ti)=="const RooAbsData&") {
        RooFactoryWSTool::as_DATA(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_DATA(%d)",i) ;
      } else if ((*ti)=="RooDataSet&" || (*ti)=="const RooDataSet&") {
        RooFactoryWSTool::as_DSET(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_DSET(%d)",i) ;
      } else if ((*ti)=="RooDataHist&" || (*ti)=="const RooDataHist&") {
        RooFactoryWSTool::as_DHIST(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_DHIST(%d)",i) ;
      } else if ((*ti)=="const RooArgSet&") {
        RooFactoryWSTool::as_SET(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_SET(%d)",i) ;
      } else if ((*ti)=="const RooArgList&") {
        RooFactoryWSTool::as_LIST(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_LIST(%d)",i) ;
      } else if ((*ti)=="const char*") {
        RooFactoryWSTool::as_STRING(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_STRING(%d)",i) ;        
      } else if ((*ti)=="Int_t" || (*ti)=="int" || (*ti)=="Bool_t" || (*ti)=="bool") {
        RooFactoryWSTool::as_INT(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_INT(%d)",i) ;   
      } else if ((*ti)=="Double_t") {
        RooFactoryWSTool::as_DOUBLE(i) ;
        cintExpr += Form(",RooFactoryWSTool::as_DOUBLE(%d)",i) ;        
      } else if (RooCintUtils::isEnum(ti->c_str())) {             
        string qualvalue ;
        if (_args[i].find(Form("%s::",className)) != string::npos) {                
          qualvalue = _args[i].c_str() ;
        } else {        
          qualvalue =  Form("%s::%s",className,_args[i].c_str()) ;          
        }
        if (RooCintUtils::isValidEnumValue(ti->c_str(),qualvalue.c_str())) {
          cintExpr += Form(",(%s)%s",ti->c_str(),qualvalue.c_str()) ;
        } else {
          throw string(Form("Supplied argument %s does not represent a valid state of enum %s",_args[i].c_str(),ti->c_str())) ;
          }
      } else {
        // Check if generic object store has argument of given name and type
        TObject& obj = RooFactoryWSTool::as_OBJ(i) ;

        // Strip argument type to bare type (i.e. const X& -> X)
        string btype ;
        if (ti->find("const ")==0) {
          btype = ti->c_str()+6 ;
        } else {
          btype = *ti ;
        }
        if (btype.find("&")) {
          btype.erase(btype.size()-1,btype.size()) ;
        }

        // If btype if a typedef, substitute it by the true type name
        btype = RooCintUtils::trueName(btype.c_str()) ;

        if (obj.InheritsFrom(btype.c_str())) {
          cintExpr += Form(",(%s&)RooFactoryWSTool::as_OBJ(%d)",ti->c_str(),i) ;
        } else {
          throw string(Form("Required argument with name %s of type '%s' is not in the workspace",_args[i].c_str(),ti->c_str())) ;
        }
      }
    }
    cintExpr += ") ;" ;
  } catch (string err) {
    coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR constructing " << className << "::" << objName << ": " << err << endl ;
    logError() ;
    return 0 ;
  }
  
  cxcoutD(ObjectHandling) << "RooFactoryWSTool::createArg() Construct expression is " << cintExpr << endl ;

  // Call CINT to perform constructor call. Catch any error thrown by argument conversion method
  RooAbsArg* arg = (RooAbsArg*) gROOT->ProcessLineFast(cintExpr.c_str()) ;

  if (arg) {    
    if (string(className)=="RooGenericPdf") {
      arg->setStringAttribute("factory_tag",Form("EXPR::%s(%s)",objName,varList)) ;
    } else if (string(className)=="RooFormulaVar") {
      arg->setStringAttribute("factory_tag",Form("expr::%s(%s)",objName,varList)) ;
    } else {
      arg->setStringAttribute("factory_tag",Form("%s::%s(%s)",className,objName,varList)) ;
    }
    if (_ws->import(*arg,Silence())) logError() ;
    RooAbsArg* ret = _ws->arg(objName) ;
    delete arg ;
    return ret ;
  } else {
    coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR in CINT constructor call to create object" << endl ;
    logError() ;
    return 0 ;
  }
}



//_____________________________________________________________________________
vector<string> RooFactoryWSTool::ctorArgs(const char* /*className*/) 
{
  return vector<string>() ;
}




//_____________________________________________________________________________
RooAddPdf* RooFactoryWSTool::add(const char *objName, const char* specList, Bool_t recursiveCoefs)
{

  // Spec list is of form a*A,b*B,c*C,D [ *d]

  RooArgList pdfList ;
  RooArgList coefList ;
  RooArgList pdfList2 ;

  try {

    char buf[1024] ;
    strlcpy(buf,specList,1024) ;
    char* save ;
    char* tok = strtok_r(buf,",",&save) ;
    while(tok) {
      char* star=strchr(tok,'*') ;
      if (star) {
        *star=0 ;
        pdfList.add(asPDF(star+1)) ;
        coefList.add(asFUNC(tok)) ;
      } else {
        pdfList2.add(asPDF(tok)) ;
      }
      tok = strtok_r(0,",",&save) ;
    }
    pdfList.add(pdfList2) ;

  } catch (string err) {
    coutE(ObjectHandling) << "RooFactoryWSTool::add(" << objName << ") ERROR creating RooAddPdf: " << err << endl ;    
    logError() ;
    return 0 ;
  }
  
  RooAddPdf* pdf =  new RooAddPdf(objName,objName,pdfList,coefList,recursiveCoefs) ;
  pdf->setStringAttribute("factory_tag",Form("SUM::%s(%s)",objName,specList)) ;
  if (_ws->import(*pdf,Silence())) logError() ;
  return (RooAddPdf*) _ws->pdf(objName) ;
}


//_____________________________________________________________________________
RooRealSumPdf* RooFactoryWSTool::amplAdd(const char *objName, const char* specList)
{

  // Spec list is of form a*A,b*B,c*C,D [ *d]

  RooArgList amplList ;
  RooArgList coefList ;
  RooArgList amplList2 ;

  try {

    char buf[1024] ;
    strlcpy(buf,specList,1024) ;
    char* save ;
    char* tok = strtok_r(buf,",",&save) ;
    while(tok) {
      char* star=strchr(tok,'*') ;
      if (star) {
        *star=0 ;
        amplList.add(asFUNC(star+1)) ;
        coefList.add(asFUNC(tok)) ;
      } else {
        amplList2.add(asFUNC(tok)) ;
      }
      tok = strtok_r(0,",",&save) ;
    }
    amplList.add(amplList2) ;

  } catch (string err) {
    coutE(ObjectHandling) << "RooFactoryWSTool::add(" << objName << ") ERROR creating RooRealSumPdf: " << err << endl ;    
    logError() ;
    return 0 ;
  }
  
  RooRealSumPdf* pdf =  new RooRealSumPdf(objName,objName,amplList,coefList) ;
  pdf->setStringAttribute("factory_tag",Form("ASUM::%s(%s)",objName,specList)) ;
  if (_ws->import(*pdf,Silence())) logError() ;
  return (RooRealSumPdf*) _ws->pdf(objName) ;
}


//_____________________________________________________________________________
RooProdPdf* RooFactoryWSTool::prod(const char *objName, const char* pdfList) 
{
  _of = this ;

  // Separate conditional and non-conditional p.d.f terms
  RooLinkedList cmdList ;
  string regPdfList="{" ;
  char buf[1024] ;
  strlcpy(buf,pdfList,1024) ;
  char* save ;
  char* tok = strtok_r(buf,",",&save) ;
  while(tok) {
    char *sep = strchr(tok,'|') ;
    if (sep) {
      // Conditional term
      *sep=0 ;
      sep++ ;
      
      try {
        cmdList.Add(Conditional(asSET(tok),asSET(sep),kTRUE).Clone()) ;
      } catch (string err) {
        coutE(ObjectHandling) << "RooFactoryWSTool::prod(" << objName << ") ERROR creating RooProdPdf Conditional argument: " << err << endl ;
        logError() ;
        return 0 ;
      }
      
    } else {
      // Regular term
      if (regPdfList.size()>1) {
        regPdfList += "," ;
      }
      regPdfList += tok ;
    }
    tok = strtok_r(0,",",&save) ;
  }
  regPdfList += "}" ;
  
  RooProdPdf* pdf = 0 ;
  try {
    pdf = new RooProdPdf(objName,objName,asSET(regPdfList.c_str()),cmdList) ;
  } catch (string err) {
    coutE(ObjectHandling) << "RooFactoryWSTool::prod(" << objName << ") ERROR creating RooProdPdf input set of regular p.d.f.s: " << err << endl ;
    logError() ;
    pdf = 0 ;
  }
  cmdList.Delete() ;
  
  if (pdf) {
    pdf->setStringAttribute("factory_tag",Form("PROD::%s(%s)",objName,pdfList)) ;
    if (_ws->import(*pdf,Silence())) logError() ;
    delete pdf ;
    return (RooProdPdf*) _ws->pdf(objName) ;
  } else {
    return 0 ;
  }
}



//_____________________________________________________________________________
RooSimultaneous* RooFactoryWSTool::simul(const char* objName, const char* indexCat, const char* pdfMap)
{
  map<string,RooAbsPdf*> theMap ;
  // Add p.d.f. to index state mappings
  char buf[1024] ;
  strlcpy(buf,pdfMap,1024) ;
  char* save ;
  char* tok = strtok_r(buf,",",&save) ;
  while(tok) {
    char* eq = strchr(tok,'=') ;
    if (!eq) {
      coutE(ObjectHandling) << "RooFactoryWSTool::simul(" << objName << ") ERROR creating RooSimultaneous::" << objName 
                            << " expect mapping token of form 'state=pdfName', but found '" << tok << "'" << endl ;
      logError() ;
      return 0 ;
    } else {
      *eq = 0 ;

      try {
        theMap[tok] = &asPDF(eq+1) ;
      } catch ( string err ) {
        coutE(ObjectHandling) << "RooFactoryWSTool::simul(" << objName << ") ERROR creating RooSimultaneous: " << err << endl ;
        logError() ;
      }
    }
    tok = strtok_r(0,",",&save) ;
  }


  // Create simultaneous p.d.f.
  RooSimultaneous* pdf(0) ;
  try {
    pdf = new RooSimultaneous(objName,objName,theMap,asCATLV(indexCat)) ;
  } catch (string err) {
    coutE(ObjectHandling) << "RooFactoryWSTool::simul(" << objName << ") ERROR creating RooSimultaneous::" << objName << " " << err << endl ;
    logError() ;
  }

  // Import p.d.f into workspace
  pdf->setStringAttribute("factory_tag",Form("SIMUL::%s(%s,%s)",objName,indexCat,pdfMap)) ;
  if (_ws->import(*pdf,Silence())) logError() ;
  return (RooSimultaneous*) _ws->pdf(objName) ;
}




//_____________________________________________________________________________
RooAddition* RooFactoryWSTool::addfunc(const char *objName, const char* specList) 
{

  RooArgList sumlist1 ;
  RooArgList sumlist2 ;

  try {

    char buf[1024] ;
    strlcpy(buf,specList,1024) ;
    char* save ;
    char* tok = strtok_r(buf,",",&save) ;
    while(tok) {
      char* star=strchr(tok,'*') ;
      if (star) {
        *star=0 ;
        sumlist2.add(asFUNC(star+1)) ;
        sumlist1.add(asFUNC(tok)) ;
      } else {
        sumlist1.add(asFUNC(tok)) ;
      }
      tok = strtok_r(0,",",&save) ;
    }

  } catch (string err) {
    coutE(ObjectHandling) << "RooFactoryWSTool::addfunc(" << objName << ") ERROR creating RooAddition: " << err << endl ;
    logError() ;
    return 0 ;
  }

  if (sumlist2.getSize()>0 && (sumlist1.getSize()!=sumlist2.getSize())) {
    coutE(ObjectHandling) << "RooFactoryWSTool::addfunc(" << objName << ") ERROR creating RooAddition: syntax error: either all sum terms must be products or none" << endl ;
    logError() ;
    return 0 ;
  }


  RooAddition* sum ;
  if (sumlist2.getSize()>0) {
    sum = new RooAddition(objName,objName,sumlist1,sumlist2) ;
  } else {
    sum = new RooAddition(objName,objName,sumlist1) ;
  }

  sum->setStringAttribute("factory_tag",Form("sum::%s(%s)",objName,specList)) ;
  if (_ws->import(*sum,Silence())) logError() ;
  delete sum ;
  return (RooAddition*) _ws->pdf(objName) ;
  
}




//_____________________________________________________________________________
RooProduct* RooFactoryWSTool::prodfunc(const char *objName, const char* pdfList) 
{
  return (RooProduct*) createArg("RooProduct",objName,Form("{%s}",pdfList)) ;
}





//_____________________________________________________________________________
RooAbsArg* RooFactoryWSTool::process(const char* expr) 
{
  // Process high-level object creation syntax
  // Accepted forms of syntax are
  //
  //
  // Creating variables
  //
  // x[-10,10]             -  Create variable x with given range and put it in workspace
  // x[3,-10,10]           -  Create variable x with given range and initial value and put it in workspace
  // x[3]                  -  Create variable x with given constant value
  //
  // <numeric literal>     - Numeric literal expressions (0.5, -3 etc..) are converted to a RooConst(<numeric literal>) 
  //                         where ever a RooAbsReal or RooAbsArg argument is expected
  //
  // Creating categories
  //
  // c[lep,kao,nt1,nt2]    -  Create category c with given state names
  // tag[B0=1,B0bar=-1]    -  Create category tag with given state names and index assignments
  //
  //
  // Creating functions and p.d.f.s
  //
  // MyPdf::g(x,m,s)       - Create p.d.f or function of type MyPdf with name g with argument x,m,s
  //                         Interpretation and number of arguments are mapped to the constructor arguments of the class 
  //                         (after the name and title).
  //
  // MyPdf(x,m,s)          - As above, but with an implicitly defined (unique) object name
  // 
  //
  // Creating sets and lists (to be used as inputs above)
  //
  // {a,b,c}               - Create RooArgSet or RooArgList (as determined by context) from given contents
  //
  //
  //
  // Objects that are not created, are assumed to exist in the workspace
  // Object creation expressions as shown above can be nested, e.g. one can do
  //
  //   RooGaussian::g(x[-10,10],m[0],3) 
  //
  // to create a p.d.f and its variables in one go. This nesting can be applied recursively e.g.
  //
  //   SUM::model( f[0.5,0,1] * RooGaussian::g( x[-10,10], m[0], 3] ),
  //                            RooChebychev::c( x, {a0[0.1],a1[0.2],a2[-0.3]} ))
  //
  // creates the sum of a Gaussian and a Chebychev and all its variables
  //
  //
  // A seperate series of operator meta-type exists to simplify the construction of composite expressions
  // meta-types in all capitals (SUM) create p.d.f.s, meta types in lower case (sum) create
  // functions.
  //
  //
  // SUM::name(f1*pdf1,f2*pdf2,pdf3]  -- Create sum p.d.f name with value f1*pdf1+f2*pdf2+(1-f1-f2)*pdf3
  // RSUM::name(f1*pdf1,f2*pdf2,pdf3] -- Create recursive sum p.d.f. name with value f1*pdf1 + (1-f1)(f2*pdf2 + (1-f2)pdf3)
  // ASUM::name(f1*amp1,f2*amp2,amp3] -- Create sum p.d.f. name with value f1*amp1+f2*amp2+(1-f1-f2)*amp3 where amplX are amplitudes of type RooAbsReal
  // sum::name(a1,a2,a3]              -- Create sum function with value a1+a2+a3
  // sum::name(a1*b1,a2*b2,a3*b 3]    -- Create sum function with value a1*b1+a2*b2+a3*b3
  //
  // PROD::name(pdf1,pdf2]            -- Create product of p.d.f with 'name' with given input p.d.fs
  // PROD::name(pdf1|x,pdf2]          -- Create product of conditional p.d.f. pdf1 given x and pdf2
  // prod::name(a,b,c]                -- Create production function with value a*b*c
  //
  // SIMUL::name(cat,a=pdf1,b=pdf2]   -- Create simultaneous p.d.f index category cat. Make pdf1 to state a, pdf2 to state b
  //
  // EXPR::name('expr',var,...]       -- Create an generic p.d.f that interprets the given expression 
  // expr::name('expr',var,...]       -- Create an generic function that interprets the given expression 
  //
  //
  // The functionality of high level object creation tools like RooSimWSTool, RooCustomizer and RooClassFactory
  // is also interfaced through meta-types in the factory
  //
  //
  // Interface to RooSimWSTool
  //
  // SIMCLONE::name( modelPdf, $ParamSplit(...), 
  //                 $ParamSplitConstrained(...), $Restrict(...) ]            -- Clone-and-customize modelPdf according to ParamSplit and ParamSplitConstrained()
  //                                                                             specifications and return a RooSimultaneous p.d.f. of all built clones
  //
  // MSIMCLONE::name( masterIndex, 
  //                  $AddPdf(mstate1, modelPdf1, $ParamSplit(...)), 
  //                  $AddPdf(mstate2,modelPdf2),...) ]                       -- Clone-and-customize multiple models (modelPdf1,modelPdf2) according to ParamSplit and 
  //                                                                             ParamSplitConstrained() specifications and return a RooSimultaneous p.d.f. of all built clones,
  //                                                                             using the specified master index to map prototype p.d.f.s to master states
  // Interface to RooCustomizer
  //
  // EDIT::name( orig, substNode=origNode), ... ]                             -- Create a clone of input object orig, with the specified replacements operations executed
  //
  //
  // Interface to RooClassFactory
  //
  // CEXPR::name('expr',var,...]       -- Create an custom compiled p.d.f that evaluates the given expression 
  // cexpr::name('expr',var,...]       -- Create an custom compiled function that evaluates the given expression 
  //
  //
  // $MetaType(...)        - Meta argument that does not result in construction of an object but is used logically organize 
  //                         input arguments in certain operator p.d.f. constructions. The defined meta arguments are context dependent.
  //
  //                         The only meta argument that is defined globally is $Alias(typeName,aliasName) to  
  //                         define aliases for type names. For the definition of meta arguments in operator p.d.f.s
  //                         see the definitions below


  // First perform basic syntax check
  if (checkSyntax(expr)) {
    return 0 ;
  }

  // Allocate work buffer
  char* buf = new char[strlen(expr)+1] ;

  // Copy to buffer while absorbing white space and newlines
  char* buftmp = buf ;
  while(*expr) {
    if (!isspace(*expr)) {
      *buftmp = *expr ;
      buftmp++ ;
    }
    expr++ ;
  }
  *buftmp=0 ;


  // Clear error count and start a transaction in the workspace
  clearError() ;
  ws().startTransaction() ;

  // Process buffer
  string out ;
  try {
    out = processExpression(buf) ;
  } catch (string error) {
    coutE(ObjectHandling) << "RooFactoryWSTool::processExpression() ERROR in parsing: " << error << endl ;
    logError() ;
  }

  // If there were no errors commit the transaction, cancel it otherwise
  if (errorCount()>0) {
    coutE(ObjectHandling) << "RooFactoryWSTool::processExpression() ERRORS detected, transaction to workspace aborted, no objects committed" << endl ;
    ws().cancelTransaction() ;
  } else {
    ws().commitTransaction() ;
  }

  
  // Delete buffer 
  delete[] buf ;

  return out.size() ? ws().arg(out.c_str()) : 0 ;
}




//_____________________________________________________________________________
std::string RooFactoryWSTool::processExpression(const char* token) 
{
  // Process a single high-level expression or list of
  // expressions. The returned string a the reduced expression where
  // all inline object creations have been executed and substituted
  // with the name of the created object
  //
  // e.g. 'RooGaussian::g(x,m,s)' --> 'g'
  //      '{x(-10,10),s}          --> '{x,s}'

  // Delegate handling to list processor if token starts with {, otherwise
  // call single expression processor
  if (string(token).find("$Alias(")==0) {
    processAliasExpression(token) ;
  }

  if (token[0]=='{') {    
    // Process token as list if it starts with '{'
    return processListExpression(token) ;
  } else {
    // Process token as single item otherwise
    return processCompositeExpression(token) ;
  }
}



//_____________________________________________________________________________
std::string RooFactoryWSTool::processCompositeExpression(const char* token) 
{
  // Process a single composite expression 
  //
  // e.g. 'A=RooGaussian::g[x,m,s]' --> 'A=g'
  // e.g. 'f[0,1]*RooGaussian::g[x,m,s]' --> 'f*g'
  // e.g. 'RooGaussian::g(x,y,s)|x' --> g|x'
  // e.g. '$MetaArg(RooGaussian::g[x,m,s],blah)' --> '$MetaArg(g,blah)'

  // Allocate and fill work buffer
  char* buf_base = new char[strlen(token)+1] ;   
  char* buf = buf_base ;
  strlcpy(buf,token,strlen(token)+1) ;
  char* p = buf ;

  list<string> singleExpr ;
  list<char> separator ;
  Int_t blevel(0) ;
  Bool_t litmode(kFALSE) ;
  while(*p) {

    // Keep track of opening and closing brackets
    if (*p=='{' || *p=='(' || *p=='[') blevel++ ;
    if (*p=='}' || *p==')' || *p==']') blevel-- ;
    
    // Keep track of string literals
    if (*p=='"' || *p=='\'') litmode = !litmode ;

    // If we are zero-bracket level and encounter a |, store
    // the remainder of the string as suffix and exit loop
    if (!litmode && blevel==0 && ( (*p)=='=' || (*p) == '|' || (*p) == '*')) {
      separator.push_back(*p) ;
      *p=0 ;
      singleExpr.push_back(buf) ;
      buf = p+1 ;
    }
    p++ ;
  }
  if (*buf) {
    singleExpr.push_back(buf) ;
  }
  if (singleExpr.size()==1) {    
    string ret = processSingleExpression(token) ;
    delete[] buf_base ;
    return ret ;
  }

  string ret ;
  list<char>::iterator ic = separator.begin() ;
  for (list<string>::iterator ii = singleExpr.begin() ; ii!=singleExpr.end() ; ii++) {
    ret += processSingleExpression(ii->c_str()) ;
    if (ic != separator.end()) {
      ret += *ic ;
      ic++ ;
    }
  }

  delete[] buf_base ;
  return ret ;
}



//_____________________________________________________________________________
std::string RooFactoryWSTool::processSingleExpression(const char* arg) 
{
  // Process a single high-level expression. The returned string a the reduced
  // expression where all inline object creations have been executed and substituted
  // with the name of the created object
  //
  // e.g. 'RooGaussian::g(x,m,s)' --> 'g'
  // e.g. 'x[-10,10]' --> 'x'

  // Handle empty strings here
  if (strlen(arg)==0) {
    return string("") ;
  }

  // Handle string literal case
  if (arg[0]=='\'' || arg[0]=='"') {
    return string(arg) ;
  }

  // Allocate and fill work buffer
  char* buf = new char[strlen(arg)+1] ; 
  strlcpy(buf,arg,strlen(arg)+1) ;
  char* bufptr = buf ;

  string func,prefix ;
  vector<string> args ;

  // Process token into arguments
  char* save ;
  char* tmpx = strtok_r(buf,"([",&save) ;
  func = tmpx ? tmpx : "" ;
  char* p = strtok_r(0,"",&save) ;
  
  // Return here if token is fundamental
  if (!p) {
    delete[] buf ;
    return arg ;
  }
    

  char* tok = p ;
  Int_t blevel=0 ;
  Bool_t litmode(kFALSE) ;
  while(*p) {

    // Keep track of opening and closing brackets
    if (*p=='{' || *p=='(' || *p=='[') blevel++ ;
    if (*p=='}' || *p==')' || *p==']') blevel-- ;

    // Keep track of string literals
    if (*p=='"' || *p=='\'') litmode = !litmode ;


    // If we encounter a comma at zero bracket level
    // finalize the current token as a completed argument
    // and start the next token
    if (!litmode && blevel==0 && ((*p)==',')) {
      *p = 0 ;
      args.push_back(tok) ;
      tok = p+1 ;
    }

    p++ ;
  }

  // If the last character was a closing bracket, kill
  // it in the buffer
  if (p>bufptr && (*(p-1)==')'||*(p-1)==']')) {
    *(p-1)=0 ;
  }

  // Finalize last token as argument
  string tmp = tok ;
  
  // If there is a suffix left in the work buffer attach it to 
  // this argument
  p = strtok_r(0,"",&save) ;
  if (p) tmp += p ;
  args.push_back(tmp) ;

  // Delete the work buffer
  delete[] buf ;

  // If function contains :: then call createArg to process this arg, otherwise
  // call createVariable
  string ret ;

  // Determine type of leading bracket
  char lb = ' ' ;  
  for(const char* pp=arg ; *pp!=0 ; pp++) {
    if (*pp=='(' || *pp=='[' || *pp=='{') {
      lb = *pp ;
      break ;
    }
  }

  if (strstr(func.c_str(),"::")) {
    if (lb=='(') {
      // Create function argument with instance name
      ret= processCreateArg(func,args) ;
    } else {
      coutE(ObjectHandling) << "RooFactoryWSTool::processSingleExpression(" << arg << "): ERROR: Syntax error: Class::Instance must be followed by (...)" << endl ;
      logError() ;
    }
  } else if (func[0]!='$'){    
    if (lb=='[') {
      // Create variable argument
      ret= processCreateVar(func,args) ;
    } else if (lb=='(') {

      // Create function argument with autoname
      string autoname ;
      if (!_autoNamePrefix.empty()) {
        // If we're inside a function creation call to a higher level object, use its
        // name as base for the autoname
        autoname = (Form("%s::%s",func.c_str(),_autoNamePrefix.top().c_str())) ;
      } else {
        // Otherwise find a free global_%d name
        static Int_t globCounter = 0 ;
        while(true) {
          autoname = Form("gobj%d",globCounter) ;
          globCounter++ ;
          if (!ws().arg(autoname.c_str())) {
            break ;
          }
        }
        autoname = Form("%s::%s",func.c_str(),autoname.c_str()) ;
      }
      ret= processCreateArg(autoname,args) ;
    } else {
      coutE(ObjectHandling) << "RooFactoryWSTool::processSingleExpression(" << arg << "): ERROR: Syntax error: expect either Class(...) or Instance[...]" << endl ;
      logError() ;
    }
  } else {
    if (lb=='(') {
      // Process meta function (compile arguments, but not meta-function itself)
      ret= processMetaArg(func,args) ;
    } else {
      coutE(ObjectHandling) << "RooFactoryWSTool::processSingleExpression(" << arg << "): ERROR: Syntax error: $MetaClass must be followed by (...)" << endl ;
      logError() ;
    }
  }

  // Return reduced token with suffix
  return ret ;
}


//_____________________________________________________________________________
string RooFactoryWSTool::processListExpression(const char* arg) 
{
  // Process a list of high-level expression. The returned string a the reduced
  // expression list where all inline object creations have been executed and substituted
  // with the name of the created object
  //
  // E.g.   '{x(-10,10),s}  --> '{x,s}'

  // Allocate and fill work buffer
  char* buf = new char[strlen(arg)+1] ;
  strlcpy(buf,arg,strlen(arg)+1) ;

  vector<string> args ;

  // Start running pointer at position 1 to skip opening bracket
  char* tok = buf+1 ;
  char* p = buf+1 ;
  
  // Processing look
  Int_t level(0) ;
  while(*p) {

    // Track bracketing level
    if (*p=='{' || *p=='(' || *p=='[') level++ ;
    if (*p=='}' || *p==')' || *p==']') level-- ;

    
    // If we encounter a comma at zero bracket level
    // finalize the current token as a completed argument
    // and start the next token
    if (level==0 && ((*p)==',')) {
      *p = 0 ;
      args.push_back(tok) ;
      tok = p+1 ;
    }

    p++ ;
  }

  // Finalize token as last argument
  if (p>buf && *(p-1)=='}') {
    *(p-1)=0 ;
  }
  args.push_back(tok) ;

  // Delete work buffer
  delete[] buf ;

  // Process each argument in list and construct reduced 
  // expression to be returned
  string ret("{") ;
  vector<string>::iterator iter = args.begin() ;
  Int_t i(0) ;
  while(iter!= args.end()) {
    if (strlen(ret.c_str())>1) ret += "," ;
    if (!_autoNamePrefix.empty()) {
      _autoNamePrefix.push(Form("%s%d",_autoNamePrefix.top().c_str(),i+1)) ;
    }
    ret += processSingleExpression(iter->c_str()) ;
    if (!_autoNamePrefix.empty()) {
      _autoNamePrefix.pop() ;
    }
    iter++ ;
    i++ ;
  }
  ret += "}" ;

  return ret ;
}



//_____________________________________________________________________________
string RooFactoryWSTool::processAliasExpression(const char* token) 
{
  // Parse token
  vector<string> args = splitFunctionArgs(token) ;
  if (args.size()!=2) {
    coutE(ObjectHandling) << "RooFactorWSTool::processAliasExpression() ERROR $Alias() takes exactly two arguments, " << args.size() << " args found" << endl ;
    logError() ;
    return string() ;
  }

  // Insert alias in table
  _typeAliases[args[1]] = args[0] ;
  
  return string() ;
}




//_____________________________________________________________________________
TClass* RooFactoryWSTool::resolveClassName(const char* className) 
{

  // First do recursive alias expansion
  while (true) {
    map<string,string>::iterator item = _typeAliases.find(className) ;
    
    // If an alias is found, recurse
    if (item != _typeAliases.end()) {
      className = item->second.c_str() ;
    } else {
      break ;
    }
  }

  // Now find dealiased class in ROOT class table
  TClass* tc =  TClass::GetClass(className,kTRUE,kTRUE) ;

  // If its not there, try prefixing with Roo
  if (!tc) {
    tc = TClass::GetClass(Form("Roo%s",className)) ;
    if (!tc) {
      coutE(ObjectHandling) << "RooFactoryWSTool::createArg() ERROR class " << className << " not defined in ROOT class table" << endl ;
      logError() ;
      return 0 ;
    } 
  }
  return tc ;
}



//_____________________________________________________________________________
string RooFactoryWSTool::varTag(string& func, vector<string>& args)
{
  string ret ;
  ret += func ;
  ret += "[" ;
  for (vector<string>::iterator iter = args.begin() ; iter!=args.end() ; ++iter) {
    if (iter!=args.begin()) {
      ret += "," ;
    }
    ret += *iter ;
  }
  ret += "]" ;
  return ret ;
}




//_____________________________________________________________________________
string RooFactoryWSTool::processCreateVar(string& func, vector<string>& args)
{  
  // Glue function between high-level syntax and low-level factory call to createVariable:
  // Process a parsed call to create a variable named 'func'
  // 
  // If initial token is non-numeric, a RooCategory will be created, and the args are interpreted
  // as either state names or 'name=id' assignments. Otherwise a RooRealvar is created and the
  // arg list is interpreted as follows:
  // If list has two args, these are interpreted as xmin,xmax
  // If list has three args, these are interpreted as xinit,xmin,xmax
  // If list has one arg, this is interpreted as xinit and the variable is set as constant


  // Determine if first arg is numeric
  string first = *(args.begin()) ;
  if (isdigit(first[0]) || first[0]=='.' || first[0]=='+' || first[0]=='-') {
    
    // Create a RooRealVar
    vector<string>::iterator ai = args.begin() ;
    if (args.size()==1) {

      // One argument, create constant variable with given value
      Double_t xinit = atof((ai)->c_str()) ;
      cxcoutD(ObjectHandling) << "CREATE variable " << func << " xinit = " << xinit << endl ;
      RooRealVar tmp(func.c_str(),func.c_str(),xinit) ;
      tmp.setStringAttribute("factory_tag",varTag(func,args).c_str()) ;
      if (_ws->import(tmp,Silence())) {
        logError() ;
      }

    } else if (args.size()==2) {

      // Two arguments, create variable with given range
      Double_t xlo = atof((ai++)->c_str()) ;
      Double_t xhi = atof(ai->c_str()) ;
      cxcoutD(ObjectHandling) << "CREATE variable " << func << " xlo = " << xlo << " xhi = " << xhi << endl ;
      RooRealVar tmp(func.c_str(),func.c_str(),xlo,xhi) ;
      tmp.setStringAttribute("factory_tag",varTag(func,args).c_str()) ;
      if (_ws->import(tmp,Silence())) {
        logError() ;
      }

    } else if (args.size()==3) {

      // Three arguments, create variable with given initial value and range
      Double_t xinit = atof((ai++)->c_str()) ;
      Double_t xlo = atof((ai++)->c_str()) ;
      Double_t xhi = atof(ai->c_str()) ;
      cxcoutD(ObjectHandling) << "CREATE variable " << func << " xinit = " << xinit << " xlo = " << xlo << " xhi = " << xhi << endl ;
      RooRealVar tmp(func.c_str(),func.c_str(),xinit,xlo,xhi) ;
      tmp.setStringAttribute("factory_tag",varTag(func,args).c_str()) ;
      if (_ws->import(tmp,Silence())) {
        logError() ;
      }
    }
  } else {

    // Create a RooAbsCategory
    string allStates ;
    for (vector<string>::iterator ai = args.begin() ; ai!=args.end() ; ai++) {
      if (allStates.size()>0) {
        allStates += "," ;
      }
      allStates += *ai ;
    }
    createCategory(func.c_str(),allStates.c_str()) ;

  }
  return func ;
}
      

//_____________________________________________________________________________
string RooFactoryWSTool::processCreateArg(string& func, vector<string>& args)
{  
  // Glue function between high-level syntax and low-level factory call to createArg:
  // Process a parsed call to create a p.d.f named func
  //
  // The func arg is interpreted as ClassName::ObjectName and the arglist is passed
  // verbatim to createArg. The received arglist is expected to be fully reduced (i.e.
  // all inline object creations must have been compiled)

  // Allocate and fill work buffer
  char buf[1024] ;
  strlcpy(buf,func.c_str(),1024) ;

  // Split function part in class name and instance name
  char* save ;
  const char *className = strtok_r(buf,":",&save) ;
  const char *instName = strtok_r(0,":",&save) ;
  if (!className) className = "";
  if (!instName) instName = "" ;

  // Concatenate list of args into comma separated string
  char pargs[1024] ;
  pargs[0] = 0 ;
  vector<string>::iterator iter = args.begin() ;
  vector<string> pargv ;
  Int_t iarg(0) ;
  while(iter!=args.end()) {
    if (strlen(pargs)>0) strlcat(pargs,",",1024) ;
    _autoNamePrefix.push(Form("%s_%d",instName,iarg+1)) ;
    string tmp = processExpression(iter->c_str()) ;
    _autoNamePrefix.pop() ;
    strlcat(pargs,tmp.c_str(),1024) ;
    pargv.push_back(tmp) ;
    iter++ ;
    iarg++ ;
  }

  // Look up if func is a special
  for (map<string,IFace*>::iterator ii=hooks().begin() ; ii!=hooks().end() ; ii++) {
  }
  if (hooks().find(className) != hooks().end()) {
    IFace* iface = hooks()[className] ;
    return iface->create(*this, className,instName,pargv) ;
  }

  createArg(className,instName,pargs) ;

  return string(instName) ;
}
      


//_____________________________________________________________________________
std::string RooFactoryWSTool::processMetaArg(std::string& func, std::vector<std::string>& args) 
{
  // Concatenate list of args into comma separated string
  char pargs[1024] ;
  pargs[0] = 0 ;
  vector<string>::iterator iter = args.begin() ;
  vector<string> pargv ;
  while(iter!=args.end()) {
    if (strlen(pargs)>0) strlcat(pargs,",",1024) ;
    string tmp = processExpression(iter->c_str()) ;
    strlcat(pargs,tmp.c_str(),1024) ;
    pargv.push_back(tmp) ;
    iter++ ;
  }

  string ret = func+"("+pargs+")" ;  
  return ret ;
}




//_____________________________________________________________________________
vector<string> RooFactoryWSTool::splitFunctionArgs(const char* funcExpr) 
{
  // Allocate and fill work buffer
  char* buf = new char[strlen(funcExpr)+1] ; 
  strlcpy(buf,funcExpr,strlen(funcExpr)+1) ;
  char* bufptr = buf ;

  string func ;
  vector<string> args ;

  // Process token into arguments
  char* save ;
  char* tmpx = strtok_r(buf,"(",&save) ;  
  func = tmpx ? tmpx : "" ;
  char* p = strtok_r(0,"",&save) ;
  
  // Return here if token is fundamental
  if (!p) {
    delete[] buf ;
    return args ;
  }

  char* tok = p ;
  Int_t blevel=0 ;
  Bool_t litmode(kFALSE) ;
  while(*p) {

    // Keep track of opening and closing brackets
    if (*p=='{' || *p=='(' || *p=='[') blevel++ ;
    if (*p=='}' || *p==')' || *p==']') blevel-- ;

    // Keep track of string literals
    if (*p=='"' || *p=='\'') litmode = !litmode ;


    // If we encounter a comma at zero bracket level
    // finalize the current token as a completed argument
    // and start the next token
    if (!litmode && blevel==0 && ((*p)==',')) {
      *p = 0 ;
      args.push_back(tok) ;
      tok = p+1 ;
    }

    p++ ;
  }

  // If the last character was a closing bracket, kill
  // it in the buffer
  if (p>bufptr && *(p-1)==')') {
    *(p-1)=0 ;
  }

  // Finalize last token as argument
  string tmp = tok ;
  
  // If there is a suffix left in the work buffer attach it to 
  // this argument
  p = strtok_r(0,"",&save) ;
  if (p) tmp += p ;
  args.push_back(tmp) ;

  // Delete the work buffer
  delete[] buf ;
  
  return args ;
}





//_____________________________________________________________________________
Bool_t RooFactoryWSTool::checkSyntax(const char* arg) 
{
  // Perform basic syntax on given factory expression. If function returns
  // true syntax errors are found.
  
  // Count parentheses
  Int_t nParentheses(0), nBracket(0), nAccolade(0) ;
  const char* ptr = arg ;
  while(*ptr) {
    if (*ptr=='(') nParentheses++ ;
    if (*ptr==')') nParentheses-- ;
    if (*ptr=='[') nBracket++ ;
    if (*ptr==']') nBracket-- ;
    if (*ptr=='{') nAccolade++ ;
    if (*ptr=='}') nAccolade-- ;
    ptr++ ;
  }
  if (nParentheses!=0) {
    coutE(ObjectHandling) << "RooFactoryWSTool::checkSyntax ERROR non-matching '" << (nParentheses>0?"(":")") << "' in expression" << endl ;
    return kTRUE ;
  }
  if (nBracket!=0) {
    coutE(ObjectHandling) << "RooFactoryWSTool::checkSyntax ERROR non-matching '" << (nBracket>0?"[":"]") << "' in expression" << endl ;
    return kTRUE ;
  }
  if (nAccolade!=0) {
    coutE(ObjectHandling) << "RooFactoryWSTool::checkSyntax ERROR non-matching '" << (nAccolade>0?"{":"}") << "' in expression" << endl ;
    return kTRUE ;
  }
  return kFALSE ;
}



//_____________________________________________________________________________
void RooFactoryWSTool::checkIndex(UInt_t idx) 
{
  if (idx>_of->_args.size()-1) {
    throw string(Form("Need argument number %d, but only %d args are provided",idx,(Int_t)_of->_args.size())) ;
  }
}



//_____________________________________________________________________________
RooAbsArg& RooFactoryWSTool::asARG(const char* arg) 
  {
  // CINT constructor interface, return constructor string argument #idx as RooAbsArg reference found in workspace

  // If arg is a numeric string, make a RooConst() of it here
  if (arg[0]=='.' || arg[0]=='+' || arg[0] == '-' || isdigit(arg[0])) {
    return RooConst(atof(arg)) ;
  }

  // Otherwise look it up by name in the workspace
  RooAbsArg* rarg = ws().arg(arg) ;
  if (!rarg) {
    throw string(Form("RooAbsArg named %s not found",arg)) ;
  }
  return *rarg ;
}



//_____________________________________________________________________________
RooAbsReal& RooFactoryWSTool::asFUNC(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooAbsReal reference found in workspace

  // If arg is a numeric string, make a RooConst() of it here
  if (arg[0]=='.' || arg[0]=='+' || arg[0] == '-' || isdigit(arg[0])) {
    return RooConst(atof(arg)) ;
  }

  RooAbsArg* rarg = ws().arg(arg) ;
  if (!rarg) {
    throw string(Form("RooAbsReal named %s not found",arg)) ;
  }
  RooAbsReal* real = dynamic_cast<RooAbsReal*>(rarg) ;
  if (!real) {
    throw string(Form("Object named %s is not of type RooAbsReal",arg)) ;
  }
  return *real ;
}



//_____________________________________________________________________________
RooAbsRealLValue& RooFactoryWSTool::asVARLV(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooAbsRealLValue reference found in workspace

  // If arg is a numeric string, throw error as lvalue is required
  if (arg[0]=='.' || arg[0]=='+' || arg[0] == '-' || isdigit(arg[0])) {
    throw string(Form("Numeric literal provided for argument (%s), but lvalue is required",arg)) ;
  }

  RooAbsArg* rarg = ws().arg(arg) ;
  if (!rarg) {
    throw string(Form("RooAbsRealLValue named %s not found",arg)) ;
  }
  RooAbsRealLValue* reallv = dynamic_cast<RooAbsRealLValue*>(rarg) ;
  if (!reallv) {
    throw string(Form("Object named %s is not of type RooAbsRealLValue",arg)) ;
  }
  return *reallv ;
}



//_____________________________________________________________________________
RooRealVar& RooFactoryWSTool::asVAR(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooRealVar reference found in workspace

  RooRealVar* var = ws().var(arg) ;
  if (!var) {
    throw string(Form("RooRealVar named %s not found",arg)) ;
  }
  return *var ;
}




//_____________________________________________________________________________
RooAbsPdf& RooFactoryWSTool::asPDF(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooAbsPdf reference found in workspace

  RooAbsPdf* pdf = ws().pdf(arg) ;
  if (!pdf) {
    throw string(Form("RooAbsPdf named %s not found",arg)) ;
  }
  return *pdf ;
}




//_____________________________________________________________________________
RooResolutionModel& RooFactoryWSTool::asRMODEL(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooResolutionModel reference found in workspace

  RooAbsArg* rarg = ws().arg(arg) ;
  if (!rarg) {
    throw string(Form("RooResolutionModel named %s not found",arg)) ;
  }
  RooResolutionModel * rmodel = dynamic_cast<RooResolutionModel*>(rarg) ;
  if (!rmodel) {
    throw string(Form("Object named %s is not of type RooResolutionModel",arg)) ;
  }
  return *rmodel ;
}




//_____________________________________________________________________________
RooAbsCategory& RooFactoryWSTool::asCATFUNC(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooAbsCategory reference found in workspace

  RooAbsArg* rarg = ws().arg(arg) ;
  if (!rarg) {
    throw string(Form("RooAbsCategory named %s not found",arg)) ;
  }
  RooAbsCategory* catf = dynamic_cast<RooAbsCategory*>(rarg) ;
  if (!catf) {
    throw string(Form("Object named %s is not of type RooAbsCategory",arg)) ;
  }
  return *catf ;
}



//_____________________________________________________________________________
RooAbsCategoryLValue& RooFactoryWSTool::asCATLV(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooAbsCategoryLValue reference found in workspace

  RooAbsArg* rarg = ws().arg(arg) ;
  if (!rarg) {
    throw string(Form("RooAbsCategoryLValue named %s not found",arg)) ;
  }

  RooAbsCategoryLValue* catlv = dynamic_cast<RooAbsCategoryLValue*>(rarg) ;
  if (!catlv) {
    throw string(Form("Object named %s is not of type RooAbsCategoryLValue",arg)) ;
  }
  return *catlv ;
}



//_____________________________________________________________________________
RooCategory& RooFactoryWSTool::asCAT(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooCategory reference found in workspace

  RooCategory* cat = ws().cat(arg) ;
  if (!cat) {
    throw string(Form("RooCategory named %s not found",arg)) ;
  }
  return *cat ;
}





//_____________________________________________________________________________
RooArgSet RooFactoryWSTool::asSET(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooArgSet of objects found in workspace

  char tmp[1024] ;
  strlcpy(tmp,arg,1024) ;

  RooArgSet s ;
  
  // If given object is not of {,,,} form, interpret given string as name of defined set
  if (arg[0]!='{') {
    const RooArgSet* defSet = ws().set(arg) ;
    if (defSet) {
      s.add(*defSet) ;
      return s ;
    }
  }

  char* save ;
  char* tok = strtok_r(tmp,",{}",&save) ;
  while(tok) {

    // If arg is a numeric string, make a RooConst() of it here
    if (tok[0]=='.' || tok[0]=='+' || tok[0] == '-' || isdigit(tok[0])) {
      s.add(RooConst(atof(tok))) ;
    } else {
      RooAbsArg* aarg = ws().arg(tok) ;
      if (aarg) {
        s.add(*aarg)  ;
      } else {
        throw string(Form("RooAbsArg named %s not found",tok)) ;
      }
    }
    tok = strtok_r(0,",{}",&save) ;
  }

  return s ;
}



//_____________________________________________________________________________
RooArgList RooFactoryWSTool::asLIST(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooArgList of objects found in workspace

  char tmp[1024] ;
  strlcpy(tmp,arg,1024) ;

  RooArgList l ;
  char* save ;
  char* tok = strtok_r(tmp,",{}",&save) ;
  while(tok) {

    // If arg is a numeric string, make a RooConst() of it here
    if (tok[0]=='.' || tok[0]=='+' || tok[0] == '-' || isdigit(tok[0])) {
      l.add(RooConst(atof(tok))) ;
    } else {
      RooAbsArg* aarg = ws().arg(tok) ;
      if (aarg) {
        l.add(*aarg)  ;
      } else {
        throw string(Form("RooAbsArg named %s not found",tok)) ;
      }
    }
    tok = strtok_r(0,",{}",&save) ;
  }
  
  return l ;
}



//_____________________________________________________________________________
RooAbsData& RooFactoryWSTool::asDATA(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooAbsData object found in workspace

  RooAbsData* data = ws().data(arg) ;
  if (!data) {
      throw string(Form("RooAbsData named %s not found",arg)) ;    
  }
  return *data ;
}



//_____________________________________________________________________________
RooDataHist& RooFactoryWSTool::asDHIST(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooDataHist object found in workspace

  RooAbsData* data = ws().data(arg) ;
  if (!data) {
    throw string(Form("RooAbsData named %s not found",arg)) ;    
  }
  RooDataHist* hist = dynamic_cast<RooDataHist*>(data) ;
  if (!hist) {
    throw string(Form("Dataset named %s is not of type RooDataHist",arg)) ;    
  }
  return *hist ;
}


//_____________________________________________________________________________
RooDataSet& RooFactoryWSTool::asDSET(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as RooDataSet object found in workspace

  RooAbsData* data = ws().data(arg) ;
  if (!data) {
    throw string(Form("RooAbsData named %s not found",arg)) ;    
  }
  RooDataSet* dset = dynamic_cast<RooDataSet*>(data) ;
  if (!dset) {
    throw string(Form("Dataset named %s is not of type RooDataSet",arg)) ;    
  }
  return *dset ;
}



//_____________________________________________________________________________
TObject& RooFactoryWSTool::asOBJ(const char* arg)
{
  TObject* obj = ws().obj(arg) ;
  if (!obj) {
    throw string(Form("Object named %s not found",arg)) ;    
  }
  return *obj ;
}



//_____________________________________________________________________________
const char* RooFactoryWSTool::asSTRING(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as const char* 

  static vector<string> cbuf(10) ;
  static unsigned int cbuf_idx = 0 ;

  // Handle empty string case: return null pointer
  if (arg==0 || strlen(arg)==0) {
    return 0 ;
  }
  
  // Fill cyclical buffer entry with quotation marked stripped version of string literal
  // and return pointer to stripped buffer
  cbuf[cbuf_idx].clear() ;
  const char* p = arg+1 ;
  while(*p && (*p) != '"' && (*p) !='\'' ) {
    cbuf[cbuf_idx] += *(p++) ;
  }  
  const char* ret = cbuf[cbuf_idx].c_str() ;

  // Increment buffer pointer by one
  cbuf_idx++ ;
  if (cbuf_idx==cbuf.size()) cbuf_idx=0 ;

  return ret ;
}


//_____________________________________________________________________________
Int_t RooFactoryWSTool::asINT(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as Int_t

  return atoi(arg) ;
}


//_____________________________________________________________________________
Double_t RooFactoryWSTool::asDOUBLE(const char* arg) 
{
  // CINT constructor interface, return constructor string argument #idx as Double_t

  return atof(arg) ;
}


//_____________________________________________________________________________
void RooFactoryWSTool::registerSpecial(const char* typeName, RooFactoryWSTool::IFace* iface) 
{
  // Register foreign special objects in factory
  hooks()[typeName] = iface ;
}



//_____________________________________________________________________________
std::map<std::string,RooFactoryWSTool::IFace*>& RooFactoryWSTool::hooks() 
{
  if (_hooks) return *_hooks ;
  _hooks = new map<string,IFace*> ;
  return *_hooks ;
}



//_____________________________________________________________________________
std::string RooFactoryWSTool::SpecialsIFace::create(RooFactoryWSTool& ft, const char* typeName, const char* instName, std::vector<std::string> args)
{
  // Concatenate list of args into comma separated string
  char pargs[1024] ;
  pargs[0] = 0 ;
  vector<string>::iterator iter = args.begin() ;
  vector<string> pargv ;
  while(iter!=args.end()) {
    if (strlen(pargs)>0) strlcat(pargs,",",1024) ;
    string tmp = ft.processExpression(iter->c_str()) ;
    strlcat(pargs,tmp.c_str(),1024) ;
    pargv.push_back(tmp) ;
    iter++ ;
  }

  // Handling of special operator pdf class names
  string cl(typeName) ;
  if (cl=="SUM") {

    // SUM::name[a*A,b*B,C]
    ft.add(instName,pargs,kFALSE) ;

  } else if (cl=="RSUM") {

    // RSUM::name[a*A,b*B,C]
    ft.add(instName,pargs,kTRUE) ;

  } else if (cl=="ASUM") {

    // ASUM::name[a*A,b*B,C]
    ft.amplAdd(instName,pargs) ;

  } else if (cl=="PROD") {

    // PROD::name[A,B,C]
    ft.prod(instName,pargs) ;

  } else if (cl=="SIMUL") {

    // PROD::name[cat,state=Pdf,...]
    if (pargv.size()>1) {
      ft.simul(instName,pargv[0].c_str(),strchr(pargs,',')+1) ;
    } else {
      throw string(Form("Need at least two arguments in call to SIMUL::%s, have %d: %s",instName,(Int_t)pargv.size(),pargs)) ;
    }

  } else if (cl=="EXPR") {
    
    // EXPR::name['expr',var,var,...]
    if (args.size()<=2) {
      ft.createArg("RooGenericPdf",instName,pargs) ;
    } else {      
      char genargs[1024] ;
      strlcpy(genargs,args[0].c_str(),1024) ;      
      strlcat(genargs,",{",1024) ;
      for (UInt_t i=1 ; i<args.size() ; i++) {
        if (i!=1) strlcat(genargs,",",1024) ;
        strlcat(genargs,args[i].c_str(),1024) ;
      }
      strlcat(genargs,"}",1024) ;
      ft.createArg("RooGenericPdf",instName,genargs) ;
    }
  
  } else if (cl=="FCONV") {

    // FCONV::name[var,pdf1,pdf2]
    ft.createArg("RooFFTConvPdf",instName,pargs) ;
  
  } else if (cl=="NCONV") {

    // NCONV::name[var,pdf1,pdf2]
    ft.createArg("RooNumConvPdf",instName,pargs) ;
  
  } else if (cl=="sum") {
    
    // sum::name[a,b,c]
    ft.addfunc(instName,pargs) ;

  } else if (cl=="prod") {

    // prod::name[a,b,c]
    ft.prodfunc(instName,pargs) ;
    
  } else if (cl=="expr") {

    // expr::name['expr',var,var,...]
    if (args.size()<=2) {
      ft.createArg("RooFormulaVar",instName,pargs) ;
    } else {      
      char genargs[1024] ;
      strlcpy(genargs,args[0].c_str(),1024) ;      
      strlcat(genargs,",{",1024) ;
      for (UInt_t i=1 ; i<args.size() ; i++) {
        if (i!=1) strlcat(genargs,",",1024) ;
        strlcat(genargs,args[i].c_str(),1024) ;
      }
      strlcat(genargs,"}",1024) ;
      ft.createArg("RooFormulaVar",instName,genargs) ;
    }

  } else if (cl=="nconv") {

    // nconv::name[var,pdf1,pdf2]
    ft.createArg("RooNumConvolution",instName,pargs) ;
  
  } else if (cl=="nll") {

    // nll::name[pdf,data]
    RooNLLVar nll(instName,instName,ft.asPDF(pargv[0].c_str()),ft.asDATA(pargv[1].c_str())) ;
    if (ft.ws().import(nll,Silence())) ft.logError() ;
    
  } else if (cl=="chi2") {

    // chi2::name[pdf,data]
    RooChi2Var nll(instName,instName,ft.asPDF(pargv[0].c_str()),ft.asDHIST(pargv[1].c_str())) ;
    if (ft.ws().import(nll,Silence())) ft.logError() ;
    
  } else if (cl=="profile") {

    // profile::name[func,vars]
    ft.createArg("RooProfileLL",instName,pargs) ;

  } else if (cl=="dataobs") {

    // dataobs::name[dset,func]
    RooAbsArg* funcClone = static_cast<RooAbsArg*>(ft.asARG(pargv[1].c_str()).clone(instName)) ;
    RooAbsArg* arg = ft.asDSET(pargv[0].c_str()).addColumn(*funcClone) ;
    if (!ft.ws().fundArg(arg->GetName())) {
      if (ft.ws().import(*arg,Silence())) ft.logError() ;
    }
    delete funcClone ;

  } else if (cl=="int") {    
    
    // int::name[func,intobs]
    // int::name[func,intobs|range]
    // int::name[func,intobs,normobs]
    // int::name[func,intobs|range,normobs]

    if (pargv.size()<2 || pargv.size()>3) {
      throw string(Form("int::%s, requires 2 or 3 arguments, have %d arguments",instName,(Int_t)pargv.size())) ;
    }

    RooAbsReal& func = ft.asFUNC(pargv[0].c_str()) ;

    char buf[256] ;
    strlcpy(buf,pargv[1].c_str(),256) ;
    char* save ;
    const char* intobs = strtok_r(buf,"|",&save) ;
    if (!intobs) intobs="" ;

    const char* range = strtok_r(0,"",&save) ;
    if (!range) range="" ;

    RooAbsReal* integral = 0 ;
    if (pargv.size()==2) {
      if (range && strlen(range)) {
        integral = func.createIntegral(ft.asSET(intobs),Range(range)) ;
      } else {
        integral = func.createIntegral(ft.asSET(intobs)) ;
      }
    } else {
      if (range && strlen(range)) {
        integral = func.createIntegral(ft.asSET(intobs),Range(range),NormSet(ft.asSET(pargv[2].c_str()))) ;
      } else {
        integral = func.createIntegral(ft.asSET(intobs),NormSet(ft.asSET(pargv[2].c_str()))) ;
      }
    }

    integral->SetName(instName) ;
    if (ft.ws().import(*integral,Silence())) ft.logError() ;
    
  } else if (cl=="deriv") {
    
    // derive::name[func,obs,order]

    if (pargv.size()<2 || pargv.size()>3) {
      throw string(Form("deriv::%s, requires 2 or 3 arguments, have %d arguments",instName,(Int_t)pargv.size())) ;
    }

    RooAbsReal& func = ft.asFUNC(pargv[0].c_str()) ;

    RooAbsReal* derivative(0) ;
    if (pargv.size()==2) {
      derivative = func.derivative(ft.asVAR(pargv[1].c_str()),1) ;
    } else {
      derivative = func.derivative(ft.asVAR(pargv[1].c_str()),ft.asINT(pargv[2].c_str())) ;
    }

    derivative->SetName(instName) ;
    if (ft.ws().import(*derivative,Silence())) ft.logError() ;

  } else if (cl=="cdf") {

    // cdf::name[pdf,obs,extranormobs]

    if (pargv.size()<2 || pargv.size()>3) {
      throw string(Form("cdf::%s, requires 2 or 3 arguments, have %d arguments",instName,(Int_t)pargv.size())) ;
    }

    RooAbsPdf& pdf = ft.asPDF(pargv[0].c_str()) ;

    RooAbsReal* cdf(0) ;
    if (pargv.size()==2) {
      cdf = pdf.createCdf(ft.asSET(pargv[1].c_str())) ;
    } else {
      cdf = pdf.createCdf(ft.asSET(pargv[1].c_str()),ft.asSET(pargv[2].c_str())) ;
    }

    cdf->SetName(instName) ;
    if (ft.ws().import(*cdf,Silence())) ft.logError() ;


  } else if (cl=="PROJ") {

    // PROJ::name(pdf,intobs)
    if (pargv.size()!=2) {
      throw string(Form("PROJ::%s, requires 2 arguments, have %d arguments",instName,(Int_t)pargv.size())) ;
    }

    RooAbsPdf& pdf = ft.asPDF(pargv[0].c_str()) ;
    RooAbsPdf* projection = pdf.createProjection(ft.asSET(pargv[1].c_str())) ;
    projection->SetName(instName) ;

    if (ft.ws().import(*projection,Silence())) ft.logError() ;

  } else if (cl=="set") {

    // set::name(arg,arg,...)
    if (ft.ws().defineSet(instName,pargs)) {
      ft.logError() ;
      return string(instName) ;
    }
    
  } else {

    throw string(Form("RooFactoryWSTool::SpecialsIFace::create() ERROR: Unknown meta-type %s",typeName)) ;

  }
  return string(instName) ;    
}

---