Opened 6 years ago
Last modified 6 years ago
#1349 new How to
Calculating Missing ET after Smearing Jets
Reported by: | keanet | Owned by: | |
---|---|---|---|
Priority: | minor | Milestone: | |
Component: | Delphes code | Version: | Delphes 3 |
Keywords: | MET, Smeared Fast Jets | Cc: |
Description
Hi, I wonder how I can recalculate missing ET after incorporating the effects of smearing fast jets using Delphes.
I am using EFlowMerger/eflow as input for the FastJetFinder to calculate the fastjets before smearing the momenta of the fastjets (see the attached Delphes card).
Someone gave me the idea that I can do the following:
- remove all the constituents of the smeared jets from EFlowMerger/eflow
- sum the momenta of the smeared jets and of the objects that remained in EFlowMerger/eflow after the previous step
but I am not sure how to remove all the constituents of the smeared jets from EFlowMerger/eflow. Do I have to tweak the FastJetFinder module to do it? If so, can you give me some guidance on how to do that? Or is it possible to just modify the Delphes card that I am using to accomplish the goal?
Thank you!
Attachments (1)
Change History (5)
by , 6 years ago
Attachment: | delphes_card_CMS_smear.tcl added |
---|
comment:1 by , 6 years ago
comment:2 by , 6 years ago
Thank you pavel for your reply. After modifying the FastJetFinder module and the Delphes card, I got an error after running Delphes...
I have made the following changes to the FastJetFinder module:
starting from line 291:
fRhoOutputArray = ExportArray(GetString("RhoOutputArray", "rho")); fConstArray = ExportArray(GetString("ConstArray", "constituents"));
starting from line 451:
fConstArray->AddCandidate(constituent); candidate->AddCandidate(constituent); }
To the Delphes card I added the following changes (shown in bold):
############
# Jet finder
############
module FastJetFinder FastJetFinder {
# set InputArray Calorimeter/towers
set InputArray EFlowMerger/eflow
set OutputArray jets
set ConstArray constituents
...
#####################################################
# Find uniquely identified photons/electrons/tau/jets
#####################################################
module UniqueObjectFinder UniqueObjectFinder {
add InputArray FastJetFinder/constituents constituents
add InputArray EFlowMerger/eflow eflow
add InputArray PhotonIsolation/photons photons
add InputArray ElectronIsolation/electrons electrons
add InputArray MuonIsolation/muons muons
add InputArray JetEnergyScale/jets jets
}
###################
# Missing ET merger
###################
module Merger MissingET {
add InputArray UniqueObjectFinder/eflow
add InputArray JetsMomentumSmearing/jets
set MomentumOutputArray momentum
}
But when I ran Delphes, I get the following error:
ERROR: can't access input list 'FastJetFinder/constituents' in module 'UniqueObjectFinder'
I have also tried doing on line 451
fConstArray->Add(constituent) instead of ->AddCandidate(constituent), but I still got the same error.
comment:3 by , 6 years ago
ERROR: can't access input list 'FastJetFinder?/constituents' in module 'UniqueObjectFinder?'
Looks like something wasn't rebuilt. Could you try to completely rebuild Delphes with the following commands?
make clean rm -rf tmp make
While adding eflow and constituents to the already existing UniqueObjectFinder would probably work, I think that creating a separate UniqueObjectFinder would be a cleaner solution.
comment:4 by , 6 years ago
You were right. I needed to rebuild Delphes. Additionally, before I ran the commands
make clean
rm -rf tmp
make
I also had to add the last line below to FastJetFinder.h:
from line 152
TObjArray *fOutputArray; !
TObjArray *fRhoOutputArray; !
TObjArray *fConstArray; '''
Besides, in FastJetFinder.cc, I had to do
fConstArray->Add(constituent) instead of ->AddCandidate(constituent)
Now the code runs and seems to be doing what it's supposed to do! Thank you very much!
I don't think that it could be done by modifying the Delphes card.
Yes, it can be done by tweaking the FastJetFinder module.
There are many possible ways of doing it.
For example, you could try the following steps: