GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » [FIXED] run_pid_sttcombi.C strange message
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 previous 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
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: problem running reco macro with the july12 release
Next Topic: run_sim_sttcombi_evtgen.C
Goto Forum:
  


Current Time: Mon Dec 02 22:08:47 CET 2024

Total time taken to generate the page: 0.00841 seconds