Difference between revisions of "DV RunGroupC Moller"

From New IAC Wiki
Jump to navigation Jump to search
 
(222 intermediate revisions by 2 users not shown)
Line 2: Line 2:
  
 
Simulating the Moller scattering background for EG12
 
Simulating the Moller scattering background for EG12
 +
=Docker=
 +
==[[Set up Docker Container]]==
 +
=[[DV_XSECT|Moller Differential Cross-Section]]=
  
 
=GEANT4 Simulation of Moller Events=
 
=GEANT4 Simulation of Moller Events=
  
−
==Simulation Setup==
+
== [[LH2 target|Simulation Setup]]==
−
Determine the Moller background using an LH2 target to check the physics in GEANT4
+
Use GEANT4 via GEMC to estimate the Moller background for electron scattering experiments in JLab's Hall B.  The first step towards this goal is to use GEANT4 without the GEMC infrastructure to create event (LUND) files that will be used as input events for GEMC.
  
−
===Detector Material and Construction===
+
==[[DV_Moller_LH2 | Benchmark GEANT4's Moller scattering prediction with the theoretical cross section using LH2]]==
−
Using GEANT4, the '''ExampleN02''' file was edited to run for a NH3 target.  The file '''ExN02DetectorConstruction.cc''' was edited with
 
  
−
<pre>
 
−
//--------- Material definition ---------
 
−
 
 
−
  G4double a, z;
 
−
  G4double density, temperature, pressure;
 
  
−
//Liquid Hydrogen
 
−
G4Material* LH2 = new G4Material("Hydrogen", z=2., a=2.02*g/mole,  density=0.07*g/cm3, kStateGas,3*kelvin,1.7e5*pascal);
 
−
</pre>
 
  
−
The target is a cylinder with a 1.5 cm diameter and 1 cm thickness following dimensions listed on page 8 of [[File:PHY02-33.pdf ]]
+
===Comparison of simulation vs. the theoretical Møller differential cross section using 11 GeV electrons impinging LH2===
  
−
<pre>
+
[[Converting to barns|Converting the number of scattered electrons per scattering angle theta to a differential cross-section in barns.]]
−
//--------- Sizes of the principal geometrical components (solids)  ---------
 
  
−
  NbOfChambers = 1;
+
[[File:XSect_new_zoom.png|frame|center|alt=Experimental and Theoretical Moller Differential Cross-Section in Center of Mass Frame Frame|'''Figure 5c:''' The experimental and theoretical Moller electron differential cross-section for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.]]
−
  ChamberWidth = 1.5*cm;
 
−
  ChamberSpacing = 40*cm;
 
−
 
 
−
  fTrackerLength = (NbOfChambers+1)*ChamberSpacing; // Full length of Tracker
 
−
  fTargetLength  = 1.0 * cm;                        // Full length of Target
 
  
−
  TargetMater = NH3;
+
=Change to a NH3 Target=
−
  ChamberMater = BadVacuum;
+
[[Replacing the LH2 target with an NH3 target]]
−
 
 
−
 
 
−
  //fWorldLength= 1.2 *(fTargetLength+fTrackerLength);
 
−
  fWorldLength= 1.2 *(10+fTrackerLength)+100 *cm;
 
−
 
 
−
  G4double targetSize  = 0.5*fTargetLength;    // Half length of the Target
 
−
  G4double trackerSize = 0.5*fTrackerLength;  // Half length of the Tracker
 
−
</pre>
 
  
−
===Recording Moller Events===
+
==[[DV_Moller_NH3|Benchmark GEANT4's Moller scattering prediction with the theoretical cross section using NH3]]==
−
Moller scattering does not have a GEANT4 specific setting, so the file '''ExN02SteppingVerbose.cc''' was edited to record only events that would correspond to Moller events.  Since Moller scattering is simply electron-electron scattering, we can look for events where electron ionization is the process in the data stream.  We will limit this to only particles with a parentID of 0 and 1 representing the parent and daughter (Moller) electrons.  This eliminates second generation Moller scatterings, which only ocur about 2 times out of 1E6 incoming electrons.  This method also eliminates Delta Rays or Knock-on-Electrons.
 
  
 +
==Comparison of simulation vs. the theoretical Møller differential cross section using 11 GeV electrons impinging NH3==
  
−
The file ExN02SteppingVerbose.cc is read for each step in the GEANT4 simulation, so recording the momentum, position, and energies of the electrons before and after the collision can be found in multiple loops.
+
[[Converting to barns|Converting the number of scattered electrons per scattering angle theta to a differential cross-section in barns.]]
  
 +
[[File:XSect_NH3.png|frame|center|alt=Theoretical and Simulated Moller Differential Cross-Section in Center of Mass Frame Frame|'''Figure 5c:''' The theoretical and simulated Moller electron differential cross-section for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference for NH3 target.]]
  
−
On the first pass of SteppingVerbose, the data is recorded into a temporary variable set.  The physical process of ionization occurs after the collision of the two electrons.  This implies that the data could possible be the initial state of the incoming electron.  This is read every time to be prepared for the subsequent pass where the process of ionization (scattering) is active.
+
==LH2 Vs. NH3==
−
<pre>
+
===[[DV_Moller_NH3_LH2|Benchmark GEANT4's Moller scattering prediction for NH3 and LH2]]===
−
void ExN02SteppingVerbose::StepInfo()
 
−
{
 
−
.
 
−
.
 
−
.
 
−
if(fTrack->GetDefinition()->GetPDGEncoding()==11 && fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName()!="eIoni" && fTrack->GetParentID()==0)
 
−
    {
 
−
        Temp_Energy=fTrack->GetKineticEnergy();
 
−
        Temp_Mom_x=fTrack->GetMomentum().x();
 
−
        Temp_Mom_y=fTrack->GetMomentum().y();
 
−
        Temp_Mom_z=fTrack->GetMomentum().z();
 
−
        Temp_Pos_x=fTrack->GetPosition().x();
 
−
        Temp_Pos_y=fTrack->GetPosition().y();
 
−
        Temp_Pos_z=fTrack->GetPosition().z();
 
−
    }
 
−
</pre>
 
  
−
On a pass afterwards, the data is read into a variable for the final state when the physical state is ionization, representing the Moller scattering.  The the temporary variable is read into the initial state. 
 
  
−
<pre>
+
==Effects Due to Target Material==
−
if( fTrack->GetDefinition()->GetPDGEncoding()==11 && fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName()=="eIoni" && fTrack->GetParentID()==0)
+
===[[DV_Target_Density|Target Density]]===
−
    {
 
−
        Final_Energy= fTrack->GetKineticEnergy();
 
−
        Final_Mom_x=fTrack->GetMomentum().x();
 
−
        Final_Mom_y=fTrack->GetMomentum().y();
 
−
        Final_Mom_z=fTrack->GetMomentum().z();
 
−
        Final_Pos_x=fTrack->GetPosition().x();
 
−
        Final_Pos_y=fTrack->GetPosition().y();
 
−
        Final_Pos_z=fTrack->GetPosition().z();
 
−
     
 
−
        Init_Energy=Temp_Energy;
 
−
        Init_Mom_x=Temp_Mom_x;
 
−
        Init_Mom_y=Temp_Mom_y;
 
−
        Init_Mom_z=Temp_Mom_z;
 
−
        Init_Pos_x=Temp_Pos_x;
 
−
        Init_Pos_y=Temp_Pos_y;
 
−
        Init_Pos_z=Temp_Pos_z;
 
−
       
 
−
        }
 
−
</pre>
 
  
−
This will write the data to an external file only for a 1st generation daughter particle and an active trigger.  Afterwards, the trigger is turned off so that the recording process can start again.
 
  
−
<pre>
+
===[[DV_a_z_target|Atomic Mass and Electron Number Effects]]===
−
if(fTrack->GetDefinition()->GetPDGEncoding()==11 && fTrack->GetParentID()==1 && trigger==1 )
 
−
  {
 
−
    outfile
 
−
    //G4cout
 
−
      << Init_Energy<< "    "
 
−
      <<  Init_Mom_x << "    "
 
−
      <<  Init_Mom_y << "    "
 
−
      <<  Init_Mom_z << "    "
 
−
      <<  Init_Pos_x << "    "
 
−
      <<  Init_Pos_y << "    "
 
−
      <<  Init_Pos_z << "    "
 
−
      << Final_Energy << "    "
 
−
      <<  Final_Mom_x << "    "
 
−
      <<  Final_Mom_y << "    "
 
−
      <<  Final_Mom_z << "    "
 
−
      <<  Final_Pos_x << "    "
 
−
      <<  Final_Pos_y << "    " 
 
−
      <<  Final_Pos_z << "    "
 
−
      << Mol_Energy << "    "
 
−
      <<  Mol_Mom_x << "    "
 
−
      <<  Mol_Mom_y << "    "
 
−
      <<  Mol_Mom_z << "    "
 
−
      <<  Mol_Pos_x << "    "
 
−
      <<  Mol_Pos_y << "    "
 
−
      <<  Mol_Pos_z << "    "
 
−
      << G4endl;
 
−
  trigger=0;
 
−
  }
 
−
.
 
−
.
 
−
.
 
−
}
 
−
</pre>
 
  
−
On a later pass, the condition that the parentID no longer represents the parent.  This implies that the particle a Moller electron and the data is recorded into the Moller final state.  The trigger is activated, which on the next pass of the program will allow a printout of all Moller Scattering data.
+
==[[Calculating the differential cross-sections for the different materials, and placing them as well as the theoretical differential cross-section into a plot:|Differential Cross-Section Offset]]==
  
−
<pre>
+
[[File:Adjusted_MollerXSect_NH3.png‎|frame|center|alt=Theoretical and Simulated Moller Differential Cross-Section in Center of Mass Frame Frame|'''Figure 5c:''' The theoretical and simulated Moller electron differential cross-section for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference for NH3 target. The theoretical differential cross-section has been adjusted to account for the 10 electrons within NH3.]]
−
void ExN02SteppingVerbose::TrackingStarted()
 
−
{
 
−
.
 
−
.
 
−
.
 
−
if(fTrack->GetDefinition()->GetPDGEncoding()==11 && fTrack->GetParentID()>0)
 
−
  {
 
−
    Mol_Energy=fTrack->GetKineticEnergy();
 
−
    Mol_Mom_x=fTrack->GetMomentum().x();
 
−
    Mol_Mom_y=fTrack->GetMomentum().y();
 
−
    Mol_Mom_z=fTrack->GetMomentum().z();
 
−
    Mol_Pos_x=fTrack->GetPosition().x();
 
−
    Mol_Pos_y=fTrack->GetPosition().y();
 
−
    Mol_Pos_z=fTrack->GetPosition().z();
 
−
    trigger=1;
 
−
  }
 
−
.
 
−
.
 
−
.
 
−
}
 
−
</pre>
 
  
 +
[[ROOT Macro to read LUND files and make plots]]
  
−
===Working with Moller Data===
 
−
The event data from the GEANT4 simulation is written to a data file in the following format:
 
  
−
{| border="1"
+
<center>[[File:LUND_MolMomCM.png]][[File:LUND_MolThetaCM.png]]</center>
−
|-
 
−
! KE<sub>i</sub>
 
−
! Px<sub>i</sub>  
 
−
! Py<sub>i</sub>
 
−
! Pz<sub>i</sub>
 
−
! x<sub>i</sub>
 
−
! y<sub>i</sub>
 
−
! z<sub>1</sub>
 
−
! KE<sub>f</sub>
 
−
! Px<sub>f</sub>
 
−
! Py<sub>f</sub>
 
−
! Pz<sub>f</sub>
 
−
! x<sub>f</sub>
 
−
! y<sub>f</sub>
 
−
! z<sub>f</sub>
 
−
! KE<sub>m</sub>
 
−
! Px<sub>m</sub>
 
−
! Py<sub>m</sub>
 
−
! Pz<sub>m</sub>
 
−
! x<sub>m</sub>
 
−
! y<sub>m</sub>
 
−
! z<sub>m</sub>
 
−
|-
 
−
| 11000
 
−
| 0
 
−
| 0
 
−
| 11000.5
 
−
| 0
 
−
| 0
 
−
| -510
 
−
| 10999.1
 
−
| 0.433025
 
−
| -0.858867
 
−
| 10999.6
 
−
| 0
 
−
| 0
 
−
| -509.276
 
−
| 0.905324
 
−
| -0.433025
 
−
| 0.858867
 
−
| 0.905366
 
−
| 0
 
−
| 0
 
−
| -509.276
 
−
|}
 
  
−
Where each line represents a Moller scattering with the kinematic variables of kinetic energy, Momentum in the x, y, and z directions, as well as the x, y, and z position of the collision for the Incoming electron, the scattered state of the incoming electron, and the scattered Moller electron in that specific order.
 
  
 +
<center>[[File:LUND_MolMomLab.png]][[File:LUND_MolThetaLab.png]]</center>
  
−
Using a c++ macro, the variables are read into a Root tree, with branches for each variable.  Specific histograms are also created for the total momentum and the scattering angle theta for the scattered and Moller electron, both in the Lab frame and Center of Mass frame.
 
  
−
<pre>
+
Run the above theta and E distribution of Mollers through GEMC with and without Solenoid on. Determine the Theta and E range of Mollers that enter the detector.
−
tree->Branch("evt",&evt.event,"event/I:IntKE/F:IntPx:IntPy:IntPz:IntPosx:IntPosy:IntPosz:FnlKE:FnlPx:FnlPy:FnlPz:FnlPosx:FnlPosy:FnlPosz:
 
−
                MolKE:MolPx:MolPy:MolPz:MolPosx:MolPosy:MolPosz");
 
−
        while(in.good())
 
−
                {
 
−
//Create Tree from GEANT4 simulation data
 
−
                evt.event=nlines;
 
−
                in >> evt.IntKE >> evt.IntMom[0] >> evt.IntMom[1] >> evt.IntMom[2]  >> evt.IntPos[0]
 
−
                >> evt.IntPos[1] >> evt.IntPos[2] >> evt.FnlKE >> evt.FnlMom[0] >> evt.FnlMom[1] >> evt.FnlMom[2]
 
−
                >> evt.FnlPos[0] >> evt.FnlPos[1] >> evt.FnlPos[2] >> evt.MolKE >> evt.MolMom[0] >> evt.MolMom[1]
 
−
                >> evt.MolMom[2] >> evt.MolPos[0] >> evt.MolPos[1] >> evt.MolPos[2];
 
−
           
 
−
                nlines++;
 
  
−
                tree->Fill();
 
  
−
                FnlE=sqrt(evt.FnlMom[0]*evt.FnlMom[0]+evt.FnlMom[1]*evt.FnlMom[1]+evt.FnlMom[2]*evt.FnlMom[2]+0.511*0.511);
+
Scanning from 2MeV to 5500MeV, and 0 to 36 degrees in Theta.
−
                IntE=sqrt(evt.IntMom[0]*evt.IntMom[0]+evt.IntMom[1]*evt.IntMom[1]+evt.IntMom[2]*evt.IntMom[2]+0.511*0.511);
 
−
                MolE=sqrt(evt.MolMom[0]*evt.MolMom[0]+evt.MolMom[1]*evt.MolMom[1]+evt.MolMom[2]*evt.MolMom[2]+0.511*0.511);//Define 4Vectors
 
−
                Fnl4Mom.SetPxPyPzE(evt.FnlMom[0],evt.FnlMom[1],evt.FnlMom[2],FnlE);
 
−
                Int4Mom.SetPxPyPzE(evt.IntMom[0],evt.IntMom[1],evt.IntMom[2],IntE);
 
−
                Mol4Mom.SetPxPyPzE(evt.MolMom[0],evt.MolMom[1],evt.MolMom[2],MolE);
 
−
                 
 
−
        //Create Lab Frame Histograms
 
−
                FinalMomentum->Fill(Fnl4Mom.P());
 
−
                MollerMomentum->Fill(Mol4Mom.P());
 
−
                FinalTheta->Fill(Fnl4Mom.Theta()*180/3.14);
 
−
                MollerTheta->Fill(Mol4Mom.Theta()*180/3.14);
 
−
        //Boost to Center of Mass Frame
 
−
                CMS=Fnl4Mom+Mol4Mom;
 
−
                Fnl4Mom.Boost(-CMS.BoostVector());
 
−
             
 
−
                Mol4Mom.Boost(-CMS.BoostVector());
 
−
        //Create CM Histograms
 
−
                FinalMomentumCM->Fill(Fnl4Mom.P());
 
−
                MollerMomentumCM->Fill(Mol4Mom.P());
 
−
                FinalThetaCM->Fill(Fnl4Mom.Theta()*180/3.14);
 
−
                MollerThetaCM->Fill(Mol4Mom.Theta()*180/3.14);
 
−
                }
 
−
</pre>
 
  
−
=='''[[Momentum distributions in the Lab Frame.]]''' ==
+
For Solenoid on with 5T:
  
 +
<pre>
 +
<gcard>
  
 +
        <detector name="../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
 +
        <detector name="../../../../../clas12/dc/dc"            factory="TEXT" variation="original"/>
 +
        <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="original"/>
 +
        <detector name="../../../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
 +
        <detector name="../../../../../clas12/pcal/pcal"            factory="TEXT" variation="original"/>
 +
        <option name="BEAM_P"  value="e-, 2.800*GeV, 18.0*deg, 10*deg"/>
 +
        <option name="SPREAD_P" value="2.798*GeV, 18*deg, 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"/>
  
−
[[File:FnlMom.png|frame|center|alt=Scattered Electron Momentum in Lab Frame|Figure 1a: The scattered electron momentum distribution for an incident 11 GeV electron in the Lab frame of reference.]]
+
</gcard>
  
 +
</pre>
  
 +
Using the standard commands for gemc
  
−
[[File:MolMom.png|frame|center|alt=Moller Electron Momentum in Lab Frame|Figure 1b: The Moller electron momentum distribution for an incident 11 GeV electron in the Lab frame of reference.]]
+
<pre>
 +
gemc -USE_GUI=0 -N=1000 eg12_sol.gcard
 +
~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -i eg12.ev -config DCHB::torus=-1.0 -config DCHB::solenoid=1.0 -config DCTB::kalman=true -o eg12_rec.ev -s DCHB:DCTB:EC:FTOF:EB
 +
~/src/CLAS/coatjava-1.0/bin/run-groovy Analysis.groovy eg12_rec.0.evio
 +
</pre>
  
−
==Angular Distribution in the Lab Frame==   
 
  
 +
Checking for a reconstructed particle that undergoes a phi shift where the 1st column is energy and the 2nd the scattering angle theta:
 +
<pre>
 +
For Energy:
 +
awk 'NR == 1 {line =$0; min =$2} NR >1 && $2 < min {line =$0; min = $2} END{print line}' Energy_Phi_Shift.dat
 +
289    0.6624017631112946      33.37507824033966      -179.63657653407975    0.7254709      29.696417      -177.79057      473.36051689004984      11.978802748439216      -129.55774724876173
  
−
[[File:FnlTheta.png|frame|center|alt=Scattered Electron Angle Theta in Lab Frame|Figure 2a: The scattered electron scattering angle theta distribution for an incident 11 GeV electron in the Lab frame of reference.]]
 
  
 +
awk 'NR == 1 {line =$0; max =$2} NR >1 && $2 > max {line =$0; max = $2} END{print line}' Energy_Phi_Shift.dat
 +
978    5.584765946130235      24.90893565820338      -170.69784116534063    5.593401        24.90894        -170.65994      472.16552666546283      23.31807917059094      -165.11125601370114
  
 +
For Theta:
 +
awk 'NR == 1 {line =$0; min =$3} NR >1 && $3 < min {line =$0; min = $3} END{print line}' Energy_Phi_Shift.dat
 +
12      4.83485521395583        10.577021312967325      -50.86708775351051      4.8205996      10.617784      -50.895683      491.627055439108        5.9032309421001985      -43.82964507576243
  
−
[[File:MolTheta.png|frame|center|alt=Moller Electron Angle Theta in Lab Frame|Figure 2b: The Moller electron scattering angle theta distribution for an incident 11 GeV electron in the Lab frame of reference.]]
+
awk 'NR == 1 {line =$0; max =$3} NR >1 && $3 > max {line =$0; max = $3} END{print line}' Energy_Phi_Shift.dat
 +
187    1.3382869250001856      35.63598792463312      29.887530073265214      1.3481133      35.449417      30.009071      472.455052447561        19.034043420497742      55.20328358414997
 +
</pre>
  
−
=='''[[Momentum distributions in the Center of Mass Frame.]]'''==
 
−
===Estimated Momentum Distribution===
 
  
−
Noting that the GEANT4 data included data for the Kinetic Energy and not the total energy, we can perform the substitution using
+
Similarly for 0T:
 +
<pre>
 +
For Energy:
 +
awk 'NR == 1 {line =$0; min =$2} NR >1 && $2 < min {line =$0; min = $2} END{print line}' Energy_Phi_Shift.dat
 +
148    0.36418038435234207    33.68296595663359      -175.89942431595827    0.35383263      51.743824      -39.45801      469.5694572319308      24.353383233947522      -176.27025218702502
  
−
<center><math>E_1^2=p_1^2+m_1^2</math></center>
 
  
−
Performing [[DV_Calculations_of_4-momentum_components|the relativistic transformation]], we find for an incoming electron with momentum of 11GeV, we should find the momentum in the center of mass to be around '''53 MeV''' .
+
awk 'NR == 1 {line =$0; max =$2} NR >1 && $2 > max {line =$0; max = $2} END{print line}' Energy_Phi_Shift.dat
 +
818    5.5278794590742235      12.937136966282836      -43.35663621477178      5.539351        12.95216        -49.09295      484.8800252169308      8.27731670021694        -43.245753766203364
  
 +
For Theta:
 +
awk 'NR == 1 {line =$0; min =$3} NR >1 && $3 < min {line =$0; min = $3} END{print line}' Energy_Phi_Shift.dat
 +
559    4.909449160222199      9.918634647646794      0.06548256930681198    4.914272        9.938521        -5.861595      492.1448533455858      7.501204563534799      -0.006721611135017764
  
 +
awk 'NR == 1 {line =$0; max =$3} NR >1 && $3 > max {line =$0; max = $3} END{print line}' Energy_Phi_Shift.dat
 +
326    1.272966353429713      35.731680814303196      85.3173707349696        1.2830343      38.27654        58.368687      483.4835073848272      3.950855595379657      83.87603767806992
 +
</pre>
  
−
[[File:FnlMomCM_new.png|frame|center|alt=Scattered Electron Momentum in Center of Mass Frame|Figure 3a: The scattered electron momentum distribution for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.]]
+
=Reconstruction of Moller Events=
  
 +
==[[DV_Creating_LUND_Files|The LUND format]]==
  
 +
==[[DV_Running_LUND_for_Moller_Distribution|Writing LUND files]]==
  
−
[[File:MolMomCM_new.png|frame|center|alt=Moller Electron Momentum in Center of Mass Frame|Figure 3b: The Moller electron momentum distribution for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.]]
+
==[[DV_Running_GEMC|Running GEMC]]==
  
−
==Angular Distribution in the Center of Mass Frame==
+
==[[DV_Analyze_Recon|Phi Shift observation using DC hit Reconstruction Data]]==
  
 +
=Gcard creation defining energy and angle range of electrons=
 +
==[[Modified gcards]]==
  
−
[[File:FnlThetaCM.png|frame|center|alt=Scattered Electron Angle Theta in Center of Mass Frame|Figure 4a: The scattered electron scattering angle theta distribution for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.]]
+
=Effects of Solenoid on Electrons=
 +
==[[Results for known Moller events|Results for defined distribution in Solenoid Fields]]==
  
 +
==[[Results for Random Spread of Energy and angle theta in the Lab frame]]==
  
 +
=Cover Full Solid Angle of Detector=
  
−
[[File:MolThetaCM.png|frame|center|alt=Moller Electron Angle Theta in Center of Mass Frame Frame|Figure 4b: The Moller electron scattering angle theta distribution for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.]]
+
==Using GEANT simulation data==
 +
===[[Calculations of 4-momentum components]]===
  
−
=='''[[Comparing experimental vs. theoretical for Møller differential cross section]]''' 11GeV==
+
===[[Alter Phi Angles]]===
  
 +
===[[Check Differential Cross-Section]]===
  
−
Using the equation from [1]
 
  
−
<center><math>\frac{d\sigma}{d\Omega}=\frac{ e^4 }{8E^2}\left \{\frac{1+cos^4\frac{\theta}{2}}{sin^4\frac{\theta}{2}}+\frac{1+sin^4\frac{\theta}{2}}{cos^4\frac{\theta}{2}}+\frac{2}{sin^2\frac{\theta}{2}cos^2\frac{\theta}{2}} \right \}</math></center>
 
  
  
−
<center>where <math>\alpha=\frac{e^2}{\hbar c}\quad with\quad \hbar = c =1</math></center>
+
==Detector Occupancy==
  
 +
clas12->Draw("Detector.wire:Detector.layer>>(7,1,7,120,0,120)","Detector.superlayer<3","colz")
  
−
This can be simplified to the form
+
===[[gcard settings for daughter and procID|Gcard settings]]===
  
 +
===[[Verfication of Mother/Daughter Occupancy]]===
  
−
<center><math>\frac{d\sigma}{d\Omega}=\frac{ \alpha^2 }{4E^2}\frac{ (3+cos^2\theta)^2}{sin^4\theta}</math></center>
+
===[[Benchmark GEMC Occupancy Prediction for 11GeV Electron Beam with 0T Solenoid for LH2]]===
  
−
Plugging in the values expected for a scattering electron:
 
  
 +
===[[Setup for Batch Job With Varying Experimental CLAS12 Quantities]]===
 +
===[[Run_in_GEMC]]===
 +
===[[Center_of_Mass_for_Stationary_Target]]===
  
 +
===[[Run Occupancy for Sector 1 DC hits]]===
 +
 
 +
===[[Wire_angle_correspondance]]===
 +
===[[Isotropic Weighted Moller Distribution in GEMC]]===
  
 +
=Papers used=
  
−
<center><math>\alpha ^2=5.3279\times 10^{-5}</math></center>
+
[1]Farrukh Azfar's Derivation of Moller Scattering
  
 +
:::[[File:FarrukAzfarMollerScatter.pdf]]
  
−
<center><math>E\approx 53 MeV</math></center>
+
A polarized target for the CLAS detector
  
 +
:::[[File:PHY02-33.pdf‎]]
  
−
Using unit analysis on the term outside the parantheses, we find that the differential cross section for an electron at this momentum should be around
+
An investigation of the spin structure of the proton in deep inelastic scattering of polarized muons on polarized protons
  
−
<center><math>\frac{5.3279\times 10^{-5}}{4\times 2.81\times 10^{15}eV^2}=4.74\times 10^{-21} eV^{-2}=\frac{4.74\times 10^{-21}}{1eV^2}\times \frac{1\times 10^{18}eV^2 }{1\ GeV^2}=\frac{.0047}{GeV^2}</math></center>
+
:::[[:File:1819.pdf]]
  
−
Using the conversion of
+
==QED Radiative Corrections to Low-Energy Moller and Bhabha Scattering==
 +
http://arxiv.org/abs/1602.07609
  
 +
[[File:ChangingRates_S1_PhiThetaHits.png | 600 px]]
  
−
<center><math>\frac{1}{1GeV^2}=.3894 mb</math></center>
+
[[File:S1_50nA_PrimaryElectronSigmasWeightedRates.png | 600 px]]
  
 +
[[File:S1_PhiThetaGammaHits.png | 600 px]]
  
−
<center><math>\frac{1\ GeV^{-2}}{.3894\ mb}=\frac{.0047\ GeV^{-2}}{x\ mb}</math></center>
 
  
−
The trigonometric function part of the equation comes out to it's minimum of 9 at 90 degrees.
+
[[File:HitMakeUp.png]]
  
−
 
−
<center><math>\frac{(3+Cos^2(90))^2}{Sin^4(90)}=9</math></center>
 
  
 +
[[File:S1_50nA_PrimaryElectronSigmasWeightedRates_Full.png | 600 px]]
  
−
We find that the differential cross section scale is <math>\frac{d\sigma}{d\Omega}\approx 1.8\times 10^{-3}mb \times 9=16.2\mu b</math>
+
[[File:ComparingOppositeFields_S1_PhiThetaHits.png | 800 px]]
  
  
 +
[[File:ChangingSolenoidRates_wo_Magnets.png|800 px]]
  
−
Converting the number of electrons to barns,
 
  
−
:::::<math>L=\frac{i_{scattered}}{\sigma} \approx i_{scattered}\times \rho_{target}\times l_{target}</math>
+
[[File:ComparingDCcomponents_S1_PhiThetaHits.png| 800 px]]
  
  
−
where ρ<sub>target</sub> is the density of the target material, l<sub>target</sub> is the length of the target, and i<sub>scattered</sub> is the number of incident particles scattered.
 
  
 +
[[File:ComparingDCEndplates_S1_PhiThetaHits.png | 800 px]]
  
−
:::::<math>\rho_{target}\times l_{target}=\frac{70.85 kg}{1 m^3}\times \frac{1 mole}{2.02 g} \times \frac{1000g}{1 kg} \times \frac{6\times10^{23} atoms}{1 mole} \times \frac{1m^3}{(100 cm)^3} \times \frac{1 cm}{ } \times \frac{10^{-24} cm^{2}}{barn} =2.10\times 10^{-2} barns</math>
 
  
 +
[[File:ComparingMagnetComponents_S1_PhiThetaHits.png| 800 px]]
  
−
:::::<math>\frac{1}{\rho_{target}\times l_{target} \times 4\times 10^7}=1.19\times 10^{-6} barns</math>
+
[[File:S1_PhiThetaGammaHits_Full.png | 600 px]]
  
  
−
[[File:MolThetaCM.png|frame|center|alt=Moller Electron Angle Theta in Center of Mass Frame|Figure 5a: The Moller electron scattering angle theta distribution for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.]]
 
  
  
 +
Tomography
  
−
[[File:Theory_new.png|frame|center|alt=Theoretical Moller Differential Cross-Section in Center of Mass Frame Frame|Figure 5b: The theoretical Moller electron differential cross-section for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.]]
+
[[File:S1_PhiThetaGammaVertex_wo_MagnetComponents.png | 600 px]]
  
−
Combining these plots, and rescaling the Final Theta in the Center of Mass for micro-barns, we find
 
  
 +
=GEANT4 Simulation of Lead Conic Moller Shield=
 +
==[[Lead Shield Cone]]==
  
 +
====OLD====
 +
{| class="wikitable"
 +
|+ FTOn from Forward Vertex
 +
! 50nA
 +
! S1R1 2ndryMoller e- rate
 +
! S1R1 2ndryMoller gamma rate
 +
! S1R1 2ndryMoller particle rate
 +
! Effective Shield Rate
 +
|-
 +
! R1_36_38_R2_36_38
 +
| ?
 +
| ?
 +
| ?
 +
| ?
  
−
<center>[[File:XSect_new_zoom.png]]</center>
+
|}
 +
{| class="wikitable"
 +
|+ OLD! FTOn from 0,0,0 Vertex
 +
! 50nA
 +
! S1R1 2ndryMoller e- rate
 +
! S1R1 2ndryMoller gamma rate
 +
! S1R1 2ndryMoller particle rate
 +
! Effective Shield Rate
 +
|-
 +
! R1_36_38_R2_36_38
 +
| 405 Hz
 +
| 15480 Hz
 +
| 160 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_50_52
 +
| 470 Hz
 +
| 15227 Hz
 +
| 146 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_70_72
 +
| 461 Hz
 +
| 15045 Hz
 +
| 150 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_75_77
 +
| 372 Hz
 +
| 14916 Hz
 +
| 130 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_80_82
 +
| 376 Hz
 +
| 14995 Hz
 +
| 109 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_90_92
 +
| 413 Hz
 +
| 14580 Hz
 +
| 119 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_95_97
 +
| 383 Hz
 +
| 14186 Hz
 +
| 117 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_111_113
 +
| 447 Hz
 +
| 14196 Hz
 +
| 109 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_116_118
 +
| 420 Hz
 +
| 14167 Hz
 +
| 144 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_121_123
 +
| 389 Hz
 +
| 14251 Hz
 +
| 117 Hz
 +
| ?
 +
|-
 +
! R1_74_76_R2_151_153
 +
| 492 Hz
 +
| 14280 Hz
 +
| 110 Hz
 +
| ?
 +
|-
 +
! R1_36_38_R2_500_503
 +
| 1000 Hz
 +
| 18363 Hz
 +
| 120 Hz
 +
| ?
 +
|}
  
−
=Step 2=
+
{| class="wikitable"
 +
|+ FTOn from (0,0,0) Vertex w/o FT
 +
! 50nA
 +
! S1R1 2ndryMoller e- rate
 +
! S1R1 2ndryMoller gamma rate
 +
! S1R1 2ndryMoller particle rate
 +
! Effective Shield Rate
 +
|-
 +
! R1_0_524.0_R2_0_1034.47
 +
| ?
 +
| ?
 +
| ?
 +
| ?
  
−
Replace the LH2 target with an NH3 target and compare with LH2 target.
+
|}
  
 +
{| class="wikitable"
 +
|+ FTOn from Forward Vertex
 +
! 50nA
 +
! S1R1 2ndryMoller e- rate
 +
! S1R1 2ndryMoller gamma rate
 +
! S1R1 2ndryMoller particle rate
 +
! Effective Shield Rate
 +
|-
 +
! R1_36_38_R2_36_38
 +
| ?
 +
| ?
 +
| ?
 +
| ?
  
−
[[File:FnlMomNH3.png]][[File:MolMomNH3.png]]
+
|}
−
 
 
−
The Moller Momentum plot for higher momenta run to increase the number of values
 
−
 
 
−
[[File:Full_MolMom.png]]
 
−
 
 
−
 
 
−
 
 
−
[[File:FnlMomNH3CM.png]][[File:MolMomNH3CM.png]]
 
−
 
 
−
 
 
−
[[File:FnlThetaNH3.png]][[File:MolThetaNH3.png]]
 
−
 
 
−
 
 
−
[[File:FnlThetaNH3CM.png]][[File:MolThetaNH3CM.png]]
 
−
 
 
−
 
 
−
 
 
−
 
 
−
==LH2 Vs. NH3==
 
−
 
 
−
 
 
−
:::::<math>L=\frac{870 kg}{1 m^3}\times \frac{1 mole}{17.03 g} \times \frac{1000g}{1 kg} \times \frac{6\times10^{23} atoms}{1 mole} \times \frac{1cm}{100 cm} \times \frac{1 m}{ } \times \frac{10^{-23} m^2}{barn} =2.10\times 10^{-2} barns</math>
 
−
 
 
−
[[File:DualFnlMom.png]][[File:DualMolMom.png]]
 
−
 
 
−
 
 
−
[[File:DualFnlTheta.png]][[File:DualMolTheta.png]]
 
−
 
 
−
===Figure out the offset===
 
−
 
 
−
Rerunning the GEANT simulation for a target of solo atoms of Carbon 12, with 4<sup>10</sup> incident 11GeV electrons
 
−
 
 
−
 
 
−
[[File:3setsMom.png]][[File:3targets.png]]
 
−
 
 
−
Density of target material:
 
−
 
 
−
C=2.26 g/cm<sup>3</sup>
 
−
 
 
−
NH3=.86g/cm<sup>3</sup>
 
−
 
 
−
LH2=.07g/cm<sup>3</sup>
 
−
 
 
−
Altering the density of LH2 to be .86g/cm<sup>3</sup> we find
 
−
 
 
−
[[File:MolThetaMod4.png]]
 
−
 
 
−
 
 
−
Comparing this to the theoretical differential cross section:
 
−
 
 
−
[[DV_RunGroupC_Moller#Comparing_experimental_vs._theoretical_for_M.C3.B8ller_differential_cross_section_11GeV |As shown above ]]
 
−
 
 
−
We find that the differential cross section scale is <math>\frac{d\sigma}{d\Omega}\approx 16.2\times 10^{-2}mb=16.2\mu b</math>
 
−
 
 
−
 
 
−
 
 
−
Converting the number of electrons to barns,
 
−
 
 
−
:::::<math>L=\frac{i_{scattered}}{\sigma} \approx i_{scattered}\times \rho_{target}\times l_{target}</math>
 
−
 
 
−
 
 
−
where ρ<sub>target</sub> is the density of the target material, l<sub>target</sub> is the length of the target, and i<sub>scattered</sub> is the number of incident particles scattered.
 
−
 
 
−
For LH2:
 
−
 
 
−
:::::<math>\frac{1}{\rho_{target}\times l_{target} \times 4\times 10^8}=1.19\times 10^{-8} barns=.0119 \mu b</math>
 
−
 
 
−
 
 
−
For Carbon:
 
−
:::::<math>\rho_{target}\times l_{target}=\frac{2.26 g}{1 cm^3}\times \frac{1 mole}{12.0107 g} \times \times \frac{6\times10^{23} atoms}{1 mole} \times \frac{1 cm}{ } \times \frac{10^{-24} cm^{2}}{barn} =1.13\times 10^{-1} barns</math>
 
−
 
 
−
 
 
−
:::::<math>\frac{1}{\rho_{target}\times l_{target} \times 4\times 10^8}=2.21\times 10^{-8} barns=.0221 \mu b</math>
 
−
 
 
−
For Ammonia:
 
−
 
 
−
 
 
−
:::::<math>\rho_{target}\times l_{target}=\frac{.8 g}{1 cm^3}\times \frac{1 mole}{17 g} \times  \frac{6\times10^{23} atoms}{1 mole} \times \frac{1 cm}{ } \times \frac{10^{-24} cm^2}{barn} =2.82\times 10^{-2} barns</math>
 
−
 
 
−
 
 
−
:::::<math>\frac{1}{\rho_{target}\times l_{target} \times 4\times 10^8}=8.87\times 10^{-8} barns=.0887 \mu b</math>
 
−
 
 
−
Combing plots in Root:
 
−
 
 
−
<pre>new TBrowser();
 
−
TH1F *LH2=new TH1F("LH2","LH2",360,90,180);
 
−
LH2->Add(MollerThetaCM,1.19e-8);
 
−
LH2->Draw();
 
−
 
 
−
TH1F *C12=new TH1F("C12","C12",360,90,180);
 
−
C12->Add(MollerThetaCM,2.21e-8);
 
−
C12->Draw();
 
−
 
 
−
TH1F *NH3=new TH1F("NH3","NH3",360,90,180);
 
−
NH3->Add(MollerThetaCM,8.87e-8);
 
−
NH3->Draw();
 
−
LH2->Draw("same");
 
−
C12->Draw("same");
 
−
Theory->Draw("same");</pre>
 
−
 
 
−
 
 
−
[[File:Combo4.png]]
 
−
 
 
−
=Step 3=
 
−
 
 
−
Determine impact of Solenoid magnet on Moller events
 
−
 
 
−
=Papers used=
 
−
 
 
−
[1]Farrukh Azfar's Derivation of Moller Scattering
 
−
 
 
−
:::[[File:FarrukAzfarMollerScatter.pdf]]
 
−
 
 
−
A polarized target for the CLAS detector
 
−
 
 
−
:::[[File:PHY02-33.pdf‎]]
 
−
 
 
−
An investigation of the spin structure of the proton in deep inelastic scattering of polarized muons on polarized protons
 
−
 
 
−
:::[[:File:1819.pdf]]
 
−
 
 
−
==QED Radiative Corrections to Low-Energy Moller and Bhabha Scattering==
 
−
http://arxiv.org/abs/1602.07609
 
  
−
=Moller Track Reconstruction=
 
  
−
[[DV_MollerTrackRecon]]
+
[[VanWasshenova_Thesis]]
  
 
[[EG12]]
 
[[EG12]]

Latest revision as of 18:03, 13 January 2019

need to insert moller shielding into card after moller LUND file is created. (see clas12/beamline)

Simulating the Moller scattering background for EG12

Docker

Set up Docker Container

Moller Differential Cross-Section

GEANT4 Simulation of Moller Events

Simulation Setup

Use GEANT4 via GEMC to estimate the Moller background for electron scattering experiments in JLab's Hall B. The first step towards this goal is to use GEANT4 without the GEMC infrastructure to create event (LUND) files that will be used as input events for GEMC.

Benchmark GEANT4's Moller scattering prediction with the theoretical cross section using LH2

Comparison of simulation vs. the theoretical Møller differential cross section using 11 GeV electrons impinging LH2

Converting the number of scattered electrons per scattering angle theta to a differential cross-section in barns.

Experimental and Theoretical Moller Differential Cross-Section in Center of Mass Frame Frame
Figure 5c: The experimental and theoretical Moller electron differential cross-section for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference.

Change to a NH3 Target

Replacing the LH2 target with an NH3 target

Benchmark GEANT4's Moller scattering prediction with the theoretical cross section using NH3

Comparison of simulation vs. the theoretical Møller differential cross section using 11 GeV electrons impinging NH3

Converting the number of scattered electrons per scattering angle theta to a differential cross-section in barns.

Theoretical and Simulated Moller Differential Cross-Section in Center of Mass Frame Frame
Figure 5c: The theoretical and simulated Moller electron differential cross-section for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference for NH3 target.

LH2 Vs. NH3

Benchmark GEANT4's Moller scattering prediction for NH3 and LH2

Effects Due to Target Material

Target Density

Atomic Mass and Electron Number Effects

Differential Cross-Section Offset

Theoretical and Simulated Moller Differential Cross-Section in Center of Mass Frame Frame
Figure 5c: The theoretical and simulated Moller electron differential cross-section for an incident 11 GeV(Lab) electron in the Center of Mass frame of reference for NH3 target. The theoretical differential cross-section has been adjusted to account for the 10 electrons within NH3.

ROOT Macro to read LUND files and make plots


LUND MolMomCM.pngLUND MolThetaCM.png


LUND MolMomLab.pngLUND MolThetaLab.png


Run the above theta and E distribution of Mollers through GEMC with and without Solenoid on.  Determine the Theta and E range of Mollers that enter the detector.


Scanning from 2MeV to 5500MeV, and 0 to 36 degrees in Theta.

For Solenoid on with 5T: 

<gcard>

        <detector name="../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/dc/dc"            factory="TEXT" variation="original"/>
        <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="original"/>
        <detector name="../../../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/pcal/pcal"            factory="TEXT" variation="original"/>
        <option name="BEAM_P"   value="e-, 2.800*GeV, 18.0*deg, 10*deg"/>
        <option name="SPREAD_P" value="2.798*GeV, 18*deg, 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>

Using the standard commands for gemc 

gemc -USE_GUI=0 -N=1000 eg12_sol.gcard
~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -i eg12.ev -config DCHB::torus=-1.0 -config DCHB::solenoid=1.0 -config DCTB::kalman=true -o eg12_rec.ev -s DCHB:DCTB:EC:FTOF:EB
~/src/CLAS/coatjava-1.0/bin/run-groovy Analysis.groovy eg12_rec.0.evio


Checking for a reconstructed particle that undergoes a phi shift where the 1st column is energy and the 2nd the scattering angle theta: 

For Energy:
awk 'NR == 1 {line =$0; min =$2} NR >1 && $2 < min {line =$0; min = $2} END{print line}' Energy_Phi_Shift.dat
289     0.6624017631112946      33.37507824033966       -179.63657653407975     0.7254709       29.696417       -177.79057      473.36051689004984      11.978802748439216      -129.55774724876173


awk 'NR == 1 {line =$0; max =$2} NR >1 && $2 > max {line =$0; max = $2} END{print line}' Energy_Phi_Shift.dat
978     5.584765946130235       24.90893565820338       -170.69784116534063     5.593401        24.90894        -170.65994      472.16552666546283      23.31807917059094       -165.11125601370114

For Theta:
awk 'NR == 1 {line =$0; min =$3} NR >1 && $3 < min {line =$0; min = $3} END{print line}' Energy_Phi_Shift.dat
12      4.83485521395583        10.577021312967325      -50.86708775351051      4.8205996       10.617784       -50.895683      491.627055439108        5.9032309421001985      -43.82964507576243

awk 'NR == 1 {line =$0; max =$3} NR >1 && $3 > max {line =$0; max = $3} END{print line}' Energy_Phi_Shift.dat
187     1.3382869250001856      35.63598792463312       29.887530073265214      1.3481133       35.449417       30.009071       472.455052447561        19.034043420497742      55.20328358414997


Similarly for 0T: 

For Energy:
awk 'NR == 1 {line =$0; min =$2} NR >1 && $2 < min {line =$0; min = $2} END{print line}' Energy_Phi_Shift.dat
148     0.36418038435234207     33.68296595663359       -175.89942431595827     0.35383263      51.743824       -39.45801       469.5694572319308       24.353383233947522      -176.27025218702502


awk 'NR == 1 {line =$0; max =$2} NR >1 && $2 > max {line =$0; max = $2} END{print line}' Energy_Phi_Shift.dat
818     5.5278794590742235      12.937136966282836      -43.35663621477178      5.539351        12.95216        -49.09295       484.8800252169308       8.27731670021694        -43.245753766203364

For Theta:
awk 'NR == 1 {line =$0; min =$3} NR >1 && $3 < min {line =$0; min = $3} END{print line}' Energy_Phi_Shift.dat
559     4.909449160222199       9.918634647646794       0.06548256930681198     4.914272        9.938521        -5.861595       492.1448533455858       7.501204563534799       -0.006721611135017764

awk 'NR == 1 {line =$0; max =$3} NR >1 && $3 > max {line =$0; max = $3} END{print line}' Energy_Phi_Shift.dat
326     1.272966353429713       35.731680814303196      85.3173707349696        1.2830343       38.27654        58.368687       483.4835073848272       3.950855595379657       83.87603767806992

Reconstruction of Moller Events

The LUND format

Writing LUND files

Running GEMC

Phi Shift observation using DC hit Reconstruction Data

Gcard creation defining energy and angle range of electrons

Modified gcards

Effects of Solenoid on Electrons

Results for defined distribution in Solenoid Fields

Results for Random Spread of Energy and angle theta in the Lab frame

Cover Full Solid Angle of Detector

Using GEANT simulation data

Calculations of 4-momentum components

Alter Phi Angles

Check Differential Cross-Section

Detector Occupancy

clas12->Draw("Detector.wire:Detector.layer>>(7,1,7,120,0,120)","Detector.superlayer<3","colz")

Gcard settings

Verfication of Mother/Daughter Occupancy

Benchmark GEMC Occupancy Prediction for 11GeV Electron Beam with 0T Solenoid for LH2

Setup for Batch Job With Varying Experimental CLAS12 Quantities

Run_in_GEMC

Center_of_Mass_for_Stationary_Target

Run Occupancy for Sector 1 DC hits

Wire_angle_correspondance

Isotropic Weighted Moller Distribution in GEMC

Papers used

[1]Farrukh Azfar's Derivation of Moller Scattering

File:FarrukAzfarMollerScatter.pdf

A polarized target for the CLAS detector

File:PHY02-33.pdf

An investigation of the spin structure of the proton in deep inelastic scattering of polarized muons on polarized protons

File:1819.pdf

QED Radiative Corrections to Low-Energy Moller and Bhabha Scattering

http://arxiv.org/abs/1602.07609

ChangingRates S1 PhiThetaHits.png

S1 50nA PrimaryElectronSigmasWeightedRates.png

S1 PhiThetaGammaHits.png


HitMakeUp.png


S1 50nA PrimaryElectronSigmasWeightedRates Full.png

ComparingOppositeFields S1 PhiThetaHits.png


ChangingSolenoidRates wo Magnets.png


ComparingDCcomponents S1 PhiThetaHits.png


ComparingDCEndplates S1 PhiThetaHits.png


ComparingMagnetComponents S1 PhiThetaHits.png

S1 PhiThetaGammaHits Full.png



Tomography

S1 PhiThetaGammaVertex wo MagnetComponents.png


GEANT4 Simulation of Lead Conic Moller Shield

Lead Shield Cone

OLD

FTOn from Forward Vertex
50nA S1R1 2ndryMoller e- rate S1R1 2ndryMoller gamma rate S1R1 2ndryMoller particle rate Effective Shield Rate
R1_36_38_R2_36_38 ? ? ? ?
OLD! FTOn from 0,0,0 Vertex
50nA S1R1 2ndryMoller e- rate S1R1 2ndryMoller gamma rate S1R1 2ndryMoller particle rate Effective Shield Rate
R1_36_38_R2_36_38 405 Hz 15480 Hz 160 Hz ?
R1_36_38_R2_50_52 470 Hz 15227 Hz 146 Hz ?
R1_36_38_R2_70_72 461 Hz 15045 Hz 150 Hz ?
R1_36_38_R2_75_77 372 Hz 14916 Hz 130 Hz ?
R1_36_38_R2_80_82 376 Hz 14995 Hz 109 Hz ?
R1_36_38_R2_90_92 413 Hz 14580 Hz 119 Hz ?
R1_36_38_R2_95_97 383 Hz 14186 Hz 117 Hz ?
R1_36_38_R2_111_113 447 Hz 14196 Hz 109 Hz ?
R1_36_38_R2_116_118 420 Hz 14167 Hz 144 Hz ?
R1_36_38_R2_121_123 389 Hz 14251 Hz 117 Hz ?
R1_74_76_R2_151_153 492 Hz 14280 Hz 110 Hz ?
R1_36_38_R2_500_503 1000 Hz 18363 Hz 120 Hz ?
FTOn from (0,0,0) Vertex w/o FT
50nA S1R1 2ndryMoller e- rate S1R1 2ndryMoller gamma rate S1R1 2ndryMoller particle rate Effective Shield Rate
R1_0_524.0_R2_0_1034.47 ? ? ? ?
FTOn from Forward Vertex
50nA S1R1 2ndryMoller e- rate S1R1 2ndryMoller gamma rate S1R1 2ndryMoller particle rate Effective Shield Rate
R1_36_38_R2_36_38 ? ? ? ?


VanWasshenova_Thesis

EG12

Retrieved from "https://wiki.iac.isu.edu/index.php?title=DV_RunGroupC_Moller&oldid=126004"