Config for propagation in water

[pkalaczynski]

Hi everyone,

I have a rather general question, does the following config make sense for you? It is supposed to be used for propagation of muons through seawater. The minimal distance to be travelled is ~1900m and the maximal, well even up to tens of km (for very inclined events).

{
	"global":
	{
		"seed" : 1,
		"continous_loss_output" : false,
		"only_loss_inside_detector" : false,

		"interpolation":
		{
			"do_interpolation" : true,
			"path_to_tables" : "resources/tables",
			"path_to_tables_readonly" : "resources/tables",
			"do_binary_tables" : false,
			"just_use_readonly_path" : false,
			"max_node_energy" : 1e14,
			"nodes_cross_section" : 100,
            "nodes_continous_randomization" : 200,
            "nodes_propagate" : 1000
		},

		"exact_time" : true,
		"stopping_decay" : true,
		"scattering" : "Highland",

		"brems_multiplier" : 1,
		"photo_multiplier" : 1,
		"ioniz_multiplier" : 1,
		"epair_multiplier" : 1,
        "ioniz": "IonizBetheBlochRossi",
        "epair": "EpairKelnerKokoulinPetrukhin",
		"brems" : "BremsKelnerKokoulinPetrukhin",
		"photo" : "PhotoAbramowiczLevinLevyMaor97",
        "annihilation": "none",
		"compton" : "none",
        "photopair": "none",
        "mupair": "none",
        "weak": "none",
		"lpm" : true,
		"photo_hard_component" : true,
		"photo_shadow" : "ShadowButkevichMikhailov",

		"cuts_infront":
		{
			"e_cut": -1,
			"v_cut": 0.05,
			"cont_rand": true
		},
		"cuts_inside":
		{
			"e_cut": 500,
			"v_cut": -1,
			"cont_rand": false
		},
		"cuts_behind":
		{
			"e_cut": -1,
			"v_cut": -1,
			"cont_rand": false
		},
		"medium": "air",
		"geometry":
		{
			"shape": "sphere",
			"origin": [0, 0, -6372186],
			"outer_radius": 1000000000,
			"inner_radius": 6374134
		}
	},

	"sectors": [
		{
			"hierarchy": 1,
			"medium": "air",
			"density_correction": 0.673,

			"geometry":
			{
				"shape": "sphere",
				"origin": [0, 0, -6372186],
				"outer_radius": 1000000000,
				"inner_radius": 6374134
			},
			"cuts_inside":
			{
				"e_cut": 500,
				"v_cut": -1,
				"cont_rand": false
			},
			"cuts_infront":
			{
				"e_cut": -1,
				"v_cut": 0.05,
				"cont_rand": true
			},
			"cuts_behind":
			{
				"e_cut": -1,
				"v_cut": -1,
				"cont_rand": false
			}
		},
		{
			"hierarchy": 1,
			"medium": "ice",
			"density_correction": 0.832,

			"geometry":
			{
				"shape": "sphere",
				"origin": [0, 0, -6372186],
				"outer_radius": 6374134,
				"inner_radius": 6373934
			},
			"cuts_inside":
			{
				"e_cut": 500,
				"v_cut": -1,
				"cont_rand": false
			},
			"cuts_infront":
			{
				"e_cut": -1,
				"v_cut": 0.05,
				"cont_rand": true
			},
			"cuts_behind":
			{
				"e_cut": -1,
				"v_cut": -1,
				"cont_rand": false
			}
		},
		{
			"hierarchy": 1,
			"medium": "ice",
			"density_correction": 1.005,

			"geometry":
			{
				"shape": "sphere",
				"origin": [0, 0, -6372186],
				"outer_radius": 6373934,
				"inner_radius": 6371324
			},
			"cuts_inside":
			{
				"e_cut": 500,
				"v_cut": -1,
				"cont_rand": false
			},
			"cuts_infront":
			{
				"e_cut": -1,
				"v_cut": 0.05,
				"cont_rand": true
			},
			"cuts_behind":
			{
				"e_cut": -1,
				"v_cut": -1,
				"cont_rand": false
			}
		},
		{
			"hierarchy": 1,
			"medium": "standardrock",
			"density_correction": 1.0,

			"geometry":
			{
				"shape": "sphere",
				"origin": [0, 0, -6372186],
				"outer_radius": 6371324,
				"inner_radius": 0
			},
			"cuts_inside":
			{
				"e_cut": 500,
				"v_cut": -1,
				"cont_rand": false
			},
			"cuts_infront":
			{
				"e_cut": -1,
				"v_cut": 0.05,
				"cont_rand": true
			},
			"cuts_behind":
			{
				"e_cut": -1,
				"v_cut": -1,
				"cont_rand": false
			}
		}
	],

	"detector":
	{
		"shape": "cylinder",
		"origin" : [0, 0, 0],
		"outer_radius": 800,
		"inner_radius": 0,
		"height": 1600
	}
}

[Jean1995]

In principle, you can use the configuration file as a basis for your needs.

You just need to adapt the sectors for your specific use case, since this is currently the IceCube propagation environment.
Furthermore, you should adapt the EnergyCut settings for your use case to that you have a good performance/precision of your simulation.

One last point is that instead of Highland as a scattering model, you might want to use HighlandIntegral, especially if you want to propagate in large steps since this will give you a more accurate modeling of the multiple scattering.

[pkalaczynski]

Thank you for the suggestions! :)