LMU/musrsim
LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

26.6.2012 Kamil Sedlak

Browse Source 
1) New (GPS and others) examples added.
2) The documentation updated accordingly.
This commit is contained in:
sedlak 2012-06-26 14:42:44 +00:00
parent 6d467b024c
commit f2b9079625
21 changed files with 3149 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
doc/musrSim.pdf
View File

Binary file not shown.

24
doc/musrSim.tex
View File

@ -1542,6 +1542,30 @@ then activated in the macro file using the command \\[2ex]
%
The results of this GPD simulation are described in the \emph{musrSimAna} manual~\cite{musrSimAna}.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Example 5 -- GPS instrument}
For the details about the GPS simulation see the documents
``Simulation of the GPS μSR instrument 1--4'' (currently saved as
/afs/psi.ch/project/lmu/Facility/musr\_simulations/documentation/GPS/*.pdf).
\begin{itemize}
\item {\tt 50121.mac} -- GPS as installed in reality in 2011.
\item {\tt 50131.mac} -- GPS planned upgrade; bottle-shaped forward veto.
\item {\tt 50161.mac, 50171.mac, 50181.mac} -- GPS planned upgrade; pyramidal forward veto.
\end{itemize}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Other Examples}
During the years a lot of ``*.mac'' files were created and used.
Most of these files are stored in the file ``run\_many\_files.tar.gz''.
It is not guaranteed that all the ``*.mac'' files stored there are
compatible with the newest version of musrSim/Geant4. However,
they can be still useful, at least as a source of inspiration.
See the file ``README.TXT'' for a short description of what the purpose
of the different files was.
Sometimes the runs need additional input files (e.g.\ a field map
,TURTLE file, ...). These files are too large to be stored in the svn repository,
and typically can be found in subdirectories of the high field project:
/afs/psi.ch/project/HighFieldMuSR/.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\clearpage

BIN
doc/musrSimAna.pdf
View File

Binary file not shown.

23
doc/musrSimAna.tex
View File

@ -772,7 +772,28 @@ The purpose of the example analysis of the GPD simulation was to illustrate
the potential of \musrSim\ and \musrSimAna\ programs to investigate features
like time-independent background, sample muon fraction, prompt peak, \ldots
This information can be used in design and optimisation of \musr\ instruments.
%========================================================================================================
\section{GPS instrument}
%
It is foreseen that GPS instrument could be arranged in two geometries
after the upgrade (depending from which side the calorimeter would be
inserted).
\begin{itemize}
\item {\tt 50130hb.v1190} -- Calorimeter inserted from one side.
\item {\tt 50130hl.v1190} -- Calorimeter inserted from the other side.
\item {\tt 50130hb1.v1190 -- 50130hb6.v1190} -- All positron counters
analysed individually.
\end{itemize}
See the document about the GPS simulations saved in the directory: \\
/afs/psi.ch/project/lmu/Facility/musr\_simulations/documentation/GPS/ \\
for more details.
%========================================================================================================
\section{Other Examples}
Many different ``*.v1190'' files are stored in the file:
``run\_musrSimAna\_many\_files.tar.gz''. They could serve as additional examples.
Note that the syntax of the ``fit'' command was changed
at some point, and therefore the ``fit'' command might cause problems
(the {\tt ``option''} has to be added in the old ``*.v1190'' files).
%========================================================================================================
\begin{thebibliography}{0}

170
musrSimAna/50130hb.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
! 10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

170
musrSimAna/50130hb1.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
! 2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
! 3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
! 4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
! 5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
! 6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
! 7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
! 8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
! 9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
! 10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

170
musrSimAna/50130hb2.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
! 1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
! 3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
! 4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
! 5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
! 6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
! 7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
! 8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
! 9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
! 10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

170
musrSimAna/50130hb3.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
! 1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
! 2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
! 4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
! 5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
! 6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
! 7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
! 8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
! 9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
! 10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

170
musrSimAna/50130hb4.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
! 1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
! 2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
! 3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
! 5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
! 6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
! 7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
! 8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
! 9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
! 10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

170
musrSimAna/50130hb5.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
! 1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
! 2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
! 3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
! 4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
! 8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
! 9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
! 10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

170
musrSimAna/50130hb6.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
! 1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
! 2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
! 3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
! 4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
! 5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
! 6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
! 7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

170
musrSimAna/50130hl.v1190 Normal file
View File

@ -0,0 +1,170 @@
# TDC V1190 Set up file generated by Deltat : 16-Apr-2008 10:30:27
# bins are always in TDC channel units (195.3125 ps/channel)
INSTRUMENT=GPS
DESCRIPTION=No Veto -- 1port
TYPE=TDCV1190
RESOLUTION=195.3125
MDELAY=0
PDELAY=2000
#REWINDTIMEBINS=524288
REWINDTIMEBINS=-1000000000
MCOINCIDENCEW=50
PCOINCIDENCEW=50
VCOINCIDENCEW=100
MUONRATEFACTOR=0.0905048
#MUONRATEFACTOR=1000
#INFINITELYLOWMUONRATE
DATAWINDOWMIN=-0.5
DATAWINDOWMAX=10.0
PILEUPWINDOWMIN=-10.5
PILEUPWINDOWMAX=10.5
102; "M up"; M; 0.4; 800; -51 -52;
1; "Forw"; P; 0.1; 800; -51 -52; F; 1485; 1515; 50995;
2; "Backw"; P; 0.1; 800; -51 -52; B; 1485; 1515; 50995;
3; "Up"; P; 0.1; 800; -51 -52; U; 1485; 1515; 50995;
4; "Down"; P; 0.1; 800; -51 -52; D; 1485; 1515; 50995;
5; "R1"; P; 0.1; 800; -51 -52; R1; 1485; 1515; 50995;
6; "R2"; P; 0.1; 800; -51 -52; R2; 1485; 1515; 50995;
! 7; "R3"; P; 0.1; 800; -51 -52; R3; 1485; 1515; 50995;
8; "L1"; P; 0.1; 800; -51 -52; L1; 1485; 1515; 50995;
9; "L2"; P; 0.1; 800; -51 -52; L2; 1485; 1515; 50995;
10; "L3"; P; 0.1; 800; -51 -52; L3; 1485; 1515; 50995;
51; "Fveto"; V; 0.02; 800;
52; "Bveto"; V; 0.02; 800;
!31; Ch31; N;
counterPhaseShifts 1 0 2 180
counterPhaseShifts 3 0 4 180
counterGrouping F 1
counterGrouping B 2
counterGrouping U 3
counterGrouping D 4
counterGrouping R 5 6 7
counterGrouping L 8 9 10
sampleID 201 202
musrTH1D hmuIniPosZ "Muon initial position;z (mm);N" 200 -100000. 100000. muIniPosZ
musrTH1D hMuTargetPol_Theta "Muon polarisation in the sample;#theta (deg);N" 300 60. 120. muTargetPol_Theta
musrTH1D hMuTargetPol_Phi "Muon polarisation in the sample;#phi (deg);N" 1800 -180. 180. muTargetPol_Phi
musrTH1D hMuTargetPol_Phi360 "Muon polarisation in the sample;#phi+360 (deg);N" 1800 0. 360. muTargetPol_Phi360
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 -180. 180. pos_Phi_MINUS_muDecayPol_Phi
musrTH1D hPos_Phi_MINUS_muDecayPol_Phi360 ";#phi_{e ini} - #phi_{#mu decay spin} (deg);N" 360 0. 360. pos_Phi_MINUS_muDecayPol_Phi360
musrTH1D hdetP_phi_MINUS_pos_Phi ";#phi_{det} - #phi_{e ini} (deg);N" 360 -180. 180. detP_phi_MINUS_pos_Phi
musrTH1D hdetP_phi_MINUS_pos_Phi360 ";#phi_{det} - #phi_{e ini} (deg);N" 360 0. 360. detP_phi_MINUS_pos_Phi360
musrTH1D hdet_time1_MINUS_muDecayTime ";TOF_{e} + TDC smearing;N" 500 0. 500. det_time1_MINUS_muDecayTime
musrTH1D hdetP_time_start_MINUS_muDecayTime ";TOF_{e};N" 500 0. 500. detP_time_start_MINUS_muDecayTime
musrTH2D hMuDecayMap "Muon decay map;z [mm]; r[mm]" 100 -80. 20. 40 0. 20. muDecayPosZ muDecayPosR
musrTH1D hMuDecayPosZall "Z of muons when they decay;z[mm];N" 4002 -2001 2001 muDecayPosZ
musrTH1D hMuDecayPosZ "Penetration of muons into the sample;z[mm];N" 100 -5.0 5. muDecayPosZ
musrTH1D hMuDecayPosX "X of decayed muons;x[mm];N" 100 -5.0 5. muDecayPosX
musrTH1D hPosMom "Momentum of the produced positrons;P [MeV];N" 55 0. 55. pos_Momentum
musrTH1D hdet_m0edep "Energy deposited in the trigger;E[MeV];N" 100 0. 2 det_m0edep
musrTH1D hdet_posEdep "Energy deposited in the positron Counters;E[MeV];N" 100 0. 2 det_posEdep
#musrTH1D hdet_time0 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time0
#musrTH1D hdet_time1 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time1
#musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 220 -11. 11. det_time10
musrTH1D hdet_time10 "detected muSR spectra;time[#mus];N" 240 -1. 11. det_time10
musrTH1D hdet_time10_bgr1 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time10_bgr3 "detected muSR spectra (phase shifted);time[#mus];N" 280 -3. 11. det_time10
musrTH1D hdet_time31 "double hit time difference;time[#mus];N" 280 -3. 11. det_time31
musrTH1D hdet_time20 "detected muSR spectra (phase shifted);time[#mus];N" 240 -1. 11. det_time20
musrTH1D hgen_time10 "generated muSR spectra;time[#mus];N" 240 -11. 11. gen_time10
musrTH1D hdet_time10_MINUS_gen_time10 "detected - generated ;time_{det}-time_{gen} [ps];N" 100 -500. 500. det_time10_MINUS_gen_time10
#musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 4000 -2000. 2000. muDecayDetID
musrTH1D hMuDecayDetID "Detector ID where the muons decay;Detector ID; N" 2002 -1001. 1001. muDecayDetID
musrTH1D hMuDecayDetIDpileup "Detector ID where the pileup muons decay;Detector ID; N" 2002 -1001. 1001. pileup_muDecayDetID
humanDecayHistograms hMuDecayDetID hMuDecayDetIDpileup 201 Sample 202 Sample 102 M 1 F 2 B 3 U 4 D 5 R 6 R 7 R 8 L 9 L 10 L 51 FV 52 BV 401 Cryo 402 Cryo 403 Cryo 404 Cryo 411 Cryo 412 Cryo 413 Cryo 414 Cryo 415 Cryo 8201 pipe 8251 pipe 8301 pipe -1 World -1000 Escaped
musrTH1D hdet_time10_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time10 rotreference 66 0
musrTH1D hdet_time10_exp "hdet_time10_exp;time[#mus];N" 240 -1.0 11. det_time10 correctexpdecay
musrTH1D hdet_time20_rotref "rot ref frame;time[#mus];N" 240 -1.0 11. det_time20 rotreference 66 0
#musrTH1D hdet_time20_rotrefTwo "rot ref frame 2;time[#mus];N" 240 -1.0 11. det_time20 rotreference -1353 0
musrTH2D hMuTargetPolPhi_pInit ";#phi_{#mu spin} (deg); p_{#mu init} (MeV/c)" 100 0. 360. 100 25.7 28.9 muTargetPol_Phi360 muIniMomZ
musrTH2D hMuTargetPolPhi_RInit ";#phi_{#mu spin} (deg); R_{#mu init} (mm)" 100 0. 360. 100 0. 40. muTargetPol_Phi360 muIniPosR
musrTH2D hpos_detID_doubleHit "Double hits;positron counter ID; positron counter ID" 11 0. 11. 11 0. 11. pos_detID pos_detID_doubleHit
condition 0 alwaysTrue
condition 1 oncePerEvent
condition 2 muonDecayedInSample_gen
condition 3 muonTriggered_gen
condition 4 muonTriggered_det
condition 5 positronHit_det
condition 6 goodEvent_det
condition 7 goodEvent_gen
condition 8 goodEvent_det_AND_goodEvent_gen
condition 9 pileupEvent
condition 10 goodEvent_det_AND_muonDecayedInSample_gen
condition 11 goodEvent_F_det
condition 12 goodEvent_B_det
condition 13 goodEvent_U_det
condition 14 goodEvent_D_det
condition 15 goodEvent_R_det
condition 16 goodEvent_L_det
condition 17 goodEvent_F_det_AND_pileupEvent
condition 18 goodEvent_B_det_AND_pileupEvent
condition 19 goodEvent_U_det_AND_pileupEvent
condition 20 goodEvent_D_det_AND_pileupEvent
condition 21 goodEvent_R_det_AND_pileupEvent
condition 22 goodEvent_L_det_AND_pileupEvent
condition 23 doubleHit
fit hMuTargetPol_Phi gaus "" -180. 180. 1000. 0. 10.
fit hMuTargetPol_Phi360 gaus "" 0. 360. 1000. 0. 10.
fit hMuTargetPol_Theta gaus "" 60. 120. 1000. 0. 10.
fit hdet_time20 simpleExpoPLUSconst "" 0.1 9.9 100 1
fit hdet_time10_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
fit hdet_time20_rotref rotFrameTime20 "" 0.2 9.8 0 0 100
#fit hdet_time10 simpleExpoPLUSconst "" 0.2 9.8 100 1
fit hdet_time10 funct4 "" 0.2 9.8 0 0 0.2 600 0
fit hdet_time10_exp TFieldCosPLUSbg "" 0.2 9.8 1 0 0.3 300 1
fit hdet_time10_bgr1 pol0 "" 0.2 9.8
#fit hdet_time10_bgr3 pol0 "" -1.9 -0.1
fit hdet_time10_bgr3 pol0 "WW" -0.5 -0.01
#fit hdet_time20_rotrefTwo rotFrameTime20 "" 0.1 9.9 0 0 100
#draw hMuDecayPosZ 0
#draw hMuDecayPosZall 1
#draw hMuDecayPosZ 1
#draw hMuDecayPosZ 2
#draw hMuDecayPosZ 3
#draw hMuDecayPosZ 4
#draw hMuDecayMap 0
#draw hMuDecayMap 1
#draw hdet_m0edep 0
#draw hdet_m0edep 1
#draw hdet_m0edep 2
#draw hdet_m0edep 3
#draw hdet_m0edep 4
#draw hdet_m0edep 5
#draw hdet_m0edep 6
draw hdet_time10 6
draw hgen_time10 7
draw hdet_time20 6
draw hgen_time20 7
draw hdet_time20_rotref 6
draw hdet_time20_bgr 6
#draw hdet_time20_rotrefTwo 6
draw hdet_time10_MINUS_gen_time10 8
#draw hdet_posEdep 5
#draw hMuDecayDetID 1
#draw hdet_m0edep 9
#draw hdet_time10 9
#draw hMuDecayDetIDpileup 1
#draw hMuDecayDetIDpileup 9
#draw humanDecayHistograms 1
#draw humanDecayHistograms 6
#draw humanDecayHistograms 9
#draw humanDecayPileupHistograms 1
#draw humanDecayPileupHistograms 6
#draw humanDecayPileupHistograms 9
$

BIN
musrSimAna/run_musrSimAna_many_files.tar.gz Normal file
View File

Binary file not shown.

344
run/50121.mac Normal file
View File

@ -0,0 +1,344 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 0 0 45
/musr/command rotation matrix2 0 180 0
/musr/command rotation matrix3 0 270 0
#/musr/command rotation matrix3 90 90
/musr/command arrayDef zPlaneGPSveto 12 0. 5. 91. 105. 115. 125. 135. 145. 155. 168.49 168.5 172.5
/musr/command arrayDef rInnerGPSveto 12 10. 11. 31. 31. 28.75 24.8. 19.7 14 10.3 9. 0. 0.
/musr/command arrayDef rOuterGPSveto 12 12.5 15.11 35.11 35. 32.85 29.1 24.2 18.6 14.6 13.2 13.0 13.0
# Rotation "fieldRot 0 0.57295 0" corresponds to 10 mrad (2.5mm/25cm) tilt of the magnetic field
#/musr/command rotation fieldRot 0 0.57295 0
#/musr/command rotation fieldRot 0 1.1459 0
# WORLD VOLUME
/musr/command construct box World 300 300 2000 G4_Galactic 0 0 0 no_logical_volume norot dead -1
# ---- BEAMPIPE vers. 1 (ALC-like)
#/musr/command construct tubs beampipe 0 34 810 0 360 G4_Galactic 0 0 -840 log_World norot dead 232
#/musr/command construct tubs beampipeshell 34 36 810 0 360 G4_Al 0 0 -840 log_World norot dead 231
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 16.5 34 25 0 360 G4_Pb 0 0 540 log_beampipe norot dead 241
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 34. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 2 (SEGMENTED and BROAD, small section 30cm long)
#/musr/command construct tubs beampipeA 0 128 500 0 360 G4_Galactic 0 0 -800 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 500 0 360 G4_Al 0 0 -800 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 135 0 360 G4_Galactic 0 0 -165 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 135 0 360 G4_Al 0 0 -165 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs pbCollimator_1 34 128 25 0 360 G4_Pb 0 0 475 log_beampipeA norot dead 241
#/musr/command construct tubs pbCollimator_2 16.5 34 25 0 360 G4_Pb 0 0 65 log_beampipeB norot dead 242
# ACTIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs activeCollim 2.5 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 3 (SEGMENTED and BROAD, small section 50cm long)
#/musr/command construct tubs beampipeA 0 128 400 0 360 G4_Galactic 0 0 -900 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 400 0 360 G4_Al 0 0 -900 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 230 0 360 G4_Galactic 0 0 -270 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 230 0 360 G4_Al 0 0 -270 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 20 128 15 0 360 G4_Pb 0 0 385 log_beampipeA norot dead 241
#/musr/command construct cons pbCollimator_2 122 128 34 128 70 0 360 G4_Pb 0 0 300 log_beampipeA norot dead 242
#/musr/command construct tubs pbCollimator_11 2.5 34 15 0 360 G4_Pb 0 0 215 log_beampipeB norot dead 245
#/musr/command construct cons pbCollimator_12 32 34 2.5 34 25 0 360 G4_Pb 0 0 175 log_beampipeB norot dead 246
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 5. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 150
#/musr/command construct tubs activeCollimB 5. 10. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 151
#/musr/command construct tubs activeCollimC 10. 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 152
#
# OLD ALC BEAMPIPE
# OLD ALC DESIGN: MSR 20.05.082 Tube
/musr/command construct tubs tube 35.45 38.05 315 0 360 Steel 0 0 -613 log_World norot dead 8201
/musr/command construct tubs tube_a 0 35.45 315 0 360 G4_Galactic 0 0 -613 log_World norot dead 8251
# SECOND COLLIMATOR
/musr/command construct tubs shield_collimatorA2 10 35 25 0 360 G4_Cu 0 0 290 log_tube_a norot dead 341
# OLD ALC DESIGN: MSR 20.05.083 Flange
/musr/command construct tubs flange 38.05 90 6 0 360 Steel 0 0 -934 log_World norot dead 8301
# ----
# POSITRON COUNTERS
/musr/command construct GPSforward forwardCounter 25 25 2.5 12.51 15.12 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward backwardCounter 30 30 2.5 8.80 7.64 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -30.0 log_World norot musr/ScintSD 2
/musr/command construct box up 2.5 40 45 G4_PLASTIC_SC_VINYLTOLUENE 40 0 0. log_World norot musr/ScintSD 3
/musr/command construct box down 2.5 40 45 G4_PLASTIC_SC_VINYLTOLUENE -40 0 0. log_World norot musr/ScintSD 4
/musr/command construct box rite 45 2.5 60 G4_PLASTIC_SC_VINYLTOLUENE 0 -54.5 0 log_World norot musr/ScintSD 5
#/musr/command construct box up 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE 57.5 0 0.log_World norot musr/ScintSD 3
#/musr/command construct box down 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE -57.5 0 0.log_World norot musr/ScintSD 4
#/musr/command construct box rite1 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 -38 log_World norot musr/ScintSD 5
#/musr/command construct box rite2 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 38 log_World norot musr/ScintSD 6
#/musr/command construct box rite3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 -51.25 0 log_World norot musr/ScintSD 7
#/musr/command construct box left1 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 -48.3 log_World norot musr/ScintSD 8
#/musr/command construct box left2 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 48.3 log_World norot musr/ScintSD 9
#/musr/command construct box left3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 53.75 0 log_World norot musr/ScintSD 10
# FORWARD VETO
/musr/command construct polyconeA forwardVeto 0 360 12 zPlaneGPSveto rInnerGPSveto rOuterGPSveto G4_PLASTIC_SC_VINYLTOLUENE 0 0 24 log_World norot musr/ScintSD 51
# BACKWARD VETO
/musr/command construct GPSbackwardVeto backwardVeto 35.13 7.63 35.13 7.63 60 31 3.5 31 3.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -87.5 log_World norot musr/ScintSD 52
#---------------------
# OUTER SHIELD
/musr/command construct tubeWithWindows tube1 20 21 33.5 14.5 33.5 14.5 33.5 G4_Cu 0 0 0 log_World matrix3 dead 401
/musr/command construct tubs bottom1 0 20 0.5 0 360 G4_Cu 0 -33.0 0 log_World matrix3 dead 402
/musr/command construct tubs window1 21 21.010 15 0 360 G4_Cu 0 0 0 log_World matrix3 dead 403
/musr/command construct tubs bottom1a 0 20 0.5 0 360 G4_Cu 0 33.0 0 log_World matrix3 dead 404
# INNER SHIELD
/musr/command construct tubeWithWindows tube2 11.5 12.5 21.5 7.5 21.5 10. 21.5 G4_Cu 0 0 0 log_World matrix3 dead 411
/musr/command construct tubs bottom2 0 11.5 3.0 0 360 G4_Cu 0 -21.5 0 log_World matrix3 dead 412
/musr/command construct tubs window2 12.5 12.510 8 0 360 G4_Cu 0 0 0 log_World matrix3 dead 413
/musr/command construct tubs bottom2a 0 11.5 3.0 0 360 G4_Cu 0 21.5 0 log_World matrix3 dead 414
# HEATER
/musr/command construct tubs heater 16 16.3 7.5 0 360 G4_Cu 0 -18.5 0 log_World matrix3 dead 415
#
#---------------------
# MAGNET
#/musr/command construct tubs magnet 44.4 100 150 0 360 G4_He 0 0 0 log_World norot dead 221
# MAGNET WALL
#/musr/command construct tubs mag_wall 44.4 46.5 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 222
# MAGNET VOLUME TO DELETE TRACKS
#/musr/command construct tubs sh0 46.5 100 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 223
# SHIELD TO DELETE TRACKS
#/musr/command construct tubs shield1 36 100 5 0 360 G4_AIR 0 0 -105 log_World norot dead -2
#/musr/command construct tubs sh2 0 100 5 0 360 G4_AIR 0 0 205 log_World norot dead -3
# NOSE EXTENSION
#/musr/command construct tubs noseExtensionA 15 25 2.0 0 360 G4_Al 0 0 -47 log_World norot dead 300
#/musr/command construct tubs noseExtensionB 25 28 126.5 0 360 G4_Al 0 0 -171.5 log_World norot dead 302
#/musr/command construct tubs noseExtensionC 28 41 4.0 0 360 G4_Al 0 0 -294 log_World norot dead 304
#
# COLLIMATOR
#/musr/command construct tubs collimatorA 2.5 25 15 0 360 G4_Pb 0 0 -64 log_World norot dead 310
#/musr/command construct tubs collimatorB 3.5 25 10 0 360 G4_Pb 0 0 -89 log_World norot dead 312
#/musr/command construct tubs collimatorC 5.0 25 10 0 360 G4_Pb 0 0 -109 log_World norot dead 314
#
# M COUNTERS AND M COUNTER HOLDER
#/musr/command construct tubs M0_holder1 5.5 7.5 3.25 0 360 G4_Al 0 0 -36.31 log_World norot dead 330
#/musr/command construct tubs M0_holder2 3.5 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -37.06 log_World norot dead 333
#/musr/command construct tubs M0_holder3 4.0 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -35.56 log_World norot dead 335
/musr/command construct tubs M0 0 7.5 0.15 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -25.0 log_World norot musr/ScintSD 102
#/musr/command construct tubs M0_electronics 4.0 9 0.85 0 360 G4_POLYCARBONATE 0 0 -40.41 log_World norot dead 337
#
# TARGET SPACE
#/musr/command construct tubs targetspace 0 5 50 0 360 G4_He 0 0 -20 log_cryostat
# SCINTILLATOR BEFORE TARGET
#/musr/command construct tubs coulombM1 0 3.0 0.1 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2 log_World norot musr/ScintSD 101
# TARGET
/musr/command construct tubs target 0 4.0 0.215 0 360 G4_Ag 0 0 0 log_World norot dead 201
/musr/command construct tubs targetFieldVol 0 0.5 0.015 0 360 G4_Ag 0 0 0 log_target norot dead 202
#/musr/command construct tubs vetoTarget 0 5.0 1. 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 1.220 log_World norot musr/ScintSD 160
#/musr/command construct tubs vetoCyl 5 6.0 3.6 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -1.39 log_World norot musr/ScintSD 161
#/musr/command construct TubeWithTubeHole vetoTarget 7 9 5.1 0 360 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct TubeWithHoleAndTubeHole vetoTarget 1.5 6 5.1 0 360 4 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct tubs sampleHolder 0 4 31.39 0 360 G4_Ag 0 0 31.61 log_World norot dead 165
#/musr/command construct tubs sampleHolder2 0 7.5 108.5 0 360 G4_Ag 0 0 171.5 log_World norot dead 166
#
#/musr/command construct TubeWithHolePlusTubeHole PlexyCyl1 4 12.5 4.5 0 360 7 12.5 5 G4_PLEXIGLASS 0 0 -10. log_World matrix2 dead 45
#/musr/command construct tubs PlexyCyl2 7.5 12.5 2 0 360 G4_PLEXIGLASS 0 0 12.5 log_World norot dead 46
#/musr/command construct TubeWithHolePlusTubeHole vetoCyl 2.5 12. 5. 0 360 4 12 9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -10.5 log_World matrix2 musr/ScintSD 51
#/musr/command construct tubs vetoCylA 8.0 12.0 10 91 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 61
#/musr/command construct tubs vetoCylB 8.0 12.0 10 181 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 62
#/musr/command construct tubs vetoCylC 8.0 12.0 10 271 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 63
#/musr/command construct tubs vetoCylD 8.0 12.0 10 1 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 64
#
#---------------------------------------------------------
#/musr/command region define goulombRegion log_target
#/musr/command region define goulombRegion log_M0
#/musr/command region setProductionCut goulombRegion 0.01 0.01 0.01
#---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
#---------------------------------------------------------
#
/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
/musr/command visattributes G4_PLEXIGLASS yellow
/musr/command visattributes G4_Pb blue
#/musr/command visattributes log_beampipe invisible
#/musr/command visattributes log_beampipeAshell invisible
#/musr/command visattributes log_beampipeBshell invisible
/musr/command visattributes log_World invisible
/musr/command visattributes log_target yellow
#/musr/command visattributes log_magnet invisible
#/musr/command visattributes log_mag_wall invisible
#/musr/command visattributes log_magnet yellow
#/musr/command visattributes log_mag_wall yellow
#/musr/command visattributes log_sh0 invisible
#/musr/command visattributes log_sh2 invisible
#/musr/command visattributes log_vetoTarget green
#/musr/command visattributes log_vetoCyl green
###################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPS_6kG_yAx_xRad_cgs.table log_targetFieldVol 0.005
/musr/command globalfield setparameter SetLargestAcceptableStep 2
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 -10
/musr/command globalfield printFieldValueAtPoint 0 10 10
/musr/command globalfield printFieldValueAtPoint 10 0 10
/musr/command globalfield printFieldValueAtPoint -10 0 10
/musr/command globalfield printFieldValueAtPoint 10 10 10
/musr/command globalfield printFieldValueAtPoint 40 40 100
/musr/command globalfield printFieldValueAtPoint 40 40 900
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
#/musr/command process addModel mu+ G4MuMultipleScattering G4UrbanMscModel92 0
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
#
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
#/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 0.1
#
/musr/run/howOftenToPrintEvent 100000
/musr/command maximumRunTimeAllowed 86000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_edep_el off
/musr/command rootOutput det_edep_pos off
/musr/command rootOutput det_edep_gam off
/musr/command rootOutput det_edep_mup off
/musr/command rootOutput det_nsteps off
/musr/command rootOutput det_length off
/musr/command rootOutput det_time_end off
/musr/command rootOutput det_kine off
/musr/command rootOutput det_VrtxKine off
/musr/command rootOutput det_VrtxX off
/musr/command rootOutput det_VrtxY off
/musr/command rootOutput det_VrtxZ off
/musr/command rootOutput det_VrtxVolID off
/musr/command rootOutput det_VrtxProcID off
/musr/command rootOutput det_VrtxTrackID off
/musr/command rootOutput det_VrtxParticleID off
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeSecond off
/musr/command rootOutput odet_timeThird off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeMean off
/musr/command rootOutput odet_timeLast off
/musr/command rootOutput odet_timeCFD off
/musr/command rootOutput odet_amplCFD off
/musr/command rootOutput nOptPhot off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visVRML.mac
#/control/execute visFromToni.mac
#/control/execute visDawn50001.mac
#/vis/open VRML2FILE
#/vis/open DAWNFILE
### (if too many tracks cause core dump => storeTrajectory 0)
#/vis/scene/create
#
#/tracking/storeTrajectory 1
#/vis/viewer/set/viewpointThetaPhi 90 0
##/vis/viewer/set/globalLineWidthScale 3
#/vis/viewer/zoom 30
###/vis/scene/add/trajectories
#/vis/drawVolume
#/vis/viewer/flush
####/hits/verbose 2
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
#/gun/vertexsigma 20 20 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 128 -999999 999999 mm
#/gun/momentum 27.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
#/gun/momentumsmearing 0.81 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
/gun/turtlefilename FOR077_reggiani_Feb2010.dat
#/gun/turtleZ0position -900 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
#/gun/tiltsigma 0.2 0.2 0 deg
#/gun/pitch 0.5 deg
# Spin rotated by 50 degrees upwards:
/gun/muonPolarizVector 0.766043969 0. -0.642788174
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
#/gun/decaytimelimits 400 410 2197.03 ns
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 3000000
#/run/beamOn 10
/run/beamOn 10000000

341
run/50131.mac Normal file
View File

@ -0,0 +1,341 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 0 0 45
/musr/command rotation matrix2 0 180 0
/musr/command rotation matrix3 0 270 0
#/musr/command rotation matrix3 90 90
/musr/command arrayDef zPlaneGPSveto 12 0. 5. 91. 105. 115. 125. 135. 145. 155. 168.49 168.5 172.5
/musr/command arrayDef rInnerGPSveto 12 10. 11. 31. 31. 28.75 24.8. 19.7 14 10.3 9. 0. 0.
/musr/command arrayDef rOuterGPSveto 12 12.5 15.11 35.11 35. 32.85 29.1 24.2 18.6 14.6 13.2 13.0 13.0
# Rotation "fieldRot 0 0.57295 0" corresponds to 10 mrad (2.5mm/25cm) tilt of the magnetic field
#/musr/command rotation fieldRot 0 0.57295 0
#/musr/command rotation fieldRot 0 1.1459 0
# WORLD VOLUME
/musr/command construct box World 300 300 2000 G4_Galactic 0 0 0 no_logical_volume norot dead -1
# ---- BEAMPIPE vers. 1 (ALC-like)
#/musr/command construct tubs beampipe 0 34 810 0 360 G4_Galactic 0 0 -840 log_World norot dead 232
#/musr/command construct tubs beampipeshell 34 36 810 0 360 G4_Al 0 0 -840 log_World norot dead 231
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 16.5 34 25 0 360 G4_Pb 0 0 540 log_beampipe norot dead 241
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 34. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 2 (SEGMENTED and BROAD, small section 30cm long)
#/musr/command construct tubs beampipeA 0 128 500 0 360 G4_Galactic 0 0 -800 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 500 0 360 G4_Al 0 0 -800 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 135 0 360 G4_Galactic 0 0 -165 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 135 0 360 G4_Al 0 0 -165 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs pbCollimator_1 34 128 25 0 360 G4_Pb 0 0 475 log_beampipeA norot dead 241
#/musr/command construct tubs pbCollimator_2 16.5 34 25 0 360 G4_Pb 0 0 65 log_beampipeB norot dead 242
# ACTIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs activeCollim 2.5 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 3 (SEGMENTED and BROAD, small section 50cm long)
#/musr/command construct tubs beampipeA 0 128 400 0 360 G4_Galactic 0 0 -900 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 400 0 360 G4_Al 0 0 -900 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 230 0 360 G4_Galactic 0 0 -270 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 230 0 360 G4_Al 0 0 -270 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 20 128 15 0 360 G4_Pb 0 0 385 log_beampipeA norot dead 241
#/musr/command construct cons pbCollimator_2 122 128 34 128 70 0 360 G4_Pb 0 0 300 log_beampipeA norot dead 242
#/musr/command construct tubs pbCollimator_11 2.5 34 15 0 360 G4_Pb 0 0 215 log_beampipeB norot dead 245
#/musr/command construct cons pbCollimator_12 32 34 2.5 34 25 0 360 G4_Pb 0 0 175 log_beampipeB norot dead 246
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 5. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 150
#/musr/command construct tubs activeCollimB 5. 10. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 151
#/musr/command construct tubs activeCollimC 10. 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 152
#
# OLD ALC BEAMPIPE
# OLD ALC DESIGN: MSR 20.05.082 Tube
/musr/command construct tubs tube 35.45 38.05 315 0 360 Steel 0 0 -613 log_World norot dead 8201
/musr/command construct tubs tube_a 0 35.45 315 0 360 G4_Galactic 0 0 -613 log_World norot dead 8251
# SECOND COLLIMATOR
/musr/command construct tubs shield_collimatorA2 10 35 25 0 360 G4_Cu 0 0 290 log_tube_a norot dead 341
# OLD ALC DESIGN: MSR 20.05.083 Flange
/musr/command construct tubs flange 38.05 90 6 0 360 Steel 0 0 -934 log_World norot dead 8301
# ----
# POSITRON COUNTERS
/musr/command construct GPSforward forwardCounter 25 25 2.5 12.51 15.12 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward backwardCounter 30 30 2.5 8.80 7.64 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -30.0 log_World norot musr/ScintSD 2
/musr/command construct box up 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE 57.5 0 0.log_World norot musr/ScintSD 3
/musr/command construct box down 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE -57.5 0 0.log_World norot musr/ScintSD 4
/musr/command construct box rite1 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 -38 log_World norot musr/ScintSD 5
/musr/command construct box rite2 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 38 log_World norot musr/ScintSD 6
/musr/command construct box rite3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 -51.25 0 log_World norot musr/ScintSD 7
/musr/command construct box left1 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 -48.3 log_World norot musr/ScintSD 8
/musr/command construct box left2 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 48.3 log_World norot musr/ScintSD 9
/musr/command construct box left3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 53.75 0 log_World norot musr/ScintSD 10
# FORWARD VETO
/musr/command construct polyconeA forwardVeto 0 360 12 zPlaneGPSveto rInnerGPSveto rOuterGPSveto G4_PLASTIC_SC_VINYLTOLUENE 0 0 24 log_World norot musr/ScintSD 51
# BACKWARD VETO
/musr/command construct GPSbackwardVeto backwardVeto 35.13 7.63 35.13 7.63 60 31 3.5 31 3.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -87.5 log_World norot musr/ScintSD 52
#---------------------
# OUTER SHIELD
/musr/command construct tubeWithWindows tube1 20 21 33.5 14.5 33.5 14.5 33.5 G4_Cu 0 0 0 log_World matrix3 dead 401
/musr/command construct tubs bottom1 0 20 0.5 0 360 G4_Cu 0 -33.0 0 log_World matrix3 dead 402
/musr/command construct tubs window1 21 21.010 15 0 360 G4_Cu 0 0 0 log_World matrix3 dead 403
/musr/command construct tubs bottom1a 0 20 0.5 0 360 G4_Cu 0 33.0 0 log_World matrix3 dead 404
# INNER SHIELD
/musr/command construct tubeWithWindows tube2 11.5 12.5 21.5 7.5 21.5 10. 21.5 G4_Cu 0 0 0 log_World matrix3 dead 411
/musr/command construct tubs bottom2 0 11.5 3.0 0 360 G4_Cu 0 -21.5 0 log_World matrix3 dead 412
/musr/command construct tubs window2 12.5 12.510 8 0 360 G4_Cu 0 0 0 log_World matrix3 dead 413
/musr/command construct tubs bottom2a 0 11.5 3.0 0 360 G4_Cu 0 21.5 0 log_World matrix3 dead 414
# HEATER
/musr/command construct tubs heater 16 16.3 7.5 0 360 G4_Cu 0 -18.5 0 log_World matrix3 dead 415
#
#---------------------
# MAGNET
#/musr/command construct tubs magnet 44.4 100 150 0 360 G4_He 0 0 0 log_World norot dead 221
# MAGNET WALL
#/musr/command construct tubs mag_wall 44.4 46.5 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 222
# MAGNET VOLUME TO DELETE TRACKS
#/musr/command construct tubs sh0 46.5 100 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 223
# SHIELD TO DELETE TRACKS
#/musr/command construct tubs shield1 36 100 5 0 360 G4_AIR 0 0 -105 log_World norot dead -2
#/musr/command construct tubs sh2 0 100 5 0 360 G4_AIR 0 0 205 log_World norot dead -3
# NOSE EXTENSION
#/musr/command construct tubs noseExtensionA 15 25 2.0 0 360 G4_Al 0 0 -47 log_World norot dead 300
#/musr/command construct tubs noseExtensionB 25 28 126.5 0 360 G4_Al 0 0 -171.5 log_World norot dead 302
#/musr/command construct tubs noseExtensionC 28 41 4.0 0 360 G4_Al 0 0 -294 log_World norot dead 304
#
# COLLIMATOR
#/musr/command construct tubs collimatorA 2.5 25 15 0 360 G4_Pb 0 0 -64 log_World norot dead 310
#/musr/command construct tubs collimatorB 3.5 25 10 0 360 G4_Pb 0 0 -89 log_World norot dead 312
#/musr/command construct tubs collimatorC 5.0 25 10 0 360 G4_Pb 0 0 -109 log_World norot dead 314
#
# M COUNTERS AND M COUNTER HOLDER
#/musr/command construct tubs M0_holder1 5.5 7.5 3.25 0 360 G4_Al 0 0 -36.31 log_World norot dead 330
#/musr/command construct tubs M0_holder2 3.5 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -37.06 log_World norot dead 333
#/musr/command construct tubs M0_holder3 4.0 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -35.56 log_World norot dead 335
/musr/command construct tubs M0 0 7.5 0.15 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -25.0 log_World norot musr/ScintSD 102
#/musr/command construct tubs M0_electronics 4.0 9 0.85 0 360 G4_POLYCARBONATE 0 0 -40.41 log_World norot dead 337
#
# TARGET SPACE
#/musr/command construct tubs targetspace 0 5 50 0 360 G4_He 0 0 -20 log_cryostat
# SCINTILLATOR BEFORE TARGET
#/musr/command construct tubs coulombM1 0 3.0 0.1 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2 log_World norot musr/ScintSD 101
# TARGET
/musr/command construct tubs target 0 4.0 0.215 0 360 G4_Ag 0 0 0 log_World norot dead 201
/musr/command construct tubs targetFieldVol 0 0.5 0.015 0 360 G4_Ag 0 0 0 log_target norot dead 202
#/musr/command construct tubs vetoTarget 0 5.0 1. 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 1.220 log_World norot musr/ScintSD 160
#/musr/command construct tubs vetoCyl 5 6.0 3.6 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -1.39 log_World norot musr/ScintSD 161
#/musr/command construct TubeWithTubeHole vetoTarget 7 9 5.1 0 360 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct TubeWithHoleAndTubeHole vetoTarget 1.5 6 5.1 0 360 4 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct tubs sampleHolder 0 4 31.39 0 360 G4_Ag 0 0 31.61 log_World norot dead 165
#/musr/command construct tubs sampleHolder2 0 7.5 108.5 0 360 G4_Ag 0 0 171.5 log_World norot dead 166
#
#/musr/command construct TubeWithHolePlusTubeHole PlexyCyl1 4 12.5 4.5 0 360 7 12.5 5 G4_PLEXIGLASS 0 0 -10. log_World matrix2 dead 45
#/musr/command construct tubs PlexyCyl2 7.5 12.5 2 0 360 G4_PLEXIGLASS 0 0 12.5 log_World norot dead 46
#/musr/command construct TubeWithHolePlusTubeHole vetoCyl 2.5 12. 5. 0 360 4 12 9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -10.5 log_World matrix2 musr/ScintSD 51
#/musr/command construct tubs vetoCylA 8.0 12.0 10 91 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 61
#/musr/command construct tubs vetoCylB 8.0 12.0 10 181 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 62
#/musr/command construct tubs vetoCylC 8.0 12.0 10 271 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 63
#/musr/command construct tubs vetoCylD 8.0 12.0 10 1 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 64
#
#---------------------------------------------------------
#/musr/command region define goulombRegion log_target
#/musr/command region define goulombRegion log_M0
#/musr/command region setProductionCut goulombRegion 0.01 0.01 0.01
#---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
#---------------------------------------------------------
#
/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
/musr/command visattributes G4_PLEXIGLASS yellow
/musr/command visattributes G4_Pb blue
#/musr/command visattributes log_beampipe invisible
#/musr/command visattributes log_beampipeAshell invisible
#/musr/command visattributes log_beampipeBshell invisible
/musr/command visattributes log_World invisible
/musr/command visattributes log_target yellow
#/musr/command visattributes log_magnet invisible
#/musr/command visattributes log_mag_wall invisible
#/musr/command visattributes log_magnet yellow
#/musr/command visattributes log_mag_wall yellow
#/musr/command visattributes log_sh0 invisible
#/musr/command visattributes log_sh2 invisible
#/musr/command visattributes log_vetoTarget green
#/musr/command visattributes log_vetoCyl green
###################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPS_6kG_yAx_xRad_cgs.table log_targetFieldVol 0.005
/musr/command globalfield setparameter SetLargestAcceptableStep 2
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 -10
/musr/command globalfield printFieldValueAtPoint 0 10 10
/musr/command globalfield printFieldValueAtPoint 10 0 10
/musr/command globalfield printFieldValueAtPoint -10 0 10
/musr/command globalfield printFieldValueAtPoint 10 10 10
/musr/command globalfield printFieldValueAtPoint 40 40 100
/musr/command globalfield printFieldValueAtPoint 40 40 900
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
#/musr/command process addModel mu+ G4MuMultipleScattering G4UrbanMscModel92 0
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
#
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
#/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 0.1
#
/musr/run/howOftenToPrintEvent 100000
/musr/command maximumRunTimeAllowed 86000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_edep_el off
/musr/command rootOutput det_edep_pos off
/musr/command rootOutput det_edep_gam off
/musr/command rootOutput det_edep_mup off
/musr/command rootOutput det_nsteps off
/musr/command rootOutput det_length off
/musr/command rootOutput det_time_end off
/musr/command rootOutput det_kine off
/musr/command rootOutput det_VrtxKine off
/musr/command rootOutput det_VrtxX off
/musr/command rootOutput det_VrtxY off
/musr/command rootOutput det_VrtxZ off
/musr/command rootOutput det_VrtxVolID off
/musr/command rootOutput det_VrtxProcID off
/musr/command rootOutput det_VrtxTrackID off
/musr/command rootOutput det_VrtxParticleID off
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeSecond off
/musr/command rootOutput odet_timeThird off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeMean off
/musr/command rootOutput odet_timeLast off
/musr/command rootOutput odet_timeCFD off
/musr/command rootOutput odet_amplCFD off
/musr/command rootOutput nOptPhot off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visVRML.mac
#/control/execute visFromToni.mac
#/control/execute visDawn50001.mac
#/vis/open VRML2FILE
#/vis/open DAWNFILE
### (if too many tracks cause core dump => storeTrajectory 0)
#/vis/scene/create
#
#/tracking/storeTrajectory 1
#/vis/viewer/set/viewpointThetaPhi 90 0
##/vis/viewer/set/globalLineWidthScale 3
#/vis/viewer/zoom 30
###/vis/scene/add/trajectories
#/vis/drawVolume
#/vis/viewer/flush
####/hits/verbose 2
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
#/gun/vertexsigma 20 20 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 128 -999999 999999 mm
#/gun/momentum 27.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
#/gun/momentumsmearing 0.81 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
/gun/turtlefilename FOR077_reggiani_Feb2010.dat
#/gun/turtleZ0position -900 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
#/gun/tiltsigma 0.2 0.2 0 deg
#/gun/pitch 0.5 deg
# Spin rotated by 50 degrees upwards:
/gun/muonPolarizVector 0.766043969 0. -0.642788174
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
#/gun/decaytimelimits 400 410 2197.03 ns
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 3000000
#/run/beamOn 10
/run/beamOn 10000000

349
run/50161.mac Normal file
View File

@ -0,0 +1,349 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 0 0 45
/musr/command rotation matrix2 0 180 0
/musr/command rotation matrix3 0 270 0
#/musr/command rotation matrix3 90 90
/musr/command arrayDef zPlaneGPSveto 12 0. 5. 91. 105. 115. 125. 135. 145. 155. 168.49 168.5 172.5
/musr/command arrayDef rInnerGPSveto 12 10. 11. 31. 31. 28.75 24.8. 19.7 14 10.3 9. 0. 0.
/musr/command arrayDef rOuterGPSveto 12 12.5 15.11 35.11 35. 32.85 29.1 24.2 18.6 14.6 13.2 13.0 13.0
# Rotation "fieldRot 0 0.57295 0" corresponds to 10 mrad (2.5mm/25cm) tilt of the magnetic field
#/musr/command rotation fieldRot 0 0.57295 0
#/musr/command rotation fieldRot 0 1.1459 0
# WORLD VOLUME
/musr/command construct box World 300 300 2000 G4_Galactic 0 0 0 no_logical_volume norot dead -1
# ---- BEAMPIPE vers. 1 (ALC-like)
#/musr/command construct tubs beampipe 0 34 810 0 360 G4_Galactic 0 0 -840 log_World norot dead 232
#/musr/command construct tubs beampipeshell 34 36 810 0 360 G4_Al 0 0 -840 log_World norot dead 231
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 16.5 34 25 0 360 G4_Pb 0 0 540 log_beampipe norot dead 241
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 34. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 2 (SEGMENTED and BROAD, small section 30cm long)
#/musr/command construct tubs beampipeA 0 128 500 0 360 G4_Galactic 0 0 -800 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 500 0 360 G4_Al 0 0 -800 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 135 0 360 G4_Galactic 0 0 -165 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 135 0 360 G4_Al 0 0 -165 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs pbCollimator_1 34 128 25 0 360 G4_Pb 0 0 475 log_beampipeA norot dead 241
#/musr/command construct tubs pbCollimator_2 16.5 34 25 0 360 G4_Pb 0 0 65 log_beampipeB norot dead 242
# ACTIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs activeCollim 2.5 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 3 (SEGMENTED and BROAD, small section 50cm long)
#/musr/command construct tubs beampipeA 0 128 400 0 360 G4_Galactic 0 0 -900 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 400 0 360 G4_Al 0 0 -900 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 230 0 360 G4_Galactic 0 0 -270 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 230 0 360 G4_Al 0 0 -270 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 20 128 15 0 360 G4_Pb 0 0 385 log_beampipeA norot dead 241
#/musr/command construct cons pbCollimator_2 122 128 34 128 70 0 360 G4_Pb 0 0 300 log_beampipeA norot dead 242
#/musr/command construct tubs pbCollimator_11 2.5 34 15 0 360 G4_Pb 0 0 215 log_beampipeB norot dead 245
#/musr/command construct cons pbCollimator_12 32 34 2.5 34 25 0 360 G4_Pb 0 0 175 log_beampipeB norot dead 246
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 5. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 150
#/musr/command construct tubs activeCollimB 5. 10. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 151
#/musr/command construct tubs activeCollimC 10. 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 152
#
# OLD ALC BEAMPIPE
# OLD ALC DESIGN: MSR 20.05.082 Tube
/musr/command construct tubs tube 35.45 38.05 315 0 360 Steel 0 0 -613 log_World norot dead 8201
/musr/command construct tubs tube_a 0 35.45 315 0 360 G4_Galactic 0 0 -613 log_World norot dead 8251
# SECOND COLLIMATOR
/musr/command construct tubs shield_collimatorA2 10 35 25 0 360 G4_Cu 0 0 290 log_tube_a norot dead 341
# OLD ALC DESIGN: MSR 20.05.083 Flange
/musr/command construct tubs flange 38.05 90 6 0 360 Steel 0 0 -934 log_World norot dead 8301
# ----
# POSITRON COUNTERS
#/musr/command construct GPSforward forwardCounter 25 25 2.5 12.51 15.12 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward forwardCounter 25 25 2.5 9.1 11.9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward backwardCounter 30 30 2.5 8.80 7.64 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -30.0 log_World norot musr/ScintSD 2
/musr/command construct box up 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE 57.5 0 0.log_World norot musr/ScintSD 3
/musr/command construct box down 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE -57.5 0 0.log_World norot musr/ScintSD 4
/musr/command construct box rite1 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 -38 log_World norot musr/ScintSD 5
/musr/command construct box rite2 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 38 log_World norot musr/ScintSD 6
/musr/command construct box rite3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 -51.25 0 log_World norot musr/ScintSD 7
/musr/command construct box left1 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 -48.3 log_World norot musr/ScintSD 8
/musr/command construct box left2 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 48.3 log_World norot musr/ScintSD 9
/musr/command construct box left3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 53.75 0 log_World norot musr/ScintSD 10
# FORWARD VETO
#/musr/command construct polyconeA forwardVeto 0 360 12 zPlaneGPSveto rInnerGPSveto rOuterGPSveto G4_PLASTIC_SC_VINYLTOLUENE 0 0 24 log_World norot musr/ScintSD 51
/musr/command construct GPSbackwardVeto forwVetoA 9. 11.8 9. 11.8 2.5 6.25 6.7 6.25 6.7 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 51
/musr/command construct GPSbackwardVeto forwVetoB 11.8 23.5 11.8 23.5 37.5 6.7 18.4 6.7 18.4 G4_PLASTIC_SC_VINYLTOLUENE 0 0 66.5 log_World norot musr/ScintSD 51
/musr/command construct GPSbackwardVeto forwVetoC 23.5 16.1 23.5 16.1 17.5 18.4 11. 18.4 11. G4_PLASTIC_SC_VINYLTOLUENE 0 0 121.5 log_World norot musr/ScintSD 51
/musr/command construct box forwVetoD 16. 16. 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 141.5 log_World norot musr/ScintSD 51
#/musr/command construct box forwVetoD 23.5 23.5 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 106.5 log_World norot musr/ScintSD 51
#/musr/command construct GPSbackwardVeto forwVetoB 11.8 18.9 11.8 18.9 15 6.7 13.8 6.7 13.8 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44.0 log_World norot musr/ScintSD 51
#/musr/command construct box forwVetoD 18.9 18.9 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 61.5 log_World norot musr/ScintSD 51
# BACKWARD VETO
/musr/command construct GPSbackwardVeto backwardVeto 35.13 7.63 35.13 7.63 60 31 3.5 31 3.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -87.5 log_World norot musr/ScintSD 52
#---------------------
# OUTER SHIELD
/musr/command construct tubeWithWindows tube1 20 21 33.5 14.5 33.5 14.5 33.5 G4_Cu 0 0 0 log_World matrix3 dead 401
/musr/command construct tubs bottom1 0 20 0.5 0 360 G4_Cu 0 -33.0 0 log_World matrix3 dead 402
/musr/command construct tubs window1 21 21.010 15 0 360 G4_Cu 0 0 0 log_World matrix3 dead 403
/musr/command construct tubs bottom1a 0 20 0.5 0 360 G4_Cu 0 33.0 0 log_World matrix3 dead 404
# INNER SHIELD
/musr/command construct tubeWithWindows tube2 11.5 12.5 21.5 7.5 21.5 10. 21.5 G4_Cu 0 0 0 log_World matrix3 dead 411
/musr/command construct tubs bottom2 0 11.5 3.0 0 360 G4_Cu 0 -21.5 0 log_World matrix3 dead 412
/musr/command construct tubs window2 12.5 12.510 8 0 360 G4_Cu 0 0 0 log_World matrix3 dead 413
/musr/command construct tubs bottom2a 0 11.5 3.0 0 360 G4_Cu 0 21.5 0 log_World matrix3 dead 414
# HEATER
/musr/command construct tubs heater 16 16.3 7.5 0 360 G4_Cu 0 -18.5 0 log_World matrix3 dead 415
#
#---------------------
# MAGNET
#/musr/command construct tubs magnet 44.4 100 150 0 360 G4_He 0 0 0 log_World norot dead 221
# MAGNET WALL
#/musr/command construct tubs mag_wall 44.4 46.5 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 222
# MAGNET VOLUME TO DELETE TRACKS
#/musr/command construct tubs sh0 46.5 100 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 223
# SHIELD TO DELETE TRACKS
#/musr/command construct tubs shield1 36 100 5 0 360 G4_AIR 0 0 -105 log_World norot dead -2
#/musr/command construct tubs sh2 0 100 5 0 360 G4_AIR 0 0 205 log_World norot dead -3
# NOSE EXTENSION
#/musr/command construct tubs noseExtensionA 15 25 2.0 0 360 G4_Al 0 0 -47 log_World norot dead 300
#/musr/command construct tubs noseExtensionB 25 28 126.5 0 360 G4_Al 0 0 -171.5 log_World norot dead 302
#/musr/command construct tubs noseExtensionC 28 41 4.0 0 360 G4_Al 0 0 -294 log_World norot dead 304
#
# COLLIMATOR
#/musr/command construct tubs collimatorA 2.5 25 15 0 360 G4_Pb 0 0 -64 log_World norot dead 310
#/musr/command construct tubs collimatorB 3.5 25 10 0 360 G4_Pb 0 0 -89 log_World norot dead 312
#/musr/command construct tubs collimatorC 5.0 25 10 0 360 G4_Pb 0 0 -109 log_World norot dead 314
#
# M COUNTERS AND M COUNTER HOLDER
#/musr/command construct tubs M0_holder1 5.5 7.5 3.25 0 360 G4_Al 0 0 -36.31 log_World norot dead 330
#/musr/command construct tubs M0_holder2 3.5 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -37.06 log_World norot dead 333
#/musr/command construct tubs M0_holder3 4.0 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -35.56 log_World norot dead 335
/musr/command construct tubs M0 0 7.5 0.15 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -25.0 log_World norot musr/ScintSD 102
#/musr/command construct tubs M0_electronics 4.0 9 0.85 0 360 G4_POLYCARBONATE 0 0 -40.41 log_World norot dead 337
#
# TARGET SPACE
#/musr/command construct tubs targetspace 0 5 50 0 360 G4_He 0 0 -20 log_cryostat
# SCINTILLATOR BEFORE TARGET
#/musr/command construct tubs coulombM1 0 3.0 0.1 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2 log_World norot musr/ScintSD 101
# TARGET
/musr/command construct tubs target 0 4.0 0.215 0 360 G4_Ag 0 0 0 log_World norot dead 201
/musr/command construct tubs targetFieldVol 0 0.5 0.015 0 360 G4_Ag 0 0 0 log_target norot dead 202
#/musr/command construct tubs vetoTarget 0 5.0 1. 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 1.220 log_World norot musr/ScintSD 160
#/musr/command construct tubs vetoCyl 5 6.0 3.6 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -1.39 log_World norot musr/ScintSD 161
#/musr/command construct TubeWithTubeHole vetoTarget 7 9 5.1 0 360 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct TubeWithHoleAndTubeHole vetoTarget 1.5 6 5.1 0 360 4 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct tubs sampleHolder 0 4 31.39 0 360 G4_Ag 0 0 31.61 log_World norot dead 165
#/musr/command construct tubs sampleHolder2 0 7.5 108.5 0 360 G4_Ag 0 0 171.5 log_World norot dead 166
#
#/musr/command construct TubeWithHolePlusTubeHole PlexyCyl1 4 12.5 4.5 0 360 7 12.5 5 G4_PLEXIGLASS 0 0 -10. log_World matrix2 dead 45
#/musr/command construct tubs PlexyCyl2 7.5 12.5 2 0 360 G4_PLEXIGLASS 0 0 12.5 log_World norot dead 46
#/musr/command construct TubeWithHolePlusTubeHole vetoCyl 2.5 12. 5. 0 360 4 12 9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -10.5 log_World matrix2 musr/ScintSD 51
#/musr/command construct tubs vetoCylA 8.0 12.0 10 91 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 61
#/musr/command construct tubs vetoCylB 8.0 12.0 10 181 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 62
#/musr/command construct tubs vetoCylC 8.0 12.0 10 271 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 63
#/musr/command construct tubs vetoCylD 8.0 12.0 10 1 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 64
#
#---------------------------------------------------------
#/musr/command region define goulombRegion log_target
#/musr/command region define goulombRegion log_M0
#/musr/command region setProductionCut goulombRegion 0.01 0.01 0.01
#---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
#---------------------------------------------------------
#
/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
/musr/command visattributes G4_PLEXIGLASS yellow
/musr/command visattributes G4_Pb blue
#/musr/command visattributes log_beampipe invisible
#/musr/command visattributes log_beampipeAshell invisible
#/musr/command visattributes log_beampipeBshell invisible
/musr/command visattributes log_World invisible
/musr/command visattributes log_target yellow
#/musr/command visattributes log_magnet invisible
#/musr/command visattributes log_mag_wall invisible
#/musr/command visattributes log_magnet yellow
#/musr/command visattributes log_mag_wall yellow
#/musr/command visattributes log_sh0 invisible
#/musr/command visattributes log_sh2 invisible
#/musr/command visattributes log_vetoTarget green
#/musr/command visattributes log_vetoCyl green
###################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPS_6kG_yAx_xRad_cgs.table log_targetFieldVol 0.005
/musr/command globalfield setparameter SetLargestAcceptableStep 2
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 -10
/musr/command globalfield printFieldValueAtPoint 0 10 10
/musr/command globalfield printFieldValueAtPoint 10 0 10
/musr/command globalfield printFieldValueAtPoint -10 0 10
/musr/command globalfield printFieldValueAtPoint 10 10 10
/musr/command globalfield printFieldValueAtPoint 40 40 100
/musr/command globalfield printFieldValueAtPoint 40 40 900
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
#/musr/command process addModel mu+ G4MuMultipleScattering G4UrbanMscModel92 0
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
#
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
#/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 0.1
#
/musr/run/howOftenToPrintEvent 100000
/musr/command maximumRunTimeAllowed 86000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_edep_el off
/musr/command rootOutput det_edep_pos off
/musr/command rootOutput det_edep_gam off
/musr/command rootOutput det_edep_mup off
/musr/command rootOutput det_nsteps off
/musr/command rootOutput det_length off
/musr/command rootOutput det_time_end off
/musr/command rootOutput det_kine off
/musr/command rootOutput det_VrtxKine off
/musr/command rootOutput det_VrtxX off
/musr/command rootOutput det_VrtxY off
/musr/command rootOutput det_VrtxZ off
/musr/command rootOutput det_VrtxVolID off
/musr/command rootOutput det_VrtxProcID off
/musr/command rootOutput det_VrtxTrackID off
/musr/command rootOutput det_VrtxParticleID off
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeSecond off
/musr/command rootOutput odet_timeThird off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeMean off
/musr/command rootOutput odet_timeLast off
/musr/command rootOutput odet_timeCFD off
/musr/command rootOutput odet_amplCFD off
/musr/command rootOutput nOptPhot off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visVRML.mac
#/control/execute visFromToni.mac
#/control/execute visDawn50001.mac
#/vis/open VRML2FILE
#/vis/open DAWNFILE
### (if too many tracks cause core dump => storeTrajectory 0)
#/vis/scene/create
#
#/tracking/storeTrajectory 1
#/vis/viewer/set/viewpointThetaPhi 90 0
##/vis/viewer/set/globalLineWidthScale 3
#/vis/viewer/zoom 30
###/vis/scene/add/trajectories
#/vis/drawVolume
#/vis/viewer/flush
####/hits/verbose 2
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
#/gun/vertexsigma 20 20 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 128 -999999 999999 mm
#/gun/momentum 27.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
#/gun/momentumsmearing 0.81 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
/gun/turtlefilename FOR077_reggiani_Feb2010.dat
#/gun/turtleZ0position -900 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
#/gun/tiltsigma 0.2 0.2 0 deg
#/gun/pitch 0.5 deg
# Spin rotated by 50 degrees upwards:
/gun/muonPolarizVector 0.766043969 0. -0.642788174
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
#/gun/decaytimelimits 400 410 2197.03 ns
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 3000000
#/run/beamOn 10
/run/beamOn 10000000

349
run/50171.mac Normal file
View File

@ -0,0 +1,349 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 0 0 45
/musr/command rotation matrix2 0 180 0
/musr/command rotation matrix3 0 270 0
#/musr/command rotation matrix3 90 90
/musr/command arrayDef zPlaneGPSveto 12 0. 5. 91. 105. 115. 125. 135. 145. 155. 168.49 168.5 172.5
/musr/command arrayDef rInnerGPSveto 12 10. 11. 31. 31. 28.75 24.8. 19.7 14 10.3 9. 0. 0.
/musr/command arrayDef rOuterGPSveto 12 12.5 15.11 35.11 35. 32.85 29.1 24.2 18.6 14.6 13.2 13.0 13.0
# Rotation "fieldRot 0 0.57295 0" corresponds to 10 mrad (2.5mm/25cm) tilt of the magnetic field
#/musr/command rotation fieldRot 0 0.57295 0
#/musr/command rotation fieldRot 0 1.1459 0
# WORLD VOLUME
/musr/command construct box World 300 300 2000 G4_Galactic 0 0 0 no_logical_volume norot dead -1
# ---- BEAMPIPE vers. 1 (ALC-like)
#/musr/command construct tubs beampipe 0 34 810 0 360 G4_Galactic 0 0 -840 log_World norot dead 232
#/musr/command construct tubs beampipeshell 34 36 810 0 360 G4_Al 0 0 -840 log_World norot dead 231
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 16.5 34 25 0 360 G4_Pb 0 0 540 log_beampipe norot dead 241
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 34. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 2 (SEGMENTED and BROAD, small section 30cm long)
#/musr/command construct tubs beampipeA 0 128 500 0 360 G4_Galactic 0 0 -800 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 500 0 360 G4_Al 0 0 -800 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 135 0 360 G4_Galactic 0 0 -165 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 135 0 360 G4_Al 0 0 -165 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs pbCollimator_1 34 128 25 0 360 G4_Pb 0 0 475 log_beampipeA norot dead 241
#/musr/command construct tubs pbCollimator_2 16.5 34 25 0 360 G4_Pb 0 0 65 log_beampipeB norot dead 242
# ACTIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs activeCollim 2.5 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 3 (SEGMENTED and BROAD, small section 50cm long)
#/musr/command construct tubs beampipeA 0 128 400 0 360 G4_Galactic 0 0 -900 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 400 0 360 G4_Al 0 0 -900 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 230 0 360 G4_Galactic 0 0 -270 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 230 0 360 G4_Al 0 0 -270 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 20 128 15 0 360 G4_Pb 0 0 385 log_beampipeA norot dead 241
#/musr/command construct cons pbCollimator_2 122 128 34 128 70 0 360 G4_Pb 0 0 300 log_beampipeA norot dead 242
#/musr/command construct tubs pbCollimator_11 2.5 34 15 0 360 G4_Pb 0 0 215 log_beampipeB norot dead 245
#/musr/command construct cons pbCollimator_12 32 34 2.5 34 25 0 360 G4_Pb 0 0 175 log_beampipeB norot dead 246
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 5. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 150
#/musr/command construct tubs activeCollimB 5. 10. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 151
#/musr/command construct tubs activeCollimC 10. 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 152
#
# OLD ALC BEAMPIPE
# OLD ALC DESIGN: MSR 20.05.082 Tube
/musr/command construct tubs tube 35.45 38.05 315 0 360 Steel 0 0 -613 log_World norot dead 8201
/musr/command construct tubs tube_a 0 35.45 315 0 360 G4_Galactic 0 0 -613 log_World norot dead 8251
# SECOND COLLIMATOR
/musr/command construct tubs shield_collimatorA2 10 35 25 0 360 G4_Cu 0 0 290 log_tube_a norot dead 341
# OLD ALC DESIGN: MSR 20.05.083 Flange
/musr/command construct tubs flange 38.05 90 6 0 360 Steel 0 0 -934 log_World norot dead 8301
# ----
# POSITRON COUNTERS
#/musr/command construct GPSforward forwardCounter 25 25 2.5 12.51 15.12 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward forwardCounter 25 25 2.5 9.1 11.9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward backwardCounter 30 30 2.5 8.80 7.64 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -30.0 log_World norot musr/ScintSD 2
/musr/command construct box up 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE 57.5 0 0.log_World norot musr/ScintSD 3
/musr/command construct box down 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE -57.5 0 0.log_World norot musr/ScintSD 4
/musr/command construct box rite1 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 -38 log_World norot musr/ScintSD 5
/musr/command construct box rite2 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 38 log_World norot musr/ScintSD 6
/musr/command construct box rite3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 -51.25 0 log_World norot musr/ScintSD 7
/musr/command construct box left1 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 -48.3 log_World norot musr/ScintSD 8
/musr/command construct box left2 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 48.3 log_World norot musr/ScintSD 9
/musr/command construct box left3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 53.75 0 log_World norot musr/ScintSD 10
# FORWARD VETO
#/musr/command construct polyconeA forwardVeto 0 360 12 zPlaneGPSveto rInnerGPSveto rOuterGPSveto G4_PLASTIC_SC_VINYLTOLUENE 0 0 24 log_World norot musr/ScintSD 51
/musr/command construct GPSbackwardVeto forwVetoA 9. 11.8 9. 11.8 2.5 6.25 6.7 6.25 6.7 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 51
/musr/command construct GPSbackwardVeto forwVetoB 11.8 23.5 11.8 23.5 37.5 6.7 18.4 6.7 18.4 G4_PLASTIC_SC_VINYLTOLUENE 0 0 66.5 log_World norot musr/ScintSD 51
#/musr/command construct GPSbackwardVeto forwVetoC 23.5 16.1 23.5 16.1 17.5 18.4 11. 18.4 11. G4_PLASTIC_SC_VINYLTOLUENE 0 0 121.5 log_World norot musr/ScintSD 51
#/musr/command construct box forwVetoD 16. 16. 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 141.5 log_World norot musr/ScintSD 51
/musr/command construct box forwVetoD 23.5 23.5 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 106.5 log_World norot musr/ScintSD 51
#/musr/command construct GPSbackwardVeto forwVetoB 11.8 18.9 11.8 18.9 15 6.7 13.8 6.7 13.8 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44.0 log_World norot musr/ScintSD 51
#/musr/command construct box forwVetoD 18.9 18.9 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 61.5 log_World norot musr/ScintSD 51
# BACKWARD VETO
/musr/command construct GPSbackwardVeto backwardVeto 35.13 7.63 35.13 7.63 60 31 3.5 31 3.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -87.5 log_World norot musr/ScintSD 52
#---------------------
# OUTER SHIELD
/musr/command construct tubeWithWindows tube1 20 21 33.5 14.5 33.5 14.5 33.5 G4_Cu 0 0 0 log_World matrix3 dead 401
/musr/command construct tubs bottom1 0 20 0.5 0 360 G4_Cu 0 -33.0 0 log_World matrix3 dead 402
/musr/command construct tubs window1 21 21.010 15 0 360 G4_Cu 0 0 0 log_World matrix3 dead 403
/musr/command construct tubs bottom1a 0 20 0.5 0 360 G4_Cu 0 33.0 0 log_World matrix3 dead 404
# INNER SHIELD
/musr/command construct tubeWithWindows tube2 11.5 12.5 21.5 7.5 21.5 10. 21.5 G4_Cu 0 0 0 log_World matrix3 dead 411
/musr/command construct tubs bottom2 0 11.5 3.0 0 360 G4_Cu 0 -21.5 0 log_World matrix3 dead 412
/musr/command construct tubs window2 12.5 12.510 8 0 360 G4_Cu 0 0 0 log_World matrix3 dead 413
/musr/command construct tubs bottom2a 0 11.5 3.0 0 360 G4_Cu 0 21.5 0 log_World matrix3 dead 414
# HEATER
/musr/command construct tubs heater 16 16.3 7.5 0 360 G4_Cu 0 -18.5 0 log_World matrix3 dead 415
#
#---------------------
# MAGNET
#/musr/command construct tubs magnet 44.4 100 150 0 360 G4_He 0 0 0 log_World norot dead 221
# MAGNET WALL
#/musr/command construct tubs mag_wall 44.4 46.5 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 222
# MAGNET VOLUME TO DELETE TRACKS
#/musr/command construct tubs sh0 46.5 100 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 223
# SHIELD TO DELETE TRACKS
#/musr/command construct tubs shield1 36 100 5 0 360 G4_AIR 0 0 -105 log_World norot dead -2
#/musr/command construct tubs sh2 0 100 5 0 360 G4_AIR 0 0 205 log_World norot dead -3
# NOSE EXTENSION
#/musr/command construct tubs noseExtensionA 15 25 2.0 0 360 G4_Al 0 0 -47 log_World norot dead 300
#/musr/command construct tubs noseExtensionB 25 28 126.5 0 360 G4_Al 0 0 -171.5 log_World norot dead 302
#/musr/command construct tubs noseExtensionC 28 41 4.0 0 360 G4_Al 0 0 -294 log_World norot dead 304
#
# COLLIMATOR
#/musr/command construct tubs collimatorA 2.5 25 15 0 360 G4_Pb 0 0 -64 log_World norot dead 310
#/musr/command construct tubs collimatorB 3.5 25 10 0 360 G4_Pb 0 0 -89 log_World norot dead 312
#/musr/command construct tubs collimatorC 5.0 25 10 0 360 G4_Pb 0 0 -109 log_World norot dead 314
#
# M COUNTERS AND M COUNTER HOLDER
#/musr/command construct tubs M0_holder1 5.5 7.5 3.25 0 360 G4_Al 0 0 -36.31 log_World norot dead 330
#/musr/command construct tubs M0_holder2 3.5 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -37.06 log_World norot dead 333
#/musr/command construct tubs M0_holder3 4.0 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -35.56 log_World norot dead 335
/musr/command construct tubs M0 0 7.5 0.15 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -25.0 log_World norot musr/ScintSD 102
#/musr/command construct tubs M0_electronics 4.0 9 0.85 0 360 G4_POLYCARBONATE 0 0 -40.41 log_World norot dead 337
#
# TARGET SPACE
#/musr/command construct tubs targetspace 0 5 50 0 360 G4_He 0 0 -20 log_cryostat
# SCINTILLATOR BEFORE TARGET
#/musr/command construct tubs coulombM1 0 3.0 0.1 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2 log_World norot musr/ScintSD 101
# TARGET
/musr/command construct tubs target 0 4.0 0.215 0 360 G4_Ag 0 0 0 log_World norot dead 201
/musr/command construct tubs targetFieldVol 0 0.5 0.015 0 360 G4_Ag 0 0 0 log_target norot dead 202
#/musr/command construct tubs vetoTarget 0 5.0 1. 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 1.220 log_World norot musr/ScintSD 160
#/musr/command construct tubs vetoCyl 5 6.0 3.6 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -1.39 log_World norot musr/ScintSD 161
#/musr/command construct TubeWithTubeHole vetoTarget 7 9 5.1 0 360 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct TubeWithHoleAndTubeHole vetoTarget 1.5 6 5.1 0 360 4 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct tubs sampleHolder 0 4 31.39 0 360 G4_Ag 0 0 31.61 log_World norot dead 165
#/musr/command construct tubs sampleHolder2 0 7.5 108.5 0 360 G4_Ag 0 0 171.5 log_World norot dead 166
#
#/musr/command construct TubeWithHolePlusTubeHole PlexyCyl1 4 12.5 4.5 0 360 7 12.5 5 G4_PLEXIGLASS 0 0 -10. log_World matrix2 dead 45
#/musr/command construct tubs PlexyCyl2 7.5 12.5 2 0 360 G4_PLEXIGLASS 0 0 12.5 log_World norot dead 46
#/musr/command construct TubeWithHolePlusTubeHole vetoCyl 2.5 12. 5. 0 360 4 12 9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -10.5 log_World matrix2 musr/ScintSD 51
#/musr/command construct tubs vetoCylA 8.0 12.0 10 91 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 61
#/musr/command construct tubs vetoCylB 8.0 12.0 10 181 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 62
#/musr/command construct tubs vetoCylC 8.0 12.0 10 271 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 63
#/musr/command construct tubs vetoCylD 8.0 12.0 10 1 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 64
#
#---------------------------------------------------------
#/musr/command region define goulombRegion log_target
#/musr/command region define goulombRegion log_M0
#/musr/command region setProductionCut goulombRegion 0.01 0.01 0.01
#---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
#---------------------------------------------------------
#
/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
/musr/command visattributes G4_PLEXIGLASS yellow
/musr/command visattributes G4_Pb blue
#/musr/command visattributes log_beampipe invisible
#/musr/command visattributes log_beampipeAshell invisible
#/musr/command visattributes log_beampipeBshell invisible
/musr/command visattributes log_World invisible
/musr/command visattributes log_target yellow
#/musr/command visattributes log_magnet invisible
#/musr/command visattributes log_mag_wall invisible
#/musr/command visattributes log_magnet yellow
#/musr/command visattributes log_mag_wall yellow
#/musr/command visattributes log_sh0 invisible
#/musr/command visattributes log_sh2 invisible
#/musr/command visattributes log_vetoTarget green
#/musr/command visattributes log_vetoCyl green
###################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPS_6kG_yAx_xRad_cgs.table log_targetFieldVol 0.005
/musr/command globalfield setparameter SetLargestAcceptableStep 2
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 -10
/musr/command globalfield printFieldValueAtPoint 0 10 10
/musr/command globalfield printFieldValueAtPoint 10 0 10
/musr/command globalfield printFieldValueAtPoint -10 0 10
/musr/command globalfield printFieldValueAtPoint 10 10 10
/musr/command globalfield printFieldValueAtPoint 40 40 100
/musr/command globalfield printFieldValueAtPoint 40 40 900
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
#/musr/command process addModel mu+ G4MuMultipleScattering G4UrbanMscModel92 0
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
#
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
#/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 0.1
#
/musr/run/howOftenToPrintEvent 100000
/musr/command maximumRunTimeAllowed 86000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_edep_el off
/musr/command rootOutput det_edep_pos off
/musr/command rootOutput det_edep_gam off
/musr/command rootOutput det_edep_mup off
/musr/command rootOutput det_nsteps off
/musr/command rootOutput det_length off
/musr/command rootOutput det_time_end off
/musr/command rootOutput det_kine off
/musr/command rootOutput det_VrtxKine off
/musr/command rootOutput det_VrtxX off
/musr/command rootOutput det_VrtxY off
/musr/command rootOutput det_VrtxZ off
/musr/command rootOutput det_VrtxVolID off
/musr/command rootOutput det_VrtxProcID off
/musr/command rootOutput det_VrtxTrackID off
/musr/command rootOutput det_VrtxParticleID off
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeSecond off
/musr/command rootOutput odet_timeThird off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeMean off
/musr/command rootOutput odet_timeLast off
/musr/command rootOutput odet_timeCFD off
/musr/command rootOutput odet_amplCFD off
/musr/command rootOutput nOptPhot off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visVRML.mac
#/control/execute visFromToni.mac
#/control/execute visDawn50001.mac
#/vis/open VRML2FILE
#/vis/open DAWNFILE
### (if too many tracks cause core dump => storeTrajectory 0)
#/vis/scene/create
#
#/tracking/storeTrajectory 1
#/vis/viewer/set/viewpointThetaPhi 90 0
##/vis/viewer/set/globalLineWidthScale 3
#/vis/viewer/zoom 30
###/vis/scene/add/trajectories
#/vis/drawVolume
#/vis/viewer/flush
####/hits/verbose 2
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
#/gun/vertexsigma 20 20 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 128 -999999 999999 mm
#/gun/momentum 27.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
#/gun/momentumsmearing 0.81 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
/gun/turtlefilename FOR077_reggiani_Feb2010.dat
#/gun/turtleZ0position -900 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
#/gun/tiltsigma 0.2 0.2 0 deg
#/gun/pitch 0.5 deg
# Spin rotated by 50 degrees upwards:
/gun/muonPolarizVector 0.766043969 0. -0.642788174
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
#/gun/decaytimelimits 400 410 2197.03 ns
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 3000000
#/run/beamOn 10
/run/beamOn 10000000

349
run/50181.mac Normal file
View File

@ -0,0 +1,349 @@
#-----------------------------------------------------------------------
# Macro file for muSR instrument simulations
# Unless specified otherwises, the default units are mm, ns, MeV, MeV/c.
# Lines starting with star "#" are comments.
###################################################################################
############################# G E O M E T R Y ###################################
###################################################################################
# ROTATION MATRIXES:
/musr/command rotation matrix1 0 0 45
/musr/command rotation matrix2 0 180 0
/musr/command rotation matrix3 0 270 0
#/musr/command rotation matrix3 90 90
/musr/command arrayDef zPlaneGPSveto 12 0. 5. 91. 105. 115. 125. 135. 145. 155. 168.49 168.5 172.5
/musr/command arrayDef rInnerGPSveto 12 10. 11. 31. 31. 28.75 24.8. 19.7 14 10.3 9. 0. 0.
/musr/command arrayDef rOuterGPSveto 12 12.5 15.11 35.11 35. 32.85 29.1 24.2 18.6 14.6 13.2 13.0 13.0
# Rotation "fieldRot 0 0.57295 0" corresponds to 10 mrad (2.5mm/25cm) tilt of the magnetic field
#/musr/command rotation fieldRot 0 0.57295 0
#/musr/command rotation fieldRot 0 1.1459 0
# WORLD VOLUME
/musr/command construct box World 300 300 2000 G4_Galactic 0 0 0 no_logical_volume norot dead -1
# ---- BEAMPIPE vers. 1 (ALC-like)
#/musr/command construct tubs beampipe 0 34 810 0 360 G4_Galactic 0 0 -840 log_World norot dead 232
#/musr/command construct tubs beampipeshell 34 36 810 0 360 G4_Al 0 0 -840 log_World norot dead 231
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 16.5 34 25 0 360 G4_Pb 0 0 540 log_beampipe norot dead 241
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 34. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 2 (SEGMENTED and BROAD, small section 30cm long)
#/musr/command construct tubs beampipeA 0 128 500 0 360 G4_Galactic 0 0 -800 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 500 0 360 G4_Al 0 0 -800 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 135 0 360 G4_Galactic 0 0 -165 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 135 0 360 G4_Al 0 0 -165 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs pbCollimator_1 34 128 25 0 360 G4_Pb 0 0 475 log_beampipeA norot dead 241
#/musr/command construct tubs pbCollimator_2 16.5 34 25 0 360 G4_Pb 0 0 65 log_beampipeB norot dead 242
# ACTIVE COLLIMATOR (for BEAMPIPE vers.2)
#/musr/command construct tubs activeCollim 2.5 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -27 log_World norot musr/ScintSD 150
#
# ---- BEAMPIPE vers. 3 (SEGMENTED and BROAD, small section 50cm long)
#/musr/command construct tubs beampipeA 0 128 400 0 360 G4_Galactic 0 0 -900 log_World norot dead 232
#/musr/command construct tubs beampipeAshell 128 130 400 0 360 G4_Al 0 0 -900 log_World norot dead 231
#/musr/command construct tubs beampipeB 0 34 230 0 360 G4_Galactic 0 0 -270 log_World norot dead 234
#/musr/command construct tubs beampipeBshell 34 36 230 0 360 G4_Al 0 0 -270 log_World norot dead 233
# PASSIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs pbCollimator_1 20 128 15 0 360 G4_Pb 0 0 385 log_beampipeA norot dead 241
#/musr/command construct cons pbCollimator_2 122 128 34 128 70 0 360 G4_Pb 0 0 300 log_beampipeA norot dead 242
#/musr/command construct tubs pbCollimator_11 2.5 34 15 0 360 G4_Pb 0 0 215 log_beampipeB norot dead 245
#/musr/command construct cons pbCollimator_12 32 34 2.5 34 25 0 360 G4_Pb 0 0 175 log_beampipeB norot dead 246
# ACTIVE COLLIMATOR (for BEAMPIPE vers.1)
#/musr/command construct tubs activeCollim 2.5 5. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 150
#/musr/command construct tubs activeCollimB 5. 10. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 151
#/musr/command construct tubs activeCollimC 10. 36. 2.5 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -31 log_World norot musr/ScintSD 152
#
# OLD ALC BEAMPIPE
# OLD ALC DESIGN: MSR 20.05.082 Tube
/musr/command construct tubs tube 35.45 38.05 315 0 360 Steel 0 0 -613 log_World norot dead 8201
/musr/command construct tubs tube_a 0 35.45 315 0 360 G4_Galactic 0 0 -613 log_World norot dead 8251
# SECOND COLLIMATOR
/musr/command construct tubs shield_collimatorA2 10 35 25 0 360 G4_Cu 0 0 290 log_tube_a norot dead 341
# OLD ALC DESIGN: MSR 20.05.083 Flange
/musr/command construct tubs flange 38.05 90 6 0 360 Steel 0 0 -934 log_World norot dead 8301
# ----
# POSITRON COUNTERS
#/musr/command construct GPSforward forwardCounter 25 25 2.5 12.51 15.12 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward forwardCounter 25 25 2.5 9.1 11.9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 1
/musr/command construct GPSbackward backwardCounter 30 30 2.5 8.80 7.64 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -30.0 log_World norot musr/ScintSD 2
/musr/command construct box up 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE 57.5 0 0.log_World norot musr/ScintSD 3
/musr/command construct box down 2.5 53 53 G4_PLASTIC_SC_VINYLTOLUENE -57.5 0 0.log_World norot musr/ScintSD 4
/musr/command construct box rite1 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 -38 log_World norot musr/ScintSD 5
/musr/command construct box rite2 53 2.5 14 G4_PLASTIC_SC_VINYLTOLUENE 0 -56.25 38 log_World norot musr/ScintSD 6
/musr/command construct box rite3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 -51.25 0 log_World norot musr/ScintSD 7
/musr/command construct box left1 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 -48.3 log_World norot musr/ScintSD 8
/musr/command construct box left2 53 2.5 19.3 G4_PLASTIC_SC_VINYLTOLUENE 0 72.5 48.3 log_World norot musr/ScintSD 9
/musr/command construct box left3 35 2.5 20 G4_PLASTIC_SC_VINYLTOLUENE 0 53.75 0 log_World norot musr/ScintSD 10
# FORWARD VETO
#/musr/command construct polyconeA forwardVeto 0 360 12 zPlaneGPSveto rInnerGPSveto rOuterGPSveto G4_PLASTIC_SC_VINYLTOLUENE 0 0 24 log_World norot musr/ScintSD 51
/musr/command construct GPSbackwardVeto forwVetoA 9. 11.8 9. 11.8 2.5 6.25 6.7 6.25 6.7 G4_PLASTIC_SC_VINYLTOLUENE 0 0 26.5 log_World norot musr/ScintSD 51
#/musr/command construct GPSbackwardVeto forwVetoB 11.8 23.5 11.8 23.5 37.5 6.7 18.4 6.7 18.4 G4_PLASTIC_SC_VINYLTOLUENE 0 0 66.5 log_World norot musr/ScintSD 51
#/musr/command construct GPSbackwardVeto forwVetoC 23.5 16.1 23.5 16.1 17.5 18.4 11. 18.4 11. G4_PLASTIC_SC_VINYLTOLUENE 0 0 121.5 log_World norot musr/ScintSD 51
#/musr/command construct box forwVetoD 16. 16. 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 141.5 log_World norot musr/ScintSD 51
#/musr/command construct box forwVetoD 23.5 23.5 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 106.5 log_World norot musr/ScintSD 51
/musr/command construct GPSbackwardVeto forwVetoB 11.8 18.9 11.8 18.9 15 6.7 13.8 6.7 13.8 G4_PLASTIC_SC_VINYLTOLUENE 0 0 44.0 log_World norot musr/ScintSD 51
/musr/command construct box forwVetoD 18.9 18.9 2.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 61.5 log_World norot musr/ScintSD 51
# BACKWARD VETO
/musr/command construct GPSbackwardVeto backwardVeto 35.13 7.63 35.13 7.63 60 31 3.5 31 3.5 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -87.5 log_World norot musr/ScintSD 52
#---------------------
# OUTER SHIELD
/musr/command construct tubeWithWindows tube1 20 21 33.5 14.5 33.5 14.5 33.5 G4_Cu 0 0 0 log_World matrix3 dead 401
/musr/command construct tubs bottom1 0 20 0.5 0 360 G4_Cu 0 -33.0 0 log_World matrix3 dead 402
/musr/command construct tubs window1 21 21.010 15 0 360 G4_Cu 0 0 0 log_World matrix3 dead 403
/musr/command construct tubs bottom1a 0 20 0.5 0 360 G4_Cu 0 33.0 0 log_World matrix3 dead 404
# INNER SHIELD
/musr/command construct tubeWithWindows tube2 11.5 12.5 21.5 7.5 21.5 10. 21.5 G4_Cu 0 0 0 log_World matrix3 dead 411
/musr/command construct tubs bottom2 0 11.5 3.0 0 360 G4_Cu 0 -21.5 0 log_World matrix3 dead 412
/musr/command construct tubs window2 12.5 12.510 8 0 360 G4_Cu 0 0 0 log_World matrix3 dead 413
/musr/command construct tubs bottom2a 0 11.5 3.0 0 360 G4_Cu 0 21.5 0 log_World matrix3 dead 414
# HEATER
/musr/command construct tubs heater 16 16.3 7.5 0 360 G4_Cu 0 -18.5 0 log_World matrix3 dead 415
#
#---------------------
# MAGNET
#/musr/command construct tubs magnet 44.4 100 150 0 360 G4_He 0 0 0 log_World norot dead 221
# MAGNET WALL
#/musr/command construct tubs mag_wall 44.4 46.5 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 222
# MAGNET VOLUME TO DELETE TRACKS
#/musr/command construct tubs sh0 46.5 100 150 0 360 G4_Cu 0 0 0 log_magnet norot dead 223
# SHIELD TO DELETE TRACKS
#/musr/command construct tubs shield1 36 100 5 0 360 G4_AIR 0 0 -105 log_World norot dead -2
#/musr/command construct tubs sh2 0 100 5 0 360 G4_AIR 0 0 205 log_World norot dead -3
# NOSE EXTENSION
#/musr/command construct tubs noseExtensionA 15 25 2.0 0 360 G4_Al 0 0 -47 log_World norot dead 300
#/musr/command construct tubs noseExtensionB 25 28 126.5 0 360 G4_Al 0 0 -171.5 log_World norot dead 302
#/musr/command construct tubs noseExtensionC 28 41 4.0 0 360 G4_Al 0 0 -294 log_World norot dead 304
#
# COLLIMATOR
#/musr/command construct tubs collimatorA 2.5 25 15 0 360 G4_Pb 0 0 -64 log_World norot dead 310
#/musr/command construct tubs collimatorB 3.5 25 10 0 360 G4_Pb 0 0 -89 log_World norot dead 312
#/musr/command construct tubs collimatorC 5.0 25 10 0 360 G4_Pb 0 0 -109 log_World norot dead 314
#
# M COUNTERS AND M COUNTER HOLDER
#/musr/command construct tubs M0_holder1 5.5 7.5 3.25 0 360 G4_Al 0 0 -36.31 log_World norot dead 330
#/musr/command construct tubs M0_holder2 3.5 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -37.06 log_World norot dead 333
#/musr/command construct tubs M0_holder3 4.0 5.5 0.75 0 360 G4_PLEXIGLASS 0 0 -35.56 log_World norot dead 335
/musr/command construct tubs M0 0 7.5 0.15 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -25.0 log_World norot musr/ScintSD 102
#/musr/command construct tubs M0_electronics 4.0 9 0.85 0 360 G4_POLYCARBONATE 0 0 -40.41 log_World norot dead 337
#
# TARGET SPACE
#/musr/command construct tubs targetspace 0 5 50 0 360 G4_He 0 0 -20 log_cryostat
# SCINTILLATOR BEFORE TARGET
#/musr/command construct tubs coulombM1 0 3.0 0.1 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2 log_World norot musr/ScintSD 101
# TARGET
/musr/command construct tubs target 0 4.0 0.215 0 360 G4_Ag 0 0 0 log_World norot dead 201
/musr/command construct tubs targetFieldVol 0 0.5 0.015 0 360 G4_Ag 0 0 0 log_target norot dead 202
#/musr/command construct tubs vetoTarget 0 5.0 1. 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 1.220 log_World norot musr/ScintSD 160
#/musr/command construct tubs vetoCyl 5 6.0 3.6 0 360 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -1.39 log_World norot musr/ScintSD 161
#/musr/command construct TubeWithTubeHole vetoTarget 7 9 5.1 0 360 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct TubeWithHoleAndTubeHole vetoTarget 1.5 6 5.1 0 360 4 8.2 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -2.89 log_World norot musr/ScintSD 160
#/musr/command construct tubs sampleHolder 0 4 31.39 0 360 G4_Ag 0 0 31.61 log_World norot dead 165
#/musr/command construct tubs sampleHolder2 0 7.5 108.5 0 360 G4_Ag 0 0 171.5 log_World norot dead 166
#
#/musr/command construct TubeWithHolePlusTubeHole PlexyCyl1 4 12.5 4.5 0 360 7 12.5 5 G4_PLEXIGLASS 0 0 -10. log_World matrix2 dead 45
#/musr/command construct tubs PlexyCyl2 7.5 12.5 2 0 360 G4_PLEXIGLASS 0 0 12.5 log_World norot dead 46
#/musr/command construct TubeWithHolePlusTubeHole vetoCyl 2.5 12. 5. 0 360 4 12 9 G4_PLASTIC_SC_VINYLTOLUENE 0 0 -10.5 log_World matrix2 musr/ScintSD 51
#/musr/command construct tubs vetoCylA 8.0 12.0 10 91 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 61
#/musr/command construct tubs vetoCylB 8.0 12.0 10 181 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 62
#/musr/command construct tubs vetoCylC 8.0 12.0 10 271 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 63
#/musr/command construct tubs vetoCylD 8.0 12.0 10 1 88 G4_PLASTIC_SC_VINYLTOLUENE 0 0 0 log_World norot musr/ScintSD 64
#
#---------------------------------------------------------
#/musr/command region define goulombRegion log_target
#/musr/command region define goulombRegion log_M0
#/musr/command region setProductionCut goulombRegion 0.01 0.01 0.01
#---------------------------------------------------------
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_1 6
#/musr/command logicalVolumeToBeReweighted mu log_pbCollimator_2 5
#/musr/command logicalVolumeToBeReweighted mu log_activeCollim 4
#---------------------------------------------------------
#
/musr/command visattributes G4_PLASTIC_SC_VINYLTOLUENE red
/musr/command visattributes G4_PLEXIGLASS yellow
/musr/command visattributes G4_Pb blue
#/musr/command visattributes log_beampipe invisible
#/musr/command visattributes log_beampipeAshell invisible
#/musr/command visattributes log_beampipeBshell invisible
/musr/command visattributes log_World invisible
/musr/command visattributes log_target yellow
#/musr/command visattributes log_magnet invisible
#/musr/command visattributes log_mag_wall invisible
#/musr/command visattributes log_magnet yellow
#/musr/command visattributes log_mag_wall yellow
#/musr/command visattributes log_sh0 invisible
#/musr/command visattributes log_sh2 invisible
#/musr/command visattributes log_vetoTarget green
#/musr/command visattributes log_vetoCyl green
###################################################################################
######################### M A G N E T I C F I E L D #########################
###################################################################################
# Set magnetic field (set field intensity in T and sigma in mm)
# syntax for magneticfield: fromfile filename fieldValue
# uniform fieldValue
# gaussian fieldValue sigma
#
/musr/command globalfield centralSolenoidField 0. 0. 0. fromfile 2D GPS_6kG_yAx_xRad_cgs.table log_targetFieldVol 0.005
/musr/command globalfield setparameter SetLargestAcceptableStep 2
/musr/command globalfield setparameter SetMinimumEpsilonStep 5e-5
/musr/command globalfield setparameter SetMaximumEpsilonStep 0.001
/musr/command globalfield setparameter SetDeltaOneStep 0.1
/musr/command globalfield setparameter SetDeltaIntersection 0.01
/musr/command globalfield printparameters
/musr/command globalfield printFieldValueAtPoint 0 0 10
/musr/command globalfield printFieldValueAtPoint 0 0 -10
/musr/command globalfield printFieldValueAtPoint 0 10 10
/musr/command globalfield printFieldValueAtPoint 10 0 10
/musr/command globalfield printFieldValueAtPoint -10 0 10
/musr/command globalfield printFieldValueAtPoint 10 10 10
/musr/command globalfield printFieldValueAtPoint 40 40 100
/musr/command globalfield printFieldValueAtPoint 40 40 900
###################################################################################
######################### P H Y S I C S P R O C E S S E S ##################
###################################################################################
# Geant 4.9.3
/musr/command process addDiscreteProcess gamma G4PhotoElectricEffect
/musr/command process addDiscreteProcess gamma G4ComptonScattering
/musr/command process addDiscreteProcess gamma G4GammaConversion
/musr/command process addDiscreteProcess gamma G4RayleighScattering
/musr/command process addProcess e- G4eMultipleScattering -1 1 1
/musr/command process addProcess e- G4eIonisation -1 2 2
/musr/command process addProcess e- G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eMultipleScattering -1 1 1
/musr/command process addProcess e+ G4eIonisation -1 2 2
/musr/command process addProcess e+ G4eBremsstrahlung -1 3 3
/musr/command process addProcess e+ G4eplusAnnihilation 0 -1 4
/musr/command process addProcess mu- G4MuMultipleScattering -1 1 1
/musr/command process addProcess mu- G4MuIonisation -1 2 2
/musr/command process addProcess mu- G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu- G4MuPairProduction -1 4 4
/musr/command process addProcess mu+ G4MuMultipleScattering -1 1 1
#/musr/command process addModel mu+ G4MuMultipleScattering G4UrbanMscModel92 0
/musr/command process addProcess mu+ G4MuIonisation -1 2 2
/musr/command process addProcess mu+ G4MuBremsstrahlung -1 3 3
/musr/command process addProcess mu+ G4MuPairProduction -1 4 4
#
#
#/home/install/geant4.9.3/source/processes/electromagnetic/utils/src/G4EnergyLossMessenger.cc
######## /process/msc/StepLimit Minimal | UseDistanceToBoundary | UseSafety
/process/msc/StepLimit UseSafety
#/process/msc/LateralDisplacement
/process/msc/RangeFactor 0.04
/process/msc/GeomFactor 2.5
#/process/msc/FactorForAngleLimit
/process/msc/Skin 3.0
#/process/msc/ThetaLimit 0.2 rad
#
###################################################################################
################## S O M E O T H E R P A R A M E T E R S ##################
###################################################################################
# Set the overall range cut
#/run/setCut 1 mm
# Set the range cut on a particular volumes (in mm)
/musr/command SetUserLimits log_target 0.01
#/musr/command SetUserLimits log_targetscint 0.01
/musr/command SetUserLimits log_M0 0.01
# Store all events into the ROOT tree or just the interesting ones ? (true is default)
#/musr/command storeOnlyEventsWithHits false
# Set the minimum time separation between two subsequent signals in the same detector (in ns)
/musr/command signalSeparationTime 0.1
#
/musr/run/howOftenToPrintEvent 100000
/musr/command maximumRunTimeAllowed 86000
/musr/run/randomOption 2
###################################################################################
######################### R O O T O U T P U T ##############################
###################################################################################
/musr/command rootOutput det_edep_el off
/musr/command rootOutput det_edep_pos off
/musr/command rootOutput det_edep_gam off
/musr/command rootOutput det_edep_mup off
/musr/command rootOutput det_nsteps off
/musr/command rootOutput det_length off
/musr/command rootOutput det_time_end off
/musr/command rootOutput det_kine off
/musr/command rootOutput det_VrtxKine off
/musr/command rootOutput det_VrtxX off
/musr/command rootOutput det_VrtxY off
/musr/command rootOutput det_VrtxZ off
/musr/command rootOutput det_VrtxVolID off
/musr/command rootOutput det_VrtxProcID off
/musr/command rootOutput det_VrtxTrackID off
/musr/command rootOutput det_VrtxParticleID off
/musr/command rootOutput det_VvvKine off
/musr/command rootOutput det_VvvX off
/musr/command rootOutput det_VvvY off
/musr/command rootOutput det_VvvZ off
/musr/command rootOutput det_VvvVolID off
/musr/command rootOutput det_VvvProcID off
/musr/command rootOutput det_VvvTrackID off
/musr/command rootOutput det_VvvParticleID off
/musr/command rootOutput odet_ID off
/musr/command rootOutput odet_nPhot off
/musr/command rootOutput odet_timeFirst off
/musr/command rootOutput odet_timeSecond off
/musr/command rootOutput odet_timeThird off
/musr/command rootOutput odet_timeA off
/musr/command rootOutput odet_timeB off
/musr/command rootOutput odet_timeC off
/musr/command rootOutput odet_timeD off
/musr/command rootOutput odet_timeMean off
/musr/command rootOutput odet_timeLast off
/musr/command rootOutput odet_timeCFD off
/musr/command rootOutput odet_amplCFD off
/musr/command rootOutput nOptPhot off
###################################################################################
######################### V I S U A L I S A T I O N ##############################
###################################################################################
/vis/disable
#/control/execute visVRML.mac
#/control/execute visFromToni.mac
#/control/execute visDawn50001.mac
#/vis/open VRML2FILE
#/vis/open DAWNFILE
### (if too many tracks cause core dump => storeTrajectory 0)
#/vis/scene/create
#
#/tracking/storeTrajectory 1
#/vis/viewer/set/viewpointThetaPhi 90 0
##/vis/viewer/set/globalLineWidthScale 3
#/vis/viewer/zoom 30
###/vis/scene/add/trajectories
#/vis/drawVolume
#/vis/viewer/flush
####/hits/verbose 2
###################################################################################
######################### P A R T I C L E G U N #################################
###################################################################################
/gun/vertex 0 0 -1000 mm
# FWHM 10mm ==> sigma = 10/2.354 = 4.2481mm
#/gun/vertexsigma 20 20 0 mm
#---/gun/vertexboundary: rMaxAllowed, zMinAllowed, zMaxAllowed
/gun/vertexboundary 128 -999999 999999 mm
#/gun/momentum 27.0 MeV
# sigma = 3% ==> sigma 27*0.03 = 0.81
#/gun/momentumsmearing 0.81 MeV
#---/gun/momentumboundary: pMinAllowed, pMaxAllowed, dummy
#/gun/momentumboundary 20 40 0 MeV
# TURTLE
#/gun/turtlefilename FOR061_2008_04_22.DAT
#/gun/turtlefilename FOR070_2008_10_17_XXII.DAT
#/gun/turtlefilename FOR077_pie3_HiField_d05_x30.dat
#/gun/turtlefilename FOR077_reggiani_Jan2010_NEW.dat
/gun/turtlefilename FOR077_reggiani_Feb2010.dat
#/gun/turtleZ0position -900 mm
#---/gun/tilt: xangle, yangle, dummy
#/gun/tilt 0 0.5 0 deg
#---/gun/tiltsigma: xangleSigma, yangleSigma, dummy (1 degree on 1 meter ~ 17mm)
#/gun/tiltsigma 0.2 0.2 0 deg
#/gun/pitch 0.5 deg
# Spin rotated by 50 degrees upwards:
/gun/muonPolarizVector 0.766043969 0. -0.642788174
#---/gun/decaytimelimits: decayMin, decayMax, decayTime
#/gun/decaytimelimits 400 410 2197.03 ns
/gun/decaytimelimits -1 -1 2197.03 ns
###################################################################################
######################## B E A M O N #########################################
###################################################################################
#/run/beamOn 3000000
#/run/beamOn 10
/run/beamOn 10000000

1
run/README.TXT Symbolic link
View File

@ -0,0 +1 @@
/afs/psi.ch/user/s/sedlak/geant4/musr/run/README.TXT

10
run/g495.sh Normal file
View File

@ -0,0 +1,10 @@
export G4LEVELGAMMADATA=/home/install/geant4.9.5/data/PhotonEvaporation2.2;
export G4RADIOACTIVEDATA=/home/install/geant4.9.5/data/RadioactiveDecay3.4;
export G4LEDATA=/home/install/geant4.9.5/data/G4EMLOW6.23;
export G4NEUTRONHPDATA=/home/install/geant4.9.5/data/G4NDL4.0;
export G4ABLADATA=/home/install/geant4.9.5/data/G4ABLA3.0;
export G4REALSURFACEDATA=/home/install/geant4.9.5/data/RealSurface1.0;
export G4NEUTRONXSDATA=/home/install/geant4.9.5/data/G4NEUTRONXS1.1;
export G4PIIDATA=/home/install/geant4.9.5/data/G4PII1.3;
export G4VRMLFILE_VIEWER="vrmlview";
echo "On this machine the G4VRMLFILE_VIEWER=$G4VRMLFILE_VIEWER"

BIN
run/run_many_files.tar.gz Normal file
View File

Binary file not shown.