Difference between revisions of "HRRL Pos Test 27Jul2012 Data Ana 3 MeV"

From New IAC Wiki
Jump to navigation Jump to search
Line 1: Line 1:
−
= 3 MeV =
 
−
== 2 ways to calculate electron charge ==
 
−
   
 
−
    NE(Int_t ri, Int_t Reprate, Int_t calc_meth){
 
−
    printf("%d %d \n",ri, Reprate);
 
−
    if(ri==1) TFile *f = new TFile("C:/Users/sadiq/Documents/phd/Runs_Jul_2012/rootfiles/7-27-12/r3735.root");
 
−
    else if(ri==2) TFile *f = new TFile("C:/Users/sadiq/Documents/phd/Runs_Jul_2012/rootfiles/7-27-12/r3737.root");
 
−
    TTree *ntuple = (TTree*)f->Get("ntuple");
 
−
    TH1F *f9=new TH1F("f9","ADC9",4000,0,4000);
 
−
        TCanvas * c2 = new TCanvas("c2", "c2",600 ,500);
 
−
        c2->cd();
 
−
        Long64_t ne;
 
−
        printf("%d \n", ne);
 
−
      //method 1 is reprate times single pulse charge
 
  
−
=== method 1===
+
 
−
+
= Hand calculation on electron rate =
−
        Double_t Qp1=0,Qp2=0;
+
 
−
        if(calc_meth == 1){
+
At full current on Faraday cup 1 we have 1201 nVs/ADC channel.
−
        //ntuple->Draw("ADC9*(1.22*1e-9)/50>>f9","1");
+
 
−
        ntuple->Draw("ADC9>>f9");
+
<math>\frac{1201~nVs}{50 \Omega} = 24~nAs = 24~nC</math>
−
        //Qp2 = (f9->Integral())*(1.22*1e-9)/50; //charge calculate with method 2.
+
 
−
        //f9->Fit("gaus");
+
charge rate under 300 Reprate: <math>300 \times 24~nAs = 7200~nC = 7.2 \frac{\mu C }{ s}</math>
−
        ADC9_fit_mean = f9->GetMean(1);
+
 
−
        //ADC9_fit_mean = ADC9->GetMean(1);
+
 
−
        printf("ADC9_fit_mean = %d \n",ADC9_fit_mean);
+
<math> \frac{7.2 \times 10^{-6}\frac{ C }{ s}} {1.6 \times 10^{-19} C} = 4.5 \times 10^{13}</math>
−
        Qp1 = ADC9_fit_mean*(1.22*1e-9)/50; //Qp: charge per pulse. The calibration: 1.22 (9nV s)/(ADC channel), 50 ohm Terminal.
+
 
−
        printf("QP1=%g\n\n",Qp1);
 
−
        printf("e- charge per pulse with method 1: %g\n",Qp1);
 
−
        //printf("Reprate: %d\n",Reprate);
 
−
        Qs=Reprate*Qp1; //Qs: Charge per second.
 
−
        //printf("e- charge per second: %g\n",Qs);
 
−
        ne=Qs/Q_e; //Ne: number of electrons per second.
 
−
                }
 
  
  
−
=== method 2===
 
−
 
−
        if(calc_meth == 2){
 
−
        //ntuple->Draw("ADC9*(1.22*1e-9)/50>>f9","1");
 
−
        ntuple->Draw("ADC9>>f9");
 
−
        printf("%bin# Bin Content%d\n");
 
−
        for(i=1; i<=4000; i++){
 
−
        //for(i=0; i<=600; i++){
 
−
        bin_content[i] = f9->GetBinContent(i);
 
−
        //Qp2+=(i*(1.22*1e-9)/50)*bin_content[i];
 
−
        Qp2+=i*bin_content[i]*(1.22*1e-9)/50;
 
−
        //printf("%d  %d  %g \n",i,bin_content[i],Qp2);
 
−
        }
 
−
 
−
        //printf("Integral = %d \n",f9->Integral());
 
−
        //Qp2 = (f9->Integral())*(1.22*1e-9)/50; //charge calculate with method 2.   
 
−
        printf("QP2=%g\n\n",Qp2);
 
−
        printf("total e- charge in this run with method 2: %g\n",Qp2);
 
−
        //printf("Reprate: %d\n",Reprate);
 
−
        Qs=Qp2; //Qs: Charge per second.
 
−
        printf("e- charge per second: %g\n",Qs);
 
−
        printf("e- charge: %g\n",Q_e);
 
−
        ne=Qs/Q_e; //Ne: number of electrons per second.
 
−
        printf("1. total number of electrons in this run: %g\n",ne);
 
−
                }
 
−
        printf("2. total number of electrons in this run: %g\n",ne);
 
−
        c2->Close();
 
−
        f->Close();
 
−
        return(ne);
 
−
        }
 
  
−
== results in 3 MeV ==
+
= Data =
  
−
=== results on 3 MeV ===
 
  
−
{| border="1"
 
 
|-
 
|-
−
|run in||run out || NaI Left: in sub out || NaI Right: in sub out|| NaI Left: overlap in and out || NaI Right: overlap in and out ||e- rate (Hz)||e+ rate (Hz)||e+/e+ ratop || e- current (A)
+
|3735||3736||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_sub_r3736_DL.png|100 px]]
 +
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_sub_r3736_DR.png|100 px]]
 +
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_ove_r3736_DL.png|100 px]]
 +
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_ove_r3736_DR.png|100 px]]
 
