float bdtjetvar[81]; //variables for boosted decision tree bdtjetvar[i] = 0.; bdtjetvar[1] = (*PtGen)[i]; // Pt of anti-top (MC) bdtjetvar[2] = (*EtaGen)[i]; // eta of anti-top (MC) bdtjetvar[3] = (*PhiGen)[i]; // phi of anti-top (MC) bdtjetvar[4] = (*MGen)[i]; // mass of anti-top (MC) bdtjetvar[5] = (*RvGen)[i]; // Rvertex of anti-top (MC) bdtjetvar[6] = (*PhiVGen)[i];// Phi vertex of anti-top (MC) bdtjetvar[7] = (*ZVGen)[i]; // Z vertex of anti-top (MC) bdtjetvar[8] = (*PtGen)[i]; // Pt of W- (MC) bdtjetvar[9] = (*EtaGen)[i]; // eta of W- (MC) bdtjetvar[10] = (*PhiGen)[i]; // phi of W- (MC) bdtjetvar[11] = (*MGen)[i]; // mass of W- (MC) bdtjetvar[12] = (*RvGen)[i]; // Rvertex of W- (MC) bdtjetvar[13] = (*PhiVGen)[i];// Phi vertex of W- (MC) bdtjetvar[14] = (*ZVGen)[i]; // Z vertex of W- (MC) bdtjetvar[15] = dWpt; // Pt of W bdtjetvar[16] = dWeta; // eta of W bdtjetvar[17] = dWphi; // phi of W bdtjetvar[18] = dWMass; // mass of W bdtjetvar[19] = dWeng; // Energy of W bdtjetvar[20] = dphi; // delta phi of W->2jets bdtjetvar[21] = dWangle; // solid angle of W->2jets bdtjetvar[22] = dTpt; // Pt of top bdtjetvar[23] = dTeta; // eta of top bdtjetvar[24] = dTphi; // phi of top bdtjetvar[25] = dTMass; // rec top mass bdtjetvar[26] = Njet20; // number of jets with Et>20 GeV bdtjetvar[27] = Njet30; // number of jets with Et>30 GeV bdtjetvar[28] = Njet40; // number of jets with Et>40 GeV bdtjetvar[29] = Njet120; // number of jets with Et>120 GeV bdtjetvar[30] = MET_EtMiss; // Missing Et bdtjetvar[31] = sqrt(yjet_px[njet1]**2 + yjet_py[njet1]**2); // pt of W jet1 bdtjetvar[32] = yjet_eta[njet1]; // eta of W jet1 bdtjetvar[33] = yjet_phi[njet1]; // phi of W jet1 bdtjetvar[34] = yjet_m[njet1]; // mass of W jet1 bdtjetvar[35] = yjet_e[njet1]; // energy of W jet1 bdtjetvar[36] = yjet_size[njet1];// size of W jet1 bdtjetvar[37] = yjet_eem[njet1]; // EM energy of W jet1 bdtjetvar[38] = sqrt(yjet_px[njet2]**2 + yjet_py[njet2]**2); // pt of W jet2 bdtjetvar[39] = yjet_eta[njet2]; // eta of W jet2 bdtjetvar[40] = yjet_phi[njet2]; // phi of W jet2 bdtjetvar[41] = yjet_m[njet2]; // mass of W jet2 bdtjetvar[42] = yjet_e[njet2]; // energy of W jet2 bdtjetvar[43] = yjet_size[njet2];// size of W jet2 bdtjetvar[44] = yjet_eem[njet2]; // EM energy of W jet2 bdtjetvar[45] = sqrt(yjet_px[njet3]**2 + yjet_py[njet3]**2); // pt of b jet bdtjetvar[46] = yjet_eta[njet3]; // eta of b jet bdtjetvar[47] = yjet_phi[njet3]; // phi of b jet bdtjetvar[48] = yjet_m[njet3]; // mass of b jet bdtjetvar[49] = yjet_e[njet3]; // energy of b jet bdtjetvar[50] = yjet_size[njet3];// size of b jet bdtjetvar[51] = yjet_eem[njet3]; // EM energy of b jet bdtjetvar[52] = yjet_et1; // Et of the 1st energetic jet bdtjetvar[53] = yjet_et2; // Et of the 2nd energetic jet bdtjetvar[54] = yjet_et3; // Et of the 3rd energetic jet bdtjetvar[55] = yjet_et4; // Et of the 4th energetic jet bdtjetvar[56] = yjet_et_sum; // Et of all jets bdtjetvar[57] = dphiE1E2; // dphi between jet1 and jet2 bdtjetvar[58] = dphiE1E3; // dphi between jet1 and jet3 bdtjetvar[59] = dphiE2E3; // dphi between jet2 and jet3 bdtjetvar[60] = drE1E2; // dR between jet1 and jet2 bdtjetvar[61] = drE1E3; // dR between jet1 and jet3 bdtjetvar[62] = drE2E3; // dR between jet2 and jet3 bdtjetvar[63] = iflag_w_lep; // number of leptons from W decay bdtjetvar[64] = YJet_sum.M(); // mass of all jets bdtjetvar[65] = nblep; // number of selected lepton bdtjetvar[66] = lepton_id; // lepton id number 1-e-, 2-e+, 3-m-, 4-m+ if(index_muon>=0) bdtjetvar[67] = muflag[index_muon]; //muon flag, >0 true, <0 fake if(index_electron>=0) bdtjetvar[68] = eflag[index_electron]; //e flag, >0 true, <0 fake bdtjetvar[69] = Run; // Run number bdtjetvar[70] = Event; // Event number bdtjetvar[71] = IEvent; // Event squential number bdtjetvar[72] = Time; // Time stamp fprintf(foutjet, "%14.6e", bdtjetvar[iv]);