62 weighted(true),divideByArea(true)
73 weighted(true),divideByArea(true)
105 if ( fDirection ==
fFlux_InOut || fDirection == dirFlag ){
120 G4double localdirL2 = localdir.
x()*localdir.
x()
121 +localdir.
y()*localdir.
y()
122 +localdir.
z()*localdir.
z();
126 G4double localR2 = localpos1.
x()*localpos1.
x()
127 +localpos1.
y()*localpos1.
y()
128 +localpos1.
z()*localpos1.
z();
129 G4double anglefactor = (localdir.
x()*localpos1.
x()
130 +localdir.
y()*localpos1.
y()
131 +localdir.
z()*localpos1.
z())
132 /std::sqrt(localdirL2)/std::sqrt(localR2);
138 G4double square = radi*radi*dph*( -std::cos(enth) + std::cos(stth) );
142 if ( divideByArea ) current = current/square;
144 current /= anglefactor;
147 EvtMap->
add(index,current);
165 G4double localR2 = localpos1.
x()*localpos1.
x()
166 +localpos1.
y()*localpos1.
y()
167 +localpos1.
z()*localpos1.
z();
172 if ( localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
173 &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
183 G4double localR2 = localpos2.
x()*localpos2.
x()
184 +localpos2.
y()*localpos2.
y()
185 +localpos2.
z()*localpos2.
z();
190 if ( localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
191 &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
221 std::map<G4int,G4double*>::iterator itr = EvtMap->
GetMap()->begin();
222 for(; itr != EvtMap->
GetMap()->end(); itr++) {
223 G4cout <<
" copy no.: " << itr->first
232 if ( divideByArea ) {