missing energy pdg code MA4

Hi,
I have unfortunately no clue about madanalysis 4, so that the only way to get an answer to your question is to ask Rikkert.
Cheers
Benjamin

-------- Original message --------
From: Rafael Marcelino <email address hidden>
Date: 21/11/2015 01:52 (GMT+01:00)
To: <email address hidden>
Subject: [Question #275132]: missing energy pdg code MA4

New question #275132 on MadAnalysis 5:
https://answers.launchpad.net/madanalysis5/+question/275132

Hi!

I know this question it's not about madana5, but i can't find a way to ask Rikkert.

I have an ufo model of Little Higgs with T symmetry. In this model i did set the dark matter candidate with pdg code 5100002. After pythia6 and pgs i got some good plots of missing et, in the index.html of madgraph. My problem is when i did the plots with my MadAnalysis4, i can't get correct distributions for missing et, even setting mET class as 5100002.

When I hack the ufo code generation and change the pdg code for the dark matter candidate to light neutralino like, 1000022, in pythia run i got an error. So it's important to leave the pdg code as 5100002.

I just have to include 5100002 in the mET class on ma_card or there's something else?

Appreciate your attention and sorry for the madanalysis4 question.

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

Dear Rafael,

I don't know how your particle is identified after the PGS run. What is the 'type' in the .lhco file for your particle with PDG code 5100002 ?

In .lhco file reader in MadAnalysis, (which is at the end of the plot_events.f file), there are only 6 types allowed. The relevant code is:

         if (type(i).eq.0)then
            ic(1,i)=22 !photon
  elseif (type(i).eq.1.and.ntrk(i).eq.-1.0)then
            ic(1,i)=11 !electron
         elseif (type(i).eq.1.and.ntrk(i).eq.1.0)then
            ic(1,i)=-11 !positron
         elseif (type(i).eq.2.and.ntrk(i).eq.-1.0)then
            ic(1,i)=13 !mu-
         elseif (type(i).eq.2.and.ntrk(i).eq.1.0)then
            ic(1,i)=-13 !mu+
         elseif (type(i).eq.3.and.ntrk(i).le.-1.0)then
            ic(1,i)=15 !tau-
         elseif (type(i).eq.3.and.ntrk(i).ge.1.0)then
            ic(1,i)=-15 !tau+
         elseif (type(i).eq.4)then
            if (btag(i).ge.1) then
               ic(1,i)=5 !bjet
            else
               ic(1,i)=21 !non-bjet (or gluon)
            endif
         elseif (type(i).eq.6)then
            ic(1,i)=12 !Missing E_t (or neutrino)
         endif

Probably your particle doesn't fit in any of the types. By checking what the type is in the .lhco event file, it should be straight-forward to add it to the "Missing E_t" type.

Let me know if this helps fixing it.

Best regards,
Rikkert

Thank you Rikkert for your fast answer.

I don't know how your particle is identified after the PGS run. What is the 'type' in the .lhco file for your particle with PDG code 5100002 ?

I'm not sure how to check this, besides the pdg code 5100002, i didn't choose any "type" number for particles.

I uploaded one of my .lhco files in a dropbox url. Could you check this lhco?

https://www.dropbox.com/s/e55tw83m5bbdozo/qhqh_1.4_1.5_0.3_pgs_events.lhco.gz?dl=0

Dear Rafael,

Nearly all the particles in your LHCO events are of type "4" (which are jets) or type "6" (which is missing ET). The latter is dealt correctly in the code (I pointed you to the relevant part in my previous message).

I don't know to what Pythia/PGS does with your particle and to which type it assigns it; at least all the particles that are there should be plotted correctly. Given that you get the right plots on-line, you should get the same plots off-line. Why do you say that you do not get the correct plots when running MadAnalysis locally?

Best regards,
Rikkert

Dear Rafael,

Nearly all the particles in your LHCO events are of type "4" (which are jets) or type "6" (which is missing ET). The latter is dealt correctly in the code (I pointed you to the relevant part in my previous message).

I don't know to what Pythia/PGS does with your particle and to which type it assigns it; at least all the particles that are there should be plotted correctly. Given that you get the right plots on-line, you should get the same plots off-line. Why do you say that you do not get the correct plots when running MadAnalysis locally?

Best regards,
Rikkert

Dear Rafael,

Nearly all the particles in your LHCO events are of type "4" (which are jets) or type "6" (which is missing ET). The latter is dealt correctly in the code (I pointed you to the relevant part in my previous message).

I don't know to what Pythia/PGS does with your particle and to which type it assigns it; at least all the particles that are there should be plotted correctly. Given that you get the right plots on-line, you should get the same plots off-line. Why do you say that you do not get the correct plots when running MadAnalysis locally?

Best regards,
Rikkert

Rikkert,

I think my main problem it's the ufo model, not the mad4 routines. Anyway, thanks for your help here.