Fork me on GitHub

Opened 11 years ago

Last modified 11 years ago

#257 new Bug

Delphes3.12 is not creating leptonic events

Reported by: Arindam Das Owned by:
Priority: minor Milestone:
Component: Delphes code Version: Delphes 3
Keywords: Cc:

Description

Dear Delphes Team ,

Delphes3.12 does not record the

p p > mu+ mu- events.

It is only recording the neutrinos.

A sample of 300 events are attached.

Regards,
Arindam

Attachments (4)

tag_1_delphes_events.lhco (77.5 KB ) - added by Arindam Das 11 years ago.
tag_1_pythia_events.hep.gz (1.2 MB ) - added by Arindam Das 11 years ago.
run_01_tag_1_banner.txt (43.6 KB ) - added by Arindam Das 11 years ago.
run_01_tag_1_banner.2.txt (43.6 KB ) - added by Arindam Das 11 years ago.

Download all attachments as: .zip

Change History (12)

by Arindam Das, 11 years ago

Attachment: tag_1_delphes_events.lhco added

comment:1 by Michele Selvaggi, 11 years ago

Hi Arindam,

Are you converting the root file to lhco with the root2lhco program inside Delphes or do you let MG handle it?
I am not sure how the root to lhco conversion is done by MG (if it is done at all...).

I suggest you do the following:

  • install from scratch the last Delphes version in some dir

(dowload the tar file, untag it, and type "make")

  • take the *.hep file from MG (say it is called events.hep)
  • run the following commands from your Delphes dir:
./DelphesSTDHEP examples/delphes_card_CMS.tcl events.root events.hep 
./root2lhco events.root events.lhco

If this still outputs some strange lhco file like the one you have attached, then there might be some issue with your MG events.

Let me know if this solves your issue.

comment:2 by Arindam Das, 11 years ago

Dear Michele,

I did the same what you asked me to do,

I compiled with the standalone Delphes using the tag_1_pythia...hep file from MadGraph (This is the only .hep file I got).

Beside that I used the Madgraph commands to produce the .lhco file which comes at the time of the evnt run. In that case I had nothing to do everything was done by the MadGraph.

In both of the cases the results are same. I mean as before.

What shall I do now ? Is this a MadGraph error? Is it creating a wrong tag_1_pythia_events.hep file for Delphes ?

Regards,

Arindam

comment:3 by Michele Selvaggi, 11 years ago

Could you please attach a small hep event file (~100 events)?

in reply to:  3 comment:4 by Arindam Das, 11 years ago

Replying to mselvaggi:

Could you please attach a small hep event file (~100 events)?

The sample hep file is here. It contains 200 events.

Regards,

Arindam

by Arindam Das, 11 years ago

Attachment: tag_1_pythia_events.hep.gz added

comment:5 by Arindam Das, 11 years ago

Dear Michele,

I have attached the .hep file.

Regards,
Arindam

comment:6 by Michele Selvaggi, 11 years ago

Hi Arindam,

your muons seem to be produced with very small pt and high rapidity.
Therefore they don't pass the cuts in the Delphes card (only muons with pt>10 GeV and eta<2.5 are stored).

You might want to put some cuts in the MG run_card to veto these (unobservable) events during generation.
Or, if you really intended to look at such events, you might want to change the efficiency formula in the delphes card (around line 358)

Michele

Last edited 11 years ago by Michele Selvaggi (previous) (diff)

by Arindam Das, 11 years ago

Attachment: run_01_tag_1_banner.txt added

comment:7 by Arindam Das, 11 years ago

Dear Michele,

is that wrong if I remove cuts from the run_card.dat ?

If I want to generate events without cuts then what should I do ?

I am attached the run_card.dat . Will you please enlighten me ?

Regards,
Arindam

by Arindam Das, 11 years ago

Attachment: run_01_tag_1_banner.2.txt added

comment:8 by Arindam Das, 11 years ago

Dear Michele,

If I want to create the events without cuts and then want to put cuts from some other file on the events.lhco then what should I do ?

My choice is :

I do not want to put any cut from run_card.dat.

I want the .lhco file with al the events exactly no cut.

Then I want to put cuts from the other file on the .lhco file.

How do I change the Delphes file to do that ? Because now a days the detector cuts are different which is easy to make in other format for me than run_card.dat .

Otherwise if I want to check the cross section and the events in the detector file before putting any cut then what should I do ? How to change the Detector file ?

You wrote to change the efficiency formula around the line 358 in the delphes card . How to do that ? What should be the correct formula and how to write it ?

Please enlighten me .

Regards,
Arindam

Note: See TracTickets for help on using tickets.