<html>
<head>
</head>
<body dir="auto">
<div>Has this always been the case? If so, is it generally known? If not, this is a serious issue which should be addressed, since we need to know how many events we're simulating if we want to calculate our acceptances correctly.</div>
<div><br>
</div>
<div>John<br>
<br>
--
<div>John Price</div>
<div>Professor and Chair, CSUDH Department of Physics</div>
<div>Co-Director, Office of Undergraduate Research, Scholarship, and Creative Activity</div>
<div>(310) 243-3403</div>
</div>
<div><br>
On May 20, 2014, at 11:07, "Haiyun Lu" <<a href="mailto:hlu@jlab.org">hlu@jlab.org</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>
<div class="moz-cite-prefix">Hi Nick,<br>
<br>
I believe it has nothing to do with gsim and fsgen is the cause. Did you check how many events you got out of fsgen? I bet it's 71 events. You better check the evt file directly instead of the log file. To my knowledge, Fsgen makes 100 loops to generate events
when you ask for 100 events. However, some events are discarded due to some reason, like it's not physical. Fsgen will NOT make up for these events. Therefore you only get a portion of what you ask for.<br>
<br>
Haiyun<br>
<br>
On 05/20/2014 01:41 PM, Johann Goetz wrote:<br>
</div>
<blockquote cite="mid:CAMfAbkighta3k=xvF498m5crqNDDHn7gyrtQPJS6402V7VLSAQ@mail.gmail.com" type="cite">
<div dir="ltr"><font color="#000000" face="arial, helvetica, sans-serif">Hi All,</font>
<div><font color="#000000" face="arial, helvetica, sans-serif">Nick Compton is having what looks like a very similar issue to Markov's (see below). Was this ever resolved?<br clear="all">
</font>
<div><font color="#000000" face="arial, helvetica, sans-serif">--
<a href="https://sites.google.com/site/theodoregoetz/">
Johann</a></font></div>
<div><font color="#000000" face="arial, helvetica, sans-serif"><br>
</font></div>
<div>
<p><font color="#000000" face="arial, helvetica, sans-serif"><br>
</font></p>
<p><font color="#000000" face="arial, helvetica, sans-serif">########## begin message from Nick Compton #########</font></p>
<p><font color="#000000" face="arial, helvetica, sans-serif">when I run a Monte Carlo simulation say of 100 events, I only get 71 back in the root tree. This may be what I expect
from a detector simulation depending on the reaction, except for the fact that in the mc bank I am also only getting 71 events. This seems odd as I would not expect the number in the mc bank to be the same as the events bank. Do you know if there is an option
I have to use with GSIM or in another part of the process to have the mc bank work properly? Please let me know of any insight you might have.</font></p>
<p><font color="#000000" face="arial, helvetica, sans-serif"> </font></p>
<p><font color="#000000" face="arial, helvetica, sans-serif">I am using fsgen to generate .evt files which according to the .log file produced 100 events. The information is
stored in the directory '/volatile/clas/clasg10/ncompton'.</font></p>
<p><font color="#000000" face="arial, helvetica, sans-serif"> </font></p>
<p><font color="#000000" face="arial, helvetica, sans-serif">Thanks,</font></p>
<p><font color="#000000" face="arial, helvetica, sans-serif">Nick</font></p>
</div>
<div><font color="#000000" face="arial, helvetica, sans-serif"><br>
</font></div>
<div><font color="#000000" face="arial, helvetica, sans-serif"><br>
</font></div>
<font color="#000000" face="arial, helvetica, sans-serif">########### begin thread of emails associated with Markov's problems ########<br>
<br>
</font>
<div class="gmail_quote"><font color="#000000" face="arial, helvetica, sans-serif">---------- Forwarded message ----------<br>
From: <b class="gmail_sendername">Michael C. Kunkel</b> <<a href="mailto:mkunkel@jlab.org">mkunkel@jlab.org</a>><br>
Date: Fri, May 16, 2014 at 9:20 AM<br>
Subject: Re: [Clas_offline] Generated events in the PART bank (PART bank sector 0) after user_ana.<br>
To: nick markov <<a href="mailto:markov@jlab.org">markov@jlab.org</a>>, KiJun Park <<a href="mailto:parkkj@jlab.org">parkkj@jlab.org</a>><br>
Cc: clas offline <<a href="mailto:clas_offline@jlab.org">clas_offline@jlab.org</a>><br>
<br>
<br>
</font>
<div bgcolor="#FFFFFF" text="#000000"><font color="#000000" face="arial, helvetica, sans-serif">Greetings,<br>
<br>
From what I recall, the STABLE build of RECSIS/user_ana will not output PART 0. However, I remember Stepan having a build that does output PART 0. Maybe you can compare the source codes and figure out why.<br>
<br>
BR<br>
MK</font>
<div>
<div><font color="#000000" face="arial, helvetica, sans-serif"><br>
<br>
<br>
</font>
<div><font color="#000000" face="arial, helvetica,
sans-serif">On 5/16/14 9:13 AM, nick markov wrote:<br>
</font></div>
<blockquote type="cite"><font color="#000000" face="arial, helvetica,
sans-serif">Hi,
</font>
<div><font color="#000000" face="arial, helvetica,
sans-serif">thank you very much for your respond. The tcl file which I use is :</font></div>
<div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> source /group/clas/builds/test3/src/clas6-trunk/reconstruction/recsis/recsis_proc.tcl;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> turnoff ALL;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> global_section off;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> turnon seb trk cc tof egn lac user pid;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> inputfile gpp_out.bos;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> setc chist_filename recsis.hbook;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> setc log_file_name out.log;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> outputfile test.recsis PROC1 2047;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set dc_xvst_choice 0;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set torus_current 2250;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set mini_torus_current 6000;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set poltarget_current 0;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> setc outbanknames(1) "PARTEVNTECPBHEADDCPBTBIDTBTRTDPLTBLACCPBTRGSSCPBTGBI";</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"><br>
</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> setc prlink_file_name "prlink_NEW_60_75.bos";</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> setc bfield_file_name "bgrid_T67to33.fpk";</font></div>
<p><font color="#000000" face="arial, helvetica,
sans-serif"> <br>
</font></p>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lseb_nt_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lall_nt_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lscr_nt_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lseb_hist -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lseb_h_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lmon_hist -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set ltrk_h_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set legn_h_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set ltof_h_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lec1_h_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lfec_hist -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lfec_h_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lpart_nt_do -1;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> fpack "timestop -9999999999";</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set lscat $false;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> set ldisplay_all $false;</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> go </font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"> exit_pend;</font></div>
</div>
<div><font color="#000000" face="arial, helvetica,
sans-serif">Best regards, Nick.</font></div>
<div><font color="#000000" face="arial, helvetica,
sans-serif"><br>
</font>
<div>
<div><font color="#000000" face="arial, helvetica,
sans-serif">On May 16, 2014, at 9:08 AM, KiJun Park <<a href="mailto:parkkj@jlab.org">parkkj@jlab.org</a>>
wrote:</font></div>
<font color="#000000" face="arial, helvetica,
sans-serif"><br>
</font>
<blockquote type="cite"><font color="#000000" face="arial, helvetica,
sans-serif">Hello,<br>
<br>
Interesting. Is this matter of creating nPRT in Ntuple after RECSIS ?<br>
For example, <br>
* 1 * I*4 * * [0,20] * PART * nPRT<br>
* 2 * I*4 * * * PART * pidPART(nPRT)<br>
These are empty ?<br>
If so, I believed that it is related with set variable. <br>
I wonder whether Nick can show his tcl file.<br>
<br>
Kijun<br>
<br>
<br>
----- Original Message -----<br>
From: "Stepan Stepanyan" <<a href="mailto:stepanya@jlab.org">stepanya@jlab.org</a>><br>
To: "clas offline" <<a href="mailto:clas_offline@jlab.org">clas_offline@jlab.org</a>><br>
Sent: Friday, May 16, 2014 8:30:10 AM<br>
Subject: Re: [Clas_offline] Generated events in the PART bank (PART bank sector 0) after user_ana.<br>
<br>
<br>
Johann, <br>
<br>
I am not familiar with how PART is handled in user_ana, it is in pid package I am sure.
<br>
But it is more likely that PART 0 is not deleted but it is just is not written out.
<br>
<br>
Stepan <br>
<br>
<br>
On 5/16/14 7:47 AM, Johann Goetz wrote: <br>
<br>
<br>
<br>
Does anyone know why user_ana should actively delete the PART bank sector 0 from the data? This sounds like a bug. In principle, gsim will digitize PART just the same it does for MCTK/VX. The low-level banks produced, and used by the reconstruction packages,
will be the same. <br>
<br>
<br>
In short, I believe the solution to Nick's problem is to NOT DROP the PART bank sector 0 during reconstruction (user_ana). Can someone tell him how to do this using user_ana?
<br>
<br>
<br>
-- Johann <br>
<br>
<br>
On Tue, May 13, 2014 at 3:47 PM, Gagik Gavalian < <a href="mailto:gavalian@jlab.org">
gavalian@jlab.org</a> > wrote: <br>
<br>
<br>
<br>
Hi Nick, <br>
<br>
<br>
The PART 0, is used when you’re running simulation and then reconstructing with a1c, which keeps both PART 0 (simulated bank) and PART 1 reconstructed bank.
<br>
If you’re using user_ana to for reconstruction, you should write your simulation into MCTK and MCVX banks, then run GSIM/RECSIS chain on it.
<br>
<br>
<br>
Gagik <br>
<br>
<br>
<br>
On May 13, 2014, at 3:38 PM, nick markov < <a href="mailto:markov@jlab.org">
markov@jlab.org</a> > wrote: <br>
<br>
<br>
user_ana <br>
<br>
_______________________________________________ <br>
Clas_offline mailing list <br>
<a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a>
<br>
<a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a>
<br>
<br>
<br>
<br>
_______________________________________________<br>
Clas_offline mailing list <a href="mailto:Clas_offline@jlab.org">
Clas_offline@jlab.org</a> <a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">
https://mailman.jlab.org/mailman/listinfo/clas_offline</a> <br>
<br>
_______________________________________________<br>
Clas_offline mailing list<br>
<a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a><br>
<br>
-- <br>
---------------------------------------******<br>
Dr. Kijun Park<br>
Physics Dept. Old Dominion University<br>
Thomas Jefferson National Accelerator Facility<br>
Suit 5 Room 12_1/A102<br>
12000 Jefferson Ave.<br>
Newport News VA 23606<br>
Phone <a href="tel:757-269-6989">
757-269-6989</a><br>
<br>
<br>
_______________________________________________<br>
Clas_offline mailing list<br>
<a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a><br>
</font></blockquote>
</div>
<font color="#000000" face="arial, helvetica,
sans-serif"><br>
</font></div>
<font color="#000000" face="arial, helvetica,
sans-serif"><br>
</font><fieldset></fieldset> <font color="#000000" face="arial, helvetica,
sans-serif">
<br>
</font>
<pre><font color="#000000" face="arial, helvetica, sans-serif">_______________________________________________
Clas_offline mailing list
<a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a>
<a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a></font></pre>
</blockquote>
<font color="#000000" face="arial, helvetica, sans-serif"><br>
</font></div>
</div>
</div>
<font color="#000000" face="arial, helvetica, sans-serif"><br>
_______________________________________________<br>
Clas_offline mailing list<br>
<a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a><br>
</font></div>
<br>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset> <br>
<pre wrap="">_______________________________________________
Clas_offline mailing list
<a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a>
<a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a></pre>
</blockquote>
<br>
</div>
</blockquote>
<blockquote type="cite">
<div>_______________________________________________<br>
Clas_offline mailing list<br>
<a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a></div>
</blockquote>
</body>
</html>