Dear Elisa,
I just saw your posting and try to reproduce your problem. For me, the errors don't look healthy, since obviously no charged track candidates are found.
What may be origin of the problem is, that your analysis macro uses a quite outdated kind of data access.
Could you try to replace in your run_ana_... macro this block
TString inPidFile = "evt_pid_stt.root";
TString inSimFile = "evt_points_stt.root";
TFile *inFile = TFile::Open(inSimFile,"READ");
TTree *tree=(TTree *) inFile->Get("cbmsim") tree->AddFriend("cbmsim",inPidFile);
TClonesArray* mc_array=new TClonesArray("PndMCTrack");
tree->SetBranchAddress("MCTrack",&mc_array);
TClonesArray* cand_array=new TClonesArray("PndPidCandidate");
tree->SetBranchAddress("PidChargedCand",&cand_array);
FairMCEventHeader* evthead;
tree->SetBranchAddress("MCEventHeader.", &evthead);
TFile *out=TFile::Open("invariantmass_2pi_stt.root","RECREATE");
PndEventReader evr(inPidFile);
with that block
FairLogger::GetLogger()->SetLogToFile(kFALSE);
FairRunAna* fRun = new FairRunAna();
FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
fRun->SetInputFile("evt_points_stt.root");
fRun->AddFriend("evt_pid_stt.root");
fRun->SetOutputFile("output.root");
fRun->Init();
TFile *out=TFile::Open("invariantmass_2pi_stt.root","RECREATE");
PndAnalysis evr;
You also need to comment the line
//evthead->GetVertex(mcVertex);
for the time being in that case.
I'm also going to try myself, but it will take some time to generate and simulate some events.
Best regards,
Klaus
GSI Forum
Search
Help
Members
Login
Home
