MaranaEventLoop.C

Go to the documentation of this file.

// Class MaranaEventLoop
//
// Author     : E. Sauvan
// Created    : 27/10/2002
// Last update: $Date: 2005/07/18 12:02:03 $
//          by: $Author: blist $
//         
//        This is an exemple class for histograming and implementing
//              an analysis channel


// system 
#include <cstdio>
#include <cctype>  // for isblank etc
#include <vector>
// root 
#include <TVirtualPS.h>
#include <TPad.h>
// Marana 
#include "Marana/TMarEvent.h"
#include "Marana/TMarBody.h"
#include "SFHMarana/MaranaEventLoop.h"
#include "SFHMarana/MarSysBinningFun.h"
// jbltools 
#include "jbltools/sfh/SFH1F.h"
#include "jbltools/sfh/SFSetOfHistograms.h"
#include "jbltools/sfh/SFMatrixOfHistograms.h"
#include "jbltools/sfh/SFHProf.h"
#include "jbltools/sfh/SFSetOfHProf.h"
#include "jbltools/sfh/SFH2F.h"
#include "jbltools/sfh/SFSetOfH2F.h"
#include "jbltools/sfh/SFROList.h"
#include "jbltools/sfh/SFCutROList.h"
#include "jbltools/sfh/FloatFun.h"
#include "jbltools/sfh/BinningFun.h"
#include "jbltools/sfh/FloatFunBinning.h"
#include "jbltools/sfh/FillIterator.h"
#include "jbltools/sfh/BaseCut.h"
#include "jbltools/sfh/IntFun.h"
#include "jbltools/sfh/FloatFun.h"
#include "jbltools/sfh/Binning.h"
#include "jbltools/sfh/NamedBinning.h"
#include "jbltools/sfh/JBLParserResult.h"

ClassImp(MaranaEventLoop)

MaranaEventLoop::MaranaEventLoop(const TRunHisto& run,string Tag,string textfile)
{
  // needed for booking self-filling histograms
  sysfun = new MarSysBinningFun (run.Syst, run.output_file_root);
  syscut = &(*sysfun == 0);
  syscut->setName("syscut");
  rolist = 0;
  fillrolist = new SFROList();
  
  // copy of Marana/THisto constructor
  debug=run.debug;
  
  NumSelEve=0;
  NumEve=0;
  NumSelEvePerRun=0;
  NumSelEvePerRunOld=0;
  fp=NULL;
  hatCopy=NULL; 
  modsCopy=NULL; 

//... Open file to write selected events
  if(textfile!=""){
    
    if ((fp  = fopen(textfile.c_str(),"w+")) == NULL) {
       cout << "THisto::THisto : error in opening "<< textfile <<endl;
    }
  }
//...read the parameters and book the histograms in the / dir
  ReadParameters(run.histo_par_file,Tag,0);

//...read the parameters and book the histograms in each sub-directory of run.output_file_root
  TList* dir=run.output_file_root->GetListOfKeys();
  for(int i=0;i<dir->GetSize();i++){
    cout << "--> Create histo in  dir "<< dir->At(i)->GetName() << endl;
    run.output_file_root->cd(dir->At(i)->GetName());
    
    ReadParameters(run.histo_par_file,Tag,i+1);
  }
  run.output_file_root->cd();

//... open hat/mods file to write selected events
  
  if(run.out_file_mods_hat!=""){
    

    for(list<string>::const_iterator itag=run.TagToBeWritten.begin();itag!=run.TagToBeWritten.end();itag++)
    {
        

        if((*itag)==Tag){ 
          string dir="";
          string file="";
          unsigned int pos=run.out_file_mods_hat.rfind("/");
          if(pos!=string::npos){
             dir=run.out_file_mods_hat.substr(0,pos+1);   
             file=run.out_file_mods_hat.substr(pos+1) ;          
             
          }  else {
             file=run.out_file_mods_hat; 
          
          }   

            
          if ( gH1Tree->GetChain(3)) {
            string hatname;

            hatname=dir+"HAT."+Tag+"."+file;
            cout << "\n--> output hat: "<< hatname<<endl;      
        
            hatCopy = gH1Tree->CopyTree("HAT", hatname.c_str() );
            if (!hatCopy) throw string("THisto::THisto= could not copy HAT-tree!");
            
        
          }

          if ( gH1Tree->GetChain(2) ) {
            string modsname;
            modsname=dir+"MODS."+Tag+"."+file;
            cout << "\n--> output mods: "<< modsname<<endl;    

            modsCopy = gH1Tree->CopyTree("MODS", modsname.c_str());
            if (!modsCopy)throw string("THisto::THisto= could not copy MODS-tree!");
          }
        }

     }  
  }
  
}

