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 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