Fork me on GitHub

Opened 11 years ago

Last modified 11 years ago

#227 new How to

delphes 3 have included fake candidates ?

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

Description

Hello,

In delphes2, there are no fake candidates entering electron muon and photon collections. For delphes 3, I think it have included fake candidates in default delphes_card_ATLAS.tcl. I don't sure whether I am right ?

Cheers,
Gavin

Change History (9)

comment:1 by Gavin, 11 years ago

In delphes2, electrons and photons may create a candidate in the jet collection. We must use "EhadOverEem" to clean lepton collections from jets. But in delphes3, I don't need to do that. Do you agree with me?

comment:2 by Michele Selvaggi, 11 years ago

Hi,

the UniqueObjectFinder module already takes care of cleaning leptons from jets (as long leptons as leptons pass the isolation criteria)

Michele

comment:3 by Gavin, 11 years ago

Hi,

The paper (http://indico.cern.ch/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=266825) has used delphes 2 and a constant jet-to-photon fake rate of 10(-3) to simulate zphotonjet production. I have simulated zphotonjet process with delphes 3. But my result with delphes 3 is larger one order than that paper. (I have checked that paper using delphes 2).

Maybe I don't express clearly. I just want to know about the jet-to-photon fake situation in delphes 3.

Thanks,
Gavin


comment:4 by Michele Selvaggi, 11 years ago

Hi Gavin,

If I understand correctly you are observing more photons than you are expecting. The isolation criteria may have changed between Delphes2 and Delphes3. Could you try to tighten the photon isolation criteria in your detector card? Or check that the isolation is compatible with that of Delphes 2.

I am not sure, however, how you could set some "jet-to-photon fake rate". To my knowledge this was not possible in Delphes2. Did you apply this directly to the Delphes output? As a side comment, please be aware that it is now possible to do this with Delphes3. You can create a module that converts reconstructed jets into photons, according to some parametrised probability. Let us know if you need this.

Michele

comment:5 by Pavel Demin, 11 years ago

Just for information.

There was no photon isolation criteria in Delphes 2.

in reply to:  4 comment:6 by Gavin, 11 years ago

Replying to mselvaggi:
Hi Michele,

If I understand correctly you are observing more photons than you are expecting. The isolation criteria may have changed between Delphes2 and Delphes3. Could you try to tighten the photon isolation criteria in your detector card? Or check that the isolation is compatible with that of Delphes 2.

There was no photon isolation criteria, so I didn't consider photon isolation in delphes 2.

I am not sure, however, how you could set some "jet-to-photon fake rate". To my knowledge this was not possible in Delphes2. Did you apply this directly to the Delphes output?

I just scaled the total events number with "jet-to-photon fake rate" when I dealt with delphes2 output.

As a side comment, please be aware that it is now possible to do this with Delphes3. You can create a module that converts reconstructed jets into photons, according to some parametrised probability. Let us know if you need this.

Yes, I need jet-fake-to-photon events. But I don't know how to do in delphes3.

Cheers,
Gavin

comment:7 by Michele Selvaggi, 11 years ago

Hi Gavin,

just to be sure, could you precise the following: do you have too many photons (or too few) in Delphes 3?

You can tune the number of photons by relaxing (if you need more photons) or tighten (if you need less photons) the isolation criteria.
Moreover if want to reproduce exactly the Delphes2 results you can relax totally the isolation criteria.

To do so you can set the parameter to a large value:
set PTRatioMax 1000

Michele

comment:8 by Gavin, 11 years ago

Hi Michele,

I use the default delphes_atlas_card in delphes2 and delphes3. But I relaxed totally the isolation criteria with PTRatioMax at 1000 in delphes3. I analysed the both delphes output files with same cuts. The results is below,

delphes2 27866 (I didn't use the jet-fake-photon rate)
delphes3 3124

Cheers,
Gavin

comment:9 by Michele Selvaggi, 11 years ago

Hi Gavin,

thanks for the testing. Could you tell whether you are using the eflow algorithm in Delphes2?
In the detector card, the corresponding flag JET_Eflow should be to set to false (i.e 0)

If I remember correctly the photon collection was bugged in Delphes 2.
Could you also provide some small event sample (100 events) so that we can do some cross-check?

Thanks

Michele

Note: See TracTickets for help on using tickets.