Accessing the MC level information of a sample processed with DelphesMA5tune

Asked by Nicolas Deutschmann

Hi

I am trying to write an analysis where I process a .hep file produced by Pythia with the MA5 Delphes tune and I would like to get the cross-section of the process to normalize my results. I understand that this information is no longer available in the LHCO output from Delphes but is there a way to get the cross-section from within the MA5 analysis since I feed it the .hep file as input ?(this information is in it right ?)

Could I for example call the function that is used internally to read the .hep file and fill in the info in sample.mc() ?

Thanks for your help,
Nicolas.

Question information

Language:
English Edit question
Status:
Solved
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Solved by:
Nicolas Deutschmann
Solved:
Last query:
Last reply:
Revision history for this message
Benjamin Fuks (fuks) said :
#1

Hi Nicolas,

The only option (without hacking the code) that is so far available is to renormalize your results offline by the correct cross section and luminosity.

Cheers

Benjamin

On 10 Sep 2014, at 15:06 , Nicolas Deutschmann <email address hidden> wrote:

> New question #254303 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/254303
>
> Hi
>
> I am trying to write an analysis where I process a .hep file produced by Pythia with the MA5 Delphes tune and I would like to get the cross-section of the process to normalize my results. I understand that this information is no longer available in the LHCO output from Delphes but is there a way to get the cross-section from within the MA5 analysis since I feed it the .hep file as input ?(this information is in it right ?)
>
> Could I for example call the function that is used internally to read the .hep file and fill in the info in sample.mc() ?
>
> Thanks for your help,
> Nicolas.
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Nicolas Deutschmann (ndeutschmann) said :
#2

Hi Benjamin,

I was guessing so. Do you think a hack would be simple enough to be worth trying ? If so where in the code should I look at to get started ?

Cheers,
Nicolas.

Revision history for this message
Benjamin Fuks (fuks) said :
#3

Hi Nicolas,

I think the easiest might be to include the cross section in the filename and to ask your analysis code to decode the filename (playing with strings).

Cheers,

Benj

Revision history for this message
Nicolas Deutschmann (ndeutschmann) said :
#4

Hi Benjamin,

I'll try something along those lines.

Thank you for your help !

Nicolas.