Accessing variables via SetBranchAddress

[air_shower]

Hi everyone,

I’m having problems accessing the variables for a given tree. Usually it worked for me like in the code example, but now it just says “unknown branch”:

    TFile *f = new TFile("selected.root");
    TTree *mytree=(TTree*)f->Get("recData");

    Float_t fEnergy;
    mytree->SetBranchAddress("fEnergy",&fEnergy);

    for(Int_t i = 0; i <mytree->GetEntries(); i++){
    mytree->GetEntry(i);
    cout << fEnergy<< endl;
    }

I think that I set the address incorrectly, although I tried variations like the full pathname (event.fFDEvents.fFdRecShower.fEnergy). Please see the attached screenshot of the tree-structure.
Thanks!

[Wile_E_Coyote]

In general, one would need your root file in order to try it, but as a first step you could try to use the “full name” of the branch which is something like “event.fFDEvents.fFdRecShower.fEnergy” (if it doesn’t work, either attach your root file here or put it somewhere in a “public area”, e.g. “afs”).

[air_shower]

Ok, here is the link:
dropbox.com/s/xevxtlb1cbzc6 … .root?dl=0

[Wile_E_Coyote]

See:

Read array from TChain

You’re trying to retrieve parts of a decomposed object (which was written with splitlevel >= 1). I’d say this kind of games is reserved for “advanced users” only, so you’re asking for trouble. You’ better have a look at an “analysis skeleton”: Header files for Trees from ROOT import * import sys, glob gROOT.ProcessLine( "struct var_t {\ Int_t Event_;\ Float_t Weight[1];\ };") files = ["/afs/cern.ch/user/s/shilpi/public/B-4p-0-1-v1510_14TEV_50PileUp_101175088.…

Problems with filling a histogram from Tree File

Your “photonbeam_E” branch is a variable length array of “Double_t” values (the actual length, “treeTagger” entry number specific, is given by the “nTagged” branch as an “Int_t” value), so it is not a simple “Double_t” variable.

Dictionary for class AttributeListLayout?

If your root file contains “unknown” classes, generate appropriate dictionary using: MyFile->MakeProject(“libMyFile”, “*”, “recreate++”); See also: http://root.cern.ch/root/html/TFile.html#TFile:MakeProject If this fails, you need to talk to those who created your root files.

Try:

{
  TFile *f = TFile::Open("selected.root");
  if (!f) return; // just a precaution
  
  TTree *t; f->GetObject("recData", t);
  if (!t) { delete f; return; } // just a precaution
  
  t->SetMakeClass(1);
  t->SetBranchStatus("*", 0);
  
  Int_t fFDEvents_;
  // note: you must know kMaxevent_fFDEvents from somewhere else,
  //       e.g. use TTree::MakeClass to get it
  const Int_t kMaxevent_fFDEvents = 3;
  Double_t fEnergy[kMaxevent_fFDEvents]; //[fFDEvents_]
  
  t->SetBranchStatus("event.fFDEvents", 1);
  t->SetBranchAddress("event.fFDEvents", &fFDEvents_);
  
  t->SetBranchStatus("event.fFDEvents.fFdRecShower.fEnergy", 1);
  t->SetBranchAddress("event.fFDEvents.fFdRecShower.fEnergy", fEnergy);
  
  for(Long64_t i = 0; i < t->GetEntries(); i++) {
    t->GetEntry(i);
    
    if (fFDEvents_ > kMaxevent_fFDEvents) {
      // note: the stack is corrupt, a SIGSEV crash is possible
      std::cout << i << " : " << fFDEvents_ << " :"
                << " FATAL ERROR : fFDEvents_ > kMaxevent_fFDEvents = "
                << kMaxevent_fFDEvents << std::endl;
      delete f; // automatically deletes "t", too
      return; // try to get out of here
    }
    
    std::cout << i << " : " << fFDEvents_ << " :";
    for (Int_t j = 0; j < fFDEvents_; j++) std::cout << " " << fEnergy[j];
    std::cout << std::endl;
  }
  
  // t->ResetBranchAddresses(); // detach "t" from local variables
  delete f; // automatically deletes "t", too
}

[air_shower]

Thank you very much!