-- 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 - 2009-11-10, BeataMichalska
Copyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Foswiki Send feedback | Imprint | Privacy Policy (in German)