GSI Forum: Tracking » GEM Tracking inside LHETRACK

Home » PANDA » PandaRoot » Tracking » GEM Tracking inside LHETRACK

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

GEM Tracking inside LHETRACK [message #8895]

Mon, 22 June 2009 22:20

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *ihep.ac.cn

Dear all,
under Ralf's suggestion I have extended the lhetrack code in order to be able to use even GEM hits.
The code is inside svn and it seems to work.

I have simulated 500 muons at 2GeV/c, from 5° to 20°, using tpc + mvd + 3 stations of gem detector.

I plot the momentum distribution for tracks reconstructed by lhetrack, with and without gems, and also after genfit.

index.php?t=getfile&id=5442&private=0

You can see the blue lines coming from MVD+TPC tacking, without GEMs.
Once I add even GEMs to the tracking, with only lhetrack fit (which means helix fit) I can see a peak but not at the correct position (dashed red line). Once I give my hits to genfit, the peak (continuous red line) goes to the proper position, and the improvement is quite evident with respect to the blu line (without GEMs).
In case of MVD+TPC+GEM genfit, the resolution at 2GeV/s seems something like 1.2%.

I would suppose that there is a lot of space for improvements (i.e. ideal tracking does not work yet for gems, and it is not so clear to me why I have an average GEM hit multiplicity at 7), and that the code should be checked for different momentum ranges and also for STT (volunteers are kindly requested), but I think this is a good starting point.

If you want to try, I have put inside tutorials/lhetrack the following macros:

run_sim_tpccombigem_pgun.C
run_digi_tpccombigem.C
run_reco_tpccombigem.C

(the standard run_kalman_tpc.C is working without any changes).

Feedback and comments are welcome.

Attachment: gem_tracking.gif
(Size: 8.01KB, Downloaded 643 times)

Report message to a moderator

Re: GEM Tracking inside LHETRACK [message #8921 is a reply to message #8895]

Fri, 26 June 2009 16:56

Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt

first-grade participant

From: *to.infn.it

Hi all.

As feedback/comment I reproduced Stefanos Plots with 1000 pions each at 1 GeV/c. Theta ranges from 20 to 140 Deg. Just to check if there is no big flaw.

You need to implement the gem library with gSystem->Load("libGem");.

Simulation:

  FairDetector *Gem = new PndGemDetector("GEM", kTRUE);
  Gem->SetGeometryFileName("gem_3Stations.root");
  fRun->AddModule(Gem);

Digi & Reco:

  // -----   GEM hit producers   ---------------------------------
  PndGemDigitize* gemDigitize = new PndGemDigitize("GEM Digitizer", iVerbose);
  fRun->AddTask(gemDigitize);
  
  PndGemFindHits* gemFindHits = new PndGemFindHits("GEM Hit Finder", iVerbose);
  fRun->AddTask(gemFindHits);

Find the parameters in pandaroot/macro/params/all.par

Add to LHETrack with:

  trackMS->SetGemMode(2);  // 0 OFF, 1 GEMPoint, 2 GEMHit     // GEMPoint smearing [cm], if negative no smearing

Here the plots with STT
index.php?t=getfile&id=5452&private=0

and TPC

I will update the Torino tutorial macros in pandaroot/tutorials/analysis.

Thanks to Stefano.

Cheers, Ralf.

Attachment: 1kpionsMvdTpcGem.png
(Size: 9.24KB, Downloaded 527 times)
Attachment: 1kpionsMvdSttGem.png
(Size: 10.18KB, Downloaded 561 times)

Report message to a moderator

Previous Topic:	GEM tracking
Next Topic:	new features in GENFIT

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Dec 04 22:10:55 CET 2024

Total time taken to generate the page: 0.00809 seconds