| [FIXED] run_pid_sttcombi.C strange message [message #13876] |
Wed, 22 August 2012 11:43  |
Elisa Fioravanti
Messages: 84
Registered: January 2008
|
continuous participant |
From: *fe.infn.it
|
|
Dear all,
I'm working with july12 version.
I'm into the directory pandaroot/macro/pid
and I'm running the macros:
run_sim_sttcombi_evtgen.C
run_digi_sttcombi.C
run_reco_sttcombi.C
run_pid_stt.C
loading the file that you can read in macro/run/tdrct/npipi/pipi.dec
in order to recontruct two pions decay channel.
Everything is fine, no crashes, but when I run the macro run_pid_stt.C I have the message that you can see below.
Is it normal?
When I tried to run my analysis macro macro/run/tdrct/npipi/run_ana_invariantmass_2pi_stt.C I had empty histograms.
This is the reasons why I'm asking if this message that I got running the pid macro is normal or not.
Thanks in advance to everybody.
Elisa
Cross section calculation concluded successfully
I- FairGeane::FairGeane: Geane is Initialized
[INFO ] Branch: FtsIdealGenTrack not found in Tree
[INFO ] Branch: FtsIdealGenTrack not found in Tree
-I- PndPidCorrelator::Init: No 2nd PndTrack array!
[INFO ] Branch: PidChargedCand not found in Tree
[INFO ] Branch: PidChargedCand not found in Tree
-I- PndPidIdealAssociatorTask::Init: No PndPidCandidate array PidChargedCand there!
InitStatus PndPidMvdAssociatorTask::Init()
[INFO ] Branch: PidChargedCand not found in Tree
[INFO ] Branch: PidChargedCand not found in Tree
-I- PndPidMvdAssociatorTask::Init: No PndPidCandidate array PidChargedCand there!
InitStatus PndPidMdtHCAssociatorTask::Init()
[INFO ] Branch: PidChargedCand not found in Tree
[INFO ] Branch: PidChargedCand not found in Tree
-I- PndPidMdtHCAssociatorTask::Init: No PndPidCandidate array PidChargedCand there!
InitStatus PndPidDrcAssociatorTask::Init()
[INFO ] Branch: PidChargedCand not found in Tree
[INFO ] Branch: PidChargedCand not found in Tree
-I- PndPidDrcAssociatorTask::Init: No PndPidCandidate array PidChargedCand there!
InitStatus PndPidDiscAssociatorTask::Init()
[INFO ] Branch: PidChargedCand not found in Tree
[INFO ] Branch: PidChargedCand not found in Tree
-I- PndPidDiscAssociatorTask::Init: No PndPidCandidate array PidChargedCand there!
[INFO ] Branch: PidChargedCand not found in Tree
[INFO ] Branch: PidChargedCand not found in Tree
-E- PndPidSttAssociatorTask::Init: No PidChargedCand array!
InitStatus PndPidEmcBayesAssociatorTask::Init()
[INFO ] Branch: PidChargedCand not found in Tree
[INFO ] Branch: PidChargedCand not found in Tree
-I- PndPidEmcBayesAssociatorTask::Init: No PndPidCandidate array PidChargedCand there!
[INFO ] The number of entries in chain is 1000
[INFO ] The number of entries in chain is 1000
......
Macro finished succesfully.
Output file is pid_sttcombi.root
Parameter file is params_sttcombi.root
Real time 13.2037 s, CPU time 12.6 s
[Updated on: Tue, 16 October 2012 14:50] by Moderator Report message to a moderator |
|
|
|
| Re: run_pid_sttcombi.C strange message [message #13878 is a reply to message #13876] |
Thu, 23 August 2012 16:49  |
Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
|
first-grade participant |
From: *gsi.de
|
|
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
