xcoll 0.3.6__py3-none-any.whl → 0.4.0__py3-none-any.whl

{xcoll-0.3.6.dist-info → xcoll-0.4.0.dist-info}/RECORD

@@ -1,36 +1,41 @@
 LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
 NOTICE,sha256=6DO_E7WCdRKc42vUoVVBPGttvQi4mRt9fAcxj9u8zy8,74
-xcoll/__init__.py,sha256=XO53i_KKg6kaJfTCpmNJLy14tou4RMF9MDuXlH70iTQ,642
-xcoll/beam_elements/__init__.py,sha256=4ekQv_1Kxwqd0VSaPzMpVQ1V9RMX5tr-ycBBiTIFyRw,776
-xcoll/beam_elements/absorber.py,sha256=yxF-NhnQs_l7s-M_0AhdKzx0YVgvhujOAwQnAbVTSOs,1260
-xcoll/beam_elements/base.py,sha256=Nwgg2v1RdEEpByFG8AZgL7E6ecESbYm3mqnIfoh5SFQ,15785
-xcoll/beam_elements/collimators_src/absorber.h,sha256=ZsAzIZ3uWP15unMJyuxbuD5lVgFeWAwQsuyOj_ZSkjI,5502
-xcoll/beam_elements/collimators_src/everest_block.h,sha256=Gpr5wyFLtcaNgVj3auKgFWB-2KPP5cjbEWTfC3f9q1U,5596
-xcoll/beam_elements/collimators_src/everest_collimator.h,sha256=wsGIFBb-VMCRhxwK2G5k9SEX-6SzIgRX4XE5g_I2okg,6296
-xcoll/beam_elements/collimators_src/everest_crystal.h,sha256=pX9cmW5lfsUoszR14tXWU4Emu-0VwX4bbfg6Q6pxXR4,5164
-xcoll/beam_elements/everest.py,sha256=JxmREXln3p3c9XEcHp-DgDbXFzc9txiGrxvjOrDUacs,11369
-xcoll/colldb.py,sha256=_MywNPz9cDocGPsfXpBVkjADyAGGVFQQRUI4kxOQaOg,49476
-xcoll/collimator_settings.py,sha256=uj25cpSaM-EY7VM65OUq3Cxb8gedSmi-r9mpkP-ikrg,17015
-xcoll/general.py,sha256=Pcqa-FWVvvs3X5gBEhnNw2qmw-Rl9IdCMX2Mww33wt4,525
-xcoll/headers/checks.h,sha256=3CLIa2etrJVKBfxdsHcS8NirqzVDHG9aCbrlMwkWuHw,1586
-xcoll/headers/particle_states.h,sha256=YxpS56_ZD8mhrqKTohcyVe3hbQZ6H9vspYHSVM6a5-U,844
-xcoll/impacts/__init__.py,sha256=yKQdIeLG2CT3qo-iZB6MODm-WHzkTmulgUvOAS41w7E,39
-xcoll/impacts/impacts.py,sha256=J7FfcxIprUMHyn5oZhmjhvBewXOp7kwJtErStOsAwJQ,3953
-xcoll/impacts/impacts_src/impacts.h,sha256=0eq4Z8xDPMJzN8yPF3lTkDsld6W54tbXmynZEbXCh4c,5950
-xcoll/impacts/interaction_types.py,sha256=iHbhC4beN1vZRWeHbALqSFo8lG10R7asc08BRpIqXx0,2493
-xcoll/manager.py,sha256=A7hCuSgS7sSFKoV26Nt8xc_GXUNOm-WRqsGftRoBMzA,45052
-xcoll/rf_sweep.py,sha256=_cnDb2d0_tfVrrwuubLOSyUHpMQHMvmSzYgd_-WfJg0,8936
+xcoll/__init__.py,sha256=EN8qYpv6nXGTqJYF307Q8BivglN7YD_2JmZYQOvAceg,929
+xcoll/beam_elements/__init__.py,sha256=NiOjiV5t266rFfJAbSrNY07GDSXbD3fSrulRJe3CnwU,981
+xcoll/beam_elements/absorber.py,sha256=jDCSjaBUmflDVS4DlfimA-MJurFYaqVvUO8j0xXQzOQ,2508
+xcoll/beam_elements/base.py,sha256=fW9QmzsKohXH5VivU55NefmJ82FJiHAN1wbkniD05b4,49646
+xcoll/beam_elements/collimators_src/black_absorber.h,sha256=PBRqmr24jngZQ8UZlYsodGRT4018_k3qvWmvmMIpjFw,4478
+xcoll/beam_elements/collimators_src/black_crystal.h,sha256=l5_9QrQgsau-0PfEaYV9m394-JMOMm9MsJNkl4WFUM8,4082
+xcoll/beam_elements/collimators_src/everest_block.h,sha256=OhJ6M2qrqxCHUWRs9GgMOBr3E15pI2UyO_eHyqWGkSU,5970
+xcoll/beam_elements/collimators_src/everest_collimator.h,sha256=hnn1zOsLXcBpNtVdQWWJHxmWmt2BhSF0wuqU8s7bO7g,9857
+xcoll/beam_elements/collimators_src/everest_crystal.h,sha256=kaswVVR1a6STSaRihevk7J9XDDPevefAIpJLmCmk3k8,11822
+xcoll/beam_elements/everest.py,sha256=nuqnIn8lOERD-lbdk-21aeh_OjEswioLbtCxqosxIu0,8405
+xcoll/colldb.py,sha256=Eu-Vphi0FStulGCG9UG37dVULkybYwB9agEydjeLN-o,29822
+xcoll/general.py,sha256=r0XT40A_S2MEDtXD8L10drQoXlfvHJR3oHULY_x4wKQ,525
+xcoll/headers/checks.h,sha256=IMLhQMnh3t-t9UhVK-1hfS4GVgl2Mry69H_17By00M0,1582
+xcoll/headers/particle_states.h,sha256=WCzeHAniEK7sX7HTwhoIgewioQ7SibGn01RLrhVI98k,840
+xcoll/initial_distribution.py,sha256=Ibsyc8h6Sw4spaxjYVbKKuqQJUxXNhKci1eCmZEKK3M,7858
+xcoll/install.py,sha256=oJ2oWMJTA1GiaU6XZrimG2MFPOvGF0Hr5-SfHIZ2184,7552
+xcoll/interaction_record/__init__.py,sha256=UFoLiKa-z2oX7YoszP-7Vgdt1nM6kT382v1CaIu8_u0,50
+xcoll/interaction_record/interaction_record.py,sha256=HasJDWRVQnMbptWFg4ElKuV2shcdtiRsJGweWiuWqvU,10733
+xcoll/interaction_record/interaction_record_src/interaction_record.h,sha256=n3CSYiC_3hk52oAeNM9-Xg_aWXMNAwhi5E2ANxnXW28,5949
+xcoll/interaction_record/interaction_types.py,sha256=_UGq-PBq30hj1lp29CpNTmmwzsP7tuCng8ZIAALnzCc,2894
+xcoll/line_tools.py,sha256=vQi9fuUQOme2XkHEj47tqcFtL6jJ2-H06U4Gw7Vds64,3250
+xcoll/lossmap.py,sha256=uzZpSe0R0UzvxDxauK5X2AegJapc7Q-Eg3B8WmmR8Js,7486
+xcoll/manager.py,sha256=Uy58saX2836q_CZjNa5at611k38Za_L9jQSTqBQXooA,343
+xcoll/rf_sweep.py,sha256=5ABhou9cZcLYQ8DSPYowPpLCea8zJQyhWgUa7JUHJIg,8936
 xcoll/scattering_routines/everest/__init__.py,sha256=yuCzolZkiL7mHaBXCKHRl66VxjKFB0t9Vyw9TtFv1Bo,286
