Unable to blind with RooBlindPrecision

Roofit and RooStats
1

Dear Experts,
I am trying to check the RooUnBlindPrecision on my simulation to check how it works (I know what the unblinded value is). Here’s my working directory blinding_forum.zip (2.6 MB). Please check the macro from line# 197

However, I not able to blind the results. I am trying to do:

Bool_t isBlind(kTRUE);
  //Bool_t isBlind(kFALSE);

  RooRealVar Araw_sig("A_{sig}","Araw",0,-1,1);
  TString blind("blind"), unblind("unblind");
  RooCategory Araw_sig_blind("Araw_sig_blind","Araw_sig_blind blinding state");  
  Araw_sig_blind.defineType("Unblind",0);
  Araw_sig_blind.defineType("Blind",1);
  
  if(isBlind)
      Araw_sig_blind.setLabel(blind);
  else
      Araw_sig_blind.setLabel(unblind);

  // Unblinding transformation
 RooUnblindPrecision Araw_sig_unblind("Araw_sig_unblind","Unblinded Araw_sig","BlindString",Araw_sig.getVal(),10,Araw_sig,Araw_sig_blind,0);
 RooFormulaVar signal_yield_d0("ND0_{signal}","(@0*(1 + @1))/2",RooArgList(signal_yield,Araw_sig_unblind));
RooFormulaVar signal_yield_d0bar("ND0bar_{signal}","(@0*(1 - @1))/2",RooArgList(signal_yield,Araw_sig_unblind));

Will you please have a look ??

Regards,
Sanjeeda
2

Hi @sanjeeda,

I guess @moneta or @jonas can help here.

Cheers,
J.

3

This topic was automatically closed 14 days after the last reply. New replies are no longer allowed.