#ifndef METSTRUCT_HH #define METSTRUCT_HH //Blocks: met const Int_t MET_CORR_CMUO = 0x01; const Int_t MET_CORR_CMIO = 0x02; const Int_t MET_CORR_BMU = 0x04; const Int_t MET_CORR_EMOBJ = 0x08; const Int_t MET_CORR_CJETS = 0x10; const Int_t MET_CORR_PJETS = 0x20; const Int_t MET_CORR_ALL = 0xffff; // "del" variables = changes to met_x and met_y due to specific corrections struct mymet_t { // For the following variables, the MET is calculated using the // vertex chosen in the vertex block, unless you set the option // revertexMet to false in the tcl file Float_t metraw; // raw met from calorimeter Float_t metrawphi; // raw metphi from calorimeter Float_t sumet; // total Et in the calorimeter Float_t sume; // total energy in the calorimeter Float_t etoutoftime; // not yet filled Float_t metcorall; // met corrected with all corrections Float_t metcorallphi; // metphi corrected with all corrections // For the individual corrections, the correct way to apply them is to // ADD them to the x & y components of the raw met. // For CMUO Corrections, we check the idword of the muon and see if it // contains 0x8f1 (CMU), 0x8f2 (CMP), 0x8f4 (CMX). If so, we // add back the calorimeter deposits and subtract out the muon momenta. Float_t delcmuo_x; // CMU/P/X muons Float_t delcmuo_y; // For CMIO corrections, we check that idword & 0x8f1 == 0x8f0. // then we do the same thing... Float_t delcmio_x; // central stubless muons Float_t delcmio_y; // These corrections are not made... Float_t delbmu_x; // BMU muons Float_t delbmu_y; // These corrections are for loose high-pt isolated // electrons and photons. A good photon has 5025 < status < 6000. // A good electron has either idword & 0x8ab = 0x8ab or // idword & 0x70000 = 0x70000 and et > 12. // The correction is then -px + jet.px[jetind] where px is the appropriate // electron or photon (corrected) momentum and jetind is the jet matching to // the em object. Float_t delem_x; // EM objects Float_t delem_y; // The following are the met corrections for Central and Plug Jets // with Et > 15 GeV. The correction used is (myjet.corfd[i]-1.0)*myjet.px[i] Float_t delcj_x; // central jets (eta < 1.1) Float_t delcj_y; Float_t delpj_x; // plug jets (1.1 < eta < 2.5) Float_t delpj_y; Int_t emtowoutoftime; Int_t hadtowoutoftime; Float_t emeoutoftime; Float_t hadeoutoftime; }; //#include "UCNtuple/mymet_t.icc" #endif