+xcoll/scattering_routines/everest/amorphous.h,sha256=9iYYsak-Z2UR1tCQjkSUs0LNPiFGkZ7nQnqXpbaZxqY,10424
+xcoll/scattering_routines/everest/channeling.h,sha256=I-YUqzqihgYOgWJyGmpXMspF_Yv9M2QB1lWToyFVqYA,11313
 xcoll/scattering_routines/everest/constants.h,sha256=9Kus-RKyxsGdlxx9tfMWUz2H5Q-qy6y_F4kddIZzO30,1390
-xcoll/scattering_routines/everest/crystal.h,sha256=WlA_e9ePuVhc9fmMVdSXFBKvHrFvCQ5aG4d67iMrtuQ,54119
-xcoll/scattering_routines/everest/everest.h,sha256=p4yfzxRvh19N-V1vorf_CCIBotodV40Lcrzzy7339WM,3743
-xcoll/scattering_routines/everest/everest.py,sha256=zfKZXHvXM6ikoI2KUZ0_BWLCuxLZEvP1pgW5vFIgwQ4,1306
-xcoll/scattering_routines/everest/jaw.h,sha256=lmijUU7ih4veifVho3ehGe0DVyIYkDHMVIIclDNbZds,7308
-xcoll/scattering_routines/everest/materials.py,sha256=9ZhBd6-n2AsTbCxLhjPoU4pOcmKstl6gbeqLno1dY4c,9410
-xcoll/scattering_routines/everest/multiple_coulomb_scattering.h,sha256=0Euxsnc2BU8YqxYwfXZIiplNcnDzFq9HgWAoqGwuSAg,4518
+xcoll/scattering_routines/everest/crystal_parameters.h,sha256=NjCYWbl3Gyv0OXKdeK7N3oisbMAkg4PnTmayGlRU9YI,5416
+xcoll/scattering_routines/everest/everest.h,sha256=rfVm0q1PQwLoEadF7wiliBSWzOcto4BLwoU1DaHA-Sk,2658
+xcoll/scattering_routines/everest/everest.py,sha256=oRxiuzCBvuEyxzxjOyV-AcgpqzMBa0iVUQqTyxkNUcw,1460
+xcoll/scattering_routines/everest/jaw.h,sha256=hLofF1skic08_HTydiwkH2Ylo2pN_rHS6C_QgmWNL5c,2270
+xcoll/scattering_routines/everest/materials.py,sha256=SaSx2jUImYkiwSgWCJj88ULeSD0y_pJhCgkck8omtTY,9524
+xcoll/scattering_routines/everest/multiple_coulomb_scattering.h,sha256=HkS1Yjxqnb7mjhkfkj2zY325F1K9RWCprnZ4sEcSbLY,5400
+xcoll/scattering_routines/everest/nuclear_interaction.h,sha256=Jr3uSDwGJNnsJIDr72JWbMsg4Zbm_wXW9180K5u9duA,3113
 xcoll/scattering_routines/everest/properties.h,sha256=aRqkX5ZpUeTnbygqFVsoIKwWyoB5VXzD5vq5d-xiI74,6585
