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.18 2009/11/24 08:02:35 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);

  void debug(int); // debugging cout marker

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