Ошибка в использовании Root / C ++: нет подходящего конструктора для инициализации 'TTree' - PullRequest
Новая
       49

Ошибка в использовании Root / C ++: нет подходящего конструктора для инициализации 'TTree'

0 голосов
/ 27 января 2020

Я пытаюсь построить инвариантную массу Z-бозона из двух распавшихся мюонов. Я использую MadGraph и Root. MadGraph моделирует события (pp> Z> mu + и mu-) и создает файл. root, который содержит события. Я вызвал генерацию событий моделирования.

Чтобы использовать Root для анализа данных и построения гистограммы, я должен написать свой код на C ++. Вот код: 

#ifdef __CLING__
R__LOAD_LIBRARY(libDelphes)
#include "TFile.h"
#include "TTree.h"
#include "TCanvas.h"
#include "TH1F.h"
#include <iostream>
#include "classes/DelphesClasses.h"
#include "external/ExRootAnalysis/ExRootTreeReader.h"
#include "external/ExRootAnalysis/ExRootResult.h"
#else

class ExRootTreeReader;
class ExRootResult;
using namespace std;
#endif

void readEvents(const char *inputFile)
{

 gSystem->Load("libDelphes");

 TChain chain("Delphes");
 chain.Add(inputFile);


 ExRootTreeReader *treeReader = new ExRootTreeReader(&chain);
 Long64_t numberOfEntries = treeReader->GetEntries();

 TTree *tree = new TTree(&chain);

 TClonesArray *branchMuon = treeReader->UseBranch("Muon");
 //Long64_t numberOfMuons = branchMuon->GetEntries();

 TClonesArray *branchElectron = treeReader->UseBranch("Electron");
 Long64_t numberOfElectrons = branchElectron->GetEntries();


cout << "There are " << numberOfEntries << " entries in your ntuple" << endl;

UInt_t  Muons;
Float_t MuonEta1;
Float_t MuonPhi1;
Float_t MuonPt1;
Float_t MuonEta2;
Float_t MuonPhi2;
Float_t MuonPt2;

tree->SetBranchAddress("MuonEta1", &MuonEta1);
tree->SetBranchAddress("MuonPhi1", &MuonPhi1);
tree->SetBranchAddress("MuonPt1", &MuonPt1);
tree->SetBranchAddress("MuonEta2", &MuonEta2);
tree->SetBranchAddress("MuonPhi2", &MuonPhi2);
tree->SetBranchAddress("MuonPt2", &MuonPt2);

TH1F *histDiMuonMass = new TH1F("mass", "M_{inv}(mu+[1], mu-[2]); M_inv (GeV/c^2); Events", 50, 0.0, 1500);
for(Int_t entry = 0; entry < numberOfEntries; ++entry)
{
 treeReader->ReadEntry(entry);

 Long64_t numberOfMuons = branchMuon->GetEntries();

 //Muons += entry;



 //Muon *muon1, *muon2, *muon3, *muon4;
 TLorentzVector muon1;
 TLorentzVector muon2;

 muon1.SetPtEtaPhiM(MuonPt1, MuonEta1, MuonPhi1, 0.1);
 muon2.SetPtEtaPhi(MuonPt2, MuonEta2, MuonPhi2, 0.1);


if(branchMuon->GetEntries() == 4)
{


 TLorentzVector Z = muon1 + muon2;
 histDiMuonMass->Fill(Z.M());
}
}
cout << "There are " << Muons << " entries in your ntuple" << endl;
histDiMuonMass->Draw();
}

Затем в терминале root я набираю: \ 

.x examples/readEvents.C("../eventgeneration/Events/run_01/tag_1_delphes_events.root")

Однако я получаю следующую ошибку: 

In file included from input_line_1563:1:
/mnt/c/1/MG5_aMC_v2_6_6/Delphes/examples/readEvents.C:30:20: error: no matching constructor for initialization of 'TTree'
 TTree *tree = new TTree(&chain);
                   ^     ~~~~~~
/home/cucip/builddir/include/TTree.h:295:4: note: candidate constructor not viable: no known conversion from 'TChain *' to 'const TTree' for 1st argument; remove &
   TTree(const TTree& tt) = delete;
   ^
/home/cucip/builddir/include/TTree.h:291:4: note: candidate constructor not viable: requires 0 arguments, but 1 was provided
   TTree();
   ^
/home/cucip/builddir/include/TTree.h:292:4: note: candidate constructor not viable: requires at least 2 arguments, but 1 was provided
   TTree(const char* name, const char* title, Int_t splitlevel = 99, TDirectory* dir = gDirectory);
   ^
In file included from input_line_1563:1:
/mnt/c/1/MG5_aMC_v2_6_6/Delphes/examples/readEvents.C:72:8: error: no member named 'SetPtEtaPhi' in 'TLorentzVector'
 muon2.SetPtEtaPhi(MuonPt2, MuonEta2, MuonPhi2, 0.1);

Кто-нибудь видит, что может вызвать эту ошибку и как я могу ее исправить?

Большое спасибо!

1 Ответ

1 голос
/ 27 января 2020

Вызывая new TTree(&chain), вы пытаетесь вызвать конструктор для TTree, и вы предоставляете ему TChain * (тип chain равен TChain, поэтому &chain является указателем на a TChain).

Конструкторы для TTree https://root.cern/doc/v618/classTTree.html не принимают ни аргумента, ни двух const char* (обычно строковые литералы или .Data() из TString или .c_str() из std::string или вывод ROOT s Form…).

Похоже, вы пытаетесь вызвать удаленный конструктор TTree(const TTree&tt)=delete, который не работает для нескольких причины:

  • вы предоставляете указатель на TTree вместо TTree по ссылке (даже не уверен, с макушки головы, если приведение с TChain на TTree будет работать, если & не было.
  • конструктор удален, поэтому он на самом деле не существует и не может быть вызван.

Однако, судя по виду это значит, что вы хотите использовать TChain в качестве TTree*, для начала не требующий кода с вашей стороны, поскольку TChain уже является TTree по наследству. n используйте chain. вместо tree-> или (если переменная указателя действительно более желательна по причинам, которые я не вижу в этом примере) создайте TTree* (т.е. создайте только указатель на TTree вместо создания фактического TTree в куче и получения указателя на него от оператора new), сделав что-то вроде 

TTree* tree = &chain;
...