无法访问 3 维数组中的特定位置

Question

我正在使用ROOT库，并且有一个ROOT论坛，但我认为这不是ROOT的问题。

我有一个二维直方图 all_histos 的 3 维数组（直方图是 ROOT 中的对象），它们的大小取决于我希望将来能够更改的一些数组：

//this are the arrays
double pT_range[]         = {8, 4, 2, 1, 0.5}; //multiplicity ranges
double Aj_range[]         = {0.22, 0, 0};
double mult_range[]       = {234, 204, 188, 168, 152, 140, 128, 116, 104, 92, 76, 0}; //multiplicity ranges
//this is the first array of histograms
TH2D* all_histos[getSize(Aj_range)][getSize(pT_range)][getSize(mult_range)];

我必须初始化它们并用数据填充它们，那里没有问题。

然后我从中提取一些数据并将新数据保存在其他直方图数组 hist_final 中并保存它们，我可以看到它们是正确的。

TH1D* hist_final[getSize(Aj_range)][getSize(pT_range)][getSize(mult_range)];

  for (int a = 0; a < getSize(Aj_range); ++a)
  {
   for (int p = 0; p < getSize(pT_range); ++p)
    {
     for (int m = 0; m < getSize(mult_range); ++m)
     {
      for (int n = 0; n < 9; ++n)
       {
        all_histos[a][p][m]->GetXaxis()->SetRangeUser( 0.2*n, 0.2*(n+1) );//we first define each ring
        hist_final[a][p][m]->SetBinContent( n+1, all_histos[a][p][m]->GetMean(2) ); //then take the average in that ring
        hist_final[a][p][m]->SetBinError( n+1, all_histos[a][p][m]->GetMeanError(2) );
       }
      }
     }
    }

所有这些都有效。然后我想添加一些数据略有不同的新直方图，所以我决定将第二个坐标中的 hist_final 的维度增加 1：

TH1D* hist_final[getSize(Aj_range)][getSize(pT_range)+1][getSize(mult_range)];

然后我继续用稍微不同的数据填充这些直方图，只是为了得到“分段违规”。

for (int a = 0; a < getSize(Aj_range); ++a)
{
 for (int m = 0; m < getSize(mult_range); ++m)
 {
  for (int p = 1; p < getSize(pT_range); ++p)
  {
    all_histos[a][0][m]->Add(all_histos[a][p][m],1);
  }
  for (int n = 0; n < 9; ++n)
  {
    all_histos[a][0][m]->GetXaxis()->SetRangeUser( 0.2*n, 0.2*(n+1) );//we first define each ring
    hist_final[a][getSize(pT_range)][m]->SetBinContent( n+1, all_histos[a][0][m]->GetMean(2) ); //then take the average in that ring
    hist_final[a][getSize(pT_range)][m]->SetBinError( n+1, all_histos[a][0][m]->GetMeanError(2) );
  }
 }
}

我认为这可能是很多事情，我开始排除它们：

all_histos[a][0][m]->GetMean(2) 

是双重的，所以是

all_histos[a][0][m]->GetMeanError(2)

它们不是问题，然后我想到了垃圾箱，但如果我保存直方图而不填充它们，我可以打开它们看看它们有多少垃圾箱，他们有 9 个，这不是问题。

我最终发现它在

hist_final[a][getSize(pT_range)][m]

如果我改变了

getSize(pT_range)

对于 0 1 2 3 或 4 没有问题，尽管这会覆盖有用的数据。记得我用过：

TH1D* hist_final[getSize(Aj_range)][getSize(pT_range)+1][getSize(mult_range)];

要创建这个数组，它应该有 6 个位置，从 0 到 5。

TH1D* hist_final[getSize(Aj_range)][500][getSize(mult_range)];

谁在乎它有多大，我永远不会使用额外的空间，重要的是它可以工作。

但事实并非如此。它仍然在 5 和大于 5 的数字处失败。

我不知道为什么，但是如果我将第二个坐标更改为任何小于 5 的整数，则代码可以正常工作，我尝试手动输入坐标值的许多组合，只有第二个为 5 的那些组合会失败。

p>

正因为如此，我知道问题出在那个坐标上，但我这辈子不知道为什么，也不知道如何解决。

我的最后一个想法是指针可能是问题所在，我在集群中运行我的代码，所以我注销了应该删除我正在使用的所有临时内存，再次输入发送我的代码，然后运行它只是为了找到同样的问题。

我没有想法。谢谢。

你可以看到我的整个代码here。

Answer 1

分段错误通常是因为您试图访问内存中无法访问的位置。

在这行代码中：

 hist_final[a][getSize(pT_range)][m]

第二个字段应该是：

  hist_final[a][getSize(pT_range) - 1 ][m]

如果你想要hist_final[][] 中的最后一个元素（数组的第二维）。我不确定您想要哪个元素，但 pT_range 是 5，正如您所说，此字段中的有效输入为 0-4。

Answer 2

这个问题已经很老了，但也许我对此的回答可能对其他人有用。可能是您刚刚超出了 ROOT 的标准内存空间（您的直方图有多大？）。

据我了解，您只需创建诸如 TH2D* all_histos 之类的指针，但您不会像使用

那样为它们分配内存

TH2D**** all_histos = new TH2D***[getSize(Aj_range)];
 for (int a = 0; a < getSize(Aj_range); a++)
  { all_histos[a] = new TH2D***[getSize(pT_range)];
    for (int a = 0; a < getSize(Aj_range); p++)
       { all_histos[a][p] = new TH2D**[getSize(mult_range)];
       }
  }

事实上，ROOT 框架的预分配内存（在堆栈上）并不大（我认为是 64 MB）。如果您使用“new”命令分配内存，则可以超出此限制。 顺便说一句，我认为在某些时候您还会在索引中混淆 'getSize(Aj_range)' 和 'getSize(mult_range)'。