-- BeataMichalska - 24 Sep 2009 + Matrices

Used cuts:
  1. graphical cut on gamma convertion from the flanch around Rich
  2. graphical PID cut
  3. closest neibourgh cut /9 degree/
  4. ep_sim_parentid ==-1 /primary track/ ++ Acceptancy matrix for electron:

  • electronAcc: p vs phi:
    1

  • electronAcc: p vs th:
    2

  • electronAcc: th vs phi:
    3

++ Efficiency matrix for electron:

  • electronEff: p vs phi:
    4

  • electronEff: p vs th:
    5

  • electronEff: th vs phi:
    6

++ Acceptancy matrix for positron:

  • positronAcc: p vs phi:
    10

  • positronAcc: p vs th:
    11

  • positronAcc: th vs phi:
    12

++ Efficiency matrix for positron:

  • positronEff: p vs phi:
    7

  • positronEff: p vs th:
    8

  • positronEff: th vs phi:
    9

+ SELF CONSISTENCY CHECK

  1. The comparison of theta, momentum and phi for one file: be06130200015

red - root file generated by pluto filtered through matrices

black - pluto simulation + hgeant + PAT analysis

++ positrons

+++++ theta


  • filtering with acc & eff matrices:
    filtrowanie macierzami acc i eff bez dodatkowych warunkow

  • filtering with acc & eff matrices, additionaly momentum >100 was required:
    filtrowanie macierzami acc i eff, dodatkowo, p>100

  • filtering with acc & eff matrices, additionaly transverse momentum <100 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<100:

  • filtering with acc & eff matrices, additionaly momentum >100 and transverse momentum <100 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<100 i p>100

  • filtering with acc & eff matrices, additionaly momentum >50 and transverse momentum <600 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<600 i p>50

+++++ momentum


  • filtering with acc & eff matrices:
    filtrowanie macierzami acc i eff bez dodatkowych warunkow

  • filtering with acc & eff matrices, additionaly transverse momentum <100 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<100

  • filtering with acc & eff matrices, additionaly transverse momentum <600 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<600

+++++ phi


  • filtering with acc & eff matrices:
    filtrowanie macierzami acc i eff bez dodatkowych warunkow

  • filtering with acc & eff matrices, additionaly transverse momentum <100 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<100

  • filtering with acc & eff matrices, additionaly momentum >100 was required:
    filtrowanie macierzami acc i eff, dodatkowo, p>100

  • filtering with acc & eff matrices, additionaly momentum >100 and transverse momentum <100 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<100 i p>100

  • filtering with acc & eff matrices, additionaly momentum >50 and transverse momentum <600 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<600 i p>50

++ electrons

+++++theta


  • filtering with acc & eff matrices:
    filtrowanie macierzami acc i eff bez dodatkowych warunkow

  • filtering with acc & eff matrices, additionaly transverse momentum <100 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<100

  • filtering with acc & eff matrices, additionaly momentum >100 was required:
    filtrowanie macierzami acc i eff, dodatkowo, p>100

  • filtering with acc & eff matrices, additionaly momentum >50 and transverse momentum <600 was required:<br /> filtrowanie macierzami acc i eff, dodatkowo, pt<600 i p>50

++ momentum >50 & theta >15 & theta < 85 was required +++ positrons

  • momentum:
    momentum

  • theta:
    theta

  • phi:
    phi

++ momentum >50 & theta >15 & theta < 85 & (in each sector phi >5 && phi<60) was required +++ positrons

  • momentum:
    momentum

  • theta:
    theta

  • phi:
    phi

+ new pictures added by Radek.
  1. Eff. matrices were changed: phi axis range from 0 to 60 degree,
  2. error bars was added to histograms
  3. only for positrons
  4. momentum,theta,phi in the filtering simulation were smeared

++ no conditions

  • momentum,no conditions:
    ep_mom_acc_eff

  • phi, no conditions:
    ep_phi_acc_eff

  • theta, no conditions:
    ep_theta_acc_eff

++ pt < 100

  • momentum, pt < 100:
    ep_mom_acc_eff_ptMN100

  • phi, pt < 100:
    ep_phi_acc_eff_ptMN100

  • theta, pt < 100:
    ep_theta_acc_eff_ptMN100

++ pt < 600

  • momentum, pt < 600:
    ep_mom_acc_eff_ptMN600

  • phi, pt < 600:
    ep_phi_acc_eff_ptMN600

  • theta, pt < 600:
    ep_theta_acc_eff_ptMN600

++ p > 50

  • momentum, p > 50:
    ep_mom_acc_eff_pWI50

  • phi, p>50:
    ep_phi_acc_eff_pWI50

  • theta, p > 50:
    ep_theta_acc_eff_pWI50

