Fork me on GitHub

Opened 8 years ago

#994 new Bug

Usage of BTag in CMS PhaseII card and PyRoot example

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

Description

Dear Developers

I an using the version of Delphes updated to the head yesterday (#591)
I am trying to analyze a file from the card

delphes/cards/CMS_PhaseII/CMS_PhaseII_0PU.tcl

My intention is to use the efficiencies for BTag working points implemented there.

I have some questions on the usage of the Jets class:

1) I see three BTag flags

a) print "btag "+str(jet.BTag) == it goes from 0 to at least 7
b) print "btag algo "+str(jet.BTagAlgo)
c) print "btag algo phys "+str(jet.BTagPhys)

browsing internet and in the files I could not figure out: Where does the BitNumber is?

For example, If I want to use medium BTag which is the flag for false/true ?
Where are the weights of the efficiency map stored?

2) Even if the BTag modules appear in FatJet, those are not equivalent to a FatJet BTag, correct?

3) I am using PyRoot.
On your Example1.py I had to add the following lines:

ROOT.gInterpreter.Declare('#include "ExRootAnalysis/ExRootTreeReader.h"')
ROOT.gInterpreter.Declare('#include "classes/DelphesClasses.h"')

And in DelphesClasses.h I had to substitute, for example

#include "classes/SortableObject.h"
to
#include "SortableObject.h"

Did anybody else using PyRoot had the same issue?
I am sure I am making this wrong.

For example, I am able to access basic information of the Jet branch (M,PT,…),
but If I try to access the pruned jet as bellow, I have a crash

jet = branchFatJet.At(part)
prunedJet = jet.PrunedP4[0]
or
prunedJetM = jet.PrunedP4[0].M()

Do you think is something related with the ‘fix’ above?
A part of the log is bellow.

4) Finally: If I use this implementation for a publication, is there any reference I should use in special to recognize your implementation?
(apart from the CMS papers from where the efficiencies are from and Delphes itself)

Thank you
Alexandra

========================================
#8 0x00007f48a92fbef7 in Cppyy::GetActualClass(long, void*) () from /afs/cern.ch/cms/slc6_amd64_gcc493/lcg/root/6.06.00-ikhhed2/lib/libPyROOT.so
#9 0x00007f48a92f2afb in PyROOT::BindCppObject(void*, long, bool) () from /afs/cern.ch/cms/slc6_amd64_gcc493/lcg/root/6.06.00-ikhhed2/lib/libPyROOT.so
#10 0x00007f48a930f922 in PyROOT::TTupleOfInstances_New(void*, long, long) () from /afs/cern.ch/cms/slc6_amd64_gcc493/lcg/root/6.06.00-ikhhed2/lib/libPyROOT.so
#11 0x00007f48a92f5b22 in PyROOT::(anonymous namespace)::pp_get(PyROOT::PropertyProxy*, PyROOT::ObjectProxy*, _object*) () from /afs/cern.ch/cms/slc6_amd64_gcc493/lcg/root/6.06.00-ikhhed2/lib/libPyROOT.so

Change History (0)

Note: See TracTickets for help on using tickets.