New particles in reco-level

Asked by SAUMYEN KUNDU on 2019-07-21

Hi,

I was using two HEPMC files (Signal and Background) genereated by MG5 in the reco-level of MA5 (version 1.7.21) with fastsim package Delphes and detector CMS.
Initially, I didn't modify anything (even the Delphes card) and just generated some plots. But there was no distribution for ET(invisible), E(invisible), etc. And also there was some warning while running the Delphes like in below:

WARNING: --------------------------------------------------------------------------------
WARNING: Msg | Unknown particle by Delphes in the final state
WARNING: Details | particle with PDGID=6000005
WARNING: Where | Function = TranslateMA5toDELPHES ; File=delphes/DetectorDelphes.cpp ; Line=295
WARNING: --------------------------------------------------------------------------------
WARNING: --------------------------------------------------------------------------------
WARNING: Msg | Unknown particle by Delphes in the final state
WARNING: Details | particle with PDGID=-6000005
WARNING: Where | Function = TranslateMA5toDELPHES ; File=delphes/DetectorDelphes.cpp ; Line=295
WARNING: --------------------------------------------------------------------------------

So I again ran the script with adding the line 'define invisible = invisible +6000005 -6000005'. And also modify the Delphes card. The modifications are :
1. Added these lines under ECAL and HCAL (under 'SimpleCalorimeter' module)
                     add EnergyFraction {6000005} {0.0}
2. Added these lines under Neutrino Filter (under 'PdgCodeFilter' module)
                     add PdgCode {6000005}
                     add PdgCode {-6000005}

But there was no change in the output. And the same warning was there while running the Delphes. So my question is thet whatever modifications I made are correct or not; I mean if I should modify these thing like this or not? And what else should I modify? Basically these paticles should contribute in the missing energy.

Another question is that regarding using a particular detector card. How can I use the detector card for the ILD?

I will be very thankful for any suggestions.

Regards,
Saumyen

Question information

Language:
English Edit question
Status:
Answered
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Last query:
2019-07-22
Last reply:
2019-07-25
Benjamin Fuks (fuks) said : #1

Hi Saumyen,

The reason of the issue is that in the input hepmc file, there is a particle with a PDG code 6000005 that Delphes does not know. In order to have delphes properly running, one needs to indicate it the properties of this particle.

> 1. Added these lines under ECAL and HCAL (under 'SimpleCalorimeter' module)
> add EnergyFraction {6000005} {0.0}
This is incorrect. One needs to provide information on what the particle does both int the hcal and ecal:
  add EnergyFraction {6000005} {0.0 0.0}

> 2. Added these lines under Neutrino Filter (under 'PdgCodeFilter' module)
> add PdgCode {6000005}
> add PdgCode {-6000005}
This is correct.

> Another question is that regarding using a particular detector card. How can I use the detector card for the ILD?
MadAnalysis always asks you whether you want to modify the detector card. You can then chose any existing card (atlas or cms) and supersede them with the ILD one.

Cheers,

Benjamin

SAUMYEN KUNDU (saumyen.k) said : #2

Hi Benjamin,

Thanks a lot for your response.

>The reason of the issue is that in the input hepmc file, there is a particle with a PDG code 6000005 that Delphes does not know. In >order to have delphes properly running, one needs to indicate it the properties of this particle.

What else do I need to modify other than adding those lines under ECAL, HCAL and Neutrino Filter?

>This is incorrect. One needs to provide information on what the particle does both in the hcal and ecal:
> add EnergyFraction {6000005} {0.0 0.0}

I added that line both under ECAL as well as HCAL. I thought that I can provide the information on what to do in the ECAL and HACAL. So, in both the cases I have to enter 'add EnergyFraction {6000005} {0.0 0.0}'!

>This is correct.

OK!

>MadAnalysis always asks you whether you want to modify the detector card. You can then chose any existing card (atlas or cms) >and supersede them with the ILD one.
Yeah, this I saw. I did once too. But then the Delphes didn't run and also this is very tiring. Because starting from the tracking efficiencies, the resolution formula a lot of things are to be changed and doing that in the editor like 'vim' is tiring. So, I was looking for alternatives. Ok, let me try again.

But I had another confusion regarding the stacking method. What actually happens when we choose the method to be 'stack' or 'superimpose'. I mean, I went through the manual. It says "the different curves can be stacked (stack, default),
superimposed (superimpose) ". But what is it meant by 'stacked' and 'superimposed'. I am asking this because I generated histograms for the same datasets and same observable with choices 'stack' and 'superimpose, respectively. In both cases, the nature of the histograms was the same. But the values are not. I had one background and one signal. With the choice 'stack' I found the signal to dominate the background and with the peak value ( for a certain observable) at a little above 1 while choosing it 'superimpose' the background dominates and the signal peak decreases to below 1. In the bothe values for the background remains the same. This character is there for all the observable. So, I wanted to know what actually happens when we choose 'stack' or 'superimpose' as the stacking method?

Regards,
Saumyen

Benjamin Fuks (fuks) said : #3

Hi,

> What else do I need to modify other than adding those lines under ECAL,
> HCAL and Neutrino Filter?
Only what has been said in the previous message..

> I added that line both under ECAL as well as HCAL. I thought that I can
> provide the information on what to do in the ECAL and HACAL. So, in both
> the cases I have to enter 'add EnergyFraction {6000005} {0.0 0.0}'!
Yes

>
>> MadAnalysis always asks you whether you want to modify the detector card. You can then chose any existing card (atlas or cms) >and supersede them with the ILD one.
> Yeah, this I saw. I did once too. But then the Delphes didn't run and also this is very tiring.
This is maybe because your card is incorrect. Have you try it externally (from tools/delphes)?

> Because starting from the tracking efficiencies, the resolution formula a lot of things are to be changed and doing that in the editor like 'vim' is tiring. So, I was looking for alternatives. Ok, let me try again.
Well, I won't do your work for you :)

> But I had another confusion regarding the stacking method. What actually happens when we choose the method to be 'stack' or 'superimpose'. I mean, I went through the manual. It says "the different curves can be stacked (stack, default),
> superimposed (superimpose) ". But what is it meant by 'stacked' and 'superimposed'. I am asking this because I generated histograms for the same datasets and same observable with choices 'stack' and 'superimpose, respectively. In both cases, the nature of the histograms was the same. But the values are not. I had one background and one signal. With the choice 'stack' I found the signal to dominate the background and with the peak value ( for a certain observable) at a little above 1 while choosing it 'superimpose' the background dominates and the signal peak decreases to below 1. In the bothe values for the background remains the same. This character is there for all the observable. So, I wanted to know what actually happens when we choose 'stack' or 'superimpose' as the stacking method?
When the two histograms are stacked, they are put on top of each other. When they are superimposed… well they are superimposed (I do not know any other synonym).

Cheers,

Benjamin

Can you help with this problem?

Provide an answer of your own, or ask SAUMYEN KUNDU for more information if necessary.

To post a message you must log in.