Optimization of Positron Capturing - Revision history

Setisadi: /* Conclusions */

2012-04-12T09:02:33Z

Conclusions

2012-04-12T09:00:55Z

2012-04-12T08:58:22Z

2011-10-07T22:57:35Z

GEANT output

2011-10-07T22:56:17Z

GEANT output

2011-10-07T22:48:15Z

GEANT output

2011-10-07T22:37:47Z

SteppingVerbose

2011-10-07T22:37:36Z

SteppingVerbose

2011-10-07T22:36:46Z

GEANT output

2011-10-07T22:36:12Z

@@ Line 379: / Line 379: @@
 <math>\beta_y = \frac{\sigma_y^2}{\epsilon_y} = \frac{0.00343^2}{0.01403}  = 0.00083121~m</math>
 <math>\gamma_x = \frac{1 + \alpha_x^2}{\beta_x} </math>
 <math> \alpha_x  = \sqrt{\gamma_x \beta_x-1} = \sqrt{37.857 \times 0.00082154 -1} = \sqrt{0.0311 -1} </math>

← Older revision		Revision as of 09:00, 12 April 2012
Line 379:		Line 379:

	<math>\beta_y = \frac{\sigma_y^2}{\epsilon_y} = \frac{0.00343^2}{0.01403} = 0.00083121~m</math>		<math>\beta_y = \frac{\sigma_y^2}{\epsilon_y} = \frac{0.00343^2}{0.01403} = 0.00083121~m</math>
		+
		+	<math>\gamma_x = \frac{1 + \alpha_x^2}{\beta_x} </math>
		+
		+	<math> \alpha_x = \sqrt{\gamma_x \beta_x-1} = \sqrt{37.857 \times 0.00082154 -1} = \sqrt{0.0311 -1} </math>

← Older revision		Revision as of 08:58, 12 April 2012
Line 362:		Line 362:

	N1:N8=27.3		N1:N8=27.3
		+
		+	= Conclusions =
		+
		+	Positrons coming out of W-target:
		+
		+	<math>\epsilon_{x,y} = 14030~mm*mrad</math>
		+
		+	<math>\sigma_x = 3.41 \pm 0.01 ~mm</math>
		+
		+	<math>\sigma_y = 3.43 \pm 0.01 ~mm</math>
		+
		+	<math>\sigma_{x',y'} = 0.732 \pm 0.002 ~rad</math>
		+
		+
		+	<math>\beta_x = \frac{\sigma_x^2}{\epsilon_x} = \frac{0.00341^2}{0.01403} = 0.00082154~m</math>
		+
		+	<math>\beta_y = \frac{\sigma_y^2}{\epsilon_y} = \frac{0.00343^2}{0.01403} = 0.00083121~m</math>

← Older revision		Revision as of 22:57, 7 October 2011
Line 276:		Line 276:

	= GEANT output =		= GEANT output =
		+
		+	Root files and data files at
		+
		+	aztec.iac.isu.edu:/home/setisadi/pos_Tar_Opt/ana/dif_tar_pos

	== SteppingVerbose ==		== SteppingVerbose ==

← Older revision		Revision as of 22:56, 7 October 2011
Line 344:		Line 344:
	[[File:hrrl_pos_sim_g4_output_det8_xp.png \| 300 px]]		[[File:hrrl_pos_sim_g4_output_det8_xp.png \| 300 px]]
	[[File:hrrl_pos_sim_g4_output_det8_yp.png \| 300 px]]		[[File:hrrl_pos_sim_g4_output_det8_yp.png \| 300 px]]
		+
		+	== Ratio of the positrons on detectors ==
		+
		+	NT = 215111
		+
		+	N1 = 103506
		+
		+	N8 = 3787
		+
		+	NT:N1=2.08
		+
		+	NT:N8=56.8
		+
		+	N1:N8=27.3

@@ Line 320: / Line 320: @@
 detector.dat  detector_T.root
 == detector1 ==
@@ Line 326: / Line 330: @@
 detector1.dat  detector1_25mm.root
 == detector8 ==
@@ Line 332: / Line 340: @@
 detector8.dat  detector8_200mm.root