MaranaEventLoop::~MaranaEventLoop()
{}


void MaranaEventLoop::Fill()
{
//--- Fill self-filling stuff
  fillrolist->Fill();
  return;
} 

void MaranaEventLoop::Draw(TPad *pad, TVirtualPS *ps,
                      Option_t *option,
                      bool across,
                      Int_t nx, Int_t ny,
                      bool log, bool recursive) {
  if (fillrolist)
    fillrolist->Draw(pad, ps, option, across, nx, ny, log, recursive);
}

// add stuff to book SFHs by steering. Overwriting/ copying from Marana/MaranaEventLoop!

bool MaranaEventLoop::ParseHistos(istringstream& inParams, string&type, Int_t hoffset,
                                string Tag, list<string>& PartTags,const Short_t SysNum)
{   
// Create a new SFCutROList for all histos of this type, and register it
    rolist = new SFCutROList (fillrolist);

//     Read and parse histo definitions
    bool   stop=false;
    string item;
    string end="end_";
    char line[1024];
    while (!stop && inParams.getline(line, 1024)) {
        istringstream linestream(line);
        item = "";
        linestream >> item;
//...read and create h1f histos ;
        if(item=="h1f") {
                ReadParamsh1f(linestream,type,hoffset,SysNum,kFALSE);
        }
        else if(item=="h1f_sys"){
                ReadParamsh1f(linestream,type,hoffset,SysNum,kTRUE);
        }
        else if(item=="h2f"){
                ReadParamsh2f(linestream,type,hoffset,SysNum,kFALSE);
        }
        else if(item=="h2f_sys"){
                ReadParamsh2f(linestream,type,hoffset,SysNum,kTRUE);
        }
//...read and create h1f generic histos ;
        else if(item=="h1fgeneric") {
                ReadParamsh1fgeneric(linestream,type,hoffset,PartTags,SysNum,kFALSE);
        }
        else if(item=="h1fgeneric_sys") {
                ReadParamsh1fgeneric(linestream,type,hoffset,PartTags,SysNum,kTRUE);
        }

//... for bin histograms with fix bin width
        else if(item=="2dbin"){
                ReadParams2dbin(linestream,type,hoffset,SysNum);
        }

//... for bin histograms with fix bin width
        else if(item=="1binfh1f"){
                ReadParams1binfh1f(linestream,type,hoffset,SysNum,kFALSE);
        }
        else if(item=="1binfh1f_sys"){
                ReadParams1binfh1f(linestream,type,hoffset,SysNum,kTRUE);
        }
//... for bin histogram with variable bin width
        else if(item=="1binvbh1f"){
                ReadParams1binVarSizeh1f(linestream,type,hoffset,SysNum,kFALSE);
        }
        else if(item=="1binvbh1f_sys"){
                ReadParams1binVarSizeh1f(linestream,type,hoffset,SysNum,kTRUE);
        }
//... for bin histograms with variable bin width
        else if(item=="1binvh1f"){
                ReadParams1binvh1f(linestream,type,hoffset,SysNum,kFALSE);
        }
        else if(item=="1binvh1f_sys"){
                ReadParams1binvh1f(linestream,type,hoffset,SysNum,kTRUE);
        }
//... for 2bin histograms with fix bin width
        else if(item=="2binfh1f"){
                ReadParams2binfh1f(linestream,type,hoffset,SysNum,kFALSE);
        }
        else if(item=="2binfh1f_sys"){
                ReadParams2binfh1f(linestream,type,hoffset,SysNum,kTRUE);
        }
//... for 2bin histograms with variable bin width
        else if(item=="2binvh1f"){
                ReadParams2binvh1f(linestream,type,hoffset,SysNum,kFALSE);
        }
        else if(item=="2binvh1f_sys"){
                ReadParams2binvh1f(linestream,type,hoffset,SysNum,kTRUE);
        } 
        else if (item=="sfh1f"){
                // cout << "booking sfh" << endl;
                ReadParamssfh1f(linestream,type,hoffset,SysNum,kFALSE);
        } 
        else if (item=="sfh1f_sys"){
                // cout << "booking sfh" << endl;
                ReadParamssfh1f(linestream,type,hoffset,SysNum,kTRUE);
        } 
        else if (item=="sfhprof"){
                // cout << "booking sfhprof" << endl;
                ReadParamssfhprof(linestream,type,hoffset,SysNum,kFALSE);
        } 
        else if (item=="sfhprof_sys"){
                // cout << "booking sfhprof" << endl;
                ReadParamssfhprof(linestream,type,hoffset,SysNum,kTRUE);
        } 
        else if (item=="sfh2f"){
                // cout << "booking sfh2f" << endl;
                ReadParamssfh2f(linestream,type,hoffset,SysNum,kFALSE);
        } 
        else if (item=="sfh2f_sys"){
                // cout << "booking sfh2f" << endl;
                ReadParamssfh2f(linestream,type,hoffset,SysNum,kTRUE);
        } 
        else if (item=="sfseth1f"){
                // cout << "booking sfh set" << endl;
                ReadParamssfseth1f(linestream,type,hoffset,SysNum,kFALSE);
        } 
        else if (item=="sfseth1f_sys"){
                // cout << "booking sfh set" << endl;
                ReadParamssfseth1f(linestream,type,hoffset,SysNum,kTRUE);
        } 
        else if (item=="def"){
                // cout << "creating function object" << endl;
                ReadParamsdef(linestream,SysNum);
        } 
        else if (item=="bin"){
                // cout << "creating binning object" << endl;
                ReadParamsbin(linestream,SysNum);
        }
        else if (item=="typecut"){
                cout << "assigning cut to type" << endl;
                ReadParamstypecut(linestream,type);
        } 
        else if (item == end + Tag) {
            stop = true;
        }                       //if
          
    }                   //while 
  return stop;
}



