120 likes | 244 Views
Test of TCFIT in reco. code (update ). Test on pure D 0 ( ) Test on embedding : D 0 mixed with HIJING events Test with Real data. Mass fit. Goal : obtain a good estimation of the signal (= number of counts under the inv. Mass peak). To do this ,we have Fit of mass+background :
E N D
Test of TCFIT in reco. code(update ) Test on pure D0 () Test on embedding : D0 mixed with HIJING events Test with Real data
Mass fit • Goal : obtain a good estimation of the signal (= number of counts under the inv. Mass peak). • To do this ,we have • Fit of mass+background : • Linear + BTW (STAR paper K*0) where the fit function is : a + bMreco+ Γ/[(Mreco -MPDG)2 + (Γ/2)2 ] • similarly, there is a linear+lorentzian function (from ROOT tutorials) which the form is : • a + bMreco+cMreco2 + (c*d/2π)/Max[(Mreco -MPDG)2 + (d/2)2] • Substraction(K-π+ + K+π-) by (K+π+ + K-π-) • Refit by linear + gaus after if necessary
Using the ROOT tutorial fit(i) No cut on the #Tracks ; pT≥3 && pT<4
Using the ROOT tutorial fit(ii) #Tracks<200 ; pT≥3 && pT<4
Using the substraction • # tracks <100 • After substraction, a gaussian fit is enough to have the inv. Mass • We can also see the the width from the the gaussian fit is agree with the PDG
Using the substraction • as a function of the cut # tracks No cut #tracks<200 #tracks<100
Using the substraction • as a function of pt pT<1 pT≥1 && pT<2 pT≥1 && pT<2 pT≥3 && pT<4
Set-up / cuts • Slow simulator • 8.5 kEvents with 1 D0/events • Cuts in MuKpi : • |Zvertex| < 30 cm • NoTracks<700 • SiliconHits>0 • p >.1GeV/c • |slength(linear)|<.15 • |nK|<3, |nπ|<3
details • Took macros used in HFT soft. • Use y2007g geometry (last updated for AuAu in simu : • http://drupal.star.bnl.gov/STAR/comp/prod/MCGeometry • Follow these different steps : • Produce AuAu MB 200GeV as background • Produce 1 D0 per event • Mix fz files from background and D0 • Run the bfc to do the reco (with slow simulation) • Analyze MuDst.root file (ie all reco. tracks) (with same macro used for real data analysis)
AuAu background detp geometry y2007g make geometry make gstar … swit 2 3 vsig 0.01 20.0 gkine -1 0 0 5 -6.3 6.3 0 6.28 -40.0 40.0 gfile o outfile.fz user/input user infile.nt trig 50 • Y2007g geom. • Vertex at 0.01 with = 20 cm • AuAu collisions in 0<pT<5 • 50 events
D0 production #***************************************************** # Input for gkine set geometry=y2007g set pid=$3 set tracks_per_event=1 set ptlow=0.1 set pthigh=5.0 set ylow=-1.1 set yhigh=1.1 set philow=0.0 set phihigh=6.28 #***************************************************** • This step produces a .fz file with D0’s and used a • vertex file (containing the PV point) for correct mixing