From: Yung-Shun Yeh (adair.py94g@nctu.edu.tw)
Date: Wed Jun 27 2007 - 00:46:16 CEST
Sorry. In the inp file, I forgot to add the direction is negative.
BEAMPOS 0.0 0.0 200.0
NEGATIVE
Even the directon is negative, the result is also the same. And I use
another way ( subroutine bxdraw ) to take data when muon go to layer from
air. It still has no output.
Thanks for your attention!
Best regards,
Yung-Shun
On 6/27/07, Yung-Shun Yeh <adair.py94g@nctu.edu.tw> wrote:
>
> Dear Roesler,
>
> Thanks for your reply. Sorry for my ambiguous description.
>
> In my original simulation, I used a source subroutine to get primary
> muons.
> For a simple case, I supposed a muon beam with energy 100GeV just at +z
> direction (0,0,1000) shoot into the lab(the direction is -z), and the
> detector is at center of geometry. I wish if a muon go through all of the
> three fold layers(region RT, RM, and RB), the MuonTracker would record its
> muon id (NCASE). But I don't know how to do this, I initially recorded the
> informations occured in any three fold layer and liquid scinllator. It seems
> like every muon will be recorded in all three fold layers and liquid
> scintillator, but I only got the informations occured in liquid scinllator.
> There is no interaction between muon and three fold layers. I don't know if
> my configurations of geometry or others are wrong. Or better way can achieve
> this purpose. Thank you very much.
>
> Attached is correct files. In the mgdraw.f file, the code which is used to
> take data of muon is labeled "**** Record the Muon information in specific
> region Start ****". Thank you very much!
>
> Best regards,
> Yung-Shun
>
>
>
>
> On 6/26/07, Stefan Roesler <sroesler@mail.cern.ch> wrote:
> >
> > Dear Yung-Shun
> >
> > I am not sure that I understand your problem. After commenting out the
> > unclosed pre-processor directive just after the ASSIGNMAT cards and
> > commenting the SOURCE card (since you didnt send the source routine) I
> > get
> > entries into all three files BoundaryNeutron, MuonTracker and SRCEPART.
> >
> > Could you please explain your problem in more detail?
> >
> > Best regards
> > Stefan
> >
> >
> >
>
>
This archive was generated by hypermail 2.1.6 : Wed Jun 27 2007 - 09:23:48 CEST