Message: Re: Geting charge spectrum Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Geting charge spectrum 

Forum: Event and Track Management
Re: Question Geting charge spectrum (Oscar Larsson)
Re: Feedback Re: Geting charge spectrum (Gumplinger Peter)
Date: 23 Aug, 2010
From: Oscar Larsson <Oscar Larsson>

Tanks for your quick reply.

First, yes, you are right, charg2 should be a PostStepPoint, thanks!

However, your take on my problem was not correct. My problem is this:

When I simulate a genericIon with properties: 26 52 26 (Fe) and 100 events. charge1 says that 100 particles with track ID == 1 entering the scintillator has charge 26 but both charge2 and charge3 say that only 56 particles with charge 26 first exits the aluminium and then later hits the scintillator.

To check this I implemented the following code in my root analysis macro:

if(detID[0] == 26){cout << "detID " << event << "  ";}
for(int i = 0; i < 500; i++){
    if(hullParticle[i][0] == 26){cout << " Hull: " << i ;}
    if(scintParticle[i][0] == 26){cout << " Scint: " << i;}
}
cout << endl;

detID[0] is the value from charge1. hullParticle[i][0] is the value from charge2 and scintParticle[i][0] is the value from charge3.

detID 2

detID 3 Hull: 0 Scint: 0

detID 4

detID 5 Hull: 0 Scint: 0

when I took a closer look at event No 2 I got the following output of charge (Z) and mass (A):

Z-hull A-hull Z-scint A-scint
20 42 20 42
 2 4 2 4
 1 2 1 2
 1 1 1 1
 1 1 1 1
 0 1 1 1
 1 1 1 0
 0 1 1 1
 0 1 0 0
 0 1 0 0
 0 1 0 0
 1 1 0 0
 1 1 0 0
 1 1 0 0

So, in event No 2, charge1 (detID[0]) is set to 26. But when looking at the particle spectrum in this event exiting the hull or entering the scintillator, no particle with charge == 26 can be found. This is why I'm confused.

Why does if(volume == detector->GetScintillator_1() && point1->GetStepStatus() == fGeomBoundary && aStep->GetTrack()->GetTrackID() == 1)

set charge1 to 26 every event when both

if(volume == detector->GetphysiAlum_3() && point2->GetStepStatus() == fGeomBoundary) if(volume == detector->GetScintillator_1() && point1->GetStepStatus() == fGeomBoundary)

show that no particle with charge 26 has exited the Al or entered the scintillator in about 40% of the events??

/Oscar

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Geting charge spectrum   (Oscar Larsson - 25 Aug, 2010)
 Add Message Add Message
to: "Re: Geting charge spectrum"

 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 ]