ROOT(cern) чтение файла +TGraph работает нормально, но TGraphErrors не хочет

Question

Проблема заключается в следующем: я прочитал мои данные из файла в массив. Построение графика работает без проблем. Но TGraphErrors не хочет ничего делать. Рут говорит

Ошибка: Невозможно вызвать TGraphErrors::SetPoint(i, Данные [0][i], Данные [1][i], Данные [2][i], Данные [2][i]) в текущем графике области действия2.C:85: Возможные кандидаты... (в TGraphErrors) (в TGraph) /usr/local/Cellar/root/5.34.26/lib/root/libHist.so -1: -1 0 public: virtual void TGraph:: SetPoint (Int_t i, Double_t x, Double_t y); * Устранена ошибка интерпретатора *"

Код ниже:

    void plot2(char* fin1){
    const int VAR_N=20,NPOINTS=10000;
    Double_t Data[VAR_N][NPOINTS];
   int counts,N1;
   float E,E1,dE;
   FILE* fin12;
    printf("%s \n",fin1);
    fin12 = fopen(fin1,"r");
      if (fin12 == NULL){
        printf("Can't open input file: %s !\n",fin1);
    }
    else{
        counts=0; 
        while (!feof(fin12)) {
            fscanf(fin12, "%f %f %f \n",&E,&E1,&dE);
            printf("%f %f %f \n",E,E1,dE);
                Data[0][counts] =   E;
                Data[1][counts] =   E1;
                Data[2][counts] =   dE;  }
        N1=counts;
        printf("NP1=%d ",N1);
        fclose(fin12);
    }
   TGraph* V_Graph = new TGraph(N1);
    TGraphErrors* V_GraphErrors = new TGraphErrors(N1);
    for(int i=0;i<N1;i++){
       V_Graph->SetPoint(i,Data[0][i],Data[1][i]);
     } 
   for(int i=0;i<N1;i++){
   V_GraphErrors->SetPoint(i,Data[0][i],Data[1][i],Data[2][i],Data[2][i]);
  }

TCanvas *c1 = new TCanvas("c1", "c1", 1200, 1000);
  c1->cd(1);
   TString VName="+Errors";
   V_Graph->SetTitle(VName);
   V_Graph->SetMarkerStyle(21);
   V_Graph->SetMarkerSize(1);
   V_Graph->SetMarkerColor(kBlue);
    V_Graph->GetXaxis()->SetTitle("x");
    V_Graph->GetYaxis()->SetTitle("y");
    V_Graph->GetXaxis()->SetLimits(4.,25.);
    V_Graph->GetYaxis()->SetRangeUser(10.,80.);
    V_Graph->Draw("APL");

    V_GraphErrors->SetMarkerStyle(21);
    V_GraphErrors->SetMarkerSize(1);
    V_GraphErrors->SetMarkerColor(kGreen);
    //V_Graph->GetXaxis()->SetLimits(4.,11.);
    //V_Graph->GetYaxis()->SetRangeUser(45.,80.);
    V_GraphErrors->Draw("C*");
}

Я новичок в использовании root. Вероятно, я сделал что-то глупое и не заметил этого. Так что я сделал не так? заранее спасибо
Olka

Answer 1

Проблема в том, что нет SetPointError метод с подписью, которую вы пытаетесь использовать (int, Double_t, Double_t, Double_t, Double_t), Чтобы установить ошибки, вы должны использовать SetPointError(Double_t ex, Double_t ey) метод. Например, в приведенном выше коде вы должны сделать что-то вроде:

V_GraphErrors->SetPoint(i,Data[0][i],Data[1][i]]);
V_GraphErrors->SetPointError(i,Data[2][i],Data[2][i]);

Хотя, вероятно, было бы легче передать arrys непосредственно в конструкторе:

V_GraphErrors(N1, Data[0], Data[1], Data[2], Data[2]);

Кстати, это установит одинаковые ошибки x и y. Я не уверен, что это то, что вы хотели.

См. https://root.cern.ch/root/html/TGraphErrors.html для получения подробной информации о методах, доступных в TGraphErrors