run14 simulation

TEST :
  • look at the tracking efficiency with PXL, using code in CVSROOT (StiPxl)
  • sample : 1 muon /event, fixed Pt and phi ~ pi (even with misaligned PXL, 1 track per event should allow the PXL hits to be taken most of the time)
  • loop over global track and plot the # of PXL hits on track via 
  • gTrack->fitTraits().numberOfFitPoints(kPxlId); 
KUMAC :
 detp geometry y2014
 gexec $STAR_LIB/xgeometry.so
 make gstar
 gclose all 
 myPID = 5
 utime=$SHELL('date +%s')
 Rndm [utime] [myPID]
 user/vxyz 0.01 0.01 0.1
 user/vsig 0.01 5
 gfile o test_muon.fz
 aguser/gkine 1 5 1.0 1.0 -1.0 1.0 1.3 1.6
 trig 1000

BFC chain : 
"tpcRS ry2014 AgML MakeEvent ITTF VFMCE pxlFastSim pxlDb pxlUtil PixelIt NoSsdIt NoSvtIt Idst BAna l0 Tree 
logger Sti tpcDB TpcHitMover TpxClu bbcSim btofsim tags emcY2 EEfs evout -dstout IdTruth geantout big fzin 
MiniMcMk clearmem"
note :
the BFC chain has been updated friday (April 11th) with a dedicated option : StiLibsHft ; this option calls StiPxl.
I used PixelIT in this chain : PixelIT calls PxlIT which calls StiPxl, so it should be fine.

Results (over 100 events) :

Fig 1 : PXL hits from PxlHitCollection (all hits, no tracking here)
Fig 2 : TPC hits in tracking distribution : over 100 tracks, 92 are passed after cuts ; distribution looks ok
Fig 3 : PXL hits in tracking distribution : most of the tracks have 0 PXL hits, only 2 tracks have 1 PXL hit, 1 track with 2 PXL hits

Previous Simulation :
similar study with y2013 geometry and StiRnD code : here 

update :
backtracked to the version where StiVolumes are build per ladder :
  1. when sensor !=1, do not build the volume
  2. Set the length in Z of the Sti Volumes to 10*TGeoBBox->GetDZ()
  3. Set the position in Z at 0
  4. old numbering : inner ladders from 0 to 9, outer ladders from 0 to 30 
  5. same numbering in StiPxlHitLoader
Test over 20 muons (17 tracks after cuts)
modif 1 : with 1) and 2)
modif 2 : with 1), 2) and 3) from the above list

# of PXL hits / track 0 1 2 3
CVSROOT 16 1 0 0
modif 1 7 5 5 0
modif 2 2 3 11 2


update (2)
Simulation is :
detp geometry y2014 
 gexec $STAR_LIB/xgeometry.so 
 make gstar 
 gclose all  
 myPID = 8 
 utime=$SHELL('date +%s') 
 Rndm [utime] [myPID] 
 user/vxyz 0.01 0.01 0.1 
 user/vsig 0.01 5 
 gfile o test_pion.fz 
 aguser/gkine 20 8 0.1 5.0 -1.0 1.0 0 6.28 
 trig 100 
Fig 4 : (# of tracks with PXL >0) / (# of all tracks)  vs pt
efficiency -> .8 (both numerator and denominator have a cut TPC>15)
Fig 5 : same of as Fig 4 but both numerator and denominator have a cut |eta|<.5 
efficiency -> .85