[Frost] Can you check some part in momcor.F?

Sungkyun Park sp06k at fsu.edu
Mon May 17 13:07:12 EDT 2010


Hi Brian

As I already has mentioned, after applying eloss correction, the number of events in CH2 target is decreased and the mean of pull distribution is shifted to minus direction. Until now, I have scanned the eloss to find that reason. I think I have found that point in eloss correction. Can you check what it is?

The problem is happened in le = L(14). L(14) give us the wrong 4 vector, pin.
For test, I make histograms with and without using the following part:
###########################################################

      flag = 6
      le = L(14)
c      write(*,*)"L:",L(1),L(2),L(3),L(4),L(5),L(6),L(7),L(8),L(9)                        
c     +,L(10),L(11),L(12),L(13),L(14),L(15)                                               
c      write(*,*)"targ:",targ(1),targ(2),targ(3)                                          
c      write(*,*)"1     :",le,pout                                                        
      ierr = loss(flag, le, pmass, pout, pin) ! 20k cup end cap     

###########################################################
The result is described in the following web:

 http://hadron.physics.fsu.edu/~skpark/research/research_apr3010.html

In the web, histograms of Before are made with the upper part of eloss correction and histograms of After are made without the upper part. They are made in Topology gamma proton -> proton pi+ pi-  and polyethylene target. Our problems happened in CH2 target can be solved if we do not use the upper part of eloss. 

I think you know what L(14) is. Can you check this variable? 

Sung
Florida State University




More information about the Frost mailing list