Recreation of Plots Using GEMC 2.5 - Revision history

Brenleyt at 16:53, 3 November 2016

2016-11-03T16:53:09Z

Brenleyt at 16:43, 3 November 2016

2016-11-03T16:43:35Z

Brenleyt at 16:33, 3 November 2016

2016-11-03T16:33:06Z

Brenleyt at 18:59, 26 October 2016

2016-10-26T18:59:39Z

Brenleyt at 18:43, 26 October 2016

2016-10-26T18:43:58Z

Brenleyt at 18:43, 26 October 2016

2016-10-26T18:43:27Z

Brenleyt at 18:39, 26 October 2016

2016-10-26T18:39:32Z

Brenleyt at 16:56, 26 October 2016

2016-10-26T16:56:24Z

Brenleyt at 16:51, 26 October 2016

2016-10-26T16:51:48Z

Foretony at 16:51, 21 October 2016

2016-10-21T16:51:53Z

← Older revision		Revision as of 16:53, 3 November 2016
Line 114:		Line 114:

	The command lines for GEMC and CoatJava-2.4 are the same as above. Using the sample root commands and the fit modification within the canvas window to restrict the fit to exclude the tail end of the plot. Here is the xmgrace plot to show my data for the LUND file vs. Dr. Forest's		The command lines for GEMC and CoatJava-2.4 are the same as above. Using the sample root commands and the fit modification within the canvas window to restrict the fit to exclude the tail end of the plot. Here is the xmgrace plot to show my data for the LUND file vs. Dr. Forest's
		+
		+
		+	[[File:GEMC Results LUND Format.png\|200px]]

← Older revision		Revision as of 16:43, 3 November 2016
Line 91:		Line 91:

	Now attempt to recreate the plot using a LUND file with phi equal to 90 degrees (0 degrees if defined in the gcard).		Now attempt to recreate the plot using a LUND file with phi equal to 90 degrees (0 degrees if defined in the gcard).
		+
		+	Here is the gcard I used
		+
		+	<gcard>
		+
		+	<detector name="../../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/dc/dc" factory="TEXT" variation="ccdb"/>
		+	<detector name="../../../../../../clas12/ec/ec" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/ctof/ctof" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/ftof/ftof" factory="TEXT" variation="java"/>
		+	<detector name="../../../../../../clas12/htcc/htcc" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/magnets/solenoid" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/pcal/pcal" factory="TEXT" variation="javageom"/>
		+	<option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
		+	<option name="HALL_FIELD" value="clas12-solenoid"/>
		+	<option name="OUTPUT" value="evio,eg12_sol_75k_1cm_DiskTarg_90DegRot_GEMC_2.5_10_26_16.ev"/>
		+
		+	<!-- <option name="BEAM_P" value="e-, 6.0GeV, 25.0deg, 0*deg"/>-->
		+	<!-- <option name="SPREAD_V" value="(1.0,0.0)cm"/>-->
		+	<option name="RANDOM" value="123483273456"/>
		+	</gcard>
		+
		+	The command lines for GEMC and CoatJava-2.4 are the same as above. Using the sample root commands and the fit modification within the canvas window to restrict the fit to exclude the tail end of the plot. Here is the xmgrace plot to show my data for the LUND file vs. Dr. Forest's

← Older revision		Revision as of 16:33, 3 November 2016
Line 89:		Line 89:

	[[File:GEMC Results LBModifiedFit.png\|200px]]		[[File:GEMC Results LBModifiedFit.png\|200px]]
		+
		+	Now attempt to recreate the plot using a LUND file with phi equal to 90 degrees (0 degrees if defined in the gcard).

← Older revision		Revision as of 18:59, 26 October 2016
Line 86:		Line 86:
	[[File:VzDiff R7 ModifiedFit.png\|200px]]		[[File:VzDiff R7 ModifiedFit.png\|200px]]

−	Here sigma is 0.1861 +/- 0.0015, which agrees identically with Dr. Forest's result.	+	Here sigma is 0.1861 +/- 0.0015, which agrees identically with Dr. Forest's result. Similarly, I can do the same with my first annular section (R1) and with any plot that does not initially agree to get different values that do agree. Using this modified fit method the following plot can be created
		+
		+	[[File:GEMC Results LBModifiedFit.png\|200px]]

