ParameterLimitsReader.cc

//!
//!  \brief Reader for Parameter Limits
//!
//!  This class add user defined parameter limits
//!
//!  \author Hartmut Stadie
//!  \date  2008/12/12
//!  $Id: ParameterLimitsReader.cc,v 1.16 2010/11/01 15:47:41 stadie Exp $
//!   
#include "ParameterLimitsReader.h"

#include "CalibData.h"
#include "ConfigFile.h"
#include "Parameters.h"
#include "ParameterLimit.h"

#include <iostream>

ParameterLimitsReader::ParameterLimitsReader(const std::string& configfile, Parameters* p) :
  EventReader(configfile,p)
{
  std::vector<double> limits = bag_of<double>(config_->read<std::string>( "Jet Parameter Limits",""));
  
  // In case limits are explicitly set via config file
  if( limits.size() > 0 && limits.size() % 4 == 0 ) {
    std::cout << "Using user defined parameter limits:" << std::endl;
    for(unsigned int i = 0 ; i < limits.size() ; i += 4) {
      int index = (int)limits[i];
      for(int j = par_->numberOfTowerParameters() + index; 
	  j <  par_->numberOfParameters() ; 
	  j += par_->numberOfJetParametersPerBin()) {
	//tower parameter constraints (par_id min max error)
	par_limits.push_back(ParLimit(j,limits[i+1],limits[i+2],
					    limits[i+3]));
      }
    }
  }
  // In case default limits are to be used
  else if( limits.size() == 1 ) {
    std::string parclass = config_->read<std::string>("Parametrization Class","");
    std::cout << "Using default parameter limits for '" << parclass << "'" << std::endl;

    // For Gauss Function
    // For Gauss Function in one bin
    if( parclass == "ResolutionGaussAvePt" ) {
      for(int i = 0; i < par_->numberOfParameters(); i++) {
	double min = 2./sqrt(M_PI);   // Log-term has to be positive
	double max = 10000.;
	par_limits.push_back(ParLimit(i,min,max,limits.at(0)));
      } // End of loop over parameters
    } else if( parclass == "ResolutionGauss" ) {
      for(int i = 0; i < par_->numberOfParameters(); i++) {
	double min = 2.*sqrt(2./M_PI);   // Log-term has to be positive
	double max = 10000.;
	par_limits.push_back(ParLimit(i,min,max,limits.at(0)));
      } // End of loop over parameters
    }
  }

  // Catch wrong syntax in config file
  else if(limits.size() > 1) {
    std::cout << "wrong number of arguments for Jet Parameter Limits:" 
	      << limits.size() << '\n';
  }
  limits.clear();
  limits = bag_of<double>(config_->read<string>( "Tower Parameter Limits",""));
  
  if(limits.size() % 4 == 0) {
    for(unsigned int i = 0 ; i < limits.size() ; i += 4) {
      int index = (int)limits[i];
      for(int j = index; 
	  j <  par_->numberOfTowerParameters() ; 
	  j += par_->numberOfTowerParametersPerBin()) {
	par_limits.push_back(ParLimit(j,limits[i+1],limits[i+2],
				      limits[i+3]));
      }
    }
  } else if(limits.size() > 1) {
    std::cout << "wrong number of arguments for Tower Parameter Limits:" 
	      << limits.size() << '\n';
  }
}

ParameterLimitsReader::~ParameterLimitsReader()
{
}

int ParameterLimitsReader::readEvents(std::vector<Event*>& data)
{
  for(std::vector<ParLimit>::const_iterator pl = par_limits.begin() ;
      pl != par_limits.end() ; ++pl) {
    std::cout << " Adding limit for parameter " << std::flush;
    if     ( pl->index+1 < 10   ) std::cout << pl->index+1 << ":  " << std::flush;
    else if( pl->index+1 < 100  ) std::cout << pl->index+1 << ": "  << std::flush;
    else if( pl->index+1 < 1000 ) std::cout << pl->index+1 << ":"   << std::flush;
    std::cout << "  min: " << pl->min << "\t" << std::flush;
    std::cout <<  " max: " << pl->max << "\t" << std::flush;
    std::cout <<  " k: "   << pl->k << std::endl;
    data.push_back(new ParameterLimit(pl->index, par_, pl->min, pl->max, pl->k));
  }
  return par_limits.size();
}