void MaranaEventLoop::ReadParamssfh1f(istream &inParams, const string& type, Int_t hoffset,  Short_t SysNum, Bool_t BookforSys)
{
//     read params for bin histo h1f

   if (SysNum != 0) return;
 
   int num=0;
   string name, comment;
   string xfun, cut, wfun;

   inParams >> num >> name;
   cout << "MaranaEventLoop::ReadParamssfh1f: reading parameters\n";
   cout << "   num=" << num << "\n";
   cout << "   name='" << name << "'\n";
   num += hoffset;
   
   
   comment = ReadString (inParams);
   cout << "   comment='" << comment << "'\n";
   
   Binning *binning = ReadBinningInfo (inParams);
   if (!binning) {
     cerr << "MaranaEventLoop::ReadParamssfh1f: Cannot read binning!" << endl;
     return;
   }
   
   xfun = ReadString (inParams);
   cut  = ReadString (inParams);
   wfun = ReadString (inParams);
   
   cout << "   xfun='" << xfun  << "'\n";
   cout << "   cut='" << cut << "'\n";
   cout << "   wfun='" << wfun << "'\n";
  
     
   name="h"+type+name;
   comment=comment+" ("+type+")";
   if(SysNum==0) cout << num <<" " <<name <<" : "<< comment << endl;
   
   FloatFun     *xf;
   BaseCut      *c;
   FloatFun     *wf;
   FillIterator *it;
           
   if(SysNum==0) {
     if (BookforSys) {
       cout << "ParseHistos: Creating SFSetOfHistograms" << endl;
       SFSetOfHistograms *sh = new SFSetOfHistograms (name.c_str(), 
                               comment.c_str(), *binning, 
                               rolist, xfun, cut, wfun, *sysfun, 0); 
       xf = sh->getXFun();       
       c  = sh->getCut();         
       wf = sh->getWFun();       
       it = sh->getIter();   
     }
     else {
       BaseCut *bc = BaseCutPoR(cut).pbc;
       BaseCut *thecut;
       if (bc) thecut = &(*syscut && *bc);
       else thecut = syscut;
       cout << "ParseHistos: Creating SFH1F" << endl;
       SFH1F *h = new SFH1F (name.c_str(), 
                             comment.c_str(), *binning, 
                              rolist, xfun, thecut, 0); 
       xf = h->getXFun();       
       c  = h->getCut();         
       wf = h->getWFun();       
       it = h->getIter();   
     
     }                   
     if (xf) cout << "  ParseHistos: XFun = " << xf->getName() << endl;                    
     if (c)  cout << "  ParseHistos: Cut = "  << c->getName()  << endl;                    
     if (wf) cout << "  ParseHistos: WFun = " << wf->getName() << endl;                    
     if (it) cout << "  ParseHistos: Iter = " << it->getName() << endl; 
     delete binning;                       
   }      

   return;
}

