GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Tracking » Genfit initial detector plane
icon4.gif  Genfit initial detector plane [message #12550] Mon, 22 August 2011 16:14 Go to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *17-79-r.retail.telecomitalia.it
Dear genfit experts,
I have one big doubt about the correct way to run the kalman inside the PndRecoKalmanTask.

As you know, the standard PR algorithms are providing a PndTrack object, which should be converted into a GFTrack in order to run the Kalman. In our case we are use GeaneTrackRep representation.
The way I am using to convert the track
(GenfitTools/recotasks/PndRecoKalmanFit.cxx) is the following:


  GFAbsTrackRep* rep = 0;
  GFDetPlane start_pl(StartPos, TVector3(1.,0.,0.), TVector3(0.,1.,0.));
  GeaneTrackRep *grep = new GeaneTrackRep(fPro,
                                          start_pl,StartMom,
                                          StartPosErr,StartMomErr,
                                          fCharge,PDGCode);
  grep->setPropDir(1);
  rep = grep;
  
  GFTrack* trk= new GFTrack(rep);
  PndTrackCand trackCand = tBefore->GetTrackCand();
  trk->setCandidate(*PndTrackCand2GenfitTrackCand(&trackCand));



My main concern is about the starting point to use, i.e. start_pl StartMom StartPosErr StartMomErr. I can see two possible approaches:

a) The start detector plane is taken at the IP after a backward extrapolation (in this case XY)
b) The start detector plane is at the first point of the track

In the first case we have a back extrapolation, and maybe errors are becoming larger due to the geane propagation... not sure what is happening for secondaries.
In the second case, I am wondering if genfit uses the first hit as fit point or skip it, in theory we should have a point less, am I right? In particular, if the first hit is "inside" the first TpcCluster/SttHit, I am wondering how could the first extrapolation be done...

Another point is about the orientation of the detector plane. In general, by default, I am using the PCA to 0,0,0 and XY as plane..but I would like to understand if maybe the more correct way could to a plane perpendicular to the track direction.. or maybe this does not change the final results (or maybe this does change the way to calculate initial errors).

Could you please clarify which should be the most correct way?

Regards
Re: Genfit initial detector plane [message #12551 is a reply to message #12550] Mon, 22 August 2011 19:56 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
From: *natpool.mwn.de
Hi Stefano,

you should initialize the TrackRep at the position of the first hit.
As the Kalman makes some iterations, the first hit is then properly taken into account for the fit.


regards,

Johannes
Re: Genfit initial detector plane [message #12552 is a reply to message #12551] Mon, 22 August 2011 20:07 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *17-79-r.retail.telecomitalia.it
Ok thanks,I will try.
What about the direction of the starting plane?
Re: Genfit initial detector plane [message #12553 is a reply to message #12552] Mon, 22 August 2011 20:41 Go to previous message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
From: *natpool.mwn.de
The direction should point into the direction of the track at the first point.
Previous Topic: new DPM runs started
Next Topic: Event Mixing with TPC
Goto Forum:
  


Current Time: Wed Apr 24 03:00:37 CEST 2024

Total time taken to generate the page: 0.00949 seconds