diff --git a/resources/Detectors/densities/ND280UPGRD/ND280UPGRD-v1.dat b/resources/Detectors/densities/ND280UPGRD/ND280UPGRD-v1.dat new file mode 100644 index 00000000..5c0a7eed --- /dev/null +++ b/resources/Detectors/densities/ND280UPGRD/ND280UPGRD-v1.dat @@ -0,0 +1,81 @@ +# Detector model file +# Detector: ND280 Upgrade +# Version: v1 +# Material model file: ND280UPGRD-v1.dat +# Date: 2023-07-11 +# Authors: Ming-Shau Liu +# Notes: +# Includes old ND280 and upgrade elements +# Old: ECAL, TPC, FGD +# New: 2x TPC, 1x Super FGD, 6x TOF, ECAL + +# TDL Dirt +object box 0 0 -142 0 0 0 24 24 300 surr_earth ROCK constant 2.900 + +# Detector Hall +object box 0 0 0 0 0 0 19 19 37 det_hall AIR constant 0.001225 # 0.673atm x 1.205e-3(g/cm3 for 1atm) + +# MAGNETS +# COIL, YOKE, CHAMBER + +object box 0 0 0 0 0 0 5.6 6.1 7.6 UA1_yoke IRON constant 5.54772 +object box 0 0 0 0 0 0 3.5 4.0 7.6 UA1_chamber_1 AIR constant 0.001225 +object box 0 0 0 0 0 0 3.5 4.0 7.4 UA1_coil ALUM constant 1.70330 +object box 0 0 0 0 0 0 3.5 3.6 7.0 UA1_chamber_2 AIR constant 0.001225 + +# NEW MODULES +# 2x TPC, 1x sFGD, 6x TOF + +# 2 new TPC +object box 0 0.74 -2.205 0 0 0 1.795 0.795 2 new_TPC_shell_top TPC_padding constant 0.164 +object box 0 0.74 -2.205 0 0 0 1.700 0.700 2 new_TPC_gas_top GAS constant 1.855 + +object box 0 -0.74 -2.205 0 0 0 1.795 0.795 2 new_TPC_shell_bot TPC_padding constant 0.164 +object box 0 -0.74 -2.205 0 0 0 1.700 0.700 2 new_TPC_gas_bot GAS constant 1.855 + +# 1 new sFGD +object box 0 0 -2.205 0 0 0 2.018 0.64 2.018 sFGD_box sFGD_box constant 0.284 +object box 0 0 -2.205 0 0 0 1.920 0.56 1.920 sFGD_cubes sFGD_cubes constant 1.05 + +# 6 new TOF +object box 0 1.165 -2.205 0 0 0 2.3 0.01 2.3 TOF TOF_scint constant 1.023 +object box 0 -1.165 -2.205 0 0 0 2.3 0.01 2.3 TOF TOF_scint constant 1.023 +object box 1.195 0 -2.205 0 0 0 0.01 2.3 2.3 TOF TOF_scint constant 1.023 +object box -1.195 0 -2.205 0 0 0 0.01 2.3 2.3 TOF TOF_scint constant 1.023 +object box 0 0 -3.26 0 0 0 2.3 2.3 0.01 TOF TOF_scint constant 1.023 +object box 0 0 -1.17 0 0 0 2.3 2.3 0.01 TOF TOF_scint constant 1.023 + +# OLD MODULES +# 4x ECAL, 3x TPC, 2x FGD + +# 3 old TPC +object box 0 0 -0.375 0 0 0 2.5 2.5 1.0 old_TPC_shell old_TPC_padding constant 0.685 +object box 0 0 -0.375 0 0 0 2.47 2.47 0.97 new_TPC_gas GAS constant 1.855 + +object box 0 0 0.99 0 0 0 2.5 2.5 1.0 old_TPC_shell old_TPC_padding constant 0.685 +object box 0 0 0.99 0 0 0 2.47 2.47 0.97 new_TPC_gas GAS constant 1.855 + +object box 0 0 2.355 0 0 0 2.5 2.5 1.0 old_TPC_shell old_TPC_padding constant 0.685 +object box 0 0 2.355 0 0 0 2.47 2.47 0.97 new_TPC_gas GAS constant 1.855 + +# 4x ECAL +object box 0 0 -3.5 0 0 0 2.04 2.04 0.5 ECAL_Upstream SCINT constant 1.05 # Upstream +object box 0 0 3.250 0 0 0 2.04 2.04 0.5 ECAL_Downstream SCINT constant 1.05 # Downstream +object box 1.5 0 1.045 0 0 0 0.5 2.36 3.84 ECAL_Side_1 SCINT constant 1.05 # Side 1 +object box -1.5 0 1.045 0 0 0 0.5 2.36 3.84 ECAL_Side_2 SCINT constant 1.05 # Side 2 +object box 0 1.5 1.045 0 0 0 1.52 0.5 3.84 ECAL_Top SCINT constant 1.05 # Top +object box 0 -1.5 1.045 0 0 0 1.52 0.5 3.84 ECAL_Bottom SCINT constant 1.05 # Bottom + +# 2x FGD +object box 0 0 0.3075 0 0 0 2.300 2.400 0.365 FGDXYBox FGDbox constant 0.06 # XY MODULE +object box 0 0 0.3075 0 0 0 1.843 1.843 0.288 FGDXYScint SCINT constant 1.05 + +object box 0 0 1.6725 0 0 0 2.300 2.400 0.365 FGDWaterBox FGDbox constant 0.06 # SCINT WATER MODULE +object box 0 0 1.6725 0 0 0 1.843 1.843 0.217 FGDWaterScint H2OSCINT constant 1.0134 + + + +detector 0 0 0 + +# Fiducial volume +fiducial box 0 0 0 0 0 0 3.5 3.6 7.0 \ No newline at end of file diff --git a/resources/Detectors/materials/ND280UPGRD/ND280UPGRD-v1.dat b/resources/Detectors/materials/ND280UPGRD/ND280UPGRD-v1.dat new file mode 100644 index 00000000..4304524c --- /dev/null +++ b/resources/Detectors/materials/ND280UPGRD/ND280UPGRD-v1.dat @@ -0,0 +1,88 @@ +# Material model file +# Detector: ND280 Upgrade +# Version: v1 +# Date: 2024-07-14 +# Authors: Ming-Shau Liu +# Notes: +# + +ROCK 2 # SiO2 +1000140280 0.4674349 # Si, 33% +1000080160 0.5325651 # O16, 66% + +IRON 4 # Assume natural iron +1000260540 0.0582000 # Fe54 5.82% +1000260560 0.9166000 # Fe56 91.66% +1000260570 0.0219000 # Fe57 2.19% +1000260580 0.0033000 # Fe58 0.33% + +ALUM 1 # Coil +1000130270 1.0000000 # Al27 100%% + +WATER 2 # H20 +1000010010 0.1150000 # H1 11.5% +1000080160 0.8850000 # 016 88.5% + +AIR 2 # N2 + O2 +1000070140 0.7562326 # N2 78% in volume +1000080160 0.2437674 # O2 22% in volume + +GAS 4 # Ar-CF4-iC4H10 95 3 2 +1000010010 0.0034483 # H1 +1000060120 0.0206426 # C12 +1000090190 0.0259091 # F19 +1000180400 0.95 # Ar40 + + +SCINT 11 +1000060120 0.8627360 # C12 +1000080160 0.0368494 # O16 88.5% +1000010010 0.0741594 # H1 11.5% +1000220460 0.0013301 # Ti46 +1000220470 0.0011994 # Ti47 +1000220480 0.0118848 # Ti48 +1000220490 0.0008722 # Ti49 +1000220500 0.0008351 # Ti50 +1000140280 0.0093432 # Si28 +1000140290 0.0004763 # Si29 +1000140300 0.0003141 # Si30 + +H2OSCINT 3 +1000060120 0.1532317 # C12 15.3% +1000080160 0.7418300 # O16 74.2% +1000010010 0.1049383 # H1 10.5% + +FGDbox # Epoxy fiberglass based C21H25ClO5 +1000060120 0.6428572 # C12 +1000010010 0.0637755 # H1 +1000080160 0.2040816 # O16 +1000170350 0.0892857 # Cl35 + +old_TPC_padding 4 # Mostly ROHACELL PMI (C4H7NO), exclude G10 epoxyglass copper clad +1000060120 0.5647059 # C12 +1000010010 0.0823529 # H1 +1000080160 0.1882353 # O16 +1000070140 0.1647059 # N14 + +TOF_scint 2 # EJ200 Scintillator, Polyvinyltoluene based CH2CH(C6H4CH3) +1000060120 0.9152542 # C12 +1000010010 0.0847458 # H1 + +sFGD_box 3 +1000060120 0.6703297 # C12 +1000010010 0.0366300 # H1 +1000080160 0.2930403 # O16 + +sFGD_cubes 4 +1000060120 0.9235685 # C12 +1000010010 0.0763900 # H1 +1000080160 0.0000151 # O16 +1000070140 0.0000264 # N14 + +TPC_padding # aramid material, PPTA based +1000060120 0.69263638 # C12 +1000010010 0.04122835 # H1 +1000080160 0.13193075 # O16 +1000070140 0.11543933 # N14 +1000290630 0.01294798 # Cu63 +1000290650 0.00581721 # Cu65 \ No newline at end of file diff --git a/resources/Detectors/visuals/Det_Visual_Examples.nb b/resources/Detectors/visuals/Det_Visual_Examples.nb new file mode 100644 index 00000000..6710d405 --- /dev/null +++ b/resources/Detectors/visuals/Det_Visual_Examples.nb @@ -0,0 +1,856 @@ +(* Content-type: application/vnd.wolfram.mathematica *) + +(*** Wolfram Notebook File ***) +(* http://www.wolfram.com/nb *) + +(* CreatedBy='Mathematica 13.3' *) + +(*CacheID: 234*) +(* Internal cache information: +NotebookFileLineBreakTest +NotebookFileLineBreakTest +NotebookDataPosition[ 158, 7] +NotebookDataLength[ 41346, 848] +NotebookOptionsPosition[ 38267, 784] +NotebookOutlinePosition[ 38662, 800] +CellTagsIndexPosition[ 38619, 797] +WindowFrame->Normal*) + +(* Beginning of Notebook Content *) +Notebook[{ + +Cell[CellGroupData[{ +Cell["Detector visualisation", "Section", + CellChangeTimes->{{3.9308471076304407`*^9, + 3.930847130291616*^9}},ExpressionUUID->"ff56253e-e079-4247-8f6b-\ +78a3d0bd9d71"], + +Cell["\<\ +Visuals[ \[OpenCurlyDoubleQuote]Experiment\[CloseCurlyDoubleQuote], { \ +Excluded Modules } ] +displays the experiment as a 3D model without the specified { Excluded \ +Modules }, e.g., +\ +\>", "Text", + CellChangeTimes->{{3.930847526527903*^9, 3.9308476997779016`*^9}, + 3.930927278187214*^9},ExpressionUUID->"2b8b62a5-ee55-4e21-b7e9-\ +9cf4f37c2621"], + +Cell[CellGroupData[{ + +Cell["Set working directory to SIREN and call function ", "Subsubsection", + CellChangeTimes->{{3.9308471343658047`*^9, 3.9308471514553337`*^9}, { + 3.9308475183650503`*^9, + 3.930847522740533*^9}},ExpressionUUID->"1828e000-2b55-4658-a39d-\ +d3ec0e65319c"], + +Cell["Save this notebook /path/to/SIREN directory ", "Text", + CellChangeTimes->{{3.930918106474534*^9, 3.930918160577179*^9}, { + 3.930918306420301*^9, + 3.9309183092048483`*^9}},ExpressionUUID->"b85c7e71-3fe5-4770-92d2-\ +3b8310f025c5"], + +Cell[BoxData[ + RowBox[{ + RowBox[{"AppendTo", "[", + RowBox[{"$Path", ",", + RowBox[{"NotebookDirectory", "[", "]"}]}], "]"}], ";"}]], "Input", + CellChangeTimes->{{3.930845768512031*^9, 3.930845780337503*^9}, { + 3.93084581693408*^9, 3.930845853905582*^9}, {3.93091362763332*^9, + 3.9309136464460287`*^9}}, + CellLabel->"In[2]:=",ExpressionUUID->"01061750-99ce-4ee2-ba52-6ccdbfbf6f3b"], + +Cell[BoxData[ + RowBox[{"<<", "DetectorVisuals`"}]], "Input", + CellChangeTimes->{{3.930845856938705*^9, 3.930845892210829*^9}, { + 3.9308460624868727`*^9, 3.930846074171977*^9}, {3.9309136516632137`*^9, + 3.930913666071559*^9}}, + CellLabel->"In[35]:=",ExpressionUUID->"a6dc3d65-b39c-4f5d-8e40-9ac1c63e10c6"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[TextData[StyleBox["Examples", "Subsubsection"]], "Subsection", + CellChangeTimes->{{3.930847164943055*^9, 3.930847192667643*^9}, { + 3.930847714726037*^9, 3.930847716563344*^9}, + 3.930914921181027*^9},ExpressionUUID->"5b67d770-d2b6-44d9-871c-\ +110d0d66dc83"], + +Cell[CellGroupData[{ + +Cell["Target detectors", "Subsubsection", + CellChangeTimes->{{3.930848381829659*^9, + 3.9308483931419697`*^9}},ExpressionUUID->"b4bb631f-eef9-40d0-aa86-\ +e66af2411a5f"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", + RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], + "]"}]], "Input", + CellChangeTimes->{{3.9308461111501207`*^9, 3.930846132819652*^9}, { + 3.9308462468135443`*^9, 3.930846247147883*^9}, {3.930846995520549*^9, + 3.9308470458198357`*^9}}, + CellLabel->"In[4]:=",ExpressionUUID->"758f5770-9701-4c2e-9c28-f1655b0a6965"], + +Cell[BoxData[ + Graphics3DBox[{ + {RGBColor[0.35093217511999275`, 0.3525201340616513, 0.8276327252716125], + Opacity[0.1], CuboidBox[{-2.8, -3.05, -3.8}, {2.8, 3.05, 3.8}]}, + {RGBColor[0.8477780541436477, 0.3160906000834609, 0.219018325039134], + Opacity[0.1], CuboidBox[{-1.75, -2., -3.8}, {1.75, 2., 3.8}]}, + {RGBColor[0.22330177948027363`, 0.31121718977849944`, 0.08734595030416648], + Opacity[0.1], CuboidBox[{-1.75, -2., -3.7}, {1.75, 2., 3.7}]}, + {RGBColor[0.8477780541436477, 0.3160906000834609, 0.219018325039134], + Opacity[0.1], CuboidBox[{-1.75, -1.8, -3.5}, {1.75, 1.8, 3.5}]}, + {RGBColor[0.5014698837210463, 0.8279721286735073, 0.25805733064978975`], + Opacity[0.1], + CuboidBox[{-0.8975, 0.34249999999999997`, -3.205}, {0.8975, + 1.1375, -1.205}]}, + {RGBColor[0.6907727880028507, 0.11798104164182699`, 0.5416072831768857], + Opacity[0.1], + CuboidBox[{-0.85, 0.39, -3.205}, {0.85, 1.0899999999999999`, -1.205}]}, + {RGBColor[0.5014698837210463, 0.8279721286735073, 0.25805733064978975`], + Opacity[0.1], + CuboidBox[{-0.8975, -1.1375, -3.205}, { + 0.8975, -0.34249999999999997`, -1.205}]}, + {RGBColor[0.6907727880028507, 0.11798104164182699`, 0.5416072831768857], + Opacity[0.1], + CuboidBox[{-0.85, -1.0899999999999999`, -3.205}, { + 0.85, -0.39, -1.205}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], + CuboidBox[{-1.009, -0.32, -3.214}, {1.009, + 0.32, -1.1960000000000002`}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], CuboidBox[{-0.96, -0.28, -3.165}, {0.96, 0.28, -1.245}]}, + {RGBColor[0.25586105262190606`, 0.2633788353773121, 0.2901412060568902], + Opacity[0.1], + CuboidBox[{-1.15, 1.1600000000000001`, -3.355}, {1.15, + 1.17, -1.0550000000000002`}]}, + {RGBColor[0.25586105262190606`, 0.2633788353773121, 0.2901412060568902], + Opacity[0.1], + CuboidBox[{-1.15, -1.17, -3.355}, { + 1.15, -1.1600000000000001`, -1.0550000000000002`}]}, + {RGBColor[0.25586105262190606`, 0.2633788353773121, 0.2901412060568902], + Opacity[0.1], + CuboidBox[{1.1900000000000002`, -1.15, -3.355}, {1.2, + 1.15, -1.0550000000000002`}]}, + {RGBColor[0.25586105262190606`, 0.2633788353773121, 0.2901412060568902], + Opacity[0.1], + CuboidBox[{-1.2, -1.15, -3.355}, {-1.1900000000000002`, + 1.15, -1.0550000000000002`}]}, + {RGBColor[0.25586105262190606`, 0.2633788353773121, 0.2901412060568902], + Opacity[0.1], + CuboidBox[{-1.15, -1.15, -3.2649999999999997`}, {1.15, 1.15, -3.255}]}, + {RGBColor[0.25586105262190606`, 0.2633788353773121, 0.2901412060568902], + Opacity[0.1], + CuboidBox[{-1.15, -1.15, -1.1749999999999998`}, {1.15, 1.15, -1.165}]}, + {RGBColor[0.7458525322351903, 0.7847146789223423, 0.9181634201906055], + Opacity[0.1], CuboidBox[{-1.25, -1.25, -0.875}, {1.25, 1.25, 0.125}]}, + {RGBColor[0.6907727880028507, 0.11798104164182699`, 0.5416072831768857], + Opacity[0.1], + CuboidBox[{-1.235, -1.235, -0.86}, {1.235, 1.235, 0.10999999999999999`}]}, + {RGBColor[0.7458525322351903, 0.7847146789223423, 0.9181634201906055], + Opacity[0.1], CuboidBox[{-1.25, -1.25, 0.49}, {1.25, 1.25, 1.49}]}, + {RGBColor[0.6907727880028507, 0.11798104164182699`, 0.5416072831768857], + Opacity[0.1], + CuboidBox[{-1.235, -1.235, 0.505}, {1.235, 1.235, 1.475}]}, + {RGBColor[0.7458525322351903, 0.7847146789223423, 0.9181634201906055], + Opacity[0.1], CuboidBox[{-1.25, -1.25, 1.855}, {1.25, 1.25, 2.855}]}, + {RGBColor[0.6907727880028507, 0.11798104164182699`, 0.5416072831768857], + Opacity[0.1], CuboidBox[{-1.235, -1.235, 1.87}, {1.235, 1.235, 2.84}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], CuboidBox[{-1.02, -1.02, -3.75}, {1.02, 1.02, -3.25}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], CuboidBox[{-1.02, -1.02, 3.}, {1.02, 1.02, 3.5}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], CuboidBox[{1.25, -1.18, -0.875}, {1.75, 1.18, 2.965}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], CuboidBox[{-1.75, -1.18, -0.875}, {-1.25, 1.18, 2.965}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], CuboidBox[{-0.76, 1.25, -0.875}, {0.76, 1.75, 2.965}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], CuboidBox[{-0.76, -1.75, -0.875}, {0.76, -1.25, 2.965}]}, + {RGBColor[0.7212050307872955, 0.3081619831157636, 0.9100551807755652], + Opacity[0.1], CuboidBox[{-1.15, -1.2, 0.125}, {1.15, 1.2, 0.49}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], + CuboidBox[{-0.9215, -0.9215, 0.1635}, {0.9215, 0.9215, 0.4515}]}, + {RGBColor[0.7212050307872955, 0.3081619831157636, 0.9100551807755652], + Opacity[0.1], + CuboidBox[{-1.15, -1.2, 1.4900000000000002`}, {1.15, 1.2, 1.855}]}, + {RGBColor[0.20715118379631448`, 0.30858200457649687`, 0.3717365440678644], + Opacity[0.1], + CuboidBox[{-0.9215, -0.9215, 1.564}, {0.9215, 0.9215, + 1.7810000000000001`}]}}, + ImageSize->{326.05528778058215`, 276.80632953762915`}, + ImageSizeRaw->Automatic, + ViewPoint->{-3.267340612285325, -0.6867651866843811, 0.5503080061832035}, + ViewVertical->{-0.7653559255008713, 0.6390814102561185, + 0.07619224609995254}]], "Output", + CellChangeTimes->{{3.9308461275306063`*^9, 3.930846133194714*^9}, + 3.930846247887536*^9, 3.930846335323381*^9, 3.930846381761827*^9, + 3.930846502477841*^9, {3.930846906082329*^9, 3.9308469304566803`*^9}, { + 3.930846961387842*^9, 3.930847046949695*^9}, 3.930848431125566*^9, + 3.930913669409268*^9, 3.93091388668988*^9, 3.930995105650139*^9, { + 3.930995673845676*^9, 3.930995696466717*^9}, 3.93099581312564*^9, + 3.930995881446875*^9, 3.930995996055388*^9, 3.930996573196425*^9, + 3.9309969008185797`*^9, 3.930997477608238*^9, 3.93100477589196*^9, + 3.931005380832015*^9, {3.931005967116878*^9, 3.931005978085596*^9}, + 3.931008220329236*^9, 3.931008284734171*^9}, + CellLabel->"Out[4]=",ExpressionUUID->"8a3d7fe4-05dd-4061-b3db-44a8b527c5fa"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", + RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], + "]"}]], "Input", + CellChangeTimes->{{3.9308464581373568`*^9, 3.9308464652613783`*^9}, { + 3.930847248876378*^9, 3.9308472858524933`*^9}}, + CellLabel->"In[39]:=",ExpressionUUID->"85f2f25e-059b-47e1-b7b7-1cd9a7fb936b"], + +Cell[BoxData[ + Graphics3DBox[{ + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], + CuboidBox[{5.0001999999999995`, -2.9465, -2.}, {11.1978, + 2.9465, -1.3900000000000001`}]}, + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], + CuboidBox[{5., -0.9144, -1.3900000000000001`}, {5.660399999999999, 0.9144, + 0.464}]}, + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], + CuboidBox[{5.6604, -0.9398, -1.3900000000000001`}, {6.574800000000001, + 0.9398, 0.464}]}, + {RGBColor[0.32506606707200136`, 0.23990300858724734`, 0.486881531810764], + Opacity[0.1], + CuboidBox[{6.5748, -0.9398, -1.3900000000000001`}, {6.7780000000000005`, + 0.9398, 0.464}]}, + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], + CuboidBox[{6.778, -1.2446, -1.3900000000000001`}, {7.2352, 1.2446, + 0.464}]}, + {RGBColor[0.521412842051586, 0.4962344456641059, 0.6541153328529785], + Opacity[0.1], + CuboidBox[{7.2352, -1.8288, -1.3900000000000001`}, {7.3114, 1.8288, + 0.464}]}, + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], + CuboidBox[{8.311399999999999, -1.8288, -1.3900000000000001`}, {8.7686, + 1.8288, 0.464}]}, + {RGBColor[0.32506606707200136`, 0.23990300858724734`, 0.486881531810764], + Opacity[0.1], + CuboidBox[NCache[{18., Rational[-5, 2], -2.}, {18., -2.5, -2.}], + NCache[{19., Rational[5, 2], 1.}, {19., 2.5, 1.}]]}, + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], + CuboidBox[NCache[{19., Rational[-5, 2], -2.}, {19., -2.5, -2.}], + NCache[{21., Rational[5, 2], 1.}, {21., 2.5, 1.}]]}, + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], CuboidBox[{21., 4., -2.}, {25., 5., 1.}]}, + {RGBColor[0.41718218051431455`, 0.6517933847257689, 0.5388337863559993], + Opacity[0.1], CuboidBox[{21., -5., -2.}, {25., -4., 1.}]}, + {RGBColor[0.32506606707200136`, 0.23990300858724734`, 0.486881531810764], + Opacity[0.1], CylinderBox[{{0, 0, -2}, {0, 0, 2}}, 5]}, + {RGBColor[0.32506606707200136`, 0.23990300858724734`, 0.486881531810764], + Opacity[0.1], + CylinderBox[{{0, 0, -0.6299999999999999}, {0, 0, 1.}}, 0.83]}, + {RGBColor[0.0848019125038495, 0.8968379325252431, 0.05613755639699236], + Opacity[0.1], + CylinderBox[{{0, 0, -0.44999999999999996`}, {0, 0, 0.79}}, 0.55]}, + {RGBColor[0.7911505744190059, 0.18559734430518615`, 0.2453924328324959], + Opacity[0.1], CylinderBox[{{0, 0, -0.391}, {0, 0, -0.241}}, 0.3]}, + {RGBColor[0.7911505744190059, 0.18559734430518615`, 0.2453924328324959], + Opacity[0.1], + CylinderBox[{{0, 0, -0.241}, {0, 0, -0.09100000000000001}}, 0.3]}, + {RGBColor[0.7911505744190059, 0.18559734430518615`, 0.2453924328324959], + Opacity[0.1], + CylinderBox[{{0, 0, 0.07999999999999999}, {0, 0, 0.5}}, 0.3]}, + {RGBColor[0.3351807133351079, 0.2763153120327406, 0.7227930518683081], + Opacity[0.1], + CylinderBox[{{0, 0, 0.09200000000000001}, {0, 0, 0.183}}, 0.05]}, + {RGBColor[0.3351807133351079, 0.2763153120327406, 0.7227930518683081], + Opacity[0.1], CylinderBox[{{0, 0, -0.39}, {0, 0, -0.092}}, 0.05]}, + {RGBColor[0.32506606707200136`, 0.23990300858724734`, 0.486881531810764], + Opacity[0.1], CylinderBox[{{23, 0, -1.96}, {23, 0, 0.66}}, 1.38]}, + {RGBColor[0.009276515470650892, 0.6652402538800535, 0.2557510553491851], + Opacity[0.1], + CylinderBox[{{23, 0, -1.9100000000000001`}, {23, 0, 0.61}}, 1.35]}, + {RGBColor[0.32506606707200136`, 0.23990300858724734`, 0.486881531810764], + Opacity[0.1], + CylinderBox[{{23, 0, -1.85}, {23, 0, 0.5499999999999999}}, 1.25]}, + {RGBColor[0.6493869773926908, 0.16552659229722266`, 0.1134446897889001], + Opacity[0.1], + CylinderBox[{{23, 0, -1.7999999999999998`}, {23, 0, 0.4999999999999999}}, + 1.2]}, + {RGBColor[0.3086780992491793, 0.10261299404488455`, 0.5061987340813581], + Opacity[0.1], + CylinderBox[{{23, 0, -1.4}, {23, 0, 0.09999999999999998}}, 1.06]}, + {RGBColor[0.6493869773926908, 0.16552659229722266`, 0.1134446897889001], + Opacity[0.1], + CylinderBox[{{23, 0, -1.266}, {23, 0, -0.03400000000000003}}, + 0.96]}}]], "Output", + CellChangeTimes->{ + 3.930847286401828*^9, 3.9309138898128967`*^9, {3.930995702475547*^9, + 3.9309957057441807`*^9}, 3.930996904682377*^9, 3.9309974796076202`*^9, + 3.931005982557057*^9, 3.93100621941457*^9, 3.9310087892293873`*^9, + 3.9310089841360407`*^9, 3.93100907290653*^9, 3.931009196389038*^9, + 3.9310093752491703`*^9}, + CellLabel->"Out[39]=",ExpressionUUID->"5a905629-0c22-4620-b59c-02655db9bf79"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", + RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], + "]"}]], "Input", + CellChangeTimes->{{3.9309951152387*^9, 3.930995153174532*^9}, + 3.930997063585266*^9, 3.931005817453622*^9, {3.9310092389209433`*^9, + 3.9310092757902718`*^9}, {3.9310093189732447`*^9, + 3.931009348185121*^9}},ExpressionUUID->"45186205-4e3c-43da-8cf6-\ +942bd6843b62"], + +Cell[BoxData[ + Graphics3DBox[{ + {RGBColor[0.2298911743818579, 0.31209133471148953`, 0.7146612930358287], + Opacity[0.1], CuboidBox[{-1, -1, -1.0375}, {1, 1, -0.9625}]}, + {RGBColor[0.16452635096365742`, 0.23517942424424865`, 0.9345244778129824], + Opacity[0.1], + PolyhedronBox[{{{0., 1.23553, -0.0027999999999996916`}, {1.07, + 0.61776, -0.0027999999999996916`}, { + 1.07, -0.61776, -0.0027999999999996916`}, { + 0., -1.23553, -0.0027999999999996916`}, {-1.07, -0.61776, \ +-0.0027999999999996916`}, {-1.07, 0.61776, -0.0027999999999996916`}}, {{0., + 1.23553, 4.1372}, {1.07, 0.61776, 4.1372}, {1.07, -0.61776, 4.1372}, { + 0., -1.23553, 4.1372}, {-1.07, -0.61776, 4.1372}, {-1.07, 0.61776, + 4.1372}}}]}, + {RGBColor[0.3275193083067711, 0.5033666810361557, 0.1674300296182314], + Opacity[0.1], + PolyhedronBox[{{{0., 1.06232, 0.}, {0.92, 0.53116, 0.}, {0.92, -0.53116, + 0.}, {0., -1.06232, 0.}, {-0.92, -0.53116, 0.}, {-0.92, 0.53116, 0.}}, {{ + 0., 1.06232, 4.1344}, {0.92, 0.53116, 4.1344}, {0.92, -0.53116, + 4.1344}, {0., -1.06232, 4.1344}, {-0.92, -0.53116, 4.1344}, {-0.92, + 0.53116, 4.1344}}}]}, + {RGBColor[0.2298911743818579, 0.31209133471148953`, 0.7146612930358287], + Opacity[0.1], + PolyhedronBox[{{{-0.2825, 0.8992, 0.136}, {0.6375, -0.6942, 0.136}, { + 0., -1.06232, 0.136}, {-0.92, -0.53116, 0.136}, {-0.92, 0.53116, + 0.136}}, {{-0.2825, 0.8992, 0.16167}, {0.6375, -0.6942, 0.16167}, { + 0., -1.06232, 0.16167}, {-0.92, -0.53116, 0.16167}, {-0.92, 0.53116, + 0.16167}}}]}, + {RGBColor[0.8775095682131331, 0.5860897175157826, 0.24202862693045568`], + Opacity[0.1], + PolyhedronBox[{{{0.6375, -0.6942, 0.136}, {-0.2825, 0.8992, 0.136}, {0., + 1.06232, 0.136}, {0.92, 0.53116, 0.136}, {0.92, -0.53116, 0.136}}, {{ + 0.6375, -0.6942, 0.16178}, {-0.2825, 0.8992, 0.16178}, {0., 1.06232, + 0.16178}, {0.92, 0.53116, 0.16178}, {0.92, -0.53116, 0.16178}}}]}, + {RGBColor[0.2298911743818579, 0.31209133471148953`, 0.7146612930358287], + Opacity[0.1], + PolyhedronBox[{{{-0.6375, -0.6942, 0.313}, {0.2825, 0.8992, 0.313}, {0.92, + 0.53116, 0.313}, {0.92, -0.53116, 0.313}, {0., -1.06232, + 0.313}}, {{-0.6375, -0.6942, 0.33863}, {0.2825, 0.8992, 0.33863}, {0.92, + 0.53116, 0.33863}, {0.92, -0.53116, 0.33863}, {0., -1.06232, + 0.33863}}}]}, + {RGBColor[0.8775095682131331, 0.5860897175157826, 0.24202862693045568`], + Opacity[0.1], + PolyhedronBox[{{{0.2825, 0.8992, 0.313}, {-0.6375, -0.6942, + 0.313}, {-0.92, -0.53116, 0.313}, {-0.92, 0.53116, 0.313}, {0., 1.06232, + 0.313}}, {{0.2825, 0.8992, 0.33881}, {-0.6375, -0.6942, + 0.33881}, {-0.92, -0.53116, 0.33881}, {-0.92, 0.53116, 0.33881}, {0., + 1.06232, 0.33881}}}]}, + {RGBColor[0.2654164649210198, 0.39657426393543016`, 0.2847295993909016], + Opacity[0.1], + PolyhedronBox[{{{0., 1.06232, 0.534}, {0.92, 0.53116, + 0.534}, {-0.92, -0.53116, 0.534}, {-0.92, 0.53116, 0.534}}, {{0., + 1.06232, 0.6102000000000001}, {0.92, 0.53116, + 0.6102000000000001}, {-0.92, -0.53116, 0.6102000000000001}, {-0.92, + 0.53116, 0.6102000000000001}}}]}, + {RGBColor[0.2298911743818579, 0.31209133471148953`, 0.7146612930358287], + Opacity[0.1], + PolyhedronBox[{{{0., 0., 0.534}, {0.92, 0.53116, 0.534}, {0.92, -0.53116, + 0.534}, {0., -1.06232, 0.534}}, {{0., 0., 0.5597300000000001}, {0.92, + 0.53116, 0.5597300000000001}, {0.92, -0.53116, 0.5597300000000001}, { + 0., -1.06232, 0.5597300000000001}}}]}, + {RGBColor[0.8775095682131331, 0.5860897175157826, 0.24202862693045568`], + Opacity[0.1], + PolyhedronBox[{{{0., 0., 0.534}, {0., -1.06232, 0.534}, {-0.92, -0.53116, + 0.534}}, {{0., 0., 0.5596300000000001}, {0., -1.06232, + 0.5596300000000001}, {-0.92, -0.53116, 0.5596300000000001}}}]}, + {RGBColor[0.8332620257988352, 0.5659513355920436, 0.30336281249235686`], + Opacity[0.1], + PolyhedronBox[{{{0., 1.06232, 0.895}, {0.92, 0.53116, 0.895}, { + 0.92, -0.53116, 0.895}, {0., -1.06232, 0.895}, {-0.92, -0.53116, + 0.895}, {-0.92, 0.53116, 0.895}}, {{0., 1.06232, 1.0756000000000001`}, { + 0.92, 0.53116, 1.0756000000000001`}, {0.92, -0.53116, + 1.0756000000000001`}, {0., -1.06232, + 1.0756000000000001`}, {-0.92, -0.53116, 1.0756000000000001`}, {-0.92, + 0.53116, 1.0756000000000001`}}}]}, + {RGBColor[0.8775095682131331, 0.5860897175157826, 0.24202862693045568`], + Opacity[0.1], + PolyhedronBox[{{{0., 1.06232, 1.256}, {0.92, 0.53116, 1.256}, { + 0.92, -0.53116, 1.256}, {0., -1.06232, 1.256}, {-0.92, -0.53116, + 1.256}, {-0.92, 0.53116, 1.256}}, {{0., 1.06232, 1.26395}, {0.92, + 0.53116, 1.26395}, {0.92, -0.53116, 1.26395}, {0., -1.06232, + 1.26395}, {-0.92, -0.53116, 1.26395}, {-0.92, 0.53116, 1.26395}}}]}, + {RGBColor[0.2298911743818579, 0.31209133471148953`, 0.7146612930358287], + Opacity[0.1], + PolyhedronBox[{{{-0.2825, 0.8992, 1.389}, {0.6375, -0.6942, 1.389}, { + 0., -1.06232, 1.389}, {-0.92, -0.53116, 1.389}, {-0.92, 0.53116, + 1.389}}, {{-0.2825, 0.8992, 1.40189}, {0.6375, -0.6942, 1.40189}, { + 0., -1.06232, 1.40189}, {-0.92, -0.53116, 1.40189}, {-0.92, 0.53116, + 1.40189}}}]}, + {RGBColor[0.8775095682131331, 0.5860897175157826, 0.24202862693045568`], + Opacity[0.1], + PolyhedronBox[{{{0.6375, -0.6942, 1.389}, {-0.2825, 0.8992, 1.389}, {0., + 1.06232, 1.389}, {0.92, 0.53116, 1.389}, {0.92, -0.53116, 1.389}}, {{ + 0.6375, -0.6942, 1.40217}, {-0.2825, 0.8992, 1.40217}, {0., 1.06232, + 1.40217}, {0.92, 0.53116, 1.40217}, {0.92, -0.53116, 1.40217}}}]}}, + ImageSize->{405.47170519299516`, 208.51104838001987`}, + ImageSizeRaw->Automatic, + ViewPoint->{-0.7702959971782253, -2.5050506072020586`, 2.140412467747233}, + ViewVertical->{0.7236335114618405, -0.6677828414185077, + 0.17441507330626652`}]], "Output", + CellChangeTimes->{ + 3.930995153669259*^9, 3.930995711063602*^9, 3.930995816691572*^9, + 3.930995885371175*^9, 3.930996011606841*^9, {3.930996567202976*^9, + 3.930996579467017*^9}, 3.930996743033431*^9, 3.930996909507143*^9, { + 3.9309970583744297`*^9, 3.930997064011712*^9}, 3.9309974820353603`*^9, + 3.93099761107055*^9, 3.931004779411798*^9, 3.9310048241809387`*^9, + 3.9310053848916073`*^9, 3.931005819054318*^9, 3.931005870348734*^9, { + 3.931005905747102*^9, 3.93100590769261*^9}, 3.931005988255396*^9, + 3.931006222390541*^9, {3.9310082243420467`*^9, 3.931008227966771*^9}, + 3.931008289070985*^9, {3.931008607676939*^9, 3.931008627241765*^9}, + 3.931008698778192*^9, 3.93100874878646*^9, 3.931008989014098*^9, + 3.93100908020753*^9, 3.9310091171045313`*^9, 3.931009199983571*^9, { + 3.931009253670957*^9, 3.931009276283464*^9}}, + CellLabel->"Out[38]=",ExpressionUUID->"7aa043b8-c09a-4a90-8f3b-914c780b39b6"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", "\"\\"", "}"}]}], "]"}]], "Input", + CellChangeTimes->{{3.9309213638291893`*^9, 3.930921408318699*^9}}, + CellLabel->"In[40]:=",ExpressionUUID->"623517c9-0b0a-43b7-941c-607589cb47ec"], + +Cell[BoxData[ + Graphics3DBox[{ + {RGBColor[0.9999669336121517, 0.957806585512984, 0.28830420784733946`], + Opacity[0.1], SphereBox[{0, 0, 0}, 9.1]}, + {RGBColor[0.3248443069364404, 0.6150261091426228, 0.006485477421444896], + Opacity[0.1], SphereBox[{0, 0, 0}, 6.1]}}]], "Output", + CellChangeTimes->{ + 3.93092140928104*^9, 3.930995708508224*^9, 3.930995999878187*^9, + 3.930996576802916*^9, 3.930996907051785*^9, 3.930997612674156*^9, + 3.931005382662628*^9, {3.931008610552994*^9, 3.931008629581517*^9}, + 3.9310086961250677`*^9, 3.931008747296633*^9, 3.931008986563098*^9, + 3.9310090761618843`*^9, 3.931009114241016*^9, 3.93100919847633*^9, + 3.931009377767623*^9}, + CellLabel->"Out[40]=",ExpressionUUID->"928e908c-d246-44d0-8e3f-99106d65c374"] +}, Open ]] +}, Open ]], + +Cell[CellGroupData[{ + +Cell["Large fiducial detectors", "Subsubsection", + CellChangeTimes->{{3.930848363932369*^9, + 3.930848370955418*^9}},ExpressionUUID->"7e1653cb-c490-4521-8bec-\ +792810d94dea"], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"GraphicsRow", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", "}"}]}], "]"}], ",", "\[IndentingNewLine]", + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", + RowBox[{ + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\""}], + "}"}]}], "]"}]}], "}"}], "]"}]], "Input", + CellChangeTimes->{{3.9308473708407297`*^9, 3.9308473846505938`*^9}, { + 3.9308474813020153`*^9, 3.930847490639382*^9}, {3.930847780142088*^9, + 3.93084785043423*^9}, {3.9308478814387417`*^9, 3.9308479120262003`*^9}}, + CellLabel->"In[7]:=",ExpressionUUID->"d7a1ae22-5013-434c-9f40-0ccdc86a3ce7"], + +Cell[BoxData[ + GraphicsBox[{{}, {InsetBox[ + Graphics3DBox[{ + {RGBColor[ + 0.33604086192392857`, 0.9653758014770579, 0.8082011998173695], + Opacity[0.1], SphereBox[{0, 0, 0}, 6478000]}, + {RGBColor[ + 0.7454100290394681, 0.7665909595605325, 0.16392086034407471`], + Opacity[0.1], SphereBox[{0, 0, 0}, 6371324]}, + {RGBColor[ + 0.7454100290394681, 0.7665909595605325, 0.16392086034407471`], + Opacity[0.1], SphereBox[{0, 0, 0}, 6356000]}, + {RGBColor[ + 0.17011058031535153`, 0.08842194355669797, 0.22208792588679205`], + Opacity[0.1], SphereBox[{0, 0, 0}, 6346600]}, + {RGBColor[ + 0.17011058031535153`, 0.08842194355669797, 0.22208792588679205`], + Opacity[0.1], SphereBox[{0, 0, 0}, 6151000]}, + {RGBColor[ + 0.17011058031535153`, 0.08842194355669797, 0.22208792588679205`], + Opacity[0.1], SphereBox[{0, 0, 0}, 5971000]}, + {RGBColor[ + 0.17011058031535153`, 0.08842194355669797, 0.22208792588679205`], + Opacity[0.1], SphereBox[{0, 0, 0}, 5771000]}, + {RGBColor[ + 0.17011058031535153`, 0.08842194355669797, 0.22208792588679205`], + Opacity[0.1], SphereBox[{0, 0, 0}, 5701000]}, + {RGBColor[0.7125743490533816, 0.6497192986538929, 0.3566876197914912], + Opacity[0.1], SphereBox[{0, 0, 0}, 3480000]}, + {RGBColor[0.805721607371412, 0.9209617322598633, 0.8867540239707707], + Opacity[0.1], SphereBox[{0, 0, 0}, 1221500]}, + {RGBColor[ + 0.33604086192392857`, 0.9653758014770579, 0.8082011998173695], + Opacity[0.1], CylinderBox[{{0, 0, 6371223}, {0, 0, 6371245}}, 11]}, + {RGBColor[ + 0.3987342262275555, 0.45800648525877863`, 0.0818749710267983], + Opacity[0.1], + CylinderBox[{{0, 0, 6.37122785*^6}, {0, 0, 6.37124015*^6}}, + 3.82]}}], {192., -205.275}, ImageScaled[{0.5, 0.5}], {360, 391}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True], InsetBox[ + Graphics3DBox[{ + {RGBColor[0.5589070573347059, 0.5978054217068245, 0.7989344230742892], + Opacity[0.1], CylinderBox[{{0, 0, 6371223}, {0, 0, 6371245}}, 11]}, + {RGBColor[ + 0.48018618880404773`, 0.0759986548288627, 0.5876392247182236], + Opacity[0.1], + CylinderBox[{{0, 0, 6.37122785*^6}, {0, 0, 6.37124015*^6}}, + 3.82]}}], {576., -205.275}, ImageScaled[{0.5, 0.5}], {360, 391}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True]}, {}}, + ImageSize->{ + UpTo[600], Automatic}, + PlotRange->{{0, 768.}, {-410.55, 0}}, + PlotRangePadding->{6, 5}]], "Output", + CellChangeTimes->{{3.930847372845573*^9, 3.930847385329747*^9}, + 3.9308474913439817`*^9, 3.93084785162274*^9, {3.930847888148061*^9, + 3.9308479131871243`*^9}, 3.930913893848425*^9}, + CellLabel->"Out[7]=",ExpressionUUID->"eb22797f-f400-4d3d-8fbd-b073c6b6fe6c"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"GraphicsRow", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", "}"}]}], "]"}], ",", "\[IndentingNewLine]", + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", + RowBox[{ + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\""}], + "}"}]}], "]"}]}], "}"}], "]"}]], "Input", + CellChangeTimes->{{3.93084793508619*^9, 3.9308479993874903`*^9}}, + CellLabel->"In[8]:=",ExpressionUUID->"c2600694-ff87-45e8-aa45-7e0762877b49"], + +Cell[BoxData[ + GraphicsBox[{{}, {InsetBox[ + Graphics3DBox[{ + {RGBColor[0.3854333319682095, 0.2892742647377069, 0.8094932070709475], + Opacity[0.1], CuboidBox[{-7, -29.1, 6369838}, {7, 29.1, 6369850}]}, + {RGBColor[0.4929121197143007, 0.8219087232339506, 0.05781061434027657], + Opacity[0.1], SphereBox[{0, 0, 0}, 6478000]}, + {RGBColor[ + 0.7307987869233605, 0.8513437211581723, 0.31005230536151496`], + Opacity[0.1], SphereBox[{0, 0, 0}, 6371324]}, + {RGBColor[ + 0.7307987869233605, 0.8513437211581723, 0.31005230536151496`], + Opacity[0.1], SphereBox[{0, 0, 0}, 6356000]}, + {RGBColor[0.1964479760083393, 0.8963638500949567, 0.5827067954797789], + Opacity[0.1], SphereBox[{0, 0, 0}, 6346600]}, + {RGBColor[0.1964479760083393, 0.8963638500949567, 0.5827067954797789], + Opacity[0.1], SphereBox[{0, 0, 0}, 6151000]}, + {RGBColor[0.1964479760083393, 0.8963638500949567, 0.5827067954797789], + Opacity[0.1], SphereBox[{0, 0, 0}, 5971000]}, + {RGBColor[0.1964479760083393, 0.8963638500949567, 0.5827067954797789], + Opacity[0.1], SphereBox[{0, 0, 0}, 5771000]}, + {RGBColor[0.1964479760083393, 0.8963638500949567, 0.5827067954797789], + Opacity[0.1], SphereBox[{0, 0, 0}, 5701000]}, + {RGBColor[ + 0.7506626862180743, 0.5253325729662548, 0.44264896991070013`], + Opacity[0.1], SphereBox[{0, 0, 0}, 3480000]}, + {RGBColor[0.786628620543359, 0.6894960380159745, 0.5808208645108668], + Opacity[0.1], SphereBox[{0, 0, 0}, 1221500]}}], {192., -205.275}, + ImageScaled[{0.5, 0.5}], {360, 391}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True], InsetBox[ + Graphics3DBox[ + {RGBColor[0.7680744425973096, 0.829550618397942, 0.35629067827826955`], + Opacity[0.1], CuboidBox[{-7, -29.1, 6369838}, {7, 29.1, 6369850}]}], {576., -205.275}, + ImageScaled[{0.5, 0.5}], {360, 391}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True]}, {}}, + ImageSize->{ + UpTo[600], Automatic}, + PlotRange->{{0, 768.}, {-410.55, 0}}, + PlotRangePadding->{6, 5}]], "Output", + CellChangeTimes->{{3.930847943143111*^9, 3.930847949745967*^9}, { + 3.930847986070044*^9, 3.930848000974144*^9}, 3.930913896769115*^9}, + CellLabel->"Out[8]=",ExpressionUUID->"f81cbe95-f330-439e-bb2a-5afd2a6189ed"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"GraphicsRow", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", "}"}]}], "]"}], ",", "\[IndentingNewLine]", + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", + RowBox[{ + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\""}], + "}"}]}], "]"}]}], "}"}], "]"}]], "Input", + CellChangeTimes->{{3.930848012050405*^9, 3.930848018088457*^9}}, + CellLabel->"In[9]:=",ExpressionUUID->"bed39234-fab0-4316-b468-3b5240c92626"], + +Cell[BoxData[ + GraphicsBox[{{}, {InsetBox[ + Graphics3DBox[{ + {RGBColor[ + 0.7265281738369762, 0.23654605801110473`, 0.4513669482600129], + Opacity[0.1], SphereBox[{0, 0, 0}, 6478000]}, + {RGBColor[0.8285668626756055, 0.6123076208047464, 0.3756755558784768], + Opacity[0.1], SphereBox[{0, 0, 0}, 6371324]}, + {RGBColor[0.8285668626756055, 0.6123076208047464, 0.3756755558784768], + Opacity[0.1], SphereBox[{0, 0, 0}, 6356000]}, + {RGBColor[ + 0.9161148893717073, 0.27921300626842793`, 0.6143387394926221], + Opacity[0.1], SphereBox[{0, 0, 0}, 6346600]}, + {RGBColor[ + 0.9161148893717073, 0.27921300626842793`, 0.6143387394926221], + Opacity[0.1], SphereBox[{0, 0, 0}, 6151000]}, + {RGBColor[ + 0.9161148893717073, 0.27921300626842793`, 0.6143387394926221], + Opacity[0.1], SphereBox[{0, 0, 0}, 5971000]}, + {RGBColor[ + 0.9161148893717073, 0.27921300626842793`, 0.6143387394926221], + Opacity[0.1], SphereBox[{0, 0, 0}, 5771000]}, + {RGBColor[ + 0.9161148893717073, 0.27921300626842793`, 0.6143387394926221], + Opacity[0.1], SphereBox[{0, 0, 0}, 5701000]}, + {RGBColor[0.05484206767198696, 0.5458487029750789, 0.4295754281321804], + Opacity[0.1], SphereBox[{0, 0, 0}, 3480000]}, + {RGBColor[ + 0.5357009384874412, 0.9447267464809403, 0.19597814354442122`], + Opacity[0.1], SphereBox[{0, 0, 0}, 1221500]}, + {RGBColor[0.2516965757577794, 0.619361342638632, 0.9908314833376688], + Opacity[0.1], + CylinderBox[{{0, 0, 6370644}, {0, 0, 6370704}}, 37]}}], {192., -205.275}, + ImageScaled[{0.5, 0.5}], {360, 391}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True], InsetBox[ + Graphics3DBox[ + {RGBColor[0.3817965102765275, 0.560240567502444, 0.9641154785852963], + Opacity[0.1], CylinderBox[{{0, 0, 6370644}, {0, 0, 6370704}}, 37]}], {576., -205.275}, + ImageScaled[{0.5, 0.5}], {360, 391}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True]}, {}}, + ImageSize->{ + UpTo[600], Automatic}, + PlotRange->{{0, 768.}, {-410.55, 0}}, + PlotRangePadding->{6, 5}]], "Output", + CellChangeTimes->{3.930848021055759*^9, 3.930913899366457*^9}, + CellLabel->"Out[9]=",ExpressionUUID->"b92d23cd-00b8-4bec-934c-3505cc26fb17"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"GraphicsRow", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", "}"}]}], "]"}], ",", "\[IndentingNewLine]", + RowBox[{"Visuals", "[", + RowBox[{"\"\\"", ",", + RowBox[{"{", + RowBox[{ + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\"", ",", + "\"\\"", ",", "\"\\""}], + "}"}]}], "]"}]}], "}"}], "]"}]], "Input", + CellChangeTimes->{{3.930848034137143*^9, 3.930848038040841*^9}, { + 3.930848069237928*^9, 3.930848095894021*^9}}, + CellLabel->"In[10]:=",ExpressionUUID->"8d947ccf-7f8c-4f79-90f3-3e4c0822da4d"], + +Cell[BoxData[ + GraphicsBox[{{}, {InsetBox[ + Graphics3DBox[{ + {RGBColor[0.1452511204960758, 0.1857823034997803, 0.8210871966901028], + Opacity[0.1], SphereBox[{0, 0, 0}, 6478000]}, + {RGBColor[ + 0.009264329959300932, 0.6010366846363033, 0.9777323081774445], + Opacity[0.1], SphereBox[{0, 0, 0}, 6374134]}, + {RGBColor[ + 0.009264329959300932, 0.6010366846363033, 0.9777323081774445], + Opacity[0.1], SphereBox[{0, 0, 0}, 6373934]}, + {RGBColor[0.8885868418790182, 0.0606888172852702, 0.1352350531160138], + Opacity[0.1], SphereBox[{0, 0, 0}, 6371324]}, + {RGBColor[0.8885868418790182, 0.0606888172852702, 0.1352350531160138], + Opacity[0.1], SphereBox[{0, 0, 0}, 6356000]}, + {RGBColor[ + 0.24378268391346758`, 0.9872974256150269, 0.9000805642525402], + Opacity[0.1], SphereBox[{0, 0, 0}, 6346600]}, + {RGBColor[ + 0.24378268391346758`, 0.9872974256150269, 0.9000805642525402], + Opacity[0.1], SphereBox[{0, 0, 0}, 6151000]}, + {RGBColor[ + 0.24378268391346758`, 0.9872974256150269, 0.9000805642525402], + Opacity[0.1], SphereBox[{0, 0, 0}, 5971000]}, + {RGBColor[ + 0.24378268391346758`, 0.9872974256150269, 0.9000805642525402], + Opacity[0.1], SphereBox[{0, 0, 0}, 5771000]}, + {RGBColor[ + 0.24378268391346758`, 0.9872974256150269, 0.9000805642525402], + Opacity[0.1], SphereBox[{0, 0, 0}, 5701000]}, + {RGBColor[ + 0.14910337100201443`, 0.12436833265249714`, 0.4709167434095338], + Opacity[0.1], SphereBox[{0, 0, 0}, 3480000]}, + {RGBColor[ + 0.14682213790113297`, 0.8224387401991233, 0.15218462186818482`], + Opacity[0.1], SphereBox[{0, 0, 0}, 1221500]}, + {RGBColor[ + 0.009264329959300932, 0.6010366846363033, 0.9777323081774445], + Opacity[0.1], CylinderBox[{{0, 0, 6371684}, {0, 0, 6372684}}, 564.19], + {Opacity[0.1], CylinderBox[{{0, 0, -500}, {0, 0, 500}}, 564.19]}}}], {186.03333333333333, -226.8}, + ImageScaled[{0.5, 0.5}], {360, 432}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True], InsetBox[ + Graphics3DBox[ + {RGBColor[0.3678945716866675, 0.40177799905618894`, 0.4900159272427156], + Opacity[0.1], CylinderBox[{{0, 0, 6371684}, {0, 0, 6372684}}, 564.19], + {Opacity[0.1], CylinderBox[{{0, 0, -500}, {0, 0, 500}}, 564.19]}}], {379.1, -226.8}, + ImageScaled[{0.5, 0.5}], {2, 432}, + BaseStyle->{Graphics3DBoxOptions -> {SphericalRegion -> False}}, + ContentSelectable->True]}, {}}, + ImageSize->{ + UpTo[600], Automatic}, + PlotRange->{{0, 386.1333333333333}, {-453.6, 0}}, + PlotRangePadding->{6, 5}]], "Output", + CellChangeTimes->{3.9308480391473513`*^9, 3.930848097045261*^9, + 3.930848147809525*^9, 3.9308483503991203`*^9, 3.930913901638027*^9}, + CellLabel->"Out[10]=",ExpressionUUID->"83fa702f-558a-482f-8040-ab82e433ee6c"] +}, Open ]] +}, Open ]] +}, Open ]] +}, Open ]] +}, +WindowSize->{849, 1027}, +WindowMargins->{{0, Automatic}, {Automatic, -124}}, +FrontEndVersion->"13.3 for Mac OS X ARM (64-bit) (July 24, 2023)", +StyleDefinitions->"Default.nb", +ExpressionUUID->"180ddf86-c472-42f5-b29b-4c869749c105" +] +(* End of Notebook Content *) + +(* Internal cache information *) +(*CellTagsOutline +CellTagsIndex->{} +*) +(*CellTagsIndex +CellTagsIndex->{} +*) +(*NotebookFileOutline +Notebook[{ +Cell[CellGroupData[{ +Cell[580, 22, 169, 3, 67, "Section",ExpressionUUID->"ff56253e-e079-4247-8f6b-78a3d0bd9d71"], +Cell[752, 27, 359, 9, 81, "Text",ExpressionUUID->"2b8b62a5-ee55-4e21-b7e9-9cf4f37c2621"], +Cell[CellGroupData[{ +Cell[1136, 40, 255, 4, 45, "Subsubsection",ExpressionUUID->"1828e000-2b55-4658-a39d-d3ec0e65319c"], +Cell[1394, 46, 237, 4, 35, "Text",ExpressionUUID->"b85c7e71-3fe5-4770-92d2-3b8310f025c5"], +Cell[1634, 52, 390, 8, 30, "Input",ExpressionUUID->"01061750-99ce-4ee2-ba52-6ccdbfbf6f3b"], +Cell[2027, 62, 307, 5, 30, "Input",ExpressionUUID->"a6dc3d65-b39c-4f5d-8e40-9ac1c63e10c6"] +}, Open ]], +Cell[CellGroupData[{ +Cell[2371, 72, 266, 4, 53, "Subsection",ExpressionUUID->"5b67d770-d2b6-44d9-871c-110d0d66dc83"], +Cell[CellGroupData[{ +Cell[2662, 80, 169, 3, 45, "Subsubsection",ExpressionUUID->"b4bb631f-eef9-40d0-aa86-e66af2411a5f"], +Cell[CellGroupData[{ +Cell[2856, 87, 432, 9, 30, "Input",ExpressionUUID->"758f5770-9701-4c2e-9c28-f1655b0a6965"], +Cell[3291, 98, 6350, 106, 294, "Output",ExpressionUUID->"8a3d7fe4-05dd-4061-b3db-44a8b527c5fa"] +}, Open ]], +Cell[CellGroupData[{ +Cell[9678, 209, 372, 8, 30, "Input",ExpressionUUID->"85f2f25e-059b-47e1-b7b7-1cd9a7fb936b"], +Cell[10053, 219, 4820, 88, 212, "Output",ExpressionUUID->"5a905629-0c22-4620-b59c-02655db9bf79"] +}, Open ]], +Cell[CellGroupData[{ +Cell[14910, 312, 456, 10, 30, "Input",ExpressionUUID->"45186205-4e3c-43da-8cf6-942bd6843b62"], +Cell[15369, 324, 6834, 113, 226, "Output",ExpressionUUID->"7aa043b8-c09a-4a90-8f3b-914c780b39b6"] +}, Open ]], +Cell[CellGroupData[{ +Cell[22240, 442, 287, 5, 30, "Input",ExpressionUUID->"623517c9-0b0a-43b7-941c-607589cb47ec"], +Cell[22530, 449, 767, 13, 408, "Output",ExpressionUUID->"928e908c-d246-44d0-8e3f-99106d65c374"] +}, Open ]] +}, Open ]], +Cell[CellGroupData[{ +Cell[23346, 468, 175, 3, 45, "Subsubsection",ExpressionUUID->"7e1653cb-c490-4521-8bec-792810d94dea"], +Cell[CellGroupData[{ +Cell[23546, 475, 1015, 20, 94, "Input",ExpressionUUID->"d7a1ae22-5013-434c-9f40-0ccdc86a3ce7"], +Cell[24564, 497, 2968, 58, 341, "Output",ExpressionUUID->"eb22797f-f400-4d3d-8fbd-b073c6b6fe6c"] +}, Open ]], +Cell[CellGroupData[{ +Cell[27569, 560, 865, 18, 94, "Input",ExpressionUUID->"c2600694-ff87-45e8-aa45-7e0762877b49"], +Cell[28437, 580, 2446, 43, 341, "Output",ExpressionUUID->"f81cbe95-f330-439e-bb2a-5afd2a6189ed"] +}, Open ]], +Cell[CellGroupData[{ +Cell[30920, 628, 864, 18, 94, "Input",ExpressionUUID->"bed39234-fab0-4316-b468-3b5240c92626"], +Cell[31787, 648, 2423, 47, 341, "Output",ExpressionUUID->"b92d23cd-00b8-4bec-934c-3505cc26fb17"] +}, Open ]], +Cell[CellGroupData[{ +Cell[34247, 700, 991, 20, 115, "Input",ExpressionUUID->"8d947ccf-7f8c-4f79-90f3-3e4c0822da4d"], +Cell[35241, 722, 2974, 56, 449, "Output",ExpressionUUID->"83fa702f-558a-482f-8040-ab82e433ee6c"] +}, Open ]] +}, Open ]] +}, Open ]] +}, Open ]] +} +] +*) + diff --git a/resources/Detectors/visuals/DetectorVisuals.wl b/resources/Detectors/visuals/DetectorVisuals.wl new file mode 100644 index 00000000..e27833df --- /dev/null +++ b/resources/Detectors/visuals/DetectorVisuals.wl @@ -0,0 +1,52 @@ +(* ::Package:: *) + +BeginPackage["DetectorVisuals`"]; + +Visuals::usage = "Visuals[Experiment, Modules] displays the experiment as a 3D model without the specified Modules" + +Begin["`Private`"]; + +Visuals[Experiment_, Modules_] := +Module[{exp = Experiment, mods = Modules, dimensionFile, materialFile, keysDim, materials, color, boxes, polygons, cylinders, spheres,numd,visual}, + dimensionFile = Import[FileNames[All, FileNameJoin[{NotebookDirectory[], "/../densities", exp}]][[-1]]]; + materialFile = Import[FileNames[All, FileNameJoin[{NotebookDirectory[], "/../materials", exp}]][[-1]]]; + + dimensionFile = Delete[dimensionFile, Position[dimensionFile, _?(#[[1]] == "#" &)]]; + dimensionFile = Delete[dimensionFile, Position[dimensionFile, {}]]; + dimensionFile = Table[If[MemberQ[d,"#"],Drop[d,{Position[d,"#"]//Flatten//First,-1}],d], {d, dimensionFile}]; + keysDim = Table[If[d[[1]]=="object",d[[Position[d,_?(StringQ[#]&)]//Flatten]][[3]],d[[1]]], {d, dimensionFile}]; + dimensionFile = dimensionFile[[Position[!MemberQ[mods, #] & /@ keysDim, True] // Flatten]]; + + materialFile = Delete[materialFile, Position[materialFile, _?(#[[1]] == "#" &)]]; + materialFile = Delete[materialFile, Position[materialFile, {}]]; + materialFile = Table[If[MemberQ[d,"#"],Drop[d,{Position[d,"#"]//Flatten//First,-1}],d], {d, materialFile}]; + + materialFile = materialFile[[Position[materialFile, _?(StringQ[#[[1]]] &)] // Flatten]]; + materials = Table[m[[1]], {m, materialFile}]; + + color = RandomColor[Length[materials]]; + + boxes = dimensionFile[[Position[dimensionFile, _?(#[[2]] == "box"&&#[[1]]=="object" &)] // Flatten]]; + polygons = dimensionFile[[Position[dimensionFile, _?(#[[2]] == "extr"&&#[[1]]=="object" &)] // Flatten]]; + cylinders = dimensionFile[[Position[dimensionFile, _?(#[[2]] == "cylinder"&&#[[1]]=="object" &)] // Flatten]]; + spheres = dimensionFile[[Position[dimensionFile, _?(#[[2]] == "sphere"&&#[[1]]=="object" &)] // Flatten]]; + + numd = Table[{d[[Position[d,_?(NumberQ[#]&)]//Flatten]],color[[Position[materials,d[[-3]]]//Flatten]]}, {d, polygons}]; + visual = Evaluate[Flatten[{Table[Flatten[{{color[[Position[materials, b[[13]]] // Flatten]], Opacity[0.1]}, + Cuboid[{b[[3]], b[[4]], b[[5]]} - {b[[9]], b[[10]], b[[11]]}/2, {b[[3]], b[[4]], b[[5]]} + {b[[9]], b[[10]], b[[11]]}/2]}], {b, boxes}], + Table[Flatten[{{color[[Position[materials, s[[11]]] // Flatten]], Opacity[0.1]}, + Sphere[{s[[3]], s[[4]], s[[5]]}, s[[9]]]}], {s, spheres}], + Table[Flatten[{{color[[Position[materials, c[[13]]] // Flatten]], Opacity[0.1]}, + Cylinder[{{c[[3]], c[[4]], c[[5]]} - {0, 0, c[[11]]}/2, + {c[[3]], c[[4]], c[[5]]} + {0, 0, c[[11]]}/2}, c[[9]]]}], {c, cylinders}], + Table[Flatten[{{n[[2]], Opacity[0.1]}, + Polyhedron[Table[Table[n[[1]][[1;;3]]+AppendTo[p,z[[1]]], + {p,ArrayReshape[n[[1]][[8;;-11]],{(Length[n[[1]]]-17)/2,2}]*z[[-1]]}], + {z,ArrayReshape[n[[1]][[-9;;-2]],{2,4}]}]]}],{n,numd}] + }]]; + Graphics3D[visual] +]; + +End[]; + +EndPackage[]; diff --git a/resources/Detectors/visuals/README.md b/resources/Detectors/visuals/README.md new file mode 100644 index 00000000..94b38be0 --- /dev/null +++ b/resources/Detectors/visuals/README.md @@ -0,0 +1,25 @@ +# DetectorVisuals.wl + +`DetectorVisuals.wl` is a Wolfram Language package that takes the detector dimension files in `/densities/` and material file in `/materials/` to create a 3D model in `MATHEMATICA`. + +To call functions in the package, save the working notebook in `\visuals\` and add `\path\to\visuals` to `$Path` so `MATHEMATICA` can find the package file. + +If the notebook directory is correctly saved in `\visuals\`, run + +``` +AppendTo[$Path, NotebookDirectory[]]; +``` + +to append `\path\to\visuals\` to `MATHEMATICA` so it has access to the package. Next, run + +``` +<< DetectorVisuals` +``` + +to import the package. To create visualisation for an experiment, run + +``` +Visuals["YourExperiment", {Excluded Modules}] +``` + +to create the 3D model, excluding specified parts of the detector. For more examples, check out `Det_Visuals_Examples.nb`. diff --git a/resources/Examples/Example2/DipolePortal_ND280UPGRD.py b/resources/Examples/Example2/DipolePortal_ND280UPGRD.py new file mode 100644 index 00000000..37def1ab --- /dev/null +++ b/resources/Examples/Example2/DipolePortal_ND280UPGRD.py @@ -0,0 +1,87 @@ +import os + +import siren +from siren.SIREN_Controller import SIREN_Controller + +# Define a DarkNews model +model_kwargs = { + "m4": 0.47, # 0.140, + "mu_tr_mu4": 2.50e-6, # 1e-6, # GeV^-1 + "UD4": 0, + "Umu4": 0, + "epsilon": 0.0, + "gD": 0.0, + "decay_product": "photon", + "noHC": True, + "HNLtype": "dirac", +} + +# Number of events to inject +events_to_inject = 100000 + +# Expeirment to run +experiment = "ND280UPGRD" + +# Define the controller +controller = SIREN_Controller(events_to_inject, experiment) + +# Particle to inject +primary_type = siren.dataclasses.Particle.ParticleType.NuMu + +xs_path = siren.utilities.get_cross_section_model_path(f"DarkNewsTables-v{siren.utilities.darknews_version()}", must_exist=False) +# Define DarkNews Model +table_dir = os.path.join( + xs_path, + "Dipole_M%2.2e_mu%2.2e" % (model_kwargs["m4"], model_kwargs["mu_tr_mu4"]), +) +controller.InputDarkNewsModel(primary_type, table_dir, **model_kwargs) + +# Primary distributions +primary_injection_distributions = {} +primary_physical_distributions = {} + +# energy distribution +flux_file = siren.utilities.get_tabulated_flux_file("T2K_NEAR","PLUS_numu") +edist = siren.distributions.TabulatedFluxDistribution(flux_file, True) +edist_gen = siren.distributions.TabulatedFluxDistribution( + model_kwargs["m4"], 20, flux_file, False +) +primary_injection_distributions["energy"] = edist_gen +primary_physical_distributions["energy"] = edist + +# direction distribution +direction_distribution = siren.distributions.FixedDirection(siren.math.Vector3D(0, 0, 1.0)) +primary_injection_distributions["direction"] = direction_distribution +primary_physical_distributions["direction"] = direction_distribution + +# position distribution +decay_range_func = siren.distributions.DecayRangeFunction( + model_kwargs["m4"], controller.DN_min_decay_width, 3, 240 +) +position_distribution = siren.distributions.RangePositionDistribution( + 5.0, 9.0, decay_range_func, set(controller.GetDetectorModelTargets()[0]) +) +primary_injection_distributions["position"] = position_distribution + +# SetProcesses +controller.SetProcesses( + primary_type, primary_injection_distributions, primary_physical_distributions +) + +controller.Initialize() + +def stop(datum, i): + secondary_type = datum.record.signature.secondary_types[i] + return secondary_type != siren.dataclasses.Particle.ParticleType.N4 + +controller.injector.SetStoppingCondition(stop) + +events = controller.GenerateEvents(fill_tables_at_exit=False) + +os.makedirs("output", exist_ok=True) + +controller.SaveEvents( + "output/ND280UPGRD_Dipole_M%2.2e_mu%2.2e_example" + % (model_kwargs["m4"], model_kwargs["mu_tr_mu4"]), + fill_tables_at_exit=True +) diff --git a/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/FluxCalculator.py b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/FluxCalculator.py new file mode 100644 index 00000000..b9a040c5 --- /dev/null +++ b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/FluxCalculator.py @@ -0,0 +1,35 @@ +import os + +def MakeFluxFile(tag, abs_flux_dir): + + ''' + Accepts the following tags: + {PLUS, MINUS}_{nue,nuebar,numu,numubar} + ''' + + enhance, particle = tag.split("_") + + if enhance not in ["MINUS", "PLUS"]: + print("%s 250kA enhancement specified in tag %s is not valid"%(enhance)) + exit(0) + if particle not in ["numu", "numubar", "nue", "nuebar"]: + print("%s particle specified in tag %s is not valid"%(particle)) + exit(0) + + input_flux_file = os.path.join(abs_flux_dir, + "T2K_%s_250kA.dat"%(enhance)) + + output_flux_file = os.path.join(abs_flux_dir, + "T2KOUT_%s.dat"%(tag)) + + with open(input_flux_file,"r") as fin: + all_lines = fin.readlines() + headers = all_lines[1].strip().split() + data = [line.strip().split() for line in all_lines[3:]] + pid = headers.index(particle) + with open(output_flux_file,"w") as fout: + for row in data: + E, flux = (float(row[1])+float(row[3]))/2, float(row[pid+2]) + flux*=2e-16 # put flux in units of nu/m^2/GeV/POT + print(E, flux, file=fout) + return output_flux_file \ No newline at end of file diff --git a/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2KOUT_PLUS_nue.dat b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2KOUT_PLUS_nue.dat new file mode 100644 index 00000000..eb32d558 --- /dev/null +++ b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2KOUT_PLUS_nue.dat @@ -0,0 +1,138 @@ +0.005 4.4e-09 +0.015 1.998e-08 +0.025 5.12e-08 +0.035 8.76e-08 +0.045 1.392e-07 +0.055 1.912e-07 +0.065 2.5399999999999997e-07 +0.07500000000000001 3.2799999999999997e-07 +0.08499999999999999 4.04e-07 +0.095 4.76e-07 +0.10500000000000001 5.62e-07 +0.11499999999999999 6.36e-07 +0.125 6.96e-07 +0.135 7.78e-07 +0.14500000000000002 8.64e-07 +0.155 9.38e-07 +0.165 1.0019999999999999e-06 +0.175 1.072e-06 +0.185 1.102e-06 +0.195 1.162e-06 +0.20500000000000002 1.2299999999999999e-06 +0.215 1.306e-06 +0.225 1.33e-06 +0.235 1.386e-06 +0.245 1.4639999999999999e-06 +0.255 1.4779999999999999e-06 +0.265 1.52e-06 +0.275 1.568e-06 +0.28500000000000003 1.6219999999999999e-06 +0.295 1.6359999999999999e-06 +0.305 1.6699999999999999e-06 +0.315 1.704e-06 +0.325 1.7319999999999999e-06 +0.335 1.764e-06 +0.345 1.782e-06 +0.355 1.824e-06 +0.365 1.846e-06 +0.375 1.828e-06 +0.385 1.852e-06 +0.395 1.866e-06 +0.405 1.89e-06 +0.415 1.9e-06 +0.425 1.89e-06 +0.435 1.872e-06 +0.445 1.89e-06 +0.455 1.9079999999999998e-06 +0.46499999999999997 1.876e-06 +0.475 1.918e-06 +0.485 1.928e-06 +0.495 1.826e-06 +0.505 1.858e-06 +0.515 1.852e-06 +0.525 1.816e-06 +0.535 1.848e-06 +0.545 1.814e-06 +0.555 1.782e-06 +0.565 1.78e-06 +0.575 1.748e-06 +0.585 1.7459999999999999e-06 +0.595 1.7259999999999999e-06 +0.605 1.696e-06 +0.615 1.6839999999999999e-06 +0.625 1.6639999999999999e-06 +0.635 1.644e-06 +0.645 1.5879999999999999e-06 +0.655 1.618e-06 +0.665 1.572e-06 +0.675 1.5399999999999999e-06 +0.685 1.518e-06 +0.695 1.522e-06 +0.705 1.4979999999999999e-06 +0.715 1.474e-06 +0.725 1.462e-06 +0.735 1.4019999999999999e-06 +0.745 1.394e-06 +0.755 1.386e-06 +0.765 1.3539999999999999e-06 +0.775 1.328e-06 +0.785 1.28e-06 +0.795 1.242e-06 +0.81 1.2639999999999999e-06 +0.83 1.2019999999999999e-06 +0.85 1.172e-06 +0.87 1.09e-06 +0.89 1.084e-06 +0.925 9.78e-07 +0.975 8.86e-07 +1.025 7.88e-07 +1.0750000000000002 7.4e-07 +1.125 6.48e-07 +1.1749999999999998 6.04e-07 +1.25 5.42e-07 +1.35 4.74e-07 +1.45 4.22e-07 +1.55 3.8999999999999997e-07 +1.65 3.6e-07 +1.75 3.38e-07 +1.85 3.2e-07 +1.95 3.06e-07 +2.05 2.96e-07 +2.1500000000000004 2.8e-07 +2.25 2.6599999999999997e-07 +2.3499999999999996 2.52e-07 +2.5 2.36e-07 +2.7 2.1e-07 +2.9 1.848e-07 +3.1 1.622e-07 +3.3 1.408e-07 +3.5 1.234e-07 +3.7 1.03e-07 +3.9 8.9e-08 +4.25 6.78e-08 +4.75 4.56e-08 +5.5 2.68e-08 +6.5 1.328e-08 +7.5 6.8e-09 +8.5 3.44e-09 +9.5 1.868e-09 +10.5 8.38e-10 +11.5 5.06e-10 +12.5 1.786e-10 +13.5 1.162e-10 +14.5 5.48e-11 +15.5 1.042e-11 +16.5 7.42e-12 +17.5 1.508e-12 +18.5 3.06e-12 +19.5 1.286e-12 +20.5 0.0 +21.5 0.0 +22.5 0.0 +23.5 0.0 +24.5 0.0 +25.5 0.0 +26.5 0.0 +27.5 0.0 +28.5 0.0 +29.5 0.0 diff --git a/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2KOUT_PLUS_numu.dat b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2KOUT_PLUS_numu.dat new file mode 100644 index 00000000..b64cdf17 --- /dev/null +++ b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2KOUT_PLUS_numu.dat @@ -0,0 +1,138 @@ +0.005 2.6599999999999997e-07 +0.015 1.426e-06 +0.025 3.06e-06 +0.035 5.02e-06 +0.045 7.06e-06 +0.055 8.92e-06 +0.065 1.0899999999999999e-05 +0.07500000000000001 1.3e-05 +0.08499999999999999 1.6e-05 +0.095 2.06e-05 +0.10500000000000001 2.26e-05 +0.11499999999999999 2.4e-05 +0.125 2.66e-05 +0.135 2.9199999999999998e-05 +0.14500000000000002 3.2e-05 +0.155 3.4399999999999996e-05 +0.165 3.8e-05 +0.175 4.18e-05 +0.185 4.56e-05 +0.195 4.9999999999999996e-05 +0.20500000000000002 5.5e-05 +0.215 6.02e-05 +0.225 6.62e-05 +0.235 7.18e-05 +0.245 7.84e-05 +0.255 8.54e-05 +0.265 9.219999999999999e-05 +0.275 9.98e-05 +0.28500000000000003 0.000108 +0.295 0.00011619999999999999 +0.305 0.0001234 +0.315 0.0001312 +0.325 0.00014099999999999998 +0.335 0.0001484 +0.345 0.000158 +0.355 0.0001664 +0.365 0.0001754 +0.375 0.0001882 +0.385 0.00019559999999999998 +0.395 0.000206 +0.405 0.00021799999999999999 +0.415 0.00022799999999999999 +0.425 0.00023799999999999998 +0.435 0.00025 +0.445 0.00026 +0.455 0.000272 +0.46499999999999997 0.00028199999999999997 +0.475 0.000292 +0.485 0.00030199999999999997 +0.495 0.000314 +0.505 0.00031999999999999997 +0.515 0.00031999999999999997 +0.525 0.00033 +0.535 0.000332 +0.545 0.000338 +0.555 0.000332 +0.565 0.00033 +0.575 0.00033 +0.585 0.000328 +0.595 0.000326 +0.605 0.00031999999999999997 +0.615 0.000316 +0.625 0.000306 +0.635 0.0003 +0.645 0.000296 +0.655 0.000296 +0.665 0.000286 +0.675 0.000274 +0.685 0.000268 +0.695 0.00026199999999999997 +0.705 0.000254 +0.715 0.000246 +0.725 0.000236 +0.735 0.000224 +0.745 0.000216 +0.755 0.000206 +0.765 0.0001944 +0.775 0.000186 +0.785 0.00017659999999999998 +0.795 0.0001672 +0.81 0.0001508 +0.83 0.0001334 +0.85 0.00011619999999999999 +0.87 0.00010339999999999999 +0.89 8.88e-05 +0.925 7.32e-05 +0.975 5.54e-05 +1.025 4.56e-05 +1.0750000000000002 3.74e-05 +1.125 3.24e-05 +1.1749999999999998 2.7799999999999998e-05 +1.25 2.28e-05 +1.35 1.802e-05 +1.45 1.4779999999999999e-05 +1.55 1.258e-05 +1.65 1.092e-05 +1.75 9.42e-06 +1.85 8.22e-06 +1.95 7.2999999999999996e-06 +2.05 6.48e-06 +2.1500000000000004 5.9e-06 +2.25 5.6199999999999996e-06 +2.3499999999999996 4.98e-06 +2.5 4.34e-06 +2.7 3.84e-06 +2.9 3.58e-06 +3.1 3.26e-06 +3.3 3.0799999999999997e-06 +3.5 2.94e-06 +3.7 2.72e-06 +3.9 2.52e-06 +4.25 2.22e-06 +4.75 1.714e-06 +5.5 1.07e-06 +6.5 5.78e-07 +7.5 3.46e-07 +8.5 2.12e-07 +9.5 1.2559999999999999e-07 +10.5 7.72e-08 +11.5 4.28e-08 +12.5 2.68e-08 +13.5 1.46e-08 +14.5 9.08e-09 +15.5 4.14e-09 +16.5 2.84e-09 +17.5 1.91e-09 +18.5 6.62e-10 +19.5 3.26e-10 +20.5 8.06e-11 +21.5 4.22e-11 +22.5 0.0 +23.5 1.74e-12 +24.5 0.0 +25.5 0.0 +26.5 0.0 +27.5 0.0 +28.5 0.0 +29.5 0.0 diff --git a/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2K_MINUS_250kA.dat b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2K_MINUS_250kA.dat new file mode 100644 index 00000000..753552fe --- /dev/null +++ b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2K_MINUS_250kA.dat @@ -0,0 +1,141 @@ + Flux [/cm^2/50MeV/1E21POT] +Bin# Energy numu numubar nue nuebar +---------------------------------------------------------------- +1 0.00 - 0.01 1.35e+09 1.17e+09 5.53e+06 1.59e+07 +2 0.01 - 0.02 7.93e+09 6.65e+09 2.08e+07 9.10e+07 +3 0.02 - 0.03 1.86e+10 1.45e+10 4.67e+07 2.10e+08 +4 0.03 - 0.04 3.16e+10 2.38e+10 7.68e+07 3.94e+08 +5 0.04 - 0.05 4.59e+10 3.33e+10 1.12e+08 6.29e+08 +6 0.05 - 0.06 5.98e+10 4.23e+10 1.44e+08 8.79e+08 +7 0.06 - 0.07 7.08e+10 5.16e+10 1.84e+08 1.18e+09 +8 0.07 - 0.08 7.34e+10 6.07e+10 2.24e+08 1.47e+09 +9 0.08 - 0.09 7.11e+10 7.68e+10 2.60e+08 1.80e+09 +10 0.09 - 0.10 6.68e+10 1.00e+11 2.92e+08 2.19e+09 +11 0.10 - 0.11 6.17e+10 1.09e+11 3.25e+08 2.51e+09 +12 0.11 - 0.12 5.74e+10 1.19e+11 3.62e+08 2.85e+09 +13 0.12 - 0.13 5.18e+10 1.29e+11 3.85e+08 3.20e+09 +14 0.13 - 0.14 4.87e+10 1.41e+11 4.30e+08 3.54e+09 +15 0.14 - 0.15 4.88e+10 1.56e+11 4.58e+08 3.85e+09 +16 0.15 - 0.16 4.57e+10 1.68e+11 4.93e+08 4.16e+09 +17 0.16 - 0.17 4.59e+10 1.83e+11 5.18e+08 4.43e+09 +18 0.17 - 0.18 4.50e+10 2.02e+11 5.59e+08 4.80e+09 +19 0.18 - 0.19 4.47e+10 2.21e+11 5.77e+08 5.04e+09 +20 0.19 - 0.20 4.58e+10 2.44e+11 5.95e+08 5.22e+09 +21 0.20 - 0.21 4.44e+10 2.67e+11 6.28e+08 5.54e+09 +22 0.21 - 0.22 4.45e+10 2.93e+11 6.56e+08 5.79e+09 +23 0.22 - 0.23 4.44e+10 3.17e+11 6.90e+08 5.95e+09 +24 0.23 - 0.24 4.41e+10 3.44e+11 7.06e+08 6.12e+09 +25 0.24 - 0.25 4.46e+10 3.71e+11 7.15e+08 6.43e+09 +26 0.25 - 0.26 4.43e+10 4.04e+11 7.36e+08 6.51e+09 +27 0.26 - 0.27 4.37e+10 4.33e+11 7.59e+08 6.75e+09 +28 0.27 - 0.28 4.41e+10 4.64e+11 7.97e+08 6.85e+09 +29 0.28 - 0.29 4.36e+10 5.02e+11 7.99e+08 7.06e+09 +30 0.29 - 0.30 4.52e+10 5.36e+11 8.05e+08 7.21e+09 +31 0.30 - 0.31 4.44e+10 5.75e+11 8.48e+08 7.33e+09 +32 0.31 - 0.32 4.56e+10 6.09e+11 8.45e+08 7.47e+09 +33 0.32 - 0.33 4.50e+10 6.44e+11 8.50e+08 7.50e+09 +34 0.33 - 0.34 4.63e+10 6.84e+11 8.58e+08 7.69e+09 +35 0.34 - 0.35 4.57e+10 7.15e+11 8.77e+08 7.76e+09 +36 0.35 - 0.36 4.68e+10 7.48e+11 9.03e+08 7.76e+09 +37 0.36 - 0.37 4.82e+10 7.91e+11 8.86e+08 7.86e+09 +38 0.37 - 0.38 4.68e+10 8.23e+11 8.96e+08 7.89e+09 +39 0.38 - 0.39 4.88e+10 8.62e+11 9.11e+08 7.97e+09 +40 0.39 - 0.40 4.79e+10 9.06e+11 9.00e+08 8.04e+09 +41 0.40 - 0.41 4.92e+10 9.47e+11 9.21e+08 7.92e+09 +42 0.41 - 0.42 4.85e+10 9.96e+11 9.29e+08 7.99e+09 +43 0.42 - 0.43 4.83e+10 1.04e+12 9.28e+08 8.01e+09 +44 0.43 - 0.44 4.78e+10 1.08e+12 9.08e+08 8.06e+09 +45 0.44 - 0.45 4.97e+10 1.13e+12 9.17e+08 7.98e+09 +46 0.45 - 0.46 4.94e+10 1.17e+12 9.24e+08 7.97e+09 +47 0.46 - 0.47 4.90e+10 1.22e+12 9.47e+08 7.98e+09 +48 0.47 - 0.48 4.81e+10 1.26e+12 9.25e+08 7.92e+09 +49 0.48 - 0.49 4.87e+10 1.30e+12 9.23e+08 7.90e+09 +50 0.49 - 0.50 5.13e+10 1.34e+12 9.39e+08 7.83e+09 +51 0.50 - 0.51 4.99e+10 1.36e+12 9.40e+08 7.85e+09 +52 0.51 - 0.52 4.96e+10 1.38e+12 9.68e+08 7.70e+09 +53 0.52 - 0.53 4.98e+10 1.39e+12 9.37e+08 7.67e+09 +54 0.53 - 0.54 4.95e+10 1.41e+12 9.35e+08 7.52e+09 +55 0.54 - 0.55 5.09e+10 1.42e+12 9.20e+08 7.44e+09 +56 0.55 - 0.56 5.19e+10 1.41e+12 9.10e+08 7.41e+09 +57 0.56 - 0.57 5.00e+10 1.40e+12 9.34e+08 7.36e+09 +58 0.57 - 0.58 4.90e+10 1.39e+12 9.24e+08 7.25e+09 +59 0.58 - 0.59 4.76e+10 1.37e+12 9.17e+08 7.17e+09 +60 0.59 - 0.60 4.71e+10 1.35e+12 9.07e+08 7.10e+09 +61 0.60 - 0.61 4.83e+10 1.33e+12 9.15e+08 7.17e+09 +62 0.61 - 0.62 4.97e+10 1.31e+12 9.18e+08 7.08e+09 +63 0.62 - 0.63 4.82e+10 1.28e+12 9.06e+08 6.81e+09 +64 0.63 - 0.64 4.57e+10 1.26e+12 9.23e+08 6.71e+09 +65 0.64 - 0.65 4.84e+10 1.23e+12 9.11e+08 6.77e+09 +66 0.65 - 0.66 4.79e+10 1.20e+12 8.84e+08 6.52e+09 +67 0.66 - 0.67 4.48e+10 1.17e+12 8.96e+08 6.37e+09 +68 0.67 - 0.68 4.61e+10 1.15e+12 9.36e+08 6.31e+09 +69 0.68 - 0.69 4.49e+10 1.11e+12 8.93e+08 6.18e+09 +70 0.69 - 0.70 4.58e+10 1.07e+12 9.02e+08 6.03e+09 +71 0.70 - 0.71 4.55e+10 1.03e+12 8.70e+08 5.92e+09 +72 0.71 - 0.72 4.45e+10 9.92e+11 8.71e+08 5.95e+09 +73 0.72 - 0.73 4.48e+10 9.44e+11 8.51e+08 5.83e+09 +74 0.73 - 0.74 4.31e+10 9.08e+11 8.58e+08 5.65e+09 +75 0.74 - 0.75 4.13e+10 8.78e+11 8.50e+08 5.53e+09 +76 0.75 - 0.76 4.23e+10 8.17e+11 8.70e+08 5.43e+09 +77 0.76 - 0.77 4.14e+10 7.74e+11 8.55e+08 5.32e+09 +78 0.77 - 0.78 4.12e+10 7.41e+11 8.31e+08 5.30e+09 +79 0.78 - 0.79 4.10e+10 6.94e+11 8.44e+08 5.12e+09 +80 0.79 - 0.80 3.93e+10 6.50e+11 8.41e+08 5.05e+09 +81 0.80 - 0.82 3.80e+10 5.92e+11 8.46e+08 4.79e+09 +82 0.82 - 0.84 3.71e+10 5.23e+11 8.24e+08 4.60e+09 +83 0.84 - 0.86 3.76e+10 4.56e+11 8.23e+08 4.43e+09 +84 0.86 - 0.88 3.72e+10 3.95e+11 8.07e+08 4.21e+09 +85 0.88 - 0.90 3.56e+10 3.40e+11 7.92e+08 4.05e+09 +86 0.90 - 0.95 3.44e+10 2.76e+11 7.85e+08 3.77e+09 +87 0.95 - 1.00 3.26e+10 2.10e+11 7.63e+08 3.26e+09 +88 1.00 - 1.05 3.08e+10 1.65e+11 7.55e+08 2.88e+09 +89 1.05 - 1.10 2.85e+10 1.36e+11 7.24e+08 2.55e+09 +90 1.10 - 1.15 2.74e+10 1.14e+11 7.12e+08 2.29e+09 +91 1.15 - 1.20 2.60e+10 9.65e+10 6.89e+08 2.07e+09 +92 1.20 - 1.30 2.40e+10 8.08e+10 6.66e+08 1.72e+09 +93 1.30 - 1.40 2.13e+10 6.33e+10 6.38e+08 1.43e+09 +94 1.40 - 1.50 1.93e+10 5.12e+10 6.19e+08 1.22e+09 +95 1.50 - 1.60 1.68e+10 4.25e+10 5.88e+08 1.06e+09 +96 1.60 - 1.70 1.56e+10 3.56e+10 5.60e+08 9.55e+08 +97 1.70 - 1.80 1.39e+10 3.11e+10 5.29e+08 8.68e+08 +98 1.80 - 1.90 1.23e+10 2.70e+10 5.06e+08 8.05e+08 +99 1.90 - 2.00 1.15e+10 2.38e+10 4.78e+08 7.41e+08 +100 2.00 - 2.10 1.02e+10 2.04e+10 4.58e+08 6.86e+08 +101 2.10 - 2.20 9.65e+09 1.80e+10 4.28e+08 6.48e+08 +102 2.20 - 2.30 9.03e+09 1.59e+10 4.15e+08 6.12e+08 +103 2.30 - 2.40 8.04e+09 1.43e+10 3.91e+08 5.66e+08 +104 2.40 - 2.60 7.49e+09 1.24e+10 3.58e+08 5.12e+08 +105 2.60 - 2.80 6.30e+09 1.02e+10 3.19e+08 4.48e+08 +106 2.80 - 3.00 5.46e+09 8.62e+09 2.87e+08 3.83e+08 +107 3.00 - 3.20 4.83e+09 7.44e+09 2.58e+08 3.25e+08 +108 3.20 - 3.40 4.11e+09 6.62e+09 2.30e+08 2.84e+08 +109 3.40 - 3.60 3.53e+09 5.93e+09 2.05e+08 2.41e+08 +110 3.60 - 3.80 3.23e+09 5.18e+09 1.86e+08 2.00e+08 +111 3.80 - 4.00 2.72e+09 4.59e+09 1.59e+08 1.71e+08 +112 4.00 - 4.50 2.33e+09 3.70e+09 1.32e+08 1.27e+08 +113 4.50 - 5.00 1.81e+09 2.55e+09 9.94e+07 8.58e+07 +114 5.00 - 6.00 1.30e+09 1.36e+09 6.39e+07 4.93e+07 +115 6.00 - 7.00 8.07e+08 6.34e+08 3.64e+07 2.47e+07 +116 7.00 - 8.00 5.60e+08 3.26e+08 2.01e+07 1.22e+07 +117 8.00 - 9.00 3.87e+08 1.68e+08 1.15e+07 5.86e+06 +118 9.00 - 10.00 2.68e+08 9.00e+07 6.23e+06 2.85e+06 +119 10.00 - 11.00 1.81e+08 4.11e+07 3.37e+06 1.74e+06 +120 11.00 - 12.00 1.26e+08 1.97e+07 1.81e+06 8.04e+05 +121 12.00 - 13.00 8.56e+07 9.47e+06 1.24e+06 4.68e+05 +122 13.00 - 14.00 5.24e+07 5.40e+06 4.80e+05 1.91e+05 +123 14.00 - 15.00 3.17e+07 3.23e+06 1.99e+05 1.08e+05 +124 15.00 - 16.00 2.32e+07 2.31e+06 1.73e+05 4.54e+04 +125 16.00 - 17.00 1.34e+07 5.33e+05 3.03e+04 7.36e+02 +126 17.00 - 18.00 7.98e+06 1.39e+05 2.11e+04 1.02e+02 +127 18.00 - 19.00 4.69e+06 7.07e+04 1.39e+04 1.30e+03 +128 19.00 - 20.00 1.77e+06 0.00e+00 1.26e+04 0.00e+00 +129 20.00 - 21.00 1.50e+06 0.00e+00 0.00e+00 0.00e+00 +130 21.00 - 22.00 8.03e+05 0.00e+00 0.00e+00 0.00e+00 +131 22.00 - 23.00 1.89e+04 0.00e+00 0.00e+00 0.00e+00 +132 23.00 - 24.00 1.20e+05 0.00e+00 0.00e+00 0.00e+00 +133 24.00 - 25.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +134 25.00 - 26.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +135 26.00 - 27.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +136 27.00 - 28.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +137 28.00 - 29.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +138 29.00 - 30.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 diff --git a/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2K_PLUS_250kA.dat b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2K_PLUS_250kA.dat new file mode 100644 index 00000000..8e20d04b --- /dev/null +++ b/resources/Fluxes/T2K_NEAR/T2K_NEAR-v1.0/T2K_PLUS_250kA.dat @@ -0,0 +1,141 @@ + Flux [/cm^2/50MeV/1E21POT] +Bin# Energy numu numubar nue nuebar +---------------------------------------------------------------- +1 0.00 - 0.01 1.33e+09 1.29e+09 2.20e+07 3.77e+06 +2 0.01 - 0.02 7.13e+09 7.74e+09 9.99e+07 1.71e+07 +3 0.02 - 0.03 1.53e+10 1.82e+10 2.56e+08 3.83e+07 +4 0.03 - 0.04 2.51e+10 3.12e+10 4.38e+08 6.67e+07 +5 0.04 - 0.05 3.53e+10 4.47e+10 6.96e+08 8.79e+07 +6 0.05 - 0.06 4.46e+10 5.96e+10 9.56e+08 1.16e+08 +7 0.06 - 0.07 5.45e+10 7.05e+10 1.27e+09 1.46e+08 +8 0.07 - 0.08 6.50e+10 7.35e+10 1.64e+09 1.71e+08 +9 0.08 - 0.09 8.00e+10 7.19e+10 2.02e+09 2.01e+08 +10 0.09 - 0.10 1.03e+11 6.79e+10 2.38e+09 2.31e+08 +11 0.10 - 0.11 1.13e+11 6.26e+10 2.81e+09 2.60e+08 +12 0.11 - 0.12 1.20e+11 5.75e+10 3.18e+09 2.82e+08 +13 0.12 - 0.13 1.33e+11 5.22e+10 3.48e+09 2.91e+08 +14 0.13 - 0.14 1.46e+11 5.01e+10 3.89e+09 3.21e+08 +15 0.14 - 0.15 1.60e+11 4.81e+10 4.32e+09 3.48e+08 +16 0.15 - 0.16 1.72e+11 4.60e+10 4.69e+09 3.65e+08 +17 0.16 - 0.17 1.90e+11 4.49e+10 5.01e+09 3.87e+08 +18 0.17 - 0.18 2.09e+11 4.38e+10 5.36e+09 3.99e+08 +19 0.18 - 0.19 2.28e+11 4.29e+10 5.51e+09 4.06e+08 +20 0.19 - 0.20 2.50e+11 4.41e+10 5.81e+09 4.27e+08 +21 0.20 - 0.21 2.75e+11 4.26e+10 6.15e+09 4.57e+08 +22 0.21 - 0.22 3.01e+11 4.22e+10 6.53e+09 4.69e+08 +23 0.22 - 0.23 3.31e+11 4.33e+10 6.65e+09 4.76e+08 +24 0.23 - 0.24 3.59e+11 4.26e+10 6.93e+09 5.02e+08 +25 0.24 - 0.25 3.92e+11 4.38e+10 7.32e+09 5.25e+08 +26 0.25 - 0.26 4.27e+11 4.43e+10 7.39e+09 5.21e+08 +27 0.26 - 0.27 4.61e+11 4.23e+10 7.60e+09 5.35e+08 +28 0.27 - 0.28 4.99e+11 4.25e+10 7.84e+09 5.55e+08 +29 0.28 - 0.29 5.40e+11 4.18e+10 8.11e+09 5.65e+08 +30 0.29 - 0.30 5.81e+11 4.34e+10 8.18e+09 5.70e+08 +31 0.30 - 0.31 6.17e+11 4.07e+10 8.35e+09 5.86e+08 +32 0.31 - 0.32 6.56e+11 4.09e+10 8.52e+09 6.14e+08 +33 0.32 - 0.33 7.05e+11 4.18e+10 8.66e+09 5.96e+08 +34 0.33 - 0.34 7.42e+11 4.41e+10 8.82e+09 6.06e+08 +35 0.34 - 0.35 7.90e+11 4.31e+10 8.91e+09 6.14e+08 +36 0.35 - 0.36 8.32e+11 4.36e+10 9.12e+09 6.33e+08 +37 0.36 - 0.37 8.77e+11 4.44e+10 9.23e+09 6.55e+08 +38 0.37 - 0.38 9.41e+11 4.02e+10 9.14e+09 6.49e+08 +39 0.38 - 0.39 9.78e+11 4.40e+10 9.26e+09 6.47e+08 +40 0.39 - 0.40 1.03e+12 4.59e+10 9.33e+09 6.58e+08 +41 0.40 - 0.41 1.09e+12 4.25e+10 9.45e+09 6.74e+08 +42 0.41 - 0.42 1.14e+12 4.63e+10 9.50e+09 6.50e+08 +43 0.42 - 0.43 1.19e+12 4.56e+10 9.45e+09 6.81e+08 +44 0.43 - 0.44 1.25e+12 4.54e+10 9.36e+09 6.59e+08 +45 0.44 - 0.45 1.30e+12 4.43e+10 9.45e+09 6.77e+08 +46 0.45 - 0.46 1.36e+12 4.55e+10 9.54e+09 6.91e+08 +47 0.46 - 0.47 1.41e+12 4.62e+10 9.38e+09 6.88e+08 +48 0.47 - 0.48 1.46e+12 4.54e+10 9.59e+09 6.77e+08 +49 0.48 - 0.49 1.51e+12 4.35e+10 9.64e+09 6.62e+08 +50 0.49 - 0.50 1.57e+12 4.66e+10 9.13e+09 6.78e+08 +51 0.50 - 0.51 1.60e+12 4.38e+10 9.29e+09 6.85e+08 +52 0.51 - 0.52 1.60e+12 4.64e+10 9.26e+09 6.77e+08 +53 0.52 - 0.53 1.65e+12 4.72e+10 9.08e+09 7.28e+08 +54 0.53 - 0.54 1.66e+12 4.41e+10 9.24e+09 6.75e+08 +55 0.54 - 0.55 1.69e+12 4.35e+10 9.07e+09 6.83e+08 +56 0.55 - 0.56 1.66e+12 4.22e+10 8.91e+09 7.04e+08 +57 0.56 - 0.57 1.65e+12 4.53e+10 8.90e+09 7.06e+08 +58 0.57 - 0.58 1.65e+12 4.14e+10 8.74e+09 6.89e+08 +59 0.58 - 0.59 1.64e+12 4.29e+10 8.73e+09 6.95e+08 +60 0.59 - 0.60 1.63e+12 4.58e+10 8.63e+09 7.00e+08 +61 0.60 - 0.61 1.60e+12 4.19e+10 8.48e+09 6.74e+08 +62 0.61 - 0.62 1.58e+12 4.01e+10 8.42e+09 6.77e+08 +63 0.62 - 0.63 1.53e+12 4.12e+10 8.32e+09 6.75e+08 +64 0.63 - 0.64 1.50e+12 4.27e+10 8.22e+09 6.65e+08 +65 0.64 - 0.65 1.48e+12 4.16e+10 7.94e+09 6.66e+08 +66 0.65 - 0.66 1.48e+12 4.10e+10 8.09e+09 6.44e+08 +67 0.66 - 0.67 1.43e+12 4.28e+10 7.86e+09 6.50e+08 +68 0.67 - 0.68 1.37e+12 3.83e+10 7.70e+09 6.82e+08 +69 0.68 - 0.69 1.34e+12 3.83e+10 7.59e+09 6.46e+08 +70 0.69 - 0.70 1.31e+12 4.04e+10 7.61e+09 6.46e+08 +71 0.70 - 0.71 1.27e+12 3.86e+10 7.49e+09 6.55e+08 +72 0.71 - 0.72 1.23e+12 3.79e+10 7.37e+09 6.52e+08 +73 0.72 - 0.73 1.18e+12 3.65e+10 7.31e+09 6.51e+08 +74 0.73 - 0.74 1.12e+12 3.69e+10 7.01e+09 6.27e+08 +75 0.74 - 0.75 1.08e+12 3.60e+10 6.97e+09 6.38e+08 +76 0.75 - 0.76 1.03e+12 3.65e+10 6.93e+09 6.37e+08 +77 0.76 - 0.77 9.72e+11 3.36e+10 6.77e+09 6.33e+08 +78 0.77 - 0.78 9.30e+11 3.36e+10 6.64e+09 6.41e+08 +79 0.78 - 0.79 8.83e+11 3.40e+10 6.40e+09 6.23e+08 +80 0.79 - 0.80 8.36e+11 3.37e+10 6.21e+09 6.40e+08 +81 0.80 - 0.82 7.54e+11 3.16e+10 6.32e+09 6.31e+08 +82 0.82 - 0.84 6.67e+11 3.26e+10 6.01e+09 6.20e+08 +83 0.84 - 0.86 5.81e+11 3.11e+10 5.86e+09 6.03e+08 +84 0.86 - 0.88 5.17e+11 2.91e+10 5.45e+09 6.10e+08 +85 0.88 - 0.90 4.44e+11 2.64e+10 5.42e+09 6.01e+08 +86 0.90 - 0.95 3.66e+11 2.77e+10 4.89e+09 5.98e+08 +87 0.95 - 1.00 2.77e+11 2.55e+10 4.43e+09 5.72e+08 +88 1.00 - 1.05 2.28e+11 2.32e+10 3.94e+09 5.72e+08 +89 1.05 - 1.10 1.87e+11 2.19e+10 3.70e+09 5.48e+08 +90 1.10 - 1.15 1.62e+11 1.94e+10 3.24e+09 5.47e+08 +91 1.15 - 1.20 1.39e+11 1.76e+10 3.02e+09 5.23e+08 +92 1.20 - 1.30 1.14e+11 1.69e+10 2.71e+09 5.21e+08 +93 1.30 - 1.40 9.01e+10 1.45e+10 2.37e+09 4.87e+08 +94 1.40 - 1.50 7.39e+10 1.31e+10 2.11e+09 4.72e+08 +95 1.50 - 1.60 6.29e+10 1.09e+10 1.95e+09 4.40e+08 +96 1.60 - 1.70 5.46e+10 1.01e+10 1.80e+09 4.27e+08 +97 1.70 - 1.80 4.71e+10 8.59e+09 1.69e+09 3.97e+08 +98 1.80 - 1.90 4.11e+10 7.95e+09 1.60e+09 3.75e+08 +99 1.90 - 2.00 3.65e+10 6.75e+09 1.53e+09 3.56e+08 +100 2.00 - 2.10 3.24e+10 6.11e+09 1.48e+09 3.44e+08 +101 2.10 - 2.20 2.95e+10 5.71e+09 1.40e+09 3.28e+08 +102 2.20 - 2.30 2.81e+10 4.98e+09 1.33e+09 2.92e+08 +103 2.30 - 2.40 2.49e+10 4.45e+09 1.26e+09 2.77e+08 +104 2.40 - 2.60 2.17e+10 3.86e+09 1.18e+09 2.53e+08 +105 2.60 - 2.80 1.92e+10 3.31e+09 1.05e+09 2.27e+08 +106 2.80 - 3.00 1.79e+10 2.77e+09 9.24e+08 1.98e+08 +107 3.00 - 3.20 1.63e+10 2.34e+09 8.11e+08 1.72e+08 +108 3.20 - 3.40 1.54e+10 1.93e+09 7.04e+08 1.51e+08 +109 3.40 - 3.60 1.47e+10 1.55e+09 6.17e+08 1.33e+08 +110 3.60 - 3.80 1.36e+10 1.37e+09 5.15e+08 1.14e+08 +111 3.80 - 4.00 1.26e+10 1.24e+09 4.45e+08 9.85e+07 +112 4.00 - 4.50 1.11e+10 8.93e+08 3.39e+08 8.07e+07 +113 4.50 - 5.00 8.57e+09 6.45e+08 2.28e+08 5.93e+07 +114 5.00 - 6.00 5.35e+09 3.86e+08 1.34e+08 3.55e+07 +115 6.00 - 7.00 2.89e+09 2.04e+08 6.64e+07 1.88e+07 +116 7.00 - 8.00 1.73e+09 1.16e+08 3.40e+07 9.78e+06 +117 8.00 - 9.00 1.06e+09 7.44e+07 1.72e+07 5.23e+06 +118 9.00 - 10.00 6.28e+08 4.02e+07 9.34e+06 3.12e+06 +119 10.00 - 11.00 3.86e+08 2.21e+07 4.19e+06 1.49e+06 +120 11.00 - 12.00 2.14e+08 1.16e+07 2.53e+06 9.93e+05 +121 12.00 - 13.00 1.34e+08 6.76e+06 8.93e+05 5.74e+05 +122 13.00 - 14.00 7.30e+07 3.34e+06 5.81e+05 3.20e+05 +123 14.00 - 15.00 4.54e+07 2.46e+06 2.74e+05 4.17e+04 +124 15.00 - 16.00 2.07e+07 2.32e+06 5.21e+04 6.72e+04 +125 16.00 - 17.00 1.42e+07 1.38e+06 3.71e+04 2.21e+02 +126 17.00 - 18.00 9.55e+06 1.42e+06 7.54e+03 2.52e-08 +127 18.00 - 19.00 3.31e+06 4.45e+04 1.53e+04 0.00e+00 +128 19.00 - 20.00 1.63e+06 6.95e+03 6.43e+03 0.00e+00 +129 20.00 - 21.00 4.03e+05 4.53e+05 0.00e+00 0.00e+00 +130 21.00 - 22.00 2.11e+05 0.00e+00 0.00e+00 0.00e+00 +131 22.00 - 23.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +132 23.00 - 24.00 8.70e+03 0.00e+00 0.00e+00 0.00e+00 +133 24.00 - 25.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +134 25.00 - 26.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +135 26.00 - 27.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +136 27.00 - 28.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +137 28.00 - 29.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 +138 29.00 - 30.00 0.00e+00 0.00e+00 0.00e+00 0.00e+00