Message: Detector response more deeply! Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Detector response more deeply! 

Keywords: GEANT4 Tracking events
Forum: Event and Track Management
Date: 18 Jan, 2005
From: jamil <jamil>

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:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

 Add Message Add Message
to: "Detector response more deeply!"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]