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.17 2009/11/18 00:58:13 kukartse Exp $
//
//

#include "DataFormats/BeamSpot/interface/BeamSpot.h"

#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"

#include "DataFormats/Math/interface/LorentzVector.h"
#include "DataFormats/Common/interface/View.h"
#include "DataFormats/PatCandidates/interface/Electron.h"
#include "DataFormats/TrackReco/interface/Track.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);

  // returns pair<dist,dcot>
  pair<float, float> getConversionInfo(const pat::Electron & el,
                                       const reco::Track & track,
                                       const reco::BeamSpot::Point _bs,
                                       float bField);

  bool ConversionAlgo(const pat::Electron& e, edm::Handle<edm::View<reco::Track> > track_h);
  double DelRTrackEle(const pat::Electron& e,const reco::Track& tk);

  bool first_event;

  // _____ quantities read from the config file _________________________
  std::string _jetSource, _electronSource, _muonSource, _METSource,
    _beamSpotSource, _trigger_source, _track_source;
  std::string pdfWeightTag_;
  bool _lepton_flavor_is_muon;
  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     ;
  double _bfield;
  double _dist_min;
  double _dist_max;
  double _abs_dcot_max;
  double _e_track_dr;

  RooGKCounter eventCounter;
  RooGKCounter selectedEvents;

  TTree * _tree;

  //
  // _____ output tree variables ________________________________________
  //
  Int_t _event, ilumi, irun;
  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;
  //
  //_____ photon conversion _____________________________________________
  //
  vector<double> * _e_dist;
  vector<double> * _e_dcot;
  bool     _is_conv;
  //
  //_____ PDF weights ___________________________________________________
  //
  vector<double> * _pdf_weights;
  Double_t _weight_pdf_pos;
  Double_t _weight_pdf_neg;
  //
  //_____ trigger _______________________________________________________
  //
  std::map<std::string,bool> _hlt;
  //vector<TBranch> b_hlt;
  //
  Double_t _event_weight;
  Double_t _lepton_et;
  Double_t _lepton_pt;
  Double_t _lepton_eta;
  Double_t _lepton_phi;
  Double_t _lepton_jet_min_dr;
  Double_t _lepton_ptrel; // relative to the nearest jet
  Double_t _lepton_d0;    // relative to the beam spot
  Double_t _lepton_d0_err;    // relative to the beam spot
  Double_t _lepton_d0_significance;    // relative to the beam spot
  Double_t _lepton_energy;
  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;
  Int_t    _muon_track_nhits;
  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 _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> * _muon_pt;
  vector<double> * _muon_pz;
  vector<double> * _muon_eta;
  vector<double> * _muon_phi;
  vector<double> * _muon_et;
  vector<double> * _muon_energy;
  vector<int> *    _muon_flavor;
  //
  vector<double> * _electron_pt;
  vector<double> * _electron_pz;
  vector<double> * _electron_eta;
  vector<double> * _electron_phi;
  vector<double> * _electron_et;
  vector<double> * _electron_energy;
  vector<int> *    _electron_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_lepton_pt;
  TBranch * b_lepton_pz;
  TBranch * b_lepton_eta;
  TBranch * b_lepton_phi;
  TBranch * b_lepton_et;
  TBranch * b_lepton_energy;
  TBranch * b_lepton_flavor;
  TBranch * b_lepton_jet_min_dr;
  TBranch * b_lepton_ptrel;
  TBranch * b_lepton_d0;
  TBranch * b_lepton_d0_err;
  TBranch * b_lepton_d0_significance;
  //
  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_muon_et;
  TBranch * b_muon_pt;
  TBranch * b_muon_pz;
  TBranch * b_muon_eta;
  TBranch * b_muon_phi;
  TBranch * b_muon_energy;
  TBranch * b_muon_flavor;
  //
  TBranch * b_electron_et;
  TBranch * b_electron_pt;
  TBranch * b_electron_pz;
  TBranch * b_electron_eta;
  TBranch * b_electron_phi;
  TBranch * b_electron_energy;
  TBranch * b_electron_flavor;
  //
  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.18
