MultivariateAnalysis/interface/TtLJetsAnalyzer.h

// -*- C++ -*-
//
// Package:    TtLJetsAnalyzer
// Class:      TtLJetsAnalyzer
// 
/**\class TtLJetsAnalyzer TtLJetsAnalyzer.cc RecoBTag/TtLJetsAnalyzer/src/TtLJetsAnalyzer.cc

 Description: Likelihood variables for ttbar events

 Implementation:
     Ttbar cross section without b tagging. Likelihood quantities.
*/
//
// Original Author:  Gena Kukartsev, kukarzev@fnal.gov
//         Created:  Fri Jun 20 06:19:59 CDT 2008
// $Id: TtLJetsAnalyzer.h,v 1.7 2009/02/24 23:33:21 kukartse Exp $
//
//

#include "LJMet/MultivariateAnalysis/interface/RooGKCounter.h"

#include <memory>
#include <string>
#include "TFile.h"
#include "TTree.h"
#include "Rtypes.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

class TtLJetsAnalyzer : public edm::EDAnalyzer {
public:
  explicit TtLJetsAnalyzer(const edm::ParameterSet&);
  ~TtLJetsAnalyzer();
  
  
private:

  virtual void beginJob(const edm::EventSetup&);
  virtual void analyze(const edm::Event&, const edm::EventSetup&);
  virtual void endJob(void);

  // _____ quantities read from the config file _________________________
  std::string _jetSource, _electronSource, _muonSource, _METSource;
  std::string _outFile;
  bool _is_mc;
  int nCaloJets_min;
  int nLepton_min;
  double jet_pt_min            ;
  double jet_eta_max           ;
  double muon_pt_min           ;
  double muon_eta_max          ;
  double muon_trackIso_max     ;
  double muon_caloIso_max      ;
  double electron_pt_min       ;
  double electron_eta_max      ;
  double electron_trackIso_max ;
  double electron_caloIso_max  ;
  double met_et_min            ;
  double min_dr_jet_lepton     ;

  RooGKCounter eventCounter;
  RooGKCounter selectedEvents;

  TFile * _file;
  TTree * _tree;

