GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » [FIXED] run_pid_sttcombi.C strange message
[FIXED] run_pid_sttcombi.C strange message [message #13876] Wed, 22 August 2012 11:43 Go to next message
Elisa Fioravanti is currently offline  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 Go to previous messageGo to next message
Klaus Götzen is currently offline  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
Re: run_pid_sttcombi.C strange message [message #13881 is a reply to message #13876] Fri, 24 August 2012 09:38 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *gsi.de
Dear Elisa,


I just run the simulation with the macros in macro/pid and the analysis with run_ana_invariantmass_2pi_stt.C on this data and got non-empty histograms. Also the errors you reported for run_pid_stt.C did not appear. Of course I don't know whether the plots are reasonable, but at least not empty.

Are you sure that the adapted all file names in the macro properly and did an 'svn update' on your july12 branch (maybe something changed, although it should be fixed in principle...)?


Best regards,
Klaus
Re: run_pid_sttcombi.C strange message [message #13882 is a reply to message #13876] Fri, 24 August 2012 09:42 Go to previous messageGo to next message
Elisa Fioravanti is currently offline  Elisa Fioravanti
Messages: 84
Registered: January 2008
continuous participant
From: *fe.infn.it
Dear Klaus,

first of all that you for helping me.

I'm running just now the pid macro again, and the error in this macro did not appear anymore. It is strange but it is good.

I'm modifying also my analysis macro, according with your suggestions.

I let you know if the problem is solved or not.

Thanks again
Elisa
Re: run_pid_sttcombi.C strange message [message #13883 is a reply to message #13876] Fri, 24 August 2012 11:09 Go to previous message
Elisa Fioravanti is currently offline  Elisa Fioravanti
Messages: 84
Registered: January 2008
continuous participant
From: *fe.infn.it
Dear Klaus,

now everything works,
and my histograms are not empty!

Thanks a lot for helping me

Elisa
Previous Topic: problem running reco macro with the july12 release
Next Topic: run_sim_sttcombi_evtgen.C
Goto Forum:
  


Current Time: Thu Mar 28 12:15:06 CET 2024

Total time taken to generate the page: 0.00987 seconds