Committed:	Tue Nov 24 08:02:35 2009 UTC (15 years, 5 months ago) by kukartse
Content type:	text/plain
Branch:	MAIN
CVS Tags:	gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409
Changes since 1.17:	+2 -2 lines
Log Message:	first collisions data
#	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.18	// $Id: TtLJetsAnalyzer.h,v 1.17 2009/11/18 00:58:13 kukartse Exp $
17	kukartse	1.1	//
18			//
19
20	kukartse	1.14	#include "DataFormats/BeamSpot/interface/BeamSpot.h"
21
22	kukartse	1.8	#include "LJMet/MultivariateAnalysis/interface/RooGKCounter.h"
23
24	kukartse	1.1	#include <memory>
25			#include <string>
26	kukartse	1.10	//#include "TFile.h"
27	kukartse	1.8	#include "TTree.h"
28			#include "Rtypes.h"
29	kukartse	1.1	#include "FWCore/Framework/interface/Frameworkfwd.h"
30			#include "FWCore/Framework/interface/EDAnalyzer.h"
31			#include "FWCore/Framework/interface/Event.h"
32			#include "FWCore/Framework/interface/MakerMacros.h"
33			#include "FWCore/ParameterSet/interface/ParameterSet.h"
34
35	kukartse	1.13	#include "DataFormats/Math/interface/LorentzVector.h"
36			#include "DataFormats/Common/interface/View.h"
37			#include "DataFormats/PatCandidates/interface/Electron.h"
38			#include "DataFormats/TrackReco/interface/Track.h"
39	kukartse	1.10
40	kukartse	1.1	class TtLJetsAnalyzer : public edm::EDAnalyzer {
41			public:
42			explicit TtLJetsAnalyzer(const edm::ParameterSet&);
43			~TtLJetsAnalyzer();
44
45
46			private:
47	kukartse	1.4
48			virtual void beginJob(const edm::EventSetup&);
49			virtual void analyze(const edm::Event&, const edm::EventSetup&);
50			virtual void endJob(void);
51
52	kukartse	1.13	// returns pair<dist,dcot>
53			pair<float, float> getConversionInfo(const pat::Electron & el,
54			const reco::Track & track,
55			const reco::BeamSpot::Point _bs,
56			float bField);
57
58			bool ConversionAlgo(const pat::Electron& e, edm::Handle<edm::View<reco::Track> > track_h);
59			double DelRTrackEle(const pat::Electron& e,const reco::Track& tk);
60
61	kukartse	1.12	bool first_event;
62
63	kukartse	1.5	// _____ quantities read from the config file _________________________
64	kukartse	1.12	std::string _jetSource, _electronSource, _muonSource, _METSource,
65	kukartse	1.13	_beamSpotSource, _trigger_source, _track_source;
66	kukartse	1.15	std::string pdfWeightTag_;
67	kukartse	1.10	bool _lepton_flavor_is_muon;
68	kukartse	1.1	std::string _outFile;
69	kukartse	1.8	bool _is_mc;
70	kukartse	1.1	int nCaloJets_min;
71			int nLepton_min;
72	jindal	1.6	double jet_pt_min ;
73	kukartse	1.1	double jet_eta_max ;
74			double muon_pt_min ;
75			double muon_eta_max ;
76			double muon_trackIso_max ;
77			double muon_caloIso_max ;
78			double electron_pt_min ;
79			double electron_eta_max ;
80			double electron_trackIso_max ;
81			double electron_caloIso_max ;
82			double met_et_min ;
83	kukartse	1.7	double min_dr_jet_lepton ;
84	kukartse	1.13	double _bfield;
85			double _dist_min;
86			double _dist_max;
87			double _abs_dcot_max;
88			double _e_track_dr;
89	kukartse	1.1
90			RooGKCounter eventCounter;
91			RooGKCounter selectedEvents;
92
93			TTree * _tree;
94
95	kukartse	1.5	//
96			// _____ output tree variables ________________________________________
97	kukartse	1.12	//
98	kukartse	1.18	Int_t _event, ilumi, irun;
99	kukartse	1.1	Int_t _process_id;
100			Int_t _n_jets;
101	kukartse	1.7	Int_t _n_removed_jets;
102	kukartse	1.1	Int_t _n_met;
103			Int_t _n_muons;
104			Int_t _n_electrons;
105			Int_t _lepton_is_muon;
106	kukartse	1.12	//
107	kukartse	1.13	//_____ photon conversion _____________________________________________
108			//
109			vector<double> * _e_dist;
110			vector<double> * _e_dcot;
111			bool _is_conv;
112			//
113	kukartse	1.15	//_____ PDF weights ___________________________________________________
114			//
115	kukartse	1.16	vector<double> * _pdf_weights;
116	kukartse	1.15	Double_t _weight_pdf_pos;
117			Double_t _weight_pdf_neg;
118			//
119	kukartse	1.12	//_____ trigger _______________________________________________________
120			//
121			std::map<std::string,bool> _hlt;
122			//vector<TBranch> b_hlt;
123			//
124	kukartse	1.1	Double_t _event_weight;
125	kukartse	1.5	Double_t _lepton_et;
126	kukartse	1.1	Double_t _lepton_pt;
127			Double_t _lepton_eta;
128			Double_t _lepton_phi;
129	kukartse	1.11	Double_t _lepton_jet_min_dr;
130	kukartse	1.10	Double_t _lepton_ptrel; // relative to the nearest jet
131			Double_t _lepton_d0; // relative to the beam spot
132			Double_t _lepton_d0_err; // relative to the beam spot
133			Double_t _lepton_d0_significance; // relative to the beam spot
134	kukartse	1.1	Double_t _lepton_energy;
135	kukartse	1.10	Double_t _lepton_track_iso;
136			Double_t _lepton_calo_iso;
137			Double_t _lepton_ecal_iso;
138			Double_t _lepton_hcal_iso;
139			Double_t _lepton_user_iso;
140			Double_t _lepton_IsoDeposit_em;
141			Double_t _lepton_IsoDeposit_had;
142	kukartse	1.12	Int_t _muon_track_nhits;
143	kukartse	1.3	Double_t _muon_track_chi2;
144			Double_t _muon_track_ndof;
145			Double_t _muon_track_nChi2;
146			Double_t _muon_outerTrack_chi2;
147			Double_t _muon_outerTrack_ndof;
148			Double_t _muon_outerTrack_nChi2;
149			Double_t _muon_globalTrack_chi2;
150			Double_t _muon_globalTrack_ndof;
151			Double_t _muon_globalTrack_nChi2;
152	kukartse	1.5	Double_t _met_et;
153	kukartse	1.1	Double_t _met_pt;
154			Double_t _met_eta;
155			Double_t _met_phi;
156			Double_t _met_energy;
157	kukartse	1.5	Double_t _jet1_et;
158	kukartse	1.1	Double_t _jet1_pt;
159			Double_t _jet1_eta;
160			Double_t _jet1_phi;
161			Double_t _jet1_energy;
162	kukartse	1.5	Double_t _jet2_et;
163	kukartse	1.1	Double_t _jet2_pt;
164			Double_t _jet2_eta;
165			Double_t _jet2_phi;
166			Double_t _jet2_energy;
167	kukartse	1.5	Double_t _jet3_et;
168	kukartse	1.1	Double_t _jet3_pt;
169			Double_t _jet3_eta;
170			Double_t _jet3_phi;
171			Double_t _jet3_energy;
172	kukartse	1.5	Double_t _jet4_et;
173	kukartse	1.1	Double_t _jet4_pt;
174			Double_t _jet4_eta;
175			Double_t _jet4_phi;
176			Double_t _jet4_energy;
177			Double_t _aplanarity ;
178			Double_t _centrality ;
179			Double_t _sphericity ;
180			Double_t _ht ;
181			Double_t _htPlusLepton ;
182			Double_t _metHtPlusLepton ;
183			Double_t _h ;
184			Double_t _ktMinPrime ;
185			Double_t _dPhiLMet ;
186			Double_t _minDijetMass ;
187			Double_t _maxJetEta ;
188			Double_t _et3 ;
189			Double_t _minDiJetDeltaR ;
190			Double_t _leptonJetDeltaR ;
191			Double_t _ht2p ;
192			Double_t _jet1Jet2DeltaR ;
193			Double_t _jet1Jet2DeltaPhi ;
194			Double_t _jet1Jet2_M ;
195			Double_t _jet1Jet2_Pt ;
196			Double_t _jet1Jet2W_M ;
197			Double_t _jet1Jet2W_Pt ;
198			Double_t _hz ;
199			Double_t _HT2 ;
200			Double_t _HT2prime ;
201			Double_t _W_MT ;
202			Double_t _W_M ;
203			Double_t _W_Pt ;
204			Double_t _DphiJMET ;
205			// Ht
206			Double_t _getHt ;
207			Double_t _getHtp ;
208			Double_t _getHtpp ;
209			Double_t _getHt2 ;
210			Double_t _getHt2p ;
211			Double_t _getHt2pp ;
212			Double_t _getHt3 ;
213			Double_t _getHt3p ;
214			Double_t _getHt3pp ;
215			Double_t _getCen ;
216			Double_t _getNJW ;
217			Double_t _getJetEtaMax ;
218			Double_t _getMdijetMin ;
219			Double_t _getMtjets ;
220			Double_t _getSqrtsT ;
221			Double_t _getMtAurelio ;
222			Double_t _getPzOverHT ;
223			Double_t _getMevent ;
224			Double_t _getM123inv ;
225			Double_t _getEta2Sum ;
226			Double_t _getMwRec ;
227			Double_t _getH ;
228			// event topo
229			Double_t _getSph ;
230			Double_t _getApl ;
231			Double_t _getAplMu ;
232			// Kt
233			Double_t _getKtminp ;
234			Double_t _getKtminpReduced ;
235			Double_t _getDrMinJetJet ;
236			// mT
237			Double_t _getDphiMuMet ;
238			Double_t _getMt ;
239	kukartse	1.5	//
240			//
241			//_____ b tagging ______________________________________________
242	kukartse	1.8	vector<double> * _jet_pt;
243			vector<double> * _jet_pz;
244			vector<double> * _jet_eta;
245			vector<double> * _jet_phi;
246			vector<double> * _jet_et;
247			vector<double> * _jet_energy;
248			vector<int> * _jet_flavor;
249	kukartse	1.9	//
250			vector<double> * _muon_pt;
251			vector<double> * _muon_pz;
252			vector<double> * _muon_eta;
253			vector<double> * _muon_phi;
254			vector<double> * _muon_et;
255			vector<double> * _muon_energy;
256			vector<int> * _muon_flavor;
257			//
258			vector<double> * _electron_pt;
259			vector<double> * _electron_pz;
260			vector<double> * _electron_eta;
261			vector<double> * _electron_phi;
262			vector<double> * _electron_et;
263			vector<double> * _electron_energy;
264			vector<int> * _electron_flavor;
265	kukartse	1.8	/* FIXME: add vectors of these:
266			Double_t _lepton_track_iso;
267			Double_t _lepton_calo_iso;
268			Double_t _lepton_ecal_iso;
269			Double_t _lepton_hcal_iso;
270			Double_t _lepton_user_iso;
271			Double_t _lepton_IsoDeposit_em;
272			Double_t _lepton_IsoDeposit_had;
273			*/
274			TBranch * b_jet_pt;
275			TBranch * b_jet_pz;
276			TBranch * b_jet_eta;
277			TBranch * b_jet_phi;
278			TBranch * b_jet_et;
279			TBranch * b_jet_energy;
280			TBranch * b_jet_flavor;
281	kukartse	1.9	//
282			TBranch * b_lepton_pt;
283			TBranch * b_lepton_pz;
284			TBranch * b_lepton_eta;
285			TBranch * b_lepton_phi;
286			TBranch * b_lepton_et;
287			TBranch * b_lepton_energy;
288			TBranch * b_lepton_flavor;
289	kukartse	1.11	TBranch * b_lepton_jet_min_dr;
290			TBranch * b_lepton_ptrel;
291			TBranch * b_lepton_d0;
292			TBranch * b_lepton_d0_err;
293			TBranch * b_lepton_d0_significance;
294	kukartse	1.9	//
295	kukartse	1.5	vector<string> * _b_tagger_label;
296			vector<float> * _b_tagger_discr;
297			TBranch * b_b_tagger_label;
298			TBranch * b_b_tagger_discr;
299	kukartse	1.4	Int_t _n_tagged_jets_trackCounting_loose;
300			Int_t _n_tagged_jets_trackCounting_medium;
301			Int_t _n_tagged_jets_trackCounting_tight;
302			Int_t _n_tagged_jets_jetProb_loose;
303			Int_t _n_tagged_jets_jetProb_medium;
304			Int_t _n_tagged_jets_jetProb_tight;
305	kukartse	1.1
306	kukartse	1.4	//
307			//
308			//_____ electron quality _______________________________________
309	kukartse	1.5	Float_t _eidLoose;
310			Float_t _eidRobustHighEnergy;
311			Float_t _eidRobustLoose;
312			Float_t _eidRobustTight;
313			Float_t _eidTight;
314	kukartse	1.4	Int_t _GsfElectron_classification;
315			Double_t _GsfElectron_hadronicOverEm;
316			Double_t _GsfElectron_caloEnergyError;
317			Double_t _GsfElectron_trackMomentumError;
318			Int_t _GsfElectron_numberOfClusters;
319			Double_t _GsfElectron_caloEnergy; // the super cluster energy corrected by EnergyScaleFactor
320			Double_t _GsfElectron_eSuperClusterOverP; // the supercluster energy / track momentum at impact point
321			Double_t _GsfElectron_eSeedClusterOverPout; // the seed cluster energy / track momentum at calo from outermost state
322			Double_t _GsfElectron_deltaEtaSuperClusterTrackAtVtx; // the supercluster eta - track eta from helix extrapolation from impact point
323			Double_t _GsfElectron_deltaEtaSeedClusterTrackAtCalo; // the seed cluster eta - track eta at calo from outermost state
324			Double_t _GsfElectron_deltaPhiSuperClusterTrackAtVtx; // the supercluster phi - track phi from helix extrapolation from impact point
325			Double_t _GsfElectron_deltaPhiSeedClusterTrackAtCalo; // the seed cluster phi - track phi at calo from outermost state
326			// ____ branches
327	kukartse	1.5	TBranch * b_eidLoose;
328			TBranch * b_eidRobustHighEnergy;
329			TBranch * b_eidRobustLoose;
330			TBranch * b_eidRobustTight;
331			TBranch * b_eidTight;
332	kukartse	1.4	TBranch * b_GsfElectron_classification;
333			TBranch * b_GsfElectron_hadronicOverEm;
334			TBranch * b_GsfElectron_caloEnergyError;
335			TBranch * b_GsfElectron_trackMomentumError;
336			TBranch * b_GsfElectron_numberOfClusters;
337			TBranch * b_GsfElectron_caloEnergy; // the super cluster energy corrected by EnergyScaleFactor
338			TBranch * b_GsfElectron_eSuperClusterOverP; // the supercluster energy / track momentum at impact point
339			TBranch * b_GsfElectron_eSeedClusterOverPout; // the seed cluster energy / track momentum at calo from outermost state
340			TBranch * b_GsfElectron_deltaEtaSuperClusterTrackAtVtx; // the supercluster eta - track eta from helix extrapolation from impact point
341			TBranch * b_GsfElectron_deltaEtaSeedClusterTrackAtCalo; // the seed cluster eta - track eta at calo from outermost state
342			TBranch * b_GsfElectron_deltaPhiSuperClusterTrackAtVtx; // the supercluster phi - track phi from helix extrapolation from impact point
343			TBranch * b_GsfElectron_deltaPhiSeedClusterTrackAtCalo; // the seed cluster phi - track phi at calo from outermost state
344			//
345			//
346			//_____ muon quality _______________________________________
347			bool _muon_isGood_All; // dummy options - always true
348			bool _muon_isGood_AllGlobalMuons; // checks isGlobalMuon flag
349			bool _muon_isGood_AllStandAloneMuons; // checks isStandAloneMuon flag
350			bool _muon_isGood_AllTrackerMuons; // checks isTrackerMuon flag
351			bool _muon_isGood_TrackerMuonArbitrated; // resolve ambiguity of sharing segments
352			bool _muon_isGood_AllArbitrated; // all muons with the tracker muon arbitrated
353			bool _muon_isGood_GlobalMuonPromptTight; // global muons with tighter fit requirements
354			bool _muon_isGood_TMLastStationLoose; // penetration depth loose selector
355			bool _muon_isGood_TMLastStationTight; // penetration depth tight selector
356			bool _muon_isGood_TM2DCompatibilityLoose; // likelihood based loose selector
357			bool _muon_isGood_TM2DCompatibilityTight; // likelihood based tight selector
358			bool _muon_isGood_TMOneStationLoose; // require one well matched segment
359			bool _muon_isGood_TMOneStationTight; // require one well matched segment
360			bool _muon_isGood_TMLastStationOptimizedLowPtLoose; // combination of TMLastStation and TMOneStation
361			bool _muon_isGood_TMLastStationOptimizedLowPtTight; // combination of TMLastStation and TMOneStation
362			//____ branches
363			TBranch * b_muon_isGood_All; // dummy options - always true
364			TBranch * b_muon_isGood_AllGlobalMuons; // checks isGlobalMuon flag
365			TBranch * b_muon_isGood_AllStandAloneMuons; // checks isStandAloneMuon flag
366			TBranch * b_muon_isGood_AllTrackerMuons; // checks isTrackerMuon flag
367			TBranch * b_muon_isGood_TrackerMuonArbitrated; // resolve ambiguity of sharing segments
368			TBranch * b_muon_isGood_AllArbitrated; // all muons with the tracker muon arbitrated
369			TBranch * b_muon_isGood_GlobalMuonPromptTight; // global muons with tighter fit requirements
370			TBranch * b_muon_isGood_TMLastStationLoose; // penetration depth loose selector
371			TBranch * b_muon_isGood_TMLastStationTight; // penetration depth tight selector
372			TBranch * b_muon_isGood_TM2DCompatibilityLoose; // likelihood based loose selector
373			TBranch * b_muon_isGood_TM2DCompatibilityTight; // likelihood based tight selector
374			TBranch * b_muon_isGood_TMOneStationLoose; // require one well matched segment
375			TBranch * b_muon_isGood_TMOneStationTight; // require one well matched segment
376			TBranch * b_muon_isGood_TMLastStationOptimizedLowPtLoose; // combination of TMLastStation and TMOneStation
377			TBranch * b_muon_isGood_TMLastStationOptimizedLowPtTight; // combination of TMLastStation and TMOneStation
378			//____ output tree variables
379	kukartse	1.1	TBranch * b_event;
380			TBranch * b_process_id;
381			TBranch * b_n_jets;
382	kukartse	1.8	TBranch * b_n_removed_jets;
383	kukartse	1.1	TBranch * b_n_met;
384			TBranch * b_n_muons;
385			TBranch * b_n_electrons;
386			TBranch * b_lepton_is_muon;
387			TBranch * b_event_weight;
388	kukartse	1.9	//
389			TBranch * b_muon_et;
390			TBranch * b_muon_pt;
391			TBranch * b_muon_pz;
392			TBranch * b_muon_eta;
393			TBranch * b_muon_phi;
394			TBranch * b_muon_energy;
395			TBranch * b_muon_flavor;
396			//
397			TBranch * b_electron_et;
398			TBranch * b_electron_pt;
399			TBranch * b_electron_pz;
400			TBranch * b_electron_eta;
401			TBranch * b_electron_phi;
402			TBranch * b_electron_energy;
403			TBranch * b_electron_flavor;
404			//
405	kukartse	1.3	TBranch * b_muon_track_chi2;
406			TBranch * b_muon_track_ndof;
407			TBranch * b_muon_track_nChi2;
408			TBranch * b_muon_outerTrack_chi2;
409			TBranch * b_muon_outerTrack_ndof;
410			TBranch * b_muon_outerTrack_nChi2;
411			TBranch * b_muon_globalTrack_chi2;
412			TBranch * b_muon_globalTrack_ndof;
413			TBranch * b_muon_globalTrack_nChi2;
414	kukartse	1.1	TBranch * b_lepton_track_iso;
415			TBranch * b_lepton_calo_iso;
416	kukartse	1.2	TBranch * b_lepton_ecal_iso;
417			TBranch * b_lepton_hcal_iso;
418			TBranch * b_lepton_user_iso;
419	kukartse	1.8	TBranch * b_lepton_IsoDeposit_em;
420			TBranch * b_lepton_IsoDeposit_had;
421	kukartse	1.5	TBranch * b_met_et;
422	kukartse	1.1	TBranch * b_met_pt;
423			TBranch * b_met_eta;
424			TBranch * b_met_phi;
425			TBranch * b_met_energy;
426	kukartse	1.5	TBranch * b_jet1_et;
427	kukartse	1.1	TBranch * b_jet1_pt;
428			TBranch * b_jet1_eta;
429			TBranch * b_jet1_phi;
430			TBranch * b_jet1_energy;
431	kukartse	1.5	TBranch * b_jet2_et;
432	kukartse	1.1	TBranch * b_jet2_pt;
433			TBranch * b_jet2_eta;
434			TBranch * b_jet2_phi;
435			TBranch * b_jet2_energy;
436	kukartse	1.5	TBranch * b_jet3_et;
437	kukartse	1.1	TBranch * b_jet3_pt;
438			TBranch * b_jet3_eta;
439			TBranch * b_jet3_phi;
440			TBranch * b_jet3_energy;
441	kukartse	1.5	TBranch * b_jet4_et;
442	kukartse	1.1	TBranch * b_jet4_pt;
443			TBranch * b_jet4_eta;
444			TBranch * b_jet4_phi;
445			TBranch * b_jet4_energy;
446			TBranch * b_n_tagged_jets_trackCounting_loose;
447			TBranch * b_n_tagged_jets_trackCounting_medium;
448			TBranch * b_n_tagged_jets_trackCounting_tight;
449			TBranch * b_n_tagged_jets_jetProb_loose;
450			TBranch * b_n_tagged_jets_jetProb_medium;
451			TBranch * b_n_tagged_jets_jetProb_tight;
452	kukartse	1.4	// ____
453	kukartse	1.1	TBranch * b_aplanarity ;
454			TBranch * b_centrality ;
455			TBranch * b_sphericity ;
456			TBranch * b_ht ;
457			TBranch * b_htPlusLepton ;
458			TBranch * b_metHtPlusLepton ;
459			TBranch * b_h ;
460			TBranch * b_ktMinPrime ;
461			TBranch * b_dPhiLMet ;
462			TBranch * b_minDijetMass ;
463			TBranch * b_maxJetEta ;
464			TBranch * b_et3 ;
465			TBranch * b_minDiJetDeltaR ;
466			TBranch * b_leptonJetDeltaR ;
467			TBranch * b_ht2p ;
468			TBranch * b_jet1Jet2DeltaR ;
469			TBranch * b_jet1Jet2DeltaPhi ;
470			TBranch * b_jet1Jet2_M ;
471			TBranch * b_jet1Jet2_Pt ;
472			TBranch * b_jet1Jet2W_M ;
473			TBranch * b_jet1Jet2W_Pt ;
474			TBranch * b_hz ;
475	kukartse	1.4	TBranch * b_HT2;
476			TBranch * b_HT2prime;
477			TBranch * b_W_MT;
478	kukartse	1.1	TBranch * b_W_M ;
479			TBranch * b_W_Pt ;
480			TBranch * b_DphiJMET ;
481			// Ht
482			TBranch * b_getHt ;
483			TBranch * b_getHtp ;
484			TBranch * b_getHtpp ;
485			TBranch * b_getHt2 ;
486			TBranch * b_getHt2p ;
487			TBranch * b_getHt2pp ;
488			TBranch * b_getHt3 ;
489			TBranch * b_getHt3p ;
490			TBranch * b_getHt3pp ;
491			TBranch * b_getCen ;
492			TBranch * b_getNJW ;
493			TBranch * b_getJetEtaMax ;
494			TBranch * b_getMdijetMin ;
495			TBranch * b_getMtjets ;
496			TBranch * b_getSqrtsT ;
497			TBranch * b_getMtAurelio ;
498			TBranch * b_getPzOverHT ;
499			TBranch * b_getMevent ;
500			TBranch * b_getM123inv ;
501			TBranch * b_getEta2Sum ;
502			TBranch * b_getMwRec ;
503			TBranch * b_getH ;
504			// event topo
505			TBranch * b_getSph ;
506			TBranch * b_getApl ;
507			TBranch * b_getAplMu ;
508			// Kt
509			TBranch * b_getKtminp ;
510			TBranch * b_getKtminpReduced ;
511			TBranch * b_getDrMinJetJet ;
512			// mT
513			TBranch * b_getDphiMuMet ;
514			TBranch * b_getMt ;
515
516	kukartse	1.4
517	kukartse	1.1	};