  //
  // _____ output tree variables ________________________________________
  Int_t _event;
  Int_t _process_id;
  Int_t _n_jets;
  Int_t _n_removed_jets;
  Int_t _n_met;
  Int_t _n_muons;
  Int_t _n_electrons;
  Int_t _lepton_is_muon;
  Double_t _event_weight;
  Double_t _lepton_et;
  Double_t _lepton_pt;
  Double_t _lepton_eta;
  Double_t _lepton_phi;
  Double_t _lepton_energy;
  Double_t _muon_track_chi2;
  Double_t _muon_track_ndof;
  Double_t _muon_track_nChi2;
  Double_t _muon_outerTrack_chi2;
  Double_t _muon_outerTrack_ndof;
  Double_t _muon_outerTrack_nChi2;
  Double_t _muon_globalTrack_chi2;
  Double_t _muon_globalTrack_ndof;
  Double_t _muon_globalTrack_nChi2;
  Double_t _lepton_track_iso;
  Double_t _lepton_calo_iso;
  Double_t _lepton_ecal_iso;
  Double_t _lepton_hcal_iso;
  Double_t _lepton_user_iso;
  Double_t _lepton_IsoDeposit_em;
  Double_t _lepton_IsoDeposit_had;
  Double_t _met_et;
  Double_t _met_pt;
  Double_t _met_eta;
  Double_t _met_phi;
  Double_t _met_energy;
  Double_t _jet1_et;
  Double_t _jet1_pt;
  Double_t _jet1_eta;
  Double_t _jet1_phi;
  Double_t _jet1_energy;
  Double_t _jet2_et;
  Double_t _jet2_pt;
  Double_t _jet2_eta;
  Double_t _jet2_phi;
  Double_t _jet2_energy;
  Double_t _jet3_et;
  Double_t _jet3_pt;
  Double_t _jet3_eta;
  Double_t _jet3_phi;
  Double_t _jet3_energy;
  Double_t _jet4_et;
  Double_t _jet4_pt;
  Double_t _jet4_eta;
  Double_t _jet4_phi;
  Double_t _jet4_energy;
  Double_t _aplanarity       ;
  Double_t _centrality       ;
  Double_t _sphericity       ;
  Double_t _ht               ;
  Double_t _htPlusLepton     ;
  Double_t _metHtPlusLepton  ;
  Double_t _h                ;
  Double_t _ktMinPrime       ;
  Double_t _dPhiLMet         ;
  Double_t _minDijetMass     ;
  Double_t _maxJetEta        ;
  Double_t _et3              ;
  Double_t _minDiJetDeltaR   ;
  Double_t _leptonJetDeltaR  ;
  Double_t _ht2p             ;
  Double_t _jet1Jet2DeltaR   ;
  Double_t _jet1Jet2DeltaPhi ;
  Double_t _jet1Jet2_M       ;
  Double_t _jet1Jet2_Pt      ;
  Double_t _jet1Jet2W_M      ;
  Double_t _jet1Jet2W_Pt     ;
  Double_t _hz               ;
  Double_t _HT2              ;
  Double_t _HT2prime         ;
  Double_t _W_MT             ;
  Double_t _W_M              ;
  Double_t _W_Pt             ;
  Double_t _DphiJMET         ;
  // Ht               
  Double_t _getHt            ;
  Double_t _getHtp           ;
  Double_t _getHtpp          ;
  Double_t _getHt2           ;
  Double_t _getHt2p          ;
  Double_t _getHt2pp         ;
  Double_t _getHt3           ;
  Double_t _getHt3p          ;
  Double_t _getHt3pp         ;
  Double_t _getCen           ;
  Double_t _getNJW           ;
  Double_t _getJetEtaMax     ;
  Double_t _getMdijetMin     ;
  Double_t _getMtjets        ;
  Double_t _getSqrtsT        ;
  Double_t _getMtAurelio     ;
  Double_t _getPzOverHT      ;
  Double_t _getMevent        ;
  Double_t _getM123inv       ;
  Double_t _getEta2Sum       ;
  Double_t _getMwRec         ;
  Double_t _getH             ;
  // event topo       
  Double_t _getSph           ;
  Double_t _getApl           ;
  Double_t _getAplMu         ;
  // Kt               
  Double_t _getKtminp        ;
  Double_t _getKtminpReduced ;
  Double_t _getDrMinJetJet   ;
  // mT               
  Double_t _getDphiMuMet     ;
  Double_t _getMt            ;
  //
  //
  //_____ b tagging ______________________________________________
  vector<double> * _jet_pt;
  vector<double> * _jet_pz;
  vector<double> * _jet_eta;
  vector<double> * _jet_phi;
  vector<double> * _jet_et;
  vector<double> * _jet_energy;
  vector<int> *    _jet_flavor;
  vector<double> * _lept_pt;
  vector<double> * _lept_pz;
  vector<double> * _lept_eta;
  vector<double> * _lept_phi;
  vector<double> * _lept_et;
  vector<double> * _lept_energy;
  vector<int> *    _lept_flavor;
  /* FIXME: add vectors of these:
  Double_t _lepton_track_iso;
  Double_t _lepton_calo_iso;
  Double_t _lepton_ecal_iso;
  Double_t _lepton_hcal_iso;
  Double_t _lepton_user_iso;
  Double_t _lepton_IsoDeposit_em;
  Double_t _lepton_IsoDeposit_had;
  */
  TBranch * b_jet_pt;
  TBranch * b_jet_pz;
  TBranch * b_jet_eta;
  TBranch * b_jet_phi;
  TBranch * b_jet_et;
  TBranch * b_jet_energy;
  TBranch * b_jet_flavor;
  TBranch * b_lept_pt;
  TBranch * b_lept_pz;
  TBranch * b_lept_eta;
  TBranch * b_lept_phi;
  TBranch * b_lept_et;
  TBranch * b_lept_energy;
  TBranch * b_lept_flavor;
  vector<string> * _b_tagger_label;
  vector<float> * _b_tagger_discr;
  TBranch * b_b_tagger_label;
  TBranch * b_b_tagger_discr;
  Int_t _n_tagged_jets_trackCounting_loose;
  Int_t _n_tagged_jets_trackCounting_medium;
  Int_t _n_tagged_jets_trackCounting_tight;
  Int_t _n_tagged_jets_jetProb_loose;
  Int_t _n_tagged_jets_jetProb_medium;
  Int_t _n_tagged_jets_jetProb_tight;

