search

klendathu2k / stroot-forward-tracker

STAR Forward Tracker Integration
GNU General Public License v3.0
0 stars 0 forks source link

We are not seeing sTGC hits in the code :: issue w/ g2t_stg.F #2

Closed klendathu2k closed 5 years ago

klendathu2k commented 5 years ago

St_geant_Maker notes the presence of both inner and outer silicon hits, along with the sTCG hits. So not clear why the tracker is not finding them.

StgMaker's conversion of the g2t hits StgMaker:INFO - track_id=10 volume_id=6 x/y/z -5.62433/-15.9559/167.193 StgMaker:INFO - track_id=7 volume_id=4 x/y/z 7.7981/5.78583/141.129 StgMaker:INFO - track_id=2 volume_id=4 x/y/z -0.228333/6.43634/141.129 StgMaker:INFO - track_id=9 volume_id=4 x/y/z -9.71069/0.69635/141.129 StgMaker:INFO - track_id=6 volume_id=4 x/y/z -10.2369/-7.08618/141.129 StgMaker:INFO - track_id=5 volume_id=4 x/y/z 3.3512/-7.02087/141.129 StgMaker:INFO - track_id=8 volume_id=4 x/y/z 8.04059/-1.65814/141.129 StgMaker:INFO - track_id=4 volume_id=4 x/y/z 11.9695/-4.95349/141.129 StgMaker:INFO - track_id=1 volume_id=4 x/y/z 10.156/-4.25897/141.129 StgMaker:INFO - track_id=7 volume_id=5 x/y/z 8.5708/6.36053/155.129 StgMaker:INFO - track_id=2 volume_id=5 x/y/z -0.252869/7.0755/155.129 StgMaker:INFO - track_id=9 volume_id=5 x/y/z -10.6666/0.743896/155.129 StgMaker:INFO - track_id=6 volume_id=5 x/y/z -11.2439/-7.80023/155.129 StgMaker:INFO - track_id=5 volume_id=5 x/y/z 3.6853/-7.71704/155.129 StgMaker:INFO - track_id=8 volume_id=5 x/y/z 8.83862/-1.82129/155.129 StgMaker:INFO - track_id=4 volume_id=5 x/y/z 13.1585/-5.44147/155.129 StgMaker:INFO - track_id=1 volume_id=5 x/y/z 11.1636/-4.67896/155.129 StgMaker:INFO - track_id=7 volume_id=6 x/y/z 9.34338/6.9353/169.129 StgMaker:INFO - track_id=2 volume_id=6 x/y/z -0.27771/7.71497/169.129 StgMaker:INFO - track_id=9 volume_id=6 x/y/z -11.6188/0.780151/169.129 StgMaker:INFO - track_id=6 volume_id=6 x/y/z -12.2511/-8.51862/169.129 StgMaker:INFO - track_id=5 volume_id=6 x/y/z 4.01947/-8.41296/169.129 StgMaker:INFO - track_id=8 volume_id=6 x/y/z 9.6369/-1.98425/169.129 StgMaker:INFO - track_id=4 volume_id=6 x/y/z 14.3481/-5.92877/169.129 StgMaker:INFO - track_id=1 volume_id=6 x/y/z 12.1709/-5.09869/169.129

St_geant_Maker's end-of-run hit totals per set / detector QA :INFO - St_geant_Maker summary -QA- g2t hits (BBCH:BPOL) nhits = 27 QA :INFO - St_geant_Maker summary -QA- g2t hits (EPDH:EPDT) nhits = 117 QA :INFO - St_geant_Maker summary -QA- g2t hits (FSTH:FTIS) nhits = 140 QA :INFO - St_geant_Maker summary -QA- g2t hits (FSTH:FTOS) nhits = 50 QA :INFO - St_geant_Maker summary -QA- g2t hits (HCAH:HSCI) nhits = 200 QA :INFO - St_geant_Maker summary -QA- g2t hits (STGH:TGCG) nhits = 352 QA :INFO - St_geant_Maker summary -QA- g2t hits (WCAH:WSCI) nhits = 155

... looking at the track table ... int n_stg_hit 0 : 0 : 0 : 0 : 0 : 0 : Nhits on sTGC [new Forward small thin gap chambers]

no hits on the track from the sTGC. Strange. And, this happens when inspecting the hits using fz2root.C. So...

In St_geant_Maker::AddHits(...) -- the hit total is coming from the Gfnhit geant call. Then the table is filled using g2t_stg.F.

.. looking at g2t_stg.F (pams/sim/g2t/g2t_stg.F)

First thing I do in the loop over hits... is print out the hits in this code. I don't see anywhere the output of this printout.

So, we are not getting the copied out of geant.

klendathu2k commented 5 years ago

g2t_stg.F fails because...

 CALL g2r_get_sys('STGH','TGCG',Iprin,Idigi)      
  write (*,*) 'g2r_get_sys returns Iprin=', Iprin, ' Idigi=',Idigi

returns with Iprin < 0, indicating no hits.

klendathu2k commented 5 years ago

CALL should use the module name 'STGM', not the hit bank name 'STGH'.

klendathu2k commented 5 years ago

Issue is resolved, code committed in the CVS repo, will autobuild tonight.