void MaranaEventLoop::ReadParamssfhprof(istream &inParams, const string& type, Int_t hoffset, Short_t SysNum, Bool_t BookforSys)
{
//     read params for bin histo h1f

   if (SysNum != 0) return;
 
   int num=0;
   string name, comment;
   string xfun, yfun, cut, wfun;
   float ylow, yup;
   string option;

   inParams >> num >> name;
   cout << "MaranaEventLoop::ReadParamssfhprof: reading parameters\n";
   cout << "   num=" << num << "\n";
   cout << "   name='" << name << "'\n";
   num += hoffset;
   
   comment = ReadString (inParams);
   cout << "   comment='" << comment << "'\n";
   
   Binning *binning = ReadBinningInfo (inParams);
   if (!binning) {
     cerr << "MaranaEventLoop::ReadParamssfhprof: Cannot read binning!" << endl;
     return;
   }
   
   inParams >> ylow >> yup; 
   option = ReadString (inParams);
   
   xfun = ReadString (inParams);
   yfun = ReadString (inParams);
   cut  = ReadString (inParams);
   wfun = ReadString (inParams);
   
   cout << "   ylow=   " << ylow   << "'\n";
   cout << "   yup=    " << yup    << "'\n";
   cout << "   option='" << option << "'\n";
   cout << "   xfun=  '" << xfun   << "'\n";
   cout << "   yfun=  '" << yfun   << "'\n";
   cout << "   cut=   '" << cut    << "'\n";
   cout << "   wfun=  '" << wfun   << "'\n";
  
     
   name="h"+type+name;
   comment=comment+" ("+type+")";
   if(SysNum==0) cout << num <<" " <<name <<" : "<< comment << endl;
   
   FloatFun     *xf;
   FloatFun     *yf;
   BaseCut      *c;
   FloatFun     *wf;
   FillIterator *it;
           
   if(SysNum==0) {
     if (BookforSys) {
       cout << "ParseHistos: Creating SFSetOfHProf" << endl;
       SFSetOfHProf *sh = new SFSetOfHProf (name.c_str(), 
                               comment.c_str(), *binning, ylow, yup, option.c_str(),
                               rolist, xfun, yfun, cut, wfun, *sysfun, 0); 
       xf = sh->getXFun();       
       yf = sh->getYFun();       
       c  = sh->getCut();         
       wf = sh->getWFun();       
       it = sh->getIter();   
     }
     else {
       BaseCut *bc = BaseCutPoR(cut).pbc;
       BaseCut *thecut;
       if (bc) thecut = &(*syscut && *bc);
       else thecut = syscut;
       cout << "ParseHistos: Creating SFHProf" << endl;
       SFHProf *h = new SFHProf (name.c_str(), 
                             comment.c_str(), *binning, ylow, yup, option.c_str(),
                              rolist, xfun, yfun, thecut, wfun); 
       xf = h->getXFun();       
       yf = h->getYFun();       
       c  = h->getCut();         
       wf = h->getWFun();       
       it = h->getIter();   
     
     }                   
     if (xf) cout << "  ParseHistos: XFun = " << xf->getName() << endl;                    
     if (yf) cout << "  ParseHistos: YFun = " << yf->getName() << endl;                    
     if (c)  cout << "  ParseHistos: Cut = "  << c->getName()  << endl;                    
     if (wf) cout << "  ParseHistos: WFun = " << wf->getName() << endl;                    
     if (it) cout << "  ParseHistos: Iter = " << it->getName() << endl; 
     delete binning;                       
   }      

   return;
}