-xcoll/scattering_routines/everest/scatter.h,sha256=jl-SAIpizYCbUfy1KAhuO00sA5BBhNBCIjueQ8bCxoQ,6333
-xcoll/scattering_routines/everest/scatter_crystal.h,sha256=BHBPNiMKynq1XNxIiWJXYxMwW2bhmmbmC8j7-Dy-SNs,9455
 xcoll/scattering_routines/fluka/build_fluka_input.py,sha256=NDPg191eVTcSrX3nEWdyg0LIWf-gCJI1YWo2NqW9xtU,2270
 xcoll/scattering_routines/fluka/flukaio/.git,sha256=5ZxurpwP1waJk2Zu2AstENGOp4rXU8-C-oU4CypmIVI,73
 xcoll/scattering_routines/fluka/flukaio/.gitignore,sha256=bhrpWiAnKChkcwIaIPX_zKxyG2nCtNzNgQAHnIJa8Fw,12
@@ -104,8 +109,18 @@ xcoll/scattering_routines/fluka/flukaio/tests/fakes/FakeFlukaIOServer.h,sha256=0
 xcoll/scattering_routines/fluka/flukaio/tests/fakes/FakeFortranFlukaIO.h,sha256=WxAjG2SPlWqMG38iHBM8GQwh5yp031jAGeyPi4y9G4c,330
 xcoll/scattering_routines/fluka/flukaio/tests/fakes/FakeNetIO.c,sha256=JcgCWxbT-ZOeG0c9mz1kk3PfXjaazaUk2P_ROtoBKaY,1452
 xcoll/scattering_routines/fluka/flukaio/tests/fakes/FakeNetIO.h,sha256=JTMU9FtCBJkitsnHp7kWKUbPSN4rviRfJgVHK26AuPY,906
