sPHENIX EMCal discussion

[eric mannel <mannel AT bnl.gov>]

Thanks to Martin Frank from the mu2e collaboration, there is a tracking tracking package that can reconstruct tracks using the FTBF chambers.  Code for the package can be found in /home/nfs/ftbf_user/experiments/t1044/mannel/mfrank on ftbfbnk01.

This version of code has the third station shifted 10mm in X to account for an unexplained shift in the chamber that seems to have occurred recently.

There is a README in the area that explains how to run his package, but it straightforward. 

1) /home/nfs/ftbf_user/experiments/t1044/mannel/mfrank]$ source setup.sh
2) /home/nfs/ftbf_user/experiments/t1044/mannel/mfrank]$ ./mwpc -M input_data_file -H output_hist_file -T output_tree_file

This version of code uses the FBTF ascii format as input. 

On the T-1044 side, Martin P has setup  automatic running of the MWPC daq that generates a "prdf" file in rcdaq and for the past week we have been logging the mwpc data for each run with some constancy.  In order to analyze the data, I have modified a version of Martins prdf2text routine to produce a file that can be used as an input into Martin Frank's tracking package.  Right now it runs as a root macro:

1) root> .L mwpcprdf2text.C
2) root> mwpcprdf2text("input_file", "output_file")

Copies of the macro can be found in ftbfbnl01:/home/nfs/ftbf_user/experiments/t1044/mannel/mwpcprdf2text.C and hcaldaq:/home/phnxsa/mannel/mwpcprdf2text.C

Feel free to copy any or all of this into your own space to work with it.

From Martin Frank:

Begin forwarded message:

From: "Frank, Martin Johannes (mjf8u)" <mjf8u AT eservices.virginia.edu>
Subject: Mu2e MWPC Code Date: April 21, 2016 at 10:33:05 AM CDT
To: Mickey Chiu <chiu AT bnl.gov> Cc: Ewa E Skup <eskup AT fnal.gov>
Hi Mickey,

Here are some more details about our wire chamber tracking code.
1. You can find my code installed on the ftbf hosts in ~ftbf_user/experiments/mfrank/mwpc/.  In order to build and run the code simply look at the README file.  Here is the command that we used to run over your file:

 ./mwpc -C -m ~/experiments/t1048/data/exp_20160421063354.dat -H mickey_hist_t1048.root


2. You can browse the code in CVS here:

 http://cdcvs.fnal.gov/cgi-bin/public-cvs/cvsweb-public.cgi/mu2e/crv/mwpc/

3. In order to get the code to work for you, I would first adjust the TDC cuts in Constants.hh until the calibration histogram looks good for all TDC boards:

 http://cdcvs.fnal.gov/cgi-bin/public-cvs/cvsweb-public.cgi/mu2e/crv/mwpc/Constants.hh?rev=1.15&content-type=text/x-cvsweb-markup

4. After that, you might want to get rid of the chi^2 cut in our track fit.  Have a look at the end of Track::fit() in Track.cc:

 http://cdcvs.fnal.gov/cgi-bin/public-cvs/cvsweb-public.cgi/mu2e/crv/mwpc/Track.cc?rev=1.5&content-type=text/x-cvsweb-markup

5. Once this lets more events in, you can try to align your chambers using the mwpc::Constants::STATION_OFFSET array.  Look at the track-dx/dy-station-* histograms to minimize your residual.

I originally developed the code on my laptop, so if you have a Mac, you can also build and run it on there.

Feel free to use as much of my code as you want and let me know if anything is unclear or you are having trouble getting it to run.
Best of luck!
Martin

-- 
Eric J. Mannel, Ph.D.
PHENIX Group
Physics Dept.
Brookhaven National Laboratory
631/344-7626 (Office)
914/659-3235 (Cell)