search

SBNSoftware / icaruscode

Main/top level repository for ICARUS specific code
11 stars 34 forks source link

SCE Simulation Broken in ICARUS MC #271

Open mrmooney opened 3 years ago

mrmooney commented 3 years ago

SCE simulation is broken in ICARUS MC because the cryostat location was moved and the SCE maps assume a fixed coordinate system that is no longer correct given the change in geometry - see here (Slide 9, right): https://indico.fnal.gov/event/51213/contributions/224937/attachments/148037/190059/ICARUS_TechnicalMeeting_21_10_05.pdf

If someone can tell me the precise coordinates of the TPC active volume for each of the four TPCs, I can make new maps and modify the SpaceCharge service accordingly. I'll also want some cosmics MC sample produced to validate it (10k events at least).

mrmooney commented 3 years ago

I can't find Alessandro Menegolli on here, but he seems like a relevant person to help address this.

SFBayLaser commented 3 years ago

Shouldn't the space charge effects depend on the local coordinate system of the TPC and not the global coordinate system? This would insulate this from any future potential geometry changes... as well, perhaps this would also allow the space charge simulation to be simplified since it could be mirrored between the 4 TPCs?

wesketchum commented 3 years ago

@SFBayLaser 's comment makes sense ... is there a reason this was not done? Are their external-to-TPC considerations?

jzennamo commented 3 years ago

@mrmooney I third the comment from @SFBayLaser, but to address your immediate question here is a slide from Marta Torti: https://sbn-docdb.fnal.gov/cgi-bin/sso/RetrieveFile?docid=21693&filename=ICARUS_geometry_update_26Apr21_v3.pdf&version=3

The TPC coordinates are on slide 16 and 24

PetrilloAtWork commented 3 years ago

Active volume coordinates (setup icaruscode v09_33_00 -q e20:prof && lar -c dump_icarus_geometry.fcl && grep 'active volume box' ICARUS-geometry.txt):

C:0 T:0      active volume box: (-358.49,-181.86,-894.951) -- (-210.29,134.96,0)
C:0 T:1      active volume box: (-358.49,-181.86,0) -- (-210.29,134.96,894.951)
C:0 T:2      active volume box: (-210.14,-181.86,-894.951) -- (-61.94,134.96,0)
C:0 T:3      active volume box: (-210.14,-181.86,0) -- (-61.94,134.96,894.951)
C:1 T:0      active volume box: (61.94,-181.86,-894.951) -- (210.14,134.96,0)
C:1 T:1      active volume box: (61.94,-181.86,0) -- (210.14,134.96,894.951)
C:1 T:2      active volume box: (210.29,-181.86,-894.951) -- (358.49,134.96,0)
C:1 T:3      active volume box: (210.29,-181.86,0) -- (358.49,134.96,894.951)

Collapsed by drift volume:

C:0 S:0      active volume box: (-358.49,-181.86,-894.951) -- (-210.29,134.96,894.951)
C:0 S:1      active volume box: (-210.14,-181.86,-894.951) -- (-61.94,134.96,894.951)
C:1 S:0      active volume box: (61.94,-181.86,-894.951) -- (210.14,134.96,894.951)
C:1 S:1      active volume box: (210.29,-181.86,-894.951) -- (358.49,134.96,894.951)
brucehoward-physics commented 3 years ago

This reminds me... I have some memory of Alessandro considering the geometry components related to the 148.2 cm vs 149.04 cm question noted in doc-23235.

I don't remember - did this result in a slight change in the geometry or a plan for one, since that would then require further changes?

wesketchum commented 2 years ago

@mrmooney any update on this?

brucehoward-physics commented 2 years ago

Following up on this... Looking in a new release at the output of

DetectorPropertiesStandard_test multitpc_detsim_icarus.fcl

I still see Drift distance: 149.04 cm which suggests that there may still be some tinkering needed with the geometry? (since the original line of inquiry to what I mentioned above, as I recall, was why it was saying 149.04 cm instead of 148.2 cm). I put these messages on this thread since it could mean further changes needed for SCE map, but perhaps it should be its own ticket.