GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Tracking » Bremsstrahlung treatment in GEANE: how to go back to the original code.
Bremsstrahlung treatment in GEANE: how to go back to the original code. [message #13036] Mon, 06 February 2012 17:51
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *
Hi all,
as promised this morning during the EVO meeting, here are the instructions to go back to the original code of GEANE, removing the changes made to treat the Bremsstrahlung. A detailed and complete treatment of the Bremsstrahlung is under development and will be delivered when ready.

the modifications are in external_packages/transport/geant3/erdecks/erland.F. The following lines should be all commented out (L. 171-188):

C  inclusion of bremsstrahlung fluctuations DEDXB for electrons
C  22.09.2008 by Alberto Rotondi
      IF(XMASS.LT.0.6E-3) THEN
         LX = 1.442695*STEP/RADL
c        EMB=E/(2**LX)
         DEDXB = 0.
c     correction by A. Rotondi and L. Lavezzi to prevent FPE (15/03/2010)
         IF(LX.LT.25) THEN
            S2B=E*E*(1/3**LX - 1/4**LX)
            DEDXB = 1.2*SB
c     end

         DEDX2 = DEDX2 + DEDXB**2

After this, you just have to do "make" in external_packages/transport/geant3 to recompile geant3.

This should allow you to work with the original code, if you want to test it. If you encounter any problem, please let me know.

P.S.: in attachment you can find the file erland.F with the lines already commented out: if you prefer, you can just take it and put it in external_packages/transport/geant3/erdecks/ instead of editing the file yourself. Please, remember to recompile geant3 after this.
  • Attachment: erland.F
    (Size: 6.66KB, Downloaded 290 times)

[Updated on: Mon, 06 February 2012 18:18]

Report message to a moderator

Previous Topic: Multipion with Event Mixing
Next Topic: new task and directory
Goto Forum:

Current Time: Sun Mar 03 10:12:10 CET 2024

Total time taken to generate the page: 0.00961 seconds