• Add dependency on:
    • #include "TtjetsClassifier/McEvent.h"
      #include "TtjetsClassifier/McSemiLeptonic.h"
      #include "TtjetsClassifier/GeneratorMatch.h"
      #include "TtjetsClassifier/PartonMatching.h"
      #include "TtjetsClassifier/ParticleIndex.h"
      #include "TLorentzVector.h"

1. MiniSLTruth::clear()

2. MiniSLTruth::initialize()

3. MiniSLTruth::writeMC()

virtual int writeMC(EventData *ed, SemiLeptonic *sl, bool pass_el, bool pass_mu, McEvent *mc, McSemiLeptonic *m_partSL, bool pass_part_el, bool pass_part_mu, float initWeight=-999.0);

same stuff from MiniSL::write()


two new calls:


MiniSLTruth::writeParticles(mc, m_partSL, pass_el, pass_part_el, pass_part_mu);

4. MiniSLTruth::writePartons(McEvent *mcEvent)

Based on TtjetsPartonTtree::write(McEvent mcEvent, McSemiLeptonic)

No significant changes added.

5. int MiniSLTruth::writeParticles(McEvent *mcEvent, McSemiLeptonic *mcSl, bool pass_el, bool pass_part_el, bool pass_part_mu)

Based on TtjetsPartTree::write(McEvent* mcEvent, McSemiLeptonic* mcSl, bool pass_part_el, bool pass_part_mu)

I also pass a bool pass_el, which tells if this is electron channel or not (this was not necessary in the original, because the other group writes electrons and muons in the same ntuple).

Notice that the original code has variables for el_ and mu_, but, as both channels are independent here, these variables are call lep_ and pass_el will decide which channel must be written when the call to writeParticles() is made.

This is it ! smile

Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r4 - 2013-09-08 - PatriciaFernandez
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback