[Rivet] selecting soft leptons from HF decays

[David Bjergaard]

Hi Kristin,

The experts will probably have a better suggestion, (I'm just an
enthusiastic user).  

Is it possible to build up your own cluster-sequence along the lines of:
PseudoJets JetConstits
foreach(const Particle& p, cfs.particles()){
   //apply any selections  to particles
  JetConstits.push_back(fastjet::PseudoJet(p.momentum()));
}
fastjet::ClusterSequence cs(JetConstits,fastjet::JetDefinition(fastjet::cambridge_algorithm,0.8));
PseudoJets Jets = fastjet::sorted_by_pt(cs.inclusive_jets());

I've used this technique to apply jet clustering on different subsets of
the charged final state (cfs) particles to study the effect of
substructure variables at truth level.

Cheers,

    Dave

Kristin Lohwasser <kristin.lohwasser at cern.ch> writes:

> Dear all,
>
> I would like to come back to the following question (which needs to be
> resolved to be able to put a certain ATLAS analyses into Rivet):
>
> Is there any possibility to select the leptons (virtual W bosons) from
> semileptonic heavy flavour decays? And then feed these as veto
> particles into the Jet Finding?
>
> (Unfortunately, this data set has the feature that the soft leptons
> are removed from the jet finding), leading to noticeable differences
> (~5% level). [note: so far all further analyses using the data was
> based on the same jet definition, this is why we prefer to have an
> adapted rivet analysis ;) and not to re-publish the paper]
>
> As we would need the code also to extract the hadron level correction
> it would be soooo awesome to get hints on how to implement this
> feature.
>
> Thanks & Best
> Kristin