GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » [FIXED] use MC particle hypothesis in PndRecoKalmanTask
icon5.gif  [FIXED] use MC particle hypothesis in PndRecoKalmanTask [message #15617] Wed, 30 October 2013 14:17 Go to next message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
Dear all,

In the rec macro, when I enable the flag to use the MC particle hypothesis as
 recoKalmanFwd->SetIdealHyp(kTRUE);


my analysis macro crashes with the following error message.
Quote:


*************************************************************
initialisation for run id 1377867347
*************************************************************
-I- FairRunTimeDB::InitContainer() FairBaseParSet
-I- FairRunTimeDB::InitContainer() PndSensorNamePar
-I- FairRunTimeDB::InitContainer() PndMultiFieldPar
-I- FairRunTimeDB::InitContainer() ANAPidSelections
[ERROR ] init() ANAPidSelections not initialized
Error in <FairRuntimeDb::initContainers()>: Error occured during initialization
[INFO ] The number of entries in chain is 2000
[INFO ] Branch: PidAlgoIdealCharged not found in Tree
[INFO ] Branch: PidAlgoIdealCharged not found in Tree
Warning in <PndAnaPidCombiner::PndAnaPidCombiner::ReadTCA()>: No "PidAlgoIdealCharged" array found.

*** Break *** segmentation violation




Do any experts know how to fix this problem?
I have tried this flag with the other Kalman task "PndRecoKalmanTask* recoKalman", this problem doesn't happen. My PandaRoot is rev#20993. The 4-step simulation codes and the test analysis macro are attached.

Many thanks in advance.

Best,
Lu
  • Attachment: dspair.dec
    (Size: 0.33KB, Downloaded 147 times)
  • Attachment: sim.C
    (Size: 7.26KB, Downloaded 153 times)
  • Attachment: dig.C
    (Size: 4.59KB, Downloaded 165 times)
  • Attachment: rec.C
    (Size: 4.79KB, Downloaded 163 times)
  • Attachment: pid.C
    (Size: 3.58KB, Downloaded 163 times)
  • Attachment: mytask.C
    (Size: 8.99KB, Downloaded 177 times)

[Updated on: Thu, 31 October 2013 13:40] by Moderator

Report message to a moderator

Re: use MC particle hypothesis in PndRecoKalmanTask [message #15620 is a reply to message #15617] Wed, 30 October 2013 16:31 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
You need to add the mc track associator for the forward ideal tracker and to uncomment the correlated line in the kalman task:

  PndFtsTrackerIdeal* trackFts = new PndFtsTrackerIdeal();
  trackFts->SetRelativeMomentumSmearing(0.05);
  trackFts->SetVertexSmearing(0.05, 0.05, 0.05);
  trackFts->SetTrackingEfficiency(1.);
  trackFts->SetTrackOutput("FtsIdealTrack");
  fRun->AddTask(trackFts);


->PndMCTrackAssociator* trackMCfwd = new PndMCTrackAssociator();
->trackMCfwd->SetTrackInBranchName("FtsIdealTrack");
->trackMCfwd->SetTrackOutBranchName("FtsIdealTrackID");
->fRun->AddTask(trackMCfwd);

  PndRecoKalmanTask* recoKalmanFwd = new PndRecoKalmanTask();
  recoKalmanFwd->SetTrackInBranchName("FtsIdealTrack");
->recoKalmanFwd->SetTrackInIDBranchName("FtsIdealTrackID");
  recoKalmanFwd->SetTrackOutBranchName("FtsIdealGenTrack");
  recoKalmanFwd->SetBusyCut(50); // CHECK to be tuned
  recoKalmanFwd->SetIdealHyp(kTRUE);
  //recoKalmanFwd->SetNumIterations(3);
  fRun->AddTask(recoKalmanFwd);


Remember to switch on the idealhyp also for the barrel kalman task:

  recoKalman->SetIdealHyp(kTRUE);


The default macros in macro/run have such modification (but idealhyp set to kFALSE)
Re: use MC particle hypothesis in PndRecoKalmanTask [message #15627 is a reply to message #15620] Thu, 31 October 2013 12:46 Go to previous message
Lu Cao is currently offline  Lu Cao
Messages: 77
Registered: February 2013
continuous participant
From: *ikp.kfa-juelich.de
It works fine now.
Many thanks. Smile

Best,
Lu
Previous Topic: FIXED Access GetMicroCandidate from RhoCandList or RhoCandidate
Next Topic: loop over PndMcTracks
Goto Forum:
  


Current Time: Wed Sep 28 12:21:40 CEST 2022

Total time taken to generate the page: 0.00800 seconds