Message: Detector response more deeply!  Not Logged In (login) 
Dear GEANT4 Users! Hi ! I need some help , Can Some body help me in solving my problem! I have a code which is working for Simulations; as GEANT4 code; I am getting the detector response for gamma's.
Now In order to see more deeply detector response or
its efficiency, I need to change( add new files) at least 2 files as
G4Tracking.cc
and G4Stepping action.cc.
In order to see
the sensitivity in the Ist gas gap the sensitivity in the IInd gas gap
the sensitivity in the Ist and IInd gas gap
the last is what we are already getting from the code which is:
the sensitivity in the Ist or IInd gas gap. So dont
need to chnage this part.
In order to see the above 3 detector responses, I need
to have these files change: which are:
 subroutine guout // out put decesion making file as //G4TrackingAction .cc
// ptting flags for decesion making for different //signals
common / iprimary / iflag_prim
common / iflag_1 / iflag1,iflag2
common /iflag_2 /iflag_and,iflag_or common / itracce / nprim,isign1,isign2,isign_and,isign_or // The case of found charged partilce signal (chanrge // particle in the gas gaps) if (iflag_prim.eq.1) then
nprim = nprim + 1
iflag_prim = 0
endif
// put charged partilce in the ist gas gap
if (iflag1.eq.1) then
isign1=isign1+1 iflag1=0
endif
// If particle in the 2nd gas gap put it there if (iflag2.eq.1) then
isign2=isign2+1
iflag2=0
endif // If the particle is between the gas gap I and gas
//Gap II then put it there ...
if(iflag_and.eq.1)then isign_and=isign_and+1
iflag_and=0
endif
// It is our case we have done, particle in both gas
//gaps either Ist or 2nd if(iflag_or.eq.1)then
isign_or=isign_or+1
iflag_or=0
endif
return end
================================================================ SUBROUTINE GUSTEP // making G4stepping action.cc
//like GEANT4
// primary particle signal
common / iprimary / iflag_prim
// putting flags to make signal decesion in 1st gas
gap or II gas gap
common / iflag_1 / iflag1,iflag2
common /iflag_2 /iflag_and,iflag_or if(ngkine.GT.0) call gsking(0) // track found //as // secondary particle
if(ntmult.eq.1.and.ipart.eq.ikine) then // //! primary particle gamma
iflag_prim = 1
endif
// particlle secondarie 2nd dary particles found
if(ipart.eq.ikine) then
if(inwvol.eq.1.and.numed.eq.14) then
endif
endif
if(abs(charge).gt.0) then //! solo particelle // secondarie cariche
// particle in the I gap
if(inwvol.eq.1.and.numed.eq.1) then
iflag1=1
endif
// Particle in the 2nd gas gap , II gap
if(inwvol.eq.1.and.numed.eq.2) then
iflag2=1
endif
endif
//Particle in the First I gas gap AND in the II gas // gap
if (iflag1.eq.1.and.iflag2.eq.1) then
iflag_and=1
endif
// This is the last case what we have found in G4 by our code ...
//.Partice in the Ist gas gap or in the II gas gap
if(iflag1.eq.1) then
iflag_or = 1
endif
END
=================================================================== SUBROUTINE UGLAST
// Its similar to our code G4HistosManager.cc
//Where we have defined the sensitivity ordetector
//response before
// by changing the s //ignals out put in the Ist gas //gap,
in the II gas gap,
//in the I and II gas gap and In the last cas e // Signals in the Ist or 2nd gas gap , one can find //the sensitivity of that gas gap respectively as the // output.
common / itracce /
nprim,isign1,isign2,isign_and,isign_or
10 format(1x,'double gap sensitivity = ',
c + F8.5,1x,'%',1x,'at ene = ',f7.2,1x,'keV')
+ F8.5,1x,'%',1x,'at ene = ',f11.8,1x,'keV')
// total Triggers 100000 like in G4
write(7,*)'Total event entering' ,nprim
// signals in the Ist gas gap
write(7,*)'First gap', isign1
// signals in the IInd gas gap
write(7,*)'Second gap',isign2
// signals in the Ist and II nd gas gap
write(7,*)'first and second ',isign_and
// signals in the Ist or II nd gas gap, no need to //change // it as it is already combined our results, // what we are already getteing good response from // here.
write(7,*)'first or second ',isign_or
sens = float(isign_or)/float(nprim)*100.
write(7,10) sens, 1000.*pkine(1)
END ======================================================
Pls help me in making these files,
God bless you all!
With best Wishes, jamil HEP lab, Seoul, Korea

Inline Depth:  Outline Depth:  Add message: 
to: 