snippets/ElectronID/InclusiveSkim.h

0001 TFile* outFile = NULL;
0002 TTree* outTree = NULL;
0003
0004 const double Me = 0.511e-3; // GeV
0005 const double Mp = 0.938272; // GeV
0006 const double crossing_angle = -0.025; // rad
0007
0008 void SetInputBranchAddresses();
0009 void CreateOutputTree(TString outFileName);
0010 void ResetVariables();
0011 void CalculateElectronKinematics(double fEe, double fEh, TLorentzVector kf, float& xB, float& Q2, float& W, float& y, float& nu);
0012 TLorentzVector GetHadronBeam(double fEh);
0013
0014 using namespace std;
0015
0016 int positive_eID;
0017
0018 float mc_p;
0019 float mc_eta;
0020 float mc_phi;
0021
0022 float track_p;
0023 float track_eta;
0024 float track_phi;
0025
0026 float mc_xB;
0027 float mc_Q2;
0028 float mc_W;
0029 float mc_y;
0030 float mc_nu;
0031
0032 float e_track_xB;
0033 float e_track_Q2;
0034 float e_track_W;
0035 float e_track_y;
0036 float e_track_nu;
0037
0038 float e_clust_xB;
0039 float e_clust_Q2;
0040 float e_clust_W;
0041 float e_clust_y;
0042 float e_clust_nu;
0043