Question #699618 “missing energy ” : Questions : MadAnalysis 5

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-24:

#1

Hi

Delphes applies reconstruction efficiencies which ends up removing particles, hence the met.

—
Jack Y. Araz

> On Nov 24, 2021, at 2:15 PM, SOUAD SEMLALI <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> New question #699618 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Dear madanalaysis team,
>
> I am generating events for a signal with 2b 2j mu+ mu- in final state (Hsm -> z*z* h; h > bb~, z> mu+ mu- and z > jj). At parton level, the plots of Pt, eta and invariant mass distributions are generated without a problem; the distributions of the mass invariant showed that mbb ~ mh and m_{jj mu+ mu- bb}~ m_{Hsm}.
>
> But, the plots after hadronization and detector simulation are a bit confusing. I am getting high missing ET for a signal without missing ET. Furthermore, the peaks around mh and mH have disappeared. The delphes card already includes a NeutrinoFilter. I am not sure about the source of the problem.
>
> Many thanks for any suggestions or help!
>
> Best regards
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-24:

#2

Dear Jack,

Many thanks for your email. I am not sure I understand very well your answer. Are you suggesting that there are cuts in Delphes that can be applied to remove the met. please correct me if I am wrong?

My problem is that I am having a missing ET for a signal without missing ET, which I find really confusing. I couldn't figure out the source of the high met.

Many thanks for any help or suggestions

Best regards

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-24:

#3

Hi

Assuming that you are not doing anything wrong in delphes, it is a detector simulator a detector is not completely efficient and its job is to apply these inefficiencies to your event. But this is not a MadAnalysis related question, MadAnalysis just reads Delphes output. I suggest asking to Delphes authors for more details. Or if its too complicated to understand you can use MadAnalysis’ internal detector simulator arXiv:2006.09387.

Cheers

—
Jack Y. Araz

> On Nov 24, 2021, at 4:11 PM, SOUAD SEMLALI <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Many thanks for your email. I am not sure I understand very well your
> answer. Are you suggesting that there are cuts in Delphes that can be
> applied to remove the met. please correct me if I am wrong?
>
> My problem is that I am having a missing ET for a signal without missing
> ET, which I find really confusing. I couldn't figure out the source of
> the high met.
>
> Many thanks for any help or suggestions
>
> Best regards
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-24:

#4

Hi Jack,

Many thanks for your reply. I will write to Delphes team. In the meantime, I will resort to the simplified fast detector in madanalysis

Best regards

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-25:

#5

Dear Jack,

I have double-checked my results. The missing ET appeared after the hadronization (pythia8), which means that the problem is not related to Delphes, but I still can't figure out the source of the MET. I didn't make any changes to Pythia's card.

Many thanks for any suggestions or help!

Best regards
Souad

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-25:

#6

Hi Souad

I haven’t developed pythia... you can analyze the hadronic output through ma5 without reconstruction. Maybe some of the hadrons are radiating neutrinos but it should not be too large. You can plot neutrino distribution in hadron mode (./bin/ma5 -H) and see if this is the same/similar as your met. Additionally you can do truth level reconstructed analysis to compare your results to each other to track down the origin of the missing energy. You can also plot PID’s to see which particles are there where you shouldn’t have too many neutrinos in principle.

—
Jack Y. Araz

> On Nov 25, 2021, at 3:30 PM, SOUAD SEMLALI <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> SOUAD SEMLALI gave more information on the question:
> Dear Jack,
>
> I have double-checked my results. The missing ET appeared after the
> hadronization (pythia8), which means that the problem is not related to
> Delphes, but I still can't figure out the source of the MET. I didn't
> make any changes to Pythia's card.
>
> Many thanks for any suggestions or help!
>
> Best regards
> Souad
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-25:

#7

Hi Jack,

Many thanks for your reply. I have analyzed the hadronic output through ma5 without reconstruction. The MET is similar, and it's indeed too large. Is there a specific command to do the truth level reconstruction? or any reference?

Many thanks for your help to alleviate the confusion.

Best,
Souad

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-25:

#8

Hi

If hadron level analysis is giving you a large MET you will get the same thing from reconstructed analysis. I suggest plotting NPID to see which particles are in your final state event and/or directly read the HEPMC file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI <<email address hidden><mailto:<email address hidden>>> wrote:

Many thanks for your reply. I have analyzed the hadronic output through

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-26:

#9

Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles). I am not sure about the reason.

However, the root file (selection_0.C) associated to the plot showed this:

  stack->GetYaxis()->SetLabelSize(0.04);
  stack->GetYaxis()->SetLabelOffset(0.005);
  stack->GetYaxis()->SetTitleSize(0.06);
  stack->GetYaxis()->SetTitleFont(22);
  stack->GetYaxis()->SetTitleOffset(1);
  stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

  // X axis
  stack->GetXaxis()->SetLabelSize(0.04);
  stack->GetXaxis()->SetLabelOffset(0.005);
  stack->GetXaxis()->SetTitleSize(0.06);
  stack->GetXaxis()->SetTitleFont(22);
  stack->GetXaxis()->SetTitleOffset(1);
  stack->GetXaxis()->SetTitle("|NPID| ");
  stack->GetXaxis()->SetBinLabel(1,"e+/e-");
  stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
  stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
  stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
  stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
  stack->GetXaxis()->SetBinLabel(6,"a");
  stack->GetXaxis()->SetBinLabel(7,"K_L0");
  stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
  stack->GetXaxis()->SetBinLabel(9,"K-/K+");
  stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
  stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>>>
> wrote:
>
> Many thanks for your reply. I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles).  I am not sure about the reason.

However, the root file (selection_0.C)  associated to the plot showed this:

stack->GetYaxis()->SetLabelSize(0.04);
  stack->GetYaxis()->SetLabelOffset(0.005);
  stack->GetYaxis()->SetTitleSize(0.06);
  stack->GetYaxis()->SetTitleFont(22);
  stack->GetYaxis()->SetTitleOffset(1);
  stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
  stack->GetXaxis()->SetLabelSize(0.04);
  stack->GetXaxis()->SetLabelOffset(0.005);
  stack->GetXaxis()->SetTitleSize(0.06);
  stack->GetXaxis()->SetTitleFont(22);
  stack->GetXaxis()->SetTitleOffset(1);
  stack->GetXaxis()->SetTitle("|NPID| ");
  stack->GetXaxis()->SetBinLabel(1,"e+/e-");
  stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
  stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
  stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
  stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
  stack->GetXaxis()->SetBinLabel(6,"a");
  stack->GetXaxis()->SetBinLabel(7,"K_L0");
  stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
  stack->GetXaxis()->SetBinLabel(9,"K-/K+");
  stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
  stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>
> wrote:
>
> Many thanks for your reply.  I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-26:

#10

Hi

Portion of the code you sent just showing you x-axis labels the information is taken from histos.saf file within `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run number>/Histograms/histos.saf` if the plot that you are looking for is empty there than there is a problem. but if it has information not translated to root there might be some interference between root and ma5 so you might want to change the plotting renderer (explained how in FAQ for normal mode for a different example). It might also be the case that for some reason ma5 can not read the xsec information from your hepmc file (if you are giving root file there is no xsec information so you need to feed it by hand) so I suggest setting sample xsec as well to be on the safe side. For how to do these you can refer to normal mode reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI <<email address hidden><mailto:<email address hidden>>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles). I am not sure about the reason.

However, the root file (selection_0.C) associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>
wrote:

Many thanks for your reply. I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Hi

Portion of the code you sent just showing you x-axis labels the information is taken from histos.saf file within `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run number>/Histograms/histos.saf` if the plot that you are looking for is empty there than there is a problem. but if it has information not translated to root there might be some interference between root and ma5 so you might want to change the plotting renderer (explained how in FAQ for normal mode for a different example). It might also be the case that for some reason ma5 can not read the xsec information from your hepmc file (if you are giving root file there is no xsec information so you need to feed it by hand) so I suggest setting sample xsec as well to be on the safe side. For how to do these you can refer to normal mode reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI <question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles).  I am not sure about the reason.

However, the root file (selection_0.C)  associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
 stack->GetYaxis()->SetLabelOffset(0.005);
 stack->GetYaxis()->SetTitleSize(0.06);
 stack->GetYaxis()->SetTitleFont(22);
 stack->GetYaxis()->SetTitleOffset(1);
 stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
 stack->GetXaxis()->SetLabelSize(0.04);
 stack->GetXaxis()->SetLabelOffset(0.005);
 stack->GetXaxis()->SetTitleSize(0.06);
 stack->GetXaxis()->SetTitleFont(22);
 stack->GetXaxis()->SetTitleOffset(1);
 stack->GetXaxis()->SetTitle("|NPID| ");
 stack->GetXaxis()->SetBinLabel(1,"e+/e-");
 stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
 stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
 stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
 stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
 stack->GetXaxis()->SetBinLabel(6,"a");
 stack->GetXaxis()->SetBinLabel(7,"K_L0");
 stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
 stack->GetXaxis()->SetBinLabel(9,"K-/K+");
 stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
 stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>
