Fork me on GitHub

Opened 11 years ago

Closed 11 years ago

#196 closed Bug (fixed)

DelphesCMSFWLite output question

Reported by: dimatteo Owned by:
Priority: major Milestone:
Component: Delphes code Version: Delphes 3
Keywords: Cc:

Description

Dear Developers,

we are running DelphesCMSFWLite on Delphes Delphes-3.0.9.
We can technically produce an output out of a CMSSW RECO sample [1], but the output [2] looks not okay.
For instance, check the delphes electron Eta distribution (attached plot), the peak at zero is non-physical.
Any filtering we should apply?
What could be the reason of such behavior?

Kind regards
Leonardo

[1]
In a CMS storage
/store/mc/Summer13/DYToEE_M-20_TuneZ2star_14TeV-pythia6/GEN-SIM-RECO/UpgradePhase1Age0START_DR61SLHCx_PU140Bx25_STAR17_61_V1A-v1/10000/FE0AFAB4-A1CB-E211-AC35-002354EF3BE1.root

[2]
if you have access to cern afs
/afs/cern.ch/work/s/sixie/public/releases/analysis/CMSSW_5_3_9_patch3/src/Delphes-3.0.9/outputTest.root

Attachments (1)

Delphes_out_ele_eta.gif (4.3 KB ) - added by dimatteo 11 years ago.

Download all attachments as: .zip

Change History (8)

by dimatteo, 11 years ago

Attachment: Delphes_out_ele_eta.gif added

comment:1 by Pavel Demin, 11 years ago

Dear Leonardo,

Thank you for testing DelphesCMSFWLite.

There is a bug in DelphesCMSFWLite.cpp. The particle's pz is wrongly set as 'pz = particle.phi()'.

Here is a preliminary version with this bug fixed:
http://cp3.irmp.ucl.ac.be/downloads/Delphes-3.0.10.pre4.tar.gz

Could you, please check if this version works better?

Regards,

Pavel

in reply to:  1 comment:2 by Si Xie, 11 years ago

Hi Pavel,

I verify that the new version fixes the above issue.

Thanks.
Si

comment:3 by Si Xie, 11 years ago

Hi,

A different question related to the CMSFWLite set up.

I'm trying to run it on a CMS AOD file which has HH->bbtautau.
However, the output file:

/afs/cern.ch/work/s/sixie/public/Phase2Upgrade/Delphes/ntuples/Delphes.dihiggs-bbgg-14tev.root

does not have any jets with tau tag = 1.

Our card file is here:

/afs/cern.ch/work/s/sixie/public/releases/analysis/CMSSW_5_3_9_patch3/src/Delphes-3.0.10/delphes_card_Snowmass_NoPileUp.tcl

and the input file is here:
/afs/cern.ch/work/s/sixie/public/releases/analysis/CMSSW_5_3_9_patch3/src/Delphes-3.0.10/AOD-058201-step2.hhbbtt.root

Can it be because in the cards we do not explicitly specify "set BitNumber 0" for the tau tag?

Thanks.
Si

comment:4 by Si Xie, 11 years ago

Component: Delphes miscellaneousDelphes code
Priority: minormajor

Hi,

I have tried using the newest CMS phase2 detector cards, and I still do not get any jets with tau tag = 1. I tried adding "set BitNumber 0" but this didn't have any effect. Can somebody look into this? Is it a problem with DelphesCMSFWLite?

delphes card that I used:
/afs/cern.ch/work/s/sixie/public/releases/analysis/CMSSW_5_3_9_patch3/src/Delphes-3.0.10/phaseII_conf3_140pu_v2.sixie.tcl

input cms file:
/afs/cern.ch/work/s/sixie/public/releases/analysis/CMSSW_5_3_9_patch3/src/Delphes-3.0.10/AOD-058201-step2.hhbbtt.root

and i'm using delphes 3.0.10.

Thanks.
Si

comment:5 by Pavel Demin, 11 years ago

There is a problem in DelphesCMSFWLite.cpp: the daughters indices are not set but they are required for the tau-tagging module to select hadronic tau decays.

Here is a preliminary version with this bug fixed:

http://cp3.irmp.ucl.ac.be/downloads/Delphes-3.0.10.pre7.tar.gz

I've tried it with your input and config files and some tau jets start to appear.

comment:6 by Si Xie, 11 years ago

OK Thanks.
I tried it and also got some taus now.
We will check on the full sample now and report back if everything is ok.

Si

comment:7 by Pavel Demin, 11 years ago

Resolution: fixed
Status: newclosed

Fixed version of the DelphesCMSFWLite.cpp code is available in the new release 3.0.10.

Note: See TracTickets for help on using tickets.