#include #include #include "CMS2.h" #include "goodrun.h" #include "electronSelections.h" #include "eventSelections.h" #include #include using namespace std; using namespace tas; using namespace wp2012; int main(int argc, char** argv) { //void ElectronLooper(TString inputList) { // PassAllWP2012Cuts enum Classification { UNKNOWN=-1, GOLDEN=0, BIGBREM=1, BADTRACK=2, SHOWERING=3, GAP=4 } ; TString reader; ifstream fileReader(argv[1]); //TString dirout = argv[2]; // gSystem->Exec("mkdir -p " + dirout ); // TFile* outputFile = TFile::Open(TString::Format("%s/%s.root",argv[2],argv[3]), "recreate"); TFile* outputFile = TFile::Open("testNtuple2.root", "recreate"); TTree* redTree = new TTree ("redTree","a small Tree"); /*------------------------------ this is for the mass---------------------------------------------------------------------------------------------*/ //Float_t dielectron_mass; // Float_t dielectron_mass,dielectron_mass_lowfbrem,dielectron_mass_onecluster,dielectron_mass_golden; // RCLSA: this is a tree, not a histogram LorentzVector dielectron_p4; float dielectron_mass; float dielectron_pt; float scale1fb; int run; int event; int lumi; //-----------------------------------------------------electron0----------------------------------------------------------------------// LorentzVector electron0_p4In; LorentzVector electron0_p4Out; LorentzVector electron0_trk_p4; LorentzVector electron0_vrtx_p4; LorentzVector electron0_p4; int electron0_class; float electron0_fbrem; float electron0_eSC; float electron0_z0; float electron0_lostHits; Bool_t electron0_pass[4]; float electron0_ecalEnergy; float electron0_eSCPresh; float electron0_HoE; float electron0_dEtaIn; float electron0_dEtaOut; float electron0_dPhiIn; float electron0_dPhiOut; float electron0_r9; int electron0_charge; float electron0_sigmaEtaEta; float electron0_sigmaPhiPhi; float electron0_etaSC; float electron0_nSeed; float electron0_eOverPIn; float electron0_eOverPOut; float electron0_ip3d; float electron0_eSeed; float eletrcon0_phiSC; float electron0_trkIso; float electron0_trkIso04; float electron0_d0; float electron0_eSeedOverPOut; float electron0_eSeedOverPIn; float electron0_ecalEnergyError; float electron0_ecalIso; float electron0_ecalIso04; float electron0_etaErr; float electron0_etaSCwidth; float electron0_hcalDepth1OverEcal; float electron0_hcalDepth1TowerSumEt; float electron0_hcalDepth1TowerSumEt04; float electron0_hcalDepth2OverEcal; float electron0_hcalDepth2TowerSumEt; float electron0_hcalDepth2TowerSumEt04; float electron0_hcalIso; float electron0_hcalIso04; float electron0_ip3derr; float electron0_iso03_pf; float electron0_iso03_pf2012_ch; float electron0_iso03_pf2012_em; float electron0_iso03_pf2012_nh; float electron0_iso03_pf_ch; float electron0_iso03_pf_gamma05; float electron0_iso03_pf_nhad05; float electron0_iso04_pf; float electron0_iso04_pf2012_ch; float electron0_iso04_pf2012_em; float electron0_iso04_pf2012_nh; float electron0_iso04_pf_ch; float electron0_iso04_pf_gamma05; float electron0_iso04_pf_nhad05; float electron0_layer1_charge; float electron0_lh; float electron0_mva; float electron0_ndof; float electron0_pfChargedHadronIso; float electron0_pfNeutralHadronIso; float electron0_pfPhotonIso; float electron0_phiErr; float electron0_phiSCwidth; float electron0_ptErr; float electron0_sigmaIEtaIEta; float electron0_sigmaIEtaIEtaSC; float electron0_sigmaIEtaIPhi; float electron0_sigmaIPhiIPhi; float electron0_sigmaIPhiIPhiSC; float electron0_trackMomentumError; float electron0_trkdr; float electron0_trkshFrac; float electron0_z0Err; float electron0_z0corr; int electron0_mc3_id; int electron0_mc3idx; int electron0_mc3_motherid; int electron0_mc3_motheridx; int electron0_mc_id; int electron0_mcidx; int electron0_mc_motherid; int electron0_closestJet; int electron0_closestMuon; int electron0_pfelectron0idx; int electron0_category; int electron0_ckf_laywithmeas; int electron0_conv_delMissHits; int electron0_conv_flag; int electron0_conv_gsftkidx; int electron0_conv_old_delMissHits; int electron0_conv_old_flag; int electron0_conv_old_gsftkidx; int electron0_conv_old_tkidx; int electron0_conv_tkidx; int electron0_exp_innerlayers; int electron0_exp_outerlayers; int electron0_fiduciality; int electron0_gsftrkidx; int electron0_layer1_det; int electron0_layer1_layer; int electron0_layer1_sizerphi; int electron0_layer1_sizerz; int electron0_lost_pixelhits; int electron0_sccharge; int electron0_scindex; int electron0_trk_charge; int electron0_trkidx; int electron0_type; int electron0_validHits; int electron0_valid_pixelhits; float electron0_ecalJuraIso; float electron0_ecalJuraTowerIso; float electron0_hcalConeIso; float electron0_tkJuraIso; float electron0_jetdr; float electron0_musdr; float electron0_isoR03_chpf_radial; float electron0_isoR03_chpf_radialTight; float electron0_isoR03_chpf_radialTight_bv; float electron0_isoR03_chpf_radial_bv; float electron0_isoR03_empf_radial; float electron0_isoR03_empf_radialTight; float electron0_isoR03_empf_radialTight_bv; float electron0_isoR03_empf_radial_bv; float electron0_isoR03_nhpf_radial; float electron0_isoR03_nhpf_radialTight; float electron0_isoR03_nhpf_radialTight_bv; float electron0_isoR03_nhpf_radial_bv; float electron0_isoR03_pf2012n0p5_ch; float electron0_isoR03_pf2012n0p5_em; float electron0_isoR03_pf2012n0p5_nh; float electron0_isoR03_pf_radial; float electron0_isoR03_pf_radialTight; float electron0_isoR03_pf_radialTight_bv; float electron0_isoR03_pf_radial_bv; float electron0_chi2; float electron0_ckf_chi2; float electron0_ckf_ndof; float electron0_conv_dcot; float electron0_conv_dist; float electron0_conv_old_dcot; float electron0_conv_old_dist; float electron0_conv_old_radius; float electron0_conv_radius; float electron0_d0Err; float electron0_d0corr; float electron0_dPhiInPhiOut; float electron0_deltaEtaEleClusterTrackAtCalo; float electron0_deltaPhiEleClusterTrackAtCalo; LorentzVector electron1_p4In; LorentzVector electron1_p4Out; LorentzVector electron1_trk_p4; LorentzVector electron1_vrtx_p4; LorentzVector electron1_p4; int electron1_class; float electron1_fbrem; float electron1_eSC; float electron1_z0; float electron1_lostHits; Bool_t electron1_pass[4]; float electron1_ecalEnergy; float electron1_eSCPresh; float electron1_HoE; float electron1_dEtaIn; float electron1_dEtaOut; float electron1_dPhiIn; float electron1_dPhiOut; float electron1_r9; int electron1_charge; float electron1_sigmaEtaEta; float electron1_sigmaPhiPhi; float electron1_etaSC; float electron1_nSeed; float electron1_eOverPIn; float electron1_eOverPOut; float electron1_ip3d; float electron1_eSeed; float eletrcon1_phiSC; float electron1_trkIso; float electron1_trkIso04; float electron1_d0; float electron1_eSeedOverPOut; float electron1_eSeedOverPIn; float electron1_ecalEnergyError; float electron1_ecalIso; float electron1_ecalIso04; float electron1_etaErr; float electron1_etaSCwidth; float electron1_hcalDepth1OverEcal; float electron1_hcalDepth1TowerSumEt; float electron1_hcalDepth1TowerSumEt04; float electron1_hcalDepth2OverEcal; float electron1_hcalDepth2TowerSumEt; float electron1_hcalDepth2TowerSumEt04; float electron1_hcalIso; float electron1_hcalIso04; float electron1_ip3derr; float electron1_iso03_pf; float electron1_iso03_pf2012_ch; float electron1_iso03_pf2012_em; float electron1_iso03_pf2012_nh; float electron1_iso03_pf_ch; float electron1_iso03_pf_gamma05; float electron1_iso03_pf_nhad05; float electron1_iso04_pf; float electron1_iso04_pf2012_ch; float electron1_iso04_pf2012_em; float electron1_iso04_pf2012_nh; float electron1_iso04_pf_ch; float electron1_iso04_pf_gamma05; float electron1_iso04_pf_nhad05; float electron1_layer1_charge; float electron1_lh; float electron1_mva; float electron1_ndof; float electron1_pfChargedHadronIso; float electron1_pfNeutralHadronIso; float electron1_pfPhotonIso; float electron1_phiErr; float electron1_phiSC; float electron1_phiSCwidth; float electron1_ptErr; float electron1_sigmaIEtaIEta; float electron1_sigmaIEtaIEtaSC; float electron1_sigmaIEtaIPhi; float electron1_sigmaIPhiIPhi; float electron1_sigmaIPhiIPhiSC; float electron1_trackMomentumError; float electron1_trkdr; float electron1_trkshFrac; float electron1_z0Err; float electron1_z0corr; int electron1_mc3_id; int electron1_mc3idx; int electron1_mc3_motherid; int electron1_mc3_motheridx; int electron1_mc_id; int electron1_mcidx; int electron1_mc_motherid; int electron1_closestJet; int electron1_closestMuon; // int electron1_pfelectron0idx; int electron1_category; int electron1_ckf_laywithmeas; int electron1_conv_delMissHits; int electron1_conv_flag; int electron1_conv_gsftkidx; int electron1_conv_old_delMissHits; int electron1_conv_old_flag; // int electron1_conv_old_gsftkidx; // int electron1_conv_old_tkidx; int electron1_conv_tkidx; int electron1_exp_innerlayers; int electron1_exp_outerlayers; int electron1_fiduciality; int electron1_gsftrkidx; int electron1_layer1_det; int electron1_layer1_layer; int electron1_layer1_sizerphi; int electron1_layer1_sizerz; int electron1_lost_pixelhits; int electron1_sccharge; int electron1_scindex; int electron1_trk_charge; int electron1_trkidx; int electron1_type; int electron1_validHits; int electron1_valid_pixelhits; float electron1_ecalJuraIso; float electron1_ecalJuraTowerIso; float electron1_hcalConeIso; float electron1_tkJuraIso; float electron1_jetdr; float electron1_musdr; float electron1_isoR03_chpf_radial; float electron1_isoR03_chpf_radialTight; float electron1_isoR03_chpf_radialTight_bv; float electron1_isoR03_chpf_radial_bv; float electron1_isoR03_empf_radial; float electron1_isoR03_empf_radialTight; float electron1_isoR03_empf_radialTight_bv; float electron1_isoR03_empf_radial_bv; float electron1_isoR03_nhpf_radial; float electron1_isoR03_nhpf_radialTight; float electron1_isoR03_nhpf_radialTight_bv; float electron1_isoR03_nhpf_radial_bv; float electron1_isoR03_pf2012n0p5_ch; float electron1_isoR03_pf2012n0p5_em; float electron1_isoR03_pf2012n0p5_nh; float electron1_isoR03_pf_radial; float electron1_isoR03_pf_radialTight; float electron1_isoR03_pf_radialTight_bv; float electron1_isoR03_pf_radial_bv; float electron1_chi2; float electron1_ckf_chi2; float electron1_ckf_ndof; float electron1_conv_dcot; float electron1_conv_dist; float electron1_conv_old_dcot; float electron1_conv_old_dist; float electron1_conv_old_radius; float electron1_conv_radius; float electron1_d0Err; float electron1_d0corr; float electron1_dPhiInPhiOut; float electron1_deltaEtaEleClusterTrackAtCalo; float electron1_deltaPhiEleClusterTrackAtCalo; redTree->Branch("dielectron_p4", &dielectron_p4); redTree->Branch("dielectron_mass", &dielectron_mass); redTree->Branch("dielectron_pt", &dielectron_pt); redTree->Branch("scale1fb", &scale1fb); redTree->Branch("run", &run); redTree->Branch("event", &event); redTree->Branch("lumi", &lumi); //------------------------------------electron0-------------------------------------------------------// redTree->Branch("electron0_p4In", &electron0_p4In); redTree->Branch("electron0_p4Out", &electron0_p4Out); redTree->Branch("electron0_trk_p4", &electron0_trk_p4); redTree->Branch("electron0_vrtx_p4", &electron0_vrtx_p4); redTree->Branch("electron0_p4", &electron0_p4); redTree->Branch("electron0_class", &electron0_class); redTree->Branch("electron0_fbrem", &electron0_fbrem); redTree->Branch("electron0_eSC", &electron0_eSC); redTree->Branch("electron0_z0", &electron0_z0); redTree->Branch("electron0_lostHits", &electron0_lostHits); redTree->Branch("electron0_pass", electron0_pass, "electron0_pass[4]/O"); redTree->Branch("electron0_ecalEnergy", &electron0_ecalEnergy); redTree->Branch("electron0_eSCPresh ", &electron0_eSCPresh); redTree->Branch("electron0_HoE", &electron0_HoE); redTree->Branch("electron0_dEtaIn", &electron0_dEtaIn); redTree->Branch("electron0_dEtaOut", &electron0_dEtaOut); redTree->Branch("electron0_dPhiIn", &electron0_dPhiIn); redTree->Branch("electron0_dPhiOut", &electron0_dPhiOut); redTree->Branch("electron0_r9", &electron0_r9); redTree->Branch("electron0_charge", &electron0_charge); redTree->Branch("electron0_sigmaEtaEta", &electron0_sigmaEtaEta); redTree->Branch("electron0_sigmaPhiPhi", &electron0_sigmaPhiPhi); redTree->Branch("electron0_etaSC", &electron0_etaSC); redTree->Branch("electron0_nSeed ", &electron0_nSeed); redTree->Branch("electron0_eOverPIn", &electron0_eOverPIn); redTree->Branch("electron0_eOverPOut", &electron0_eOverPOut); redTree->Branch("electron0_ip3d ", &electron0_ip3d); redTree->Branch("electron0_eSeed ", &electron0_eSeed); redTree->Branch("eletrcon0_phiSC ", &eletrcon0_phiSC); redTree->Branch("electron0_trkIso ", &electron0_trkIso); redTree->Branch("electron0_trkIso04", &electron0_trkIso04); redTree->Branch("electron0_d0", &electron0_d0); redTree->Branch("electron0_eSeedOverPIn", &electron0_eSeedOverPIn); redTree->Branch("electron0_eSeedOverPOut", &electron0_eSeedOverPOut); redTree->Branch("electron0_ecalEnergyError", &electron0_ecalEnergyError); redTree->Branch("electron0_ecalIso", &electron0_ecalIso); redTree->Branch("electron0_ecalIso04", &electron0_ecalIso04); redTree->Branch("electron0_etaErr", &electron0_etaErr); redTree->Branch("electron0_etaSCwidth", &electron0_etaSCwidth); redTree->Branch("electron0_hcalDepth1OverEcal", &electron0_hcalDepth1OverEcal); redTree->Branch("electron0_hcalDepth1TowerSumEt", &electron0_hcalDepth1TowerSumEt); redTree->Branch("electron0_hcalDepth1TowerSumEt04", &electron0_hcalDepth1TowerSumEt04); redTree->Branch("electron0_hcalDepth2OverEcal", &electron0_hcalDepth2OverEcal); redTree->Branch("electron0_hcalDepth2TowerSumEt", &electron0_hcalDepth2TowerSumEt); redTree->Branch("electron0_hcalDepth2TowerSumEt04", &electron0_hcalDepth2TowerSumEt04); redTree->Branch("electron0_hcalIso ", &electron0_hcalIso); redTree->Branch("electron0_hcalIso04", &electron0_hcalIso04); redTree->Branch("electron0_ip3derr", &electron0_ip3derr); redTree->Branch("electron0_iso03_pf", &electron0_iso03_pf); redTree->Branch("electron0_iso03_pf2012_ch", &electron0_iso03_pf2012_ch); redTree->Branch("electron0_iso03_pf2012_em ", &electron0_iso03_pf2012_em); redTree->Branch("electron0_iso03_pf2012_nh", &electron0_iso03_pf2012_nh); redTree->Branch("electron0_iso03_pf_ch", &electron0_iso03_pf_ch); redTree->Branch("electron0_iso03_pf_gamma05", &electron0_iso03_pf_gamma05); redTree->Branch("electron0_iso03_pf_nhad05", &electron0_iso03_pf_nhad05); redTree->Branch("electron0_iso04_pf", &electron0_iso04_pf); redTree->Branch("electron0_iso04_pf2012_ch", &electron0_iso04_pf2012_ch); redTree->Branch("electron0_iso04_pf2012_em", &electron0_iso04_pf2012_em); redTree->Branch("electron0_iso04_pf2012_nh", &electron0_iso04_pf2012_nh); redTree->Branch("electron0_iso04_pf_ch", &electron0_iso04_pf_ch); redTree->Branch("electron0_iso04_pf_gamma05", &electron0_iso04_pf_gamma05); redTree->Branch("electron0_iso04_pf_nhad05", &electron0_iso04_pf_nhad05); redTree->Branch("electron0_layer1_charge", &electron0_layer1_charge); redTree->Branch("electron0_lh ", &electron0_lh); redTree->Branch("electron0_mva ", &electron0_mva); redTree->Branch("electron0_ndof", &electron0_ndof); redTree->Branch("electron0_pfChargedHadronIso", &electron0_pfChargedHadronIso); redTree->Branch("electron0_pfNeutralHadronIso", &electron0_pfNeutralHadronIso); redTree->Branch("electron0_pfPhotonIso", &electron0_pfPhotonIso); redTree->Branch("electron0_phiErr", &electron0_phiErr); redTree->Branch("electron0_phiSCwidth", &electron0_phiSCwidth); redTree->Branch("electron0_ptErr", &electron0_ptErr); redTree->Branch("electron0_sigmaIEtaIEta ", &electron0_sigmaIEtaIEta); redTree->Branch("electron0_sigmaIEtaIEtaSC", &electron0_sigmaIEtaIEtaSC); redTree->Branch("electron0_sigmaIEtaIPhi", &electron0_sigmaIEtaIPhi); redTree->Branch("electron0_sigmaIPhiIPhi", &electron0_sigmaIPhiIPhi); redTree->Branch("electron0_sigmaIPhiIPhiSC", &electron0_sigmaIPhiIPhiSC); redTree->Branch("electron0_trackMomentumError", &electron0_trackMomentumError); redTree->Branch("electron0_trkdr", &electron0_trkdr); redTree->Branch("electron0_trkshFrac", &electron0_trkshFrac); redTree->Branch("electron0_z0Err", &electron0_z0Err); redTree->Branch("electron0_z0corr", &electron0_z0corr); redTree->Branch("electron0_mc3_id", &electron0_mc3_id); redTree->Branch("electron0_mc3idx", &electron0_mc3idx); redTree->Branch("electron0_mc3_motherid", &electron0_mc3_motherid); redTree->Branch("electron0_mc3_motheridx", &electron0_mc3_motheridx); redTree->Branch("electron0_mc_id", &electron0_mc_id); redTree->Branch("electron0_mcidx", &electron0_mcidx); redTree->Branch("electron0_mc_motherid", &electron0_mc_motherid); redTree->Branch("electron0_closestJet", &electron0_closestJet); redTree->Branch("electron0_closestMuon", &electron0_closestMuon); // redTree->Branch("electron0_pfelectron0idx", &electron0_pfelectron0idx); redTree->Branch("electron0_category", &electron0_category); redTree->Branch("electron0_ckf_laywithmeas", &electron0_ckf_laywithmeas); redTree->Branch("electron0_conv_delMissHits", &electron0_conv_delMissHits); redTree->Branch("electron0_conv_flag", &electron0_conv_flag); redTree->Branch("electron0_conv_gsftkidx", &electron0_conv_gsftkidx); redTree->Branch("electron0_conv_old_delMissHits", &electron0_conv_old_delMissHits); redTree->Branch("electron0_conv_old_flag ", &electron0_conv_old_flag); // redTree->Branch("electron0_conv_old_gsftkidx ", &electron0_conv_old_gsftkidx); // redTree->Branch("electron0_conv_old_tkidx", &electron0_conv_old_tkidx ); redTree->Branch("electron0_conv_tkidx", &electron0_conv_tkidx); redTree->Branch("electron0_exp_innerlayers", &electron0_exp_innerlayers); redTree->Branch("electron0_exp_outerlayers", &electron0_exp_outerlayers); redTree->Branch("electron0_fiduciality", &electron0_fiduciality); redTree->Branch("electron0_gsftrkidx", &electron0_gsftrkidx); redTree->Branch("electron0_layer1_det", &electron0_layer1_det); redTree->Branch("electron0_layer1_layer", &electron0_layer1_layer); redTree->Branch("electron0_layer1_sizerphi", &electron0_layer1_sizerphi); redTree->Branch("electron0_layer1_sizerz", &electron0_layer1_sizerz ); redTree->Branch("electron0_lost_pixelhits", &electron0_lost_pixelhits ); redTree->Branch("electron0_sccharge", &electron0_sccharge); redTree->Branch("electron0_scindex", &electron0_scindex); redTree->Branch("electron0_trk_charge", &electron0_trk_charge); redTree->Branch("electron0_trkidx", &electron0_trkidx); redTree->Branch("electron0_type", &electron0_type); redTree->Branch("electron0_validHits", &electron0_validHits); redTree->Branch("electron0_valid_pixelhits", &electron0_valid_pixelhits); redTree->Branch("electron0_calJuraTowerIso", &electron0_ecalJuraTowerIso ); redTree->Branch("electron0_calConeIso", &electron0_hcalConeIso ); redTree->Branch("electron0_kJuraIso ", &electron0_tkJuraIso); redTree->Branch("electron0_etdr", &electron0_jetdr); redTree->Branch("electron0_usdr", &electron0_musdr); redTree->Branch("electron0_soR03_chpf_radial", &electron0_isoR03_chpf_radial); redTree->Branch("electron0_soR03_chpf_radialTight", &electron0_isoR03_chpf_radialTight); redTree->Branch("electron0_soR03_chpf_radialTight_bv", &electron0_isoR03_chpf_radialTight_bv); redTree->Branch("electron0_soR03_chpf_radial_bv", &electron0_isoR03_chpf_radial_bv); redTree->Branch("electron0_soR03_empf_radial", &electron0_isoR03_empf_radial); redTree->Branch("electron0_soR03_empf_radialTight", &electron0_isoR03_empf_radialTight); redTree->Branch("electron0_soR03_empf_radialTight_bv", &electron0_isoR03_empf_radialTight_bv); redTree->Branch("electron0_isoR03_empf_radial_bv", &electron0_isoR03_empf_radial_bv); redTree->Branch("electron0_isoR03_nhpf_radial", &electron0_isoR03_nhpf_radial); redTree->Branch("electron0_isoR03_nhpf_radialTight", &electron0_isoR03_nhpf_radialTight); redTree->Branch("electron0_isoR03_nhpf_radialTight_bv", &electron0_isoR03_nhpf_radialTight_bv); redTree->Branch("electron0_isoR03_nhpf_radial_bv", &electron0_isoR03_nhpf_radial_bv); redTree->Branch("electron0_isoR03_pf2012n0p5_ch", &electron0_isoR03_pf2012n0p5_ch); redTree->Branch("electron0_isoR03_pf2012n0p5_em", &electron0_isoR03_pf2012n0p5_em); redTree->Branch("electron0_isoR03_pf2012n0p5_nh", &electron0_isoR03_pf2012n0p5_nh); redTree->Branch("electron0_isoR03_pf_radial", &electron0_isoR03_pf_radial); redTree->Branch("electron0_isoR03_pf_radialTight", &electron0_isoR03_pf_radialTight); redTree->Branch("electron0_isoR03_pf_radialTight_bv", &electron0_isoR03_pf_radialTight_bv); redTree->Branch("electron0_isoR03_pf_radial_bv", &electron0_isoR03_pf_radial_bv); redTree->Branch("electron0_chi2", &electron0_chi2); redTree->Branch("electron0_ckf_chi2", &electron0_ckf_chi2); redTree->Branch("electron0_ckf_ndof", &electron0_ckf_ndof); redTree->Branch("electron0_conv_dcot", &electron0_conv_dcot); redTree->Branch("electron0_conv_dist", &electron0_conv_dist); redTree->Branch("electron0_conv_old_dcot", &electron0_conv_old_dcot); redTree->Branch("electron0_conv_old_dist", &electron0_conv_old_dist); redTree->Branch("electron0_conv_old_radius", &electron0_conv_old_radius); redTree->Branch("electron0_conv_radius", &electron0_conv_radius); redTree->Branch("electron0_d0Err", &electron0_d0Err); redTree->Branch("electron0_d0corr", &electron0_d0corr); redTree->Branch("electron0_dPhiInPhiOut", &electron0_dPhiInPhiOut); redTree->Branch("electron0_deltaEtaEleClusterTrackAtCal", &electron0_deltaEtaEleClusterTrackAtCalo); redTree->Branch("electron0_deltaPhiEleClusterTrackAtCal", &electron0_deltaPhiEleClusterTrackAtCalo); //--------------------------------electron1------------------------------------------------------------------// redTree->Branch("electron1_p4In", &electron1_p4In); redTree->Branch("electron1_p4Out", &electron1_p4Out); redTree->Branch("electron1_trk_p4", &electron1_trk_p4); redTree->Branch("electron1_vrtx_p4", &electron1_vrtx_p4); redTree->Branch("electron1_p4", &electron1_p4); redTree->Branch("electron1_class", &electron1_class); redTree->Branch("electron1_fbrem", &electron1_fbrem); redTree->Branch("electron1_eSC", &electron1_eSC); redTree->Branch("electron1_z0", &electron1_z0); redTree->Branch("electron1_lostHits", &electron1_lostHits); redTree->Branch("electron1_pass", electron1_pass, "electron1_pass[4]/O"); redTree->Branch("electron1_ecalEnergy", &electron1_ecalEnergy); redTree->Branch("electron1_eSCPresh ", &electron1_eSCPresh); redTree->Branch("electron1_HoE", &electron1_HoE); redTree->Branch("electron1_dEtaIn", &electron1_dEtaIn); redTree->Branch("electron1_dEtaOut", &electron1_dEtaOut); redTree->Branch("electron1_dPhiIn", &electron1_dPhiIn); redTree->Branch("electron1_dPhiOut", &electron1_dPhiOut); redTree->Branch("electron1_r9", &electron1_r9); redTree->Branch("electron1_charge", &electron1_charge); redTree->Branch("electron1_sigmaEtaEta", &electron1_sigmaEtaEta); redTree->Branch("electron1_sigmaPhiPhi", &electron1_sigmaPhiPhi); redTree->Branch("electron1_etaSC", &electron1_etaSC); redTree->Branch("electron1_nSeed ", &electron1_nSeed); redTree->Branch("electron1_eOverPIn", &electron1_eOverPIn); redTree->Branch("electron1_eOverPOut", &electron1_eOverPOut); redTree->Branch("electron1_ip3d ", &electron1_ip3d); redTree->Branch("electron1_eSeed ", &electron1_eSeed); redTree->Branch("eletrcon1_phiSC ", &eletrcon1_phiSC); redTree->Branch("electron1_trkIso ", &electron1_trkIso); redTree->Branch("electron1_trkIso04", &electron1_trkIso04); redTree->Branch("electron1_d0", &electron1_d0); redTree->Branch("electron1_eSeedOverPIn", &electron1_eSeedOverPIn); redTree->Branch("electron1_eSeedOverPOut", &electron1_eSeedOverPOut); redTree->Branch("electron1_ecalEnergyError", &electron1_ecalEnergyError); redTree->Branch("electron1_ecalIso", &electron1_ecalIso); redTree->Branch("electron1_ecalIso04", &electron1_ecalIso04); redTree->Branch("electron1_etaErr", &electron1_etaErr); redTree->Branch("electron1_etaSCwidth", &electron1_etaSCwidth); redTree->Branch("electron1_hcalDepth1OverEcal", &electron1_hcalDepth1OverEcal); redTree->Branch("electron1_hcalDepth1TowerSumEt", &electron1_hcalDepth1TowerSumEt); redTree->Branch("electron1_hcalDepth1TowerSumEt04", &electron1_hcalDepth1TowerSumEt04); redTree->Branch("electron1_hcalDepth2OverEcal", &electron1_hcalDepth2OverEcal); redTree->Branch("electron1_hcalDepth2TowerSumEt", &electron1_hcalDepth2TowerSumEt); redTree->Branch("electron1_hcalDepth2TowerSumEt04", &electron1_hcalDepth2TowerSumEt04); redTree->Branch("electron1_hcalIso ", &electron1_hcalIso); redTree->Branch("electron1_hcalIso04", &electron1_hcalIso04); redTree->Branch("electron1_ip3derr", &electron1_ip3derr); redTree->Branch("electron1_iso03_pf", &electron1_iso03_pf); redTree->Branch("electron1_iso03_pf2012_ch", &electron1_iso03_pf2012_ch); redTree->Branch("electron1_iso03_pf2012_em ", &electron1_iso03_pf2012_em); redTree->Branch("electron1_iso03_pf2012_nh", &electron1_iso03_pf2012_nh); redTree->Branch("electron1_iso03_pf_ch", &electron1_iso03_pf_ch); redTree->Branch("electron1_iso03_pf_gamma05", &electron1_iso03_pf_gamma05); redTree->Branch("electron1_iso03_pf_nhad05", &electron1_iso03_pf_nhad05); redTree->Branch("electron1_iso04_pf", &electron1_iso04_pf); redTree->Branch("electron1_iso04_pf2012_ch", &electron1_iso04_pf2012_ch); redTree->Branch("electron1_iso04_pf2012_em", &electron1_iso04_pf2012_em); redTree->Branch("electron1_iso04_pf2012_nh", &electron1_iso04_pf2012_nh); redTree->Branch("electron1_iso04_pf_ch", &electron1_iso04_pf_ch); redTree->Branch("electron1_iso04_pf_gamma05", &electron1_iso04_pf_gamma05); redTree->Branch("electron1_iso04_pf_nhad05", &electron1_iso04_pf_nhad05); redTree->Branch("electron1_layer1_charge", &electron1_layer1_charge); redTree->Branch("electron1_lh ", &electron1_lh); redTree->Branch("electron1_mva ", &electron1_mva); redTree->Branch("electron1_ndof", &electron1_ndof); redTree->Branch("electron1_pfChargedHadronIso", &electron1_pfChargedHadronIso); redTree->Branch("electron1_pfNeutralHadronIso", &electron1_pfNeutralHadronIso); redTree->Branch("electron1_pfPhotonIso", &electron1_pfPhotonIso); redTree->Branch("electron1_phiErr", &electron1_phiErr); redTree->Branch("electron1_phiSCwidth", &electron1_phiSCwidth); redTree->Branch("electron1_ptErr", &electron1_ptErr); redTree->Branch("electron1_sigmaIEtaIEta ", &electron1_sigmaIEtaIEta); redTree->Branch("electron1_sigmaIEtaIEtaSC", &electron1_sigmaIEtaIEtaSC); redTree->Branch("electron1_sigmaIEtaIPhi", &electron1_sigmaIEtaIPhi); redTree->Branch("electron1_sigmaIPhiIPhi", &electron1_sigmaIPhiIPhi); redTree->Branch("electron1_sigmaIPhiIPhiSC", &electron1_sigmaIPhiIPhiSC); redTree->Branch("electron1_trackMomentumError", &electron1_trackMomentumError); redTree->Branch("electron1_trkdr", &electron1_trkdr); redTree->Branch("electron1_trkshFrac", &electron1_trkshFrac); redTree->Branch("electron1_z0Err", &electron1_z0Err); redTree->Branch("electron1_z0corr", &electron1_z0corr); redTree->Branch("electron1_mc3_id", &electron1_mc3_id); redTree->Branch("electron1_mc3idx", &electron1_mc3idx); redTree->Branch("electron1_mc3_motherid", &electron1_mc3_motherid); redTree->Branch("electron1_mc3_motheridx", &electron1_mc3_motheridx); redTree->Branch("electron1_mc_id", &electron1_mc_id); redTree->Branch("electron1_mcidx", &electron1_mcidx); redTree->Branch("electron1_mc_motherid", &electron1_mc_motherid); redTree->Branch("electron1_closestJet", &electron1_closestJet); redTree->Branch("electron1_closestMuon", &electron1_closestMuon); // redTree->Branch("electron1_pfelectron0idx", &electron1_pfelectron0idx); redTree->Branch("electron1_category", &electron1_category); redTree->Branch("electron1_ckf_laywithmeas", &electron1_ckf_laywithmeas); redTree->Branch("electron1_conv_delMissHits", &electron1_conv_delMissHits); redTree->Branch("electron1_conv_flag", &electron1_conv_flag); redTree->Branch("electron1_conv_gsftkidx", &electron1_conv_gsftkidx); redTree->Branch("electron1_conv_old_delMissHits", &electron1_conv_old_delMissHits); redTree->Branch("electron1_conv_old_flag ", &electron1_conv_old_flag); // redTree->Branch("electron1_conv_old_gsftkidx ", &electron1_conv_old_gsftkidx); // redTree->Branch("electron1_conv_old_tkidx", &electron1_conv_old_tkidx ); redTree->Branch("electron1_conv_tkidx", &electron1_conv_tkidx); redTree->Branch("electron1_exp_innerlayers", &electron1_exp_innerlayers); redTree->Branch("electron1_exp_outerlayers", &electron1_exp_outerlayers); redTree->Branch("electron1_fiduciality", &electron1_fiduciality); redTree->Branch("electron1_gsftrkidx", &electron1_gsftrkidx); redTree->Branch("electron1_layer1_det", &electron1_layer1_det); redTree->Branch("electron1_layer1_layer", &electron1_layer1_layer); redTree->Branch("electron1_layer1_sizerphi", &electron1_layer1_sizerphi); redTree->Branch("electron1_layer1_sizerz", &electron1_layer1_sizerz ); redTree->Branch("electron1_lost_pixelhits", &electron1_lost_pixelhits ); redTree->Branch("electron1_sccharge", &electron1_sccharge); redTree->Branch("electron1_scindex", &electron1_scindex); redTree->Branch("electron1_trk_charge", &electron1_trk_charge); redTree->Branch("electron1_trkidx", &electron1_trkidx); redTree->Branch("electron1_type", &electron1_type); redTree->Branch("electron1_validHits", &electron1_validHits); redTree->Branch("electron1_valid_pixelhits", &electron1_valid_pixelhits); redTree->Branch("electron1_calJuraTowerIso", &electron1_ecalJuraTowerIso ); redTree->Branch("electron1_calConeIso", &electron1_hcalConeIso ); redTree->Branch("electron1_kJuraIso ", &electron1_tkJuraIso); redTree->Branch("electron1_etdr", &electron1_jetdr); redTree->Branch("electron1_usdr", &electron1_musdr); redTree->Branch("electron1_soR03_chpf_radial", &electron1_isoR03_chpf_radial); redTree->Branch("electron1_soR03_chpf_radialTight", &electron1_isoR03_chpf_radialTight); redTree->Branch("electron1_soR03_chpf_radialTight_bv", &electron1_isoR03_chpf_radialTight_bv); redTree->Branch("electron1_soR03_chpf_radial_bv", &electron1_isoR03_chpf_radial_bv); redTree->Branch("electron1_soR03_empf_radial", &electron1_isoR03_empf_radial); redTree->Branch("electron1_soR03_empf_radialTight", &electron1_isoR03_empf_radialTight); redTree->Branch("electron1_soR03_empf_radialTight_bv", &electron1_isoR03_empf_radialTight_bv); redTree->Branch("electron1_isoR03_empf_radial_bv", &electron1_isoR03_empf_radial_bv); redTree->Branch("electron1_isoR03_nhpf_radial", &electron1_isoR03_nhpf_radial); redTree->Branch("electron1_isoR03_nhpf_radialTight", &electron1_isoR03_nhpf_radialTight); redTree->Branch("electron1_isoR03_nhpf_radialTight_bv", &electron1_isoR03_nhpf_radialTight_bv); redTree->Branch("electron1_isoR03_nhpf_radial_bv", &electron1_isoR03_nhpf_radial_bv); redTree->Branch("electron1_isoR03_pf2012n0p5_ch", &electron1_isoR03_pf2012n0p5_ch); redTree->Branch("electron1_isoR03_pf2012n0p5_em", &electron1_isoR03_pf2012n0p5_em); redTree->Branch("electron1_isoR03_pf2012n0p5_nh", &electron1_isoR03_pf2012n0p5_nh); redTree->Branch("electron1_isoR03_pf_radial", &electron1_isoR03_pf_radial); redTree->Branch("electron1_isoR03_pf_radialTight", &electron1_isoR03_pf_radialTight); redTree->Branch("electron1_isoR03_pf_radialTight_bv", &electron1_isoR03_pf_radialTight_bv); redTree->Branch("electron1_isoR03_pf_radial_bv", &electron1_isoR03_pf_radial_bv); redTree->Branch("electron1_chi2", &electron1_chi2); redTree->Branch("electron1_ckf_chi2", &electron1_ckf_chi2); redTree->Branch("electron1_ckf_ndof", &electron1_ckf_ndof); redTree->Branch("electron1_conv_dcot", &electron1_conv_dcot); redTree->Branch("electron1_conv_dist", &electron1_conv_dist); redTree->Branch("electron1_conv_old_dcot", &electron1_conv_old_dcot); redTree->Branch("electron1_conv_old_dist", &electron1_conv_old_dist); redTree->Branch("electron1_conv_old_radius", &electron1_conv_old_radius); redTree->Branch("electron1_conv_radius", &electron1_conv_radius); redTree->Branch("electron1_d0Err", &electron1_d0Err); redTree->Branch("electron1_d0corr", &electron1_d0corr); redTree->Branch("electron1_dPhiInPhiOut", &electron1_dPhiInPhiOut); redTree->Branch("electron1_deltaEtaEleClusterTrackAtCal", &electron1_deltaEtaEleClusterTrackAtCalo); redTree->Branch("electron1_deltaPhiEleClusterTrackAtCal", &electron1_deltaPhiEleClusterTrackAtCalo); set_goodrun_file("goodruns.txt"); TChain *treeHandle = new TChain("Events"); while (!fileReader.eof()) { reader.ReadLine(fileReader); if( reader.Length() < 2 ) continue; cout << " Adding " << reader << endl; treeHandle->Add( reader ); } Long64_t nentries = treeHandle->GetEntries(); cout << " Number of Entries : " << nentries << endl; // nentries = 1000; cms2.Init(treeHandle); for (Long64_t ientry = 0; ientry < nentries; ientry++) { CMS2::progress(ientry, nentries); cms2.GetEntry(ientry); // clean up data if (!cleaning_goodVertexApril2011()) continue; if (evt_isRealData() && !goodrun(evt_run(), evt_lumiBlock())) continue; dielectron_p4 = LorentzVector(0, 0, 0, 0); dielectron_mass = -999.; dielectron_pt = -999.; scale1fb = -999.; run = -999; event = -999; lumi = -999; //--------------------------------------electron0----------------------------------------// electron0_p4In = LorentzVector(0,0,0,0); electron0_p4Out = LorentzVector(0,0,0,0); electron0_trk_p4 = LorentzVector(0,0,0,0); electron0_vrtx_p4 = LorentzVector(0,0,0,0); electron0_p4 = LorentzVector(0, 0, 0, 0); electron0_class = -999; electron0_fbrem = -999.; electron0_z0 = -999.; electron0_eSC = -999.; electron0_lostHits = -999.; electron0_ecalEnergy = -999.; electron0_eSCPresh = -999.; electron0_HoE = -999.; electron0_dEtaIn = -999.; electron0_dEtaOut = -999.; electron0_dPhiIn = -999.; electron0_dPhiOut = -999.; electron0_r9 = -999.; electron0_charge = -999; electron0_sigmaEtaEta = -999.; electron0_sigmaPhiPhi = -999.; electron0_etaSC = -999.; electron0_nSeed = -999.; electron0_eOverPIn = -999.; electron0_eOverPOut = -999.; electron0_ip3d = -999.; electron0_eSeed = -999.; eletrcon0_phiSC = -999.; electron0_trkIso = -999.; electron0_trkIso04 = -999.; electron0_d0 = -999.; electron0_eSeedOverPIn = -999.; electron0_eSeedOverPOut = -999.; electron0_ecalEnergyError = -999.; electron0_ecalIso = -999.; electron0_ecalIso04 = -999.; electron0_etaErr = -999.; electron0_etaSCwidth = -999.; electron0_hcalDepth1OverEcal = -999.; electron0_hcalDepth1TowerSumEt = -999.; electron0_hcalDepth1TowerSumEt04 = -999.; electron0_hcalDepth2OverEcal = -999.; electron0_hcalDepth2TowerSumEt = -999.; electron0_hcalDepth2TowerSumEt04 = -999.; electron0_hcalIso = -999.; electron0_hcalIso04 = -999.; electron0_ip3derr = -999.; electron0_iso03_pf = -999.; electron0_iso03_pf2012_ch = -999.; electron0_iso03_pf2012_em = -999.; electron0_iso03_pf2012_nh = -999.; electron0_iso03_pf_ch = -999.; electron0_iso03_pf_gamma05 = -999.; electron0_iso03_pf_nhad05 = -999.; electron0_iso04_pf = -999.; electron0_iso04_pf2012_ch = -999.; electron0_iso04_pf2012_em = -999.; electron0_iso04_pf2012_nh = -999.; electron0_iso04_pf_ch = -999.; electron0_iso04_pf_gamma05 = -999.; electron0_iso04_pf_nhad05 = -999.; electron0_layer1_charge = -999.; electron0_lh = -999.; electron0_mva = -999.; electron0_ndof = -999.; electron0_pfChargedHadronIso = -999.; electron0_pfNeutralHadronIso = -999.; electron0_pfPhotonIso = -999.; electron0_phiErr = -999.; electron0_phiSCwidth = -999.; electron0_ptErr = -999.; electron0_sigmaIEtaIEta = -999.; electron0_sigmaIEtaIEtaSC = -999.; electron0_sigmaIEtaIPhi = -999.; electron0_sigmaIPhiIPhi = -999.; electron0_sigmaIPhiIPhiSC = -999.; electron0_trackMomentumError = -999.; electron0_trkdr = -999.; electron0_trkshFrac = -999.; electron0_z0Err = -999.; electron0_z0corr = -999.; electron0_mc3_id =-999; electron0_mc3idx =-999; electron0_mc3_motherid =-999; electron0_mc3_motheridx =-999; electron0_mc_id =-999; electron0_mcidx =-999; electron0_mc_motherid =-999; electron0_closestJet =-999; electron0_closestMuon =-999; // electron0_pfelectron0idx =-999; electron0_category =-999; electron0_ckf_laywithmeas =-999; electron0_conv_delMissHits =-999; electron0_conv_flag =-999; electron0_conv_gsftkidx =-999; electron0_conv_old_delMissHits =-999; electron0_conv_old_flag =-999; // electron0_conv_old_gsftkidx =-999; // electron0_conv_old_tkidx =-999; electron0_conv_tkidx =-999; electron0_exp_innerlayers =-999; electron0_exp_outerlayers =-999; electron0_fiduciality =-999; electron0_gsftrkidx =-999; electron0_layer1_det =-999; electron0_layer1_layer =-999; electron0_layer1_sizerphi =-999; electron0_layer1_sizerz =-999; electron0_lost_pixelhits =-999; electron0_sccharge =-999; electron0_scindex =-999; electron0_trk_charge =-999; electron0_trkidx =-999; electron0_type =-999; electron0_validHits =-999; electron0_valid_pixelhits =-999; electron0_ecalJuraTowerIso =-999.; electron0_hcalConeIso =-999.; electron0_tkJuraIso =-999.; electron0_jetdr =-999.; electron0_musdr =-999.; electron0_isoR03_chpf_radial =-999.; electron0_isoR03_chpf_radialTight =-999.; electron0_isoR03_chpf_radialTight_bv =-999; electron0_isoR03_chpf_radial_bv =-999; electron0_isoR03_empf_radial =-999.; electron0_isoR03_empf_radialTight =-999.; electron0_isoR03_empf_radialTight_bv =-999.; electron0_isoR03_empf_radial_bv =-999.; electron0_isoR03_nhpf_radial =-999.; electron0_isoR03_nhpf_radialTight =-999.; electron0_isoR03_nhpf_radialTight_bv =-999.; electron0_isoR03_nhpf_radial_bv =-999.; electron0_isoR03_pf2012n0p5_ch =-999.; electron0_isoR03_pf2012n0p5_em =-999.; electron0_isoR03_pf2012n0p5_nh =-999.; electron0_isoR03_pf_radial =-999.; electron0_isoR03_pf_radialTight =-999.; electron0_isoR03_pf_radialTight_bv =-999.; electron0_isoR03_pf_radial_bv =-999.; electron0_chi2 =-999.; electron0_ckf_chi2 =-999.; electron0_ckf_ndof =-999.; electron0_conv_dcot =-999.; electron0_conv_dist =-999.; electron0_conv_old_dcot =-999.; electron0_conv_old_dist =-999.; electron0_conv_old_radius =-999.; electron0_conv_radius =-999.; electron0_d0Err =-999.; electron0_d0corr =-999.; electron0_dPhiInPhiOut =-999.; electron0_deltaEtaEleClusterTrackAtCalo =-999.; electron0_deltaPhiEleClusterTrackAtCalo =-999.; //--------------------------------------electron1----------------------------------------// electron1_p4In = LorentzVector(0,0,0,0); electron1_p4Out = LorentzVector(0,0,0,0); electron1_trk_p4 = LorentzVector(0,0,0,0); electron1_vrtx_p4 = LorentzVector(0,0,0,0); electron1_p4 = LorentzVector(0, 0, 0, 0); electron1_class = -999; electron1_fbrem = -999.; electron1_z0 = -999.; electron1_eSC = -999.; electron1_lostHits = -999.; electron1_ecalEnergy = -999.; electron1_eSCPresh = -999.; electron1_HoE = -999.; electron1_dEtaIn = -999.; electron1_dEtaOut = -999.; electron1_dPhiIn = -999.; electron1_dPhiOut = -999.; electron1_r9 = -999.; electron1_charge = -999; electron1_sigmaEtaEta = -999.; electron1_sigmaPhiPhi = -999.; electron1_etaSC = -999.; electron1_nSeed = -999.; electron1_eOverPIn = -999.; electron1_eOverPOut = -999.; electron1_ip3d = -999.; electron1_eSeed = -999.; eletrcon1_phiSC = -999.; electron1_trkIso = -999.; electron1_trkIso04 = -999.; electron1_d0 = -999.; electron1_eSeedOverPIn = -999.; electron1_eSeedOverPOut = -999.; electron1_ecalEnergyError = -999.; electron1_ecalIso = -999.; electron1_ecalIso04 = -999.; electron1_etaErr = -999.; electron1_etaSCwidth = -999.; electron1_hcalDepth1OverEcal = -999.; electron1_hcalDepth1TowerSumEt = -999.; electron1_hcalDepth1TowerSumEt04 = -999.; electron1_hcalDepth2OverEcal = -999.; electron1_hcalDepth2TowerSumEt = -999.; electron1_hcalDepth2TowerSumEt04 = -999.; electron1_hcalIso = -999.; electron1_hcalIso04 = -999.; electron1_ip3derr = -999.; electron1_iso03_pf = -999.; electron1_iso03_pf2012_ch = -999.; electron1_iso03_pf2012_em = -999.; electron1_iso03_pf2012_nh = -999.; electron1_iso03_pf_ch = -999.; electron1_iso03_pf_gamma05 = -999.; electron1_iso03_pf_nhad05 = -999.; electron1_iso04_pf = -999.; electron1_iso04_pf2012_ch = -999.; electron1_iso04_pf2012_em = -999.; electron1_iso04_pf2012_nh = -999.; electron1_iso04_pf_ch = -999.; electron1_iso04_pf_gamma05 = -999.; electron1_iso04_pf_nhad05 = -999.; electron1_layer1_charge = -999.; electron1_lh = -999.; electron1_mva = -999.; electron1_ndof = -999.; electron1_pfChargedHadronIso = -999.; electron1_pfNeutralHadronIso = -999.; electron1_pfPhotonIso = -999.; electron1_phiErr = -999.; electron1_phiSCwidth = -999.; electron1_ptErr = -999.; electron1_sigmaIEtaIEta = -999.; electron1_sigmaIEtaIEtaSC = -999.; electron1_sigmaIEtaIPhi = -999.; electron1_sigmaIPhiIPhi = -999.; electron1_sigmaIPhiIPhiSC = -999.; electron1_trackMomentumError = -999.; electron1_trkdr = -999.; electron1_trkshFrac = -999.; electron1_z0Err = -999.; electron1_z0corr = -999.; electron1_mc3_id =-999; electron1_mc3idx =-999; electron1_mc3_motherid =-999; electron1_mc3_motheridx =-999; electron1_mc_id =-999; electron1_mcidx =-999; electron1_mc_motherid =-999; electron1_closestJet =-999; electron1_closestMuon =-999; /// electron1_pfelectron0idx =-999; electron1_category =-999; electron1_ckf_laywithmeas =-999; electron1_conv_delMissHits =-999; electron1_conv_flag =-999; electron1_conv_gsftkidx =-999; electron1_conv_old_delMissHits =-999; electron1_conv_old_flag =-999; // electron1_conv_old_gsftkidx =-999; // electron1_conv_old_tkidx =-999; electron1_conv_tkidx =-999; electron1_exp_innerlayers =-999; electron1_exp_outerlayers =-999; electron1_fiduciality =-999; electron1_gsftrkidx =-999; electron1_layer1_det =-999; electron1_layer1_layer =-999; electron1_layer1_sizerphi =-999; electron1_layer1_sizerz =-999; electron1_lost_pixelhits =-999; electron1_sccharge =-999; electron1_scindex =-999; electron1_trk_charge =-999; electron1_trkidx =-999; electron1_type =-999; electron1_validHits =-999; electron1_valid_pixelhits =-999; electron1_ecalJuraTowerIso =-999.; electron1_hcalConeIso =-999.; electron1_tkJuraIso =-999.; electron1_jetdr =-999.; electron1_musdr =-999.; electron1_isoR03_chpf_radial =-999.; electron1_isoR03_chpf_radialTight =-999.; electron1_isoR03_chpf_radialTight_bv =-999; electron1_isoR03_chpf_radial_bv =-999; electron1_isoR03_empf_radial =-999.; electron1_isoR03_empf_radialTight =-999.; electron1_isoR03_empf_radialTight_bv =-999.; electron1_isoR03_empf_radial_bv =-999.; electron1_isoR03_nhpf_radial =-999.; electron1_isoR03_nhpf_radialTight =-999.; electron1_isoR03_nhpf_radialTight_bv =-999.; electron1_isoR03_nhpf_radial_bv =-999.; electron1_isoR03_pf2012n0p5_ch =-999.; electron1_isoR03_pf2012n0p5_em =-999.; electron1_isoR03_pf2012n0p5_nh =-999.; electron1_isoR03_pf_radial =-999.; electron1_isoR03_pf_radialTight =-999.; electron1_isoR03_pf_radialTight_bv =-999.; electron1_isoR03_pf_radial_bv =-999.; electron1_chi2 =-999.; electron1_ckf_chi2 =-999.; electron1_ckf_ndof =-999.; electron1_conv_dcot =-999.; electron1_conv_dist =-999.; electron1_conv_old_dcot =-999.; electron1_conv_old_dist =-999.; electron1_conv_old_radius =-999.; electron1_conv_radius =-999.; electron1_d0Err =-999.; electron1_d0corr =-999.; electron1_dPhiInPhiOut =-999.; electron1_deltaEtaEleClusterTrackAtCalo =-999.; electron1_deltaPhiEleClusterTrackAtCalo =-999.; int electron0_idx = -1; int electron1_idx = -1; // Analysis selection for (int ielec = 2; ielec < els_p4().size(); ielec++) { float pt0; float pt1; // cout << els_p4().size() << endl; if (electron0_idx == -1) { // first loop electron0_idx = 0; // the higher PT electron1_idx = 1; // the samller PT pt0 = els_p4().at(electron0_idx).pt(); pt1 = els_p4().at(electron1_idx).pt(); if (pt0 < pt1) swap(electron0_idx, electron0_idx); } float pt_current = els_p4().at(ielec).pt(); if (pt_current < pt1) // nothing to do continue; if (pt_current > pt1 && pt_current <= pt0) { // between the two values electron1_idx = ielec; pt1 = els_p4().at(electron1_idx).pt(); } else { // greater than both values electron1_idx = electron0_idx; pt1 = pt0; electron0_idx = ielec; pt0 = pt_current; } // check the electrons agisnt the four tests: TIGHT, LOOSE, etc ... // bool electron0_pass = false; //bool electron1_pass = false; /* if (electron0_idx != -1 && electron1_idx != -1) { electron0_pass = electronId_WP2012_v3(electron0_idx, TIGHT) == PassAllWP2012Cuts; electron1_pass = electronId_WP2012_v3(electron1_idx, TIGHT) == PassAllWP2012Cuts; } electronId_WP2012_v3(ielec, TIGHT) == PassAllWP2012Cuts; electron0_pass = true; electron1_pass = true;*/ if (electron0_idx != -1 && electron1_idx != -1 && electronId_WP2012_v3(electron0_idx, TIGHT) == PassAllWP2012Cuts && electronId_WP2012_v3(electron1_idx, TIGHT) == PassAllWP2012Cuts) { electron0_pass[0] = true; electron1_pass[0] = true; } else { electron0_pass[0] = false; electron1_pass[0] = false; } if (electron0_idx != -1 && electron1_idx != -1 && electronId_WP2012_v3(electron0_idx, MEDIUM) == PassAllWP2012Cuts && electronId_WP2012_v3(electron1_idx, MEDIUM) == PassAllWP2012Cuts) { electron0_pass[1] = true; electron1_pass[1] = true; } else { electron0_pass[1] = false; electron1_pass[1] = false; } if (electron0_idx != -1 && electron1_idx != -1 && electronId_WP2012_v3(electron0_idx, LOOSE) == PassAllWP2012Cuts && electronId_WP2012_v3(electron1_idx, LOOSE) == PassAllWP2012Cuts) { electron0_pass[2] = true; electron1_pass[2] = true; } else { electron0_pass[2] = false; electron1_pass[2] = false; } if (electron0_idx != -1 && electron1_idx != -1 && electronId_WP2012_v3(electron0_idx, VETO) == PassAllWP2012Cuts && electronId_WP2012_v3(electron1_idx, VETO) == PassAllWP2012Cuts) { electron0_pass[3] = true; electron1_pass[3] = true; } else { electron0_pass[3] = false; electron1_pass[3] = false; } } if (electron0_idx < 0 || electron1_idx < 0) continue; scale1fb = evt_scale1fb(); run = evt_run(); event = evt_event(); lumi = evt_lumiBlock(); //------------------------------------------electron0-----------------------------------------------------// electron0_p4In = els_p4In().at(electron0_idx); electron0_p4Out = els_p4Out().at(electron0_idx); electron0_trk_p4 = els_trk_p4().at(electron0_idx); electron0_vrtx_p4 = els_vertex_p4().at(electron0_idx); electron0_p4 = els_p4().at(electron0_idx); electron0_class = els_class().at(electron0_idx); electron0_fbrem = els_fbrem().at(electron0_idx); electron0_eSC = els_eSC().at(electron0_idx); electron0_z0 = els_z0().at(electron0_idx); electron0_lostHits = els_lostHits().at(electron0_idx); electron0_ecalEnergy = els_ecalEnergy().at(electron0_idx); electron0_eSCPresh = els_eSCPresh().at(electron0_idx); electron0_HoE = els_hOverE().at(electron0_idx); electron0_dEtaIn = els_dEtaIn().at(electron0_idx); electron0_dEtaOut = els_dEtaOut().at(electron0_idx); electron0_dPhiIn = els_dPhiIn().at(electron0_idx); electron0_dPhiOut = els_dPhiOut().at(electron0_idx); electron0_r9 = els_r9().at(electron0_idx); electron0_charge = els_charge().at(electron0_idx); electron0_sigmaEtaEta = els_sigmaEtaEta().at(electron0_idx); electron0_sigmaPhiPhi = els_sigmaPhiPhi().at(electron0_idx); electron0_etaSC = els_etaSC().at(electron0_idx); electron0_nSeed = els_nSeed().at(electron0_idx); electron0_eOverPIn = els_eOverPIn().at(electron0_idx); electron0_eOverPOut = els_eOverPOut().at(electron0_idx); electron0_ip3d = els_ip3d().at(electron0_idx); electron0_eSeed = els_eSeed().at(electron0_idx); eletrcon0_phiSC = els_phiSC().at(electron0_idx); electron0_trkIso = els_tkIso().at(electron0_idx); electron0_trkIso04 = els_tkIso04().at(electron0_idx); electron0_d0 = els_d0().at(electron0_idx); electron0_eSeedOverPIn = els_eSeedOverPIn().at(electron0_idx); electron0_eSeedOverPOut = els_eSeedOverPOut().at(electron0_idx); electron0_ecalEnergyError = els_ecalEnergyError().at(electron0_idx); electron0_ecalIso = els_ecalIso().at(electron0_idx); electron0_ecalIso04 = els_ecalIso04().at(electron0_idx); electron0_etaErr = els_etaErr().at(electron0_idx); electron0_etaSCwidth = els_etaSCwidth().at(electron0_idx); electron0_hcalDepth1OverEcal = els_hcalDepth1OverEcal().at(electron0_idx); electron0_hcalDepth1TowerSumEt = els_hcalDepth1TowerSumEt().at(electron0_idx); electron0_hcalDepth1TowerSumEt04 = els_hcalDepth1TowerSumEt04().at(electron0_idx); electron0_hcalDepth2OverEcal = els_hcalDepth2OverEcal().at(electron0_idx); electron0_hcalDepth2TowerSumEt = els_hcalDepth2TowerSumEt().at(electron0_idx); electron0_hcalDepth2TowerSumEt04 = els_hcalDepth2TowerSumEt04().at(electron0_idx); electron0_hcalIso = els_hcalIso().at(electron0_idx); electron0_hcalIso04 = els_hcalIso04().at(electron0_idx); electron0_ip3derr = els_ip3derr().at(electron0_idx); electron0_iso03_pf = els_iso03_pf().at(electron0_idx); electron0_iso03_pf2012_ch = els_iso03_pf2012_ch().at(electron0_idx); electron0_iso03_pf2012_em = els_iso03_pf2012_em().at(electron0_idx); electron0_iso03_pf2012_nh = els_iso03_pf2012_nh().at(electron0_idx); electron0_iso03_pf_ch = els_iso03_pf_ch().at(electron0_idx); electron0_iso03_pf_gamma05 = els_iso03_pf_gamma05().at(electron0_idx); electron0_iso03_pf_nhad05 = els_iso03_pf_nhad05().at(electron0_idx); electron0_iso04_pf = els_iso04_pf().at(electron0_idx); electron0_iso04_pf2012_ch = els_iso04_pf2012_ch().at(electron0_idx); electron0_iso04_pf2012_em = els_iso04_pf2012_em().at(electron0_idx); electron0_iso04_pf2012_nh = els_iso04_pf2012_nh().at(electron0_idx); electron0_iso04_pf_ch = els_iso04_pf_ch().at(electron0_idx); electron0_iso04_pf_gamma05 = els_iso04_pf_gamma05().at(electron0_idx); electron0_iso04_pf_nhad05 = els_iso04_pf_nhad05().at(electron0_idx); electron0_layer1_charge = els_layer1_charge().at(electron0_idx); electron0_lh = els_lh().at(electron0_idx); electron0_mva = els_mva().at(electron0_idx); electron0_ndof = els_ndof().at(electron0_idx); electron0_pfChargedHadronIso = els_pfChargedHadronIso().at(electron0_idx); electron0_pfNeutralHadronIso = els_pfNeutralHadronIso().at(electron0_idx); electron0_pfPhotonIso = els_pfPhotonIso().at(electron0_idx); electron0_phiErr = els_phiErr().at(electron0_idx); electron0_phiSCwidth = els_phiSCwidth().at(electron0_idx); electron0_ptErr = els_ptErr().at(electron0_idx); electron0_sigmaIEtaIEta = els_sigmaIEtaIEta().at(electron0_idx); electron0_sigmaIEtaIEtaSC = els_sigmaIEtaIEtaSC().at(electron0_idx); electron0_sigmaIEtaIPhi = els_sigmaIEtaIPhi().at(electron0_idx); electron0_sigmaIPhiIPhi = els_sigmaIPhiIPhi().at(electron0_idx); electron0_sigmaIPhiIPhiSC = els_sigmaIPhiIPhiSC().at(electron0_idx); electron0_trackMomentumError = els_trackMomentumError().at(electron0_idx); electron0_trkdr = els_trkdr().at(electron0_idx); electron0_trkshFrac = els_trkshFrac().at(electron0_idx); electron0_z0Err = els_z0Err().at(electron0_idx); electron0_z0corr = els_z0corr().at(electron0_idx); electron0_mc_id = els_mc_id().at(electron0_idx); electron0_mcidx = els_mcidx().at(electron0_idx); electron0_mc_motherid = els_mc_motherid().at(electron0_idx); electron0_closestJet = els_closestJet().at(electron0_idx); electron0_closestMuon = els_closestMuon().at(electron0_idx); ///electron0_pfelectron0idx = els_pfelsidx.at(electron0_idx); electron0_category = els_category().at(electron0_idx); electron0_ckf_laywithmeas = els_ckf_laywithmeas().at(electron0_idx); electron0_conv_delMissHits = els_conv_delMissHits().at(electron0_idx); electron0_conv_flag = els_conv_flag().at(electron0_idx); electron0_conv_gsftkidx = els_conv_gsftkidx().at(electron0_idx); electron0_conv_old_delMissHits = els_conv_old_delMissHits().at(electron0_idx); electron0_conv_old_flag = els_conv_old_flag().at(electron0_idx); //electron0_conv_old_gsftkidx = els_conv_old_gsftkidx().at(electron0_idx); //electron0_conv_old_tkidx = els_conv_old_tkidx().at(electron0_idx); electron0_conv_tkidx = els_conv_tkidx().at(electron0_idx); electron0_exp_innerlayers = els_exp_innerlayers().at(electron0_idx); electron0_exp_outerlayers = els_exp_outerlayers().at(electron0_idx); electron0_fiduciality = els_fiduciality().at(electron0_idx); electron0_gsftrkidx = els_gsftrkidx().at(electron0_idx); electron0_layer1_det = els_layer1_det().at(electron0_idx); electron0_layer1_layer = els_layer1_layer().at(electron0_idx); electron0_layer1_sizerphi = els_layer1_sizerphi().at(electron0_idx); electron0_layer1_sizerz = els_layer1_sizerz().at(electron0_idx); electron0_lost_pixelhits = els_lost_pixelhits().at(electron0_idx); electron0_sccharge = els_sccharge().at(electron0_idx); electron0_scindex = els_scindex().at(electron0_idx); electron0_trk_charge = els_trk_charge().at(electron0_idx); electron0_trkidx = els_trkidx().at(electron0_idx); electron0_type = els_type().at(electron0_idx); electron0_validHits = els_validHits().at(electron0_idx); electron0_valid_pixelhits = els_valid_pixelhits().at(electron0_idx); electron0_ecalJuraTowerIso = els_ecalJuraTowerIso().at(electron0_idx); electron0_hcalConeIso = els_hcalConeIso().at(electron0_idx); electron0_tkJuraIso = els_tkJuraIso().at(electron0_idx); electron0_jetdr = els_jetdr().at(electron0_idx); electron0_musdr = els_musdr().at(electron0_idx); electron0_isoR03_chpf_radial = els_isoR03_chpf_radial().at(electron0_idx); electron0_isoR03_chpf_radialTight = els_isoR03_chpf_radialTight().at(electron0_idx); electron0_isoR03_chpf_radialTight_bv = els_isoR03_chpf_radialTight_bv().at(electron0_idx); electron0_isoR03_chpf_radial_bv = els_isoR03_chpf_radial_bv().at(electron0_idx); electron0_isoR03_empf_radial = els_isoR03_empf_radial().at(electron0_idx); electron0_isoR03_empf_radialTight = els_isoR03_empf_radialTight().at(electron0_idx); electron0_isoR03_empf_radialTight_bv = els_isoR03_empf_radialTight_bv().at(electron0_idx); electron0_isoR03_empf_radial_bv = els_isoR03_empf_radial_bv().at(electron0_idx); electron0_isoR03_nhpf_radial = els_isoR03_nhpf_radial().at(electron0_idx); electron0_isoR03_nhpf_radialTight = els_isoR03_nhpf_radialTight().at(electron0_idx); electron0_isoR03_nhpf_radialTight_bv = els_isoR03_nhpf_radialTight_bv().at(electron0_idx); electron0_isoR03_nhpf_radial_bv = els_isoR03_nhpf_radial_bv().at(electron0_idx); electron0_isoR03_pf2012n0p5_ch = els_isoR03_pf2012n0p5_ch().at(electron0_idx); electron0_isoR03_pf2012n0p5_em = els_isoR03_pf2012n0p5_em().at(electron0_idx); electron0_isoR03_pf2012n0p5_nh = els_isoR03_pf2012n0p5_nh().at(electron0_idx); electron0_isoR03_pf_radial = els_isoR03_pf_radial().at(electron0_idx); electron0_isoR03_pf_radialTight = els_isoR03_pf_radialTight().at(electron0_idx); electron0_isoR03_pf_radialTight_bv = els_isoR03_pf_radialTight_bv().at(electron0_idx); electron0_isoR03_pf_radial_bv = els_isoR03_pf_radial_bv().at(electron0_idx); electron0_chi2 = els_chi2().at(electron0_idx); electron0_ckf_chi2 = els_ckf_chi2().at(electron0_idx); electron0_ckf_ndof = els_ckf_ndof().at(electron0_idx); electron0_conv_dcot = els_conv_dcot().at(electron0_idx); electron0_conv_dist = els_conv_dist().at(electron0_idx); electron0_conv_old_dcot = els_conv_old_dcot().at(electron0_idx); electron0_conv_old_dist = els_conv_old_dist().at(electron0_idx); electron0_conv_old_radius = els_conv_old_radius().at(electron0_idx); electron0_conv_radius = els_conv_radius().at(electron0_idx); electron0_d0Err = els_d0Err().at(electron0_idx); electron0_d0corr = els_d0corr().at(electron0_idx); electron0_dPhiInPhiOut = els_dPhiInPhiOut().at(electron0_idx); electron0_deltaEtaEleClusterTrackAtCalo = els_deltaEtaEleClusterTrackAtCalo().at(electron0_idx); electron0_deltaPhiEleClusterTrackAtCalo = els_deltaPhiEleClusterTrackAtCalo().at(electron0_idx); //---------------------------------------electron1-------------------------------------------------------// electron1_p4In = els_p4In().at(electron1_idx); electron1_p4Out = els_p4Out().at(electron1_idx); electron1_trk_p4 = els_trk_p4().at(electron1_idx); electron1_vrtx_p4 = els_vertex_p4().at(electron1_idx); electron1_p4 = els_p4().at(electron1_idx); electron1_class = els_class().at(electron1_idx); electron1_fbrem = els_fbrem().at(electron1_idx); electron1_eSC = els_eSC().at(electron1_idx); electron1_z0 = els_z0().at(electron1_idx); electron1_lostHits = els_lostHits().at(electron1_idx); electron1_ecalEnergy = els_ecalEnergy().at(electron1_idx); electron1_eSCPresh = els_eSCPresh().at(electron1_idx); electron1_HoE = els_hOverE().at(electron1_idx); electron1_dEtaIn = els_dEtaIn().at(electron1_idx); electron1_dEtaOut = els_dEtaOut().at(electron1_idx); electron1_dPhiIn = els_dPhiIn().at(electron1_idx); electron1_dPhiOut = els_dPhiOut().at(electron1_idx); electron1_r9 = els_r9().at(electron1_idx); electron1_charge = els_charge().at(electron1_idx); electron1_sigmaEtaEta = els_sigmaEtaEta().at(electron1_idx); electron1_sigmaPhiPhi = els_sigmaPhiPhi().at(electron1_idx); electron1_etaSC = els_etaSC().at(electron1_idx); electron1_nSeed = els_nSeed().at(electron1_idx); electron1_eOverPIn = els_eOverPIn().at(electron1_idx); electron1_eOverPOut = els_eOverPOut().at(electron1_idx); electron1_ip3d = els_ip3d().at(electron1_idx); electron1_eSeed = els_eSeed().at(electron1_idx); eletrcon1_phiSC = els_phiSC().at(electron1_idx); electron1_trkIso = els_tkIso().at(electron1_idx); electron1_trkIso04 = els_tkIso04().at(electron1_idx); electron1_d0 = els_d0().at(electron1_idx); electron1_eSeedOverPIn = els_eSeedOverPIn().at(electron1_idx); electron1_eSeedOverPOut = els_eSeedOverPOut().at(electron1_idx); electron1_ecalEnergyError = els_ecalEnergyError().at(electron1_idx); electron1_ecalIso = els_ecalIso().at(electron1_idx); electron1_ecalIso04 = els_ecalIso04().at(electron1_idx); electron1_etaErr = els_etaErr().at(electron1_idx); electron1_etaSCwidth = els_etaSCwidth().at(electron1_idx); electron1_hcalDepth1OverEcal = els_hcalDepth1OverEcal().at(electron1_idx); electron1_hcalDepth1TowerSumEt = els_hcalDepth1TowerSumEt().at(electron1_idx); electron1_hcalDepth1TowerSumEt04 = els_hcalDepth1TowerSumEt04().at(electron1_idx); electron1_hcalDepth2OverEcal = els_hcalDepth2OverEcal().at(electron1_idx); electron1_hcalDepth2TowerSumEt = els_hcalDepth2TowerSumEt().at(electron1_idx); electron1_hcalDepth2TowerSumEt04 = els_hcalDepth2TowerSumEt04().at(electron1_idx); electron1_hcalIso = els_hcalIso().at(electron1_idx); electron1_hcalIso04 = els_hcalIso04().at(electron1_idx); electron1_ip3derr = els_ip3derr().at(electron1_idx); electron1_iso03_pf = els_iso03_pf().at(electron1_idx); electron1_iso03_pf2012_ch = els_iso03_pf2012_ch().at(electron1_idx); electron1_iso03_pf2012_em = els_iso03_pf2012_em().at(electron1_idx); electron1_iso03_pf2012_nh = els_iso03_pf2012_nh().at(electron1_idx); electron1_iso03_pf_ch = els_iso03_pf_ch().at(electron1_idx); electron1_iso03_pf_gamma05 = els_iso03_pf_gamma05().at(electron1_idx); electron1_iso03_pf_nhad05 = els_iso03_pf_nhad05().at(electron1_idx); electron1_iso04_pf = els_iso04_pf().at(electron1_idx); electron1_iso04_pf2012_ch = els_iso04_pf2012_ch().at(electron1_idx); electron1_iso04_pf2012_em = els_iso04_pf2012_em().at(electron1_idx); electron1_iso04_pf2012_nh = els_iso04_pf2012_nh().at(electron1_idx); electron1_iso04_pf_ch = els_iso04_pf_ch().at(electron1_idx); electron1_iso04_pf_gamma05 = els_iso04_pf_gamma05().at(electron1_idx); electron1_iso04_pf_nhad05 = els_iso04_pf_nhad05().at(electron1_idx); electron1_layer1_charge = els_layer1_charge().at(electron1_idx); electron1_lh = els_lh().at(electron1_idx); electron1_mva = els_mva().at(electron1_idx); electron1_ndof = els_ndof().at(electron1_idx); electron1_pfChargedHadronIso = els_pfChargedHadronIso().at(electron1_idx); electron1_pfNeutralHadronIso = els_pfNeutralHadronIso().at(electron1_idx); electron1_pfPhotonIso = els_pfPhotonIso().at(electron1_idx); electron1_phiErr = els_phiErr().at(electron1_idx); electron1_phiSCwidth = els_phiSCwidth().at(electron1_idx); electron1_ptErr = els_ptErr().at(electron1_idx); electron1_sigmaIEtaIEta = els_sigmaIEtaIEta().at(electron1_idx); electron1_sigmaIEtaIEtaSC = els_sigmaIEtaIEtaSC().at(electron1_idx); electron1_sigmaIEtaIPhi = els_sigmaIEtaIPhi().at(electron1_idx); electron1_sigmaIPhiIPhi = els_sigmaIPhiIPhi().at(electron1_idx); electron1_sigmaIPhiIPhiSC = els_sigmaIPhiIPhiSC().at(electron1_idx); electron1_trackMomentumError = els_trackMomentumError().at(electron1_idx); electron1_trkdr = els_trkdr().at(electron1_idx); electron1_trkshFrac = els_trkshFrac().at(electron1_idx); electron1_z0Err = els_z0Err().at(electron1_idx); electron1_z0corr = els_z0corr().at(electron1_idx); electron1_mc_id = els_mc_id().at(electron1_idx); electron1_mcidx = els_mcidx().at(electron1_idx); electron1_mc_motherid = els_mc_motherid().at(electron1_idx); electron1_closestJet = els_closestJet().at(electron1_idx); electron1_closestMuon = els_closestMuon().at(electron1_idx); ///electron1_pfelectron1idx = els_pfelsidx.at(electron1_idx); electron1_category = els_category().at(electron1_idx); electron1_ckf_laywithmeas = els_ckf_laywithmeas().at(electron1_idx); electron1_conv_delMissHits = els_conv_delMissHits().at(electron1_idx); electron1_conv_flag = els_conv_flag().at(electron1_idx); electron1_conv_gsftkidx = els_conv_gsftkidx().at(electron1_idx); electron1_conv_old_delMissHits = els_conv_old_delMissHits().at(electron1_idx); electron1_conv_old_flag = els_conv_old_flag().at(electron1_idx); //electron1_conv_old_gsftkidx = els_conv_old_gsftkidx().at(electron1_idx); //electron1_conv_old_tkidx = els_conv_old_tkidx().at(electron1_idx); electron1_conv_tkidx = els_conv_tkidx().at(electron1_idx); electron1_exp_innerlayers = els_exp_innerlayers().at(electron1_idx); electron1_exp_outerlayers = els_exp_outerlayers().at(electron1_idx); electron1_fiduciality = els_fiduciality().at(electron1_idx); electron1_gsftrkidx = els_gsftrkidx().at(electron1_idx); electron1_layer1_det = els_layer1_det().at(electron1_idx); electron1_layer1_layer = els_layer1_layer().at(electron1_idx); electron1_layer1_sizerphi = els_layer1_sizerphi().at(electron1_idx); electron1_layer1_sizerz = els_layer1_sizerz().at(electron1_idx); electron1_lost_pixelhits = els_lost_pixelhits().at(electron1_idx); electron1_sccharge = els_sccharge().at(electron1_idx); electron1_scindex = els_scindex().at(electron1_idx); electron1_trk_charge = els_trk_charge().at(electron1_idx); electron1_trkidx = els_trkidx().at(electron1_idx); electron1_type = els_type().at(electron1_idx); electron1_validHits = els_validHits().at(electron1_idx); electron1_valid_pixelhits = els_valid_pixelhits().at(electron1_idx); electron1_ecalJuraTowerIso = els_ecalJuraTowerIso().at(electron1_idx); electron1_hcalConeIso = els_hcalConeIso().at(electron1_idx); electron1_tkJuraIso = els_tkJuraIso().at(electron1_idx); electron1_jetdr = els_jetdr().at(electron1_idx); electron1_musdr = els_musdr().at(electron1_idx); electron1_isoR03_chpf_radial = els_isoR03_chpf_radial().at(electron1_idx); electron1_isoR03_chpf_radialTight = els_isoR03_chpf_radialTight().at(electron1_idx); electron1_isoR03_chpf_radialTight_bv = els_isoR03_chpf_radialTight_bv().at(electron1_idx); electron1_isoR03_chpf_radial_bv = els_isoR03_chpf_radial_bv().at(electron1_idx); electron1_isoR03_empf_radial = els_isoR03_empf_radial().at(electron1_idx); electron1_isoR03_empf_radialTight = els_isoR03_empf_radialTight().at(electron1_idx); electron1_isoR03_empf_radialTight_bv = els_isoR03_empf_radialTight_bv().at(electron1_idx); electron1_isoR03_empf_radial_bv = els_isoR03_empf_radial_bv().at(electron1_idx); electron1_isoR03_nhpf_radial = els_isoR03_nhpf_radial().at(electron1_idx); electron1_isoR03_nhpf_radialTight = els_isoR03_nhpf_radialTight().at(electron1_idx); electron1_isoR03_nhpf_radialTight_bv = els_isoR03_nhpf_radialTight_bv().at(electron1_idx); electron1_isoR03_nhpf_radial_bv = els_isoR03_nhpf_radial_bv().at(electron1_idx); electron1_isoR03_pf2012n0p5_ch = els_isoR03_pf2012n0p5_ch().at(electron1_idx); electron1_isoR03_pf2012n0p5_em = els_isoR03_pf2012n0p5_em().at(electron1_idx); electron1_isoR03_pf2012n0p5_nh = els_isoR03_pf2012n0p5_nh().at(electron1_idx); electron1_isoR03_pf_radial = els_isoR03_pf_radial().at(electron1_idx); electron1_isoR03_pf_radialTight = els_isoR03_pf_radialTight().at(electron1_idx); electron1_isoR03_pf_radialTight_bv = els_isoR03_pf_radialTight_bv().at(electron1_idx); electron1_isoR03_pf_radial_bv = els_isoR03_pf_radial_bv().at(electron1_idx); electron1_chi2 = els_chi2().at(electron1_idx); electron1_ckf_chi2 = els_ckf_chi2().at(electron1_idx); electron1_ckf_ndof = els_ckf_ndof().at(electron1_idx); electron1_conv_dcot = els_conv_dcot().at(electron1_idx); electron1_conv_dist = els_conv_dist().at(electron1_idx); electron1_conv_old_dcot = els_conv_old_dcot().at(electron1_idx); electron1_conv_old_dist = els_conv_old_dist().at(electron1_idx); electron1_conv_old_radius = els_conv_old_radius().at(electron1_idx); electron1_conv_radius = els_conv_radius().at(electron1_idx); electron1_d0Err = els_d0Err().at(electron1_idx); electron1_d0corr = els_d0corr().at(electron1_idx); electron1_dPhiInPhiOut = els_dPhiInPhiOut().at(electron1_idx); electron1_deltaEtaEleClusterTrackAtCalo = els_deltaEtaEleClusterTrackAtCalo().at(electron1_idx); electron1_deltaPhiEleClusterTrackAtCalo = els_deltaPhiEleClusterTrackAtCalo().at(electron1_idx); // cout << electron0_pass << endl; dielectron_p4 = electron0_p4+electron1_p4; dielectron_mass =dielectron_p4.mass(); dielectron_pt = dielectron_p4.pt(); redTree->Fill(); } outputFile->cd(); redTree->Write(); outputFile->Close(); }