GSI Helmholtzzentrum für Schwerionenforschung
[FIXED] run_pid_sttcombi.C strange message
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
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
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
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
Dear Klaus, now everything works, and my histograms are not empty! Thanks a lot for helping me Elisa