opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

RE: [Opal] Electric field map in BANDRF cyclotron

From: Daniel Winklehner <winklehn AT mit.edu>
To: Zhang Hui <hui.zhang AT psi.ch>, Dirk Vandeplassche <dirkvdpla AT skynet.be>, "opal AT lists.psi.ch" <Opal AT lists.psi.ch>
Cc: "opal-developers AT lists.psi.ch" <opal-developers AT lists.psi.ch>
Subject: RE: [Opal] Electric field map in BANDRF cyclotron
Date: Mon, 16 Jun 2025 19:06:41 +0000
Accept-language: en-US
Arc-authentication-results: i=3; mx.microsoft.com 1; spf=pass (sender ip is 2a01:111:f403:c112::3) smtp.rcpttodomain=lists.psi.ch smtp.mailfrom=mit.edu; dmarc=pass (p=none sp=none pct=100) action=none header.from=mit.edu; dkim=pass (signature was verified) header.d=mit.edu; arc=pass (0 oda=1 ltdi=1 spf=[1,1,smtp.mailfrom=mit.edu] dkim=[1,1,header.d=mit.edu] dmarc=[1,1,header.from=mit.edu])
Arc-authentication-results: i=2; mx.microsoft.com 1; spf=pass (sender ip is 2a01:111:f403:c112::3) smtp.rcpttodomain=lists.psi.ch smtp.mailfrom=mit.edu; dmarc=pass (p=none sp=none pct=100) action=none header.from=mit.edu; dkim=pass (signature was verified) header.d=mit.edu; arc=pass (0 oda=1 ltdi=1 spf=[1,1,smtp.mailfrom=mit.edu] dkim=[1,1,header.d=mit.edu] dmarc=[1,1,header.from=mit.edu])
Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass smtp.mailfrom=mit.edu; dmarc=pass action=none header.from=mit.edu; dkim=pass header.d=mit.edu; arc=none
Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector10001; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-AntiSpam-MessageData-ChunkCount:X-MS-Exchange-AntiSpam-MessageData-0:X-MS-Exchange-AntiSpam-MessageData-1; bh=VfbRZKdTXoLziMeObT6srpE4mcSZIAv16303JT34r94=; b=JE7h+wwkVWG9qK3Zft9VkzQpBut6IPwyhCwh/1796V1C2J9xxGCW/UJykNeOhD0TVNJbAsqLCoJi0pYjMKEMejtcIxk2TTML8VH+AJgEQpqDPVPvokzEwXkSQYJ84OVqis2bh/XdOicqSFGBI5R22DZhSNKAotMuR71evF5lg7bJ1q2UwPzdMm5mRRMluMUuuePqRfJsLhDtRRNyZUg9tcK/gclAS7f4zr1oWFvcyA7j576LeNMhbPKwysYVYITmXP1uOXZpIn4zhVZpw5x1mS28FQNyqzLYNd1CikvbQEM42frbl4Z3Vp4atOG0yTwzfEj4gniGCHHlzWgORlgevA==
Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector10001; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-AntiSpam-MessageData-ChunkCount:X-MS-Exchange-AntiSpam-MessageData-0:X-MS-Exchange-AntiSpam-MessageData-1; bh=VfbRZKdTXoLziMeObT6srpE4mcSZIAv16303JT34r94=; b=ZPE/VoqcxxVieoKQPDDsJGHjPFHG9BuJGv8DxAJRh7bhk+pYeoVJkecBncWk1u0vogT36W3xyqBJzD7GhSX7xNQDvWSKqpVJngwFCf2wrYsi3Ck3fEH5dPvayCaqKjUkkgln9EGaGChk3dEnu1D4RnLzbhvjFbgwJW9buCW9S2F3+Bfy1HiAgx97CRwXm+Xy8v01NOt29FgZelmXGhSlZrE3bDNOwslk8Z0dFhfifO1WoXopqzwbdyNSBhufSGDsZZv7lP2r+klHq3BWGLpW5GJcov/+5sxxvYXldPR4nGQYkj04arIiG2djuTvFVWqk40Z2wZcTvQotu3cwPiVFrA==
Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector10001; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-AntiSpam-MessageData-ChunkCount:X-MS-Exchange-AntiSpam-MessageData-0:X-MS-Exchange-AntiSpam-MessageData-1; bh=VfbRZKdTXoLziMeObT6srpE4mcSZIAv16303JT34r94=; b=QoeFPA9/duQ327P9TS//59qt//z145w8v3nCYyrM/rKg78mzW6AnHSnP0qqn4eOoZly33XxaMlYw7yZ0Zuuc12m1NcdVcZXeAX232KIfy9XPiXwpaFb5XwcYsbaHbiAzJlqcota1Czl2AZereZN6qMHKYFciKgUZg9Coehfy5RCRwSkF7z7CZ9Cciamw04R4+arj0KUaPElRX2OiVPpsrfvkAUaOSPMm0X4ElZKCT3TXPeOGIsbkd7dJbdI6GUDlMcgJZ0SGSS2P5AQbx3R79IdhnxY+FGOnX7TpIM8qBbCbm0MOU5rR1YGn4dxnj5JpH7IUGUtFbAI1954mJ9nd1g==
Arc-seal: i=3; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=pass; b=KfKfKdBByEwlCe9xS7owlJr80w7raFnSykQ0QshQyaCI+jgZ72JD+VO8j/EDbJFZf5jxOTLUR9KEg7JMNkH4icdg0qkuzcy1QE2fg88KcvQCTYVdPmbCkXaGBFFoFM+IdJ+5walJpJiimWfA+LNbWXxoLL9HK7NSqWPqGzhZWDONFkcFPiK4OuLvt2rbfK5+JkV5aEJVieo1Ah2Bxdf6GjlbxoPLYPaV+rTHZnneZUZar2eQIoyGip2xGQwls7F6WnmCop9izkn7tSbv2lwy2BNnbjOR3v87mS40LNhKtNa9dw5As64VimXUQsTWb16VQrS+yLDyfhNBzQf2rQRs7w==
Arc-seal: i=2; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=pass; b=u8yH0/UgKLS9owJIkqKFNLpUp0btX/auYx8mkl3ZW1PwUaWXN3gaVwlwJ4LgOAYqtODU4p2aMJPlNQnr+P02tNZDkSyBHqlwBi7bVJ+EwT/RYBwng1a++L2E4W3HzqqWzNingMGZpQzeRicKbpKLal+lMKqQSHtwu0VKBp5BuIhybqyzZmd6xcNX8hRB4DT3Ayz8DZPoJlPDiEEjKqYpvsp16SBzOBELnwetxFFWMrsP0s2XZSoVTl9QZYCBGvDJ1YlYhZiS8bIOxotOBisBVzX0nItX4h9NO2BTJQwQ9CVFEEN0EderCZmiBxauSsh3eHymYbaOBueLqJTCi/31CQ==
Arc-seal: i=1; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=none; b=GT5SOzMwa0W1ilB+JTsv5/eY+FrX3FpMpjBnT3rFJNiMGjWZ8uP8Acfe4WZmLzk+ffsBHBrv4ot1xJ3ZrajQCZxiHRKT/YkvOOYUAfScQHSQ1JUVlWf+OGEEqbFdo+Li01En8kvPNWYJViA2AdMVPCPqeSbJjZfw/1g4v2rhCy1RiGNwWOLUAF9QUBRLDbU7y90dtzFhkD5VO7VDsio7r8MiDvEUZqs9cEliKCHF7lgSYqz4YquJEaa1nM+HwTLevPDWgnHVjRmBdP0phR0bGylJk3Ncq92IOpNUHx7gdeB+clIJs6ecGYVjY9ieix36ZU3flXu9JaOuCCPgDiOS2Q==
Authentication-results: spf=pass (sender IP is 2a01:111:f403:c112::3) smtp.mailfrom=mit.edu; dkim=pass (signature was verified) header.d=mit.edu;dmarc=pass action=none header.from=mit.edu;compauth=pass reason=100
Authentication-results-original: dkim=none (message not signed) header.d=none;dmarc=none action=none header.from=mit.edu;

Hi Dirk,

I have attached a python script that I wrote to replace the C++ program. It may not cover all pertinent cases, but it should give you a starting point. I zipped it because some email providers do not allow sending .py files.

Best,

Daniel

From: opal-request AT lists.psi.ch <opal-request AT lists.psi.ch> On Behalf Of Zhang Hui
Sent: Monday, June 16, 2025 16:12
To: Dirk Vandeplassche <dirkvdpla AT skynet.be>; opal AT lists.psi.ch
Cc: opal-developers AT lists.psi.ch
Subject: RE: [Opal] Electric field map in BANDRF cyclotron

Dear Dirk,

I was transforming ascii file to *.h5part map four years ago. The format of my ascii file looks like the following

nx ny nz

x1 y1 z1 Ex1 Ey1 Ez1

x1 y1 z2 Ex2 Ey2 Ez2

x1 y1 z3 Ex3 Ey3 Ez3

…

The original c++ code doing the transformation must be from Andreas. But I am not sure about that.

I would think, that the head could be different. However, the numbers of grids in x, y, y-directions have to be defined in head. And this information has to be correctly transferred into .h5part file.

The .h5part map can be visualized by a program called VisIt. If the .h5part file cannot be opened, or it looks strange, then something must be wrong, in the ascii data, and/or during the transformation.

Best regards,

Hui

__________________________________________
Paul Scherrer Institut
Hui Zhang
WBGA/C34
Forschungsstrasse 111
5232 Villigen PSI
Schweiz

Telefon: +41563105398
E-Mail: hui.zhang AT psi.ch

From: opal-request AT lists.psi.ch <opal-request AT lists.psi.ch> On Behalf Of Dirk Vandeplassche
Sent: Monday, June 16, 2025 1:05 PM
To: opal AT lists.psi.ch
Cc: opal-developers AT lists.psi.ch
Subject: [Opal] Electric field map in BANDRF cyclotron

Dear all,

It would really be very helpful for me to obtain a detailed description of the ASCII-version of the electric field map to be used in a BANDRF-type cyclotron. I am especially interested in the header of this file, hoping that it may solve the error in OPAL "Couldn't determine type of fieldmap in file “efield.h5part””. The conversion from the ASCII file to the binary file with the “ascii2h5block” tool seems to work well — anyway it does not give any error message and returns a correct number of data with a header like

101

0.0005 0.0005 0.0001

-0.025 -0.025 -0.003 502.1066965 342.215468 -1294.568941

…

Its output is

Combine efield2.txt and to efield2.h5part

Frequency 7.2615e+07 [Hz]

H5block grid

H5gridPx 101

H5gridPy 101

H5gridPz 61

number Edata 622261

gridPx 101 stepEx 0.0005

gridPy 101 stepEy 0.0005

gridPz 61 stepEz 0.0001

number Bdata 0

Done bye ...

Any hint that may solve the problem is welcome !

Many thanks in advance,

Dirk Vandeplassche

Attachment: ascii2h5block.zip
Description: ascii2h5block.zip

[Opal] Electric field map in BANDRF cyclotron, Dirk Vandeplassche, 06/16/2025
- RE: [Opal] Electric field map in BANDRF cyclotron, Zhang Hui, 06/16/2025
  - RE: [Opal] Electric field map in BANDRF cyclotron, Daniel Winklehner, 06/16/2025