void MaranaEventLoop::ReadParamssfh2f(istream &inParams, const string& type, Int_t hoffset, Short_t SysNum, Bool_t BookforSys)
{
//     read params for bin histo h1f

   if (SysNum != 0) return;
 
   int num=0;
   string name, comment;
   string xfun, yfun, cut, wfun;

   inParams >> num >> name;
   cout << "MaranaEventLoop::ReadParamssfh2f: reading parameters\n";
   cout << "   num=" << num << "\n";
   cout << "   name='" << name << "'\n";
   num += hoffset;
   
   comment = ReadString (inParams);
   cout << "   comment='" << comment << "'\n";
   
   Binning *xbinning = ReadBinningInfo (inParams);
   if (!xbinning) {
     cerr << "MaranaEventLoop::ReadParamssfh2f: Cannot read x binning!" << endl;
     return;
   }
   
   Binning *ybinning = ReadBinningInfo (inParams);
   if (!ybinning) {
     cerr << "MaranaEventLoop::ReadParamssfh2f: Cannot read y binning!" << endl;
     return;
   }
   
   xfun = ReadString (inParams);
   yfun = ReadString (inParams);
   cut  = ReadString (inParams);
   wfun = ReadString (inParams);
   
   cout << "   xfun=  '" << xfun   << "'\n";
   cout << "   yfun=  '" << xfun   << "'\n";
   cout << "   cut=   '" << cut    << "'\n";
   cout << "   wfun=  '" << wfun   << "'\n";
  
     
   name="h"+type+name;
   comment=comment+" ("+type+")";
   if(SysNum==0) cout << num <<" " <<name <<" : "<< comment << endl;
   
   FloatFun     *xf;
   FloatFun     *yf;
   BaseCut      *c;
   FloatFun     *wf;
   FillIterator *it;
           
   if(SysNum==0) {
     if (BookforSys) {
       cout << "ParseHistos: Creating SFSetOfH2F" << endl;
       SFSetOfH2F *sh = new SFSetOfH2F (name.c_str(), 
                               comment.c_str(), *xbinning, *ybinning, 
                               rolist, xfun, yfun, cut, wfun, *sysfun); 
       xf = sh->getXFun();       
       yf = sh->getYFun();       
       c  = sh->getCut();         
       wf = sh->getWFun();       
       it = sh->getIter();   
     }
     else {
       BaseCut *bc = BaseCutPoR(cut).pbc;
       BaseCut *thecut;
       if (bc) thecut = &(*syscut && *bc);
       else thecut = syscut;
       cout << "ParseHistos: Creating SFH2F" << endl;
       SFH2F *h = new SFH2F (name.c_str(), 
                             comment.c_str(), *xbinning, *ybinning, 
                              rolist, xfun, yfun, thecut, wfun); 
       xf = h->getXFun();       
       yf = h->getYFun();       
       c  = h->getCut();         
       wf = h->getWFun();       
       it = h->getIter();   
     
     }                   
     if (xf) cout << "  ParseHistos: XFun = " << xf->getName() << endl;                    
     if (yf) cout << "  ParseHistos: YFun = " << yf->getName() << endl;                    
     if (c)  cout << "  ParseHistos: Cut = "  << c->getName()  << endl;                    
     if (wf) cout << "  ParseHistos: WFun = " << wf->getName() << endl;                    
     if (it) cout << "  ParseHistos: Iter = " << it->getName() << endl; 
     delete xbinning;                       
     delete ybinning;                       
   }      

   return;
}