-xcoll-0.3.6.dist-info/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
-xcoll-0.3.6.dist-info/METADATA,sha256=-HXnJmvM23LxNmSOOxUuLhk49SuI8DnG7uIWl7DtHPg,2709
-xcoll-0.3.6.dist-info/NOTICE,sha256=6DO_E7WCdRKc42vUoVVBPGttvQi4mRt9fAcxj9u8zy8,74
-xcoll-0.3.6.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
-xcoll-0.3.6.dist-info/RECORD,,
+xcoll/scattering_routines/geometry/__init__.py,sha256=dnq9E5n8TUIIAjooE3zXHFAh_QTlSzpI2pCiuJ17NAs,221
+xcoll/scattering_routines/geometry/collimator_geometry.h,sha256=tfyZl9eVQKqJHGTveErCTpn-FO8EigCAgy3-NkWri2Q,7887
+xcoll/scattering_routines/geometry/crystal_geometry.h,sha256=SVsUHz3AgERmYycVbkAHoNDAUNzAEiBkYHwMMI0SGUc,4887
+xcoll/scattering_routines/geometry/geometry.py,sha256=wvvTWhF3mbZEowZu4oEi7QDxFwQxDrV88N3jkdXAjJ8,1102
+xcoll/scattering_routines/geometry/get_s.h,sha256=1gss-BVzVzD5E4ebPtRUeq4t8eHIBofWA6a84EGTnt4,2947
+xcoll/scattering_routines/geometry/methods.h,sha256=SUo6N7XWsGN9S6gg5RoGNEJ8mmbHuHKzaH9R4R_UDPQ,3747
+xcoll/scattering_routines/geometry/objects.h,sha256=lNFEfqrSDnJuB6sKh59bbg52bGqOHqjt44jl4YlazmI,5499
+xcoll/scattering_routines/geometry/rotation.h,sha256=AvNwcqAAlhjEUtSbTqi53BTvzR256gtfUYcM5jlXUco,701
+xcoll/scattering_routines/geometry/segments.h,sha256=IAz4IWvarIIpX7Waph-Hh2j6o6n_2uCPmHxxQp8Q3yg,7667
+xcoll/scattering_routines/geometry/sort.h,sha256=wNdQLwEn58jgeSG4B9x-jfm7vOjN4hnAlxZ3SG-SkWg,5763
+xcoll-0.4.0.dist-info/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
+xcoll-0.4.0.dist-info/METADATA,sha256=7--PbyNsEcm2q1tjmu4HHEBdv6MzA9DLmyK_VlTqsxA,2709
+xcoll-0.4.0.dist-info/NOTICE,sha256=6DO_E7WCdRKc42vUoVVBPGttvQi4mRt9fAcxj9u8zy8,74
+xcoll-0.4.0.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
+xcoll-0.4.0.dist-info/RECORD,,

xcoll/beam_elements/collimators_src/absorber.h

