DV MollerTrackRecon - Revision history

Vanwdani: /* Phi shifts */

2016-03-31T14:54:38Z

Phi shifts

Vanwdani: /* Transformation Matrix */

2016-03-30T22:28:50Z

Transformation Matrix

Vanwdani: /* Moller Lund */

2016-03-30T17:14:30Z

Moller Lund

Vanwdani: /* Euler Angles */

2016-03-30T16:49:45Z

Euler Angles

Vanwdani: /* =Commands */

2016-03-29T17:19:22Z

=Commands

Vanwdani: /* Commands */

2016-03-29T17:15:21Z

Commands

Vanwdani: /* Phi shifts */

2016-03-29T16:56:53Z

Phi shifts

Vanwdani: /* Alter Phi Angles */

2016-03-27T23:18:56Z

Alter Phi Angles

Vanwdani: /* Alter Phi Angles */

2016-03-26T02:34:49Z

Alter Phi Angles

Vanwdani: /* Alter Phi Angles */

2016-03-25T19:36:46Z

Alter Phi Angles

← Older revision		Revision as of 14:54, 31 March 2016
Line 23:		Line 23:

	</pre>		</pre>
−	~~==Commands==~~
−	~~<pre>gemc -USE_GUI=0 -Hall_Material="Vacuum" -N=1000000 eg12_sol.gcard~~
−	~~</pre>~~
−
−	~~==Results==~~
−	~~gcard to generate electrons.~~
−
−
−	~~<option name="BEAM_P" value="e-, 6.0GeV, 30.0deg, 10*deg">~~
−	~~<option name="SPREAD_P" value="5.5GeV, 25deg, 180*deg">~~
−
−
−	~~[[File:Composite_Fields.png]]~~
−
−
−
−	~~[[File:GeV_graph.png]]~~
−
−
−	~~[[File:MeV_graph.png]]~~
−
−
−	~~[[File:Total_graph.png]]~~

	=Cross-section =		=Cross-section =

← Older revision		Revision as of 22:28, 30 March 2016
Line 1:		Line 1:
−	~~=Transformation Matrix=~~
−
−	~~The Euler angles can be applied using a transformation matrix~~
−
−	~~<math>\left(~~
−	~~\begin{array}{ccc}~~
−	~~\cos (\theta ) & 0 & -\sin (\theta ) \\~~
−	~~0 & 1 & 0 \\~~
−	~~\sin (\theta ) & 0 & \cos (\theta ) \\~~
−	~~\end{array}~~
−	~~\right).\left(~~
−	~~\begin{array}{c}~~
−	~~x \\~~
−	~~y \\~~
−	~~z \\~~
−	~~\end{array}~~
−	~~\right)</math>~~
−
−
−	~~<math>=\left(~~
−	~~\begin{array}{c}~~
−	~~x \cos (\theta )-z \sin (\theta ) \\~~
−	~~y \\~~
−	~~z \cos (\theta )+x \sin (\theta ) \\~~
−	~~\end{array}~~
−	~~\right)</math>~~
−
−
−
−
−	~~For event #29, in sector 3, the location of the first interaction is given by~~
−
−	~~[[File:conversions.png]]~~
−
−
−	~~Converting -25 degrees to radians,~~
−	~~<math>\theta =-0.436332</math>~~
−	~~which is the rotation the detectors are rotated from the y axis.~~
−
−	~~<math>\left(~~
−	~~\begin{array}{ccc}~~
−	~~\cos (\theta ) & 0 & -\sin (\theta ) \\~~
−	~~0 & 1 & 0 \\~~
−	~~\sin (\theta ) & 0 & \cos (\theta ) \\~~
−	~~\end{array}~~
−	~~\right).\left(~~
−	~~\begin{array}{c}~~
−	~~-15.76 \\~~
−	~~0 \\~~
−	~~237.43 \\~~
−	~~\end{array}~~
−	~~\right)</math>~~
−
−	~~<math>=\left(~~
−	~~\begin{array}{c}~~
−	~~86.0588 \\~~
−	~~0. \\~~
−	~~221.845 \\~~
−	~~\end{array}~~
−	~~\right)</math>~~
−
−	~~Finding <math>\phi =\frac{120\ 2 \pi }{360};</math> since "sector -1" =3-1=2*60=120 degrees~~
−
−	~~<math>\left(~~
−	~~\begin{array}{ccc}~~
−	~~\cos (\phi ) & -\sin (\phi ) & 0 \\~~
−	~~\sin (\phi ) & \cos (\phi ) & 0 \\~~
−	~~0 & 0 & 1 \\~~
−	~~\end{array}~~
−	~~\right).\left(~~
−	~~\begin{array}{c}~~
−	~~86.0588 \\~~
−	~~0. \\~~
−	~~221.845 \\~~
−	~~\end{array}~~
−	~~\right)</math>~~
−
−	~~<math>\left(~~
−	~~\begin{array}{c}~~
−	~~-43.0294 \\~~
−	~~74.5291 \\~~
−	~~221.845 \\~~
−	~~\end{array}~~
−	~~\right)</math>~~
−
−	~~This shows how the coordinates are transformed and explains the validity of using the TBTracking information to obtain a phi angle in the lab frame.~~
−
−
−
	=Phi shifts=		=Phi shifts=

