Мне нужно суммировать тысячи гистограмм из одного каталога - PullRequest
Новая
       39

Мне нужно суммировать тысячи гистограмм из одного каталога

0 голосов
/ 02 июля 2019

У меня есть каталог Processed_Data с тысячами hists*****_blinded.root файлов.Каждый hists*****_blinded.root содержит около 15 графиков и гистограмм.Моя цель - просто перекрыть 1 конкретную гистограмму sc***** из каждого файла, чтобы получить окончательную гистограмму finalhists_blinded.root, которая будет представлять все перекрывающиеся вместе.

Я пробовал следующий макрос: 

void final()
{
   TCanvas *time = new TCanvas("c1","overlap" ,600,1000);
   time ->Divide(1,1);
   time ->cd(1);

   TH1F *h1 = new TH1F("h1","time" ,4096,0,4096);
   ifstream in;
   Float_t t;
   Int_t nlines= 0;
   in.open("Processed_Data", ios::in);
   while (1) {
      in >> t;
      if (!in.good()) break;
      h1->Fill(t);
      nlines++;
   }
   in.close();

Но я получаю чистый холст в конце.Идея состоит в том, чтобы пропустить каждый файл hists через код и добавить каждый по одному.

В результате я хочу видеть, что все эти гистограммы sc***** перекрываются, так что пики в каждой из них будут создавать шаблон в файле finalhists_blinded.root.

1 Ответ

1 голос
/ 03 июля 2019

Не должно быть так сложно, попробуйте это: 

void overlap()
{
        TCanvas *time = new TCanvas("c1", "overlap", 0, 0, 800, 600);

        const char* histoname = "sc";

        const int NFiles = 100000;
        for (int fileNumber = 0; fileNumber < NFiles; fileNumber++)
        {
                TFile* myFile = TFile::Open(Form("Processed_Data/hists%i_blinded.root", fileNumber));
                if (!myFile)
                {
                        printf("Nope, no such file!\n");
                        return;
                }
                TH1* h1 = (TH1*)myFile->Get(histoname);
                if (!h1)
                {
                        printf("Nope, no such histogram!\n");
                        return;
                }
                h1->SetDirectory(gROOT);
                h1->Draw("same");
                myFile->Close();
        }
}

Он зацикливается на всех Processed_Data/histsXXXXXi_blinded.root файлах (учитывая их имена Processed_Data/hists0_blinded.root, Processed_Data/hists1_blinded.root, Processed_Data/hists2_blinded.root, ..., Processed_Data/hists99998_blinded.root, Processed_Data/hists99999_blinded.root), открывает каждый из них, получает гистограмму 1D sc, добавляет ее на холст, закрывает файл и переходит к следующему файлу.

...