82 |
|
|
83 |
|
Int_t numjets; |
84 |
|
//Jet Kinematic |
85 |
< |
/*Double_t jetpt[NMAX]; |
86 |
< |
Double_t jeteta[NMAX]; |
85 |
> |
Double_t jet1pt; |
86 |
> |
Double_t jet2pt; |
87 |
> |
Double_t jet3pt; |
88 |
> |
Double_t jet4pt; |
89 |
> |
Double_t deltaRdijet; |
90 |
> |
/*Double_t jeteta[NMAX]; |
91 |
|
Double_t jetphi[NMAX]; |
92 |
< |
Double_t jete[NMAX]; |
92 |
> |
Double_t jete[NMAX];*/ |
93 |
|
|
94 |
|
//mc |
95 |
+ |
/* |
96 |
|
Int_t jetparton[NMAX]; |
97 |
|
|
98 |
|
//b tag |
100 |
|
Double_t jettche[NMAX]; |
101 |
|
Double_t jetsvhe[NMAX]; |
102 |
|
Double_t jetsvhp[NMAX]; |
103 |
< |
Double_t jetcsv[NMAX]; |
104 |
< |
|
103 |
> |
Double_t jetcsv[NMAX]; */ |
104 |
> |
Int_t numBjets_csvl; |
105 |
> |
Int_t numBjets_csvm; |
106 |
> |
Int_t numBjets_csvt; |
107 |
> |
|
108 |
> |
Double_t jet1CSVLpt; |
109 |
> |
Double_t jet2CSVLpt; |
110 |
> |
Double_t jet3CSVLpt; |
111 |
> |
Double_t jet4CSVLpt; |
112 |
> |
Double_t jet1CSVMpt; |
113 |
> |
Double_t jet2CSVMpt; |
114 |
> |
Double_t jet3CSVMpt; |
115 |
> |
Double_t jet4CSVMpt; |
116 |
|
//Jet Track number |
117 |
< |
Int_t jettracksize[NMAX]; |
118 |
< |
*/ |
103 |
< |
|
104 |
< |
// B discriminator |
105 |
< |
Double_t bdisc_1st; |
106 |
< |
Double_t bdisc_2nd; |
107 |
< |
Double_t bdisc_3rd; |
108 |
< |
Double_t bdisc_4th; |
109 |
< |
|
117 |
> |
//Int_t jettracksize[NMAX]; |
118 |
> |
|
119 |
|
}; |
120 |
|
|
121 |
|
|
135 |
|
|
136 |
|
struct MetVariables{ |
137 |
|
|
138 |
< |
Double_t met; |
139 |
< |
Double_t metphi; |
138 |
> |
Double_t MET; |
139 |
> |
/* |
140 |
> |
Double_t metphi; |
141 |
|
Double_t mht; |
142 |
|
Double_t neupt; |
143 |
|
Double_t rho; |
160 |
|
Int_t csvbtagjetnum_up; |
161 |
|
Int_t csvbtagjetnum_down; |
162 |
|
Int_t flavorhistory; |
163 |
+ |
*/ |
164 |
+ |
Double_t Ht; |
165 |
+ |
Double_t Stlep; |
166 |
+ |
Double_t Stjet; |
167 |
+ |
Double_t DeltaPhiWW; |
168 |
+ |
Double_t DeltaRWW; |
169 |
+ |
Double_t DeltaPhiMETWlep; |
170 |
+ |
Double_t DeltaPhiMETlep; |
171 |
+ |
Double_t DeltaPhiNulep; |
172 |
+ |
Double_t LepWPt; |
173 |
|
}; |
174 |
|
|
175 |
|
struct GeneralVariables { |
308 |
|
cout<<"We are initial Jet Branch"<<endl; |
309 |
|
|
310 |
|
Top->Branch("jet_number",&VStoreJet.numjets,"numjets/I"); |
311 |
< |
/* |
312 |
< |
Top->Branch("jet_pt",VStoreJet.jetpt,"jetpt[numjets]/D"); |
313 |
< |
Top->Branch("jet_eta",VStoreJet.jeteta,"jeteta[numjets]/D"); |
311 |
> |
|
312 |
> |
Top->Branch("jet1_pt",&VStoreJet.jet1pt,"jet1pt/D"); |
313 |
> |
Top->Branch("jet2_pt",&VStoreJet.jet2pt,"jet2pt/D"); |
314 |
> |
Top->Branch("jet3_pt",&VStoreJet.jet3pt,"jet3pt/D"); |
315 |
> |
Top->Branch("jet4_pt",&VStoreJet.jet4pt,"jet4pt/D"); |
316 |
> |
Top->Branch("deltaRdijet",&VStoreJet.deltaRdijet,"deltaRdijet/D"); |
317 |
> |
/*Top->Branch("jet_eta",VStoreJet.jeteta,"jeteta[numjets]/D"); |
318 |
|
Top->Branch("jet_phi",VStoreJet.jetphi,"jetphi[numjets]/D"); |
319 |
|
Top->Branch("jet_energy",VStoreJet.jete,"jete[numjets]/D"); |
320 |
|
Top->Branch("jet_partonflavor",VStoreJet.jetparton,"jetparton[numjets]/I"); |
325 |
|
Top->Branch("jet_btag_SVhighputy",VStoreJet.jetsvhp,"jetsvhp[numjets]/D"); |
326 |
|
Top->Branch("jet_btag_csv",VStoreJet.jetcsv,"jetcsv[numjets]/D"); |
327 |
|
*/ |
328 |
< |
//b discriminant |
329 |
< |
Top->Branch("jet_bdisc_1st",&VStoreJet.bdisc_1st,"bdisc_1st/D"); |
330 |
< |
Top->Branch("jet_bdisc_2nd",&VStoreJet.bdisc_2nd,"bdisc_2nd/D"); |
331 |
< |
Top->Branch("jet_bdisc_3rd",&VStoreJet.bdisc_3rd,"bdisc_3rd/D"); |
332 |
< |
Top->Branch("jet_bdisc_4th",&VStoreJet.bdisc_4th,"bdisc_4th/D"); |
333 |
< |
|
328 |
> |
Top->Branch("numBjets_csvl",&VStoreJet.numBjets_csvl,"numBjets_csvl/I"); |
329 |
> |
Top->Branch("numBjets_csvm",&VStoreJet.numBjets_csvm,"numBjets_csvm/I"); |
330 |
> |
Top->Branch("numBjets_csvt",&VStoreJet.numBjets_csvt,"numBjets_csvt/I"); |
331 |
> |
|
332 |
> |
Top->Branch("jet1CSVLpt",&VStoreJet.jet1CSVLpt,"jet1CSVLpt/D"); |
333 |
> |
Top->Branch("jet2CSVLpt",&VStoreJet.jet2CSVLpt,"jet2CSVLpt/D"); |
334 |
> |
Top->Branch("jet3CSVLpt",&VStoreJet.jet3CSVLpt,"jet3CSVLpt/D"); |
335 |
> |
Top->Branch("jet4CSVLpt",&VStoreJet.jet4CSVLpt,"jet4CSVLpt/D"); |
336 |
> |
Top->Branch("jet1CSVMpt",&VStoreJet.jet1CSVMpt,"jet1CSVMpt/D"); |
337 |
> |
Top->Branch("jet2CSVMpt",&VStoreJet.jet2CSVMpt,"jet2CSVMpt/D"); |
338 |
> |
Top->Branch("jet3CSVMpt",&VStoreJet.jet3CSVMpt,"jet3CSVMpt/D"); |
339 |
> |
Top->Branch("jet4CSVMpt",&VStoreJet.jet4CSVMpt,"jet4CSVMpt/D"); |
340 |
|
//Jet track size |
341 |
|
//Top->Branch("jet_tracksize",VStoreJet.jettracksize,"jettracksize[numjets]/I"); |
312 |
– |
|
342 |
|
} |
343 |
|
|
344 |
|
} |
348 |
|
if(StoreMet){ |
349 |
|
cout<<"We are inital Met Branch"<<endl; |
350 |
|
|
351 |
< |
Top->Branch("MET_energy",&VStoreMet.met,"met/D"); |
352 |
< |
Top->Branch("MET_phi",&VStoreMet.metphi,"metphi/D"); |
351 |
> |
Top->Branch("MET",&VStoreMet.MET,"MET/D"); |
352 |
> |
/*Top->Branch("MET_phi",&VStoreMet.metphi,"metphi/D"); |
353 |
|
Top->Branch("MHT_energy",&VStoreMet.mht,"mht/D"); |
354 |
|
Top->Branch("Neutrinuo_pt",&VStoreMet.neupt,"neupt/D"); |
355 |
|
Top->Branch("pileup_rho",&VStoreMet.rho,"rho/D"); |
371 |
|
Top->Branch("csvbtagjetnum_up",&VStoreMet.csvbtagjetnum_up,"csvbtagjetnum_up/I"); |
372 |
|
Top->Branch("csvbtagjetnum_down",&VStoreMet.csvbtagjetnum_down,"csvbtagjetnum_down/I"); |
373 |
|
Top->Branch("flavorhistory",&VStoreMet.flavorhistory,"flavorhistory/I"); |
374 |
< |
|
374 |
> |
*/ |
375 |
> |
Top->Branch("Ht",&VStoreMet.Ht,"Ht/D"); |
376 |
> |
Top->Branch("Stlep",&VStoreMet.Stlep,"Stlep/D"); |
377 |
> |
Top->Branch("Stjet",&VStoreMet.Stjet,"Stjet/D"); |
378 |
> |
Top->Branch("DeltaPhiWW",&VStoreMet.DeltaPhiWW,"DeltaPhiWW/D"); |
379 |
> |
Top->Branch("DeltaRWW",&VStoreMet.DeltaRWW,"DeltaRWW/D"); |
380 |
> |
Top->Branch("DeltaPhiMETWlep",&VStoreMet.DeltaPhiMETWlep,"DeltaPhiMETWlep/D"); |
381 |
> |
Top->Branch("DeltaPhiMETlep",&VStoreMet.DeltaPhiMETlep,"DeltaPhiMETlep/D"); |
382 |
> |
Top->Branch("DeltaPhiNulep",&VStoreMet.DeltaPhiNulep,"DeltaPhiNulep/D"); |
383 |
> |
Top->Branch("LepWPt",&VStoreMet.LepWPt,"LepWPt/D"); |
384 |
|
} |
385 |
|
} |
386 |
|
|