为什么这个循环不会加载正确的事件?

时间:2016-11-16 14:15:16

标签: c++ loops root-framework

这是我写的代码,它运行正常,但它似乎并没有实际加载事件......相反,它似乎用一个任意的微小数字填充直方图和cout。

include "TTree.h"
include "TFile.h"
include "TRandom.h"
include "TTree.h"
include "TCanvas.h"
include "TStyle.h"
include <iostream>
include "TH1.h"
include "TChain.h"

using namespace std;

void readtestsimple(){

  TFile *file = TFile::Open("/file/");

  TTree* comp_tree = (TTree*)file->Get("_compress_tree");

  TChain* compchain = new TChain("_compress_tree");
  compchain->Add("/file/");

  TH1D *comp = new TH1D("comp","_compression",100,0,0.05);

  Double_t _compression, _compressionU, _compressionV, _compressionY, _ch_compression;
  Int_t _ch;
  comp_tree->SetBranchAddress("_compression",&_compression);
  comp_tree->SetBranchAddress("_compressionU",&_compressionU);
  comp_tree->SetBranchAddress("_compressionV",&_compressionV);
  comp_tree->SetBranchAddress("_compressionY",&_compressionY);

      Long64_t nentries = compchain->GetEntries();
      Long64_t nb = 0, nbytes = 0 ;

      for (Long64_t i=0; i<nentries; i++) {
              Long64_t entry = compchain->LoadTree(i);
              if (entry <0 ) break;

             nb = compchain->GetEntry(i); nbytes += nb;

              cout << "entry = " << nb << endl;
              cout << "compression = " << _compression << endl;

              comp->Fill(_compression);
      }

   cout << "entries = " << nentries << endl;

}

我得到的输出是长期重复:

entry = 36
compression = 1.54768e-316

entry = 36
compression = 1.54768e-316

entry = 36
compression = 1.54768e-316

(最后,正确的条目数)。

如何让它读出正确的条目? 请帮忙!

1 个答案:

答案 0 :(得分:0)

看起来你根本不需要使用comp_tree指针。 compchain指针应该足以访问TTree中的数据,因为TChainTTree。因此,您应该能够将变量“链接”到链/树分支:

compchain->SetBranchAddress("_compression",&_compression);

而不是

comp_tree->SetBranchAddress("_compression",&_compression);