++ p > 50 && pt < 600

  • momentum, p > 50 && pt < 600:
    ep_mom_acc_eff_pWI50_ptMN600

  • phi, p > 50 && pt < 600:
    ep_phi_acc_eff_pWI50_ptMN600

  • theta, p > 100:
    ep_theta_acc_eff_pWI100

++ p > 100

  • momentum, p > 100:
    ep_mom_acc_eff_pWI100

  • phi, p > 100:
    ep_phi_acc_eff_pWI100

  • theta, p > 50 && pt < 600:
    ep_theta_acc_eff_pWI50_ptMN600

+ 01.10.09
  1. rebin histograms in order to reach to compatibility with the bining of efficiency matrices(acc. matr. have different bining!!!)
  2. again eff matrices have 0-360 phi range
  3. these matrices are using reconstruct variable
  4. smearing in filtered ntuple is ON

++ no conditions

  • mom: with ACC and EFF:
    ep_mom_acc_eff

  • phi:
    ep_phi_acc_eff

  • theta:
    ep_theta_acc_eff

++ pt < 100

  • mom: pt < 100:
    ep_mom_acc_eff_ptMN100

  • phi: pt < 100:
    ep_phi_acc_eff_ptMN100

  • theta: pt < 100:
    ep_theta_acc_eff_ptMN100

++ pt < 600

  • mom: pt < 600:
    ep_mom_acc_eff_ptMN600

  • phi pt < 600:
    ep_phi_acc_eff_ptMN600

  • theta pt < 600:
    ep_theta_acc_eff_ptMN600

++ p > 50

  • mom: p > 50:
    ep_mom_acc_eff_pWI50

  • phi: p > 50:
    ep_phi_acc_eff_pWI50

  • theta: p > 50:
    ep_theta_acc_eff_pWI50

++ p > 50 && pt < 600

  • mom: p > 50 && pt < 600:
    ep_mom_acc_eff_pWI50_ptMN600

  • phi: p > 50 && pt < 600:
    ep_phi_acc_eff_pWI50_ptMN600

  • theta: p > 50 && pt < 600:
    ep_theta_acc_eff_pWI50_ptMN600

++ p > 100

  • mom: p > 100:
    ep_mom_acc_eff_pWI100

  • phi: p > 100:
    ep_phi_acc_eff_pWI100

  • theta: p > 100:
    ep_theta_acc_eff_pWI100

++ no conditions, only ACC && EFF filtered

* p vs pt, filtr:
filt_ep_p_pt_acc_eff

  • p vs th, filtr:
    filt_ep_p_th_acc_eff

  • th vs pt, filtr:
    filt_ep_th_pt_acc_eff

  • p vs pt, pat:
    pat_ep_p_pt_acc_eff

  • p vs th, pat:
    pat_ep_p_th_acc_eff

  • th vs pt, pat:
    pat_ep_th_pt_acc_eff

  • (pR-pI)/pR:
    pat_ep_pR-pI_pR_acc_eff

  • (pR-pI)/pR zoom:
    pat_ep_pR-pI_pR_acc_eff3

++ no conditions, ONLY ACC filtered, w/o efficiency

  • mom:
  • PAT draw from /lustre/hades/user/beatam/Embedding/out_ntupleMianownik/be06130200015_eff_ntuple_ep.root
    and /lustre/hades/user/beatam/Embedding/out_ntupleMianownik_z1/be06130200015_eff_ntuple_ep_1.root
    without conditions (isBest==1)&&(ep_sim_iscommon>=76)&&(ep_sim_parentid==-1)
    ep_mom_acc3

  • theta:
    ep_th_acc

  • phi:
    ep_ph_acc

+ systematic shift of reconstracted momentum

  • (reconstr mom - ideal mom) / ideal mom:
    pat_ep_pR-pI_pR_acc_eff2

  • slice_momentum 30 GeV:
    slice_ped30

  • slice_momentum 50 GeV:
    slice_ped50

  • slice_momentum 70 GeV:
    slice_ped70

  • slice_momentum 150 GeV:
    slice_ped150

  • slice_momentum 590 GeV:
    slice_ped590

  • slice_momentum 1130 GeV:
    slice_ped1130

  • slice_momentum 1690 GeV:
    slice_ped1690

  • slice_momentum 1950 GeV:
    slice_ped1950

  • Pideal - P recons / P ideal:
    mean_Pideal_precons_acc_eff

  • (reconstr mom - ideal mom):

  • pat_ep_pR-pI_acc_eff:
    pat_ep_pR-pI_acc_eff

  • slice_ped_230:
    slice_ped_230

  • slice_ped_250:
    slice_ped_250

  • slice_ped_2490:
    slice_ped_2490

+ projections of apr07 acc matrix

  • apr07_phi_p:
    apr07_phi_p

  • apr07_ph_th:
    apr07_ph_th

  • apr07_th_p:
    apr07_th_p

