Ticket #997: mononew_zero_jet_left_subjet_1.C

#define mononew_zero_jet_left_subjet_1_cxx
#include "mononew_zero_jet_left_subjet_1.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <TLorentzVector.h>

void mononew_zero_jet_left_subjet_1::Loop()
{
//   In a ROOT session, you can do:
//      Root > .L mononew_zero_jet_left_subjet_1.C
//      Root > mononew_zero_jet_left_subjet_1 t
//      Root > t.GetEntry(12); // Fill t data members with entry number 12
//      Root > t.Show();       // Show values of entry 12
//      Root > t.Show(16);     // Read and show values of entry 16
//      Root > t.Loop();       // Loop on all entries
//

//     This is the loop skeleton where:
//    jentry is the global entry number in the chain
//    ientry is the entry number in the current Tree
//  Note that the argument to GetEntry must be:
//    jentry for TChain::GetEntry
//    ientry for TTree::GetEntry and TBranch::GetEntry
//
//       To read only selected branches, Insert statements like:
// METHOD1:
//    fChain->SetBranchStatus("*",0);  // disable all branches
//    fChain->SetBranchStatus("branchname",1);  // activate branchname
// METHOD2: replace line
//    fChain->GetEntry(jentry);       //read all branches
//by  b_branchname->GetEntry(ientry); //read only this branch
   if (fChain == 0) return;

   Long64_t nentries = fChain->GetEntriesFast();

   Long64_t nbytes = 0, nb = 0;
   for (Long64_t jentry=0; jentry<nentries;jentry++) {
      Long64_t ientry = LoadTree(jentry);
      if (ientry < 0) break;
      nb = fChain->GetEntry(jentry);   nbytes += nb;
      // if (Cut(ientry) < 0) continue;
   }
}

void mononew_zero_jet_left_subjet_1::Preprocess()
{
  int evEntry = 0;

  ofstream out_file("mononew_zero_jet_left_subjet_1.ppdat");

  // while (GetEntry(evEntry) != 0)
  // GetEntry(evEntry);
  while (evEntry < 100)
  {
    GetEntry(evEntry);
 
    for (int i = 0; i < Jet_; i++)
    {
      out_file << Jet_NSubJetsTrimmed[i] << endl;

      // This is just an example: 
      out_file << (Jet_TrimmedP4[0][i])[0] << " " 
               << (Jet_TrimmedP4[1][i])[0] << " " 
               << (Jet_TrimmedP4[2][i])[0] << " " 
               << (Jet_TrimmedP4[3][i])[0] << " " 
               << (Jet_TrimmedP4[4][i])[0] << endl;
    }

    out_file << Jet_ << endl;

    evEntry++;

  } // end while on events

  out_file.close();

}