double countTotal(process &process) { // Create temporary histogram TH1D* tmp = new TH1D("tmp", "tmp", 30, 0, 1000); double Total = 0; process.GetChain()->Draw("Njet>>tmp", "XS*5000/Nevt", "goff"); Total = int(tmp->Integral()); delete tmp; process.SetNevt_nosel(Total); return Total; }
double countEvt(process &process, TCut Cut) { // Create temporary histogram TH1D* tmp = new TH1D("tmp", "tmp", 30, 0, 1000); // Get number of events after a cut double Sel = 0; process.GetChain()->Draw("Njet>>tmp","XS*5000/Nevt"*Cut, "goff"); Sel = tmp->Integral(); delete tmp; process.SetSelection(Cut.GetName()); process.SetNevt_sel(Sel); return Sel; }