wrote:

Many thanks for your reply.  I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-26:

#11

Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
16). I have tried to select events with MET < 10.0, but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles). I am not sure about the reason.
>
> However, the root file (selection_0.C) associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>
> wrote:
>
> Many thanks for your reply. I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
16). I have tried to select events with MET < 10.0,  but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles).  I am not sure about the reason.
>
> However, the root file (selection_0.C)  associated to the plot showed
> this:
>
>  stack->GetYaxis()->SetLabelSize(0.04);
>  stack->GetYaxis()->SetLabelOffset(0.005);
>  stack->GetYaxis()->SetTitleSize(0.06);
>  stack->GetYaxis()->SetTitleFont(22);
>  stack->GetYaxis()->SetTitleOffset(1);
>  stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
>  // X axis
>  stack->GetXaxis()->SetLabelSize(0.04);
>  stack->GetXaxis()->SetLabelOffset(0.005);
>  stack->GetXaxis()->SetTitleSize(0.06);
>  stack->GetXaxis()->SetTitleFont(22);
>  stack->GetXaxis()->SetTitleOffset(1);
>  stack->GetXaxis()->SetTitle("|NPID| ");
>  stack->GetXaxis()->SetBinLabel(1,"e+/e-");
>  stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
>  stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
>  stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
>  stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
>  stack->GetXaxis()->SetBinLabel(6,"a");
>  stack->GetXaxis()->SetBinLabel(7,"K_L0");
>  stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
>  stack->GetXaxis()->SetBinLabel(9,"K-/K+");
>  stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
>  stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>
> wrote:
>
> Many thanks for your reply.  I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-26:

#12

Hi

In reco mode there is no such thing as neutrino you can not reconstruct a neutrino. MET is define via energy conservation -> visible spectrum should balance the invisible so the MET is the neg. vector sum of visible spectra. Hence if you type `select MET < 10` ma5 will only choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI <<email address hidden><mailto:<email address hidden>>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
16). I have tried to select events with MET < 10.0, but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles). I am not sure about the reason.

However, the root file (selection_0.C) associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>
wrote:

Many thanks for your reply. I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Hi

In reco mode there is no such thing as neutrino you can not reconstruct a neutrino. MET is define via energy conservation -> visible spectrum should balance the invisible so the MET is the neg. vector sum of visible spectra. Hence if you type `select MET < 10` ma5 will only choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI <question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
16). I have tried to select events with MET < 10.0,  but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles).  I am not sure about the reason.

However, the root file (selection_0.C)  associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>
wrote:

Many thanks for your reply.  I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-29:

#13

Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj. I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
> 16). I have tried to select events with MET < 10.0, but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles). I am not sure about the reason.
>
> However, the root file (selection_0.C) associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>
> wrote:
>
> Many thanks for your reply. I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj.  I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
> 16). I have tried to select events with MET < 10.0,  but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles).  I am not sure about the reason.
>
> However, the root file (selection_0.C)  associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>
> wrote:
>
> Many thanks for your reply.  I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-29:

#14

Hi

How many events left after the cut?

—
Jack Y. Araz

On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI <<email address hidden><mailto:<email address hidden>>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj. I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

In reco mode there is no such thing as neutrino you can not reconstruct
a neutrino. MET is define via energy conservation -> visible spectrum
should balance the invisible so the MET is the neg. vector sum of
visible spectra. Hence if you type `select MET < 10` ma5 will only
choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
<<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
16). I have tried to select events with MET < 10.0, but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles). I am not sure about the reason.

However, the root file (selection_0.C) associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>>
wrote:

Many thanks for your reply. I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Hi

How many events left after the cut?

—
Jack Y. Araz

On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI <question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj.  I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

In reco mode there is no such thing as neutrino you can not reconstruct
a neutrino. MET is define via energy conservation -> visible spectrum
should balance the invisible so the MET is the neg. vector sum of
visible spectra. Hence if you type `select MET < 10` ma5 will only
choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
16). I have tried to select events with MET < 10.0,  but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles).  I am not sure about the reason.

However, the root file (selection_0.C)  associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>>
wrote:

Many thanks for your reply.  I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-29:

#15

Hi Jack,

Many thanks for your reply!
After the cut : MET < 10.0
events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
efficiency(K/ (K+R)): 0.926+/-0.119
Ps: the initial number of the events is 10000.

Thank you very much for any help or suggestion!

Best
Souad

Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> How many events left after the cut?
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Many thanks for your clarification. Apologies, this problem has dragged on
> for so long.
> Select MET< 10GeV is working fine, and I am only getting events
> with missing ET less than 10GeV. But I am still getting empty plots, at
> hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
> mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
> the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj. I got a peak
> at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .
>
> I couldn't identify the reason behind empty plots.
>
> Many thanks for any suggestions or help to alleviate the confusion.
>
> Best regards
> Souad
>
> Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
> 16). I have tried to select events with MET < 10.0, but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles). I am not sure about the reason.
>
> However, the root file (selection_0.C) associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>>
> wrote:
>
> Many thanks for your reply. I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Hi Jack,

Many thanks for your reply!
After the cut : MET < 10.0
 events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
efficiency(K/ (K+R)): 0.926+/-0.119
 Ps: the initial number of the events is 10000.

Thank you very much for any help or suggestion!

Best
Souad

Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> How many events left after the cut?
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Many thanks for your clarification. Apologies, this problem has dragged on
> for so long.
> Select MET< 10GeV is working fine, and I am only getting events
> with missing ET less than 10GeV. But I am still getting empty plots, at
> hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
> mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
> the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj.  I got a peak
> at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .
>
> I couldn't identify the reason behind empty plots.
>
> Many thanks for any suggestions or help to alleviate the confusion.
>
> Best regards
> Souad
>
> Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
> 16). I have tried to select events with MET < 10.0,  but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>  Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles).  I am not sure about the reason.
>
> However, the root file (selection_0.C)  associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>  Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>>
> wrote:
>
> Many thanks for your reply.  I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-29:

#16

Hi

This value is scaled with cross section and luminosity you need to check how many MC events left in your sample after the cut if that value is too low that will explain you why your plots are empty.

—
Jack Y. Araz

On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI <<email address hidden><mailto:<email address hidden>>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Many thanks for your reply!
After the cut : MET < 10.0
events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
efficiency(K/ (K+R)): 0.926+/-0.119
Ps: the initial number of the events is 10000.

Thank you very much for any help or suggestion!

Best
Souad

Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

How many events left after the cut?

—
Jack Y. Araz

On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
<<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj. I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

In reco mode there is no such thing as neutrino you can not reconstruct
a neutrino. MET is define via energy conservation -> visible spectrum
should balance the invisible so the MET is the neg. vector sum of
visible spectra. Hence if you type `select MET < 10` ma5 will only
choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
16). I have tried to select events with MET < 10.0, but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles). I am not sure about the reason.

However, the root file (selection_0.C) associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:
<email address hidden>><mailto:
<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>>>
wrote:

Many thanks for your reply. I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Hi

This value is scaled with cross section and luminosity you need to check how many MC events left in your sample after the cut if that value is too low that will explain you why your plots are empty.

—
Jack Y. Araz

On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI <question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Many thanks for your reply!
After the cut : MET < 10.0
events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
efficiency(K/ (K+R)): 0.926+/-0.119
Ps: the initial number of the events is 10000.

Thank you very much for any help or suggestion!

Best
Souad

Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

How many events left after the cut?

—
Jack Y. Araz

On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj.  I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

In reco mode there is no such thing as neutrino you can not reconstruct
a neutrino. MET is define via energy conservation -> visible spectrum
should balance the invisible so the MET is the neg. vector sum of
visible spectra. Hence if you type `select MET < 10` ma5 will only
choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
16). I have tried to select events with MET < 10.0,  but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles).  I am not sure about the reason.

However, the root file (selection_0.C)  associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>>>
wrote:

Many thanks for your reply.  I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-29:

#17

Hi Jack,

I have checked Cutflows/myregion_SAF, and the number of events kept is
9261. I don't think this number is too small?
Please correct me if I am wrong.
<Counter>
"1_MET < 10.0" # 1st cut
9261 0 # nentries
4.471067e+00 0.000000e+00 # sum of weights
2.158562e-03 0.000000e+00 # sum of weights^2
</Counter>

Many thanks for your help!

Best
Souad

