diff --git a/examples/gtrk_e/comp_gtrk_e.C b/examples/gtrk_e/comp_gtrk_e.C deleted file mode 100644 index a3e1a0675f9e5f2f9a135a3c6f56c2f3684b3447..0000000000000000000000000000000000000000 --- a/examples/gtrk_e/comp_gtrk_e.C +++ /dev/null @@ -1,200 +0,0 @@ -#define hcana_class_cxx -#define engine_class_cxx - -#include <TTree.h> -#include <TFile.h> -#include <TROOT.h> -#include <iostream> -#include <TMath.h> -#include <TLegend.h> -#include "hcana_class.h" -#include "engine_class.h" - -using namespace std; - -void comp_gtrk_e() { - - // Compare Ecalo of the golden track from Engine and hcana. - - // hcana tree - - hcana_class* hcana_tree = new hcana_class(0); - - Long64_t hcana_nentries = hcana_tree->fChain->GetEntriesFast(); - cout << "hcana_nentries = " << hcana_nentries << endl; - - // hcana_tree->fChain->SetBranchStatus("*",0); // disable all branches - // hcana_tree->fChain->SetBranchStatus("*H.gold.e",1); // activate branchname - // hcana_tree->fChain->SetBranchStatus("*H.gold.p",1); - // hcana_tree->fChain->SetBranchStatus("*H.gold.index",1); - // hcana_tree->fChain->SetBranchStatus("*H.tr.n",1); - // hcana_tree->fChain->SetBranchStatus("*g.evtyp",1); - // hcana_tree->fChain->SetBranchStatus("*fEvtHdr.fEvtType",1); - // hcana_tree->fChain->SetBranchStatus("*fEvtHdr.fEvtNum",1); - - // Engine tree - - engine_class* engine_tree = new engine_class(0); - - Long64_t engine_nentries = engine_tree->fChain->GetEntriesFast(); - cout << "engine_nentries = " << engine_nentries << endl; - - // engine_tree->fChain->SetBranchStatus("*",0); - // engine_tree->fChain->SetBranchStatus("*hstrk_et",1); - // engine_tree->fChain->SetBranchStatus("*hsp",1); - // engine_tree->fChain->SetBranchStatus("*ev_type",1); - // engine_tree->fChain->SetBranchStatus("*eventid",1); - // engine_tree->fChain->SetBranchStatus("*hstrk_in",1); - // engine_tree->fChain->SetBranchStatus("*hntracks",1); - - // Histograms. - - TH1F* he_hcana = new TH1F("e_hcana", "Golden track Ecalo", - 100, 0.011, 1.8); - he_hcana->SetLineColor(kRed); - - TH1F* he_engine = (TH1F*) he_hcana->Clone("e_engine"); - he_engine->SetLineColor(kBlue); - - TH1F* ee_diff = new TH1F("ee_diff", - "Event by event hcana - Engine Ecalo difference", 200,-1.,1.); - TH1F* pp_diff = new TH1F("pp_diff","Event by event P difference", - 200,-1.,1.); - TH2F* pvse_diff = new TH2F("pvse_diff","Event by event P vs Ecalo difference", - 200,-1.,1., 200,-100.,100.); - pvse_diff->SetMarkerStyle(3); - TH2F* evse = new TH2F("evse","Event by event hcana vs Engine Ecalo", - 210,-0.1,2., 210,-0.1,2.); - evse->SetMarkerStyle(3); - - // Loop over entries. - - Long64_t nentries = TMath::Min(hcana_nentries, engine_nentries); - cout << "nentries = " << nentries << endl; - - Long64_t jentry_hcana = 0; - Long64_t jentry_engine = 0; - for (Long64_t jentry=0; jentry<nentries; jentry++) { - - hcana_tree->fChain->GetEntry(jentry_hcana++); - engine_tree->fChain->GetEntry(jentry_engine++); - - while (hcana_tree->fEvtHdr_fEvtNum > engine_tree->eventid && - jentry_engine < engine_nentries) { - engine_tree->fChain->GetEntry(jentry_engine++); - //cout<< hcana_tree->fEvtHdr_fEvtNum << " " << engine_tree->eventid << endl; - } - - while (hcana_tree->fEvtHdr_fEvtNum < engine_tree->eventid && - jentry_hcana < hcana_nentries) { - hcana_tree->fChain->GetEntry(jentry_hcana++); - //cout<< hcana_tree->fEvtHdr_fEvtNum << " " << engine_tree->eventid << endl; - } - - if (hcana_tree->fEvtHdr_fEvtType==1 && engine_tree->ev_type==1 - && hcana_tree->H_gold_index>-1 && engine_tree->hstrk_in>0 - && hcana_tree->H_gold_index == engine_tree->hstrk_in - 1) { - - he_hcana->Fill(hcana_tree->H_gold_e, 1.); - //cout << hcana_tree->H_gold_e << " " << hcana_tree->g_evtyp << " " - // << hcana_tree->fEvtHdr_fEvtType << endl; - - he_engine->Fill(engine_tree->hstrk_et, 1.); - // cout << engine_tree->hsshtrk << " " << engine_tree->ev_type << " " - // << engine_tree->eventid << endl; - - ee_diff->Fill(hcana_tree->H_gold_e - engine_tree->hstrk_et, 1.); - pp_diff->Fill(hcana_tree->H_gold_p - engine_tree->hsp, 1.); - pvse_diff->Fill(hcana_tree->H_gold_e - engine_tree->hstrk_et, - hcana_tree->H_gold_p - engine_tree->hsp, 1.); - evse->Fill(engine_tree->hstrk_et, hcana_tree->H_gold_e, 1.); - - //Debug print out. - - // if (hcana_tree->H_gold_e>0.6 && hcana_tree->H_gold_e<0.95 && - // engine_tree->hstrk_et<0.1) { - // if (hcana_tree->H_gold_e>0.6 && hcana_tree->H_gold_e<0.95 && - // TMath::Abs(engine_tree->hstrk_et - hcana_tree->H_gold_e) < 0.0001 ) { - // if (engine_tree->hstrk_et>0.6 && engine_tree->hstrk_et<0.95 && - // hcana_tree->H_gold_e<0.1) { - if (TMath::Abs(engine_tree->hstrk_et - hcana_tree->H_gold_e) < 0.65 && - TMath::Abs(engine_tree->hstrk_et - hcana_tree->H_gold_e) > 0.05 ) { - - cout << "===========================================================\n"; - cout << "Engine event " << engine_tree->eventid - << ": hstrk_et = " << engine_tree->hstrk_et - << ", hstrk_index = " << engine_tree->hstrk_in - << ", hntracks = " << engine_tree->hntracks << endl; - cout << " hsxfp = " << engine_tree->hsxfp - << ", hsyfp = " << engine_tree->hsyfp << endl; - cout << " hsxpfp= " << engine_tree->hsxpfp - << ", hsypfp= " << engine_tree->hsypfp << endl; - cout << " hsp = " << engine_tree->hsp << endl; - cout << " At calorimeter:" - << " x = " << engine_tree->hsxfp+338.69*engine_tree->hsxpfp - << " y = " << engine_tree->hsyfp+338.69*engine_tree->hsypfp - << endl; - cout << "-----------------------------------------------------------\n"; - cout << "hcana event " << hcana_tree->fEvtHdr_fEvtNum - << ": H_gold_e = " << hcana_tree->H_gold_e - << " H_gold_index = " << hcana_tree->H_gold_index - << " H_tr_n = " << hcana_tree->H_tr_n << endl; - cout << " H_gold_fp_x = " << hcana_tree->H_gold_fp_x - << " H_gold_fp_y = " << hcana_tree->H_gold_fp_y << endl; - cout << " H_gold_fp_th= " << hcana_tree->H_gold_fp_th - << " H_gold_fp_ph= " << hcana_tree->H_gold_fp_ph << endl; - cout << " H_gold_p = " << hcana_tree->H_gold_p << endl; - cout << " At calorimeter:" - << " x = " << hcana_tree->H_gold_fp_x+338.69*hcana_tree->H_gold_fp_th - << " y = " << hcana_tree->H_gold_fp_y+338.69*hcana_tree->H_gold_fp_ph - << endl; - - } - - } - - } //over entries - - // Plot. - - TCanvas* c1 = new TCanvas("c1","Golden Track E from hcana",600,750); - - c1->Divide(1,2); - c1->cd(1); - he_hcana->Draw(); - he_engine->Draw("same"); - - TLegend* leg = new TLegend(0.25, 0.70, 0.40, 0.82); - leg->AddEntry(he_hcana,"hcana","l"); - leg->AddEntry(he_engine,"Engine","l"); - leg->Draw(); - - // Difference between the histograms. - - c1->cd(2); - - TH1F* dif = he_hcana->Clone("difference"); - dif->Reset(); - dif->Add(he_hcana,he_engine,1.,-1.); - - dif->SetTitle("hcana -- Engine difference"); - dif->SetFillColor(kRed); - dif->SetLineColor(kRed); - dif->SetLineWidth(1); - dif->SetFillStyle(1111); - dif->Draw(); - - TCanvas* c2 = new TCanvas("c2","Event by event hcana - Engine Ecalo diff.", - 600,750); - - c2->Divide(1,2); - c2->cd(1); - gPad->SetLogy(); - ee_diff->Draw(); - c2->cd(2); - // gPad->SetLogy(); - // pp_diff->Draw(); - // pvse_diff->Draw(); - evse->Draw(); - -} diff --git a/examples/gtrk_e/engine_class.C b/examples/gtrk_e/engine_class.C deleted file mode 100644 index b437035942eff550c1172237ffbe2728edd6c0b8..0000000000000000000000000000000000000000 --- a/examples/gtrk_e/engine_class.C +++ /dev/null @@ -1,43 +0,0 @@ -#define engine_class_cxx -#include "engine_class.h" -#include <TH2.h> -#include <TStyle.h> -#include <TCanvas.h> - -void engine_class::Loop() -{ -// In a ROOT session, you can do: -// Root > .L engine_class.C -// Root > engine_class 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; - } -} diff --git a/examples/gtrk_e/engine_class.h b/examples/gtrk_e/engine_class.h deleted file mode 100644 index 3a6d2e5f8321532e4682d742e46fc6aa54524a53..0000000000000000000000000000000000000000 --- a/examples/gtrk_e/engine_class.h +++ /dev/null @@ -1,255 +0,0 @@ -////////////////////////////////////////////////////////// -// This class has been automatically generated on -// Mon Sep 29 17:54:32 2014 by ROOT version 5.34/06 -// from TTree h9010/ -// found on file: hms52949_rz.root -////////////////////////////////////////////////////////// - -#ifndef engine_class_h -#define engine_class_h - -#include <TROOT.h> -#include <TChain.h> -#include <TFile.h> - -// Header file for the classes stored in the TTree if any. - -// Fixed size dimensions of array or collections stored in the TTree if any. - -class engine_class { -public : - TTree *fChain; //!pointer to the analyzed TTree or TChain - Int_t fCurrent; //!current Tree number in a TChain - - // Declaration of leaf types - Float_t hcer_npe; - Float_t hsp; - Float_t hse; - Float_t charge; - Float_t hsdelta; - Float_t hstheta; - Float_t hsphi; - Float_t w; - Float_t hszbeam; - Float_t hsdedx1; - Float_t hsbeta; - Float_t hsshtrk; - Float_t hsprtrk; - Float_t hsxfp; - Float_t hsyfp; - Float_t hsxpfp; - Float_t hsypfp; - Float_t hsytar; - Float_t hsxptar; - Float_t hsyptar; - Float_t hstart; - Float_t eventid; - Float_t ev_type; - Float_t gfrx_raw; - Float_t gfry_raw; - Float_t gbeam_x; - Float_t gbeam_y; - Float_t bpma_x; - Float_t bpma_y; - Float_t bpmb_x; - Float_t bpmb_y; - Float_t bpmc_x; - Float_t bpmc_y; - Float_t hseloss; - Float_t hntracks; - Float_t hcal_et; - Float_t hgoodsc; - Float_t hcal_e1; - Float_t hcal_e2; - Float_t hcal_e3; - Float_t hcal_e4; - Float_t hstrk_et; - Float_t hstrk_in; - - // List of branches - TBranch *b_hcer_npe; //! - TBranch *b_hsp; //! - TBranch *b_hse; //! - TBranch *b_charge; //! - TBranch *b_hsdelta; //! - TBranch *b_hstheta; //! - TBranch *b_hsphi; //! - TBranch *b_w; //! - TBranch *b_hszbeam; //! - TBranch *b_hsdedx1; //! - TBranch *b_hsbeta; //! - TBranch *b_hsshtrk; //! - TBranch *b_hsprtrk; //! - TBranch *b_hsxfp; //! - TBranch *b_hsyfp; //! - TBranch *b_hsxpfp; //! - TBranch *b_hsypfp; //! - TBranch *b_hsytar; //! - TBranch *b_hsxptar; //! - TBranch *b_hsyptar; //! - TBranch *b_hstart; //! - TBranch *b_eventid; //! - TBranch *b_ev_type; //! - TBranch *b_gfrx_raw; //! - TBranch *b_gfry_raw; //! - TBranch *b_gbeam_x; //! - TBranch *b_gbeam_y; //! - TBranch *b_bpma_x; //! - TBranch *b_bpma_y; //! - TBranch *b_bpmb_x; //! - TBranch *b_bpmb_y; //! - TBranch *b_bpmc_x; //! - TBranch *b_bpmc_y; //! - TBranch *b_hseloss; //! - TBranch *b_hntracks; //! - TBranch *b_hcal_et; //! - TBranch *b_hgoodsc; //! - TBranch *b_hcal_e1; //! - TBranch *b_hcal_e2; //! - TBranch *b_hcal_e3; //! - TBranch *b_hcal_e4; //! - TBranch *b_hstrk_et; //! - TBranch *b_hstrk_in; //! - - engine_class(TTree *tree=0); - virtual ~engine_class(); - virtual Int_t Cut(Long64_t entry); - virtual Int_t GetEntry(Long64_t entry); - virtual Long64_t LoadTree(Long64_t entry); - virtual void Init(TTree *tree); - virtual void Loop(); - virtual Bool_t Notify(); - virtual void Show(Long64_t entry = -1); -}; - -#endif - -#ifdef engine_class_cxx -engine_class::engine_class(TTree *tree) : fChain(0) -{ -// if parameter tree is not specified (or zero), connect the file -// used to generate this class and read the Tree. - if (tree == 0) { - TFile *f = (TFile*)gROOT->GetListOfFiles()->FindObject("hms52949_rz.root"); - if (!f || !f->IsOpen()) { - f = new TFile("hms52949_rz.root"); - } - f->GetObject("h9010",tree); - - } - Init(tree); -} - -engine_class::~engine_class() -{ - if (!fChain) return; - delete fChain->GetCurrentFile(); -} - -Int_t engine_class::GetEntry(Long64_t entry) -{ -// Read contents of entry. - if (!fChain) return 0; - return fChain->GetEntry(entry); -} -Long64_t engine_class::LoadTree(Long64_t entry) -{ -// Set the environment to read one entry - if (!fChain) return -5; - Long64_t centry = fChain->LoadTree(entry); - if (centry < 0) return centry; - if (fChain->GetTreeNumber() != fCurrent) { - fCurrent = fChain->GetTreeNumber(); - Notify(); - } - return centry; -} - -void engine_class::Init(TTree *tree) -{ - // The Init() function is called when the selector needs to initialize - // a new tree or chain. Typically here the branch addresses and branch - // pointers of the tree will be set. - // It is normally not necessary to make changes to the generated - // code, but the routine can be extended by the user if needed. - // Init() will be called many times when running on PROOF - // (once per file to be processed). - - // Set branch addresses and branch pointers - if (!tree) return; - fChain = tree; - fCurrent = -1; - fChain->SetMakeClass(1); - - fChain->SetBranchAddress("hcer_npe", &hcer_npe, &b_hcer_npe); - fChain->SetBranchAddress("hsp", &hsp, &b_hsp); - fChain->SetBranchAddress("hse", &hse, &b_hse); - fChain->SetBranchAddress("charge", &charge, &b_charge); - fChain->SetBranchAddress("hsdelta", &hsdelta, &b_hsdelta); - fChain->SetBranchAddress("hstheta", &hstheta, &b_hstheta); - fChain->SetBranchAddress("hsphi", &hsphi, &b_hsphi); - fChain->SetBranchAddress("w", &w, &b_w); - fChain->SetBranchAddress("hszbeam", &hszbeam, &b_hszbeam); - fChain->SetBranchAddress("hsdedx1", &hsdedx1, &b_hsdedx1); - fChain->SetBranchAddress("hsbeta", &hsbeta, &b_hsbeta); - fChain->SetBranchAddress("hsshtrk", &hsshtrk, &b_hsshtrk); - fChain->SetBranchAddress("hsprtrk", &hsprtrk, &b_hsprtrk); - fChain->SetBranchAddress("hsxfp", &hsxfp, &b_hsxfp); - fChain->SetBranchAddress("hsyfp", &hsyfp, &b_hsyfp); - fChain->SetBranchAddress("hsxpfp", &hsxpfp, &b_hsxpfp); - fChain->SetBranchAddress("hsypfp", &hsypfp, &b_hsypfp); - fChain->SetBranchAddress("hsytar", &hsytar, &b_hsytar); - fChain->SetBranchAddress("hsxptar", &hsxptar, &b_hsxptar); - fChain->SetBranchAddress("hsyptar", &hsyptar, &b_hsyptar); - fChain->SetBranchAddress("hstart", &hstart, &b_hstart); - fChain->SetBranchAddress("eventid", &eventid, &b_eventid); - fChain->SetBranchAddress("ev_type", &ev_type, &b_ev_type); - fChain->SetBranchAddress("gfrx_raw", &gfrx_raw, &b_gfrx_raw); - fChain->SetBranchAddress("gfry_raw", &gfry_raw, &b_gfry_raw); - fChain->SetBranchAddress("gbeam_x", &gbeam_x, &b_gbeam_x); - fChain->SetBranchAddress("gbeam_y", &gbeam_y, &b_gbeam_y); - fChain->SetBranchAddress("bpma_x", &bpma_x, &b_bpma_x); - fChain->SetBranchAddress("bpma_y", &bpma_y, &b_bpma_y); - fChain->SetBranchAddress("bpmb_x", &bpmb_x, &b_bpmb_x); - fChain->SetBranchAddress("bpmb_y", &bpmb_y, &b_bpmb_y); - fChain->SetBranchAddress("bpmc_x", &bpmc_x, &b_bpmc_x); - fChain->SetBranchAddress("bpmc_y", &bpmc_y, &b_bpmc_y); - fChain->SetBranchAddress("hseloss", &hseloss, &b_hseloss); - fChain->SetBranchAddress("hntracks", &hntracks, &b_hntracks); - fChain->SetBranchAddress("hcal_et", &hcal_et, &b_hcal_et); - fChain->SetBranchAddress("hgoodsc", &hgoodsc, &b_hgoodsc); - fChain->SetBranchAddress("hcal_e1", &hcal_e1, &b_hcal_e1); - fChain->SetBranchAddress("hcal_e2", &hcal_e2, &b_hcal_e2); - fChain->SetBranchAddress("hcal_e3", &hcal_e3, &b_hcal_e3); - fChain->SetBranchAddress("hcal_e4", &hcal_e4, &b_hcal_e4); - fChain->SetBranchAddress("hstrk_et", &hstrk_et, &b_hstrk_et); - fChain->SetBranchAddress("hstrk_in", &hstrk_in, &b_hstrk_in); - Notify(); -} - -Bool_t engine_class::Notify() -{ - // The Notify() function is called when a new file is opened. This - // can be either for a new TTree in a TChain or when when a new TTree - // is started when using PROOF. It is normally not necessary to make changes - // to the generated code, but the routine can be extended by the - // user if needed. The return value is currently not used. - - return kTRUE; -} - -void engine_class::Show(Long64_t entry) -{ -// Print contents of entry. -// If entry is not specified, print current entry - if (!fChain) return; - fChain->Show(entry); -} -Int_t engine_class::Cut(Long64_t entry) -{ -// This function may be called from Loop. -// returns 1 if entry is accepted. -// returns -1 otherwise. - return 1; -} -#endif // #ifdef engine_class_cxx diff --git a/examples/gtrk_e/hcana_class.C b/examples/gtrk_e/hcana_class.C deleted file mode 100644 index 61a1bfb8a575837022217429ec3287ed1b6241a9..0000000000000000000000000000000000000000 --- a/examples/gtrk_e/hcana_class.C +++ /dev/null @@ -1,43 +0,0 @@ -#define hcana_class_cxx -#include "hcana_class.h" -#include <TH2.h> -#include <TStyle.h> -#include <TCanvas.h> - -void hcana_class::Loop() -{ -// In a ROOT session, you can do: -// Root > .L hcana_class.C -// Root > hcana_class 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; - } -} diff --git a/examples/gtrk_e/hcana_class.h b/examples/gtrk_e/hcana_class.h deleted file mode 100644 index bd15df66bbf3678993020df98850959c9f742841..0000000000000000000000000000000000000000 --- a/examples/gtrk_e/hcana_class.h +++ /dev/null @@ -1,649 +0,0 @@ -////////////////////////////////////////////////////////// -// This class has been automatically generated on -// Wed Oct 22 09:01:11 2014 by ROOT version 5.34/06 -// from TTree T/Hall A Analyzer Output DST -// found on file: hodtest_52949.root -////////////////////////////////////////////////////////// - -#ifndef hcana_class_h -#define hcana_class_h - -#include <TROOT.h> -#include <TChain.h> -#include <TFile.h> - -// Header file for the classes stored in the TTree if any. - -// Fixed size dimensions of array or collections stored in the TTree if any. - -class hcana_class { -public : - TTree *fChain; //!pointer to the analyzed TTree or TChain - Int_t fCurrent; //!current Tree number in a TChain - - // Declaration of leaf types - Int_t Ndata_H_cal_1pr_aneg; - Double_t H_cal_1pr_aneg[13]; //[Ndata.H.cal.1pr.aneg] - Int_t Ndata_H_cal_1pr_aneg_p; - Double_t H_cal_1pr_aneg_p[13]; //[Ndata.H.cal.1pr.aneg_p] - Int_t Ndata_H_cal_1pr_apos; - Double_t H_cal_1pr_apos[13]; //[Ndata.H.cal.1pr.apos] - Int_t Ndata_H_cal_1pr_apos_p; - Double_t H_cal_1pr_apos_p[13]; //[Ndata.H.cal.1pr.apos_p] - Int_t Ndata_H_cal_1pr_emean; - Double_t H_cal_1pr_emean[13]; //[Ndata.H.cal.1pr.emean] - Int_t Ndata_H_cal_1pr_eneg; - Double_t H_cal_1pr_eneg[13]; //[Ndata.H.cal.1pr.eneg] - Int_t Ndata_H_cal_1pr_epos; - Double_t H_cal_1pr_epos[13]; //[Ndata.H.cal.1pr.epos] - Int_t Ndata_H_cal_1pr_negadchits; - Double_t H_cal_1pr_negadchits[10]; //[Ndata.H.cal.1pr.negadchits] - Int_t Ndata_H_cal_1pr_posadchits; - Double_t H_cal_1pr_posadchits[8]; //[Ndata.H.cal.1pr.posadchits] - Int_t Ndata_H_cal_2ta_aneg; - Double_t H_cal_2ta_aneg[13]; //[Ndata.H.cal.2ta.aneg] - Int_t Ndata_H_cal_2ta_aneg_p; - Double_t H_cal_2ta_aneg_p[13]; //[Ndata.H.cal.2ta.aneg_p] - Int_t Ndata_H_cal_2ta_apos; - Double_t H_cal_2ta_apos[13]; //[Ndata.H.cal.2ta.apos] - Int_t Ndata_H_cal_2ta_apos_p; - Double_t H_cal_2ta_apos_p[13]; //[Ndata.H.cal.2ta.apos_p] - Int_t Ndata_H_cal_2ta_emean; - Double_t H_cal_2ta_emean[13]; //[Ndata.H.cal.2ta.emean] - Int_t Ndata_H_cal_2ta_eneg; - Double_t H_cal_2ta_eneg[13]; //[Ndata.H.cal.2ta.eneg] - Int_t Ndata_H_cal_2ta_epos; - Double_t H_cal_2ta_epos[13]; //[Ndata.H.cal.2ta.epos] - Int_t Ndata_H_cal_2ta_negadchits; - Double_t H_cal_2ta_negadchits[9]; //[Ndata.H.cal.2ta.negadchits] - Int_t Ndata_H_cal_2ta_posadchits; - Double_t H_cal_2ta_posadchits[10]; //[Ndata.H.cal.2ta.posadchits] - Int_t Ndata_H_cal_3ta_aneg; - Double_t H_cal_3ta_aneg[13]; //[Ndata.H.cal.3ta.aneg] - Int_t Ndata_H_cal_3ta_aneg_p; - Double_t H_cal_3ta_aneg_p[13]; //[Ndata.H.cal.3ta.aneg_p] - Int_t Ndata_H_cal_3ta_apos; - Double_t H_cal_3ta_apos[13]; //[Ndata.H.cal.3ta.apos] - Int_t Ndata_H_cal_3ta_apos_p; - Double_t H_cal_3ta_apos_p[13]; //[Ndata.H.cal.3ta.apos_p] - Int_t Ndata_H_cal_3ta_emean; - Double_t H_cal_3ta_emean[13]; //[Ndata.H.cal.3ta.emean] - Int_t Ndata_H_cal_3ta_eneg; - Double_t H_cal_3ta_eneg[13]; //[Ndata.H.cal.3ta.eneg] - Int_t Ndata_H_cal_3ta_epos; - Double_t H_cal_3ta_epos[13]; //[Ndata.H.cal.3ta.epos] - Int_t Ndata_H_cal_3ta_negadchits; - Double_t H_cal_3ta_negadchits[1]; //[Ndata.H.cal.3ta.negadchits] - Int_t Ndata_H_cal_3ta_posadchits; - Double_t H_cal_3ta_posadchits[9]; //[Ndata.H.cal.3ta.posadchits] - Int_t Ndata_H_cal_4ta_aneg; - Double_t H_cal_4ta_aneg[13]; //[Ndata.H.cal.4ta.aneg] - Int_t Ndata_H_cal_4ta_aneg_p; - Double_t H_cal_4ta_aneg_p[13]; //[Ndata.H.cal.4ta.aneg_p] - Int_t Ndata_H_cal_4ta_apos; - Double_t H_cal_4ta_apos[13]; //[Ndata.H.cal.4ta.apos] - Int_t Ndata_H_cal_4ta_apos_p; - Double_t H_cal_4ta_apos_p[13]; //[Ndata.H.cal.4ta.apos_p] - Int_t Ndata_H_cal_4ta_emean; - Double_t H_cal_4ta_emean[13]; //[Ndata.H.cal.4ta.emean] - Int_t Ndata_H_cal_4ta_eneg; - Double_t H_cal_4ta_eneg[13]; //[Ndata.H.cal.4ta.eneg] - Int_t Ndata_H_cal_4ta_epos; - Double_t H_cal_4ta_epos[13]; //[Ndata.H.cal.4ta.epos] - Int_t Ndata_H_cal_4ta_negadchits; - Double_t H_cal_4ta_negadchits[1]; //[Ndata.H.cal.4ta.negadchits] - Int_t Ndata_H_cal_4ta_posadchits; - Double_t H_cal_4ta_posadchits[9]; //[Ndata.H.cal.4ta.posadchits] - Int_t Ndata_H_tr_beta; - Double_t H_tr_beta[10]; //[Ndata.H.tr.beta] - Int_t Ndata_H_tr_chi2; - Double_t H_tr_chi2[10]; //[Ndata.H.tr.chi2] - Int_t Ndata_H_tr_d_ph; - Double_t H_tr_d_ph[10]; //[Ndata.H.tr.d_ph] - Int_t Ndata_H_tr_d_th; - Double_t H_tr_d_th[10]; //[Ndata.H.tr.d_th] - Int_t Ndata_H_tr_d_x; - Double_t H_tr_d_x[10]; //[Ndata.H.tr.d_x] - Int_t Ndata_H_tr_d_y; - Double_t H_tr_d_y[10]; //[Ndata.H.tr.d_y] - Int_t Ndata_H_tr_dbeta; - Double_t H_tr_dbeta[10]; //[Ndata.H.tr.dbeta] - Int_t Ndata_H_tr_dtime; - Double_t H_tr_dtime[10]; //[Ndata.H.tr.dtime] - Int_t Ndata_H_tr_flag; - Double_t H_tr_flag[10]; //[Ndata.H.tr.flag] - Int_t Ndata_H_tr_ndof; - Double_t H_tr_ndof[10]; //[Ndata.H.tr.ndof] - Int_t Ndata_H_tr_p; - Double_t H_tr_p[10]; //[Ndata.H.tr.p] - Int_t Ndata_H_tr_pathl; - Double_t H_tr_pathl[10]; //[Ndata.H.tr.pathl] - Int_t Ndata_H_tr_ph; - Double_t H_tr_ph[10]; //[Ndata.H.tr.ph] - Int_t Ndata_H_tr_px; - Double_t H_tr_px[10]; //[Ndata.H.tr.px] - Int_t Ndata_H_tr_py; - Double_t H_tr_py[10]; //[Ndata.H.tr.py] - Int_t Ndata_H_tr_pz; - Double_t H_tr_pz[10]; //[Ndata.H.tr.pz] - Int_t Ndata_H_tr_r_ph; - Double_t H_tr_r_ph[10]; //[Ndata.H.tr.r_ph] - Int_t Ndata_H_tr_r_th; - Double_t H_tr_r_th[10]; //[Ndata.H.tr.r_th] - Int_t Ndata_H_tr_r_x; - Double_t H_tr_r_x[10]; //[Ndata.H.tr.r_x] - Int_t Ndata_H_tr_r_y; - Double_t H_tr_r_y[10]; //[Ndata.H.tr.r_y] - Int_t Ndata_H_tr_tg_dp; - Double_t H_tr_tg_dp[10]; //[Ndata.H.tr.tg_dp] - Int_t Ndata_H_tr_tg_ph; - Double_t H_tr_tg_ph[10]; //[Ndata.H.tr.tg_ph] - Int_t Ndata_H_tr_tg_th; - Double_t H_tr_tg_th[10]; //[Ndata.H.tr.tg_th] - Int_t Ndata_H_tr_tg_y; - Double_t H_tr_tg_y[10]; //[Ndata.H.tr.tg_y] - Int_t Ndata_H_tr_th; - Double_t H_tr_th[10]; //[Ndata.H.tr.th] - Int_t Ndata_H_tr_time; - Double_t H_tr_time[10]; //[Ndata.H.tr.time] - Int_t Ndata_H_tr_vx; - Double_t H_tr_vx[10]; //[Ndata.H.tr.vx] - Int_t Ndata_H_tr_vy; - Double_t H_tr_vy[10]; //[Ndata.H.tr.vy] - Int_t Ndata_H_tr_vz; - Double_t H_tr_vz[10]; //[Ndata.H.tr.vz] - Int_t Ndata_H_tr_x; - Double_t H_tr_x[10]; //[Ndata.H.tr.x] - Int_t Ndata_H_tr_y; - Double_t H_tr_y[10]; //[Ndata.H.tr.y] - Double_t H_cal_1pr_eplane; - Double_t H_cal_1pr_eplane_neg; - Double_t H_cal_1pr_eplane_pos; - Double_t H_cal_2ta_eplane; - Double_t H_cal_2ta_eplane_neg; - Double_t H_cal_2ta_eplane_pos; - Double_t H_cal_3ta_eplane; - Double_t H_cal_3ta_eplane_neg; - Double_t H_cal_3ta_eplane_pos; - Double_t H_cal_4ta_eplane; - Double_t H_cal_4ta_eplane_neg; - Double_t H_cal_4ta_eplane_pos; - Double_t H_cal_nclust; - Double_t H_cal_nhits; - Double_t H_cal_ntracks; - Double_t H_gold_dp; - Double_t H_gold_e; - Double_t H_gold_fp_ph; - Double_t H_gold_fp_th; - Double_t H_gold_fp_x; - Double_t H_gold_fp_y; - Double_t H_gold_index; - Double_t H_gold_ok; - Double_t H_gold_p; - Double_t H_gold_ph; - Double_t H_gold_px; - Double_t H_gold_py; - Double_t H_gold_pz; - Double_t H_gold_th; - Double_t H_gold_x; - Double_t H_gold_y; - Double_t H_tr_n; - Double_t g_evtyp; - //THaEvent *Event_Branch; - ULong64_t fEvtHdr_fEvtTime; - UInt_t fEvtHdr_fEvtNum; - Int_t fEvtHdr_fEvtType; - Int_t fEvtHdr_fEvtLen; - Int_t fEvtHdr_fHelicity; - Int_t fEvtHdr_fTargetPol; - Int_t fEvtHdr_fRun; - - // List of branches - TBranch *b_Ndata_H_cal_1pr_aneg; //! - TBranch *b_H_cal_1pr_aneg; //! - TBranch *b_Ndata_H_cal_1pr_aneg_p; //! - TBranch *b_H_cal_1pr_aneg_p; //! - TBranch *b_Ndata_H_cal_1pr_apos; //! - TBranch *b_H_cal_1pr_apos; //! - TBranch *b_Ndata_H_cal_1pr_apos_p; //! - TBranch *b_H_cal_1pr_apos_p; //! - TBranch *b_Ndata_H_cal_1pr_emean; //! - TBranch *b_H_cal_1pr_emean; //! - TBranch *b_Ndata_H_cal_1pr_eneg; //! - TBranch *b_H_cal_1pr_eneg; //! - TBranch *b_Ndata_H_cal_1pr_epos; //! - TBranch *b_H_cal_1pr_epos; //! - TBranch *b_Ndata_H_cal_1pr_negadchits; //! - TBranch *b_H_cal_1pr_negadchits; //! - TBranch *b_Ndata_H_cal_1pr_posadchits; //! - TBranch *b_H_cal_1pr_posadchits; //! - TBranch *b_Ndata_H_cal_2ta_aneg; //! - TBranch *b_H_cal_2ta_aneg; //! - TBranch *b_Ndata_H_cal_2ta_aneg_p; //! - TBranch *b_H_cal_2ta_aneg_p; //! - TBranch *b_Ndata_H_cal_2ta_apos; //! - TBranch *b_H_cal_2ta_apos; //! - TBranch *b_Ndata_H_cal_2ta_apos_p; //! - TBranch *b_H_cal_2ta_apos_p; //! - TBranch *b_Ndata_H_cal_2ta_emean; //! - TBranch *b_H_cal_2ta_emean; //! - TBranch *b_Ndata_H_cal_2ta_eneg; //! - TBranch *b_H_cal_2ta_eneg; //! - TBranch *b_Ndata_H_cal_2ta_epos; //! - TBranch *b_H_cal_2ta_epos; //! - TBranch *b_Ndata_H_cal_2ta_negadchits; //! - TBranch *b_H_cal_2ta_negadchits; //! - TBranch *b_Ndata_H_cal_2ta_posadchits; //! - TBranch *b_H_cal_2ta_posadchits; //! - TBranch *b_Ndata_H_cal_3ta_aneg; //! - TBranch *b_H_cal_3ta_aneg; //! - TBranch *b_Ndata_H_cal_3ta_aneg_p; //! - TBranch *b_H_cal_3ta_aneg_p; //! - TBranch *b_Ndata_H_cal_3ta_apos; //! - TBranch *b_H_cal_3ta_apos; //! - TBranch *b_Ndata_H_cal_3ta_apos_p; //! - TBranch *b_H_cal_3ta_apos_p; //! - TBranch *b_Ndata_H_cal_3ta_emean; //! - TBranch *b_H_cal_3ta_emean; //! - TBranch *b_Ndata_H_cal_3ta_eneg; //! - TBranch *b_H_cal_3ta_eneg; //! - TBranch *b_Ndata_H_cal_3ta_epos; //! - TBranch *b_H_cal_3ta_epos; //! - TBranch *b_Ndata_H_cal_3ta_negadchits; //! - TBranch *b_H_cal_3ta_negadchits; //! - TBranch *b_Ndata_H_cal_3ta_posadchits; //! - TBranch *b_H_cal_3ta_posadchits; //! - TBranch *b_Ndata_H_cal_4ta_aneg; //! - TBranch *b_H_cal_4ta_aneg; //! - TBranch *b_Ndata_H_cal_4ta_aneg_p; //! - TBranch *b_H_cal_4ta_aneg_p; //! - TBranch *b_Ndata_H_cal_4ta_apos; //! - TBranch *b_H_cal_4ta_apos; //! - TBranch *b_Ndata_H_cal_4ta_apos_p; //! - TBranch *b_H_cal_4ta_apos_p; //! - TBranch *b_Ndata_H_cal_4ta_emean; //! - TBranch *b_H_cal_4ta_emean; //! - TBranch *b_Ndata_H_cal_4ta_eneg; //! - TBranch *b_H_cal_4ta_eneg; //! - TBranch *b_Ndata_H_cal_4ta_epos; //! - TBranch *b_H_cal_4ta_epos; //! - TBranch *b_Ndata_H_cal_4ta_negadchits; //! - TBranch *b_H_cal_4ta_negadchits; //! - TBranch *b_Ndata_H_cal_4ta_posadchits; //! - TBranch *b_H_cal_4ta_posadchits; //! - TBranch *b_Ndata_H_tr_beta; //! - TBranch *b_H_tr_beta; //! - TBranch *b_Ndata_H_tr_chi2; //! - TBranch *b_H_tr_chi2; //! - TBranch *b_Ndata_H_tr_d_ph; //! - TBranch *b_H_tr_d_ph; //! - TBranch *b_Ndata_H_tr_d_th; //! - TBranch *b_H_tr_d_th; //! - TBranch *b_Ndata_H_tr_d_x; //! - TBranch *b_H_tr_d_x; //! - TBranch *b_Ndata_H_tr_d_y; //! - TBranch *b_H_tr_d_y; //! - TBranch *b_Ndata_H_tr_dbeta; //! - TBranch *b_H_tr_dbeta; //! - TBranch *b_Ndata_H_tr_dtime; //! - TBranch *b_H_tr_dtime; //! - TBranch *b_Ndata_H_tr_flag; //! - TBranch *b_H_tr_flag; //! - TBranch *b_Ndata_H_tr_ndof; //! - TBranch *b_H_tr_ndof; //! - TBranch *b_Ndata_H_tr_p; //! - TBranch *b_H_tr_p; //! - TBranch *b_Ndata_H_tr_pathl; //! - TBranch *b_H_tr_pathl; //! - TBranch *b_Ndata_H_tr_ph; //! - TBranch *b_H_tr_ph; //! - TBranch *b_Ndata_H_tr_px; //! - TBranch *b_H_tr_px; //! - TBranch *b_Ndata_H_tr_py; //! - TBranch *b_H_tr_py; //! - TBranch *b_Ndata_H_tr_pz; //! - TBranch *b_H_tr_pz; //! - TBranch *b_Ndata_H_tr_r_ph; //! - TBranch *b_H_tr_r_ph; //! - TBranch *b_Ndata_H_tr_r_th; //! - TBranch *b_H_tr_r_th; //! - TBranch *b_Ndata_H_tr_r_x; //! - TBranch *b_H_tr_r_x; //! - TBranch *b_Ndata_H_tr_r_y; //! - TBranch *b_H_tr_r_y; //! - TBranch *b_Ndata_H_tr_tg_dp; //! - TBranch *b_H_tr_tg_dp; //! - TBranch *b_Ndata_H_tr_tg_ph; //! - TBranch *b_H_tr_tg_ph; //! - TBranch *b_Ndata_H_tr_tg_th; //! - TBranch *b_H_tr_tg_th; //! - TBranch *b_Ndata_H_tr_tg_y; //! - TBranch *b_H_tr_tg_y; //! - TBranch *b_Ndata_H_tr_th; //! - TBranch *b_H_tr_th; //! - TBranch *b_Ndata_H_tr_time; //! - TBranch *b_H_tr_time; //! - TBranch *b_Ndata_H_tr_vx; //! - TBranch *b_H_tr_vx; //! - TBranch *b_Ndata_H_tr_vy; //! - TBranch *b_H_tr_vy; //! - TBranch *b_Ndata_H_tr_vz; //! - TBranch *b_H_tr_vz; //! - TBranch *b_Ndata_H_tr_x; //! - TBranch *b_H_tr_x; //! - TBranch *b_Ndata_H_tr_y; //! - TBranch *b_H_tr_y; //! - TBranch *b_H_cal_1pr_eplane; //! - TBranch *b_H_cal_1pr_eplane_neg; //! - TBranch *b_H_cal_1pr_eplane_pos; //! - TBranch *b_H_cal_2ta_eplane; //! - TBranch *b_H_cal_2ta_eplane_neg; //! - TBranch *b_H_cal_2ta_eplane_pos; //! - TBranch *b_H_cal_3ta_eplane; //! - TBranch *b_H_cal_3ta_eplane_neg; //! - TBranch *b_H_cal_3ta_eplane_pos; //! - TBranch *b_H_cal_4ta_eplane; //! - TBranch *b_H_cal_4ta_eplane_neg; //! - TBranch *b_H_cal_4ta_eplane_pos; //! - TBranch *b_H_cal_nclust; //! - TBranch *b_H_cal_nhits; //! - TBranch *b_H_cal_ntracks; //! - TBranch *b_H_gold_dp; //! - TBranch *b_H_gold_e; //! - TBranch *b_H_gold_fp_ph; //! - TBranch *b_H_gold_fp_th; //! - TBranch *b_H_gold_fp_x; //! - TBranch *b_H_gold_fp_y; //! - TBranch *b_H_gold_index; //! - TBranch *b_H_gold_ok; //! - TBranch *b_H_gold_p; //! - TBranch *b_H_gold_ph; //! - TBranch *b_H_gold_px; //! - TBranch *b_H_gold_py; //! - TBranch *b_H_gold_pz; //! - TBranch *b_H_gold_th; //! - TBranch *b_H_gold_x; //! - TBranch *b_H_gold_y; //! - TBranch *b_H_tr_n; //! - TBranch *b_g_evtyp; //! - TBranch *b_Event_Branch_fEvtHdr_fEvtTime; //! - TBranch *b_Event_Branch_fEvtHdr_fEvtNum; //! - TBranch *b_Event_Branch_fEvtHdr_fEvtType; //! - TBranch *b_Event_Branch_fEvtHdr_fEvtLen; //! - TBranch *b_Event_Branch_fEvtHdr_fHelicity; //! - TBranch *b_Event_Branch_fEvtHdr_fTargetPol; //! - TBranch *b_Event_Branch_fEvtHdr_fRun; //! - - hcana_class(TTree *tree=0); - virtual ~hcana_class(); - virtual Int_t Cut(Long64_t entry); - virtual Int_t GetEntry(Long64_t entry); - virtual Long64_t LoadTree(Long64_t entry); - virtual void Init(TTree *tree); - virtual void Loop(); - virtual Bool_t Notify(); - virtual void Show(Long64_t entry = -1); -}; - -#endif - -#ifdef hcana_class_cxx -hcana_class::hcana_class(TTree *tree) : fChain(0) -{ -// if parameter tree is not specified (or zero), connect the file -// used to generate this class and read the Tree. - if (tree == 0) { - TFile *f = (TFile*)gROOT->GetListOfFiles()->FindObject("hodtest_52949.root"); - if (!f || !f->IsOpen()) { - f = new TFile("hodtest_52949.root"); - } - f->GetObject("T",tree); - - } - Init(tree); -} - -hcana_class::~hcana_class() -{ - if (!fChain) return; - delete fChain->GetCurrentFile(); -} - -Int_t hcana_class::GetEntry(Long64_t entry) -{ -// Read contents of entry. - if (!fChain) return 0; - return fChain->GetEntry(entry); -} -Long64_t hcana_class::LoadTree(Long64_t entry) -{ -// Set the environment to read one entry - if (!fChain) return -5; - Long64_t centry = fChain->LoadTree(entry); - if (centry < 0) return centry; - if (fChain->GetTreeNumber() != fCurrent) { - fCurrent = fChain->GetTreeNumber(); - Notify(); - } - return centry; -} - -void hcana_class::Init(TTree *tree) -{ - // The Init() function is called when the selector needs to initialize - // a new tree or chain. Typically here the branch addresses and branch - // pointers of the tree will be set. - // It is normally not necessary to make changes to the generated - // code, but the routine can be extended by the user if needed. - // Init() will be called many times when running on PROOF - // (once per file to be processed). - - // Set branch addresses and branch pointers - if (!tree) return; - fChain = tree; - fCurrent = -1; - fChain->SetMakeClass(1); - - fChain->SetBranchAddress("Ndata.H.cal.1pr.aneg", &Ndata_H_cal_1pr_aneg, &b_Ndata_H_cal_1pr_aneg); - fChain->SetBranchAddress("H.cal.1pr.aneg", H_cal_1pr_aneg, &b_H_cal_1pr_aneg); - fChain->SetBranchAddress("Ndata.H.cal.1pr.aneg_p", &Ndata_H_cal_1pr_aneg_p, &b_Ndata_H_cal_1pr_aneg_p); - fChain->SetBranchAddress("H.cal.1pr.aneg_p", H_cal_1pr_aneg_p, &b_H_cal_1pr_aneg_p); - fChain->SetBranchAddress("Ndata.H.cal.1pr.apos", &Ndata_H_cal_1pr_apos, &b_Ndata_H_cal_1pr_apos); - fChain->SetBranchAddress("H.cal.1pr.apos", H_cal_1pr_apos, &b_H_cal_1pr_apos); - fChain->SetBranchAddress("Ndata.H.cal.1pr.apos_p", &Ndata_H_cal_1pr_apos_p, &b_Ndata_H_cal_1pr_apos_p); - fChain->SetBranchAddress("H.cal.1pr.apos_p", H_cal_1pr_apos_p, &b_H_cal_1pr_apos_p); - fChain->SetBranchAddress("Ndata.H.cal.1pr.emean", &Ndata_H_cal_1pr_emean, &b_Ndata_H_cal_1pr_emean); - fChain->SetBranchAddress("H.cal.1pr.emean", H_cal_1pr_emean, &b_H_cal_1pr_emean); - fChain->SetBranchAddress("Ndata.H.cal.1pr.eneg", &Ndata_H_cal_1pr_eneg, &b_Ndata_H_cal_1pr_eneg); - fChain->SetBranchAddress("H.cal.1pr.eneg", H_cal_1pr_eneg, &b_H_cal_1pr_eneg); - fChain->SetBranchAddress("Ndata.H.cal.1pr.epos", &Ndata_H_cal_1pr_epos, &b_Ndata_H_cal_1pr_epos); - fChain->SetBranchAddress("H.cal.1pr.epos", H_cal_1pr_epos, &b_H_cal_1pr_epos); - fChain->SetBranchAddress("Ndata.H.cal.1pr.negadchits", &Ndata_H_cal_1pr_negadchits, &b_Ndata_H_cal_1pr_negadchits); - fChain->SetBranchAddress("H.cal.1pr.negadchits", H_cal_1pr_negadchits, &b_H_cal_1pr_negadchits); - fChain->SetBranchAddress("Ndata.H.cal.1pr.posadchits", &Ndata_H_cal_1pr_posadchits, &b_Ndata_H_cal_1pr_posadchits); - fChain->SetBranchAddress("H.cal.1pr.posadchits", H_cal_1pr_posadchits, &b_H_cal_1pr_posadchits); - fChain->SetBranchAddress("Ndata.H.cal.2ta.aneg", &Ndata_H_cal_2ta_aneg, &b_Ndata_H_cal_2ta_aneg); - fChain->SetBranchAddress("H.cal.2ta.aneg", H_cal_2ta_aneg, &b_H_cal_2ta_aneg); - fChain->SetBranchAddress("Ndata.H.cal.2ta.aneg_p", &Ndata_H_cal_2ta_aneg_p, &b_Ndata_H_cal_2ta_aneg_p); - fChain->SetBranchAddress("H.cal.2ta.aneg_p", H_cal_2ta_aneg_p, &b_H_cal_2ta_aneg_p); - fChain->SetBranchAddress("Ndata.H.cal.2ta.apos", &Ndata_H_cal_2ta_apos, &b_Ndata_H_cal_2ta_apos); - fChain->SetBranchAddress("H.cal.2ta.apos", H_cal_2ta_apos, &b_H_cal_2ta_apos); - fChain->SetBranchAddress("Ndata.H.cal.2ta.apos_p", &Ndata_H_cal_2ta_apos_p, &b_Ndata_H_cal_2ta_apos_p); - fChain->SetBranchAddress("H.cal.2ta.apos_p", H_cal_2ta_apos_p, &b_H_cal_2ta_apos_p); - fChain->SetBranchAddress("Ndata.H.cal.2ta.emean", &Ndata_H_cal_2ta_emean, &b_Ndata_H_cal_2ta_emean); - fChain->SetBranchAddress("H.cal.2ta.emean", H_cal_2ta_emean, &b_H_cal_2ta_emean); - fChain->SetBranchAddress("Ndata.H.cal.2ta.eneg", &Ndata_H_cal_2ta_eneg, &b_Ndata_H_cal_2ta_eneg); - fChain->SetBranchAddress("H.cal.2ta.eneg", H_cal_2ta_eneg, &b_H_cal_2ta_eneg); - fChain->SetBranchAddress("Ndata.H.cal.2ta.epos", &Ndata_H_cal_2ta_epos, &b_Ndata_H_cal_2ta_epos); - fChain->SetBranchAddress("H.cal.2ta.epos", H_cal_2ta_epos, &b_H_cal_2ta_epos); - fChain->SetBranchAddress("Ndata.H.cal.2ta.negadchits", &Ndata_H_cal_2ta_negadchits, &b_Ndata_H_cal_2ta_negadchits); - fChain->SetBranchAddress("H.cal.2ta.negadchits", H_cal_2ta_negadchits, &b_H_cal_2ta_negadchits); - fChain->SetBranchAddress("Ndata.H.cal.2ta.posadchits", &Ndata_H_cal_2ta_posadchits, &b_Ndata_H_cal_2ta_posadchits); - fChain->SetBranchAddress("H.cal.2ta.posadchits", H_cal_2ta_posadchits, &b_H_cal_2ta_posadchits); - fChain->SetBranchAddress("Ndata.H.cal.3ta.aneg", &Ndata_H_cal_3ta_aneg, &b_Ndata_H_cal_3ta_aneg); - fChain->SetBranchAddress("H.cal.3ta.aneg", H_cal_3ta_aneg, &b_H_cal_3ta_aneg); - fChain->SetBranchAddress("Ndata.H.cal.3ta.aneg_p", &Ndata_H_cal_3ta_aneg_p, &b_Ndata_H_cal_3ta_aneg_p); - fChain->SetBranchAddress("H.cal.3ta.aneg_p", H_cal_3ta_aneg_p, &b_H_cal_3ta_aneg_p); - fChain->SetBranchAddress("Ndata.H.cal.3ta.apos", &Ndata_H_cal_3ta_apos, &b_Ndata_H_cal_3ta_apos); - fChain->SetBranchAddress("H.cal.3ta.apos", H_cal_3ta_apos, &b_H_cal_3ta_apos); - fChain->SetBranchAddress("Ndata.H.cal.3ta.apos_p", &Ndata_H_cal_3ta_apos_p, &b_Ndata_H_cal_3ta_apos_p); - fChain->SetBranchAddress("H.cal.3ta.apos_p", H_cal_3ta_apos_p, &b_H_cal_3ta_apos_p); - fChain->SetBranchAddress("Ndata.H.cal.3ta.emean", &Ndata_H_cal_3ta_emean, &b_Ndata_H_cal_3ta_emean); - fChain->SetBranchAddress("H.cal.3ta.emean", H_cal_3ta_emean, &b_H_cal_3ta_emean); - fChain->SetBranchAddress("Ndata.H.cal.3ta.eneg", &Ndata_H_cal_3ta_eneg, &b_Ndata_H_cal_3ta_eneg); - fChain->SetBranchAddress("H.cal.3ta.eneg", H_cal_3ta_eneg, &b_H_cal_3ta_eneg); - fChain->SetBranchAddress("Ndata.H.cal.3ta.epos", &Ndata_H_cal_3ta_epos, &b_Ndata_H_cal_3ta_epos); - fChain->SetBranchAddress("H.cal.3ta.epos", H_cal_3ta_epos, &b_H_cal_3ta_epos); - fChain->SetBranchAddress("Ndata.H.cal.3ta.negadchits", &Ndata_H_cal_3ta_negadchits, &b_Ndata_H_cal_3ta_negadchits); - fChain->SetBranchAddress("H.cal.3ta.negadchits", &H_cal_3ta_negadchits, &b_H_cal_3ta_negadchits); - fChain->SetBranchAddress("Ndata.H.cal.3ta.posadchits", &Ndata_H_cal_3ta_posadchits, &b_Ndata_H_cal_3ta_posadchits); - fChain->SetBranchAddress("H.cal.3ta.posadchits", H_cal_3ta_posadchits, &b_H_cal_3ta_posadchits); - fChain->SetBranchAddress("Ndata.H.cal.4ta.aneg", &Ndata_H_cal_4ta_aneg, &b_Ndata_H_cal_4ta_aneg); - fChain->SetBranchAddress("H.cal.4ta.aneg", H_cal_4ta_aneg, &b_H_cal_4ta_aneg); - fChain->SetBranchAddress("Ndata.H.cal.4ta.aneg_p", &Ndata_H_cal_4ta_aneg_p, &b_Ndata_H_cal_4ta_aneg_p); - fChain->SetBranchAddress("H.cal.4ta.aneg_p", H_cal_4ta_aneg_p, &b_H_cal_4ta_aneg_p); - fChain->SetBranchAddress("Ndata.H.cal.4ta.apos", &Ndata_H_cal_4ta_apos, &b_Ndata_H_cal_4ta_apos); - fChain->SetBranchAddress("H.cal.4ta.apos", H_cal_4ta_apos, &b_H_cal_4ta_apos); - fChain->SetBranchAddress("Ndata.H.cal.4ta.apos_p", &Ndata_H_cal_4ta_apos_p, &b_Ndata_H_cal_4ta_apos_p); - fChain->SetBranchAddress("H.cal.4ta.apos_p", H_cal_4ta_apos_p, &b_H_cal_4ta_apos_p); - fChain->SetBranchAddress("Ndata.H.cal.4ta.emean", &Ndata_H_cal_4ta_emean, &b_Ndata_H_cal_4ta_emean); - fChain->SetBranchAddress("H.cal.4ta.emean", H_cal_4ta_emean, &b_H_cal_4ta_emean); - fChain->SetBranchAddress("Ndata.H.cal.4ta.eneg", &Ndata_H_cal_4ta_eneg, &b_Ndata_H_cal_4ta_eneg); - fChain->SetBranchAddress("H.cal.4ta.eneg", H_cal_4ta_eneg, &b_H_cal_4ta_eneg); - fChain->SetBranchAddress("Ndata.H.cal.4ta.epos", &Ndata_H_cal_4ta_epos, &b_Ndata_H_cal_4ta_epos); - fChain->SetBranchAddress("H.cal.4ta.epos", H_cal_4ta_epos, &b_H_cal_4ta_epos); - fChain->SetBranchAddress("Ndata.H.cal.4ta.negadchits", &Ndata_H_cal_4ta_negadchits, &b_Ndata_H_cal_4ta_negadchits); - fChain->SetBranchAddress("H.cal.4ta.negadchits", &H_cal_4ta_negadchits, &b_H_cal_4ta_negadchits); - fChain->SetBranchAddress("Ndata.H.cal.4ta.posadchits", &Ndata_H_cal_4ta_posadchits, &b_Ndata_H_cal_4ta_posadchits); - fChain->SetBranchAddress("H.cal.4ta.posadchits", H_cal_4ta_posadchits, &b_H_cal_4ta_posadchits); - fChain->SetBranchAddress("Ndata.H.tr.beta", &Ndata_H_tr_beta, &b_Ndata_H_tr_beta); - fChain->SetBranchAddress("H.tr.beta", H_tr_beta, &b_H_tr_beta); - fChain->SetBranchAddress("Ndata.H.tr.chi2", &Ndata_H_tr_chi2, &b_Ndata_H_tr_chi2); - fChain->SetBranchAddress("H.tr.chi2", H_tr_chi2, &b_H_tr_chi2); - fChain->SetBranchAddress("Ndata.H.tr.d_ph", &Ndata_H_tr_d_ph, &b_Ndata_H_tr_d_ph); - fChain->SetBranchAddress("H.tr.d_ph", H_tr_d_ph, &b_H_tr_d_ph); - fChain->SetBranchAddress("Ndata.H.tr.d_th", &Ndata_H_tr_d_th, &b_Ndata_H_tr_d_th); - fChain->SetBranchAddress("H.tr.d_th", H_tr_d_th, &b_H_tr_d_th); - fChain->SetBranchAddress("Ndata.H.tr.d_x", &Ndata_H_tr_d_x, &b_Ndata_H_tr_d_x); - fChain->SetBranchAddress("H.tr.d_x", H_tr_d_x, &b_H_tr_d_x); - fChain->SetBranchAddress("Ndata.H.tr.d_y", &Ndata_H_tr_d_y, &b_Ndata_H_tr_d_y); - fChain->SetBranchAddress("H.tr.d_y", H_tr_d_y, &b_H_tr_d_y); - fChain->SetBranchAddress("Ndata.H.tr.dbeta", &Ndata_H_tr_dbeta, &b_Ndata_H_tr_dbeta); - fChain->SetBranchAddress("H.tr.dbeta", H_tr_dbeta, &b_H_tr_dbeta); - fChain->SetBranchAddress("Ndata.H.tr.dtime", &Ndata_H_tr_dtime, &b_Ndata_H_tr_dtime); - fChain->SetBranchAddress("H.tr.dtime", H_tr_dtime, &b_H_tr_dtime); - fChain->SetBranchAddress("Ndata.H.tr.flag", &Ndata_H_tr_flag, &b_Ndata_H_tr_flag); - fChain->SetBranchAddress("H.tr.flag", H_tr_flag, &b_H_tr_flag); - fChain->SetBranchAddress("Ndata.H.tr.ndof", &Ndata_H_tr_ndof, &b_Ndata_H_tr_ndof); - fChain->SetBranchAddress("H.tr.ndof", H_tr_ndof, &b_H_tr_ndof); - fChain->SetBranchAddress("Ndata.H.tr.p", &Ndata_H_tr_p, &b_Ndata_H_tr_p); - fChain->SetBranchAddress("H.tr.p", H_tr_p, &b_H_tr_p); - fChain->SetBranchAddress("Ndata.H.tr.pathl", &Ndata_H_tr_pathl, &b_Ndata_H_tr_pathl); - fChain->SetBranchAddress("H.tr.pathl", H_tr_pathl, &b_H_tr_pathl); - fChain->SetBranchAddress("Ndata.H.tr.ph", &Ndata_H_tr_ph, &b_Ndata_H_tr_ph); - fChain->SetBranchAddress("H.tr.ph", H_tr_ph, &b_H_tr_ph); - fChain->SetBranchAddress("Ndata.H.tr.px", &Ndata_H_tr_px, &b_Ndata_H_tr_px); - fChain->SetBranchAddress("H.tr.px", H_tr_px, &b_H_tr_px); - fChain->SetBranchAddress("Ndata.H.tr.py", &Ndata_H_tr_py, &b_Ndata_H_tr_py); - fChain->SetBranchAddress("H.tr.py", H_tr_py, &b_H_tr_py); - fChain->SetBranchAddress("Ndata.H.tr.pz", &Ndata_H_tr_pz, &b_Ndata_H_tr_pz); - fChain->SetBranchAddress("H.tr.pz", H_tr_pz, &b_H_tr_pz); - fChain->SetBranchAddress("Ndata.H.tr.r_ph", &Ndata_H_tr_r_ph, &b_Ndata_H_tr_r_ph); - fChain->SetBranchAddress("H.tr.r_ph", H_tr_r_ph, &b_H_tr_r_ph); - fChain->SetBranchAddress("Ndata.H.tr.r_th", &Ndata_H_tr_r_th, &b_Ndata_H_tr_r_th); - fChain->SetBranchAddress("H.tr.r_th", H_tr_r_th, &b_H_tr_r_th); - fChain->SetBranchAddress("Ndata.H.tr.r_x", &Ndata_H_tr_r_x, &b_Ndata_H_tr_r_x); - fChain->SetBranchAddress("H.tr.r_x", H_tr_r_x, &b_H_tr_r_x); - fChain->SetBranchAddress("Ndata.H.tr.r_y", &Ndata_H_tr_r_y, &b_Ndata_H_tr_r_y); - fChain->SetBranchAddress("H.tr.r_y", H_tr_r_y, &b_H_tr_r_y); - fChain->SetBranchAddress("Ndata.H.tr.tg_dp", &Ndata_H_tr_tg_dp, &b_Ndata_H_tr_tg_dp); - fChain->SetBranchAddress("H.tr.tg_dp", H_tr_tg_dp, &b_H_tr_tg_dp); - fChain->SetBranchAddress("Ndata.H.tr.tg_ph", &Ndata_H_tr_tg_ph, &b_Ndata_H_tr_tg_ph); - fChain->SetBranchAddress("H.tr.tg_ph", H_tr_tg_ph, &b_H_tr_tg_ph); - fChain->SetBranchAddress("Ndata.H.tr.tg_th", &Ndata_H_tr_tg_th, &b_Ndata_H_tr_tg_th); - fChain->SetBranchAddress("H.tr.tg_th", H_tr_tg_th, &b_H_tr_tg_th); - fChain->SetBranchAddress("Ndata.H.tr.tg_y", &Ndata_H_tr_tg_y, &b_Ndata_H_tr_tg_y); - fChain->SetBranchAddress("H.tr.tg_y", H_tr_tg_y, &b_H_tr_tg_y); - fChain->SetBranchAddress("Ndata.H.tr.th", &Ndata_H_tr_th, &b_Ndata_H_tr_th); - fChain->SetBranchAddress("H.tr.th", H_tr_th, &b_H_tr_th); - fChain->SetBranchAddress("Ndata.H.tr.time", &Ndata_H_tr_time, &b_Ndata_H_tr_time); - fChain->SetBranchAddress("H.tr.time", H_tr_time, &b_H_tr_time); - fChain->SetBranchAddress("Ndata.H.tr.vx", &Ndata_H_tr_vx, &b_Ndata_H_tr_vx); - fChain->SetBranchAddress("H.tr.vx", H_tr_vx, &b_H_tr_vx); - fChain->SetBranchAddress("Ndata.H.tr.vy", &Ndata_H_tr_vy, &b_Ndata_H_tr_vy); - fChain->SetBranchAddress("H.tr.vy", H_tr_vy, &b_H_tr_vy); - fChain->SetBranchAddress("Ndata.H.tr.vz", &Ndata_H_tr_vz, &b_Ndata_H_tr_vz); - fChain->SetBranchAddress("H.tr.vz", H_tr_vz, &b_H_tr_vz); - fChain->SetBranchAddress("Ndata.H.tr.x", &Ndata_H_tr_x, &b_Ndata_H_tr_x); - fChain->SetBranchAddress("H.tr.x", H_tr_x, &b_H_tr_x); - fChain->SetBranchAddress("Ndata.H.tr.y", &Ndata_H_tr_y, &b_Ndata_H_tr_y); - fChain->SetBranchAddress("H.tr.y", H_tr_y, &b_H_tr_y); - fChain->SetBranchAddress("H.cal.1pr.eplane", &H_cal_1pr_eplane, &b_H_cal_1pr_eplane); - fChain->SetBranchAddress("H.cal.1pr.eplane_neg", &H_cal_1pr_eplane_neg, &b_H_cal_1pr_eplane_neg); - fChain->SetBranchAddress("H.cal.1pr.eplane_pos", &H_cal_1pr_eplane_pos, &b_H_cal_1pr_eplane_pos); - fChain->SetBranchAddress("H.cal.2ta.eplane", &H_cal_2ta_eplane, &b_H_cal_2ta_eplane); - fChain->SetBranchAddress("H.cal.2ta.eplane_neg", &H_cal_2ta_eplane_neg, &b_H_cal_2ta_eplane_neg); - fChain->SetBranchAddress("H.cal.2ta.eplane_pos", &H_cal_2ta_eplane_pos, &b_H_cal_2ta_eplane_pos); - fChain->SetBranchAddress("H.cal.3ta.eplane", &H_cal_3ta_eplane, &b_H_cal_3ta_eplane); - fChain->SetBranchAddress("H.cal.3ta.eplane_neg", &H_cal_3ta_eplane_neg, &b_H_cal_3ta_eplane_neg); - fChain->SetBranchAddress("H.cal.3ta.eplane_pos", &H_cal_3ta_eplane_pos, &b_H_cal_3ta_eplane_pos); - fChain->SetBranchAddress("H.cal.4ta.eplane", &H_cal_4ta_eplane, &b_H_cal_4ta_eplane); - fChain->SetBranchAddress("H.cal.4ta.eplane_neg", &H_cal_4ta_eplane_neg, &b_H_cal_4ta_eplane_neg); - fChain->SetBranchAddress("H.cal.4ta.eplane_pos", &H_cal_4ta_eplane_pos, &b_H_cal_4ta_eplane_pos); - fChain->SetBranchAddress("H.cal.nclust", &H_cal_nclust, &b_H_cal_nclust); - fChain->SetBranchAddress("H.cal.nhits", &H_cal_nhits, &b_H_cal_nhits); - fChain->SetBranchAddress("H.cal.ntracks", &H_cal_ntracks, &b_H_cal_ntracks); - fChain->SetBranchAddress("H.gold.dp", &H_gold_dp, &b_H_gold_dp); - fChain->SetBranchAddress("H.gold.e", &H_gold_e, &b_H_gold_e); - fChain->SetBranchAddress("H.gold.fp_ph", &H_gold_fp_ph, &b_H_gold_fp_ph); - fChain->SetBranchAddress("H.gold.fp_th", &H_gold_fp_th, &b_H_gold_fp_th); - fChain->SetBranchAddress("H.gold.fp_x", &H_gold_fp_x, &b_H_gold_fp_x); - fChain->SetBranchAddress("H.gold.fp_y", &H_gold_fp_y, &b_H_gold_fp_y); - fChain->SetBranchAddress("H.gold.index", &H_gold_index, &b_H_gold_index); - fChain->SetBranchAddress("H.gold.ok", &H_gold_ok, &b_H_gold_ok); - fChain->SetBranchAddress("H.gold.p", &H_gold_p, &b_H_gold_p); - fChain->SetBranchAddress("H.gold.ph", &H_gold_ph, &b_H_gold_ph); - fChain->SetBranchAddress("H.gold.px", &H_gold_px, &b_H_gold_px); - fChain->SetBranchAddress("H.gold.py", &H_gold_py, &b_H_gold_py); - fChain->SetBranchAddress("H.gold.pz", &H_gold_pz, &b_H_gold_pz); - fChain->SetBranchAddress("H.gold.th", &H_gold_th, &b_H_gold_th); - fChain->SetBranchAddress("H.gold.x", &H_gold_x, &b_H_gold_x); - fChain->SetBranchAddress("H.gold.y", &H_gold_y, &b_H_gold_y); - fChain->SetBranchAddress("H.tr.n", &H_tr_n, &b_H_tr_n); - fChain->SetBranchAddress("g.evtyp", &g_evtyp, &b_g_evtyp); - fChain->SetBranchAddress("fEvtHdr.fEvtTime", &fEvtHdr_fEvtTime, &b_Event_Branch_fEvtHdr_fEvtTime); - fChain->SetBranchAddress("fEvtHdr.fEvtNum", &fEvtHdr_fEvtNum, &b_Event_Branch_fEvtHdr_fEvtNum); - fChain->SetBranchAddress("fEvtHdr.fEvtType", &fEvtHdr_fEvtType, &b_Event_Branch_fEvtHdr_fEvtType); - fChain->SetBranchAddress("fEvtHdr.fEvtLen", &fEvtHdr_fEvtLen, &b_Event_Branch_fEvtHdr_fEvtLen); - fChain->SetBranchAddress("fEvtHdr.fHelicity", &fEvtHdr_fHelicity, &b_Event_Branch_fEvtHdr_fHelicity); - fChain->SetBranchAddress("fEvtHdr.fTargetPol", &fEvtHdr_fTargetPol, &b_Event_Branch_fEvtHdr_fTargetPol); - fChain->SetBranchAddress("fEvtHdr.fRun", &fEvtHdr_fRun, &b_Event_Branch_fEvtHdr_fRun); - Notify(); -} - -Bool_t hcana_class::Notify() -{ - // The Notify() function is called when a new file is opened. This - // can be either for a new TTree in a TChain or when when a new TTree - // is started when using PROOF. It is normally not necessary to make changes - // to the generated code, but the routine can be extended by the - // user if needed. The return value is currently not used. - - return kTRUE; -} - -void hcana_class::Show(Long64_t entry) -{ -// Print contents of entry. -// If entry is not specified, print current entry - if (!fChain) return; - fChain->Show(entry); -} -Int_t hcana_class::Cut(Long64_t entry) -{ -// This function may be called from Loop. -// returns 1 if entry is accepted. -// returns -1 otherwise. - return 1; -} -#endif // #ifdef hcana_class_cxx diff --git a/examples/gtrk_e/make_class.C b/examples/gtrk_e/make_class.C deleted file mode 100644 index a22c5e9cad888e3f01b38a3f18f70ae384ee195f..0000000000000000000000000000000000000000 --- a/examples/gtrk_e/make_class.C +++ /dev/null @@ -1,17 +0,0 @@ -#include <TTree.h> -#include <TFile.h> -#include <TROOT.h> -#include <iostream> - -using namespace std; - -void make_class(string root_file="hodtest_52949.root", - string tree_name="T;1", - string class_name="hcana_class") { - - TFile* r_file = new TFile(root_file.c_str()); - // TTree *T = (TTree*)gROOT->FindObject("T;1"); - TTree *T = (TTree*)gROOT->FindObject(tree_name.c_str()); - T->MakeClass(class_name.c_str()); - -}