+ New ACC matr. with new bining/ compatible with EFF matr./ were used

  • positron: phi:
    ep_ph_accNB_M3

  • positron: mom:
    ep_mom_accNB_M

  • positron: th:
    ep_th_accNB_M

+ EFF added

  • ep_mom_acc_eff:
    ep_mom_acc_eff

  • ep_phi_acc_eff:
    ep_phi_acc_eff

  • ep_theta_acc_eff:
    ep_theta_acc_eff

+ new bining for ACC MATRIX - positron

  • ep_acc_p_phi:
    ep_acc_p_phi

  • ep_acc_p_th:
    ep_acc_p_th

  • ep_acc_th_phi:
    ep_acc_th_phi

+ after EFF, no condition

  • ep_mom_acc_eff3:
    ep_mom_acc_eff3

  • ep_phi_acc_eff3:
    ep_phi_acc_eff3

  • ep_theta_acc_eff3:
    ep_theta_acc_eff3

+ ACC + EFF

  1. acc,eff matrices based on 100 files
  2. ideal params
  3. ep_sim_iscommon>=76
  4. NoHadronE1 && NoLeptonE1
  5. ep_sim_parentid ==-1
  6. isBest==1
  7. no Positron condition

  • ep_mom_acc_eff7:
    ep_mom_acc_eff7

  • ep_theta_acc_eff7:
    ep_theta_acc_eff7

+++ lets check with smearing

+ after Smearing the Eff matrix and data from PAT

  • ep_mom_acc_effS:
    ep_mom_acc_effS

  • ep_phi_acc_effS:
    ep_phi_acc_effS

  • ep_theta_acc_effS:
    ep_theta_acc_effS

+ momentum

  • pI_iR_th_2D:
    pI_iR_th_2D

  • pI_iR_pI_2D:
    pI_iR_pI_2D

  • pI_iR_pI_2D_slices:
    pI_iR_pI_2D_slices

  • pI_iR_th2:
    pI_iR_th2

+ theta [degree]

  • thI_thR_pI_2D:
    thI_thR_pI_2D

  • thI_thR_pI_2D_slices:
    thI_thR_pI_2D_slices

+ phi [degree]

  • phI_phR_pI_2D:
    phI_phR_pI_2D

  • phI_phR_pI_2D_slices2:
    phI_phR_pI_2D_slices2

+ after smearing with the above histograms

  • ep_mom_acc_effSm2:
    ep_mom_acc_effSm2

  • ep_phi_acc_effSm2:
    ep_phi_acc_effSm2

  • ep_theta_acc_effSm2:
    ep_theta_acc_effSm2

+ after smearing with the original histograms + pI_iR_pI_2D, phI_phR_pI_2D and thI_thR_pI_2D. +++Discrepancy is about 1%.

  • ep_mom_acc_eff9.gif:
    ep_mom_acc_eff9.gif

  • ep_theta_acc_eff9.gif:
    ep_theta_acc_eff9.gif

  • ep_phi_acc_eff9.gif:
    ep_phi_acc_eff9.gif

+ no shift to the direction of negative momentum [smearing: hist->GetRandom()-hist->GetMean()]

  • ep_mom_acc_eff6:
    ep_mom_acc_eff6

  • ep_phi_acc_eff6:
    ep_phi_acc_eff6

  • ep_theta_acc_eff6:
    ep_theta_acc_eff6

+ smearing: hist->GetRandom(); eff taken from 1Dim-EFF(p)

  • ep_mom_acc_eff1d:
    ep_mom_acc_eff1d

  • ep_phi_acc_eff1d:
    ep_phi_acc_eff1d

  • ep_theta_acc_eff1d:
    ep_theta_acc_eff1d

+ standard 3D matrices; smearing only for momentum

  • ep_mom_acc_effOM:
    ep_mom_acc_effOM

  • ep_phi_acc_effOM:
    ep_phi_acc_effOM

  • ep_theta_acc_effOM:
    ep_theta_acc_effOM

+ standard 3D matrices; smearing momentum in the standard way, theta/phi from gaus(0.023)

  • ep_mom_acc_effG:
    ep_mom_acc_effG

  • ep_phi_acc_effG:
    ep_phi_acc_effG

  • ep_theta_acc_effG:
    ep_theta_acc_effG

+ PAT mom ideal vs reconstructed. Yields are identical

  • pat momentum ideal vs zrek:
    pat momentum ideal vs zrek

+ filtr: 1. eff_ideal, 2. smear; PAT: recons == 1. smear; 2. eff_ideal
  • p_mom_acc_effs.:
    p_mom_acc_effs.