@@ Line 313: / Line 313: @@
 So, the data file contains columns of
- KE x y z Px Py Pz xp yp En DelEn
+KE x y z Px Py Pz xp yp En DelEn
 == detector ==

@@ Line 280: / Line 280: @@
    if( fTrack->GetVolume()->GetName() =="Detector4"  &&  fTrack->GetDefinition()->GetPDGEncoding()==-11 )//e+
+       {
            float Px, Py, Pz, Px2, Py2, Pz2;
            Px = fTrack->GetMomentum().x();
            Py = fTrack->GetMomentum().y();
            Pz = fTrack->GetMomentum().z();
            Px2 = pow(Px, 2);
            Py2 = pow(Py, 2);
            Pz2 = pow(Pz, 2);
 	  float  xp, yp, P2, En, DelEn;
            xp = atan(Px/Pz);
            yp = atan(Py/Pz);
 	  P2 = Px2 + Py2 + Pz2;
 	  En = sqrt ( P2 + 0.511*0.511 );
 	  DelEn = 10 - En;
            if(Px2>0 || Py2>0 || Pz2>0 ){
            <<fTrack->GetKineticEnergy() << " "
            << fTrack->GetPosition().x()<< " "
            << fTrack->GetPosition().y()<< " "
            << fTrack->GetPosition().z()<< " "
            << fTrack->GetMomentum().x() << " "
            << fTrack->GetMomentum().y() << " "
            << fTrack->GetMomentum().z() << " "
 	  << xp << " "
 	  << yp << " "
 	  << En << " "
 	  << DelEn << " "
            << G4endl;
+		}

← Older revision		Revision as of 22:36, 7 October 2011
Line 272:		Line 272:

	[[File:HRRL_e+_pos_opt_sim_20cm_from_q4_P.png \| 300]]		[[File:HRRL_e+_pos_opt_sim_20cm_from_q4_P.png \| 300]]
		+
		+
		+
		+	= GEANT output =
		+
		+	== SteppingVerbose ==
		+
		+	if( fTrack->GetVolume()->GetName() =="Detector4" && fTrack->GetDefinition()->GetPDGEncoding()==-11 )//e+
		+
		+	{
		+	float Px, Py, Pz, Px2, Py2, Pz2;
		+
		+	Px = fTrack->GetMomentum().x();
		+	Py = fTrack->GetMomentum().y();
		+	Pz = fTrack->GetMomentum().z();
		+
		+	Px2 = pow(Px, 2);
		+	Py2 = pow(Py, 2);
		+	Pz2 = pow(Pz, 2);
		+
		+	float xp, yp, P2, En, DelEn;
		+
		+	xp = atan(Px/Pz);
		+	yp = atan(Py/Pz);
		+
		+	P2 = Px2 + Py2 + Pz2;
		+	En = sqrt ( P2 + 0.511*0.511 );
		+	DelEn = 10 - En;
		+
		+	if(Px2>0 \|\| Py2>0 \|\| Pz2>0 ){
		+
		+	outfile_detector4
		+	<<fTrack->GetKineticEnergy() << " "
		+	<< fTrack->GetPosition().x()<< " "
		+	<< fTrack->GetPosition().y()<< " "
		+	<< fTrack->GetPosition().z()<< " "
		+
		+	<< fTrack->GetMomentum().x() << " "
		+	<< fTrack->GetMomentum().y() << " "
		+	<< fTrack->GetMomentum().z() << " "
		+
		+	<< xp << " "
		+	<< yp << " "
		+
		+	<< En << " "
		+	<< DelEn << " "
		+
		+	<< G4endl;
		+	}
		+
		+	So, the data file contains columns of
		+
		+	KE x y z Px Py Pz xp yp En DelEn
		+
		+	== detector ==
		+
		+	Detector right after the W target.
		+
		+	detector.dat detector_T.root
		+
		+	== detector1 ==
		+
		+	Detector 25 mm from W target.
		+
		+	detector1.dat detector1_25mm.root
		+
		+	== detector8 ==
		+
		+	Detector 200 mm from W target.
		+
		+	detector8.dat detector8_200mm.root