Fork me on GitHub

Opened 8 years ago

Last modified 8 years ago

#1061 new Bug

Distribution of transvere W boson mass not modelled correctly by Delphes 3.4.0

Reported by: Soeren Stamm Owned by:
Priority: minor Milestone:
Component: Delphes code Version: Delphes 3
Keywords: ATLAS Cc:

Description

Hi,

I have generated single-top s-channel events (pp > t b~) with MadGraph5 2.4.2 and stored them in LHE files. These files are provided to Delphes 3.4.0 and showered/decayed by Pythia 8.219 using the executable "DelphesPythia8".

I found that the distribution of the transverse mass of the W boson (mtW) from the top decay is not modelled as expected. I did some research and was able to get a good modelling of the mtW distribution with Delphes 3.3.3 (and lower).

I attached a plot of mtW for Delphes 3.4.0 (red line) compared to Delphes 3.3.3 (blue line). The blue line is in agreement with my particle level study using Rivet. I also attached distributions for lepton pt (e/mu) and missing transverse energy and jet pt (highest pt jet only).

The content of the attached pdf is: page 1 - lepton pt, p.2 - jet pt, p.3 - mtW, p.4 - missing transverse energy.

What was changed between version 3.3.3 and 3.4.0 that could cause this odd looking shape of the mtW distribution?

Cheers,
Soeren

Attachments (4)

DELPHES.pdf (23.0 KB ) - added by Soeren Stamm 8 years ago.
Distributions for lepton pt, jet pt, mtW and Et miss
DelphesFastAnalysis.C (18.0 KB ) - added by Soeren Stamm 8 years ago.
Script for creating mtW plot
events_run1.lhe.gz (211.6 KB ) - added by Soeren Stamm 8 years ago.
2000 events MG5 pp > tb~ + pp > t~b
events_run2.lhe.gz (211.6 KB ) - added by Soeren Stamm 8 years ago.
2000 events MG5 pp > tb~ + pp > t~b

Download all attachments as: .zip

Change History (9)

by Soeren Stamm, 8 years ago

Attachment: DELPHES.pdf added

Distributions for lepton pt, jet pt, mtW and Et miss

comment:1 by Soeren Stamm, 8 years ago

I am using the ATLAS card "delphes_card_ATLAS.tcl".

comment:2 by Michele Selvaggi, 8 years ago

Hi Soren,

thanks for spotting this. Can you please attach a small (~2k evts) LHE event sample and the root macro that produces the Mtw plot?
I can then have a look at when things started to go wrong.

Thanks
Michele

by Soeren Stamm, 8 years ago

Attachment: DelphesFastAnalysis.C added

Script for creating mtW plot

by Soeren Stamm, 8 years ago

Attachment: events_run1.lhe.gz added

2000 events MG5 pp > tb~ + pp > t~b

by Soeren Stamm, 8 years ago

Attachment: events_run2.lhe.gz added

2000 events MG5 pp > tb~ + pp > t~b

in reply to:  2 comment:3 by Soeren Stamm, 8 years ago

Hi Michele,

thanks for looking into this.

I attached a root script "DelphesFastAnalysis.C" that you can use to analyse Delphes root files. I usually compile it using the ROOT cint interpreter via ".L DelphesFastAnalysis.C+" and then run it via DelphesFastAnalysis("input.root", "output.root").
Please let me know if you have trouble to run it. I am using ROOT 5.34.

I also attached 2x 2000 evts for pp > tb~ + pp > t~b.

Cheers,
Soeren

Replying to mselvaggi:

Hi Soren,

thanks for spotting this. Can you please attach a small (~2k evts) LHE event sample and the root macro that produces the Mtw plot?
I can then have a look at when things started to go wrong.

Thanks
Michele

comment:4 by Michele Selvaggi, 8 years ago

Hi Soeren,

Given your final state the MT distribution is going to be very sensitive to which MET and isolation definiton are adopted and both these things have changed from the card in 3.3.3 and 3.4.0. Unfortunately the signal selection is very small so I can't really see if the changes I propose are really effective.
You'll have to try the following with larger stats:

Replace the inputArray to MissingET by "set InputArray EFlowMerger/eflow"

This will put back leptons in the MissingET definition (with EFlowMerger/eflow they are excluded)

Another thing you can try is putting back leptons in the isolation definition. So in every module XXXIsolation (where XXX=Electron/Photon/Muon) you can try to replace the IsolationInputArray EFlowFilter/eflow by EFlowMerger/eflow.

These changes should get you close to the card in 3.3.3.

Cheers
Michele

in reply to:  4 comment:5 by Soeren Stamm, 8 years ago

Hi Michele,

I tried two different setups:

1st setup) I replaced the input array Calorimeter/towers with EFlowMerger/eflow in

module Merger MissingET {
# add InputArray InputArray
# add InputArray Calorimeter/towers

set InputArray EFlowMerger/eflow
set MomentumOutputArray momentum

}

I left the isolation as it was (i.e. so still using "set IsolationInputArray EFlowFilter/eflow")

This gives the expected mtW distribution. Very good!

However, I don't understand why this works... (I just started a week ago to look into Delphes)

  • The EFlowMerger/eflow are made from HCal/eflowTracks, ECal/eflowPhotons, HCal/eflowNeutralHadrons.

I don't see electrons here?

Are they included via HCal/eflowTracks, i.e. the four-momentum of all electrons are included through their tracking information contained in the tracks? HCal/eflowTracks would then also include all charged particle track momenta (i.e. pions, etc.)?

  • The Calorimeter/towers are made from ECal/ecalTowers, HCal/hcalTowers.

Shouldn't this include contributions from electrons and photons to the ecalTowers?

Is there any documentation on how the modules work?
I would like to optimize the ATLAS card to get better results with my analysis. For this I would like to understand the code a bit better.

2nd setup) I replaced in all XXXIsolationModules the lines:

set IsolationInputArray EFlowFilter/eflow
by
set IsolationInputArray EFlowMerger/eflow

but kept the Calorimeter/towers in the MissingET computation. Only changing the isolation, does not fix the problem.

Thank you, cheers,
Soeren

Replying to mselvaggi:

Hi Soeren,

Given your final state the MT distribution is going to be very sensitive to which MET and isolation definiton are adopted and both these things have changed from the card in 3.3.3 and 3.4.0. Unfortunately the signal selection is very small so I can't really see if the changes I propose are really effective.
You'll have to try the following with larger stats:

Replace the inputArray to MissingET by "set InputArray EFlowMerger/eflow"

This will put back leptons in the MissingET definition (with EFlowMerger/eflow they are excluded)

Another thing you can try is putting back leptons in the isolation definition. So in every module XXXIsolation (where XXX=Electron/Photon/Muon) you can try to replace the IsolationInputArray EFlowFilter/eflow by EFlowMerger/eflow.

These changes should get you close to the card in 3.3.3.

Cheers
Michele

Note: See TracTickets for help on using tickets.