[G12] G12 efficiency script

[Michael C. Kunkel]

Greetings,

I have written a script for this. This script is in beta, I will add 
more to it later, but the output is my efficiency correction for single 
photon selection.

Please look at my svn
https://jlabsvn.jlab.org/svnroot/clas/users/mkunkel/clas/g12_corrections/g12_trackEfficiency.hpp

you will need the histograms located in
https://jlabsvn.jlab.org/svnroot/clas/users/mkunkel/clas/g12_corrections/EFFICIENCY_HISTS/

To use,
     double tot_eff = clas::g12::g12_trackEfficiency(Vz, P_Ptot_fit, 
P_Phi_fit, P_Theta_fit, Ep_Ptot_fit, Ep_Phi_fit, Ep_Theta_fit, 
Em_Ptot_fit, Em_Phi_fit, Em_Theta_fit, TString DEG_or_RAD ,TString dir);

I have only tested this with lab angles that are in degrees, but there 
is a conversion in the script, not tested though.

I use it as
     Double_t xSection_scale = 
(1.e-30*Flux_norm*density_of_target*dTheta/tot_eff);

Then I weight my event by this number.

BR
MK
----------------------------------------
Michael C. Kunkel, PhD
Forschungszentrum Jülich
Nuclear Physics Institute and Juelich Center for Hadron Physics
Experimental Hadron Structure (IKP-1)
www.fz-juelich.de/ikp

On 7/2/15 4:41 PM, Akbar, Zulkaida wrote:
>
> Sure, big thanks MK.
>
>
> Best Regard,
>
>
> Zulkaida Akbar
>
>
>
> ------------------------------------------------------------------------
> *From:* Michael C. Kunkel <mkunkel at jlab.org>
> *Sent:* Thursday, July 2, 2015 8:39 PM
> *To:* Akbar, Zulkaida
> *Subject:* Re: omega cross sections:
> I do not have it in script form.
> What I can do is get a package together that you can mimic, but you 
> will still have to do some coding on your own.
> BR
> MK
> ----------------------------------------
> Michael C. Kunkel, PhD
> Forschungszentrum Jülich
> Nuclear Physics Institute and Juelich Center for Hadron Physics
> Experimental Hadron Structure (IKP-1)
> www.fz-juelich.de/ikp
> On 7/2/15 4:37 PM, Akbar, Zulkaida wrote:
>>
>> Hi MK,
>>
>>
>> I select single photon for my event.
>>
>>
>>
>> Best Regard,
>>
>>
>> Zulkaida Akbar
>>
>>
>>
>> ------------------------------------------------------------------------
>> *From:* Michael C. Kunkel <mkunkel at jlab.org>
>> *Sent:* Thursday, July 2, 2015 8:36 PM
>> *To:* Akbar, Zulkaida; Lei Guo
>> *Cc:* Johann Goetz
>> *Subject:* Re: omega cross sections:
>> This depends, do you use all "in time" photons in the event, or do 
>> you select single photons?
>> BR
>> MK
>> ----------------------------------------
>> Michael C. Kunkel, PhD
>> Forschungszentrum Jülich
>> Nuclear Physics Institute and Juelich Center for Hadron Physics
>> Experimental Hadron Structure (IKP-1)
>> www.fz-juelich.de/ikp
>> On 7/2/15 4:33 PM, Akbar, Zulkaida wrote:
>>>
>>> Dear Lei and MK,
>>>
>>>
>>> I am doing trigger inefficiency studies now, but it is very good 
>>> idea to use MK inefficiency map.  Where can I get the trigger 
>>> inefficiency map/Routine? Thanks a Lot.
>>>
>>>
>>>
>>> Best Regard
>>>
>>>
>>> Zulkaida Akbar
>>>
>>>
>>>
>>> ------------------------------------------------------------------------
>>> *From:* Lei Guo <lguo at jlab.org>
>>> *Sent:* Thursday, July 2, 2015 8:00 PM
>>> *To:* Akbar, Zulkaida
>>> *Cc:* Michael C. Kunkel; Johann Goetz
>>> *Subject:* omega cross sections:
>>> HI, Zulkaida,
>>>
>>> MK, Johann and I had a meeting today going over the g12 review 
>>> comments, specifically related to this comments:
>>>
>>> "When you compare cross sections to previous measurements please
>>> compare in the same kinematical bins. A strong forward peak pushes 
>>> other
>>> point to the axis and it is not possible to judge if they agree or 
>>> not. In addition
>>> to overlaying two data sets, show their ratio. Even without it we 
>>> see that in
>>> some bins the agreement is not very good. Make something similar to 
>>> pull
>>> distributions as well, to see if there is a bias."
>>>
>>> Could you do the following two
>>> 1. Apply MK’s track-dependent correction for the normalization. It 
>>> should be similar to your overall 15% (or is it 18% correction). MK 
>>> has some routine you can just call. If you want to do it quickly, at 
>>> least your should apply separate corrections (globally) for p pi+ pi-
>>> 2. Plot the ratio of your xsections over g11 and/or pull 
>>> distributions. If you have questions about details, please contact 
>>> MK/JOhann. We are looking to submit the response next week, so 
>>> hopefully we are not ruining your long weekend.
>>>
>>> Lei
>>>
>>>
>>>
>>> Lei Guo
>>> Assistant Professor
>>> Physics Department
>>> Florida International University
>>> Miami, FL
>>>
>>> email: leguo at fiu.edu <mailto:leguo at fiu.edu> or lguo at jlab.org 
>>> <mailto:lguo at jlab.org>
>>> Office:305-348-0234
>>>
>>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.jlab.org/pipermail/g12/attachments/20150706/d5b76f6d/attachment-0001.html>