|-
 
|-
 +
|}
 +
 +
Reprate = 300 Hz
 +
 +
Run time = 1002 s
 +
 +
Pulses = 301462
  
−
|3735||3736
+
Events = 9045
  
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_sub_r3736_DL_1.png|100 px]]
+
e- rate = 7.80042e+013 +- 2.25903e+012 (Hz)
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_sub_r3736_DR_1.png|100 px]]
 
  
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_over_r3736_DL_1.png|100 px]]
+
e+ rate = 1.41748 +- 0.0376119 (Hz).
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3735_over_r3736_DR_1.png|100 px]]
 
  
−
|| 1.17925e+012||0.255007||2.16244e-013||1.88937e-007
+
e+/e-=1.81719e-014+-4.8238e-016.
  
 +
e- current:I_e- =6.26271e-006+-3.61937e-007 (A).
 +
 +
e+/I_e-=1.4 Hz e+ / 6.3 \muA e-
 
|-
 
|-
 +
|3737||3736||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_sub_r3736_DL.png|100 px]]
 +
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_sub_r3736_DR.png|100 px]]
 +
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_ove_r3736_DL.png|100 px]]
 +
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_ove_r3736_DR.png|100 px]]
 +
|-
 +
|}
 +
 +
Reprate = 300 Hz
 +
 +
Run time = 1094 s
 +
 +
Pulses = 329368
 +
 +
Events = 15361
 +
 +
e- rate = 1.36517e+014 +- 2.39286e+012 (Hz)
  
−
|3737||3736
+
e+ rate = 1.60462 +- 0.0382982 (Hz).
  
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_sub_r3736_DL_1.png|100 px]]
+
e+/e-=1.1754e-014+-2.8058e-016.
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_sub_r3736_DR_1.png|100 px]]
+
 
 +
e- current:I_e- =1.52232e-005+-3.83378e-007 (A).
 +
 
 +
e+/I_e-=1.6 Hz e+ / 15 \muA e-
  
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_over_r3736_DL_1.png|100 px]]
 
−
||[[File:hrrl_pos_27jul2012_data_ana_with_Cuts_r3737_over_r3736_DR_1.png|100 px]]
 
  
−
||1.97625e+012||0.251467||1.27244e-013||3.1663e-007
 
  
−
|-
 
−
|}
 
  
−
=== Hand calculation on electron rate ===
 
  
−
At full current on Faraday cup 1 we have 1201 nVs/ADC channel.
 
  
−
<math>\frac{1201~nVs}{50 \Omega} = 24~nAs = 24~nC</math>
 
  
−
charge rate under 300 Reprate: <math>300 \times 24~nAs = 7200~nC = 7.2 \frac{\mu C }{ s}</math>
 
  
  
−
<math> \frac{7.2 \times 10^{-6}\frac{ C }{ s}} {1.6 \times 10^{-19} C} = 4.5 \times 10^{13}</math>
 
  
  
  
 
[[https://wiki.iac.isu.edu/index.php/Beam_Test_Jul_27_2012 go back]]
 
[[https://wiki.iac.isu.edu/index.php/Beam_Test_Jul_27_2012 go back]]

Revision as of 05:39, 20 September 2012


Hand calculation on electron rate

At full current on Faraday cup 1 we have 1201 nVs/ADC channel.

[math]\frac{1201~nVs}{50 \Omega} = 24~nAs = 24~nC[/math]

charge rate under 300 Reprate: [math]300 \times 24~nAs = 7200~nC = 7.2 \frac{\mu C }{ s}[/math]


[math] \frac{7.2 \times 10^{-6}\frac{ C }{ s}} {1.6 \times 10^{-19} C} = 4.5 \times 10^{13}[/math]



Data

|- |3735||3736||Hrrl pos 27jul2012 data ana with Cuts r3735 sub r3736 DL.png ||Hrrl pos 27jul2012 data ana with Cuts r3735 sub r3736 DR.png ||100 px ||100 px |- |}

Reprate = 300 Hz

Run time = 1002 s

Pulses = 301462

Events = 9045

e- rate = 7.80042e+013 +- 2.25903e+012 (Hz)

e+ rate = 1.41748 +- 0.0376119 (Hz).

e+/e-=1.81719e-014+-4.8238e-016.

e- current:I_e- =6.26271e-006+-3.61937e-007 (A).

e+/I_e-=1.4 Hz e+ / 6.3 \muA e- |- |3737||3736||Hrrl pos 27jul2012 data ana with Cuts r3737 sub r3736 DL.png ||Hrrl pos 27jul2012 data ana with Cuts r3737 sub r3736 DR.png ||100 px ||100 px |- |}

Reprate = 300 Hz

Run time = 1094 s

Pulses = 329368

Events = 15361

e- rate = 1.36517e+014 +- 2.39286e+012 (Hz)

e+ rate = 1.60462 +- 0.0382982 (Hz).

e+/e-=1.1754e-014+-2.8058e-016.

e- current:I_e- =1.52232e-005+-3.83378e-007 (A).

e+/I_e-=1.6 Hz e+ / 15 \muA e-







[go back]

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