FWD Software

["Brandenburg, Daniel" <brandenburg.89 AT osu.edu>]

Hi Akio, Lets discuss this today, Want to make sure I know what to look at Daniel Brandenburg Assistant Professor College of Arts and Sciences Department of Physics Ohio State University From: star-fwd-software-l-request@ lists. bnl. gov <star-fwd-software-l-request@ lists. bnl. gov>

ZjQcmQRYFpfptBannerStart

This Message Is From an External Sender

This message came from outside your organization.

 

ZjQcmQRYFpfptBannerEnd

Hi Akio,

Lets discuss this today, Want to make sure I know what to look at 

Daniel Brandenburg

Assistant Professor

College of Arts and Sciences

Department of Physics

Ohio State University

From: star-fwd-software-l-request AT lists.bnl.gov <star-fwd-software-l-request AT lists.bnl.gov> on behalf of Ogawa, Akio <akio AT bnl.gov>
Date: Wednesday, June 3, 2026 at 2:35 PM
To: "Brandenburg, Daniel" <star-fwd-software-l AT lists.bnl.gov>
Subject: [[Star-fwd-software-l] ] Pythia Jpsi sample

This Message Is From an External Sender

This message came from outside your organization.

Report Suspicious

 

Hello

I did produce some Pythia Jpsi samples:

/gpfs01/star/pwg_tasks/FwdCalib/akio/pythia/pythia.JPsi.vz0/pythia.JPsi.vz0.runXXX.fzd

and .MuDst.root and .picoDst.root inn the same dir, where XXX is 1~100 and each contains 200events.

I do get nice Jpsi peak from FCS alone, but no good track matchings.

https://www.star.bnl.gov/protected/spin/akio/fcs/jpsi/20260603/dilep2...1.png

https://www.star.bnl.gov/protected/spin/akio/fcs/trackmatch/20260603/index.html

Now I went to single electron gun (pt=1.5GeV, eta=2.5~4, full phi) events:

cd  ~/fcstrk11/star-sw-fwd/

root4star -b -q runSimFlat.C'(1000)’

root4star -b -q sim.C > log_sim

root4star -b -q pico3.C'(0,11,0)' > log_pico

Here is a part of log_pico on tracks and Ecal clusters:

ev=   1 trk=  0 Type=0 NHit=  8 Cvg=1 Chi2=     9.3 pxyz=  -0.20    0.26    3.53 E=   3.55 ProjE=   1.04   -0.12 ( -43.01   55.74)

ev=   1 trk=  1 Type=2 NHit= -9 Cvg=1 Chi2=    14.3 pxyz=  -0.79    1.17   15.36 E=  15.42 ProjE=   0.07    0.04 ( -38.81   56.99)

ev=   1 trk=  2 Type=1 NHit=  9 Cvg=1 Chi2=    11.5 pxyz= -10.14   14.16  187.51 E= 188.31 ProjE=   0.35   -0.17 ( -40.56   56.63)

ev=   1 clu=  0 det=0 ehp=0 xyz= -38.25   56.02 E=15.396 SigMax/Min= 0.515  0.380 Ctg=1

ev=   3 trk=  0 Type=0 NHit=  8 Cvg=1 Chi2=     3.2 pxyz=  -0.33   -0.64    9.75 E=   9.77 ProjE=  -0.18    0.25    0.00 ( -25.30  -49.61)

ev=   3 trk=  1 Type=2 NHit= -9 Cvg=1 Chi2=     4.1 pxyz=  -1.01   -1.86   28.63 E=  28.71 ProjE=  -0.03    0.02    0.00 ( -26.46  -48.67)

ev=   3 trk=  2 Type=1 NHit= -9 Cvg=1 Chi2=     3.6 pxyz= -10.50  -19.85  302.82 E= 303.65 ProjE=  -0.03    0.13    0.00 ( -26.01  -49.17)

ev=   3 clu=  0 det=0 ehp=0 xyz= -25.25  -48.80 E=19.674 SigMax/Min= 0.466  0.390 Ctg=1

Last 2 numbers for track in () are 

   trk->momentum().X()/trk->momentum().Z()*750.0

   trk->momentum().Y()/trk->momentum().Z()*750.0

as poor man's xy projection to Ecal based on momentum.

I’m using Danie’s dev code and his sim.C macro as well. I’m not sure what I am doing wrong with Ecal projection. 

As far as I can see… pico code looks correct about z position to project to, x100 to store as short, and /100.

It seems like something wrong upstream...

-Akio