  //
  //
  //_____ electron quality _______________________________________
  Float_t _eidLoose;
  Float_t _eidRobustHighEnergy;
  Float_t _eidRobustLoose;
  Float_t _eidRobustTight;
  Float_t _eidTight;
  Int_t    _GsfElectron_classification;
  Double_t _GsfElectron_hadronicOverEm;
  Double_t _GsfElectron_caloEnergyError;
  Double_t _GsfElectron_trackMomentumError;
  Int_t    _GsfElectron_numberOfClusters;
  Double_t _GsfElectron_caloEnergy;  // the super cluster energy corrected by EnergyScaleFactor
  Double_t _GsfElectron_eSuperClusterOverP;  // the supercluster energy / track momentum at impact point
  Double_t _GsfElectron_eSeedClusterOverPout;  // the seed cluster energy / track momentum at calo from outermost state
  Double_t _GsfElectron_deltaEtaSuperClusterTrackAtVtx;  // the supercluster eta - track eta from helix extrapolation from impact point
  Double_t _GsfElectron_deltaEtaSeedClusterTrackAtCalo;  // the seed cluster eta - track eta at calo from outermost state
  Double_t _GsfElectron_deltaPhiSuperClusterTrackAtVtx;  // the supercluster phi - track phi from helix extrapolation from impact point
  Double_t _GsfElectron_deltaPhiSeedClusterTrackAtCalo;  // the seed cluster phi - track phi at calo from outermost state
  // ____ branches
  TBranch * b_eidLoose;
  TBranch * b_eidRobustHighEnergy;
  TBranch * b_eidRobustLoose;
  TBranch * b_eidRobustTight;
  TBranch * b_eidTight;
  TBranch * b_GsfElectron_classification;
  TBranch * b_GsfElectron_hadronicOverEm;
  TBranch * b_GsfElectron_caloEnergyError;
  TBranch * b_GsfElectron_trackMomentumError;
  TBranch * b_GsfElectron_numberOfClusters;
  TBranch * b_GsfElectron_caloEnergy;  // the super cluster energy corrected by EnergyScaleFactor
  TBranch * b_GsfElectron_eSuperClusterOverP;  // the supercluster energy / track momentum at impact point
  TBranch * b_GsfElectron_eSeedClusterOverPout;  // the seed cluster energy / track momentum at calo from outermost state
  TBranch * b_GsfElectron_deltaEtaSuperClusterTrackAtVtx;  // the supercluster eta - track eta from helix extrapolation from impact point
  TBranch * b_GsfElectron_deltaEtaSeedClusterTrackAtCalo;  // the seed cluster eta - track eta at calo from outermost state
  TBranch * b_GsfElectron_deltaPhiSuperClusterTrackAtVtx;  // the supercluster phi - track phi from helix extrapolation from impact point
  TBranch * b_GsfElectron_deltaPhiSeedClusterTrackAtCalo;  // the seed cluster phi - track phi at calo from outermost state
  //
  //
  //_____ muon quality _______________________________________
  bool _muon_isGood_All;                      // dummy options - always true
  bool _muon_isGood_AllGlobalMuons;           // checks isGlobalMuon flag
  bool _muon_isGood_AllStandAloneMuons;       // checks isStandAloneMuon flag
  bool _muon_isGood_AllTrackerMuons;          // checks isTrackerMuon flag
  bool _muon_isGood_TrackerMuonArbitrated;    // resolve ambiguity of sharing segments
  bool _muon_isGood_AllArbitrated;            // all muons with the tracker muon arbitrated
  bool _muon_isGood_GlobalMuonPromptTight;    // global muons with tighter fit requirements
  bool _muon_isGood_TMLastStationLoose;       // penetration depth loose selector
  bool _muon_isGood_TMLastStationTight;       // penetration depth tight selector
  bool _muon_isGood_TM2DCompatibilityLoose;   // likelihood based loose selector
  bool _muon_isGood_TM2DCompatibilityTight;   // likelihood based tight selector
  bool _muon_isGood_TMOneStationLoose;        // require one well matched segment
  bool _muon_isGood_TMOneStationTight;        // require one well matched segment
  bool _muon_isGood_TMLastStationOptimizedLowPtLoose; // combination of TMLastStation and TMOneStation
  bool _muon_isGood_TMLastStationOptimizedLowPtTight;  // combination of TMLastStation and TMOneStation
  //____ branches
  TBranch * b_muon_isGood_All;                      // dummy options - always true
  TBranch * b_muon_isGood_AllGlobalMuons;           // checks isGlobalMuon flag
  TBranch * b_muon_isGood_AllStandAloneMuons;       // checks isStandAloneMuon flag
  TBranch * b_muon_isGood_AllTrackerMuons;          // checks isTrackerMuon flag
  TBranch * b_muon_isGood_TrackerMuonArbitrated;    // resolve ambiguity of sharing segments
  TBranch * b_muon_isGood_AllArbitrated;            // all muons with the tracker muon arbitrated
  TBranch * b_muon_isGood_GlobalMuonPromptTight;    // global muons with tighter fit requirements
  TBranch * b_muon_isGood_TMLastStationLoose;       // penetration depth loose selector
  TBranch * b_muon_isGood_TMLastStationTight;       // penetration depth tight selector
  TBranch * b_muon_isGood_TM2DCompatibilityLoose;   // likelihood based loose selector
  TBranch * b_muon_isGood_TM2DCompatibilityTight;   // likelihood based tight selector
  TBranch * b_muon_isGood_TMOneStationLoose;        // require one well matched segment
  TBranch * b_muon_isGood_TMOneStationTight;        // require one well matched segment
  TBranch * b_muon_isGood_TMLastStationOptimizedLowPtLoose; // combination of TMLastStation and TMOneStation
  TBranch * b_muon_isGood_TMLastStationOptimizedLowPtTight;  // combination of TMLastStation and TMOneStation
  //____ output tree variables
  TBranch * b_event;
  TBranch * b_process_id;
  TBranch * b_n_jets;
  TBranch * b_n_removed_jets;
  TBranch * b_n_met;
  TBranch * b_n_muons;
  TBranch * b_n_electrons;
  TBranch * b_lepton_is_muon;
  TBranch * b_event_weight;
  TBranch * b_lepton_et;
  TBranch * b_lepton_pt;
  TBranch * b_lepton_eta;
  TBranch * b_lepton_phi;
  TBranch * b_lepton_energy;
  TBranch * b_muon_track_chi2;
  TBranch * b_muon_track_ndof;
  TBranch * b_muon_track_nChi2;
  TBranch * b_muon_outerTrack_chi2;
  TBranch * b_muon_outerTrack_ndof;
  TBranch * b_muon_outerTrack_nChi2;
  TBranch * b_muon_globalTrack_chi2;
  TBranch * b_muon_globalTrack_ndof;
  TBranch * b_muon_globalTrack_nChi2;
  TBranch * b_lepton_track_iso;
  TBranch * b_lepton_calo_iso;
  TBranch * b_lepton_ecal_iso;
  TBranch * b_lepton_hcal_iso;
  TBranch * b_lepton_user_iso;
  TBranch * b_lepton_IsoDeposit_em;
  TBranch * b_lepton_IsoDeposit_had;
  TBranch * b_met_et;
  TBranch * b_met_pt;
  TBranch * b_met_eta;
  TBranch * b_met_phi;
  TBranch * b_met_energy;
  TBranch * b_jet1_et;
  TBranch * b_jet1_pt;
  TBranch * b_jet1_eta;
  TBranch * b_jet1_phi;
  TBranch * b_jet1_energy;
  TBranch * b_jet2_et;
  TBranch * b_jet2_pt;
  TBranch * b_jet2_eta;
  TBranch * b_jet2_phi;
  TBranch * b_jet2_energy;
  TBranch * b_jet3_et;
  TBranch * b_jet3_pt;
  TBranch * b_jet3_eta;
  TBranch * b_jet3_phi;
  TBranch * b_jet3_energy;
  TBranch * b_jet4_et;
  TBranch * b_jet4_pt;
  TBranch * b_jet4_eta;
  TBranch * b_jet4_phi;
  TBranch * b_jet4_energy;
  TBranch * b_n_tagged_jets_trackCounting_loose;
  TBranch * b_n_tagged_jets_trackCounting_medium;
  TBranch * b_n_tagged_jets_trackCounting_tight;
  TBranch * b_n_tagged_jets_jetProb_loose;
  TBranch * b_n_tagged_jets_jetProb_medium;
  TBranch * b_n_tagged_jets_jetProb_tight;       
  // ____
  TBranch * b_aplanarity       ;
  TBranch * b_centrality       ;
  TBranch * b_sphericity       ;
  TBranch * b_ht               ;
  TBranch * b_htPlusLepton     ;
  TBranch * b_metHtPlusLepton  ;
  TBranch * b_h                ;
  TBranch * b_ktMinPrime       ;
  TBranch * b_dPhiLMet         ;
  TBranch * b_minDijetMass     ;
  TBranch * b_maxJetEta        ;
  TBranch * b_et3              ;
  TBranch * b_minDiJetDeltaR   ;
  TBranch * b_leptonJetDeltaR  ;
  TBranch * b_ht2p             ;
  TBranch * b_jet1Jet2DeltaR   ;
  TBranch * b_jet1Jet2DeltaPhi ;
  TBranch * b_jet1Jet2_M       ;
  TBranch * b_jet1Jet2_Pt      ;
  TBranch * b_jet1Jet2W_M      ;
  TBranch * b_jet1Jet2W_Pt     ;
  TBranch * b_hz               ;
  TBranch * b_HT2;
  TBranch * b_HT2prime;
  TBranch * b_W_MT;
  TBranch * b_W_M              ;
  TBranch * b_W_Pt             ;
  TBranch * b_DphiJMET         ;
  // Ht               
  TBranch * b_getHt            ;
  TBranch * b_getHtp           ;
  TBranch * b_getHtpp          ;
  TBranch * b_getHt2           ;
  TBranch * b_getHt2p          ;
  TBranch * b_getHt2pp         ;
  TBranch * b_getHt3           ;
  TBranch * b_getHt3p          ;
  TBranch * b_getHt3pp         ;
  TBranch * b_getCen           ;
  TBranch * b_getNJW           ;
  TBranch * b_getJetEtaMax     ;
  TBranch * b_getMdijetMin     ;
  TBranch * b_getMtjets        ;
  TBranch * b_getSqrtsT        ;
  TBranch * b_getMtAurelio     ;
  TBranch * b_getPzOverHT      ;
  TBranch * b_getMevent        ;
  TBranch * b_getM123inv       ;
  TBranch * b_getEta2Sum       ;
  TBranch * b_getMwRec         ;
  TBranch * b_getH             ;
  // event topo       
  TBranch * b_getSph           ;
  TBranch * b_getApl           ;
  TBranch * b_getAplMu         ;
  // Kt               
  TBranch * b_getKtminp        ;
  TBranch * b_getKtminpReduced ;
  TBranch * b_getDrMinJetJet   ;
  // mT               
  TBranch * b_getDphiMuMet     ;
  TBranch * b_getMt            ;

  
};
Revision:	1.8
Committed:	Thu Mar 5 17:30:20 2009 UTC (16 years, 2 months ago) by kukartse
Content type:	text/plain
Branch:	MAIN
CVS Tags:	V00-01-10, gak031009, gak030509
Changes since 1.7:	+49 -7 lines
Log Message:	* empty log message *
#	User	Rev	Content
1	kukartse	1.1	// -- C++ --
2			//
3			// Package: TtLJetsAnalyzer
4			// Class: TtLJetsAnalyzer
5			//
6			/**\class TtLJetsAnalyzer TtLJetsAnalyzer.cc RecoBTag/TtLJetsAnalyzer/src/TtLJetsAnalyzer.cc
7
8			Description: Likelihood variables for ttbar events
9
10			Implementation:
11			Ttbar cross section without b tagging. Likelihood quantities.
12			*/
13			//
14			// Original Author: Gena Kukartsev, kukarzev@fnal.gov
15			// Created: Fri Jun 20 06:19:59 CDT 2008
16	kukartse	1.8	// $Id: TtLJetsAnalyzer.h,v 1.7 2009/02/24 23:33:21 kukartse Exp $
17	kukartse	1.1	//
18			//
19
20	kukartse	1.8	#include "LJMet/MultivariateAnalysis/interface/RooGKCounter.h"
21
22	kukartse	1.1	#include <memory>
23			#include <string>
24	kukartse	1.8	#include "TFile.h"
25			#include "TTree.h"
26			#include "Rtypes.h"
27	kukartse	1.1	#include "FWCore/Framework/interface/Frameworkfwd.h"
28			#include "FWCore/Framework/interface/EDAnalyzer.h"
29			#include "FWCore/Framework/interface/Event.h"
30			#include "FWCore/Framework/interface/MakerMacros.h"
31			#include "FWCore/ParameterSet/interface/ParameterSet.h"
32
33			class TtLJetsAnalyzer : public edm::EDAnalyzer {
34			public:
35			explicit TtLJetsAnalyzer(const edm::ParameterSet&);
36			~TtLJetsAnalyzer();
37
38
39			private:
40	kukartse	1.4
41			virtual void beginJob(const edm::EventSetup&);
42			virtual void analyze(const edm::Event&, const edm::EventSetup&);
43			virtual void endJob(void);
44
45	kukartse	1.5	// _____ quantities read from the config file _________________________
46	kukartse	1.1	std::string _jetSource, _electronSource, _muonSource, _METSource;
47			std::string _outFile;
48	kukartse	1.8	bool _is_mc;
49	kukartse	1.1	int nCaloJets_min;
50			int nLepton_min;
51	jindal	1.6	double jet_pt_min ;
52	kukartse	1.1	double jet_eta_max ;
53			double muon_pt_min ;
54			double muon_eta_max ;
55			double muon_trackIso_max ;
56			double muon_caloIso_max ;
57			double electron_pt_min ;
58			double electron_eta_max ;
59			double electron_trackIso_max ;
60			double electron_caloIso_max ;
61			double met_et_min ;
62	kukartse	1.7	double min_dr_jet_lepton ;
63	kukartse	1.1
64			RooGKCounter eventCounter;
65			RooGKCounter selectedEvents;
66
67			TFile * _file;
68			TTree * _tree;
69
70	kukartse	1.5	//
71			// _____ output tree variables ________________________________________
72	kukartse	1.1	Int_t _event;
73			Int_t _process_id;
74			Int_t _n_jets;
75	kukartse	1.7	Int_t _n_removed_jets;
76	kukartse	1.1	Int_t _n_met;
77			Int_t _n_muons;
78			Int_t _n_electrons;
79			Int_t _lepton_is_muon;
80			Double_t _event_weight;
81	kukartse	1.5	Double_t _lepton_et;
82	kukartse	1.1	Double_t _lepton_pt;
83			Double_t _lepton_eta;
84			Double_t _lepton_phi;
85			Double_t _lepton_energy;
86	kukartse	1.3	Double_t _muon_track_chi2;
87			Double_t _muon_track_ndof;
88			Double_t _muon_track_nChi2;
89			Double_t _muon_outerTrack_chi2;
90			Double_t _muon_outerTrack_ndof;
91			Double_t _muon_outerTrack_nChi2;
92			Double_t _muon_globalTrack_chi2;
93			Double_t _muon_globalTrack_ndof;
94			Double_t _muon_globalTrack_nChi2;
95	kukartse	1.1	Double_t _lepton_track_iso;
96			Double_t _lepton_calo_iso;
97	kukartse	1.2	Double_t _lepton_ecal_iso;
98			Double_t _lepton_hcal_iso;
99			Double_t _lepton_user_iso;
100	kukartse	1.8	Double_t _lepton_IsoDeposit_em;
101			Double_t _lepton_IsoDeposit_had;
102	kukartse	1.5	Double_t _met_et;
103	kukartse	1.1	Double_t _met_pt;
104			Double_t _met_eta;
105			Double_t _met_phi;
106			Double_t _met_energy;
107	kukartse	1.5	Double_t _jet1_et;
108	kukartse	1.1	Double_t _jet1_pt;
109			Double_t _jet1_eta;
110			Double_t _jet1_phi;
111			Double_t _jet1_energy;
112	kukartse	1.5	Double_t _jet2_et;
113	kukartse	1.1	Double_t _jet2_pt;
114			Double_t _jet2_eta;
115			Double_t _jet2_phi;
116			Double_t _jet2_energy;
117	kukartse	1.5	Double_t _jet3_et;
118	kukartse	1.1	Double_t _jet3_pt;
119			Double_t _jet3_eta;
120			Double_t _jet3_phi;
121			Double_t _jet3_energy;
122	kukartse	1.5	Double_t _jet4_et;
123	kukartse	1.1	Double_t _jet4_pt;
124			Double_t _jet4_eta;
125			Double_t _jet4_phi;
126			Double_t _jet4_energy;
127			Double_t _aplanarity ;
128			Double_t _centrality ;
129			Double_t _sphericity ;
130			Double_t _ht ;
131			Double_t _htPlusLepton ;
132			Double_t _metHtPlusLepton ;
133			Double_t _h ;
134			Double_t _ktMinPrime ;
135			Double_t _dPhiLMet ;
136			Double_t _minDijetMass ;
137			Double_t _maxJetEta ;
138			Double_t _et3 ;
139			Double_t _minDiJetDeltaR ;
140			Double_t _leptonJetDeltaR ;
141			Double_t _ht2p ;
142			Double_t _jet1Jet2DeltaR ;
143			Double_t _jet1Jet2DeltaPhi ;
144			Double_t _jet1Jet2_M ;
145			Double_t _jet1Jet2_Pt ;
146			Double_t _jet1Jet2W_M ;
147			Double_t _jet1Jet2W_Pt ;
148			Double_t _hz ;
149			Double_t _HT2 ;
150			Double_t _HT2prime ;
151			Double_t _W_MT ;
152			Double_t _W_M ;
153			Double_t _W_Pt ;
154			Double_t _DphiJMET ;
155			// Ht
156			Double_t _getHt ;
157			Double_t _getHtp ;
158			Double_t _getHtpp ;
159			Double_t _getHt2 ;
160			Double_t _getHt2p ;
161			Double_t _getHt2pp ;
162			Double_t _getHt3 ;
163			Double_t _getHt3p ;
164			Double_t _getHt3pp ;
165			Double_t _getCen ;
166			Double_t _getNJW ;
167			Double_t _getJetEtaMax ;
168			Double_t _getMdijetMin ;
169			Double_t _getMtjets ;
170			Double_t _getSqrtsT ;
171			Double_t _getMtAurelio ;
172			Double_t _getPzOverHT ;
173			Double_t _getMevent ;
174			Double_t _getM123inv ;
175			Double_t _getEta2Sum ;
176			Double_t _getMwRec ;
177			Double_t _getH ;
178			// event topo
179			Double_t _getSph ;
180			Double_t _getApl ;
181			Double_t _getAplMu ;
182			// Kt
183			Double_t _getKtminp ;
184			Double_t _getKtminpReduced ;
185			Double_t _getDrMinJetJet ;
186			// mT
187			Double_t _getDphiMuMet ;
188			Double_t _getMt ;
189	kukartse	1.5	//
190			//
191			//_____ b tagging ______________________________________________
192	kukartse	1.8	vector<double> * _jet_pt;
193			vector<double> * _jet_pz;
194			vector<double> * _jet_eta;
195			vector<double> * _jet_phi;
196			vector<double> * _jet_et;
197			vector<double> * _jet_energy;
198			vector<int> * _jet_flavor;
199			vector<double> * _lept_pt;
200			vector<double> * _lept_pz;
201			vector<double> * _lept_eta;
202			vector<double> * _lept_phi;
203			vector<double> * _lept_et;
204			vector<double> * _lept_energy;
205			vector<int> * _lept_flavor;
206			/* FIXME: add vectors of these:
207			Double_t _lepton_track_iso;
208			Double_t _lepton_calo_iso;
209			Double_t _lepton_ecal_iso;
210			Double_t _lepton_hcal_iso;
211			Double_t _lepton_user_iso;
212			Double_t _lepton_IsoDeposit_em;
213			Double_t _lepton_IsoDeposit_had;
214			*/
215			TBranch * b_jet_pt;
216			TBranch * b_jet_pz;
217			TBranch * b_jet_eta;
218			TBranch * b_jet_phi;
219			TBranch * b_jet_et;
220			TBranch * b_jet_energy;
221			TBranch * b_jet_flavor;
222			TBranch * b_lept_pt;
223			TBranch * b_lept_pz;
224			TBranch * b_lept_eta;
225			TBranch * b_lept_phi;
226			TBranch * b_lept_et;
227			TBranch * b_lept_energy;
228			TBranch * b_lept_flavor;
229	kukartse	1.5	vector<string> * _b_tagger_label;
230			vector<float> * _b_tagger_discr;
231			TBranch * b_b_tagger_label;
232			TBranch * b_b_tagger_discr;
233	kukartse	1.4	Int_t _n_tagged_jets_trackCounting_loose;
234			Int_t _n_tagged_jets_trackCounting_medium;
235			Int_t _n_tagged_jets_trackCounting_tight;
236			Int_t _n_tagged_jets_jetProb_loose;
237			Int_t _n_tagged_jets_jetProb_medium;
238			Int_t _n_tagged_jets_jetProb_tight;
239	kukartse	1.1
240	kukartse	1.4	//
241			//
242			//_____ electron quality _______________________________________
243	kukartse	1.5	Float_t _eidLoose;
244			Float_t _eidRobustHighEnergy;
245			Float_t _eidRobustLoose;
246			Float_t _eidRobustTight;
247			Float_t _eidTight;
248	kukartse	1.4	Int_t _GsfElectron_classification;
249			Double_t _GsfElectron_hadronicOverEm;
250			Double_t _GsfElectron_caloEnergyError;
251			Double_t _GsfElectron_trackMomentumError;
252			Int_t _GsfElectron_numberOfClusters;
253			Double_t _GsfElectron_caloEnergy; // the super cluster energy corrected by EnergyScaleFactor
254			Double_t _GsfElectron_eSuperClusterOverP; // the supercluster energy / track momentum at impact point
255			Double_t _GsfElectron_eSeedClusterOverPout; // the seed cluster energy / track momentum at calo from outermost state
256			Double_t _GsfElectron_deltaEtaSuperClusterTrackAtVtx; // the supercluster eta - track eta from helix extrapolation from impact point
257			Double_t _GsfElectron_deltaEtaSeedClusterTrackAtCalo; // the seed cluster eta - track eta at calo from outermost state
258			Double_t _GsfElectron_deltaPhiSuperClusterTrackAtVtx; // the supercluster phi - track phi from helix extrapolation from impact point
259			Double_t _GsfElectron_deltaPhiSeedClusterTrackAtCalo; // the seed cluster phi - track phi at calo from outermost state
260			// ____ branches
261	kukartse	1.5	TBranch * b_eidLoose;
262			TBranch * b_eidRobustHighEnergy;
263			TBranch * b_eidRobustLoose;
264			TBranch * b_eidRobustTight;
265			TBranch * b_eidTight;
266	kukartse	1.4	TBranch * b_GsfElectron_classification;
267			TBranch * b_GsfElectron_hadronicOverEm;
268			TBranch * b_GsfElectron_caloEnergyError;
269			TBranch * b_GsfElectron_trackMomentumError;
270			TBranch * b_GsfElectron_numberOfClusters;
271			TBranch * b_GsfElectron_caloEnergy; // the super cluster energy corrected by EnergyScaleFactor
272			TBranch * b_GsfElectron_eSuperClusterOverP; // the supercluster energy / track momentum at impact point
273			TBranch * b_GsfElectron_eSeedClusterOverPout; // the seed cluster energy / track momentum at calo from outermost state
274			TBranch * b_GsfElectron_deltaEtaSuperClusterTrackAtVtx; // the supercluster eta - track eta from helix extrapolation from impact point
275			TBranch * b_GsfElectron_deltaEtaSeedClusterTrackAtCalo; // the seed cluster eta - track eta at calo from outermost state
276			TBranch * b_GsfElectron_deltaPhiSuperClusterTrackAtVtx; // the supercluster phi - track phi from helix extrapolation from impact point
277			TBranch * b_GsfElectron_deltaPhiSeedClusterTrackAtCalo; // the seed cluster phi - track phi at calo from outermost state
278			//
279			//
280			//_____ muon quality _______________________________________
281			bool _muon_isGood_All; // dummy options - always true
282			bool _muon_isGood_AllGlobalMuons; // checks isGlobalMuon flag
283			bool _muon_isGood_AllStandAloneMuons; // checks isStandAloneMuon flag
284			bool _muon_isGood_AllTrackerMuons; // checks isTrackerMuon flag
285			bool _muon_isGood_TrackerMuonArbitrated; // resolve ambiguity of sharing segments
286			bool _muon_isGood_AllArbitrated; // all muons with the tracker muon arbitrated
287			bool _muon_isGood_GlobalMuonPromptTight; // global muons with tighter fit requirements
288			bool _muon_isGood_TMLastStationLoose; // penetration depth loose selector
289			bool _muon_isGood_TMLastStationTight; // penetration depth tight selector
290			bool _muon_isGood_TM2DCompatibilityLoose; // likelihood based loose selector
291			bool _muon_isGood_TM2DCompatibilityTight; // likelihood based tight selector
292			bool _muon_isGood_TMOneStationLoose; // require one well matched segment
293			bool _muon_isGood_TMOneStationTight; // require one well matched segment
294			bool _muon_isGood_TMLastStationOptimizedLowPtLoose; // combination of TMLastStation and TMOneStation
295			bool _muon_isGood_TMLastStationOptimizedLowPtTight; // combination of TMLastStation and TMOneStation
296			//____ branches
297			TBranch * b_muon_isGood_All; // dummy options - always true
298			TBranch * b_muon_isGood_AllGlobalMuons; // checks isGlobalMuon flag
299			TBranch * b_muon_isGood_AllStandAloneMuons; // checks isStandAloneMuon flag
300			TBranch * b_muon_isGood_AllTrackerMuons; // checks isTrackerMuon flag
301			TBranch * b_muon_isGood_TrackerMuonArbitrated; // resolve ambiguity of sharing segments
302			TBranch * b_muon_isGood_AllArbitrated; // all muons with the tracker muon arbitrated
303			TBranch * b_muon_isGood_GlobalMuonPromptTight; // global muons with tighter fit requirements
304			TBranch * b_muon_isGood_TMLastStationLoose; // penetration depth loose selector
305			TBranch * b_muon_isGood_TMLastStationTight; // penetration depth tight selector
306			TBranch * b_muon_isGood_TM2DCompatibilityLoose; // likelihood based loose selector
307			TBranch * b_muon_isGood_TM2DCompatibilityTight; // likelihood based tight selector
308			TBranch * b_muon_isGood_TMOneStationLoose; // require one well matched segment
309			TBranch * b_muon_isGood_TMOneStationTight; // require one well matched segment
310			TBranch * b_muon_isGood_TMLastStationOptimizedLowPtLoose; // combination of TMLastStation and TMOneStation
311			TBranch * b_muon_isGood_TMLastStationOptimizedLowPtTight; // combination of TMLastStation and TMOneStation
312			//____ output tree variables
313	kukartse	1.1	TBranch * b_event;
314			TBranch * b_process_id;
315			TBranch * b_n_jets;
316	kukartse	1.8	TBranch * b_n_removed_jets;
317	kukartse	1.1	TBranch * b_n_met;
318			TBranch * b_n_muons;
319			TBranch * b_n_electrons;
320			TBranch * b_lepton_is_muon;
321			TBranch * b_event_weight;
322	kukartse	1.5	TBranch * b_lepton_et;
323	kukartse	1.1	TBranch * b_lepton_pt;
324			TBranch * b_lepton_eta;
325			TBranch * b_lepton_phi;
326			TBranch * b_lepton_energy;
327	kukartse	1.3	TBranch * b_muon_track_chi2;
328			TBranch * b_muon_track_ndof;
329			TBranch * b_muon_track_nChi2;
330			TBranch * b_muon_outerTrack_chi2;
331			TBranch * b_muon_outerTrack_ndof;
332			TBranch * b_muon_outerTrack_nChi2;
333			TBranch * b_muon_globalTrack_chi2;
334			TBranch * b_muon_globalTrack_ndof;
335			TBranch * b_muon_globalTrack_nChi2;
336	kukartse	1.1	TBranch * b_lepton_track_iso;
337			TBranch * b_lepton_calo_iso;
338	kukartse	1.2	TBranch * b_lepton_ecal_iso;
339			TBranch * b_lepton_hcal_iso;
340			TBranch * b_lepton_user_iso;
341	kukartse	1.8	TBranch * b_lepton_IsoDeposit_em;
342			TBranch * b_lepton_IsoDeposit_had;
343	kukartse	1.5	TBranch * b_met_et;
344	kukartse	1.1	TBranch * b_met_pt;
345			TBranch * b_met_eta;
346			TBranch * b_met_phi;
347			TBranch * b_met_energy;
348	kukartse	1.5	TBranch * b_jet1_et;
349	kukartse	1.1	TBranch * b_jet1_pt;
350			TBranch * b_jet1_eta;
351			TBranch * b_jet1_phi;
352			TBranch * b_jet1_energy;
353	kukartse	1.5	TBranch * b_jet2_et;
354	kukartse	1.1	TBranch * b_jet2_pt;
355			TBranch * b_jet2_eta;
356			TBranch * b_jet2_phi;
357			TBranch * b_jet2_energy;
358	kukartse	1.5	TBranch * b_jet3_et;
359	kukartse	1.1	TBranch * b_jet3_pt;
360			TBranch * b_jet3_eta;
361			TBranch * b_jet3_phi;
362			TBranch * b_jet3_energy;
363	kukartse	1.5	TBranch * b_jet4_et;
364	kukartse	1.1	TBranch * b_jet4_pt;
365			TBranch * b_jet4_eta;
366			TBranch * b_jet4_phi;
367			TBranch * b_jet4_energy;
368			TBranch * b_n_tagged_jets_trackCounting_loose;
369			TBranch * b_n_tagged_jets_trackCounting_medium;
370			TBranch * b_n_tagged_jets_trackCounting_tight;
371			TBranch * b_n_tagged_jets_jetProb_loose;
372			TBranch * b_n_tagged_jets_jetProb_medium;
373			TBranch * b_n_tagged_jets_jetProb_tight;
374	kukartse	1.4	// ____
375	kukartse	1.1	TBranch * b_aplanarity ;
376			TBranch * b_centrality ;
377			TBranch * b_sphericity ;
378			TBranch * b_ht ;
379			TBranch * b_htPlusLepton ;
380			TBranch * b_metHtPlusLepton ;
381			TBranch * b_h ;
382			TBranch * b_ktMinPrime ;
383			TBranch * b_dPhiLMet ;
384			TBranch * b_minDijetMass ;
385			TBranch * b_maxJetEta ;
386			TBranch * b_et3 ;
387			TBranch * b_minDiJetDeltaR ;
388			TBranch * b_leptonJetDeltaR ;
389			TBranch * b_ht2p ;
390			TBranch * b_jet1Jet2DeltaR ;
391			TBranch * b_jet1Jet2DeltaPhi ;
392			TBranch * b_jet1Jet2_M ;
393			TBranch * b_jet1Jet2_Pt ;
394			TBranch * b_jet1Jet2W_M ;
395			TBranch * b_jet1Jet2W_Pt ;
396			TBranch * b_hz ;
397	kukartse	1.4	TBranch * b_HT2;
398			TBranch * b_HT2prime;
399			TBranch * b_W_MT;
400	kukartse	1.1	TBranch * b_W_M ;
401			TBranch * b_W_Pt ;
402			TBranch * b_DphiJMET ;
403			// Ht
404			TBranch * b_getHt ;
405			TBranch * b_getHtp ;
406			TBranch * b_getHtpp ;
407			TBranch * b_getHt2 ;
408			TBranch * b_getHt2p ;
409			TBranch * b_getHt2pp ;
410			TBranch * b_getHt3 ;
411			TBranch * b_getHt3p ;
412			TBranch * b_getHt3pp ;
413			TBranch * b_getCen ;
414			TBranch * b_getNJW ;
415			TBranch * b_getJetEtaMax ;
416			TBranch * b_getMdijetMin ;
417			TBranch * b_getMtjets ;
418			TBranch * b_getSqrtsT ;
419			TBranch * b_getMtAurelio ;
420			TBranch * b_getPzOverHT ;
421			TBranch * b_getMevent ;
422			TBranch * b_getM123inv ;
423			TBranch * b_getEta2Sum ;
424			TBranch * b_getMwRec ;
425			TBranch * b_getH ;
426			// event topo
427			TBranch * b_getSph ;
428			TBranch * b_getApl ;
429			TBranch * b_getAplMu ;
430			// Kt
431			TBranch * b_getKtminp ;
432			TBranch * b_getKtminpReduced ;
433			TBranch * b_getDrMinJetJet ;
434			// mT
435			TBranch * b_getDphiMuMet ;
436			TBranch * b_getMt ;
437
438	kukartse	1.4
439	kukartse	1.1	};