<div dir="ltr"><div><div><div><div>Hi Stepan, Hi all,<br><br></div>I managed to run successfully user_ana in a batch farm after implementing some <br>changes in one of recsis routines. There is still some issues with for example a Fortran <br>
"IF" statement which is not acting as it should be, and other things which need more investigation. However, I am thinking if at least the RTPC response is similar to what <br>we have in Pass1_v1, and if that what only matters in the drift path calibration, can I <br>
go ahead with processing the data for this purpose and don't hold this work while doing <br>the other tests? I may need the RTPC experts to give me their words on that after <br>having a look to the outputs of this batch cooking located at:<br>
/volatile/clas/claseg6/tgemRHL6_v1/1p2gev/<br></div>I processed in this test 10 files of a run 61448 with the RTPC calibration bosbanks. <br>Please ignore the outputs of a file 61448.A03 which crashed for a different issue!<br>
<br></div>Thank you,<br><br></div>Lamiaa<br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Dec 21, 2013 at 7:55 AM, Stepan Stepanyan <span dir="ltr"><<a href="mailto:stepanya@jlab.org" target="_blank">stepanya@jlab.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
Hello Nathan,<br>
<br>
Thanks for looking into this. Although we now managed to run the
code<br>
on 64-bin machines, there is still this issue with FORTRAN "IF"
statement.<br>
That could be one of sources of difference.<br>
<br>
But I agree with you, more cooking and comparison has to be done.<br>
<br>
Regards, Stepan <br><div><div class="h5">
<br>
On 12/20/13 8:10 PM, Nathan Baltzell wrote:
<blockquote type="cite">
<div>Dear EG6ers,<br>
</div>
<div><br>
</div>
<div>Here's an h10 comparison of this new 64bit-compiled cooking
versus pass1v1.</div>
<div>The only significant differences I noticed are trigger time
and neutral yields. </div>
<div><br>
</div>
<div>Similarities:</div>
<div>- same number of reconstructed events </div>
<div>- same number of good TPC tracks</div>
<div>- same number of IC clusters</div>
<div><br>
</div>
<div>Differences: (only in CLAS)</div>
<div>- few-% more charged tracks for 64bit</div>
<div>- few-% more electrons and protons </div>
<div>- 40% more neutrons, but same number of gammas (!?)</div>
<div><br>
</div>
<div>Looking a little closer, event-by-event:</div>
<div>- trigger time is ALWAYS different (by hundreds of ps up to
~1 ns) (!?)</div>
<div>- beta for neutrals differs due to trigger time</div>
<div>
<div>
<div>- reconstructed SC time, EC energy, EC time, CC nphe are
identical</div>
<div>- DC momentum is very similar (discrepancy is usually
sub-MeV)</div>
<div>- reconstructed TPC and IC banks are practically
identical</div>
<div><br>
</div>
<div>Note this is based on an early 1.2 GeV file, maybe we
should check a 6 GeV run too.</div>
<div><br>
</div>
<div>Regards,</div>
</div>
</div>
<div>Nathan</div>
<div><br>
</div>
<div>
<div><br>
</div>
<div><br>
</div>
<div>P.S. The files I compared are here:</div>
<div>
<div>
<div><br>
</div>
<div>new 64bit:</div>
<div>/w/hallb/claseg6/top_dir/test_64compil/ntp_61448_05.hbook</div>
<div>/w/hallb/claseg6/baltzell/test_64compil/hroot_61448_05.root<br>
</div>
</div>
</div>
<div><br>
</div>
<div>pass1v1:</div>
<div>/cache/mss/clas/eg6/production/pass1_v1/1p2gev/PAW/ntp_61448_05_pass1v1.hbook</div>
<div>/cache/mss/clas/eg6/production/pass1_v1/1p2gev/HROOT/hroot_61448_05_pass1v1.root</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
</div>
On Thu, 19 Dec 2013 20:57:19 -0600, Lamiaa El Fassi
<a href="mailto:lamiomar@gmail.com" target="_blank"><lamiomar@gmail.com></a> wrote:<br>
<br>
<blockquote style="margin:0 0 0.80ex;border-left:#0000ff 2px solid;padding-left:1ex">
<div dir="ltr">
<div>
<div>
<div>
<div>Hi Stepan, Hi All,<br>
<br>
</div>
I succeeded today to build the appropriate claseg6 32
bit libraries and binaries in the <br>
new 64 operating machine. Many thanks to Hovanes, who
showed me how to do that!<br>
<br>
</div>
The new user_ana executable under ~claseg6/bin/LinuxRHEL6/
runs successfully in the interactive mode. From the
printout below, it seems that all HBT and TBT conditions
are fulfilled:<br>
........<br>
TRK_STATI L: Cut# Description Value
Passed % <br>
TRK_STATI L: 1 trk_evnt: Ntrig (Nevt) >
0.0 50000 100.0 <br>
TRK_STATI L: 2 ana_segm: Nevt w. segm = 5.0
49645 99.3 <br>
TRK_STATI L: 3 trk_link: Nevt w. link = 5.0
47862 95.7 100.0 <br>
TRK_STATI L: 4 ana_prfit: Nevt PR chi2 <
70.0 47119 94.2 98.4 <br>
TRK_STATI L: 5 SEB_eid: Nevt SEB = 0.0
46394 92.8 96.9 <br>
TRK_STATI L: 6 trk_dcam: Nevt w. LRseg= 4.0
45968 91.9 96.0 <br>
TRK_STATI L: 7 trk_dcam: Nevt w. LRseg= 5.0
44345 88.7 92.7 <br>
TRK_STATI L: 8 ana_fit: Nevt FF1 chi2 <
70.0 43367 86.7 90.6 <br>
TRK_STATI L: 9 ana_fit: Nevt in TBT >
0.0 43367 86.7 90.6 <br>
..............<br>
</div>
<div>The outputs of this test run, in case anyone want to
have a look at them, are saved <br>
under this directory:
/w/hallb/claseg6/top_dir/test_64compil. Files produced
with RTPC bosbanks and 50K events.<br>
</div>
<div><br>
On the other hand, running user_ana on batch mode is still
giving the I/O errors while reading the tcl flags on
"retcl_yield.F" routine. I am attaching an example of the
two farm <br>
output files; example.err & example.txt. <br>
</div>
<br>
In order to know if this error is happening at the end of a
tcl file or while reading it, I <br>
put some write/print statements around the line causing this
issue, as follow:<br>
<br>
</div>
<div>-----------------------------------retcl_yield.F----------------------------------------------<br>
</div>
<div>...........................<br>
if (nbyt .gt. 0) then<br>
write(6,*) "RETCL YIELD READING FILE" ! Added 1st
write<br>
read(*,13, end=999) cline
! line 101<br>
print *,
cline ! Added
print<br>
write(6,*) "END RETCL YIELD READING FILE" !
Added 2nd write<br>
13 format(a200)<br>
cline = spaces(cline,1)<br>
iend = lenocc(cline)<br>
.......<br>
</div>
999 continue<br>
............<br>
<div> --------------------------------------------------------------------------------------------------------<br>
<br>
As you can see in .txt file, only the 1st write statement
was printed out (twice), and <br>
nothing after reading the tcl flag!<br>
<br>
</div>
<div>Any comments/suggestions on how to deal with this issue
will be appreciated!<br>
<br>
</div>
<div>Best regards,<br>
<br>
</div>
<div>Lamiaa<br>
<br>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Thu, Dec 19, 2013 at 5:34 PM,
Stepan Stepanyan <span dir="ltr"><<a href="mailto:stepanya@jlab.org" target="_blank">stepanya@jlab.org</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello
everyone,<br>
<br>
With Hovanes's help, I managed to compile user_ana on
ifarm1102 with<br>
-m32 flag.<br>
More, it runs and produced reasonable results. My build
with all libs is in<br>
/home/stepanya/eg6/test<br>
Appropriate files in cms directory have been modified to
take into<br>
account -m32<br>
and other related changes.<br>
<br>
Have a nice holidays, Stepan<br>
<br>
_______________________________________________<br>
Eg6_analysis mailing list<br>
<a href="mailto:Eg6_analysis@jlab.org" target="_blank">Eg6_analysis@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/eg6_analysis" target="_blank">https://mailman.jlab.org/mailman/listinfo/eg6_analysis</a><br>
</blockquote>
</div>
<br>
</div>
</div>
</blockquote>
<br>
<br>
<br>
</blockquote>
</div></div></div>
</blockquote></div><br></div>