RooBlindTools.cxx

Go to the documentation of this file.

/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitModels                                                     *
 * @(#)root/roofit:$Id$
 * Authors:                                                                  *
 *   AR, Aaron Roodman, Stanford University, roodman@slac.stanford.edu       *
 *   WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.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 RooBlindTools
    \ingroup Roofit
 
**/
 
#include "RooBlindTools.h"
 
#include "RooErrorHandler.h"
#include "strlcpy.h"
 
#include <iostream>
#include <fstream>
#include <cmath>
#include <cstring>
#include <cctype>
 
 
////////////////////////////////////////////////////////////////////////////////
/// Constructor
 
RooBlindTools::RooBlindTools(const char *stSeedIn, blindMode Mode,
              double centralValue, double sigmaOffset, bool s2bMode) :
 
  _PrecisionOffsetScale(sigmaOffset),
  _PrecisionCentralValue(centralValue),
  _mode(Mode),
  _s2bMode(s2bMode)
{
  setup(stSeedIn);
}
 
////////////////////////////////////////////////////////////////////////////////
/// Constructor
 
RooBlindTools::RooBlindTools(const RooBlindTools& blindTool):
  _PrecisionOffsetScale(blindTool.getPrecisionOffsetScale()),
  _PrecisionCentralValue(blindTool.getPrecisionCentralValue()),
  _mode(blindTool.mode()),
  _s2bMode(blindTool._s2bMode)
{
  setup(blindTool.stSeed());
}
 
////////////////////////////////////////////////////////////////////////////////
 
void RooBlindTools::setup(const char *stSeedIn)
{
  _stSeed = stSeedIn;
 
  _DeltaZScale = 1.56;
 
  _DeltaZOffset = _DeltaZScale*MakeOffset("abcdefghijklmnopqrstuvwxyz");
 
  _DeltaZSignFlip = MakeSignFlip("ijklmnopqrstuvwxyzabcdefgh");
 
  _AsymOffset = MakeGaussianOffset("opqrstuvwxyzabcdefghijklmn");
 
  _AsymSignFlip = MakeSignFlip("zyxwvutsrqponmlkjihgfedcba");
 
  _DeltaMScale = 0.1;
 
  _DeltaMOffset = _DeltaMScale*MakeOffset("opqrstuvwxyzabcdefghijklmn");
 
  _MysteryPhase = 3.14159 *
                  MakeOffset("wxyzabcdefghijklmnopqrstuv");
 
  if (_s2bMode) {
    _PrecisionSignFlip = MakeSignFlip("zyxwvutsrqponmlkjihgfedcba");
  } else {
    _PrecisionSignFlip = MakeSignFlip("klmnopqrstuvwxyzabcdefghij");
  }
 
  _PrecisionOffset = _PrecisionOffsetScale*MakeGaussianOffset("opqrstuvwxyzabcdefghijklmn");
 
  _PrecisionUniform = _PrecisionOffsetScale*MakeOffset("jihgfedcbazyxwvutsrqponmlk");
 
  _STagConstant = Randomizer("fghijklmnopqrstuvwxyzabcde");
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HideDeltaZ(double DeltaZ, double STag)const{
  Int_t sTag = SignOfTag(STag);
  double DeltaZPrime = _DeltaZSignFlip*DeltaZ*sTag + _DeltaZOffset;
 
  return DeltaZPrime;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HiDelZPdG(double DeltaZ, double STag, double PdG) const{
  Int_t sTag = SignOfTag(STag);
  double DeltaZPrime = _DeltaZSignFlip*(DeltaZ - PdG)*sTag + _DeltaZOffset;
 
  return DeltaZPrime;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHideDeltaZ(double DeltaZPrime, double STag) const{
  Int_t sTag = SignOfTag(STag);
  double DeltaZ = (DeltaZPrime - _DeltaZOffset)/(sTag*_DeltaZSignFlip);
 
  return DeltaZ;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHiDelZPdG(double DeltaZPrime, double STag, double PdG) const{
  Int_t sTag = SignOfTag(STag);
  double DeltaZ = PdG + (DeltaZPrime - _DeltaZOffset)/(sTag*_DeltaZSignFlip);
 
  return DeltaZ;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHideAsym(double AsymPrime) const{
  if(mode()==dataonly) return AsymPrime;
 
  double Asym = (AsymPrime - _AsymOffset)/_AsymSignFlip;
 
  return Asym;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HideAsym(double Asym) const{
  if(mode()==dataonly) return Asym;
 
  double AsymPrime = Asym*_AsymSignFlip + _AsymOffset;
 
  return AsymPrime;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHideDeltaM(double DeltaMPrime) const{
  if(mode()==dataonly) return DeltaMPrime;
 
  double DeltaM = DeltaMPrime - _DeltaMOffset;
 
  return DeltaM;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HideDeltaM(double DeltaM) const{
  if(mode()==dataonly) return DeltaM;
 
  double DeltaMPrime = DeltaM + _DeltaMOffset;
 
  return DeltaMPrime;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHiAsPdG(double AsymPrime, double PdG) const{
  if(mode()==dataonly) return AsymPrime;
 
  double Asym = PdG + (AsymPrime - _AsymOffset)/_AsymSignFlip;
 
  return Asym;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::MysteryPhase() const{
  if(mode()==dataonly) return 0.0;
 
  return _MysteryPhase;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HiAsPdG(double Asym, double PdG) const{
  if(mode()==dataonly) return Asym;
 
  double AsymPrime = (Asym - PdG)*_AsymSignFlip + _AsymOffset;
 
  return AsymPrime;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHidePrecision(double PrecisionPrime) const{
  if(mode()==dataonly) return PrecisionPrime;
 
  double Precision(0.);
 
  if (_PrecisionSignFlip>0) {
    Precision = PrecisionPrime - _PrecisionOffset;
  }
  else {
    Precision = 2.0*_PrecisionCentralValue - PrecisionPrime + _PrecisionOffset;
  }
 
 
  return Precision;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HidePrecision(double Precision) const{
  if(mode()==dataonly) return Precision;
 
  double PrecisionPrime(0.);
 
  if (_PrecisionSignFlip>0) {
    PrecisionPrime = Precision + _PrecisionOffset;
  }
  else {
    PrecisionPrime = 2.0*_PrecisionCentralValue - Precision + _PrecisionOffset;
  }
 
  return PrecisionPrime;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHideOffset(double PrecisionPrime) const{
  if(mode()==dataonly) return PrecisionPrime;
 
  return PrecisionPrime - _PrecisionOffset;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HideOffset(double Precision) const{
  if(mode()==dataonly) return Precision;
 
  return Precision + _PrecisionOffset;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::UnHideUniform(double PrecisionPrime) const{
  if(mode()==dataonly) return PrecisionPrime;
 
  return PrecisionPrime - _PrecisionUniform;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::HideUniform(double Precision) const{
  if(mode()==dataonly) return Precision;
 
  return Precision + _PrecisionUniform;
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::RandomizeTag(double STag, Int_t EventNumber) const{
  Int_t Seed = EventNumber % 7997 + 2;
  double r = PseudoRandom(Seed);
  double STagPrime(0.0);
 
  if (r < _STagConstant){
    STagPrime = STag;
  } else {
    STagPrime = -1.0 * STag ;
  }
 
  return STagPrime;
 
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::Randomizer(const char *StringAlphabet) const{
  char lowerseed[1024] ;
  strlcpy(lowerseed,_stSeed,1024) ;
 
  Int_t lengthSeed = strlen(lowerseed);
 
  for (Int_t j=0; j<lengthSeed; j++){
    lowerseed[j] =tolower(_stSeed[j]);
  }
  Int_t sumSeed = 0;
  for (Int_t i=0; i<lengthSeed; i++){
    for (Int_t iAlphabet=0; iAlphabet<26; iAlphabet++){
      if ( lowerseed[i] == StringAlphabet[iAlphabet] ){
   if (_s2bMode) {
     sumSeed =  (iAlphabet<<(5*(i%3)))^sumSeed;
   } else {
     sumSeed = sumSeed + iAlphabet ;
   }
      }
    }
  }
 
  if (lengthSeed<5 || ((sumSeed<1 || sumSeed>8000)&&!_s2bMode)) {
    std::cout<< "RooBlindTools::Randomizer: Your String Seed is Bad: '" << _stSeed << "'" << std::endl ;
    RooErrorHandler::softAbort() ;
  }
 
  Int_t ia = 8121;
  Int_t ic = 28411;
  Int_t im = 134456;
  UInt_t jRan = (sumSeed*ia + ic) % im;
 
  jRan = (jRan*ia + ic) % im;
  jRan = (jRan*ia + ic) % im;
  jRan = (jRan*ia + ic) % im;
 
  double theRan = (float) jRan / (float) im;
 
  return theRan;    //theRan is between 0.0 - 1.0
 
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::PseudoRandom(Int_t Seed) const{
  if (Seed<1 || Seed>8000 ) {
    std::cout<< "RooBlindTools::PseudoRandom: Your integer Seed is Bad" << std::endl;
  }
 
  Int_t ia = 8121;
  Int_t ic = 28411;
  Int_t im = 134456;
  UInt_t jRan = (Seed*ia + ic) % im;
 
  jRan = (jRan*ia + ic) % im;
  jRan = (jRan*ia + ic) % im;
  jRan = (jRan*ia + ic) % im;
 
  double theRan = (float) jRan / (float) im;
 
  return theRan;    //theRan is between 0.0 - 1.0
 
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::MakeOffset(const char *StringAlphabet) const{
  double theRan = Randomizer(StringAlphabet);
 
  double theOffset = (2.0)*theRan - (1.0);
 
  return theOffset;   //theOffset lies between -1.0 and 1.0
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::MakeGaussianOffset(const char *StringAlphabet) const{
  double theRan1 = Randomizer(StringAlphabet);
  double theRan2 = Randomizer("cdefghijklmnopqrstuvwxyzab");
 
  if (theRan1==0.0 || theRan1==1.0){
    theRan1 = 0.5;
  }
  if (theRan2==0.0 || theRan2==1.0){
    theRan2 = 0.5;
  }
 
  double theOffset = sin(2.0*3.14159*theRan1)*sqrt(-2.0*log(theRan2));
 
  return theOffset;   //theOffset is Gaussian with mean 0, sigma 1
}
 
////////////////////////////////////////////////////////////////////////////////
 
double RooBlindTools::MakeSignFlip(const char *StringAlphabet) const{
  double theRan = Randomizer(StringAlphabet);
 
  double theSignFlip = 1.0;
  if (theRan>0.5){
    theSignFlip = 1.0;
  } else {
    theSignFlip = -1.0;
  }
 
  return theSignFlip;  //theSignFlip is = +1 or -1
}
 
////////////////////////////////////////////////////////////////////////////////
 
Int_t RooBlindTools::SignOfTag(double STag) const{
  Int_t sTag;
  if (STag < 0.0){
    sTag = -1;
  }
  else if (STag > 0.0) {
    sTag = 1;
  }
  else {
    sTag = 1;
  }
 
  return sTag;
 
}