← Older revision		Revision as of 18:43, 26 October 2016
Line 86:		Line 86:
	[[File:VzDiff R7 ModifiedFit.png\|200px]]		[[File:VzDiff R7 ModifiedFit.png\|200px]]

−	Here sigma is 0.1861 +/- 0.0015	+	Here sigma is 0.1861 +/- 0.0015, which agrees identically with Dr. Forest's result.

← Older revision		Revision as of 18:43, 26 October 2016
Line 81:		Line 81:

	I haven't perfected setting stat box titles in my macro yet, but the order is that black is 0.0cm-0.2cm, then the stat boxes rotate around clockwise where magenta is 0.8cm-1.0cm.		I haven't perfected setting stat box titles in my macro yet, but the order is that black is 0.0cm-0.2cm, then the stat boxes rotate around clockwise where magenta is 0.8cm-1.0cm.
		+
		+	Despite my reconstruction failing from 0.4cm<Z<0.6cm, I tried to plot only that region and use the fit panel to restrict the region of interest for the fit. The plot is below
		+
		+	[[File:VzDiff R7 ModifiedFit.png\|200px]]
		+
		+	Here sigma is 0.1861 +/- 0.0015

← Older revision		Revision as of 18:39, 26 October 2016
Line 71:		Line 71:

	~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -config GEOM::new=true -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true -config DCTB::useRaster=true -i eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_new.ev -o eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_rec_new.evio		~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -config GEOM::new=true -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true -config DCTB::useRaster=true -i eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_new.ev -o eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_rec_new.evio
		+
		+	Here are the histograms to be compared to mine created on 10/26/16
		+
		+	https://clasweb.jlab.org/wiki/index.php/File:EG12_SweetSpotVertex_ZdiffDist_10-26-16.png

	and here are the plots using the sample cuts above provided by Dr. Forest:		and here are the plots using the sample cuts above provided by Dr. Forest:

@@ Line 57: / Line 57: @@
         <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
         <option name="HALL_FIELD"  value="clas12-solenoid"/>
-        <option name="OUTPUT" value="evio,eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5.ev"/>
+        <option name="OUTPUT" value="evio,eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_new.ev"/>
          <option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>
@@ Line 65: / Line 65: @@
 Here is the GEMC command line:

← Older revision		Revision as of 16:51, 26 October 2016
Line 43:		Line 43:

	[[File:VzDiff Z=0 DeltaR=0.2.png\|200px]]		[[File:VzDiff Z=0 DeltaR=0.2.png\|200px]]
		+
		+	Since these 2 plots were statistically different, I will try to use the same random number generator. Here is the gcard I used:
		+	<gcard>
		+
		+	<detector name="../../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/dc/dc" factory="TEXT" variation="ccdb"/>
		+	<detector name="../../../../../../clas12/ec/ec" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/ctof/ctof" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/ftof/ftof" factory="TEXT" variation="java"/>
		+	<detector name="../../../../../../clas12/htcc/htcc" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/magnets/solenoid" factory="TEXT" variation="original"/>
		+	<detector name="../../../../../../clas12/pcal/pcal" factory="TEXT" variation="javageom"/>
		+	<option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
		+	<option name="HALL_FIELD" value="clas12-solenoid"/>
		+	<option name="OUTPUT" value="evio,eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5.ev"/>
		+
		+	<option name="BEAM_P" value="e-, 6.0GeV, 25.0deg, 0*deg"/>
		+	<option name="SPREAD_V" value="(1.0,0.0)cm"/>
		+	<option name="RANDOM" value="123483273456"/>
		+	</gcard>
		+
		+	Here is the GEMC command line:

@@ Line 3: / Line 3: @@
 https://clasweb.jlab.org/wiki/index.php/File:EG12_SweetSpotVertex_ZdiffDist_10-6-16.png
-using the command
+using the commands
   clas12->Draw("GenPart.z-Event.z >> R1(1000,-5,5)","abs(GenPart.x-Event.x)<1 && abs(GenPart.y-Event.y)<1 && sqrt(GenPart.x*GenPart.x+GenPart.y*GenPart.y)<0.2");
 Here is the gcard used