invariant mass distribution plot

Asked by niujj

hello
I am facing the following problem with MG5_aMC_v2.2.3.
I have root (ROOT 5.32) installed.
I installed MG5, ExRootAnalysis, MadAnalysis and Pythia-pgs.

I runned a process p p > a a. The invariant mass distribution plots in the delphes are not what i want. What I want is m_aa invariant mass distribution.
 How can i modify the plot_card.dat?

And i want to konw what(jet1,jet2, jet 3,jet 4,a1 ,a2 ,b1,b2 ... ) they mean in the plot?

Please help me ,
Thanks in advance for all your help,
Best Regards,
niujj

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Rikkert Frederix Edit question
Solved by:
niujj
Solved:
Last query:
Last reply:
Revision history for this message
Rikkert Frederix (frederix) said :
#1

Dear niujj,

To plot the M_aa, you have to create a separate class for the photons and then specify that you want to plot the invariant mass of the two hardest particles in that class.

The jet1...jet4 are the four hardest jets. 'a' are photons and 'b' are the b-jets.

Best,
Rikkert

Revision history for this message
niujj (501852601-j) said :
#2

Hi Rikkert,

I did just as you said. I create a separate class for the photons,just as
 a 6 7 # Class number 3
and
mij 3 2 # use the first two particles in the 3d class to plot inv. mass
but it made no sense.
And then i set pythia and delphes off, There is plots of a'pt and y .but no invariant mass distribution. I wonder why and how.
There is a other question.I generate p p > h z ,then i set the mass of mh: range(170,180,2) ,but report a bug,could not convert string to float:range(120,130,2). I want to know if the format is right?

thank you very much ,
niujj

Revision history for this message
Rikkert Frederix (frederix) said :
#3

Dear niujj,

The PDG code for a photon is 22, (6 is the top quark and 7 doesn't exist). Hence it's

a 22 # Class number 3

to add the photons to the class with name 'a'.

For the parameter scan, see here:
https://answers.launchpad.net/mg5amcnlo/+faq/2735

Best,
Rikkert

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#4

Dear niujj,

I guess that you are missing the
scan: flag before the range(170,180,2)

the line should be
   25 scan:range(170,180,2)

Cheers,

Olivier
> On Apr 6, 2016, at 08:27, Rikkert Frederix <email address hidden> wrote:
>
> Question #289929 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/289929
>
> Status: Open => Answered
>
> Rikkert Frederix proposed the following answer:
> Dear niujj,
>
> The PDG code for a photon is 22, (6 is the top quark and 7 doesn't
> exist). Hence it's
>
> a 22 # Class number 3
>
> to add the photons to the class with name 'a'.
>
> For the parameter scan, see here:
> https://answers.launchpad.net/mg5amcnlo/+faq/2735
>
> Best,
> Rikkert
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
niujj (501852601-j) said :
#5

Dear Rikkert and Olivier,

Thank you a lot for the problem of scan. I did miss the scan.
 It really help me save more time and energe by using the scan.

The first ploblem is still there. I add the photons as you said,
a 22 # Class number 3
However there is still not m_aa plot,
Anywhere eles need to be modified?

Thanks again,
niujj

Revision history for this message
Rikkert Frederix (frederix) said :
#6

Dear Niujj,

I don't know what you are doing wrong. If I do

generate p p > a a j
output
launch --laststep=parton

Then in the plot_card, I add the following (line 84):

a 22

So that 'a' becomes the third available class (the 1st class remains the jets, and 2nd class the b-quarks).
Around line 152 I add:

mij 3 2

to add the invariant mass plot for the two hardest particles in the 3rd class to the list of plots.

And I start the run and get the invariant mass plot for the photon pair. Could you please check that this works in your case?

best regards,
Rikkert

Revision history for this message
niujj (501852601-j) said :
#7

Dear Rikkert,

I miss the modification of line 152,So this time i get the expected plot.

Thanhs again for you legible explanation.

best wishes,
niujjj

Revision history for this message
Nico (nico818) said :
#8

Hello,
I worked on MG5_aMC_2_5_4 and root 6.08, but MG5 plot card is not open.
Actually, plot card is not change.
I just generate p p > zh, so how I do invariant mass plot in root data analysis.
How can I work in root?

Please help me.

Best,
Bulgaa

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#9

Hi,

I'm not sure to understand what you ask for.
The situation is also highly different depending which kind of generation you are trying to do.

For LO,
You can either use MadAnalysis4 or MadAnalysis5 for having the plot.
You can install those tools by typing (inside MG5aMC command interface): "install MadAnalysis" or "install MadAnalysis5"
MadAnalysis5 is much more performant/flexible than MA4.

You can also decide to convert your output to root file.
For this you have to install another external package called ExRootAnalysis.
This will just convert the output file in root format. We do not do any plot, but you can obviously use that file to do any plot that you want.

Note that installing is not enough. When you run your code, you will be asked which one to use inside the first question.
You can answer analysis=MADANALYSIS_5, analysis=MADANALYSIS_4 , analysis=EXROOTANALYSIS
Depending of the one that you want to use.

For NLO, it also depends if you run NLO+PS or fixed order NLO
At fix-order you have two choice:
1) you can ask to generate events (following a lhe format), there we do not do plot automatically but you can use any tools externally (include those above) to do your plot
2) you can ask to generate the plot on the flight by writing your analysis (either in root or in fortran)

For NLO+PS, you have the MadAnalysis_5 option available (MA4 will also do the plot but those are not using a NLO save jet algorithm so those are certainly not correct for publication)

Cheers,

Olivier

> On 14 Sep 2017, at 11:33, bulgaa <email address hidden> wrote:
>
> Question #289929 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/289929
>
> bulgaa posted a new comment:
> Hello,
> I worked on MG5_aMC_2_5_4 and root 6.08, but MG5 plot card is not open.
> Actually, plot card is not change.
> I just generate p p > zh, so how I do invariant mass plot in root data analysis.
> How can I work in root?
>
>
> Please help me.
>
>
> Best,
> Bulgaa
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.