-
Notifications
You must be signed in to change notification settings - Fork 37
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
latest changes on gep trigger simulations, sanity-checking PYTHIA6 cr…
…oss sections, etc
- Loading branch information
Showing
10 changed files
with
3,751 additions
and
662 deletions.
There are no files selected for viewing
Large diffs are not rendered by default.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,43 @@ | ||
#define gep_tree_pythia_cxx | ||
#include "gep_tree_pythia.h" | ||
#include <TH2.h> | ||
#include <TStyle.h> | ||
#include <TCanvas.h> | ||
|
||
void gep_tree_pythia::Loop() | ||
{ | ||
// In a ROOT session, you can do: | ||
// root> .L gep_tree_pythia.C | ||
// root> gep_tree_pythia 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; | ||
} | ||
} |
Oops, something went wrong.