Plot dei genJet relativi ai tag. Il matching e' fatto guardando il quark primario del genJet e assicurandosi che corrisponda in flavor e deltaR con uno dei due quark tag generati.

JetFlavour jetFlavour = m_jfi.identifyBasedOnPartons(*(jet));
         int myflav = jetFlavour.flavour ();
         m_genJetFlavour -> push_back (myflav);
   if (m_evtFlag == 123 || m_evtFlag == 124)
   {
    TLorentzVector match;
    match.SetPxPyPzE (jetFlavour.vec4MainParton ().px (), jetFlavour.vec4MainParton ().py (),
          jetFlavour.vec4MainParton ().pz (), jetFlavour.vec4MainParton ().energy ()) ;
         
    double mydeltaR_F = m_genqTagF->DeltaR (match);
    double mydeltaE_F = m_genqTagF->E() - match.E() ;
    double mydeltaR_B = m_genqTagB->DeltaR (match);
    double mydeltaE_B = m_genqTagB->E() - match.E() ;
    
    if (mydeltaR_F < 0.9 && myflav == abs (m_genqTagF_Flavour))
    *m_genJetCloneTagF = match;
   
    else if (mydeltaR_B < 1.5 && myflav == abs (m_genqTagB_Flavour))
    *m_genJetCloneTagB = match;
         
   }