void MaranaEventLoop::ReadParamsdef(istream &inParams, Short_t SysNum)
{
//     read params for bin histo h1f

   if (SysNum != 0) return;
 
   string name;

   inParams >> name;
   // cout << "MaranaEventLoop::ReadParamsdef: reading parameters\n";
   // cout << "   name='" << name << "'\n";
   
   JBLParserResult definition;
   inParams >> definition;
   // cout << "   definition='" << definition << "'\n";
   
   if (definition.isIntFun()) {
     IntFun& fun = *definition.getIntFunValue();
     fun.setName(name);
     cout << "Definition of new FloatFun '" << fun.getName() <<  "', definition " << fun.getDefinition() << endl;
   } else if (definition.isFloatFun()) {
     FloatFun& fun = *definition.getFloatFunValue();
     assert (&fun);
     fun.setName(name);
     cout << "Definition of new FloatFun '" << fun.getName() <<  "', definition " << fun.getDefinition() << endl;
   } else if (definition.isBaseCut()) {
     BaseCut& cut = *definition.getBaseCutValue();
     assert (&cut);
     cut.setName(name);
     cout << "Definition of new BaseCut '" << cut.getName() <<  "', definition " << cut.getDefinition() << endl;
   } else if (definition.isInt()) {
     IntFun& fun = *new ConstIntFun (definition.getIntValue(), name);
     cout << "Definition of  new Integer '" << fun.getName() <<  "', definition " << fun.getDefinition() << endl;
   } else if (definition.isFloat()) {
     FloatFun& fun = *new ConstFun (definition.getFloatValue(), name);
     cout << "Definition of new Float '" << fun.getName() <<  "', definition " << fun.getDefinition() << endl;
   } else if (definition.isBool()) {
     BaseCut& fun = *new ConstBaseCut (definition.getBoolValue(), name);
     cout << "Definition of new Bool '" << fun.getName() <<  "', definition " << fun.getDefinition() << endl;
   }
   else cerr << "MaranaEventLoop::ReadParamsdef: Cannot define '" << name << "'\n";
   return;
}

void MaranaEventLoop::ReadParamssfseth1f(istream &inParams, const string& type, Int_t hoffset,  Short_t SysNum, Bool_t BookforSys)
{
//     read params for bin histo h1f

   if (SysNum != 0) return;
 
   int num=0;
   string name, comment;
   string xfun, cut, wfun;

   inParams >> num >> name;
   cout << "MaranaEventLoop::ReadParamssfseth1f: reading parameters\n";
   cout << "   num=" << num << "\n";
   cout << "   name='" << name << "'\n";
   num += hoffset;
   
   comment = ReadString (inParams);
   cout << "   comment='" << comment << "'\n";
   
   Binning *axisbinning = ReadBinningInfo (inParams);
   if (!axisbinning) {
     cerr << "MaranaEventLoop::ReadParamssfseth1f: Cannot read binning!" << endl;
     return;
   }
   
   xfun = ReadString (inParams);
   cut  = ReadString (inParams);
   wfun = ReadString (inParams);
   
   cout << "   xfun='"  << xfun   << "'\n";
   cout << "   cut="    << cut    << "\n";
   cout << "   wfun=" << wfun << "\n";
  
   BinningFun *binningfun = ReadBinningFunInfo (inParams);
   if (!binningfun) {
     cerr << "MaranaEventLoop::ReadParamssfseth1f: Cannot read binning function!" << endl;
     return;
   }
     
   name="h"+type+name;
   comment=comment+" ("+type+")";
   if(SysNum==0) cout << num <<" " <<name <<" : "<< comment << endl;
   FloatFun     *xf;
   BaseCut      *c;
   FloatFun     *wf;
   FillIterator *it;
           
   if(SysNum==0) {
     if (BookforSys) {
       cout << "ParseHistos: Creating SFMatrixOfHistograms" << endl;
       SFMatrixOfHistograms *mh = new SFMatrixOfHistograms (name.c_str(), 
                               comment.c_str(), *axisbinning, 
                               rolist, xfun, cut, wfun, *sysfun, *binningfun); 
       xf = mh->getXFun();       
       c  = mh->getCut();         
       wf = mh->getWFun();       
       it = mh->getIter();   
     }
     else {
       BaseCut *bc = BaseCutPoR(cut).pbc;
       BaseCut *thecut;
       if (bc) thecut = &(*syscut && *bc);
       else thecut = syscut;
       cout << "ParseHistos: Creating SFSetOfHistograms" << endl;
       SFSetOfHistograms *sh = new SFSetOfHistograms (name.c_str(), 
                             comment.c_str(), *axisbinning, 
                              rolist, xfun, thecut, wfun, *binningfun); 
     
       xf = sh->getXFun();       
       c  = sh->getCut();         
       wf = sh->getWFun();       
       it = sh->getIter();   
     }                   
     if (xf) cout << "  ParseHistos: XFun = " << xf->getName() << endl;                    
     if (c)  cout << "  ParseHistos: Cut = "  << c->getName()  << endl;                    
     if (wf) cout << "  ParseHistos: WFun = " << wf->getName() << endl;                    
     if (it) cout << "  ParseHistos: Iter = " << it->getName() << endl; 
     delete axisbinning;                       
   }      

   return;
}


void MaranaEventLoop::ReadParamsbin(istream &inParams, Short_t SysNum)
{
//     read params for bin histo h1f

   if (SysNum != 0) return;
 
   string name;

   inParams >> name;
   cout << "MaranaEventLoop::ReadParamsbin: reading parameters\n";
   cout << "   name='" << name << "'\n";
   
   Binning *binning = ReadBinningInfo (inParams);
   if (!binning) {
     cerr << "MaranaEventLoop::ReadParamsbin: Cannot create binning!";
   }
   
   NamedBinning& result = *new NamedBinning(*binning, name);
   delete binning;
   std::string definition = "";
   if (result.isEquidistant()) {
     definition = NamedO::str('(')+NamedO::str(result.getNBins())+", "
                  +NamedO::str(result.getLowerEdge())+", "+NamedO::str(result.getUpperEdge())+')';
   }
   else {
     definition = NamedO::str('{');
     for (int i = 0; i < result.getNBins(); i++) 
       definition += NamedO::str(result.getLowerBinEdge(i))+", ";
     definition += NamedO::str(result.getUpperEdge())+'}';
   }
   result.setDefinition(definition);
   cout << "created binning " << result << endl;
   return;
}

void MaranaEventLoop::ReadParamstypecut(istream &inParams, const string& type) {
  cout << "ReadParamstypecut, type='" << type << "'" << endl;
  bool doesapply = false;
  string thistype;
  string cut;
  do {
    thistype = ReadString (inParams, ',');
    cout << "  ...type = '" << thistype << "', peek=" << inParams.peek() << endl;
    if (thistype == type) doesapply = true;
    if (inParams.peek()!=',') break;
    inParams.get(); // swallow comma
  } while (inParams);
  cut = ReadString (inParams);
  if (doesapply && rolist) {
    cout << "  cut '" << cut << "' applies to type '" << type << "'" << endl;
    BaseCutPoR newcut(cut);
    BaseCut *oldcut = rolist->getCut();
    if (newcut.pbc) {
      cout << "    Additional cut is '" << newcut.pbc->getName()  << "'" << endl;
      if (oldcut) {
        cout << "    Old cut is '" << oldcut->getName()  << "'" << endl;
        rolist->setCut(*oldcut && *newcut.pbc);
      }
      else {
        rolist->setCut(newcut.pbc);
      }
    }
    if (rolist->getCut())
      cout << "    Cut is now '" << rolist->getCut()->getName() << "'." << endl;
  }
}