Le lun. 29 nov. 2021 à 12:15, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> This value is scaled with cross section and luminosity you need to check
> how many MC events left in your sample after the cut if that value is
> too low that will explain you why your plots are empty.
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Many thanks for your reply!
> After the cut : MET < 10.0
> events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
> efficiency(K/ (K+R)): 0.926+/-0.119
> Ps: the initial number of the events is 10000.
>
> Thank you very much for any help or suggestion!
>
> Best
> Souad
>
> Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> How many events left after the cut?
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Many thanks for your clarification. Apologies, this problem has dragged on
> for so long.
> Select MET< 10GeV is working fine, and I am only getting events
> with missing ET less than 10GeV. But I am still getting empty plots, at
> hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
> mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
> the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj. I got a peak
> at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .
>
> I couldn't identify the reason behind empty plots.
>
> Many thanks for any suggestions or help to alleviate the confusion.
>
> Best regards
> Souad
>
> Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
> 16). I have tried to select events with MET < 10.0, but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles). I am not sure about the reason.
>
> However, the root file (selection_0.C) associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>>>
> wrote:
>
> Many thanks for your reply. I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=15
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Hi Jack,

I have checked Cutflows/myregion_SAF, and the number of events kept is
9261.  I don't think this number is too small?
  Please correct me if I am wrong.
 <Counter>
"1_MET < 10.0"                  # 1st cut
9261            0               # nentries
4.471067e+00    0.000000e+00    # sum of weights
2.158562e-03    0.000000e+00    # sum of weights^2
</Counter>

Many thanks for your help!

Best
Souad

Le lun. 29 nov. 2021 à 12:15, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> This value is scaled with cross section and luminosity you need to check
> how many MC events left in your sample after the cut if that value is
> too low that will explain you why your plots are empty.
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Many thanks for your reply!
> After the cut : MET < 10.0
> events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
> efficiency(K/ (K+R)): 0.926+/-0.119
> Ps: the initial number of the events is 10000.
>
> Thank you very much for any help or suggestion!
>
> Best
> Souad
>
> Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> How many events left after the cut?
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Many thanks for your clarification. Apologies, this problem has dragged on
> for so long.
> Select MET< 10GeV is working fine, and I am only getting events
> with missing ET less than 10GeV. But I am still getting empty plots, at
> hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
> mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
> the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj.  I got a peak
> at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .
>
> I couldn't identify the reason behind empty plots.
>
> Many thanks for any suggestions or help to alleviate the confusion.
>
> Best regards
> Souad
>
> Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>  Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
> 16). I have tried to select events with MET < 10.0,  but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>  Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles).  I am not sure about the reason.
>
> However, the root file (selection_0.C)  associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>>>
> wrote:
>
> Many thanks for your reply.  I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=15
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-29:

#18

No it is not, can you also check histos.saf to see if they are filled. If not you might be plotting something wrong or maybe using the wrong syntax

—
Jack Y. Araz

On Nov 29, 2021, at 12:30 PM, SOUAD SEMLALI <<email address hidden><mailto:<email address hidden>>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

I have checked Cutflows/myregion_SAF, and the number of events kept is
9261. I don't think this number is too small?
Please correct me if I am wrong.
<Counter>
"1_MET < 10.0" # 1st cut
9261 0 # nentries
4.471067e+00 0.000000e+00 # sum of weights
2.158562e-03 0.000000e+00 # sum of weights^2
</Counter>

Many thanks for your help!

Best
Souad

Le lun. 29 nov. 2021 à 12:15, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

This value is scaled with cross section and luminosity you need to check
how many MC events left in your sample after the cut if that value is
too low that will explain you why your plots are empty.

—
Jack Y. Araz

On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI
<<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Many thanks for your reply!
After the cut : MET < 10.0
events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
efficiency(K/ (K+R)): 0.926+/-0.119
Ps: the initial number of the events is 10000.

Thank you very much for any help or suggestion!

Best
Souad

Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

How many events left after the cut?

—
Jack Y. Araz

On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj. I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

In reco mode there is no such thing as neutrino you can not reconstruct
a neutrino. MET is define via energy conservation -> visible spectrum
should balance the invisible so the MET is the neg. vector sum of
visible spectra. Hence if you type `select MET < 10` ma5 will only
choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
16). I have tried to select events with MET < 10.0, but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:
<email address hidden>><mailto:
<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles). I am not sure about the reason.

However, the root file (selection_0.C) associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:
<email address hidden>><mailto:
<email address hidden><mailto:
<email address hidden><mailto:<email address hidden>>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>><mailto:
<email address hidden><mailto:
<email address hidden>><mailto:
<email address hidden><mailto:
<email address hidden><mailto:
<email address hidden>>>><mailto:
<email address hidden><mailto:
<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>><mailto:
<email address hidden><mailto:<email address hidden>><mailto:
<email address hidden><mailto:<email address hidden>>>>>>>>
wrote:

Many thanks for your reply. I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=15

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

No it is not, can you also check histos.saf to see if they are filled. If not you might be plotting something wrong or maybe using the wrong syntax

—
Jack Y. Araz

On Nov 29, 2021, at 12:30 PM, SOUAD SEMLALI <question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

I have checked Cutflows/myregion_SAF, and the number of events kept is
9261.  I don't think this number is too small?
 Please correct me if I am wrong.
<Counter>
"1_MET < 10.0"                  # 1st cut
9261            0               # nentries
4.471067e+00    0.000000e+00    # sum of weights
2.158562e-03    0.000000e+00    # sum of weights^2
</Counter>

Many thanks for your help!

Best
Souad

Le lun. 29 nov. 2021 à 12:15, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

This value is scaled with cross section and luminosity you need to check
how many MC events left in your sample after the cut if that value is
too low that will explain you why your plots are empty.

—
Jack Y. Araz

On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Many thanks for your reply!
After the cut : MET < 10.0
events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
efficiency(K/ (K+R)): 0.926+/-0.119
Ps: the initial number of the events is 10000.

Thank you very much for any help or suggestion!

Best
Souad

Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

How many events left after the cut?

—
Jack Y. Araz

On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Dear Jack,

Many thanks for your clarification. Apologies, this problem has dragged on
for so long.
Select MET< 10GeV is working fine, and I am only getting events
with missing ET less than 10GeV. But I am still getting empty plots, at
hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj.  I got a peak
at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .

I couldn't identify the reason behind empty plots.

Many thanks for any suggestions or help to alleviate the confusion.

Best regards
Souad

Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

In reco mode there is no such thing as neutrino you can not reconstruct
a neutrino. MET is define via energy conservation -> visible spectrum
should balance the invisible so the MET is the neg. vector sum of
visible spectra. Hence if you type `select MET < 10` ma5 will only
choose the events with met < 10 GeV etc etc.

—
Jack Y. Araz

On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your help! I have checked the hepmc file. The MET
is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
16). I have tried to select events with MET < 10.0,  but I am getting an
error.

define nu = 12 -12 14 -14 16 -16
reject (nu) MET > 50.0

Is there another way to remove these neutrinos after hadronization?
Many thanks for any help you can provide.

Best regards
Souad

Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

Portion of the code you sent just showing you x-axis labels the
information is taken from histos.saf file within
`<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
number>/Histograms/histos.saf` if the plot that you are looking for is
empty there than there is a problem. but if it has information not
translated to root there might be some interference between root and ma5
so you might want to change the plotting renderer (explained how in FAQ
for normal mode for a different example). It might also be the case that
for some reason ma5 can not read the xsec information from your hepmc
file (if you are giving root file there is no xsec information so you
need to feed it by hand) so I suggest setting sample xsec as well to be
on the safe side. For how to do these you can refer to normal mode
reference card in your local madanalysis/doc folder.

—
Jack Y. Araz

On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>>>
wrote:

[EXTERNAL EMAIL]

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

Status: Answered => Open

SOUAD SEMLALI is still having a problem:
Hi Jack,

Thank you very much for your response! I tried to plot NPID (plot NPID),
but I am getting an empty plot (no labels, no information about the
particles).  I am not sure about the reason.

However, the root file (selection_0.C)  associated to the plot showed
this:

stack->GetYaxis()->SetLabelSize(0.04);
stack->GetYaxis()->SetLabelOffset(0.005);
stack->GetYaxis()->SetTitleSize(0.06);
stack->GetYaxis()->SetTitleFont(22);
stack->GetYaxis()->SetTitleOffset(1);
stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");