← Older revision		Revision as of 17:14, 30 March 2016
Line 1:		Line 1:
−	~~=Moller Lund=~~
−
−
−
−	~~LUND file with Moller events (with origin of coordinates occurring at each event)~~
−	~~<pre>~~
−	~~2 1 1 1 1 0 0.000563654 3.53715 0 6.2002~~
−	~~1 -1 1 11 0 0 0.69 -2.4999 10993.7998 10993.80 0.000511 0 0 0~~
−	~~2 -1 1 11 0 0 -0.69 2.4999 6.5852 7.08 0.000511 0 0 0~~
−	~~</pre>~~
−
−
−	~~From a GEMC run WITH the Solenoid ced is used to obtain the information from the eg12_rec.ev file.~~
−
−	~~[[File:Event29.png]]~~
−
−
−	~~We take the phi angle from the Generated Event momentum as the initial phi angle. The obtain the final phi angle, we can look at the final position of the electron with in the drift chambers.~~
−
−	~~[[File:Detector_position.png]]~~
−
−	~~Examining the position from Timer Based Tracking, we can see that after rotations about first the y-axis, then the z-axis transforms from the detector frame of reference to the lab frame of reference.~~
−
	=Transformation Matrix=		=Transformation Matrix=

← Older revision		Revision as of 16:49, 30 March 2016
Line 21:		Line 21:

	Examining the position from Timer Based Tracking, we can see that after rotations about first the y-axis, then the z-axis transforms from the detector frame of reference to the lab frame of reference.		Examining the position from Timer Based Tracking, we can see that after rotations about first the y-axis, then the z-axis transforms from the detector frame of reference to the lab frame of reference.
−
−	~~=Euler Angles=~~
−
−	~~We can use the Euler angles to perform the rotations.~~
−
−	~~For the rotation about the y axis.~~
−
−	~~[[File:Euler1.png]]~~
−
−	~~And the rotation about the z axis.~~
−
−	~~[[File:Euler2.png]]~~

	=Transformation Matrix=		=Transformation Matrix=

@@ Line 147: / Line 147: @@
 </pre>
-==Commands=
+==Commands==
 <pre>gemc -USE_GUI=0 -Hall_Material="Vacuum" -N=1000000 eg12_sol.gcard
 </pre>

← Older revision		Revision as of 16:56, 29 March 2016
Line 125:		Line 125:
	=Phi shifts=		=Phi shifts=

		+	==gcard==

		+	<pre>
		+	<gcard>
		+
		+	<detector name="../../../../../clas12/fc/forwardCarriage" factory="TEXT$
		+	<detector name="../../../../../clas12/dc/dc" factory="TEXT" $
		+	<detector name="../../../../../clas12/ec/ec" factory="TEXT" $
		+	<detector name="../../../../../clas12/ctof/ctof" factory="TE$
		+	<detector name="../../../../../clas12/ftof/ftof" factory="TE$
		+	<detector name="../../../../../clas12/htcc/htcc" factory="TE$
		+	<detector name="../../../../../clas12/pcal/pcal" factory="TE$
		+	<option name="BEAM_P" value="e-, 6.0GeV, 30.0deg, 10*deg"/>
		+	<option name="SPREAD_P" value="5.5GeV, 25deg, 180*deg"/>
		+	<option name="SCALE_FIELD" value="clas12-torus-big, -1.0"/>
		+	<option name="HALL_FIELD" value="clas12-solenoid"/>
		+	<option name="SCALE_FIELD" value="clas12-solenoid, 1.0"/>
		+	<option name="OUTPUT" value="evio,eg12.ev"/>
		+
		+	</gcard>
		+
		+	</pre>
		+	==Commands==
		+
		+
		+	==Results==
	gcard to generate electrons.		gcard to generate electrons.

@@ Line 274: / Line 274: @@
 <pre>cat set1/Extra_Phi.LUND set2/Extra_Phi.LUND set3/Extra_Phi.LUND set4/Extra_Phi.LUND set5/Extra_Phi.LUND set6/Extra_Phi.LUND set7/Extra_Phi.LUND set8/Extra_Phi.LUND set9/Extra_Phi.LUND set10/Extra_Phi.LUND >Total_Extra_Phi.LUND</pre>
-resulting in a LUND file that was 13294719 lines in length, which equates to 4431573 entries.  This was divided into 177 file parts of 75000 each.  The first set from the original data set is shown below.
+resulting in a LUND file that was 13309755 lines in length, which equates to 4436585 entries.  This was divided into 177 file parts of 75000 each.  The first set from the original data set is shown below.
 <pre>split -a 4 -d -l 75000 Total_Extra_Phi.LUND Phi_Parts_</pre>

@@ Line 326: / Line 326: @@
 Combo->Add(MolThetaCM,8.87e-10);
 Combo->Draw();
-Double_t Bins[15]={90,110,120,130,135,140,142,144,146,148,150,152,154,156,180};
+Double_t Bins[16]={90,100,110,120,130,135,140,142,144,146,148,150,152,154,156,180};
-hnew=Combo->Rebin(14,"hnew",Bins);
+hnew=Combo->Rebin(15,"hnew",Bins);
 hnew->Draw();
 Theory->Draw("same");</pre>

← Older revision		Revision as of 19:36, 25 March 2016
Line 320:		Line 320:

	:::::<math>\frac{1}{\rho_{target}\times l_{target} \times 4\times 10^{10}}=8.87\times 10^{-10} barns=.000887 \mu b</math>		:::::<math>\frac{1}{\rho_{target}\times l_{target} \times 4\times 10^{10}}=8.87\times 10^{-10} barns=.000887 \mu b</math>
		+
		+	Rebining the histogram to account for the unequal weighting of the bins outlined in the table above
		+
		+	<pre>TH1F *Combo=new TH1F("TheoryExperiment","Theoretical and Experimental Differential Cross-Section CM Frame",360,90,180);
		+	Combo->Add(MolThetaCM,8.87e-10);
		+	Combo->Draw();
		+	Double_t Bins[15]={90,110,120,130,135,140,142,144,146,148,150,152,154,156,180};
		+	hnew=Combo->Rebin(14,"hnew",Bins);
		+	hnew->Draw();
		+	Theory->Draw("same");</pre>