DV Plotting XSect - Revision history

Vanwdani at 18:17, 28 June 2016

2016-06-28T18:17:59Z

Vanwdani at 17:58, 28 June 2016

2016-06-28T17:58:35Z

Vanwdani at 17:56, 28 June 2016

2016-06-28T17:56:48Z

Vanwdani at 17:56, 28 June 2016

2016-06-28T17:56:04Z

Vanwdani at 17:51, 28 June 2016

2016-06-28T17:51:36Z

Vanwdani at 17:42, 28 June 2016

2016-06-28T17:42:53Z

Vanwdani at 17:36, 28 June 2016

2016-06-28T17:36:56Z

Vanwdani at 17:36, 28 June 2016

2016-06-28T17:36:27Z

Vanwdani at 17:35, 28 June 2016

2016-06-28T17:35:02Z

Show changes

Vanwdani at 17:26, 28 June 2016

2016-06-28T17:26:44Z

← Older revision		Revision as of 18:17, 28 June 2016
Line 285:		Line 285:

	<center>[[File:XSect_Theory_NH3_11GeV.png]]</center>		<center>[[File:XSect_Theory_NH3_11GeV.png]]</center>
		+
		+
		+	<center>[[File:XSect_LH2_NH3overlay_11GeV.png]]</center>

← Older revision		Revision as of 17:58, 28 June 2016
Line 280:		Line 280:
	}		}
	</pre>		</pre>
		+
		+	<center>[[File:XSect_Theory_LH2_11GeV.png]]</center>
		+
		+
		+	<center>[[File:XSect_Theory_NH3_11GeV.png]]</center>

@@ Line 271: / Line 271: @@
                  tree->Fill();
                  Th->Fill(evt.Angle,evt.dSigma_dOmega);
-−
+         }
   tree->Print();
   tree->Write();

@@ Line 255: / Line 255: @@
    Int_t nlines=0;
-  int count=0;
    TFile *f = new TFile("DiffXSect.root","RECREATE");
    TTree *tree = new TTree("DiffXSect","Moller Diff Xsect");
-   TH1F *Th=new TH1F("Theory","Theoretical Differential Cross Section",360,90,180);
+   TH1F *Th=new TH1F("Theory","Theoretical Differential Cross Section",181,90,180);
    tree->Branch("evt",&evt.event,"event/F:Angle/F:dSigma_dOmega");
@@ Line 267: / Line 268: @@
                  evt.event=nlines;
                  in >> evt.Angle >> evt.dSigma_dOmega;
                   nlines++;
                  tree->Fill();
@@ Line 283: / Line 279: @@
   delete f;
+}

@@ Line 49: / Line 49: @@
 This creates a data file with the scattering angle theta and the corresponding theoretical Moller differential cross-section for said angle.
-</pre>
 A data file is created with the differential cross-section for the scattering angle theta, by increments of .5 degrees
@@ Line 238: / Line 238: @@
 A c++ macro can be written to input this into root
 <pre>
 void MollerDiffXSect2Root()

@@ Line 48: / Line 48: @@
 </pre>
-This creates a data file with the scattering angle theta and the corresponding theoretical Moller differential cross-section for said angle.  A c++ macro can be written to input this into root
+This creates a data file with the scattering angle theta and the corresponding theoretical Moller differential cross-section for said angle.
 </pre>
@@ Line 277: / Line 235: @@
 .000000     inf
 </pre>

@@ Line 93: / Line 93: @@
 </pre>
-A data file is created with the differential cross-section for the scattering angle theta, by increments of 1
+A data file is created with the differential cross-section for the scattering angle theta, by increments of .5 degrees
 <pre>
 .000000      0.000008

@@ Line 95: / Line 95: @@
 A data file is created with the differential cross-section for the scattering angle theta, by increments of 1
 <pre>
 .000000      0.000008
 .500000      0.000008
 .000000      0.000008
@@ Line 135: / Line 135: @@
 .000000     0.000011
 .500000     0.000011
 .000000     0.000011
 .500000     0.000012
 .000000     0.000012
@@ Line 173: / Line 173: @@
 .000000     0.000027
 .500000     0.000028
 .000000     0.000029
 .500000     0.000030
 .000000     0.000031
@@ Line 211: / Line 211: @@
 .000000     0.000144
 .500000     0.000152
 .000000     0.000162
 .500000     0.000172
 .000000     0.000183
@@ Line 249: / Line 249: @@
 .000000     0.004183
 .500000     0.004834
 .000000     0.005618
 .500000     0.006568
 .000000     0.007727

@@ Line 8: / Line 8: @@
 float angle,dSigma_dOmega;
 double alpha=7.2973525664e-3;
-double E_CM=53e6;
+double E_CM=53e6;                                               //CM E for 11GeV incident electrons
                  /* using equation from Landau-Lifshitz & Azfar */
                  /* units in MeV */
-                 for(angle=90;angle<=180; angle=angle+.25)
+                 for(angle=90;angle<=180; angle=angle+.5)
+                         {
-                         dSigma_dOmega=(alpha*(3+cos(angle*3.14/180)*cos(angle*3.14/180)));
+                         dSigma_dOmega=(alpha*(3+cos(angle*3.14159265359/180)*cos(angle*3.14159265359/180)));
                          dSigma_dOmega=dSigma_dOmega*dSigma_dOmega;
-dSigma_dOmega=dSigma_dOmega/(4*E_CM*E_CM*sin(angle*3.14/180)*sin(angle*3.14/180)*sin(angle*3.14/180)*sin(angle*3.14/180));
+                        dSigma_dOmega=dSigma_dOmega/(4*E_CM*E_CM*sin(angle*3.14159265359/180)*sin(angle*3.14159265359/180)
                          dSigma_dOmega=dSigma_dOmega*1e18;
-                         dSigma_dOmega=dSigma_dOmega*.3892e-3;//barns
+                         dSigma_dOmega=dSigma_dOmega*.3892e-3;   //barns
-                         dSigma_dOmega=dSigma_dOmega/1e-6;//micro barns
+                         dSigma_dOmega=dSigma_dOmega/2;          //2 electrons for LH2
                          fprintf(stdout," %f\t%f\n", angle, dSigma_dOmega);
+                         }