// X axis
stack->GetXaxis()->SetLabelSize(0.04);
stack->GetXaxis()->SetLabelOffset(0.005);
stack->GetXaxis()->SetTitleSize(0.06);
stack->GetXaxis()->SetTitleFont(22);
stack->GetXaxis()->SetTitleOffset(1);
stack->GetXaxis()->SetTitle("|NPID| ");
stack->GetXaxis()->SetBinLabel(1,"e+/e-");
stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
stack->GetXaxis()->SetBinLabel(6,"a");
stack->GetXaxis()->SetBinLabel(7,"K_L0");
stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
stack->GetXaxis()->SetBinLabel(9,"K-/K+");
stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");

Many thanks for any suggestions or help.

Best regards
Souad

Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>>> a écrit :

Your question #699618 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/699618

Status: Open => Answered

Jack Y. Araz proposed the following answer:
Hi

If hadron level analysis is giving you a large MET you will get the same
thing from reconstructed analysis. I suggest plotting NPID to see which
particles are in your final state event and/or directly read the HEPMC
file and try to understand what is producing missing energy.

For truth level reco analysis; please see the description given in
arXiv:2006.09387.

—
Jack Y. Araz

On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
<question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net>>><mailto:
question699618@answers.launchpad.net<mailto:
question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net><mailto:
question699618@answers.launchpad.net<mailto:question699618@answers.launchpad.net>>>>>>>
wrote:

Many thanks for your reply.  I have analyzed the hadronic output through

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

--
If this answers your question, please go to the following page to let us
know that it is solved:

https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=15

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madanalysis5/+question/699618

You received this question notification because you asked the question.

--
Souad SEMLALI

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

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-29:

#19

Hi Jack,

I have checked histos.saf. The values are set to 0. This is for the
invariant mass distribution of b b~.
( at parton level, it is not ste to 0. I got mbb~ ~ mh)
This is from histos.saf:
<Histo>
  <Description>
    "3_M"
    # nbins xmin xmax
      40 0.000000e+00 2.000000e+02
    # Defined regions
      myregion # Region nr. 1
  </Description>
  <Statistics>
      0 0 # nevents
      0.000000e+00 0.000000e+00 # sum of event-weights over events
      0 0 # nentries
      0.000000e+00 0.000000e+00 # sum of event-weights over entries
      0.000000e+00 0.000000e+00 # sum weights^2
      0.000000e+00 0.000000e+00 # sum value*weight
      0.000000e+00 0.000000e+00 # sum value^2*weight
  </Statistics>
  <Data>
      0.000000e+00 0.000000e+00 # underflow
      0.000000e+00 0.000000e+00 # bin 1 / 40
      0.000000e+00 0.000000e+00 # bin 2 / 40
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00 # bin 39 / 40
      0.000000e+00 0.000000e+00 # bin 40 / 40
      0.000000e+00 0.000000e+00 # overflow
  </Data>
</Histo>

Many thanks for any help!

Best regards

Le lun. 29 nov. 2021 à 12:41, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> No it is not, can you also check histos.saf to see if they are filled.
> If not you might be plotting something wrong or maybe using the wrong
> syntax
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 12:30 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> I have checked Cutflows/myregion_SAF, and the number of events kept is
> 9261. I don't think this number is too small?
> Please correct me if I am wrong.
> <Counter>
> "1_MET < 10.0" # 1st cut
> 9261 0 # nentries
> 4.471067e+00 0.000000e+00 # sum of weights
> 2.158562e-03 0.000000e+00 # sum of weights^2
> </Counter>
>
> Many thanks for your help!
>
> Best
> Souad
>
>
> Le lun. 29 nov. 2021 à 12:15, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> This value is scaled with cross section and luminosity you need to check
> how many MC events left in your sample after the cut if that value is
> too low that will explain you why your plots are empty.
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Many thanks for your reply!
> After the cut : MET < 10.0
> events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
> efficiency(K/ (K+R)): 0.926+/-0.119
> Ps: the initial number of the events is 10000.
>
> Thank you very much for any help or suggestion!
>
> Best
> Souad
>
> Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> How many events left after the cut?
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Many thanks for your clarification. Apologies, this problem has dragged on
> for so long.
> Select MET< 10GeV is working fine, and I am only getting events
> with missing ET less than 10GeV. But I am still getting empty plots, at
> hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
> mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
> the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj. I got a peak
> at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .
>
> I couldn't identify the reason behind empty plots.
>
> Many thanks for any suggestions or help to alleviate the confusion.
>
> Best regards
> Souad
>
> Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12, -16, -14, 12, 14,
> 16). I have tried to select events with MET < 10.0, but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles). I am not sure about the reason.
>
> However, the root file (selection_0.C) associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <<email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden><mailto:
> <email address hidden>>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>><mailto:
> <email address hidden><mailto:
> <email address hidden>><mailto:
> <email address hidden><mailto:
> <email address hidden>>>>>>>>
> wrote:
>
> Many thanks for your reply. I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=15
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=17
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Hi Jack,

I have checked histos.saf. The values are set to 0.  This is for the
invariant mass distribution of b b~.
( at parton level, it is not ste to 0.  I got mbb~ ~ mh)
This is from histos.saf:
<Histo>
  <Description>
    "3_M"
    # nbins   xmin           xmax
      40      0.000000e+00   2.000000e+02
    # Defined regions
      myregion    # Region nr. 1
  </Description>
  <Statistics>
      0              0               # nevents
      0.000000e+00   0.000000e+00    # sum of event-weights over events
      0              0               # nentries
      0.000000e+00   0.000000e+00    # sum of event-weights over entries
      0.000000e+00   0.000000e+00    # sum weights^2
      0.000000e+00   0.000000e+00    # sum value*weight
      0.000000e+00   0.000000e+00    # sum value^2*weight
  </Statistics>
  <Data>
      0.000000e+00   0.000000e+00    # underflow
      0.000000e+00   0.000000e+00    # bin 1 / 40
      0.000000e+00   0.000000e+00    # bin 2 / 40
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00
      0.000000e+00   0.000000e+00    # bin 39 / 40
      0.000000e+00   0.000000e+00    # bin 40 / 40
      0.000000e+00   0.000000e+00    # overflow
  </Data>
</Histo>

Many thanks for any help!

Best regards

