Ploting the Differential cross section

Hi

This is not a default usage of the module, an easy workaround is setting luminosity to one which will give you dsigma/dpt in fb but it wont be written like that in the y axis label. You can then go to the outputs folder and modify the code that is prepared for plotting to get the output that you want.

cheers

—
Jack Y. Araz

> On Mar 25, 2021, at 2:50 PM, Jossif <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> New question #696245 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/696245
>
>
> Dear MadAnalysis team,
>
> I need to plot a differential cross section . I would like to have a differential cross section dsigma/d PT on the y-axis and PT on the x-axis .
>
> Thanks a lot.
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Hi Jack ,
can you give me an exmeple for that .how to set the luminosity
Thank you .

Hi

From madanalysis interface type
set main.lumi = 1

PS: MadAnalysis does have a manual. See 1808.00480 and 2006.09387.

Cheers

—
Jack Y. Araz

> On Mar 25, 2021, at 3:25 PM, Jossif <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #696245 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/696245
>
> Status: Answered => Open
>
> Jossif is still having a problem:
> Hi Jack ,
> can you give me an exmeple for that .how to set the luminosity
> Thank you .
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Thank you so much for your answer!
For this example:
N=100000
set un_evt.xsection = 483.2
set main.lumi =( 1/1000)*12.5
plot PT(b[1]) 40 0 500 [logY]
To obtain the differential cross section in y-axis in units of pb/GeV
right?

no.

set main.lumi = 1e-3

—
Jack Y. Araz

> On Mar 25, 2021, at 3:41 PM, Jossif <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #696245 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/696245
>
> Status: Answered => Open
>
> Jossif is still having a problem:
> Thank you so much for your answer!
> For this example:
> N=100000
> set un_evt.xsection = 483.2
> set main.lumi =( 1/1000)*12.5
> plot PT(b[1]) 40 0 500 [logY]
> To obtain the differential cross section in y-axis in units of pb/GeV
> right?
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

But for the bin size 12.5 GeV.

Luminosity is just scaling the distribution does not change the bin size. Please read the manual.

—
Jack Y. Araz

> On Mar 25, 2021, at 3:50 PM, Jossif <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #696245 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/696245
>
> Status: Answered => Open
>
> Jossif is still having a problem:
>
> But for the bin size 12.5 GeV.
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Yes but how we get the differential cross section

If you set luminosity to 1e-3 that will give you the differential cross section automatically

number of events = sigma X lumi X 1000

if lumi is 1 you get the sigma [fb]. Is this clear?

A histogram means sigma is divided into dXXX bins hence when you plot number of events its number of events per bin so when you set lumi to one you get sigma per bin. These are all written in the manual.

—
Jack Y. Araz

> On Mar 25, 2021, at 4:05 PM, Jossif <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #696245 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/696245
>
> Status: Answered => Open
>
> Jossif is still having a problem:
> Yes but how we get the differential cross section
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Yes but how we get the dsigma/d PT

I might have misunderstood the problem that you are having I was giving you the recipe to get pb/XX GeV. So yes if you use the lumi=( 1/1000)*12.5 you will get pb/GeV

jack

so si on use lumi=( 1/1000)*12.5 we get dsigma/d PT (pb/GeV) on the y-axis and PT on the x-axis ???

Hi

It won't be written on the label of the plot but the content will be dsigma/dPT. Ma5 writes the code to create the plots automatically in your output folder so you can modify them to your liking.

Cheers
Jack

Hi
Yes i know the writing is modified by the output.
but i want to get a distribution of the dsigma/dPT on the y-axis and PT on the x-axis.
Can I send you a image of what I want in your mail?
Thank you.

Hi

As I said the simple script that generates those plots in your report is written under your Output folder. You can go in and modify that particular plot that you want to change, run the code and generate the plot.

Cheers

--
Jack Y. Araz

> On 18 Sep 2021, at 9:11 am, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #696245 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/696245
>
> Ayoub Hmissou requested more information:
>
> Hi
> Yes i know the writing is modified by the output.
> but i want to get a distribution of the dsigma/dPT on the y-axis and PT on the x-axis.
> Can I send you a image of what I want in your mail?
> Thank you.
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

For this example:
N=1000
set un_evt.xsection = 483.2
set main.lumi =( 1/1000)*12.5
plot PT(b[1]) 40 0 500 [logY]
We get a distribution of dsigma/dPT

And
         si
          set main.lumi = 1e-3 we get a distribution
           of dsigma
right?

Yes

—
Jack Y. Araz

> On Sep 18, 2021, at 2:01 PM, Ayoub Hmissou <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #696245 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/696245
>
> Ayoub Hmissou requested more information:
> For this example:
> N=1000
> set un_evt.xsection = 483.2
> set main.lumi =( 1/1000)*12.5
> plot PT(b[1]) 40 0 500 [logY]
> We get a distribution of dsigma/dPT
>
> Andi si
> set main.lumi = 1e-3 we get a distribution
> of dsigma
> right?
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.