@@ -1,141 +0,0 @@
-// copyright ############################### #
-// This file is part of the Xcoll Package.   #
-// Copyright (c) CERN, 2024.                 #
-// ######################################### #
-
-#ifndef XCOLL_ABSORBER_H
-#define XCOLL_ABSORBER_H
-
-/*gpufun*/
-int64_t is_within_aperture(LocalParticle* part, double jaw_L, double jaw_R){
-    double const x = LocalParticle_get_x(part);
-    return (int64_t)((x > jaw_R) && (x < jaw_L) );
-}
-
-/*gpufun*/
-void BlackAbsorber_track_local_particle(BlackAbsorberData el, LocalParticle* part0){
-
-    // Collimator active and length
-    int8_t is_active = BlackAbsorberData_get_active(el);
-    is_active       *= BlackAbsorberData_get__tracking(el);
-    double const inactive_front = BlackAbsorberData_get_inactive_front(el);
-    double const inactive_back = BlackAbsorberData_get_inactive_back(el);
-    double const active_length = BlackAbsorberData_get_active_length(el);
-    // Collimator jaws
-    double const jaw_L  = BlackAbsorberData_get_jaw_L(el);
-    double const jaw_R  = BlackAbsorberData_get_jaw_R(el);
-    double const sin_yL = BlackAbsorberData_get_sin_yL(el);
-    double const sin_yR = BlackAbsorberData_get_sin_yR(el);
-    double const jaw_LU = jaw_L - sin_yL*active_length/2.;
-    double const jaw_LD = jaw_L + sin_yL*active_length/2.;
-    double const jaw_RU = jaw_R - sin_yR*active_length/2.;
-    double const jaw_RD = jaw_R + sin_yR*active_length/2.;
-    // TODO: need shortening of active length!
-    // Collimator reference frame
-    double const sin_zL = BlackAbsorberData_get_sin_zL(el);
-    double const cos_zL = BlackAbsorberData_get_cos_zL(el);
-    double const sin_zR = BlackAbsorberData_get_sin_zR(el);
-    double const cos_zR = BlackAbsorberData_get_cos_zR(el);
-    if (fabs(sin_zL-sin_zR) > 1.e-10 || fabs(cos_zL-cos_zR) > 1.e-10 ){
-        kill_all_particles(part0, XC_ERR_NOT_IMPLEMENTED);
-    };
-    double const dx = BlackAbsorberData_get_ref_x(el);
-    double const dy = BlackAbsorberData_get_ref_y(el);
-    // Impact table
-    CollimatorImpactsData record = BlackAbsorberData_getp_internal_record(el, part0);
-    RecordIndex record_index = NULL;
-    if (record){
-        record_index = CollimatorImpactsData_getp__index(record);
-    }
-
-    //start_per_particle_block (part0->part)
-    
-        double s_collimator = LocalParticle_get_s(part);
-
-        // Go to collimator reference system (centered around orbit)
-        XYShift_single_particle(part, dx, dy);
-        SRotation_single_particle(part, sin_zL, cos_zL);
-
-        int64_t is_alive = 1;
-
-        if (!is_active){
-
-            // If collimator not active, replace with drift
-            Drift_single_particle(part, inactive_front+active_length+inactive_back);
-
-        } else {
-
-            int8_t is_tracking = assert_tracking(part, XC_ERR_INVALID_TRACK);
-            if (is_tracking) {
-           
-            // Drift inactive length before jaw
-            Drift_single_particle(part, inactive_front);
-
-            // Store transversal coordinates for potential backtracking later
-            double x_F = LocalParticle_get_x(part);
-//             double y_F = LocalParticle_get_y(part);
-
-            // Check if hit on the collimator jaw at the front
-            is_alive = is_within_aperture(part, jaw_LU, jaw_RU);
-
-            // Continue if the particle didn't hit the collimator
-            if (is_alive){
-
-                // Drift the jaw length
-                Drift_single_particle(part, active_length);
-
-                // Check if hit on the collimator jaw at the back
-                is_alive = is_within_aperture(part, jaw_LD, jaw_RD);
-
-                // TODO: is there a performance difference with nesting the ifs or not?
-                // Continue if the particle didn't hit the collimator
-                if (is_alive){
-
-                    // Drift inactive length after jaw
-                    Drift_single_particle(part, inactive_back);
-
-                } else {
-
-                    // Backtrack to the particle position of impact
-                    // This should only be done if the particle did NOT hit the front jaw
-                    double x_B = LocalParticle_get_x(part);
-//                     double y_B = LocalParticle_get_y(part);
-                    double length;
-
-                    if (x_B > 0){        // Left jaw
-                        length = (jaw_LD - x_B) / (jaw_LD - jaw_LU - x_B + x_F) * active_length;
-                    } else if (x_B < 0){ // Right jaw
-                        length = (jaw_RD - x_B) / (jaw_RD - jaw_RU - x_B + x_F) * active_length;
-                    // TODO: check this
-//                     } else if (y_B > 0){ // Upper jaw
-//                         length = (y_B - jaw_U) / (y_B - y_F) * active_length;
-//                     } else if (y_B < 0){ // Lower jaw
-//                         length = (y_B - jaw_D) / (y_B - y_F) * active_length;
-                    } else {
-                        length = 0;
-                    }
-                    Drift_single_particle(part, -length);
-
-                }
-            }
-            }
-        }
-
-        // Update dead particles
-        if (!is_alive){
-
-            LocalParticle_set_state(part, XC_LOST_ON_ABSORBER);
-
-            // Record impact data
-            CollimatorImpactsData_log(record, record_index, part, XC_ABSORBED);
-        }
-
-        // Move back from collimator reference system
-        SRotation_single_particle(part, -sin_zL, cos_zL);
-        XYShift_single_particle(part, -dx, -dy);
-
-    //end_per_particle_block
-
-}
-
-#endif