-
Notifications
You must be signed in to change notification settings - Fork 0
/
looper.h
81 lines (65 loc) · 1.77 KB
/
looper.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
#include "EventTree.h"
#include "GenParticleTree.h"
#include "LeptonTree.h"
#include "JetTree.h"
#include "TauTree.h"
#include "IsoTracksTree.h"
#include "TFile.h"
#include "TTree.h"
#include "TRandom.h"
#include "TF1.h"
#include "Math/VectorUtil.h"
#include "TChain.h"
#include "Math/LorentzVector.h"
// // typedefs
typedef ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<float> > LorentzVector;
#ifndef LOOPER_H
#define LOOPER_H
#pragma GCC diagnostic ignored "-Wwrite-strings"
using namespace std;
class babyMaker {
public:
//constructor/destructor
babyMaker ();
babyMaker (const std::string &prefix);
virtual ~babyMaker (){}
void MakeBabyNtuple(const char* output_name);
void InitBabyNtuple();
int looper(TChain* chain, char* output_name, int nEvents = -1, char* path = "./");
char* babypath;
//vaariables for baby skim
int skim_nvtx;
float skim_met;
int skim_nlep;
float skim_leppt;
float skim_lepeta;
int skim_njets;
float skim_jetpt;
float skim_jeteta;
void setSkimVariables(int nvtx, float met, int nlep, float leppt, float lepeta, int njets, int jetpt, int jeteta);
protected:
TFile* BabyFile;
TTree* BabyTree;
private:
//for tree
EventTree StopEvt;
LeptonTree lep1;
LeptonTree lep2;
JetTree jets;
GenParticleTree gen_els;
GenParticleTree gen_mus;
GenParticleTree gen_leptau_els;
GenParticleTree gen_leptau_mus;
GenParticleTree gen_taus;
GenParticleTree gen_nus;
GenParticleTree gen_bs;
GenParticleTree gen_tops;
GenParticleTree gen_qs;
GenParticleTree gen_lsp;
GenParticleTree gen_stop;
TauTree Taus;
IsoTracksTree Tracks;
// = new EventTree(); //dont make a pointer
};
struct val_err_t { float value; float error; };
#endif