def makedist(coeff_barl, coeff_endc):
ebmap = TH2F("EB", "EB", 360, 1, 261, 171, -85, 86)
eePmap = TH2F("EE", "EE", 100, 1, 101, 100, 1, 101)
eeMmap = TH2F("EEminus", "EEminus", 100, 1, 101, 100, 1, 101)
ebdist = TH1F("EBdist", "EBdist", 100, -2, 2)
ebBorderdist = TH1F("EBBorderdist", "EBBorderdist", 100, -2, 2)
ebeta = TH2F("ebeta", "ebeta", 171, -85, 86, 100, -2, 2)
ebphi = TH2F("ebphi", "ebphi", 360, 1, 361, 100, -2, 2)
eePL = TH2F("EEPL", "EEPlus Left", 50, 10, 55, 100, -2, 2)
eePR = TH2F("EEPR", "EEPlus Right", 50, 10, 55, 100, -2, 2)
eeML = TH2F("EEML", "EEMinus Left", 50, 10, 55, 100, -2, 2)
eeMR = TH2F("EEMR", "EEMinus Right", 50, 10, 55, 100, -2, 2)
for i, c in enumerate(coeff_barl):
ieta, iphi = EcalPyUtils.unhashEBIndex(i)
ebmap.Fill(iphi, ieta, c)
ebdist.Fill(c)
ebeta.Fill(ieta, c)
ebphi.Fill(iphi, c)
if (abs(ieta) == 85 or abs(ieta) == 65 or abs(ieta) == 64
or abs(ieta) == 45 or abs(ieta) == 44 or abs(ieta) == 25
or abs(ieta) == 24 or abs(ieta) == 1 or iphi % 20 == 1
or iphi % 20 == 0):
ebBorderdist.Fill(c)
for i, c in enumerate(coeff_endc):
ix, iy, iz = EcalPyUtils.unhashEEIndex(i)
R = sqrt((ix - 50) * (ix - 50) + (iy - 50) * (iy - 50))
if iz > 0:
eePmap.Fill(ix, iy, c)
if ix < 50:
eePL.Fill(R, c, 1)
if ix > 50:
eePR.Fill(R, c, 1)
if iz < 0:
eeMmap.Fill(ix, iy, c)
if ix < 50:
eeML.Fill(R, c, 1)
if ix > 50:
eeMR.Fill(R, c, 1)
prof_eePL = eePL.ProfileX()
prof_eePR = eePR.ProfileX()
prof_eeML = eeML.ProfileX()
prof_eeMR = eeMR.ProfileX()
return ebmap, ebeta, ebphi, eePmap, ebdist, eeMmap, prof_eePL, prof_eePR, prof_eeML, prof_eeMR, ebBorderdist
def makedist(coeff_barl, coeff_endc) :
ebmap = TH2F("EB","EB",360,1,261,171, -85,86)
eePmap = TH2F("EE","EE",100, 1,101,100,1,101)
eeMmap = TH2F("EEminus","EEminus",100,1,101,100,1,101)
ebdist = TH1F("EBdist","EBdist",100,-2,2)
ebBorderdist = TH1F("EBBorderdist","EBBorderdist",100,-2,2)
ebeta = TH2F("ebeta","ebeta",171,-85,86,100,-2,2)
ebphi = TH2F("ebphi","ebphi",360,1,361,100,-2,2)
eePL = TH2F("EEPL","EEPlus Left",50,10,55,100,-2,2)
eePR = TH2F("EEPR","EEPlus Right",50,10,55,100,-2,2)
eeML = TH2F("EEML","EEMinus Left",50,10,55,100,-2,2)
eeMR = TH2F("EEMR","EEMinus Right",50,10,55,100,-2,2)
for i,c in enumerate(coeff_barl):
ieta,iphi = EcalPyUtils.unhashEBIndex(i)
ebmap.Fill(iphi,ieta,c)
ebdist.Fill(c)
ebeta.Fill(ieta,c)
ebphi.Fill(iphi,c)
if (abs(ieta)==85 or abs(ieta)==65 or abs(ieta)==64 or abs(ieta)==45 or abs(ieta)==44 or abs(ieta)==25 or abs(ieta)==24 or abs(ieta)==1 or iphi%20==1 or iphi%20==0):
ebBorderdist.Fill(c)
for i,c in enumerate(coeff_endc):
ix,iy,iz = EcalPyUtils.unhashEEIndex(i)
R = sqrt((ix-50)*(ix-50)+(iy-50)*(iy-50))
if iz>0:
eePmap.Fill(ix,iy,c)
if ix<50:
eePL.Fill(R,c,1)
if ix>50:
eePR.Fill(R,c,1)
if iz<0:
eeMmap.Fill(ix,iy,c)
if ix<50:
eeML.Fill(R,c,1)
if ix>50:
eeMR.Fill(R,c,1)
prof_eePL = eePL.ProfileX()
prof_eePR = eePR.ProfileX()
prof_eeML = eeML.ProfileX()
prof_eeMR = eeMR.ProfileX()
return ebmap, ebeta, ebphi, eePmap, ebdist, eeMmap, prof_eePL, prof_eePR, prof_eeML, prof_eeMR, ebBorderdist
def compareBarrelBorder(coeff_barl_1,coeff_barl_2) :
'''Return an histogram and a map of the differences '''
diff_distro_border_h = TH1F("diffborderh","diffh",100,-2,2)
diff_distro_border_m = TH2F("diffborderm","diffm",360,1,361,171,-85,86)
diff_distro_border_m.SetStats(0)
for i,c in enumerate(coeff_barl_1):
diff = c - coeff_barl_2[i]
ieta,iphi= EcalPyUtils.unhashEBIndex(i)
if (abs(ieta)==85 or abs(ieta)==65 or abs(ieta)==64 or abs(ieta)==45 or abs(ieta)==44 or abs(ieta)==25 or abs(ieta)==24 or abs(ieta)==1 or iphi%20==1 or iphi%20==0):
diff_distro_border_h.Fill(diff)
if (abs(ieta)==85 or abs(ieta)==65 or abs(ieta)==64 or abs(ieta)==45 or abs(ieta)==44 or abs(ieta)==25 or abs(ieta)==24 or abs(ieta)==1 or iphi%20==0 or iphi%20==1):
diff_distro_border_m.Fill(iphi,ieta,diff)
return diff_distro_border_h, diff_distro_border_m
def compareBarrelBorder(coeff_barl_1,coeff_barl_2) :
'''Return an histogram and a map of the differences '''
diff_distro_border_h = TH1F("diffborderh","diffh",100,-2,2)
diff_distro_border_m = TH2F("diffborderm","diffm",360,1,361,171,-85,86)
diff_distro_border_m.SetStats(0)
for i,c in enumerate(coeff_barl_1):
diff = c - coeff_barl_2[i]
ieta,iphi= EcalPyUtils.unhashEBIndex(i)
if (abs(ieta)==85 or abs(ieta)==65 or abs(ieta)==64 or abs(ieta)==45 or abs(ieta)==44 or abs(ieta)==25 or abs(ieta)==24 or abs(ieta)==1 or iphi%20==1 or iphi%20==0):
diff_distro_border_h.Fill(diff)
if (abs(ieta)==85 or abs(ieta)==65 or abs(ieta)==64 or abs(ieta)==45 or abs(ieta)==44 or abs(ieta)==25 or abs(ieta)==24 or abs(ieta)==1 or iphi%20==0 or iphi%20==1):
diff_distro_border_m.Fill(iphi,ieta,diff)
return diff_distro_border_h, diff_distro_border_m
def compareBarrel(coeff_barl_1, coeff_barl_2):
'''Return an histogram and a map of the differences '''
diff_distro_h = TH1F("diffh", "diffh", 100, -2, 2)
diff_distro_m = TH2F("diffm", "diffm", 360, 1, 361, 171, -85, 86)
diff_distro_m.SetStats(0)
ebeta = TH2F("ebeta", "ebeta", 171, -85, 86, 100, -2, 2)
ebphi = TH2F("ebphi", "ebphi", 360, 1, 361, 100, -2, 2)
for i, c in enumerate(coeff_barl_1):
diff = c - coeff_barl_2[i]
ieta, iphi = EcalPyUtils.unhashEBIndex(i)
diff_distro_h.Fill(diff)
diff_distro_m.Fill(iphi, ieta, diff)
ebeta.Fill(ieta, diff)
ebphi.Fill(iphi, diff)
prof_x_h = ebeta.ProfileX()
prof_y_h = ebphi.ProfileX()
return diff_distro_h, diff_distro_m, prof_x_h, prof_y_h
def compareBarrel(coeff_barl_1,coeff_barl_2) :
'''Return an histogram and a map of the differences '''
diff_distro_h = TH1F("diffh","diffh",100,-2,2)
diff_distro_m = TH2F("diffm","diffm",360,1,361,171,-85,86)
diff_distro_m.SetStats(0)
ebeta = TH2F("ebeta","ebeta",171,-85,86,100,-2,2)
ebphi = TH2F("ebphi","ebphi",360,1,361,100,-2,2)
for i,c in enumerate(coeff_barl_1):
diff = c - coeff_barl_2[i]
ieta,iphi= EcalPyUtils.unhashEBIndex(i)
diff_distro_h.Fill(diff)
diff_distro_m.Fill(iphi,ieta,diff)
ebeta.Fill(ieta,diff)
ebphi.Fill(iphi,diff)
prof_x_h = ebeta.ProfileX()
prof_y_h = ebphi.ProfileX()
return diff_distro_h, diff_distro_m, prof_x_h, prof_y_h