Binning *MaranaEventLoop::ReadBinningInfo (istream &is) {
  Binning *result = 0;
  cout << "MaranaEventLoop::ReadBinningInfo starting..." << endl;
  std::string word;
  
  
  // eat leading whitespace
  while (isblank(is.peek())) is.get();
  
  if (is.peek() == '{') {
    cout << "parsing non-equidistant binning" << endl;
    // { a, b, c...} defines a non-equidistant binning
    is.get();
    double d = 0;
    std::vector<double> boundaries;
    while (is && is.peek() != '}') {
      is >> d;
      cout << "   boundary " << d << endl;
      boundaries.push_back(d);
      while (isblank(is.peek())) is.get();
      if (is.peek() == ',') is.get();
      while (isblank(is.peek())) is.get();
    }
    if (!is) {
      cerr << "Error in MaranaEventLoop::ReadBinningInfo: " 
           << "encountered input end before '}' was found!" << endl;
    }
    else {
      // eat last '}'
      is.get();
    }
    if (boundaries.size() > 1) {
      result = new Binning (boundaries.size()-1, &(boundaries[0]));
    }
    else {
      cerr << "Error in MaranaEventLoop::ReadBinningInfo: " 
           << "Less than 2 bin boundaries found!" << endl;
    }
  }
  else if (isdigit(is.peek()) || is.peek() == '.') {
    cout << "parsing equidistant binning" << endl;
    int nbins;
    double xlow, xup;
    is >> nbins >> xlow >> xup;
    cout << "  nbins = " << nbins 
         << ", xlow = " << xlow << ", xup " << xup << endl;
    result = new Binning (nbins, xlow, xup);
  }
  else  {
    cout << "parsing named binning" << endl;
    std::string identifier;
    is >> identifier;
    NamedBinning *def = NamedBinning::getObject(identifier);
    if (def) {
      result = new Binning (*def);
      cout << "  found " << *def << endl;
    }
    else {
      cerr << "Error in MaranaEventLoop::ReadBinningInfo: " 
           << "Cannot find named binning '" << identifier << "'!" << endl;
    }
  }
  return result;
}

BinningFun *MaranaEventLoop::ReadBinningFunInfo (istream &is) {
  BinningFun *result = 0;
  FloatFun *fun = 0;
  cout << "MaranaEventLoop::ReadBinningFunInfo starting..." << endl;
  
  // eat leading whitespace
  while (isblank(is.peek())) is.get();
  
  if (is.peek() == '[') {
    is.get();
    cout << "   parsing FloatFun definition\n";
    
    std::string word;
    is >> word;
  
    JBLParserResult definition = JBLParserResult::parseExpression (word);
    cout << "   definition='" << definition << "'\n";
  
    if (definition.isFloatFun()) {
       fun = definition.getFloatFunValue();
       cout << "  Found FloatFun '" << fun->getName() <<  "', definition " << fun->getDefinition() << endl;
     }
     else {
       cerr << "MaranaEventLoop::ReadBinningFunInfo: Expected a '" << definition.getName() 
            << "' to be a FloatFun definition!\n";
       return 0;
     }
     
     Binning *binning = ReadBinningInfo(is);
     
     if (!binning) {
       cerr << "MaranaEventLoop::ReadBinningFunInfo: Expected a binning here!\n";
       return 0;
     }
     
     result = new FloatFunBinning (*fun, *binning, fun->getName());
     
    // eat trailing whitespace
    while (isblank(is.peek())) is.get();
  
    if (is.peek() == ']') {
      is.get();
    }
    else {
      cerr << "MaranaEventLoop::ReadBinningFunInfo: Expected a ']' instead of " << is.peek() << "!\n";
      return result;
    }
  }
  else {
    cerr << "MaranaEventLoop::ReadBinningFunInfo: Cannot handle named FloatFunBinnings yet!\n";
    return 0;
  }
  return result;
}


std::string MaranaEventLoop::ReadString (istream &is, char delim) {
  std::string result;

  // eat leading whitespace
  while (is && isblank(is.peek())) is.get();
  
  if (is && is.peek() == '\047') {
    // single quotation marks
    is.get();
    while (is && isblank(is.peek())) is.get();
    while (is && is.peek() != '\047' && is.peek()!='\377') result += is.get();
    if (is && is.peek()!='\377') is.get();
  }
  else if (is && is.peek() == '"') {
    // double quotation marks
    is.get();
    while (is && isblank(is.peek())) is.get();
    while (is && is.peek() != '"' && is.peek()!='\377') result += is.get();
    if (is && is.peek()!='\377') is.get();
  }
  else {
    // no quotation marks
    while (is && !isblank(is.peek()) && is.peek()!=delim && is.peek()!='\377') result += is.get();
  }
  // eat trailing whitespace
  while (is && isblank(is.peek())) is.get();
  
  // remove trailing whitespace from result
  string::size_type pos = result.find_last_not_of (' ');
  if (pos != string::npos) result.erase (pos+1);
  return result;
}

---