Plot the diffrentail cross section

Asked by Ayoub Hmissou

Dear MadAnalysis team

I Need to plot the diffrentail with unit attobarn not picobarn in the axe of y

Question information

Language:
English Edit question
Status:
Answered
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Jack Y. Araz (jackaraz) said :
#1

Hi

You can include the unit conversion within luminosity value and calculate plot accordingly. Unit conversion is a universal transformation you can just multiply with the existing cross section and calculate the result in the way you want.

Cheers


Jack Y. Araz

On Nov 9, 2021, at 5:30 PM, Ayoub Hmissou <<email address hidden><mailto:<email address hidden>>> wrote:

[EXTERNAL EMAIL]

Question #699415 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699415

Description changed to:

Dear MadAnalysis team

I Need to plot the diffrentail with unit attobarn not picobarn in the
axe of y

--
You received this question notification because you are an answer
contact for MadAnalysis 5.

Revision history for this message
Ayoub Hmissou (ayb+) said :
#2

Please give me a example

Revision history for this message
Jack Y. Araz (jackaraz) said :
#3

Number of events per X GeV = xsec [pb] * 1000 * Lumi [1/fb] / X GeV

X is your bin size assuming you are plotting against mass or momentum etc.


Jack Y. Araz

> On Nov 9, 2021, at 5:45 PM, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699415 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699415
>
> Status: Answered => Open
>
> Ayoub Hmissou is still having a problem:
> Please give me a example
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Ayoub Hmissou (ayb+) said :
#4

Thank you
That give the diffrentail cross-section?

Revision history for this message
Jack Y. Araz (jackaraz) said :
#5

The differential cross-section is cross-section per bin.

--
Jack Y. Araz

> On 9 Nov 2021, at 6:25 pm, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699415 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699415
>
> Status: Answered => Open
>
> Ayoub Hmissou is still having a problem:
> Thank you
> That give the diffrentail cross-section?
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Ayoub Hmissou (ayb+) said :
#6

How to get it

Revision history for this message
Jack Y. Araz (jackaraz) said :
#7

I gave you everything you need, i think you should read a bit more and think about it.

--
Jack Y. Araz

> On 9 Nov 2021, at 9:01 pm, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699415 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699415
>
> Status: Answered => Open
>
> Ayoub Hmissou is still having a problem:
> How to get it
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Ayoub Hmissou (ayb+) said :
#8

Hi

for my example

set un_evt.xsection = 0.07092
set main.lumi = (1/1000)*0.25
plot Y(z ) 12 0 3

Number of events per X GeV = xsec [pb] * 1000 * Lumi [1/fb] / X GeV

with X=0.25 GeV in my example

so Number of evants per X = xsec [pb] * 1000*(1/1000)*0.25/0.25
                                                    = xec[pb]

that don't give dsigma/dY (pb/GeV) ?

Revision history for this message
Jack Y. Araz (jackaraz) said :
#9

Y is rapidity it does not have units of GeV please read the manual and the reference card for the normal mode which can be found under docs folder. All your questions are answered in the documentations available online.

--
Jack Y. Araz

> On 10 Nov 2021, at 8:30 pm, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699415 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699415
>
> Status: Answered => Open
>
> Ayoub Hmissou is still having a problem:
> Hi
>
> for my example
>
> set un_evt.xsection = 0.07092
> set main.lumi = (1/1000)*0.25
> plot Y(z ) 12 0 3
>
> Number of events per X GeV = xsec [pb] * 1000 * Lumi [1/fb] / X GeV
>
> with X=0.25 GeV in my example
>
> so Number of evants per X = xsec [pb] * 1000*(1/1000)*0.25/0.25
> = xec[pb]
>
> that don't give dsigma/dY (pb/GeV) ?
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Ayoub Hmissou (ayb+) said :
#10

i don't found that in documentation please cite the reference

Revision history for this message
Jack Y. Araz (jackaraz) said :
#11

As i said docs folder in your local madanalysis folder or ma5 home page. There are ton of tutorials as well.

--
Jack Y. Araz

> On 10 Nov 2021, at 9:01 pm, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699415 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699415
>
> Status: Answered => Open
>
> Ayoub Hmissou is still having a problem:
> i don't found that in documentation please cite the reference
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Jack Y. Araz (jackaraz) said :
#12
Revision history for this message
Ayoub Hmissou (ayb+) said :
#13

Thank you
In your example you set the lumi=1e-3*10 (10 is the bin size)
But you give me those relation
      Number of events per 10 GeV = xsec [pb] * 1000 * Lumi [1/fb] / 10 GeV
 So = xec[pb] isn't dsigma/dpT

Wath the difference between Number of events per bin size And Number of events per bin??
I don't understand

Revision history for this message
Jack Y. Araz (jackaraz) said :
#14

A bin has a size which the case of that example 10 GeV. Please read about histograms and how they are plotted. Your questions are not about MadAnalysis but about general knowledge of physics and can be found in any textbook or internet.

--
Jack Y. Araz

> On 13 Nov 2021, at 10:31 pm, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699415 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699415
>
> Status: Answered => Open
>
> Ayoub Hmissou is still having a problem:
> Thank you
> In your example you set the lumi=1e-3*10 (10 is the bin size)
> But you give me those relation
> Number of events per 10 GeV = xsec [pb] * 1000 * Lumi [1/fb] / 10 GeV
> So = xec[pb] isn't dsigma/dpT
>
> Wath the difference between Number of events per bin size And Number of events per bin??
> I don't understand
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Ayoub Hmissou (ayb+) said :
#15

Yes I know Wath you said

I have in my histogram number of Z (particle) =number Of events

But Which relation we can use

number of events = sigma X lumi X 1000

Or

Number of events per X GeV = xsec [pb] * 1000 * Lumi [1/fb] / X GeV

Revision history for this message
Jack Y. Araz (jackaraz) said :
#16

Hi

Please learn what histogram is.

--
Jack Y. Araz

> On 14 Nov 2021, at 4:15 am, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699415 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699415
>
> Status: Answered => Open
>
> Ayoub Hmissou is still having a problem:
> Yes I know Wath you said
>
> I have in my histogram number of Z (particle) =number Of events
>
> But Which relation we can use
>
> number of events = sigma X lumi X 1000
>
> Or
>
> Number of events per X GeV = xsec [pb] * 1000 * Lumi [1/fb] / X GeV
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Can you help with this problem?

Provide an answer of your own, or ask Ayoub Hmissou for more information if necessary.

To post a message you must log in.