Le lun. 29 nov. 2021 à 12:41, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> No it is not, can you also check histos.saf to see if they are filled.
> If not you might be plotting something wrong or maybe using the wrong
> syntax
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 12:30 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> I have checked Cutflows/myregion_SAF, and the number of events kept is
> 9261.  I don't think this number is too small?
>  Please correct me if I am wrong.
> <Counter>
> "1_MET < 10.0"                  # 1st cut
> 9261            0               # nentries
> 4.471067e+00    0.000000e+00    # sum of weights
> 2.158562e-03    0.000000e+00    # sum of weights^2
> </Counter>
>
> Many thanks for your help!
>
> Best
> Souad
>
>
> Le lun. 29 nov. 2021 à 12:15, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>    Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> This value is scaled with cross section and luminosity you need to check
> how many MC events left in your sample after the cut if that value is
> too low that will explain you why your plots are empty.
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 12:11 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Many thanks for your reply!
> After the cut : MET < 10.0
> events kept (K): 4.471+/-0.575, rejected events(R): 0.3758+/- 0.575,
> efficiency(K/ (K+R)): 0.926+/-0.119
> Ps: the initial number of the events is 10000.
>
> Thank you very much for any help or suggestion!
>
> Best
> Souad
>
> Le lun. 29 nov. 2021 à 11:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>   Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> How many events left after the cut?
>
> —
> Jack Y. Araz
>
> On Nov 29, 2021, at 11:50 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>  Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Many thanks for your clarification. Apologies, this problem has dragged on
> for so long.
> Select MET< 10GeV is working fine, and I am only getting events
> with missing ET less than 10GeV. But I am still getting empty plots, at
> hadron level (./bin/ma5 -H script), for invariant mass distributions (mbb~,
> mjj, m{bbjj}, m_{bbjjmu+mu-}), which wasn't the case at Parton level (for
> the signal Hsm -> z*z* h with h > bb~, z> mu+ mu- and z > jj.  I got a peak
> at mbb~mh, peak at m{bb~jjmu+ mu-}~mHsm at Parton level) .
>
> I couldn't identify the reason behind empty plots.
>
> Many thanks for any suggestions or help to alleviate the confusion.
>
> Best regards
> Souad
>
> Le ven. 26 nov. 2021 à 18:35, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>  Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> In reco mode there is no such thing as neutrino you can not reconstruct
> a neutrino. MET is define via energy conservation -> visible spectrum
> should balance the invisible so the MET is the neg. vector sum of
> visible spectra. Hence if you type `select MET < 10` ma5 will only
> choose the events with met < 10 GeV etc etc.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 6:25 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your help! I have checked the hepmc file. The MET
> is indeed due to the presence of the neutrinos (-12,  -16,  -14, 12, 14,
> 16). I have tried to select events with MET < 10.0,  but I am getting an
> error.
>
> define nu = 12 -12 14 -14 16 -16
> reject (nu) MET > 50.0
>
> Is there another way to remove these neutrinos after hadronization?
> Many thanks for any help you can provide.
>
> Best regards
> Souad
>
>
> Le ven. 26 nov. 2021 à 11:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Portion of the code you sent just showing you x-axis labels the
> information is taken from histos.saf file within
> `<my_ana>/Outputs/SAF/<my_smp>/MadAnalysisJob_<latest run
> number>/Histograms/histos.saf` if the plot that you are looking for is
> empty there than there is a problem. but if it has information not
> translated to root there might be some interference between root and ma5
> so you might want to change the plotting renderer (explained how in FAQ
> for normal mode for a different example). It might also be the case that
> for some reason ma5 can not read the xsec information from your hepmc
> file (if you are giving root file there is no xsec information so you
> need to feed it by hand) so I suggest setting sample xsec as well to be
> on the safe side. For how to do these you can refer to normal mode
> reference card in your local madanalysis/doc folder.
>
> —
> Jack Y. Araz
>
> On Nov 26, 2021, at 11:41 AM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>>>
> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Jack,
>
> Thank you very much for your response! I tried to plot NPID (plot NPID),
> but I am getting an empty plot (no labels, no information about the
> particles).  I am not sure about the reason.
>
> However, the root file (selection_0.C)  associated to the plot showed
> this:
>
> stack->GetYaxis()->SetLabelSize(0.04);
> stack->GetYaxis()->SetLabelOffset(0.005);
> stack->GetYaxis()->SetTitleSize(0.06);
> stack->GetYaxis()->SetTitleFont(22);
> stack->GetYaxis()->SetTitleOffset(1);
> stack->GetYaxis()->SetTitle("N. of particles ( scaled to one )");
>
> // X axis
> stack->GetXaxis()->SetLabelSize(0.04);
> stack->GetXaxis()->SetLabelOffset(0.005);
> stack->GetXaxis()->SetTitleSize(0.06);
> stack->GetXaxis()->SetTitleFont(22);
> stack->GetXaxis()->SetTitleOffset(1);
> stack->GetXaxis()->SetTitle("|NPID| ");
> stack->GetXaxis()->SetBinLabel(1,"e+/e-");
> stack->GetXaxis()->SetBinLabel(2,"ve~/ve");
> stack->GetXaxis()->SetBinLabel(3,"mu+/mu-");
> stack->GetXaxis()->SetBinLabel(4,"vm~/vm");
> stack->GetXaxis()->SetBinLabel(5,"vt~/vt");
> stack->GetXaxis()->SetBinLabel(6,"a");
> stack->GetXaxis()->SetBinLabel(7,"K_L0");
> stack->GetXaxis()->SetBinLabel(8,"pi-/pi+");
> stack->GetXaxis()->SetBinLabel(9,"K-/K+");
> stack->GetXaxis()->SetBinLabel(10,"nbar0/n0");
> stack->GetXaxis()->SetBinLabel(11,"pbar-/p+");
>
> Many thanks for any suggestions or help.
>
> Best regards
> Souad
>
> Le jeu. 25 nov. 2021 à 17:55, Jack Y. Araz <
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>>> a écrit :
>
> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> If hadron level analysis is giving you a large MET you will get the same
> thing from reconstructed analysis. I suggest plotting NPID to see which
> particles are in your final state event and/or directly read the HEPMC
> file and try to understand what is producing missing energy.
>
> For truth level reco analysis; please see the description given in
> arXiv:2006.09387.
>
> —
> Jack Y. Araz
>
> On Nov 25, 2021, at 5:45 PM, SOUAD SEMLALI
> <question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net><mailto:
> question699618@answers.launchpad.net<mailto:
> question699618@answers.launchpad.net>>>>>>>
> wrote:
>
> Many thanks for your reply.  I have analyzed the hadronic output through
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=7
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=9
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=11
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=13
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=15
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=17
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-29:

#20

Hi Jack,

Following the previous message, I have checked the histos.saf. The values are set to 0.

This is from histos.saf:
<Histo>
  <Description>
    "3_M"
    # nbins xmin xmax
      40 0.000000e+00 2.000000e+02
    # Defined regions
      myregion # Region nr. 1
  </Description>
  <Statistics>
      0 0 # nevents
      0.000000e+00 0.000000e+00 # sum of event-weights over events
      0 0 # nentries
      0.000000e+00 0.000000e+00 # sum of event-weights over entries
      0.000000e+00 0.000000e+00 # sum weights^2
      0.000000e+00 0.000000e+00 # sum value*weight
      0.000000e+00 0.000000e+00 # sum value^2*weight
  </Statistics>
  <Data>
      0.000000e+00 0.000000e+00 # underflow
      0.000000e+00 0.000000e+00 # bin 1 / 40
      0.000000e+00 0.000000e+00 # bin 2 / 40
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00 # bin 39 / 40
      0.000000e+00 0.000000e+00 # bin 40 / 40
      0.000000e+00 0.000000e+00 # overflow
  </Data>
</Histo>

I am using the following syntax:
define q = 5 -5
plot M(q[1] q[2]) 40 0 200 [logY ]

Am I missing something?

Many thanks for any help!

Best
Souad

Hi Jack,

Following the previous message, I have checked the histos.saf. The values are set to 0.

This is from histos.saf:
<Histo>
  <Description>
    "3_M"
    # nbins xmin xmax
      40 0.000000e+00 2.000000e+02
    # Defined regions
      myregion # Region nr. 1
  </Description>
  <Statistics>
      0 0 # nevents
      0.000000e+00 0.000000e+00 # sum of event-weights over events
      0 0 # nentries
      0.000000e+00 0.000000e+00 # sum of event-weights over entries
      0.000000e+00 0.000000e+00 # sum weights^2
      0.000000e+00 0.000000e+00 # sum value*weight
      0.000000e+00 0.000000e+00 # sum value^2*weight
  </Statistics>
  <Data>
      0.000000e+00 0.000000e+00 # underflow
      0.000000e+00 0.000000e+00 # bin 1 / 40
      0.000000e+00 0.000000e+00 # bin 2 / 40
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00
      0.000000e+00 0.000000e+00 # bin 39 / 40
      0.000000e+00 0.000000e+00 # bin 40 / 40
      0.000000e+00 0.000000e+00 # overflow
  </Data>
</Histo>

I am using the following syntax:
 define q = 5 -5
 plot M(q[1] q[2]) 40 0  200 [logY ]

Am I missing something?

Many thanks for any help!

Best
Souad

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-29:

#21

you can not use PDGID in reco mode this is not defined. reco mode consists of reconstructed objects not particles. please try the following

$ ./bin/ma5 -R
ma5> set main.fastsim.package = fastjet
ma5>import <your HEPMC sample>
ma5>plot M(b[1] b[2])
ma5>submit

—
Jack Y. Araz

> On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> SOUAD SEMLALI gave more information on the question:
> Hi Jack,
>
> Following the previous message, I have checked the histos.saf. The
> values are set to 0.
>
> This is from histos.saf:
> <Histo>
> <Description>
> "3_M"
> # nbins xmin xmax
> 40 0.000000e+00 2.000000e+02
> # Defined regions
> myregion # Region nr. 1
> </Description>
> <Statistics>
> 0 0 # nevents
> 0.000000e+00 0.000000e+00 # sum of event-weights over events
> 0 0 # nentries
> 0.000000e+00 0.000000e+00 # sum of event-weights over entries
> 0.000000e+00 0.000000e+00 # sum weights^2
> 0.000000e+00 0.000000e+00 # sum value*weight
> 0.000000e+00 0.000000e+00 # sum value^2*weight
> </Statistics>
> <Data>
> 0.000000e+00 0.000000e+00 # underflow
> 0.000000e+00 0.000000e+00 # bin 1 / 40
> 0.000000e+00 0.000000e+00 # bin 2 / 40
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00
> 0.000000e+00 0.000000e+00 # bin 39 / 40
> 0.000000e+00 0.000000e+00 # bin 40 / 40
> 0.000000e+00 0.000000e+00 # overflow
> </Data>
> </Histo>
>
>
> I am using the following syntax:
> define q = 5 -5
> plot M(q[1] q[2]) 40 0 200 [logY ]
>
> Am I missing something?
>
> Many thanks for any help!
>
> Best
> Souad
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

you can not use PDGID in reco mode this is not defined. reco mode consists of reconstructed objects not particles. please try the following

$ ./bin/ma5 -R
ma5> set main.fastsim.package = fastjet
ma5>import <your HEPMC sample>
ma5>plot M(b[1] b[2])
ma5>submit

—
Jack Y. Araz

