c++ - 无法访问 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) 

是一个double，也是

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

它们不是问题，然后我想到了 bin，但如果我保存直方图而不填充它们，我可以打开它们并查看它们有多少个 bin，它们有 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 的组合才会失败。

正因为如此，我知道问题出在那个坐标上，但我终究无法弄清楚为什么，也不知道如何解决。

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

我没主意了。谢谢。

你可以看到我的全部代码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。