[Frost] About momcor.F in eloss correction
Sungkyun Park
sp06k at fsu.edu
Fri May 14 16:52:29 EDT 2010
Hi Eugene,
I have a question in momcor.F of eloss correction.
The following part is the energy loss correction for the air and Start Counter in FROST:
(There is the momcor.F code part in ###.)
###########################################################
mate = 8 ! air
ierr = loss(mate, gap, pmass, po, p)
do i = 1,3
pi(i) = p(i)
enddo
mate = 2 ! Start Counter scintillator
ierr = loss(mate, dist, pmass, pi, p)
do i = 1,3
pi(i) = p(i)
enddo
#############################################################
I know pout is the input 4 vector and pin is the output 4 vector.
The last variable for the 4 vector in the upper part is pi. That is, eloss calculates the energy loss in the air and Start Counter and save the new 4 vector in the pi. But In the part for FROST,
###############################################################
if (icell .eq. 9) then
call FROST(vertex, cdir, L) ! target cell
call targets(vertex, cdir, targ) !targets
flag = 2
le = L(17) / 10
ierr = loss(flag, le, pmass, pout, pin) ! scattering chamber end
do i=1,3
pout(i) = pin(i)
enddo
###########################################################
eloss use pout as the input 4 vector, not pi.
I check this pout is the initial 4 vector we give an eloss program.
We do not need to use the eloss part for air and Start Counter so do we use the pout as the initial variable in the FROST part? or is it incorrect?
Sung
More information about the Frost
mailing list