+ filtr: 1. eff. reconst, 2. smear; PAT: recons

  • ep_mom_acc_eff_effR_Sm:
    ep_mom_acc_eff_effR_Sm

smearing mom hists extend to -2000,2000. First getEff[ideal], next smearing

Procedure: Acc matrix - based on ideal parameter Eff matrix - based on ideal parameter conditions:
  1. (((ep_system==0&&insideEpS0==0)||(ep_system==1&&insideEpS1==0)))
  2. ep_sim_iscommon>=76
  3. !(ep_oa_lept< 9 &&ep_oa_lept>-9 )
  4. !(ep_oa_hadr< 9 &&ep_oa_hadr>-9 )
  5. ep_sim_parentid ==-1
  6. isBest==1
  7. x smearing histograms for momentum(without any limits), theta, phi takes as reconstructed_parameter - ideal_parameter. Smearing is based on the sampling of the value from these histograms, which is subtracted from real value taken from simulation. Procedure:
  8. filtering with ACC matrix
  9. filtering with EFF matrix
  10. Smearing

- compared with PAT reconstracted

  • ep_mom_acc_effM:
    ep_mom_acc_effM

  • ep_phi_acc_effM:
    ep_phi_acc_effM

  • ep_theta_acc_effM:
    ep_theta_acc_effM

+ electrons

  • em_mom_acc_eff:
    em_mom_acc_eff

  • em_phi_acc_eff:
    em_phi_acc_eff

  • em_theta_acc_eff:
    em_theta_acc_eff

+ protons
condition
p_sim_iscommon>=74 instead of p_sim_iscommon>=76
  • p_mom_acc_eff:
    p_mom_acc_eff

  • p_theta_acc_eff:
    p_theta_acc_eff

  • p_phi_acc_eff2:
    p_phi_acc_eff2

+ 22.10

  • ep_pat_pseudoideal_cor_eff_from_pseudoideal_effG0.01:
    ep_pat_pseudoideal_cor_eff_from_pseudoideal_effG0.01

  • test_SIGNAL_exp_sim:
    test_SIGNAL_exp_sim

  • test_SIGNAL_fiducial_exp_sim:
    test_SIGNAL_fiducial_exp_sim + resmearing histograms to [0,3000]

  • filtr: fitrACC ; PAT: EFFcorr for pseudoideal; draw pseudoideal:
    new_ideal.

  • filtr: filtACC, smearing ; PAT: EFFcorr for pseudoideal; draw recons:
    new_rekon

++ EXP & SIM 4prong analysis
  • SIGNIFICANCE_novert_Fmala_mass:
    SIGNIFICANCE_novert_Fmala_mass

  • test_SIGNAL_exp_sim:
    test_SIGNAL_exp_sim

  • test_SIGNAL_fiducial_exp_sim:
    test_SIGNAL_fiducial_exp_sim

  • EXP, helicity after correction; fit in range (-0.8,0.8)
    Anizotropic factor B = 0.92 ± 0.17; 1+B*x2
    08

  • helicity EXP(black)_SIM(red):
    helicity_poprawionyExp_sim

  • helicity after correction; fit in range (-0.6,0.6)
    Anizotropic factor B = 1.26 ± 0.29; 1+B*x2
    :
    helicity after correction fit in range (-0.6,0.6) <br />Anizotropic factor B = 1.26 ± 0.29<br />

  • cos(theta of pi0 in CM); pi0 was boosted to CM:
    pi0 distribution

  • cos(theta of pi0 in CM) in fiducial area: em/ep_mom >100MeV && em/ep_theta>20MeV:
    pi0 distribution in fiducial area: MOM100_TH20

  • Delta1_invMass: delta1 = em + ep + first proton
    Delta1_invMass

  • helicity correction from Hades acc to 4pi:
    helicity_poprawka

Files on the tape: tsmcli archive "*root" hadesuser michalsk/eff_mat/out_geant_pip tsmcli archive "*root" hadesuser michalsk/eff_mat/out_geant_p tsmcli archive "*root" hadesuser michalsk/eff_mat/out_geant_ep tsmcli archive "*root" hadesuser michalsk/eff_mat/out_geant_em

tsmcli archive "*root" hadesuser michalsk/eff_mat/out_embedded tsmcli archive "*root" hadesuser michalsk/eff_mat/out_embedded_em tsmcli archive "*root" hadesuser michalsk/eff_mat/out_embedded_p tsmcli archive "*root" hadesuser michalsk/eff_mat/out_embedded_pip

=====================
beatam@yeti:/net/disk2/beatam/out_PAT4_bind was vertex ntuple created thay were copied to /lustre/hades/user/beatam/Embedding/out_vertex_from_dst

Topic revision: r50 - 10 Nov 2009, BeataMichalska
 
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Hades Wiki? Send feedback
Imprint (in German)
Privacy Policy (in German)