Single pion simulation - 05/14/2014 

1) Simulation setup

a) Single primary pion fz files:

/star/u/yfzhang/lbl/hft/HFT_simu/single/pion/events/fz
file 0-9: 5 event per file, 10 pions per event.
file 10-49: 5 event per file, 100 pions per event. (used for this simulation)
file 50-59: 5 event per file, 300 pions per event.

Geometry: y2014
flat pT 0-10 GeV/c, flat eta (-1,1), flat full phi
Vertex:
With sigma = 40 cm for Vz, 0.01 cm for (Vx,Vy). 
|Vz| < 5cm
 

b) Bfc chain:

/star/u/yfzhang/lbl/hft/HFT_simu/single/bfc

StRoot: follow offline/hft/README.md
StarDb from Amilkar, but using tuned hit errors: 

for IST:
  Double_t sigma_x = 0.017; 
  Double_t sigma_z = 0.17;
for PXL:
  Double_t sigma_x = 0.0012; 
  Double_t sigma_z = 0.0012; 

Chain option used:

tpcRS,y2014,StiPulls,AgML,MakeEvent,tables,SimuDb,VFMinuit,pxlFastSim,istFastSim,Idst,BAna,l0,Tree,logger,Sti,StiHftC,
-SstIT,-ssdDb,TpcHitMover,TpxClu,bbcSim,btofsim,-tags,emcY2,EEfs,evout,geantout,IdTruth,geant,big,fzin,clearmem

SSD geometry is in but no SSD hits loaded for tracking.
 

No pile-ups so far.
 

c) Analysis ntuples:

/star/u/yfzhang/lbl/hft/HFT_simu/single/bfc/ntuple/out

2) Efficiencies

a) Mc efficiency:
Black: dN_dpt(pions with at least 2 PXL MC hits) / dN_dpt(MC pions)
Red:   dN_dpt(pions with at least 1 IST MC hits) / dN_dpt(MC pions)
Blue:  dN_dpt(pions with at least 1 IST MC hits and 2 PXL MC hits) / dN_dpt(MC pions)

b) Mc efficiency with secondary pions removed
Same as a), but require pions are all from primary vertex.

c) MC efficiency as a function of eta (eta acceptance), definition is the same as b)

d) MC efficiency as a function of phi (phi acceptance), definition is the same as b) with |eta|<0.5

e) Reconstructed efficiency:

Black: dN_dpt(pions with nTPCFitPts > 20) / dN_dpt(MC pions)
Red:   dN_dpt(pions with nTPCFitPts > 20 & IST fired) / dN_dpt(MC pions)
Blue:  dN_dpt(pions with nTPCFitPts > 20 & IST fired & at least 2 PXL hits) / dN_dpt(MC pions)