higgs_csv.C

#define higgs_csv_cxx
#include <iostream>
#include "TMath.h"
#include "higgs_csv.h"
#include "treeChecker.h"

// Class for doing quick counting experiments on the treeChecker trees. This one checks a csv cut
// John Hakala, 5/11/2016

int higgs_csv::Loop(string workingpoint, float lowerMassBound, float upperMassBound)
{
   if (fChain == 0) return -1;
   Long64_t nentries = fChain->GetEntriesFast();
   Long64_t nbytes = 0, nb = 0;

   uint nHiggsJets = 0;
   for (Long64_t jentry=0; jentry<nentries;jentry++) {
      Long64_t ientry = LoadTree(jentry);
      if (ientry < 0) break;
      nb = fChain->GetEntry(jentry);   nbytes += nb;
      if (  Cut(ientry) < 0 
         || phJetInvMass_pruned_higgs < lowerMassBound 
         || phJetInvMass_pruned_higgs > upperMassBound ) continue;
      if (workingpoint == "nonenone" && leading < -999) { 
          cout << "The subjets csv values are: " << leading << ", " << subleading << endl;
          continue;
      }
      if (workingpoint == "loosenone" && leading < 0.605) {
         //cout << "The subjets csv values are: " << leading << ", " << subleading << endl;
         continue;
      }
      if (workingpoint == "mediumnone" && leading < 0.89) {
         //cout << "The subjets csv values are: " << leading << ", " << subleading << endl;
         continue;
      }
      if (workingpoint == "tightnone" && leading < 0.97) {
         //cout << "The subjets csv values are: " << leading << ", " << subleading << endl;
         continue;
      }
      if (workingpoint == "looseloose"
          && !loose_loose) continue;
      if (workingpoint == "mediumloose"
          && !medium_loose) continue;
      if (workingpoint == "tightloose"
          && !tight_loose) continue;
      if (workingpoint == "mediummedium"
          && !medium_medium) continue;
      if (workingpoint == "tightmedium"
          && !tight_medium) continue;
      if (workingpoint == "tighttight"
          && !tight_tight) continue;
      else if (   workingpoint != "nonenone"
               && workingpoint != "loosenone"
               && workingpoint != "mediumnone"
               && workingpoint != "tightnone"
               && workingpoint != "looseloose"
               && workingpoint != "mediumloose"
               && workingpoint != "tightloose"
               && workingpoint != "mediummedium"
               && workingpoint != "tightmedium"
               && workingpoint != "tighttight" 
              ) {
                  cout << "forgot to pick a valid csv working point!" << endl;
                  return 0;
                }
          
      ++nHiggsJets;
   }
   return nHiggsJets;
}

Int_t higgs_csv::Cut(Long64_t entry) {
  if (   cosThetaStar<-0.627*TMath::ATan((-0.005938*phJetInvMass_pruned_higgs)+3.427)
     &&  leadingPhAbsEta          < 1.4442
     &&  higgsJet_pruned_abseta   < 2.2
     &&  phJetDeltaR_higgs        > 1.1 ) return 1;

  else return -1;
}