> On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <question699618@answers.launchpad.net> wrote:
> 
> [EXTERNAL EMAIL]
> 
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
> 
> SOUAD SEMLALI gave more information on the question:
> Hi Jack,
> 
> Following the previous message, I have checked the histos.saf. The
> values are set to 0.
> 
> This is from histos.saf:
> <Histo>
>  <Description>
>    "3_M"
>    # nbins xmin xmax
>      40 0.000000e+00 2.000000e+02
>    # Defined regions
>      myregion # Region nr. 1
>  </Description>
>  <Statistics>
>      0 0 # nevents
>      0.000000e+00 0.000000e+00 # sum of event-weights over events
>      0 0 # nentries
>      0.000000e+00 0.000000e+00 # sum of event-weights over entries
>      0.000000e+00 0.000000e+00 # sum weights^2
>      0.000000e+00 0.000000e+00 # sum value*weight
>      0.000000e+00 0.000000e+00 # sum value^2*weight
>  </Statistics>
>  <Data>
>      0.000000e+00 0.000000e+00 # underflow
>      0.000000e+00 0.000000e+00 # bin 1 / 40
>      0.000000e+00 0.000000e+00 # bin 2 / 40
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00
>      0.000000e+00 0.000000e+00 # bin 39 / 40
>      0.000000e+00 0.000000e+00 # bin 40 / 40
>      0.000000e+00 0.000000e+00 # overflow
>  </Data>
> </Histo>
> 
> 
> I am using the following syntax:
> define q = 5 -5
> plot M(q[1] q[2]) 40 0  200 [logY ]
> 
> Am I missing something?
> 
> Many thanks for any help!
> 
> Best
> Souad
> 
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-29:

#22

Dear Jack,

Thank you for your reply and your help all this time. The plots are
produced, but I am afraid this brings me back to my starting point. The
distributions at reco-mode are completely different from the ones at the
partonic level and confusing. I thought that this is due to the presence of
large MET, but suppressing MET doesn't solve the problem. As I mentioned
before, I am working on signal Hsm > z* z* h, with z* > mu+ mu-, z* > j j
and h > b b~. But The invariant mass distributions are completely different
from this; distributions with several peaks, and m_{bb~} != m_h and m_{mu+
mu- jj bb~}!= m_{Hsm}.
I would be grateful for any suggestions or help!

Best
Souad

Le lun. 29 nov. 2021 à 14:11, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> you can not use PDGID in reco mode this is not defined. reco mode
> consists of reconstructed objects not particles. please try the
> following
>
> $ ./bin/ma5 -R
> ma5> set main.fastsim.package = fastjet
> ma5>import <your HEPMC sample>
> ma5>plot M(b[1] b[2])
> ma5>submit
>
> —
> Jack Y. Araz
>
> > On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <
> <email address hidden>> wrote:
> >
> > [EXTERNAL EMAIL]
> >
> > Question #699618 on MadAnalysis 5 changed:
> > https://answers.launchpad.net/madanalysis5/+question/699618
> >
> > SOUAD SEMLALI gave more information on the question:
> > Hi Jack,
> >
> > Following the previous message, I have checked the histos.saf. The
> > values are set to 0.
> >
> > This is from histos.saf:
> > <Histo>
> > <Description>
> > "3_M"
> > # nbins xmin xmax
> > 40 0.000000e+00 2.000000e+02
> > # Defined regions
> > myregion # Region nr. 1
> > </Description>
> > <Statistics>
> > 0 0 # nevents
> > 0.000000e+00 0.000000e+00 # sum of event-weights over events
> > 0 0 # nentries
> > 0.000000e+00 0.000000e+00 # sum of event-weights over entries
> > 0.000000e+00 0.000000e+00 # sum weights^2
> > 0.000000e+00 0.000000e+00 # sum value*weight
> > 0.000000e+00 0.000000e+00 # sum value^2*weight
> > </Statistics>
> > <Data>
> > 0.000000e+00 0.000000e+00 # underflow
> > 0.000000e+00 0.000000e+00 # bin 1 / 40
> > 0.000000e+00 0.000000e+00 # bin 2 / 40
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00
> > 0.000000e+00 0.000000e+00 # bin 39 / 40
> > 0.000000e+00 0.000000e+00 # bin 40 / 40
> > 0.000000e+00 0.000000e+00 # overflow
> > </Data>
> > </Histo>
> >
> >
> > I am using the following syntax:
> > define q = 5 -5
> > plot M(q[1] q[2]) 40 0 200 [logY ]
> >
> > Am I missing something?
> >
> > Many thanks for any help!
> >
> > Best
> > Souad
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=20
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Dear Jack,

Thank you for your reply and your help all this time. The plots are
produced, but I am afraid this brings me back to my starting point. The
distributions at reco-mode are completely different from the ones at the
partonic level and confusing. I thought that this is due to the presence of
large MET, but suppressing MET doesn't solve the problem. As I mentioned
before, I am working on signal Hsm > z* z* h, with z* > mu+ mu-, z* > j j
and h > b b~. But The invariant mass distributions are completely different
from this; distributions with several peaks, and m_{bb~} != m_h and m_{mu+
mu- jj bb~}!= m_{Hsm}.
I would be grateful for any suggestions or help!

Best
Souad

Le lun. 29 nov. 2021 à 14:11, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> you can not use PDGID in reco mode this is not defined. reco mode
> consists of reconstructed objects not particles. please try the
> following
>
> $ ./bin/ma5 -R
> ma5> set main.fastsim.package = fastjet
> ma5>import <your HEPMC sample>
> ma5>plot M(b[1] b[2])
> ma5>submit
>
> —
> Jack Y. Araz
>
> > On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <
> question699618@answers.launchpad.net> wrote:
> >
> > [EXTERNAL EMAIL]
> >
> > Question #699618 on MadAnalysis 5 changed:
> > https://answers.launchpad.net/madanalysis5/+question/699618
> >
> > SOUAD SEMLALI gave more information on the question:
> > Hi Jack,
> >
> > Following the previous message, I have checked the histos.saf. The
> > values are set to 0.
> >
> > This is from histos.saf:
> > <Histo>
> >  <Description>
> >    "3_M"
> >    # nbins xmin xmax
> >      40 0.000000e+00 2.000000e+02
> >    # Defined regions
> >      myregion # Region nr. 1
> >  </Description>
> >  <Statistics>
> >      0 0 # nevents
> >      0.000000e+00 0.000000e+00 # sum of event-weights over events
> >      0 0 # nentries
> >      0.000000e+00 0.000000e+00 # sum of event-weights over entries
> >      0.000000e+00 0.000000e+00 # sum weights^2
> >      0.000000e+00 0.000000e+00 # sum value*weight
> >      0.000000e+00 0.000000e+00 # sum value^2*weight
> >  </Statistics>
> >  <Data>
> >      0.000000e+00 0.000000e+00 # underflow
> >      0.000000e+00 0.000000e+00 # bin 1 / 40
> >      0.000000e+00 0.000000e+00 # bin 2 / 40
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00
> >      0.000000e+00 0.000000e+00 # bin 39 / 40
> >      0.000000e+00 0.000000e+00 # bin 40 / 40
> >      0.000000e+00 0.000000e+00 # overflow
> >  </Data>
> > </Histo>
> >
> >
> > I am using the following syntax:
> > define q = 5 -5
> > plot M(q[1] q[2]) 40 0  200 [logY ]
> >
> > Am I missing something?
> >
> > Many thanks for any help!
> >
> > Best
> > Souad
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=20
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

Jack Y. Araz (jackaraz) said on 2021-11-29:

#23

Hi

I’m gonna leave the physics part to you, can’t do your project for you. I’m sure you will figure in no time.

cheers

—
Jack Y. Araz

> On Nov 29, 2021, at 3:45 PM, SOUAD SEMLALI <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
>
> Thank you for your reply and your help all this time. The plots are
> produced, but I am afraid this brings me back to my starting point. The
> distributions at reco-mode are completely different from the ones at the
> partonic level and confusing. I thought that this is due to the presence of
> large MET, but suppressing MET doesn't solve the problem. As I mentioned
> before, I am working on signal Hsm > z* z* h, with z* > mu+ mu-, z* > j j
> and h > b b~. But The invariant mass distributions are completely different
> from this; distributions with several peaks, and m_{bb~} != m_h and m_{mu+
> mu- jj bb~}!= m_{Hsm}.
> I would be grateful for any suggestions or help!
>
> Best
> Souad
>
>
> Le lun. 29 nov. 2021 à 14:11, Jack Y. Araz <
> <email address hidden>> a écrit :
>
>> Your question #699618 on MadAnalysis 5 changed:
>> https://answers.launchpad.net/madanalysis5/+question/699618
>>
>> Status: Open => Answered
>>
>> Jack Y. Araz proposed the following answer:
>> you can not use PDGID in reco mode this is not defined. reco mode
>> consists of reconstructed objects not particles. please try the
>> following
>>
>> $ ./bin/ma5 -R
>> ma5> set main.fastsim.package = fastjet
>> ma5>import <your HEPMC sample>
>> ma5>plot M(b[1] b[2])
>> ma5>submit
>>
>> —
>> Jack Y. Araz
>>
>>> On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <
>> <email address hidden>> wrote:
>>>
>>> [EXTERNAL EMAIL]
>>>
>>> Question #699618 on MadAnalysis 5 changed:
>>> https://answers.launchpad.net/madanalysis5/+question/699618
>>>
>>> SOUAD SEMLALI gave more information on the question:
>>> Hi Jack,
>>>
>>> Following the previous message, I have checked the histos.saf. The
>>> values are set to 0.
>>>
>>> This is from histos.saf:
>>> <Histo>
>>> <Description>
>>> "3_M"
>>> # nbins xmin xmax
>>> 40 0.000000e+00 2.000000e+02
>>> # Defined regions
>>> myregion # Region nr. 1
>>> </Description>
>>> <Statistics>
>>> 0 0 # nevents
>>> 0.000000e+00 0.000000e+00 # sum of event-weights over events
>>> 0 0 # nentries
>>> 0.000000e+00 0.000000e+00 # sum of event-weights over entries
>>> 0.000000e+00 0.000000e+00 # sum weights^2
>>> 0.000000e+00 0.000000e+00 # sum value*weight
>>> 0.000000e+00 0.000000e+00 # sum value^2*weight
>>> </Statistics>
>>> <Data>
>>> 0.000000e+00 0.000000e+00 # underflow
>>> 0.000000e+00 0.000000e+00 # bin 1 / 40
>>> 0.000000e+00 0.000000e+00 # bin 2 / 40
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00
>>> 0.000000e+00 0.000000e+00 # bin 39 / 40
>>> 0.000000e+00 0.000000e+00 # bin 40 / 40
>>> 0.000000e+00 0.000000e+00 # overflow
>>> </Data>
>>> </Histo>
>>>
>>>
>>> I am using the following syntax:
>>> define q = 5 -5
>>> plot M(q[1] q[2]) 40 0 200 [logY ]
>>>
>>> Am I missing something?
>>>
>>> Many thanks for any help!
>>>
>>> Best
>>> Souad
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadAnalysis 5.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=20
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/madanalysis5/+question/699618
>>
>> You received this question notification because you asked the question.
>>
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Hi

I’m gonna leave the physics part to you, can’t do your project for you. I’m sure you will figure in no time.

cheers

—
Jack Y. Araz

> On Nov 29, 2021, at 3:45 PM, SOUAD SEMLALI <question699618@answers.launchpad.net> wrote:
> 
> [EXTERNAL EMAIL]
> 
> Question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
> 
>    Status: Answered => Open
> 
> SOUAD SEMLALI is still having a problem:
> Dear Jack,
> 
> Thank you for your reply and your help all this time. The plots are
> produced, but I am afraid this brings me back to my starting point. The
> distributions at reco-mode are completely different from the ones at the
> partonic level and confusing. I thought that this is due to the presence of
> large MET, but suppressing MET doesn't solve the problem. As I mentioned
> before, I am working on signal Hsm > z* z* h, with z* > mu+ mu-, z* > j j
> and h > b b~. But The invariant mass distributions are completely different
> from this; distributions with several peaks, and m_{bb~} != m_h and m_{mu+
> mu- jj bb~}!= m_{Hsm}.
> I would be grateful for any suggestions or help!
> 
> Best
> Souad
> 
> 
> Le lun. 29 nov. 2021 à 14:11, Jack Y. Araz <
> question699618@answers.launchpad.net> a écrit :
> 
>> Your question #699618 on MadAnalysis 5 changed:
>> https://answers.launchpad.net/madanalysis5/+question/699618
>> 
>>    Status: Open => Answered
>> 
>> Jack Y. Araz proposed the following answer:
>> you can not use PDGID in reco mode this is not defined. reco mode
>> consists of reconstructed objects not particles. please try the
>> following
>> 
>> $ ./bin/ma5 -R
>> ma5> set main.fastsim.package = fastjet
>> ma5>import <your HEPMC sample>
>> ma5>plot M(b[1] b[2])
>> ma5>submit
>> 
>> —
>> Jack Y. Araz
>> 
>>> On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <
>> question699618@answers.launchpad.net> wrote:
>>> 
>>> [EXTERNAL EMAIL]
>>> 
>>> Question #699618 on MadAnalysis 5 changed:
>>> https://answers.launchpad.net/madanalysis5/+question/699618
>>> 
>>> SOUAD SEMLALI gave more information on the question:
>>> Hi Jack,
>>> 
>>> Following the previous message, I have checked the histos.saf. The
>>> values are set to 0.
>>> 
>>> This is from histos.saf:
>>> <Histo>
>>> <Description>
>>>   "3_M"
>>>   # nbins xmin xmax
>>>     40 0.000000e+00 2.000000e+02
>>>   # Defined regions
>>>     myregion # Region nr. 1
>>> </Description>
>>> <Statistics>
>>>     0 0 # nevents
>>>     0.000000e+00 0.000000e+00 # sum of event-weights over events
>>>     0 0 # nentries
>>>     0.000000e+00 0.000000e+00 # sum of event-weights over entries
>>>     0.000000e+00 0.000000e+00 # sum weights^2
>>>     0.000000e+00 0.000000e+00 # sum value*weight
>>>     0.000000e+00 0.000000e+00 # sum value^2*weight
>>> </Statistics>
>>> <Data>
>>>     0.000000e+00 0.000000e+00 # underflow
>>>     0.000000e+00 0.000000e+00 # bin 1 / 40
>>>     0.000000e+00 0.000000e+00 # bin 2 / 40
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00
>>>     0.000000e+00 0.000000e+00 # bin 39 / 40
>>>     0.000000e+00 0.000000e+00 # bin 40 / 40
>>>     0.000000e+00 0.000000e+00 # overflow
>>> </Data>
>>> </Histo>
>>> 
>>> 
>>> I am using the following syntax:
>>> define q = 5 -5
>>> plot M(q[1] q[2]) 40 0  200 [logY ]
>>> 
>>> Am I missing something?
>>> 
>>> Many thanks for any help!
>>> 
>>> Best
>>> Souad
>>> 
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadAnalysis 5.
>> 
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> 
>> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=20
>> 
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/madanalysis5/+question/699618
>> 
>> You received this question notification because you asked the question.
>> 
> 
> 
> --
> Souad SEMLALI
> 
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message

SOUAD SEMLALI (souad-sem) said on 2021-11-29:

#24

Hi Jack,

Thank you very much for your help all this time!

Best,
Souad

Le lun. 29 nov. 2021 à 15:56, Jack Y. Araz <
<email address hidden>> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> I’m gonna leave the physics part to you, can’t do your project for you.
> I’m sure you will figure in no time.
>
> cheers
>
> —
> Jack Y. Araz
>
> > On Nov 29, 2021, at 3:45 PM, SOUAD SEMLALI <
> <email address hidden>> wrote:
> >
> > [EXTERNAL EMAIL]
> >
> > Question #699618 on MadAnalysis 5 changed:
> > https://answers.launchpad.net/madanalysis5/+question/699618
> >
> > Status: Answered => Open
> >
> > SOUAD SEMLALI is still having a problem:
> > Dear Jack,
> >
> > Thank you for your reply and your help all this time. The plots are
> > produced, but I am afraid this brings me back to my starting point. The
> > distributions at reco-mode are completely different from the ones at the
> > partonic level and confusing. I thought that this is due to the presence
> of
> > large MET, but suppressing MET doesn't solve the problem. As I mentioned
> > before, I am working on signal Hsm > z* z* h, with z* > mu+ mu-, z* > j j
> > and h > b b~. But The invariant mass distributions are completely
> different
> > from this; distributions with several peaks, and m_{bb~} != m_h and
> m_{mu+
> > mu- jj bb~}!= m_{Hsm}.
> > I would be grateful for any suggestions or help!
> >
> > Best
> > Souad
> >
> >
> > Le lun. 29 nov. 2021 à 14:11, Jack Y. Araz <
> > <email address hidden>> a écrit :
> >
> >> Your question #699618 on MadAnalysis 5 changed:
> >> https://answers.launchpad.net/madanalysis5/+question/699618
> >>
> >> Status: Open => Answered
> >>
> >> Jack Y. Araz proposed the following answer:
> >> you can not use PDGID in reco mode this is not defined. reco mode
> >> consists of reconstructed objects not particles. please try the
> >> following
> >>
> >> $ ./bin/ma5 -R
> >> ma5> set main.fastsim.package = fastjet
> >> ma5>import <your HEPMC sample>
> >> ma5>plot M(b[1] b[2])
> >> ma5>submit
> >>
> >> —
> >> Jack Y. Araz
> >>
> >>> On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <
> >> <email address hidden>> wrote:
> >>>
> >>> [EXTERNAL EMAIL]
> >>>
> >>> Question #699618 on MadAnalysis 5 changed:
> >>> https://answers.launchpad.net/madanalysis5/+question/699618
> >>>
> >>> SOUAD SEMLALI gave more information on the question:
> >>> Hi Jack,
> >>>
> >>> Following the previous message, I have checked the histos.saf. The
> >>> values are set to 0.
> >>>
> >>> This is from histos.saf:
> >>> <Histo>
> >>> <Description>
> >>> "3_M"
> >>> # nbins xmin xmax
> >>> 40 0.000000e+00 2.000000e+02
> >>> # Defined regions
> >>> myregion # Region nr. 1
> >>> </Description>
> >>> <Statistics>
> >>> 0 0 # nevents
> >>> 0.000000e+00 0.000000e+00 # sum of event-weights over events
> >>> 0 0 # nentries
> >>> 0.000000e+00 0.000000e+00 # sum of event-weights over entries
> >>> 0.000000e+00 0.000000e+00 # sum weights^2
> >>> 0.000000e+00 0.000000e+00 # sum value*weight
> >>> 0.000000e+00 0.000000e+00 # sum value^2*weight
> >>> </Statistics>
> >>> <Data>
> >>> 0.000000e+00 0.000000e+00 # underflow
> >>> 0.000000e+00 0.000000e+00 # bin 1 / 40
> >>> 0.000000e+00 0.000000e+00 # bin 2 / 40
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00
> >>> 0.000000e+00 0.000000e+00 # bin 39 / 40
> >>> 0.000000e+00 0.000000e+00 # bin 40 / 40
> >>> 0.000000e+00 0.000000e+00 # overflow
> >>> </Data>
> >>> </Histo>
> >>>
> >>>
> >>> I am using the following syntax:
> >>> define q = 5 -5
> >>> plot M(q[1] q[2]) 40 0 200 [logY ]
> >>>
> >>> Am I missing something?
> >>>
> >>> Many thanks for any help!
> >>>
> >>> Best
> >>> Souad
> >>>
> >>> --
> >>> You received this question notification because you are an answer
> >>> contact for MadAnalysis 5.
> >>
> >> --
> >> If this answers your question, please go to the following page to let us
> >> know that it is solved:
> >>
> >>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=20
> >>
> >> If you still need help, you can reply to this email or go to the
> >> following page to enter your feedback:
> >> https://answers.launchpad.net/madanalysis5/+question/699618
> >>
> >> You received this question notification because you asked the question.
> >>
> >
> >
> > --
> > Souad SEMLALI
> >
> > You received this question notification because you are an answer
> > contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=22
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Hi Jack,

Thank you very much for your help all this time!

Best,
 Souad

Le lun. 29 nov. 2021 à 15:56, Jack Y. Araz <
question699618@answers.launchpad.net> a écrit :

> Your question #699618 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
>     Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> I’m gonna leave the physics part to you, can’t do your project for you.
> I’m sure you will figure in no time.
>
> cheers
>
> —
> Jack Y. Araz
>
> > On Nov 29, 2021, at 3:45 PM, SOUAD SEMLALI <
> question699618@answers.launchpad.net> wrote:
> >
> > [EXTERNAL EMAIL]
> >
> > Question #699618 on MadAnalysis 5 changed:
> > https://answers.launchpad.net/madanalysis5/+question/699618
> >
> >    Status: Answered => Open
> >
> > SOUAD SEMLALI is still having a problem:
> > Dear Jack,
> >
> > Thank you for your reply and your help all this time. The plots are
> > produced, but I am afraid this brings me back to my starting point. The
> > distributions at reco-mode are completely different from the ones at the
> > partonic level and confusing. I thought that this is due to the presence
> of
> > large MET, but suppressing MET doesn't solve the problem. As I mentioned
> > before, I am working on signal Hsm > z* z* h, with z* > mu+ mu-, z* > j j
> > and h > b b~. But The invariant mass distributions are completely
> different
> > from this; distributions with several peaks, and m_{bb~} != m_h and
> m_{mu+
> > mu- jj bb~}!= m_{Hsm}.
> > I would be grateful for any suggestions or help!
> >
> > Best
> > Souad
> >
> >
> > Le lun. 29 nov. 2021 à 14:11, Jack Y. Araz <
> > question699618@answers.launchpad.net> a écrit :
> >
> >> Your question #699618 on MadAnalysis 5 changed:
> >> https://answers.launchpad.net/madanalysis5/+question/699618
> >>
> >>    Status: Open => Answered
> >>
> >> Jack Y. Araz proposed the following answer:
> >> you can not use PDGID in reco mode this is not defined. reco mode
> >> consists of reconstructed objects not particles. please try the
> >> following
> >>
> >> $ ./bin/ma5 -R
> >> ma5> set main.fastsim.package = fastjet
> >> ma5>import <your HEPMC sample>
> >> ma5>plot M(b[1] b[2])
> >> ma5>submit
> >>
> >> —
> >> Jack Y. Araz
> >>
> >>> On Nov 29, 2021, at 1:55 PM, SOUAD SEMLALI <
> >> question699618@answers.launchpad.net> wrote:
> >>>
> >>> [EXTERNAL EMAIL]
> >>>
> >>> Question #699618 on MadAnalysis 5 changed:
> >>> https://answers.launchpad.net/madanalysis5/+question/699618
> >>>
> >>> SOUAD SEMLALI gave more information on the question:
> >>> Hi Jack,
> >>>
> >>> Following the previous message, I have checked the histos.saf. The
> >>> values are set to 0.
> >>>
> >>> This is from histos.saf:
> >>> <Histo>
> >>> <Description>
> >>>   "3_M"
> >>>   # nbins xmin xmax
> >>>     40 0.000000e+00 2.000000e+02
> >>>   # Defined regions
> >>>     myregion # Region nr. 1
> >>> </Description>
> >>> <Statistics>
> >>>     0 0 # nevents
> >>>     0.000000e+00 0.000000e+00 # sum of event-weights over events
> >>>     0 0 # nentries
> >>>     0.000000e+00 0.000000e+00 # sum of event-weights over entries
> >>>     0.000000e+00 0.000000e+00 # sum weights^2
> >>>     0.000000e+00 0.000000e+00 # sum value*weight
> >>>     0.000000e+00 0.000000e+00 # sum value^2*weight
> >>> </Statistics>
> >>> <Data>
> >>>     0.000000e+00 0.000000e+00 # underflow
> >>>     0.000000e+00 0.000000e+00 # bin 1 / 40
> >>>     0.000000e+00 0.000000e+00 # bin 2 / 40
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00
> >>>     0.000000e+00 0.000000e+00 # bin 39 / 40
> >>>     0.000000e+00 0.000000e+00 # bin 40 / 40
> >>>     0.000000e+00 0.000000e+00 # overflow
> >>> </Data>
> >>> </Histo>
> >>>
> >>>
> >>> I am using the following syntax:
> >>> define q = 5 -5
> >>> plot M(q[1] q[2]) 40 0  200 [logY ]
> >>>
> >>> Am I missing something?
> >>>
> >>> Many thanks for any help!
> >>>
> >>> Best
> >>> Souad
> >>>
> >>> --
> >>> You received this question notification because you are an answer
> >>> contact for MadAnalysis 5.
> >>
> >> --
> >> If this answers your question, please go to the following page to let us
> >> know that it is solved:
> >>
> >>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=20
> >>
> >> If you still need help, you can reply to this email or go to the
> >> following page to enter your feedback:
> >> https://answers.launchpad.net/madanalysis5/+question/699618
> >>
> >> You received this question notification because you asked the question.
> >>
> >
> >
> > --
> > Souad SEMLALI
> >
> > You received this question notification because you are an answer
> > contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/699618/+confirm?answer_id=22
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/699618
>
> You received this question notification because you asked the question.
>

-- 
Souad SEMLALI

Revision history for this message

Launchpad Janitor (janitor) said on 2021-12-15:

#25

This question was expired because it remained in the 'Open' state without activity for the last 15 days.

Revision history for this message

Benjamin Fuks (fuks) said on 2021-12-15:

#26

The flag was no updated. This is now done.

MadAnalysis 5

missing energy

Question information

Subscribers

MadAnalysis 5

missing energy

Question information

Related bugs

Related FAQ:

Subscribers