当您将ROOT对象作为参数传递给C++中的函数时，您可以访问该对象中的方法吗

已经指出，不能返回局部变量的地址，当函数结束时，该地址将被销毁。

你的另一个问题是：

histo->GetNbinsX()

不起作用。我试着在主脚本和脚本的子例程中调用它：它非常适合我，使用当前的ROOT版本。在这个问题中，您拼错了GetNBinsX(是的，这将与骆驼式策略更加一致(。可能

无论如何，我相信你会很高兴知道ROOT有一个非常智能的1D峰值搜索算法可以使用：寻找TSspectrum类。

您的代码存在各种问题。

最引人注目的是这个：

int *peak_counter1d(TH1F *histogram,int peak_num,int threshold = 5,int display = 0)
{
    //...
    double peak_bins[peak_num];
    //...
    return peak_bins;
}

您正在返回一个指向局部变量的指针。返回指向局部变量的指针是未定义的行为。

下一个问题是：

  int *check_array; // to put previously found peak bins
  while(counter < peak_num)
  {
      //...
      check_array = new int[counter];
  }
  delete[] check_array;

您有潜在的内存泄漏，因为您正在循环而从未释放check_array。另一个问题是，如果该循环从未执行，那么您正在对未初始化的变量调用delete []。

下一个问题是：

int * peak_counter1d(...)
{
    double peak_bins[peak_num];
    //...
    return peak_bins;
}

即使可以安全地返回指向局部变量的指针，函数也会返回int*，但返回的是double *。

下一个问题是：

TCanvas *look = new TCanvas("look","look",500,400);

您正在分配look，但从未解除分配，甚至从未使用过它

您也可以在main:中执行相同的操作

 TH1F *histo = new TH1F("histo","try",100,0,10);
  TF1 *f1 = new TF1("f1","exp(-x/10)*sin(x)*sin(x)",0,10);

C++不是Java。您不必使用new创建对象。

TH1F histo("histo","try",100,0,10);
  TF1 f1("f1","exp(-x/10)*sin(x)*sin(x)",0,10);

除了最后一个问题外，如果您使用std::vector而不使用new[]来创建动态阵列，那么总体而言，这些问题都可以得到解决。

应用这些更改，代码应该看起来像这样(未编译(：

#include <iostream>
#include <algorithm>
#include <iterator>
#include <vector>
#include "TROOT.h"
#include "TCanvas.h"
#include "TH1.h"
#include "TF1.h"
using namespace std;
vector<int> peak_counter1d(TH1F *histogram,int peak_num,int threshold = 5,int display = 0)
{
  if(display == 1)
  {
    // TCanvas *look = new TCanvas("look","look",500,400);
     histogram->Draw();
  }
  int total_bins = histogram->GetNBinsX();
  vector<int> peak_bins(peak_num);
  peak_bins[0] = histogram->GetMaximumBin();
  int counter = 1;
  vector<int> check_array; // to put previously found peak bins
  while(counter < peak_num){
    double peak = threshold;
    double peak_loc = -500;
    check_array.resize(counter);
    for(int i=0; i<counter; i++){
      check_array[i] = peak_bins[i]; // fills the array with previously found peak bins
    }
    for(int i=0; i<total_bins; i++){
      if(peak < histogram->GetBinContent(i)){ 
        bool exists = find(begin(check_array),end(checkarray),i); // makes sure this is a peak we haven't already found
        if(!exists){
          peak = histogram->GetBinContent(i);
          peak_loc = i;
        }
      }
    }
    peak_bins[counter] = peak_loc;
    counter ++;
  }
  return peak_bins;
}
void testing(){
  gROOT->Reset();
  TH1F histo("histo","try",100,0,10);
  TF1 f1("f1","exp(-x/10)*sin(x)*sin(x)",0,10);
  double val;
  for(int i=0; i<100; i++){
    val = f1.Eval(i/10.0);
    //cout << i << "t" << i/100.0 << "t" << val << endl;
    histo.SetBinContent(i,val);
  }
  vector<int> peak_bins = peak_counter1d(&histo,3,5,1);
  for(int i=0; i<3; i++){
    cout << i << "t" << peak_bins[i] << endl;
  }
}