exclusion_CLs

Dear Ahmed,

Could you please be more specific. What is the question exactly? I do not see what you need.

Regards,

Benjamin

On 16 May 2015, at 24:12 , ahh <email address hidden> wrote:

> New question #266966 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/266966
>
> Dear MA authors
> ​Is there any example of XML.info file to mimic it in my analysis ????
>
> Thanks
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Dear dr Benjamin
Iam sorry if i made it fuzzy,
I just used ./exclusion_CLs.py
After, the message is
Usage: ./exclusion_CLs.py analysis_name benchmark_point [run_number] [cross section in pb]
I just ask about the file which i put my benchmarks "erros, number of observation , the number of bkg,...." is this file generated automatically and i just put this variables in it or shall i construct it from scratch ?????

Many thanks

Dear Ahmed,

Please have a look to http://madanalysis.irmp.ucl.ac.be/wiki/PhysicsAnalysisDatabase where you will be able to get more information and links to examples of files such as the one you are asking about.

Regards,

Benjamin

On 16 May 2015, at 12:21 , ahh <email address hidden> wrote:

> Question #266966 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/266966
>
> Status: Answered => Open
>
> ahh is still having a problem:
> Dear dr Benjamin
> Iam sorry if i made it fuzzy,
> I just used ./exclusion_CLs.py
> After, the message is
> Usage: ./exclusion_CLs.py analysis_name benchmark_point [run_number] [cross section in pb]
> I just ask about the file which i put my benchmarks "erros, number of observation , the number of bkg,...." is this file generated automatically and i just put this variables in it or shall i construct it from scratch ?????
>
> Many thanks
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Dear Ahmed,

You can find the answer to your question more easily on pages 14-15 in http://arxiv.org/abs/1407.3278
Also, if you download any of the existing analysis codes (for example http://inspirehep.net/record/1301484/ ) you will find an .info file that you can use as a template.

Best,
Sabine

Dear Authors
Thank you very much, but what if i tried to calculate the exclusion limit at 13 TeV where no data observed.
Shall i use signal instead or signal+background or i can still use the data collected at 8 TeV ???
Many thanks

Dear Ahmed,

I am not sure to fully understand the question. Could you please be a little bit more explicit? Why do you want to ignore the background?

Cheers,

Benjamin

Dear dr Benjamin
The code takes as an input:
1- lumi 2- Number of background <nb> 3- number of observed events <nobs> 5- Uncertainty
And the question is :
Now iam simulating at 13 TeV, how can i set number of observed <nobs> ???

 Many thanks

nobs means number of observed events. there is no nobs for 13 TEV yet. and
also nb will change will change at the higher energy. you can do a rough
estimate by rescaling the 8 TeV numbers by the 13TeV cross section, but as
said this is rough.

best, sabi e

On Saturday, June 20, 2015, ahh <email address hidden>
wrote:

> Question #266966 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/266966
>
> Status: Needs information => Open
>
> ahh gave more information on the question:
> Dear dr Benjamin
> The code takes as an input:
> 1- lumi 2- Number of background <nb> 3- number of observed events <nobs>
> 5- Uncertainty
> And the question is :
> Now iam simulating at 13 TeV, how can i set number of observed <nobs> ???
>
> Many thanks
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>

--
_______________________________________________________
Sabine Kraml - <email address hidden> - www.kraml.net

LPSC
53 Av des Martyrs
38026 Grenoble
France
(+33)(0)4 76 28 40 52

dear prof Sabine
Thank you very much, indeed that what i was asking for.
Just one last question:
In case to calculate the expected mean shall i replace "weak" by "mean" in
scipy.stats.percentileofscore(ToyBGs, NumObserved, kind='weak')*.01
for both p_b and p_SplusB functions. then the return is expected median limit ????

Many thanks again for your help

Hi Ahmed,

You do not have to do this. See here:
  http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.stats.percentileofscore.html

Cheers,

Benjamin