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.15 2009/06/11 19:18:34 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 "FWCore/ServiceRegistry/interface/Service.h"
#include "PhysicsTools/UtilAlgos/interface/TFileService.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);

  edm::Service<TFileService> fs;

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