[Frost] About MVRT bank

Eugene Pasyuk pasyuk at jlab.org
Thu Apr 22 13:32:33 EDT 2010


Hi Sung,

For you first point, NO. You don't compare between GPID and MVRT. Ypou 
compare z-vertexes in GPID for protons an pions.

2 and 3 will be unnecessary if you require mass_ref=0 for all particles 
you are interested in and they all have the same tagrid.

So, for pi+pi-p events you do:
1. mass_ref=0 for all particles you need.
2. targrid is the same for all particles
3. check if difference between z-vertexes of particles you need is less 
than some number.

This would be the initial criteria for the good event.

You use GPID vertex to find that tracks are coming from the same vertex. 
Once you know that this is the case, you use MVRT vertex for everything 
else.

-Eugene

On 4/22/10 1:16 PM, Sungkyun Park wrote:
> Thanks Eugene,
> I understand.
>
> I will add more cuts in my code.
>    1. compare the vertex position between GPID.vertex and MVRT.vertex
>    2. compare the sc_time with st_time
>    3. check MVRT.chi2
>    4. compare mass_ref of the each particle.
>
> I have one more question about vertex.
> Do I use the MVRT vertex information as the input parameter of ELoss Correction?
> Until now, I used the GPID vertex information.
> If so, Do I need to use the MVRT vertex information when I need to use the vertex?
> Thanks very much.
>
> Best Wishes
>
> Sung
>
>
> ----- Original Message -----
> From: Eugene Pasyuk<pasyuk at jlab.org>
> Date: Wednesday, April 21, 2010 10:05 pm
> Subject: Re: About MVRT bank
> To: Sungkyun Park<sp06k at fsu.edu>
> Cc: FROST<frost at jlab.org>
>
>> Hi Sung,
>>
>> You found very good example of accidental tracks.
>> First of all, as you noticed the vertexes are different. This is an
>> indication that they may have nothing to do with each other.
>> Furthemore, if you look at the timing information you see that for
>> the
>> "proton" ST time is way off, it is even larger than  SC time.
>> Also you see that mass_ref is not good for it. All together makes
>> it a
>> bad track. chi2 for the fitted vertex in MVRT bank is also
>> terrible, 45!
>> When you select multi-track events you have to make sure that they
>> all
>> belong to the same interaction, apply cuts on vertex time and
>> vertex
>> position differences. Also check that they have good pid status
>> (mass_ref amand they matched to the same photon).
>>
>> -Eugene
>>
>>
>>
>> On 4/21/10 2:29 PM, Sungkyun Park wrote:
>>> Hi Eugene,
>>>
>>> In today meeting, you said that I need to use the vertex
>> information in MVRT bank.
>>> Until now, I use the vertex in GPID bank.
>>>
>>> I have a question about MVRT bank.
>>>
>>> In my code, the part to get the vertex information of the proton
>> is the following:
>>> 	
>>>           if( GPID[row].pid == 14 ){
>>>
>> My_pro_vertex.SetXYZ(GPID[row].x,GPID[row].y,GPID[row].z);>
>> }
>>>
>>> 1. Do the MVRT bank not have any information for the particle ID?
>>>
>>> In the event with a proton and pi-,
>>> the GPID vertex for proton is   x: -0.000       y: 0.019
>> z: 2.180   ,
>>> the GPID vertex for pi- is  x: 0.000        y: 2.320        z: -
>> 9.450    , and
>>> the MVRT vertex is  x: 3.328  y: 1.341  z: 2.930
>>> 2. I do not understand why this event has two GPID vertexes and
>> one MVRT vertex.
>>> Who is this MVRT vertex the vertex for? How do I use the MVRT
>> vertex in this case.
>>>
>>> The following is the outputs of the GPID and MVRT bank in the
>> event with a proton and a pi-
>>> ------------------------------------------------------------------
>> -------------------------------------------------------
>>> Group:  GPID    Sector: 1       Nhits:  2       Next ind:       0
>>> pid:     14
>>>      vert->         x: 0.000        y: 2.320        z: -9.450
>>>      p->            E: 1.106        px: 0.151       py: 0.010
>> pz: 0.565
>>>      q:   1       trkid:   1      sec:   1        paddle:  16
>>>      dedx: 33.069 beta: 0.529
>>>      sc_stat:   0 sc_time: 40.619         sc_len: 526.997
>>>      st_stat:   1 st_time: 57.143         st_len: 34.069
>>>      mass: 0.882  mass_ref:   2   betam: 0.553
>>>      epho: 0.726  tpho: 8.749     tagrid:   2     ngrf:   1
>> ppid:  14
>>> pid:      9
>>>      vert->         x: -0.000       y: 0.019        z: 2.180
>>>      p->            E: 0.311        px: -0.270      py: -0.023
>> pz: -0.064
>>>      q:  -1       trkid:   2      sec:   4        paddle:  41
>>>      dedx: 11.189 beta: 0.894
>>>      sc_stat:   2 sc_time: 24.078         sc_len: 414.331
>>>      st_stat:   1 st_time:  9.746         st_len: 10.415
>>>      mass: 0.130  mass_ref:   0   betam: 0.906
>>>      epho: 0.726  tpho: 8.749     tagrid:   2     ngrf:   1
>> ppid:   9
>>>
>>> Group:  MVRT    Sector: 1       Nhits:  1       Next ind:       0
>>> v_id:      6  ntrk:      2  x: 3.328  y: 1.341  z: 2.930  chi2:
>> 45.419  cxx: 0.003  cxy: -0.001  cxz: 0.001  cyy: 0.027  cyz: -
>> 0.001  czz: 0.002  stat:      0
>>>
>>> Best Wishes
>>> Sung
>>
> _______________________________________________
> Frost mailing list
> Frost at jlab.org
> https://mailman.jlab.org/mailman/listinfo/frost


More information about the Frost mailing list