[G12] G12 efficiency script

[Michael C. Kunkel]

Greetings,

To clarify. I am not placing this script into g12 svn until I am 
finished with the multiple photon efficiency.

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/6/15 4:40 PM, Michael C. Kunkel wrote:
> 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
>>>>
>>>
>>
>
>
>
> _______________________________________________
> G12 mailing list
> G12 at jlab.org
> https://mailman.jlab.org/mailman/listinfo/g12

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