4

简短的问题： 如果您必须有许多（甚至超过 200 个）成员变量（它们中的每一个都将是我们分析中感兴趣的一些物理量的图）。声明这些变量的最佳位置是什么？

长解释： 我用于分析的框架创建了一个循环事件的类，它可以缩小到

constructor()
initialize()
execute()
finalize()

在标题中，您可以声明如下指针（这显然是ROOT我们必须使用的包的要求）：

std::vector<double> *m_jet_pt;

和指向直方图类的指针：

TH1F *h_jet_pt;

然后constructor必须初始化指向某个明确内存地址的指针（据我所知，这是为了以后从文件中读取数据）

constructor()
{
  // this is data, will be associated to a TTree later
  m_jet_pt = 0;

  // this is a histogram
  h_jet_pt = new TH1F("name", "title", nbins, min_bin, max_bin); 
}

然后在initialize函数中打开一个包含先前存储的数据的文件，并将类成员的地址设置为指向文件中包含的对象（我不知道这是如何工作的，或者我的陈述是否准确）：

TFile *file = new TFile("filename.root");

// The tree is where the data is (for some reason ROOT uses C style casts)
TTree *tree = (TTree*)file->Get("MyTree");

// this is how to set the address of the member class pointer
// to point to the data in the tree:
tree->SetBranchAddress("m_jet_pt", &jet_pt);

然后在execute循环中调用函数，每个事件一次，你会做你想要的物理选择

execute()
{
  // I omit here where the `i` index comes from as not relevant
  if(m_jet_pt->at(i) > 30)

    h_jet_pt->Fill();
}

最后在finalize函数中，您将在选择完成后的每个事件中执行所有您需要做的事情，例如将直方图存储在另一个文件中

finalize()
{
  h_jet_pt->Write();
}

现在想象一下，我有 20 个 quantites jets，30 个electrons，30 个muons等等，你可以看到成员变量的数量是如何变得巨大的！很快代码就变得一团糟，那么您的专业程序员将如何处理这种情况？希望这足够清楚！

4

1 回答 1

0

ROOT 具有相当不错的 I/O 功能，并支持从/向 TTree 读取和写入任意数据结构。通常，您需要做的就是编写一个具有与您的单个“测量”或“事件”相对应的数据结构的类。如何在此类中组织多个物理量取决于您。该结构可以是平面结构，具有 200 个同一级别的数据成员，也可以包含封装的其他类类型的非 POD 数据成员。好消息是 ROOT 可以很好地与所有不同的类型一起工作，包括标量、数组、指针和 stl 容器作为数据成员。

创建类的对象后，它可以与 TTree 上的顶级分支相关联，如下所示：

tree->Branch(branchname, className, &p_object, bufsize, splitlevel)

有关详细信息，请参见案例 B：https ://root.cern.ch/doc/master/classTTree.html

如果您为所有数据结构生成字典和流媒体，ROOT 将自动创建子分支，从而为您节省大量手动工作。之后，您可以使用 TTree 的方法，例如 Fill() 和 GetEntry()，通过将对象保存在树中或从树中读取值到对象来创建条目。

于 2016-06-17T05:14:25.967 回答