Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Work up until 74X MC samples #164

Closed
wants to merge 26 commits into from

Conversation

EmyrClement
Copy link
Contributor

This pull request contains a large amount of work, up until the point when 74X MC samples became available.

Top reconstruction

Doug has been working on the Rochester reco method. I've almost implemented the hit fit reco as done in Top-12-028.

74X MC

Mainly update sample names and cross sections, but also some changes for running and recognising TTJet samples.

Lepton pt threshold

I loosened the pt threshold at ntuple production, so we need to apply the tighter final selection here.

General tidy up/debugging

Everything else...

Remove obsolete function from MCParticle.h/.cc
Store B-quark information
Store single lepton information
Store neutrino information (GEN levels)
…bles that don't require and pseudo top recosntruction.
@@ -59,9 +60,11 @@ void Analysis::analyse() {
bjetWeights = BjetWeights(jets, numberOfBJets);

ttbar_plus_X_analyser_->analyse(currentEvent);
if ( ( currentEvent->getDataType() == DataType::TTJets || currentEvent->getDataType() == DataType::TT_Pythia8 ) && Globals::treePrefix_ == "" ) {
if ( ( currentEvent->getDataType() == DataType::TTJets_amcatnloFXFX || currentEvent->getDataType() == DataType::TTJets_madgraphMLM ) && Globals::treePrefix_ == "" ) {
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

To make this a bit easier in the future I created a feature request in issue #166

double ht = Event::HT(jets);
return ht + met->et() + lepton->pt();
double ht = 0;
if ( jets.size() > 0 ) ht = Event::HT(jets);
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If you are really concerned about having 0 jets, this should got into Event::HT and not here!

@kreczko
Copy link
Member

kreczko commented Jun 17, 2015

Except for the block of commented lines in the configs, I am happy with the PR.
There are some improvements necessary, but they are of architectural nature and will be followed up in separate issues

Before merging, could you please have a look at the travis errors? The code fails to compile.

return ht + lepton->pt();
else
return ht + met->et() + lepton->pt();

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The alternative is:

    // ST = HT + MET + lepton pt
    double ht = Event::HT(jets); #HT is already safe for 0 jets
    double MET = 0 ? met == 0 : met->et();
    double lpt = 0 ? lepton == 0 : lepton->pt();

    return ht + MET + lpt;

@kreczko
Copy link
Member

kreczko commented Jul 1, 2015

This should be now superseded by #169. @EmyrClement if you are happy with the changes, please merge #169 .

@kreczko kreczko closed this Jul 1, 2015
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants