Changeset 490 in Sophya for trunk/SophyaLib/NTools/histos2.cc
- Timestamp:
- Oct 21, 1999, 5:25:53 PM (26 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/NTools/histos2.cc
r308 r490 45 45 , hprojx(NULL), hprojy(NULL) 46 46 { 47 DBASSERT(nxBin>0 && nyBin>0 && xMin<xMax && yMin<yMax);47 ASSERT(nxBin>0 && nyBin>0 && xMin<xMax && yMin<yMax); 48 48 for(int i=0;i<3;i++) for(int j=0;j<3;j++) over[i][j]=0.; 49 49 Zero(); … … 663 663 int i,j; 664 664 v.Realloc(nx,ny); 665 if(!err2) for(i=0;i<nx;i++) 666 for(j=0;j<ny;j++) { v(i,j) = 0.; return;} 667 for(i=0;i<nx;i++) 668 for(j=0;j<ny;j++) v(i,j) = Error2(i,j); 665 if(!err2) 666 {for(i=0;i<nx;i++) for(j=0;j<ny;j++) v(i,j) = 0.; return;} 667 for(i=0;i<nx;i++) for(j=0;j<ny;j++) v(i,j) = Error2(i,j); 669 668 return; 670 669 } … … 678 677 int i,j; 679 678 v.Realloc(nx,ny); 680 if(!err2) for(i=0;i<nx;i++) 681 for(j=0;j<ny;j++) { v(i,j) = 0.; return;} 682 for(i=0;i<nx;i++) 683 for(j=0;j<ny;j++) v(i,j) = Error(i,j); 679 if(!err2) 680 {for(i=0;i<nx;i++) for(j=0;j<ny;j++) v(i,j) = 0.; return;} 681 for(i=0;i<nx;i++) for(j=0;j<ny;j++) v(i,j) = Error(i,j); 684 682 return; 685 683 } … … 1786 1784 /////////////////////////////////////////////////////////// 1787 1785 1788 FIO_Histo2D::FIO_Histo2D() 1789 { 1790 dobj=new Histo2D; 1791 ownobj=true; 1792 } 1793 1794 FIO_Histo2D::FIO_Histo2D(string const & filename) 1795 { 1796 dobj=new Histo2D; 1797 ownobj=true; 1798 Read(filename); 1799 } 1800 1801 FIO_Histo2D::FIO_Histo2D(const Histo2D & obj) 1802 { 1803 dobj = new Histo2D(obj); 1804 ownobj=true; 1805 } 1806 1807 FIO_Histo2D::FIO_Histo2D(Histo2D * obj) 1808 { 1809 dobj = obj; 1810 ownobj=false; 1811 } 1812 1813 FIO_Histo2D::~FIO_Histo2D() 1814 { 1815 if (ownobj && dobj) delete dobj; 1816 } 1817 1818 AnyDataObj* FIO_Histo2D::DataObj() 1819 { 1820 return(dobj); 1821 } 1822 1823 void FIO_Histo2D::ReadSelf(PInPersist& is) 1786 1787 void ObjFileIO<Histo2D>::ReadSelf(PInPersist& is) 1824 1788 { 1825 1789 char strg[256]; … … 1879 1843 is.GetLine(strg, 255); 1880 1844 dobj->SetProjX(); 1881 FIO_Histofio_h(dobj->hprojx);1845 ObjFileIO<Histo> fio_h(dobj->hprojx); 1882 1846 fio_h.Read(is); 1883 1847 } … … 1885 1849 is.GetLine(strg, 255); 1886 1850 dobj->SetProjY(); 1887 FIO_Histofio_h(dobj->hprojy);1851 ObjFileIO<Histo> fio_h(dobj->hprojy); 1888 1852 fio_h.Read(is); 1889 1853 } … … 1893 1857 is.GetLine(strg, 255); 1894 1858 dobj->SetSliX(nslix); 1895 DBASSERT (nslix==dobj->NSliX());1859 ASSERT (nslix==dobj->NSliX()); 1896 1860 for(int j=0;j<dobj->NSliX();j++) 1897 { FIO_Histofio_h(dobj->HSliX(j)); fio_h.Read(is);}1861 {ObjFileIO<Histo> fio_h(dobj->HSliX(j)); fio_h.Read(is);} 1898 1862 } 1899 1863 if(nsliy>0) { 1900 1864 is.GetLine(strg, 255); 1901 1865 dobj->SetSliY(nsliy); 1902 DBASSERT (nsliy==dobj->NSliY());1866 ASSERT (nsliy==dobj->NSliY()); 1903 1867 for(int j=0;j<dobj->NSliY();j++) 1904 { FIO_Histofio_h(dobj->HSliY(j)); fio_h.Read(is);}1868 {ObjFileIO<Histo> fio_h(dobj->HSliY(j)); fio_h.Read(is);} 1905 1869 } 1906 1870 … … 1912 1876 dobj->SetBandX(min,max); 1913 1877 }} 1914 DBASSERT (nbanx==dobj->NBandX());1878 ASSERT (nbanx==dobj->NBandX()); 1915 1879 {for(int j=0; j<dobj->NBandX(); j++) { 1916 FIO_Histofio_h(dobj->HBandX(j));1880 ObjFileIO<Histo> fio_h(dobj->HBandX(j)); 1917 1881 fio_h.Read(is); 1918 1882 }} … … 1924 1888 dobj->SetBandY(min,max); 1925 1889 }} 1926 DBASSERT (nbany==dobj->NBandY());1890 ASSERT (nbany==dobj->NBandY()); 1927 1891 {for(int j=0; j<dobj->NBandY(); j++) { 1928 FIO_Histofio_h(dobj->HBandY(j));1892 ObjFileIO<Histo> fio_h(dobj->HBandY(j)); 1929 1893 fio_h.Read(is); 1930 1894 }} … … 1934 1898 } 1935 1899 1936 void FIO_Histo2D::WriteSelf(POutPersist& os) const1900 void ObjFileIO<Histo2D>::WriteSelf(POutPersist& os) const 1937 1901 { 1938 1902 if (dobj == NULL) return; … … 2010 1974 sprintf(strg,"Histo2D: Projection X"); 2011 1975 os.PutLine(strg); 2012 FIO_Histofio_h(dobj->hprojx); fio_h.Write(os);1976 ObjFileIO<Histo> fio_h(dobj->hprojx); fio_h.Write(os); 2013 1977 } 2014 1978 if(projy) { 2015 1979 sprintf(strg,"Histo2D: Projection Y"); 2016 1980 os.PutLine(strg); 2017 FIO_Histofio_h(dobj->hprojy); fio_h.Write(os);1981 ObjFileIO<Histo> fio_h(dobj->hprojy); fio_h.Write(os); 2018 1982 } 2019 1983 … … 2024 1988 for(int j=0;j<nslix;j++) { 2025 1989 Histo* h = dobj->HSliX(j); 2026 FIO_Histofio_h(h); fio_h.Write(os);1990 ObjFileIO<Histo> fio_h(h); fio_h.Write(os); 2027 1991 } 2028 1992 } … … 2032 1996 for(int j=0;j<nsliy;j++) { 2033 1997 Histo* h = dobj->HSliY(j); 2034 FIO_Histofio_h(h); fio_h.Write(os);1998 ObjFileIO<Histo> fio_h(h); fio_h.Write(os); 2035 1999 } 2036 2000 } … … 2047 2011 for(it = dobj->lbandx.begin(); it != dobj->lbandx.end(); it++) { 2048 2012 Histo* h = (*it).H; 2049 FIO_Histofio_h(h); fio_h.Write(os);2013 ObjFileIO<Histo> fio_h(h); fio_h.Write(os); 2050 2014 } 2051 2015 } … … 2060 2024 for(it = dobj->lbandy.begin(); it != dobj->lbandy.end(); it++) { 2061 2025 Histo* h = (*it).H; 2062 FIO_Histofio_h(h); fio_h.Write(os);2026 ObjFileIO<Histo> fio_h(h); fio_h.Write(os); 2063 2027 } 2064 2028 } … … 2066 2030 return; 2067 2031 } 2032 2033 2034 #ifdef __CXX_PRAGMA_TEMPLATES__ 2035 #pragma define_template ObjFileIO<Histo2D> 2036 #endif 2037 2038 #if defined(ANSI_TEMPLATES) || defined(GNU_TEMPLATES) 2039 template class ObjFileIO<Histo2D>; 2040 #endif
Note:
See TracChangeset
for help on using the changeset viewer.