[Clas12_verystrange] GEMC 1.3.1
Veronique Ziegler
ziegler at jlab.org
Wed Jun 22 14:14:59 EDT 2016
Hi,
it looks like it was reconstructed without the CVT service. In the event bank the momentum is in GeV/c and the vertex in cm.
The law efficiency is due to missing the central detector and also to event where the tracks bend outside of the fiducial detector. You could look at the reconstructed file in ced to see what is happening in terms of track lost.
Veronique.
----- Original Message -----
From: "Nicholas Zachariou" <nickzachariou at gmail.com>
To: "Veronique Ziegler" <ziegler at jlab.org>, "Nathan Harrison" <nathanh at jlab.org>, "Maurizio Ungaro" <ungaro at jlab.org>
Cc: "Vicente Pais" <dav30 at hw.ac.uk>, "clas12 verystrange" <clas12_verystrange at jlab.org>
Sent: Tuesday, June 21, 2016 7:27:26 AM
Subject: Re: GEMC 1.3.1
Hi Veronique,
Judging from the GEMC output it looks like it reads the LUND file correctly now (run with version 2.4). There seems to be issues though with the reconstruction.
We generated 5000 events (gamma d —> K+ Lambda; Lambda decayed in the generator) processed through GEMC and reconstructed. Out of the 5000 events there is no single event with 3 charged track reconstructed and only 39 events with a proton and pion. It looks like the EVIO information for the GenPart::true bank has different units from the EVENTHB::particle bank (for both momenta and vertices) and the vertex reconstruction in EVENTHB seem way off (PID also fails quite often). See attached for two examples.
Files are located at
They are in this directory:
/volatile/clas/clasg13/duartep/test
The files are the following:
outFile.txt: Generated Lund File
out.ev: GEMC evio
GEMCout.dat: GEMC output
myoutput.0.evio: Reconstructed evio
Let me know what you think,
Nick
Example 1:
Generated event (LUND format):
4 1 1 0 0 0 0 1 0 8.87824
1 0 11 3122 0 -1 -1.62786 -0.831217 6.00681 6.3771 1.11568 -0.15508 0.0800463 -1.23106
2 0 1 321 0 -1 1.62786 0.831217 2.87143 3.43941 0.4936 0.0169529 0.0315987 -1.26968
3 0 1 2212 1 -1 -1.28489 -0.601995 4.94704 5.23133 0.93827 -10.0265 -5.22436 36.1426
4 0 1 -211 1 -1 -0.342966 -0.229222 1.05978 1.14577 0.13957 -9.97781 -5.12122 36.3364
GEMC output:
Beam Settings >> Particle Number: 2, id=321 (kaon+) Vertex=(0.0169529,0.0315987,-1.26968)cm, momentum=(1.62786,0.831217,2.87143) GeV
Beam Settings >> Particle Number: 3, id=2212 (proton) Vertex=(-10.0265,-5.22436,36.1426)cm, momentum=(-1.28489,-0.601995,4.94704) GeV
Beam Settings >> Particle Number: 4, id=-211 (pi-) Vertex=(-9.97781,-5.12122,36.3364)cm, momentum=(-0.342966,-0.229222,1.05978) GeV
Reconstructed EVIO:
Press Enter for Next Event or Bank Name: GenPart::true
*****>>>>> BANK GenPart::true >>>> SIZE = 7
pid (int) : 321 2212 -211
px (double) : 1627.86000 -1284.89000 -342.96600
py (double) : 831.21700 -601.99500 -229.22200
pz (double) : 2871.43000 4947.04000 1059.78000
vx (double) : 0.16953 -100.26500 -99.77810
vy (double) : 0.31599 -52.24360 -51.21220
vz (double) : -12.69680 361.42600 363.36400
Press Enter for Next Event or Bank Name: EVENTHB::particle
*****>>>>> BANK EVENTHB::particle >>>> SIZE = 12
status (int) : 100
charge (int) : 1
pid (int) : 0
chi2pid (float) : 0.00000
mass (float) : 0.00000
beta (float) : 0.00000
px (float) : 1.61035
py (float) : 0.81359
pz (float) : 2.84154
vx (float) : 0.90305
vy (float) : 0.23706
vz (float) : -0.42210
Example 2:
Generated event:
4 1 1 0 0 0 0 1 0 7.88471
1 0 11 3122 0 -1 1.19836 -1.32457 4.23444 4.72924 1.11568 0.00658124 0.092797 3.10574
2 0 1 321 0 -1 -1.19836 1.32457 3.65027 4.09373 0.4936 -0.154459 0.0929773 3.27535
3 0 1 2212 1 -1 1.09197 -1.07384 3.51083 3.94358 0.93827 5.59435 -6.13771 23.27
4 0 1 -211 1 -1 0.10639 -0.250726 0.72361 0.785668 0.13957 5.83073 -6.27152 23.2612
GEMC output:
Beam Settings >> Particle Number: 2, id=321 (kaon+) Vertex=(-0.154459,0.0929773,3.27535)cm, momentum=(-1.19836,1.32457,3.65027) GeV
Beam Settings >> Particle Number: 3, id=2212 (proton) Vertex=(5.59435,-6.13771,23.27)cm, momentum=(1.09197,-1.07384,3.51083) GeV
Beam Settings >> Particle Number: 4, id=-211 (pi-) Vertex=(5.83073,-6.27152,23.2612)cm, momentum=(0.10639,-0.250726,0.72361) GeV
Reconstructed EVIO:
Press Enter for Next Event or Bank Name: GenPart::true
*****>>>>> BANK GenPart::true >>>> SIZE = 7
pid (int) : 321 2212 -211
px (double) : -1198.36000 1091.97000 106.39000
py (double) : 1324.57000 -1073.84000 -250.72600
pz (double) : 3650.27000 3510.83000 723.61000
vx (double) : -1.54459 55.94350 58.30730
vy (double) : 0.92977 -61.37710 -62.71520
vz (double) : 32.75350 232.70000 232.61200
*****>>>>> BANK EVENTHB::particle >>>> SIZE = 12
status (int) : 100 100
charge (int) : 1 -1
pid (int) : 0 0
chi2pid (float) : 0.00000 0.00000
mass (float) : 0.00000 0.00000
beta (float) : 0.00000 0.00000
px (float) : -1.20575 0.01975
py (float) : 1.29930 -0.26616
pz (float) : 3.63421 0.76682
vx (float) : 0.74334 7.58224
vy (float) : -0.49951 -1.89821
vz (float) : 0.37457 -2.53571
> On 21 Jun 2016, at 11:15, Vicente Pais, Duarte <dav30 at hw.ac.uk> wrote:
>
> They are in this directory:
> /lustre/expphy/volatile/clas/clasg13/duartep/test
>
> The files are the following:
>
> Generated Lund File -> outFile.txt
>
> GEMC evio -> out.ev
>
> GEMC output -> GEMCout.dat
>
> Reconstructed evio -> myoutput.0.evio
>
> If you want to look at the generated and reconstructed events, the file is output groovy.txt. I went through the generated ones and the vertex positions are more or less multiplied by ten, again…
>
> And those 39 entries give respect to the events with pions and protons. Still no event is completely reconstructed…
>
> Duarte
>
>> On 21/06/2016, at 10:44, Nicholas Zachariou <nickzachariou at gmail.com <mailto:nickzachariou at gmail.com>> wrote:
>>
>> Can you put the files in a directory at volatile and send me the link?
>> I want the Generated Lund file, GEMC evio, GEMC output, and Reconstructed EVIO. Is the Generated vertex in the reconstructed evio OK?
>>
>> Do we reconstruct 39 events now or is it just events with proton and pions?
>>
>> -N
>>
>>
>>> On 21 Jun 2016, at 10:42, Vicente Pais, Duarte <dav30 at hw.ac.uk <mailto:dav30 at hw.ac.uk>> wrote:
>>>
>>> Nick,
>>>
>>> the reconstructed file still has the vertex positions with zeroes… This is the same plot for the invariant mass
>>>
>>> <Screen Shot 2016-06-21 at 10.40.19.png>
>>> Founded in 1821, Heriot-Watt is a leader in ideas and solutions. With campuses and students across the entire globe we span the world, delivering innovation and educational excellence in business, engineering, design and the physical, social and life sciences.
>>>
>>> The contents of this e-mail (including any attachments) are confidential. If you are not the intended recipient of this e-mail, any disclosure, copying, distribution or use of its contents is strictly prohibited, and you should please notify the sender immediately and then delete it (including any attachments) from your system.
>>>
>>
>
More information about the Clas12_verystrange
mailing list