GSI Forum: Analysis » [FIXED] use MC particle hypothesis in PndRecoKalmanTask

Home » PANDA » PandaRoot » Analysis » [FIXED] use MC particle hypothesis in PndRecoKalmanTask

Show: Today's Messages :: Polls :: Message Navigator

[FIXED] use MC particle hypothesis in PndRecoKalmanTask [message #15617]

Wed, 30 October 2013 14:17

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 224 times)
Attachment: sim.C
(Size: 7.26KB, Downloaded 223 times)
Attachment: dig.C
(Size: 4.59KB, Downloaded 245 times)
Attachment: rec.C
(Size: 4.79KB, Downloaded 247 times)
Attachment: pid.C
(Size: 3.58KB, Downloaded 236 times)
Attachment: mytask.C
(Size: 8.99KB, Downloaded 253 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

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: