JeffersonLab / remoll

Simulations for the MOLLER Experiment at Jefferson Lab, http://moller.jlab.org
http://jeffersonlab.github.io/remoll/
11 stars 55 forks source link

TM2 and TM3 radial move-new default maps #574

Closed rahmans1 closed 1 year ago

rahmans1 commented 1 year ago

The subcoils in TM2 and TM3 torus magnets in the downstream spectrometer is now shifted outward by radius 3mm to allow more room for belly plate shielding and better clearance with central beam. The new default field map is subcoil_2_3_3mm_full.txt and the realistic asymmetric map for this case is subcoil_2_3_3mm_real_asymmetric.txt. Also adding the realistic asymmetric map for the TM0 torus magnet in upstream region.

Didn't notice any significant deterioration in deconvolution with new default maps. https://moller.jlab.org/DocDB/0009/000981/001/OffsetCoil2and3.pdf

cipriangal commented 1 year ago

@rahmans1 which coils were offset? were the old ones in the wrong place and the b-field produced by TOSCA had a mistake? does this affect the actual coil geometry we have in the gdml?

rahmans1 commented 1 year ago

@rahmans1 which coils were offset? were the old ones in the wrong place and the b-field produced by TOSCA had a mistake? does this affect the actual coil geometry we have in the gdml?

There was no mistake. Subcoils in segment 2 and 3 (TM2 and TM3) in the downstream spectrometer region have been consciously moved outward by 3 mm to allow more room for belly plate shielding if needed. This field map corresponds to this new default configuration. The physical positioning of the coils will be updated in a separate PR. Shouldn't be a noticeable effect as the move is negligible to the coil cross-section. image

cipriangal commented 1 year ago

@rahmans1 thanks for clarifying.

Does this mean that collimator 5 needs to be moved as well?

From the overall results of the deconvolution you ran, I'd like to take a look at the big spreadsheet we usually make. Could you point me to that? The ~1sigma in some of the asymmetries is a bit concerning. We may have to live with it, but I just want to understand it.

rahmans1 commented 1 year ago

@rahmans1 thanks for clarifying.

Does this mean that collimator 5 needs to be moved as well?

From the overall results of the deconvolution you ran, I'd like to take a look at the big spreadsheet we usually make. Could you point me to that? The ~1sigma in some of the asymmetries is a bit concerning. We may have to live with it, but I just want to understand it.

Collimator 5 stays where it is in the global coordinate system. Just the relative position of that element with the coils change since the coils are being moved outward. I am running an a version of deconvolution code that i had from a while ago. So, it would be good if @zdemirog double checked my numbers too. Here is the detailed table for the new default map

zdemirog commented 1 year ago

@rahmans1 thanks for clarifying. Does this mean that collimator 5 needs to be moved as well? From the overall results of the deconvolution you ran, I'd like to take a look at the big spreadsheet we usually make. Could you point me to that? The ~1sigma in some of the asymmetries is a bit concerning. We may have to live with it, but I just want to understand it.

Collimator 5 stays where it is in the global coordinate system. Just the relative position of that element with the coils change since the coils are being moved outward. I am running an a version of deconvolution code that i had from a while ago. So, it would be good if @zdemirog double checked my numbers too. Here is the detailed table for the new default map

I ran the simulation for the shielding configuration 36 by using new default field maps. You can find the results of the deconvolution study (along with ratio plots for the rate weighted radial distribution) here. The link of the detailed tables can be found on page 5.

cipriangal commented 1 year ago

@rahmans1 : I am still looking at the deconvolution output (do you by any chance have a similar table like Zuhal pointed us to?).

However, it seems that if we move TM2 and TM3 we should see some impact on coll6a/b with asymmetric maps (https://moller.jlab.org/DocDB/0009/000974/004/NewCollimatorDesign.pdf). Did you look at this?

cipriangal commented 1 year ago

@rahmans1 : Sorry about the table. Zuhal pointed out that you already sent it to me and I compared it to some recent results.

The first thing to notice is that the e-Al elastic is very different for Ring 2 closed: Ring 2C A f Sakib 74.3459 0.059782 Zuhal 985 0.00245102 Do you have any idea why this is?

The other thing to note is that in terms of f*A there are some big >10% differences for epElastic in ring 4C and epInelasticW1 for Ring 1O, Ring 4C.

How about we discuss this at the detector meeting on Th? I am pretty sure these won't be too big of an issue but we should discuss them so we understand things better.

If you manage to put together coll6a/b plots that would be something we can discuss carefully in the sim meeting as well.

rahmans1 commented 1 year ago

@rahmans1 : Sorry about the table. Zuhal pointed out that you already sent it to me and I compared it to some recent results.

The first thing to notice is that the e-Al elastic is very different for Ring 2 closed: Ring 2C A f Sakib 74.3459 0.059782 Zuhal 985 0.00245102 Do you have any idea why this is?

The other thing to note is that in terms of f*A there are some big >10% differences for epElastic in ring 4C and epInelasticW1 for Ring 1O, Ring 4C.

How about we discuss this at the detector meeting on Th? I am pretty sure these won't be too big of an issue but we should discuss them so we understand things better.

If you manage to put together coll6a/b plots that would be something we can discuss carefully in the sim meeting as well.

The geometry I used to run my numbers are https://github.com/JeffersonLab/remoll/commits/develop head 83e4f37. Since then there have been multiple shielding updates. Some of the difference may be down to that. To truly compare apples to apples, I would probably need to rerun it with the latest develop. I am ok with discussing this on Thursday.