Line: 1 to 1  

<! 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  
Changed:  
< <  Then, one should check three things:  
> >  Also define:  
Changed:  
< < 
 
> >  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
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  
