Hi,
I am trying to reconstruct jets for my self-generated signal samples in which I require a Z boson to decay into a quark pair. I am using the durham kt algorithm as follows:
.Define(“RP_noMu”, “FCCAnalyses::ReconstructedParticle::remove(ReconstructedParticles, Leptons)”)
.Define("pseudo_jets_noMu", #Jet clustering
"FCCAnalyses::JetClusteringUtils::set_pseudoJets("
"ReconstructedParticle::get_px(RP_noMu),"
"ReconstructedParticle::get_py(RP_noMu),"
"ReconstructedParticle::get_pz(RP_noMu),"
"ReconstructedParticle::get_e (RP_noMu))")
.Define("clustered_durham2_noMu", "JetClustering::clustering_ee_kt(1, 20, 0, 0)(pseudo_jets_noMu)")
.Define("jets_durham2_noMu", "FCCAnalyses::JetClusteringUtils::get_pseudoJets(clustered_durham2_noMu)")
.Define("zjj_e", "FCCAnalyses::JetClusteringUtils::get_e (jets_durham2_noMu)")
.Define("zjj_px", "FCCAnalyses::JetClusteringUtils::get_px(jets_durham2_noMu)")
.Define("zjj_py", "FCCAnalyses::JetClusteringUtils::get_py(jets_durham2_noMu)")
.Define("zjj_pz", "FCCAnalyses::JetClusteringUtils::get_pz(jets_durham2_noMu)")
.Define("zjj_pt", "return sqrt(zjj_px*zjj_px + zjj_py*zjj_py);")
I would expect my signal events to have mostly 2 jets, because that is how I generated them, but this is what I observe:
I do not understand why the algorithm is reconstructing so many jets. Am I implementing this algorithm correctly, or should I use a different set of parameters?
Many thanks!