One of the primary goals of my research this summer was to incorporate a 'C' encoded program for the analysis of lambda (1520) particles into the Fortran encoded RECON (RECSIS) software. The 'C' program was written by Steve Barrow, and was used mainly in an independent fashion for analysis of lambda (1520) particles. Incorporating this program into RECON, then, involved isolating the appropriate sections of the 'C' code and calling them from the appropriate files in the Fortran code. Below is where the primary event processing code for lambdas is called: File : user_evnt.F (RECON)
The above code handles the processing of the lambdas by calling the procedure pe_sbarrow(). This procedure name is essentially an abbreviation for "process event (Steve Barrow's code)." Involved with this code is the procedure "pl_1520_sbarrow()" which is more specific code for the lambdas. Again, the abbreviation as such is for "process lambda 1520 particles (Steve Barrow's code)." The processed lambda events are also entered into histograms, which are 'booked' in a process called bh_sbarrow() in the same code as processevent. The abbreviation in this case stands for "book histograms (Steve Barrow's software)." This section of Barrow's code is called in the user_init.F file of RECON.
File : user_init.F (RECON)
The links in the previous section of my text both link to the full 'C' code used (authored by Steve Barrow) in the incorporation. Once this was accomplished, and all the bugs were worked out of the code, the software was run with the following input files:
Track (trkmon) Histograms:
Histogram trkmon 100 gives a typical output of the number of tracks seen in a particular event. As can be observed, the software was modified to the point that approximately 50% of the events contained tracks (~500/1000.) Histogram trkmon 105 gives some information about the W (GeV) values detected in the events with tracks, and a distinct peak at ~2 GeV can be clearly seen. This peak corresponds to the lambda (1520) particle (actually a lambda and a kaon.) This histogram represents the first sector only of the CLAS detector, but is representative of the data obtained in the other sectors as well. The last histogram, trkmon 109, is a Q^2 versus W relationship for the lambda events, again, for the first sector only. However, this is also fairly representative of the data in the other sectors as well.
The actual events can be viewed with CED, and several events can be observed with the links below. The curvature of the tracks, therein, is the result of the toroidal magnets in the CLAS that bend the tracks of charged particles, which assists in their identification.
|