RooUnblindPrecision.cxx

Go to the documentation of this file.

/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitModels                                                     *
 * @(#)root/roofit:$Id$
 * 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)             *
 *****************************************************************************/
 
/** \class RooUnblindPrecision
    \ingroup Roofit
 
Implementation of BlindTools' precision blinding method
A RooUnblindPrecision object is a real valued function
object, constructed from a blind value holder and a
set of unblinding parameters. When supplied to a PDF
in lieu of a regular parameter, the blind value holder
supplied to the unblinded objects will in a fit be minimized
to blind value corresponding to the actual minimum of the
parameter. The transformation is chosen such that the
the error on the blind parameters is identical to that
of the unblind parameter
**/
 
#include "RooFit.h"
 
#include "RooArgSet.h"
#include "RooUnblindPrecision.h"
 
 
using namespace std;
 
ClassImp(RooUnblindPrecision);
 
////////////////////////////////////////////////////////////////////////////////
/// Default constructor
 
RooUnblindPrecision::RooUnblindPrecision()
{
}
 
////////////////////////////////////////////////////////////////////////////////
/// Constructor from a given RooAbsReal (to hold the blind value) and a set of blinding parameters
 
RooUnblindPrecision::RooUnblindPrecision(const char *name, const char *title,
                const char *blindString, Double_t centralValue,
                Double_t scale, RooAbsReal& value,
                Bool_t sin2betaMode)
  : RooAbsHiddenReal(name,title),
  _value("value","Precision blinded value",this,value),
  _blindEngine(blindString,RooBlindTools::full,centralValue,scale,sin2betaMode)
{
}
 
////////////////////////////////////////////////////////////////////////////////
/// Constructor from a given RooAbsReal (to hold the blind value) and a set of blinding parameters
 
RooUnblindPrecision::RooUnblindPrecision(const char *name, const char *title,
                const char *blindString, Double_t centralValue,
                Double_t scale, RooAbsReal& value, RooAbsCategory& blindState,
                Bool_t sin2betaMode)
  : RooAbsHiddenReal(name,title,blindState),
  _value("value","Precision blinded value",this,value),
  _blindEngine(blindString,RooBlindTools::full,centralValue,scale,sin2betaMode)
{
}
 
////////////////////////////////////////////////////////////////////////////////
/// Copy constructor
 
RooUnblindPrecision::RooUnblindPrecision(const RooUnblindPrecision& other, const char* name) :
  RooAbsHiddenReal(other, name),
  _value("asym",this,other._value),
  _blindEngine(other._blindEngine)
{
}
 
////////////////////////////////////////////////////////////////////////////////
/// Destructor
 
RooUnblindPrecision::~RooUnblindPrecision()
{
}
 
////////////////////////////////////////////////////////////////////////////////
/// Evaluate RooBlindTools unhide-precision method on blind value
 
Double_t RooUnblindPrecision::evaluate() const
{
  if (isHidden()) {
    // Blinding active for this event
    return _blindEngine.UnHidePrecision(_value);
  } else {
    // Blinding not active for this event
    return _value ;
  }
}