CERN根-根据时间筛选子树的事件时出错

CERN root - error when filtering events for subtree based on time

本文关键字：事件 出错 筛选 -根 时间 CERN      更新时间：2023-10-16

我之前在官方根目录（CERN）论坛上问过这个问题，但到目前为止这个问题仍未解决。也许这里的任何人都可以帮忙，要么指出我的错误，要么提出一种替代方法？

我有一个事件树；TTree有一个分支，其中包含每个事件的UNIX时间，还有一些其他分支。我想根据时间间隔选择事件的子集，这样我就可以分别分析这些事件。为了进行选择，我创建了另一个树，并在正确的时间间隔内复制所有条目。

以下代码工作得非常好，并将所有事件从树复制到子树：

void tree_time_filter(TTree* tree, TTree** subtree, Int_t time_i, Int_t time_f){
     *subtree = tree->CloneTree(0);
     Int_t t;
     tree->SetBranchAddress("UNIX time", &t);
     for(Long64_t i = 0; i<tree->GetEntries(); i++){
        tree->GetEntry(i);
        if (true)
        {
            (*subtree)->Fill();
        }
     }
     cout<<"Tree filtered. "<<(*subtree)->GetEntries()<< " entries were selected.n";
     return;
}

当我用实际条件替换if（true）时，问题发生了：

void tree_time_filter(TTree* tree, TTree** subtree, Int_t time_i, Int_t time_f){
     *subtree = tree->CloneTree(0);
     Int_t t;
     tree->SetBranchAddress("UNIX time", &t);
     for(Long64_t i = 0; i<tree->GetEntries(); i++){
        tree->GetEntry(i);
        if (t > time_i && t < time_f) //-> the condition
        {
            (*subtree)->Fill(); //-> this line now gives an error
        }
     }
     cout<<"Tree filtered. "<<(*subtree)->GetEntries()<< " entries were selected.n";
     return;
}

我收到错误："错误：指向类对象子树0x0 3084 c:/….的非法指针。*解释器错误已恢复*"

错误行引用了（*subtree）->Fill（），这与第一个示例中运行良好的代码相同。对于任何不涉及t的条件或任何不引用子树的if体，代码都有效。有人能解释一下这里出了什么问题吗？

谢谢！

（参考原始问题链接：http://root.cern.ch/phpBB3/posting.php?mode=edit&f＝3&p=79722）

我认为这是解释器的问题，在我看来是一个bug。

当宏第一次经过(*subtree)->Fill();时，它没有填充任何内容时，我可以重现您的错误。这很奇怪，因为它没有进入那里，但是。。。

我可以通过做两个循环来解决它：第一个循环循环，直到找到第一个通过剪切的事件，然后停止。第二个，(*subtree)->Fill();从这个事件开始，确保处理的第一个事件通过了剪切并被填充。

感谢Asen Christov，这里是该函数的工作版本。我把两个for循环放在while循环中，以防事件不按时间顺序排列；这应该有效，但没有经过广泛的测试。我真的不敢相信这是完成这样一项共同任务的最佳方式，所以任何"标准"的方式仍然受到欢迎。

void tree_time_filter(TTree* tree, TTree** subtree, Int_t time_i, Int_t time_f){
     *subtree = tree->CloneTree(0);
     Int_t t;
     tree->SetBranchAddress("UNIX time", &t);
     Long64_t i = 0;
     while (i < tree->GetEntries())
     {
         for(i; i<tree->GetEntries(); i++){
            tree->GetEntry(i);
            if (t > time_i && t < time_f) break;
         }
         for (i; i<tree->GetEntries(); i++){
            tree->GetEntry(i);
            (*subtree)->Fill();
            if (!(t > time_i && t < time_f)) break;
         }
     }
     cout<<"Tree filtered. "<<(*subtree)->GetEntries()<< " entries were selected.n";
     return;
}
相关文章：