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