Difference: BottomUpUncerts (7 vs. 8)

Revision 82018-04-19 - BenjaminNachman

Line: 1 to 1
 
META TOPICPARENT name="BenjaminNachman"
<!-- This is the default ATLAS template. 
Please modify it in the sections indicated to create your topic. In particular, notice that at the bottom there are some sections that must be filled for publicly accessible pages.
Line: 191 to 191
 This is clearly something that could be important and we have no measurements of this so we need to make sure that we aren't sensitive to these effects. First, define sigma to be the effective cluster size (Eq. 20 here):
Changed:
<
<
sigma = atan(sqrt(SECOND_R)/CENTER_MAG)*cosh(eta)
>
>
sigma = atan(sqrt(cluster->auxdata("SECOND_R"))/cluster->auxdata("CENTER_MAG"))*cosh(eta)
 
Changed:
<
<
Then, one should check three things:
>
>
Also define:
 
Changed:
<
<
  1. test splitting between HAD and EM; for any pair of clusters where dR(i,j)<sigma_{eta/phi}^{i}+sigma_{eta/phi}^{j} && ((E_EMB3^i =0 && E_EMB3^j =0) || (E_Tile0^i =0 && E_Tile0^j =0) || (E_EME3^i =0 && E_EME3^j =0) || (E_HEC0^i =0 && E_HEC0^j =0)) - this condition basically asks if the clusters are overlapping and they both have energy in the first HAD layer or last EM layer (so are touching each other). These should be added together to form a single cluster with eta=(eta_i*E_i+eta_j*E_j)/(E_i+E_j), phi=(phi_i*E_i+phi_j*E_j)/(E_i+E_j), E = E_i+E_j TODO BPN> check if variables are present
>
>
E_EMB3 = cluster->auxdata<vector<float> >("EnergyPerSampling")[3]; 
E_Tile0 = cluster->auxdata<vector<float> >("EnergyPerSampling")[12];
E_EME3 = cluster->auxdata<vector<float> >("EnergyPerSampling")[7]; 
E_HEC0 = cluster->auxdata<vector<float> >("EnergyPerSampling")[8];

Then, we have the following procedure

  1. Splitting between HAD and EM:

vector merged_constituents;
vector already_merged;
for i in constituents do
      already_merged[i] = false
end do

for i in constituents do
      if (already_merged[i])
          continue
      end if
      for j in constituents do
          if i == j
              continue
          if (DeltaR(constituent[i], constituent[j]) > sigma[i]+sigma[j]) #close to each other transversely
              continue
          if  ((E_EMB3[i] !=0 && E_EMB3[j] !=0) || (E_Tile0[i] !=0 && E_Tile0[j] !=0) || (E_EME3[i] !=0 && E_EME3[j] !=0) || (E_HEC0[i] !=0 && E_HEC0[j] !=0)) then #This condition basically asks if the clusters are overlapping and they both have energy in the first HAD layer or last EM layer (so are touching each other).
              already_merged[i] = true
              already_merged[j] = true
              #merged_constituent has eta = (eta[i]*E[i]+eta[j]*E[j])/(E[i]+E[j]), phi = (phi[i]*E[i]+phi[j]*E[j])/(E[i]+E[j]), E = E[i]+E[j]
              merged_constituents.push_back(merged_constituent)
          end if
     end do
     if (!already_merged[i])
        merged_constituents.push_back(constituents[i]);
     end if
end do
 
  1. splitting in the EM calorimeter of hadronic clusters; for any pair of clusters with the maximum energy in EMB2 or EME2, and EM_PROBABILITY<0.9 (not extremely likely to be a photon), and dR(i,j)<sigma_{eta/phi}^{i}+sigma_{eta/phi}^{j} then we combine them into a single cluster with eta=(eta_i*E_i+eta_j*E_j)/(E_i+E_j), phi=(phi_i*E_i+phi_j*E_j)/(E_i+E_j), E = E_i+E_j
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback