Hello, 'residual' kinematic corrections here.
I did not parameterize them as a polynomial because I
didn't like the extreme variation on the ends. If anyone
would like to parameterize the points themselves,
set up some TGraphs and see if you can get better looking
pol8, pol9 TFunctions
Just initialize the offsets for your version and:
TGraph* eLeft = new TGraph(20,eTheta,eL);
eLeft->Draw("AP*");
The algorithm for corrections...
float eCorrectedMomentum(int sector, float Pe, float THe) {
if (!sector) { //assume zero == eLeft
int Bin = (THe-24.0)/2.7; // electron theta binning (good from 24-78)
if (Bin>=0 && Bin<20) Pe *= eL[Bin]; // Multiplicative offset
else if (Bin<0) Pe *= eL[0];
else Pe *= eL[19];
}
else if (sector) {
int Bin = (THe-24.0)/2.7; // electron theta binning
if (Bin>=0 && Bin<20) Pe *= eR[Bin]; // Multiplicative offset
else if (Bin<0) Pe *= eR[0];
else Pe *= eR[19];
}
return Pe;
}
float pCorrectedMomentum(int sector, float Pp, float THp) {
if (!sector) { //assume zero == eLeft
int Bin = (THp-24.0)/2.7; // proton theta binning (good from 33-69)
if (Bin>=0 && Bin<20) Pp *= pR[Bin]; // Multiplicative offset
else if (Bin<0) Pp *= pR[0];
else Pp *= pR[19];
}
else if (sector) {
int Bin = (THp-33.0)/1.8; // proton theta binning (good from 33-69)
if (Bin>=0 && Bin<20) Pp *= pL[Bin]; // Multiplicative offset
else if (Bin<0) Pp *= pL[0];
else Pp *= pL[19];
}
return Pp;
}
//**///////////////////////////////////////////////**
//**//// V17 ELoss subtracted ////**
//**///////////////////////////////////////////////**
//Multiplicative offsets
//**************************
//Electron Theta Binning
Float_t eTheta[20] =
{25.350000,28.050000,30.750000,33.450000,36.150000,38.850000,41.550000,44.250000,46.95000
0,49.650000,52.350001,55.050001,57.750001,60.450001,63.150001,65.850001,68.550001,71.2500
01,73.950001,76.650000};
//Left Sector Electron
//--------------------
Float_t eL[20] =
{1.015636,1.017280,1.020067,1.022152,1.024380,1.024887,1.022006,1.028023,1.025257,1.031380
,1.029298,1.027893,1.027475,1.024166,1.030168,1.029777,1.023533,1.029012,1.021426,1.022334
};
//Right Sector Electron
//---------------------
Float_t eR[20] =
{1.011404,1.022541,1.029819,1.028645,1.018602,1.017281,1.018398,1.021183,1.021764,1.019833
,1.019112,1.016707,1.013193,1.011465,1.011629,1.009358,1.010387,1.014550,1.018511,1.020436
};
//Proton Theta Binning
Float_t pTheta[20] =
{33.900000,35.700000,37.500000,39.300000,41.100000,42.900000,44.700000,46.500000,48.30000
0,50.100000,51.899999,53.699999,55.499999,57.299999,59.099999,60.899999,62.699999,64.4999
99,66.299999,68.100000};
//Left Sector Proton
//--------------------
Float_t pL[20] =
{1.001029,0.996832,0.993442,0.992234,0.996468,0.998291,0.995417,0.992569,0.995615,0.997916
,0.998761,0.998967,0.997358,1.003048,1.005920,1.009334,1.013055,1.014365,1.016129,0.998989
};
//Right Sector Proton
//---------------------
Float_t pR[20] =
{1.001661,1.007988,1.004757,1.001283,0.997066,0.996414,0.997819,0.996660,0.999607,1.005860
,1.010342,1.008562,1.007439,1.012586,1.012559,1.012973,1.012452,1.013697,1.017158,0.982159
};
//**///////////////////////////////////////////////**
//**//// V18 ELoss subtracted ////**
//**///////////////////////////////////////////////**
//Multiplicative offsets
//**************************
Electron Theta Binning
Float_t eTheta[20] =
{25.350000,28.050000,30.750000,33.450000,36.150000,38.850000,41.550000,44.250000,46.95000
0,49.650000,52.350001,55.050001,57.750001,60.450001,63.150001,65.850001,68.550001,71.2500
01,73.950001,76.650000};
//Left Sector Electron
//--------------------
Float_t eL[20] =
{1.016035,1.018352,1.020069,1.021894,1.024036,1.024766,1.021387,1.027320,1.024610,1.030527
,1.028730,1.027321,1.026770,1.023156,1.029372,1.027028,1.023978,1.027151,1.020025,1.020711
};
//Right Sector Electron
//---------------------
Float_t eR[20] =
{1.012006,1.022475,1.029579,1.023681,1.011731,1.016225,1.018426,1.022030,1.021525,1.019610
,1.018735,1.016122,1.012492,1.010626,1.010086,1.007577,1.002128,1.009946,1.017456,1.023062
};
//Proton Theta Binning
Float_t pTheta[20] =
{33.900000,35.700000,37.500000,39.300000,41.100000,42.900000,44.700000,46.500000,48.30000
0,50.100000,51.899999,53.699999,55.499999,57.299999,59.099999,60.899999,62.699999,64.4999
99,66.299999,68.100000};
//Left Sector Proton
//--------------------
Float_t pL[20] =
{0.999504,0.997195,0.993397,0.991293,0.996002,0.998386,0.995234,0.992361,0.995301,0.997454
,0.998314,0.998575,0.997084,1.001981,1.004527,1.006676,1.012316,1.013670,1.017440,0.999979
};
//Right Sector Proton
//---------------------
Float_t pR[20] =
{1.004651,1.007496,1.005751,1.005151,0.997538,0.996472,0.997084,0.995883,0.999052,1.005195
,1.009478,1.007644,1.006487,1.012024,1.012160,1.013766,1.012421,1.014731,1.018490,0.976884
};
Quoting Michael Kohl <kohlm@mit.edu>:
> Sorry,
>
> forgot to mention: The hydrogen v3_4_18 recrunch results were written
> to
> /net/data/5/Analysis/data/v3_4_18/lrd
>
> Regards,
>
> Michael
>
>
> On Sun, 2 Apr 2006, Yuan Xiao wrote:
>
> > Hi, Michael,
> >
> > Where are the new files?
> >
> > Yuan
> >
> > Quoting Michael Kohl <kohlm@MIT.EDU>:
> >
> >> Hi,
> >>
> >> the hydrogen v3_4_18 recrunch is complete.
> >> Major difference to v3_4_17 is that the holding field has been turned
> on
> >> for reconstruction, and the bugfix in pathlength.
> >>
> >> Eugene, please evaluate kinematic corrections for v3_4_18 (and
> v3_4_17!)
> >> while separately accounting for proton eloss and send the
> parametrizations
> >> to blast_anaware. How much has the coplanarity peak changed compared
> to
> >> v3_4_17?
> >>
> >> Yuan, please look at the pathlength variables, both the original ones
> and
> >> the recently added ones.
> >>
> >> Regards,
> >>
> >> Michael
> >>
> >>
> >> On Sat, 1 Apr 2006, Michael Kohl wrote:
> >>
> >>> Hi,
> >>>
> >>> I have compiled, tagged and linked v3_4_18
> libBlast.so/lrn/lrd/lrs/nsed.
> >>> I started a recrunch of the hydrogen data 12144-13278 with
> lrd_v3.4.18,
> >>> expect completion by Saturday 4/1/2006 at night. I'll notify.
> >>>
> >>> Changes from v3_4_17 to v3_4_18 are:
> >>> + Bugfix for pathlength
> >>> + Holding field turned on in tracking + Chris' extension of geometry
> class
> >>> to introduce shifted/rotated
> >>> wire chambers; code to calculate derivatives of p,th,ph,z with
> respect
> >>> to shifts and rotations
> >>>
> >>>
> >>> The holding field is expected to affect the coplanarity offset.
> >>>
> >>> Eugene, please evaluate kinematic corrections for v3_4_18 (and
> v3_4_17!)
> >>> and send the parametrizations to blast_anaware.
> >>>
> >>> Yuan, please look at the pathlength variables, both the original
> ones and
> >>> the recently added ones.
> >>>
> >>> Regards,
> >>>
> >>> Michael
> >>>
> >>>
> >>>
> >>> +-------------------------------------+--------------------------+
> >>> | Office: | Home: |
> >>> |-------------------------------------|--------------------------|
> >>> | Dr. Michael Kohl | Michael Kohl |
> >>> | Laboratory for Nuclear Science | 5 Ibbetson Street |
> >>> | MIT-Bates Linear Accelerator Center | Somerville, MA 02143 |
> >>> | Middleton, MA 01949 | U.S.A. |
> >>> | U.S.A. | |
> >>> | - - - - - - - - - - - - | - - - - - - - - -|
> >>> | Email: kohlm@mit.edu | K.Michael.Kohl@gmx.de |
> >>> | Work: +1-617-253-9207 | Home: +1-617-629-3147 |
> >>> | Fax: +1-617-253-9599 | Mobile: +1-978-580-4190 |
> >>> | http://blast.lns.mit.edu | |
> >>> +-------------------------------------+--------------------------+
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>>
> >>
> >>
> >> +-------------------------------------+--------------------------+
> >> | Office: | Home: |
> >> |-------------------------------------|--------------------------|
> >> | Dr. Michael Kohl | Michael Kohl |
> >> | Laboratory for Nuclear Science | 5 Ibbetson Street |
> >> | MIT-Bates Linear Accelerator Center | Somerville, MA 02143 |
> >> | Middleton, MA 01949 | U.S.A. |
> >> | U.S.A. | |
> >> | - - - - - - - - - - - - | - - - - - - - - -|
> >> | Email: kohlm@mit.edu | K.Michael.Kohl@gmx.de |
> >> | Work: +1-617-253-9207 | Home: +1-617-629-3147 |
> >> | Fax: +1-617-253-9599 | Mobile: +1-978-580-4190 |
> >> | http://blast.lns.mit.edu | |
> >> +-------------------------------------+--------------------------+
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >
> >
> >
>
>
> +-------------------------------------+--------------------------+
> | Office: | Home: |
> |-------------------------------------|--------------------------|
> | Dr. Michael Kohl | Michael Kohl |
> | Laboratory for Nuclear Science | 5 Ibbetson Street |
> | MIT-Bates Linear Accelerator Center | Somerville, MA 02143 |
> | Middleton, MA 01949 | U.S.A. |
> | U.S.A. | |
> | - - - - - - - - - - - - | - - - - - - - - -|
> | Email: kohlm@mit.edu | K.Michael.Kohl@gmx.de |
> | Work: +1-617-253-9207 | Home: +1-617-629-3147 |
> | Fax: +1-617-253-9599 | Mobile: +1-978-580-4190 |
> | http://blast.lns.mit.edu | |
> +-------------------------------------+--------------------------+
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
--------------------------------------------------------------------------
Eugene Geis
PhD Student, Physics Department, ASU
Research Affiliate, MIT-Bates Laboratory of Nuclear Science
eugene.geis@asu.edu
--------------------------------------------------------------------------
http://quickreaction.blogspot.com
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:33 EST