From j.voigt at fz-juelich.de Mon Jan 23 09:19:51 2006 From: j.voigt at fz-juelich.de (=?iso-8859-1?q?J=F6rg_Voigt?=) Date: Mon, 23 Jan 2006 09:19:51 +0100 Subject: [neutron-mc] Re: TbMnO3 In-Reply-To: <54718.193.49.43.124.1137577700.squirrel@195.83.126.3> References: <200601171508.19212.j.voigt@fz-juelich.de> <54718.193.49.43.124.1137577700.squirrel@195.83.126.3> Message-ID: <200601230919.52147.j.voigt@fz-juelich.de> Dear Tapan, could you please send me a copy of the proposal for my files? Thank you in advance J?rg Am Wednesday 18 January 2006 10:48 schrieb chatt at ill.fr: > Dear Joerg, > You are very welcome to our TbMnO3 project. I have included > your name as a proposer. It would be nice if you take care of > the crystal. We would like to have a plate in the a-b plane (Pbnm) > with c axis perpendicular to the plane of the plate. The dimensions > can be of the order of 5x5x2 mm^3 if possible. > I am now in India but shall read and respond to messages time to > time if not regularly. > Best regards. > Tapan > > > Dear Tapan, > > Thomas told me, that you and Jane Brown proposed a polarimetry study of > > TbMnO3. I would be happy, if I could participate in this experiment. > > We could contribute a single crystal of TbMnO3, our colleague J?rg Person > > has > > grown a first test crystal and will soon make a larger specimen. We would > > also take care of preparation of the crystal. > > Perhaps we can discuss details via mail. > > > > regards > > > > J?rg > > -- > > =====================Dr.J?rg Voigt > > Institut f?r Streumethoden > > Forschungszentrum J?lich > > 52424 J?lich > > Deutschland > > Tel: +49 2461 616020 > > Fax: +49 2461 612610 > > Mobil: +49 160 8479827 > > Email: j.voigt at fz-juelich.de > > ===================== -- ====================== Dr.J?rg Voigt Institut f?r Streumethoden Forschungszentrum J?lich 52424 J?lich Deutschland Tel: +49 2461 616020 Fax: +49 2461 612610 Mobil: +49 160 8479827 Email: j.voigt at fz-juelich.de ====================== From j.voigt at fz-juelich.de Mon Jan 23 09:38:56 2006 From: j.voigt at fz-juelich.de (=?iso-8859-1?q?J=F6rg_Voigt?=) Date: Mon, 23 Jan 2006 09:38:56 +0100 Subject: [neutron-mc] Re: TbMnO3 In-Reply-To: <200601230919.52147.j.voigt@fz-juelich.de> References: <200601171508.19212.j.voigt@fz-juelich.de> <54718.193.49.43.124.1137577700.squirrel@195.83.126.3> <200601230919.52147.j.voigt@fz-juelich.de> Message-ID: <200601230938.56276.j.voigt@fz-juelich.de> Sorry, I don' know why the mailing list appeared in the address of my last mail Regards J?rg Am Monday 23 January 2006 09:19 schrieb J?rg Voigt: > Dear Tapan, > could you please send me a copy of the proposal for my files? > > Thank you in advance > > J?rg > > Am Wednesday 18 January 2006 10:48 schrieb chatt at ill.fr: > > Dear Joerg, > > You are very welcome to our TbMnO3 project. I have included > > your name as a proposer. It would be nice if you take care of > > the crystal. We would like to have a plate in the a-b plane (Pbnm) > > with c axis perpendicular to the plane of the plate. The dimensions > > can be of the order of 5x5x2 mm^3 if possible. > > I am now in India but shall read and respond to messages time to > > time if not regularly. > > Best regards. > > Tapan > > > > > Dear Tapan, > > > Thomas told me, that you and Jane Brown proposed a polarimetry study of > > > TbMnO3. I would be happy, if I could participate in this experiment. > > > We could contribute a single crystal of TbMnO3, our colleague J?rg > > > Person has > > > grown a first test crystal and will soon make a larger specimen. We > > > would also take care of preparation of the crystal. > > > Perhaps we can discuss details via mail. > > > > > > regards > > > > > > J?rg > > > -- > > > =====================Dr.J?rg Voigt > > > Institut f?r Streumethoden > > > Forschungszentrum J?lich > > > 52424 J?lich > > > Deutschland > > > Tel: +49 2461 616020 > > > Fax: +49 2461 612610 > > > Mobil: +49 160 8479827 > > > Email: j.voigt at fz-juelich.de > > > ===================== -- ====================== Dr.J?rg Voigt Institut f?r Streumethoden Forschungszentrum J?lich 52424 J?lich Deutschland Tel: +49 2461 616020 Fax: +49 2461 612610 Mobil: +49 160 8479827 Email: j.voigt at fz-juelich.de ====================== From zilly at hmi.de Tue Feb 14 16:12:08 2006 From: zilly at hmi.de (Frederik Zilly) Date: Tue, 14 Feb 2006 16:12:08 +0100 Subject: [neutron-mc] Error in component Source_simple.comp ? Message-ID: <43F1F348.4010205@hmi.de> Dear Mcstas experts, The point which emits the neutron is calculated by: chi=2*PI*rand01(); /* Choose point on source */ r=sqrt(rand01())*radius; /* with uniform distribution. */ x=r*cos(chi); y=r*sin(chi); In my humble opinion, this will not lead to an uniform distribution over the sources surface. Points near the center have greater probility than points with a higher radius. Thanks, Frederik. From peter.willendrup at risoe.dk Wed Feb 15 13:41:26 2006 From: peter.willendrup at risoe.dk (Peter Willendrup) Date: Wed, 15 Feb 2006 13:41:26 +0100 Subject: [neutron-mc] Error in component Source_simple.comp ? In-Reply-To: <43F1F348.4010205@hmi.de> References: <43F1F348.4010205@hmi.de> Message-ID: <1140007287.9137.33.camel@localhost.localdomain> Hello Frederik, On Tue, 2006-02-14 at 16:12 +0100, Frederik Zilly wrote: > Points near the center have greater probility than points with a higher > radius. Actually, this is not the case. The square root actually gives zero probability for points exactly at the origin; more events must be distributed over a circular area element at the circle perimeter than at the center to achieve a uniform distribution. I wrote up a few lines of matlab code to show this in graphics (you need the attached hist2d code to reproduce the attached graphics): chi=2*pi*rand(1e6,1); r=sqrt(rand(1e6,1)); % Source radius assumed one metre x=r.*cos(chi); y=r.*sin(chi); mXY=[x y]; vXEdge = linspace(-1,1,128); vYEdge = linspace(-1,1,128); img=hist2d(mXY,vXEdge,vYEdge); imagesc(img); As you can see in the jpg graphic, the result of the above code (which is equvalent to the c-code in Source_simple) is completely uniform across the source area. Hope this helps? Cheers, Peter -- Peter Willendrup Ris? National Laboratory -------------- next part -------------- A non-text attachment was scrubbed... Name: hist2d.m Type: text/x-objcsrc Size: 1652 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: uniform.jpg Type: image/jpeg Size: 64596 bytes Desc: not available URL: From zilly at hmi.de Wed Feb 15 16:01:04 2006 From: zilly at hmi.de (Frederik Zilly) Date: Wed, 15 Feb 2006 16:01:04 +0100 Subject: [neutron-mc] Error in component Source_simple.comp ? In-Reply-To: <1140007287.9137.33.camel@localhost.localdomain> References: <43F1F348.4010205@hmi.de> <1140007287.9137.33.camel@localhost.localdomain> Message-ID: <43F34230.5040100@hmi.de> Hello Peter ! Thank you very much for your quick answer. Indeed, I must have been blind not seeing the square root. With this square root everything is alright as your jpg shows. Sorry for this "false alarm". Cheers, Frederik. Peter Willendrup wrote: > Hello Frederik, > > On Tue, 2006-02-14 at 16:12 +0100, Frederik Zilly wrote: > >>Points near the center have greater probility than points with a higher >>radius. > > > Actually, this is not the case. The square root actually gives zero > probability for points exactly at the origin; more events must be > distributed over a circular area element at the circle perimeter than > at the center to achieve a uniform distribution. > > I wrote up a few lines of matlab code to show this in graphics (you need > the attached hist2d code to reproduce the attached graphics): > > chi=2*pi*rand(1e6,1); > r=sqrt(rand(1e6,1)); % Source radius assumed one metre > x=r.*cos(chi); > y=r.*sin(chi); > mXY=[x y]; > vXEdge = linspace(-1,1,128); > vYEdge = linspace(-1,1,128); > img=hist2d(mXY,vXEdge,vYEdge); > imagesc(img); > > As you can see in the jpg graphic, the result of the above code (which > is equvalent to the c-code in Source_simple) is completely uniform > across the source area. > > Hope this helps? > > > Cheers, > > Peter > > > ------------------------------------------------------------------------ > > %function mHist = hist2d ([vY, vX], vYEdge, vXEdge) > %2 Dimensional Histogram > %Counts number of points in the bins defined by vYEdge, vXEdge. > %size(vX) == size(vY) == [n,1] > %size(mHist) == [length(vYEdge) -1, length(vXEdge) -1] > % > %EXAMPLE > % mYX = rand(100,2); > % vXEdge = linspace(0,1,10); > % vYEdge = linspace(0,1,20); > % mHist2d = hist2d(mYX,vYEdge,vXEdge); > % > % nXBins = length(vXEdge); > % nYBins = length(vYEdge); > % vXLabel = 0.5*(vXEdge(1:(nXBins-1))+vXEdge(2:nXBins)); > % vYLabel = 0.5*(vYEdge(1:(nYBins-1))+vYEdge(2:nYBins)); > % pcolor(vXLabel, vYLabel,mHist2d); colorbar > function mHist = hist2d (mX, vYEdge, vXEdge) > [nRow, nCol] = size(mX); > if nCol < 2 > error ('mX has less than two columns') > end > > nRow = length (vYEdge)-1; > nCol = length (vXEdge)-1; > > vRow = mX(:,1); > vCol = mX(:,2); > > mHist = zeros(nRow,nCol); > > for iRow = 1:nRow > rRowLB = vYEdge(iRow); > rRowUB = vYEdge(iRow+1); > > [mIdxRow] = find (vRow > rRowLB & vRow <= rRowUB); > vColFound = vCol(mIdxRow); > > if (~isempty(vColFound)) > > > vFound = histc (vColFound, vXEdge); > > nFound = (length(vFound)-1); > > if (nFound ~= nCol) > [nFound nCol] > error ('hist2d error: Size Error') > end > > [nRowFound, nColFound] = size (vFound); > > nRowFound = nRowFound - 1; > nColFound = nColFound - 1; > > if nRowFound == nCol > mHist(iRow, :)= vFound(1:nFound)'; > elseif nColFound == nCol > mHist(iRow, :)= vFound(1:nFound); > else > error ('hist2d error: Size Error') > end > end > > end > > > > > ------------------------------------------------------------------------ > > > ------------------------------------------------------------------------ > > _______________________________________________ > neutron-mc mailing list > neutron-mc at risoe.dk > http://mailman.risoe.dk/mailman/listinfo/neutron-mc From peter.willendrup at risoe.dk Sun Mar 12 19:36:11 2006 From: peter.willendrup at risoe.dk (Peter Willendrup) Date: Sun, 12 Mar 2006 19:36:11 +0100 Subject: [neutron-mc] Test message, please ignore Message-ID: <1142188571.27571.15.camel@localhost.localdomain> Hello everyone? The mailing list server at RIS? has been upgraded, this is a test message to check all is working again. Regards, Peter -- Peter Willendrup Ris? National Laboratory From peter at willendrup.org Mon Mar 27 15:34:13 2006 From: peter at willendrup.org (Peter Willendrup) Date: Mon, 27 Mar 2006 15:34:13 +0200 Subject: [neutron-mc] Bug and fix: mcgui interface crashing with Gtk2 apps running Message-ID: <1143466454.8497.83.camel@localhost.localdomain> Hello everyone, NOTE: The problem below applies to Unix systems with gnome only. Recent updates to Gtk (gnome library) has made the mcgui graphical interface crash on startup/use of certain gnome applications. If you experience this problem, please either 1) Install our locally patched version of perl-Tk: http://www.mcstas.org/download/Tk-804.027_gtk2_patch.tar.gz (Source code package, requires you to do perl Makefile.PL && make && make install) 2) Apply the patch yourself from the patchfile available at http://www.mcstas.org/download/gtk2-patch Regards, Peter Willendrup Ps.: mcstas-1.9.1 with minor updates and bugfixes will be released very soon. -- Peter Willendrup From peter.willendrup at risoe.dk Wed Mar 29 16:11:58 2006 From: peter.willendrup at risoe.dk (Peter Willendrup) Date: Wed, 29 Mar 2006 16:11:58 +0200 Subject: [neutron-mc] Announcing McStas 1.9.1 Message-ID: <1143641519.1509.39.camel@localhost.localdomain> Hello everyone, The McStas (http://www.mcstas.org) team is happy to announce that an update release of the general neutron ray-trace simulation package is now available. The new release can be downloaded from the project website at http://www.mcstas.org/download For a list of the most important new bugfixes, additions and features, see the attached CHANGES.txt document. Regards from the McStas team: Peter Willendrup - peter.willendrup at risoe.dk Peter Christiansen - peter.christiansen at risoe.dk Kim Lefmann - kim.lefmann at risoe.dk Emmanuel Farhi - farhi at ill.fr Klaus Lieutenant - lieutena at ill.fr -- Peter Willendrup Ris? National Laboratory -------------- next part -------------- /******************************************************************************* * * McStas, neutron ray-tracing package * Copyright (C) 1997-2006, All rights reserved * Risoe National Laboratory, Roskilde, Denmark * Institut Laue Langevin, Grenoble, France * * Documentation: CHANGES * * %Identification * Written by: KN, KL, PEO, EF, PW * Date: 1997 * Origin: Risoe, ILL, ANSTO * Release: McStas 1.8 * Version: 1.8 * Modified by: KN, October 26, 1998 : initial release 1.0 * Modified by: KN, March 31, 1999 : release 1.1 * Modified by: KN, January 31, 2000 : release 1.2 * Modified by: KN, May 18, 2000 : release 1.3 * Modified by: KN, July 28, 2000 : release 1.4 * Modified by: KN, PEO, March 16, 2001: release 1.4.1 * Modified by: PEO, EF, October, 10th, 2001: release 1.5 * Modified by: PW, EF, May 19th 2003: version 1.7 * Modified by: PW, EF, March 4th 2004: version 1.8 * Modified by: PW, EF, November 16th 2005: version 1.9 * * This file is part of McStas version 1.9, released November, 16th, 2005. * It gives a 'changes' list from the beginning of the project * *******************************************************************************/ Changes in McStas v1.9.1, Mar 29, 2006 General - Update release, mainly bugfixes plus a few other things Documentation - Component manual slightly updated according to adding/modification of components. Tools - mcgui: A few bugfixes. Note also that if you experience crashes when starting/running GNOME/Gtk2 applications, update your perl-Tk installation. Two ways to do this: 1) Get the tarball from http://www.mcstas.org/download/Tk-804.027_gtk2_patch.tar.gz and install using 'perl Makefile.PL && make && make install'. 2) Do a 'make tk' in the McStas unpacked McStas source directory. This will download and install the tarball into your McStas system directory. Components - DiskChopper.comp - new disc chopper component. Chopper.comp made obsolete. See mcdoc page / component manual for further details. - Guide.comp - reads reflectivity tables from disk. - Isotropic_Sqw - minor bugfixes, validated against PowerN, V_sample, etc. - PowderN - validated against Isotropix_Sqw, Squires, experimental data. WARNINGS - Because of rapid changes in the support software for McStas, the current release has support for only some releases of: * Scilab (Win32 and Unix) - pick release 4.0 if available or pick release 3.0 from , e.g. from http://www.scilab.org/download/index_download.php?page=oldReleases.html * Perl (Win32) - on Windows we only support Perl 5.6 - pick http://downloads.activestate.com/ActivePerl/Windows/5.6/ActivePerl-5.6.1.638-MSWin32-x86.msi * If you are running Unix with GNOME/Gtk2 applications, you may experience problems with mcgui crashing, please install updated perl-Tk as noted above.