1 #ifndef GAUDISVC_GENERIC2D_H 2 #define GAUDISVC_GENERIC2D_H 1 4 #include "AIDA/IProfile2D.h" 15 # pragma clang diagnostic push 16 # pragma clang diagnostic ignored "-Winconsistent-missing-override" 17 #elif defined( __GNUC__ ) && __GNUC__ >= 5 18 # pragma GCC diagnostic push 19 # pragma GCC diagnostic ignored "-Wsuggest-override" 34 template <
class INTERFACE,
class IMPLEMENTATION>
49 void adoptRepresentation( TObject* rep )
override;
59 AIDA::IAnnotation&
annotation()
override {
return m_annotation; }
61 const AIDA::IAnnotation&
annotation()
const override {
return m_annotation; }
64 const AIDA::IAxis&
xAxis()
const override {
return m_xAxis; }
66 const AIDA::IAxis&
yAxis()
const override {
return m_yAxis; }
68 virtual int rIndexX(
int index )
const {
return m_xAxis.rIndex( index ); }
70 virtual int rIndexY(
int index )
const {
return m_yAxis.rIndex( index ); }
73 int entries()
const override;
75 int allEntries()
const override;
77 int extraEntries()
const override;
79 double sumBinHeights()
const override;
81 double sumAllBinHeights()
const override;
85 double minBinHeight()
const override;
87 double maxBinHeight()
const override;
90 double binMeanX(
int indexX,
int indexY )
const override;
92 double binMeanY(
int indexX,
int indexY )
const override;
94 int binEntries(
int indexX,
int indexY )
const override;
96 int binEntriesX(
int indexX )
const override;
98 int binEntriesY(
int indexY )
const override;
100 double binHeight(
int indexX,
int indexY )
const override;
102 double binHeightX(
int indexX )
const override;
104 double binHeightY(
int indexY )
const override;
106 double binError(
int indexX,
int indexY )
const override;
108 virtual double binRms(
int indexX,
int indexY )
const;
110 double meanX()
const override;
112 double meanY()
const override;
114 double rmsX()
const override;
116 double rmsY()
const override;
118 int coordToIndexX(
double coordX )
const override;
120 int coordToIndexY(
double coordY )
const override;
122 virtual double equivalentBinEntries()
const;
125 virtual bool scale(
double scaleFactor );
127 bool add(
const INTERFACE& h )
override;
129 bool reset()
override;
131 void* cast(
const std::string& className )
const override;
141 int write(
const char* file_name )
const override;
155 int m_sumEntries = 0;
158 template <
class INTERFACE,
class IMPLEMENTATION>
160 m_rep->SetTitle( title.
c_str() );
161 if ( !annotation().addItem(
"Title", title ) ) m_annotation.setValue(
"Title", title );
162 if ( !annotation().addItem(
"title", title ) ) annotation().setValue(
"title", title );
166 template <
class INTERFACE,
class IMPLEMENTATION>
168 m_rep->SetName( newName.
c_str() );
169 m_annotation.setValue(
"Name", newName );
173 template <
class INTERFACE,
class IMPLEMENTATION>
175 return m_rep->GetEntries();
178 template <
class INTERFACE,
class IMPLEMENTATION>
180 return m_rep->GetEntries();
183 template <
class INTERFACE,
class IMPLEMENTATION>
185 return m_rep->GetMinimum();
188 template <
class INTERFACE,
class IMPLEMENTATION>
190 return m_rep->GetMaximum();
193 template <
class INTERFACE,
class IMPLEMENTATION>
195 return m_rep->GetSumOfWeights();
198 template <
class INTERFACE,
class IMPLEMENTATION>
200 return m_rep->GetSum();
203 template <
class INTERFACE,
class IMPLEMENTATION>
205 return m_rep->GetBinError( rIndexX( indexX ), rIndexY( indexY ) );
208 template <
class INTERFACE,
class IMPLEMENTATION>
210 return m_rep->GetXaxis()->GetBinCenter( rIndexX( indexX ) );
213 template <
class INTERFACE,
class IMPLEMENTATION>
215 return m_rep->GetYaxis()->GetBinCenter( rIndexY( indexY ) );
218 template <
class INTERFACE,
class IMPLEMENTATION>
221 for (
int iY = -2; iY < yAxis().bins(); ++iY ) n += binEntries( index, iY );
225 template <
class INTERFACE,
class IMPLEMENTATION>
228 for (
int iX = -2; iX < xAxis().bins(); ++iX ) n += binEntries( iX, index );
232 template <
class INTERFACE,
class IMPLEMENTATION>
234 return m_rep->GetBinContent( rIndexX( indexX ), rIndexY( indexY ) );
237 template <
class INTERFACE,
class IMPLEMENTATION>
240 for (
int iY = -2; iY < yAxis().bins(); ++iY ) { s += binHeight( index, iY ); }
244 template <
class INTERFACE,
class IMPLEMENTATION>
247 for (
int iX = -2; iX < xAxis().bins(); ++iX ) s += binHeight( iX, index );
251 template <
class INTERFACE,
class IMPLEMENTATION>
253 return m_rep->GetBinError( rIndexX( indexX ), rIndexY( indexY ) );
256 template <
class INTERFACE,
class IMPLEMENTATION>
258 return m_rep->GetMean( 1 );
261 template <
class INTERFACE,
class IMPLEMENTATION>
263 return m_rep->GetMean( 2 );
266 template <
class INTERFACE,
class IMPLEMENTATION>
268 return m_rep->GetRMS( 1 );
271 template <
class INTERFACE,
class IMPLEMENTATION>
273 return m_rep->GetRMS( 2 );
276 template <
class INTERFACE,
class IMPLEMENTATION>
278 return xAxis().coordToIndex( coord );
281 template <
class INTERFACE,
class IMPLEMENTATION>
283 return yAxis().coordToIndex( coord );
286 template <
class INTERFACE,
class IMPLEMENTATION>
288 const Base* p =
dynamic_cast<const Base*
>( &hist );
289 if ( !p )
throw std::runtime_error(
"Cannot add profile histograms of different implementations." );
294 template <
class INTERFACE,
class IMPLEMENTATION>
296 return binEntries( AIDA::IAxis::UNDERFLOW_BIN, AIDA::IAxis::UNDERFLOW_BIN ) +
297 binEntries( AIDA::IAxis::UNDERFLOW_BIN, AIDA::IAxis::OVERFLOW_BIN ) +
298 binEntries( AIDA::IAxis::OVERFLOW_BIN, AIDA::IAxis::UNDERFLOW_BIN ) +
299 binEntries( AIDA::IAxis::OVERFLOW_BIN, AIDA::IAxis::OVERFLOW_BIN );
302 template <
class INTERFACE,
class IMPLEMENTATION>
304 if ( sumBinHeights() <= 0 )
return 0;
306 m_rep->GetStats( stats );
307 return stats[0] * stats[0] / stats[1];
310 template <
class INTERFACE,
class IMPLEMENTATION>
312 m_rep->Scale( scaleFactor );
316 template <
class INTERFACE,
class IMPLEMENTATION>
323 template <
class INTERFACE,
class IMPLEMENTATION>
326 m_rep->Print(
"all" );
331 template <
class INTERFACE,
class IMPLEMENTATION>
334 s <<
"BinX, BinY, Height, Error " <<
std::endl;
335 for (
int i = 0; i < xAxis().bins(); ++i ) {
336 for (
int j = 0; j < yAxis().bins(); ++j ) {
337 s << binMeanX( i, j ) <<
", " << binMeanY( i, j ) <<
", " << binHeight( i, j ) <<
", " << binError( i, j )
346 template <
class INTERFACE,
class IMPLEMENTATION>
348 TFile* f = TFile::Open( file_name,
"RECREATE" );
349 Int_t nbytes = m_rep->Write();
356 # pragma clang diagnostic pop 357 #elif defined( __GNUC__ ) && __GNUC__ >= 5 358 # pragma GCC diagnostic pop 361 #endif // GAUDIPI_GENERIC2D_H int coordToIndexX(double coordX) const override
Convenience method, equivalent to xAxis().coordToIndex(coord).
virtual double equivalentBinEntries() const
Number of equivalent entries, i.e. SUM[ weight ] ^ 2 / SUM[ weight^2 ]
double rmsX() const override
Returns the rms of the profile as calculated on filling-time projected on the X axis.
int extraEntries() const override
Get the number of entries in the underflow and overflow bins.
int dimension() const override
Get the Histogram's dimension.
Generic2D(IMPLEMENTATION *p)
constructor
double binMeanY(int indexX, int indexY) const override
The weighted mean along y of a given bin.
virtual bool scale(double scaleFactor)
Scale the weights and the errors of all the IHistogram's bins (in-range and out-of-range ones) by a g...
double rmsY() const override
Returns the rms of the profile as calculated on filling-time projected on the Y axis.
double maxBinHeight() const override
Get the maximum height of the in-range bins.
std::string name() const
object name
bool setName(const std::string &newName)
Set the name of the object.
virtual int rIndexX(int index) const
operator methods
int binEntriesX(int indexX) const override
Equivalent to projectionX().binEntries(indexX).
double sumBinHeights() const override
Get the sum of in range bin heights in the IProfile.
double binHeightY(int indexY) const override
Equivalent to projectionY().binHeight(indexY).
double binHeightX(int indexX) const override
Equivalent to projectionX().binHeight(indexX).
std::string m_classType
class type
Axis m_xAxis
X axis member.
Generic2D< INTERFACE, IMPLEMENTATION > Base
int allEntries() const override
Get the number or all the entries, both in range and underflow/overflow bins of the IProfile...
Implementation of the AIDA IAnnotation interface class.
int coordToIndexY(double coordY) const override
Convenience method, equivalent to yAxis().coordToIndex(coord).
double sumExtraBinHeights() const override
Get the sum of the underflow and overflow bin height.
std::ostream & print(std::ostream &s) const override
Print (ASCII) the histogram into the output stream.
const AIDA::IAnnotation & annotation() const override
Access annotation object (cons)
const std::string & userLevelClassType() const
The AIDA user-level unterface leaf class type.
AIDA::IAnnotation & annotation() override
Access annotation object.
double binError(int indexX, int indexY) const override
The error on this bin.
Common base class for all histograms Use is solely functional to minimize dynamic_casts inside Histog...
AIDA::Annotation m_annotation
Object annotations.
bool setTitle(const std::string &title) override
Set the title of the object.
TObject * representation() const override
ROOT object implementation.
double meanX() const override
Returns the mean of the profile, as calculated on filling-time projected on the X axis...
Axis m_yAxis
Y axis member.
std::ostream & write(std::ostream &s) const override
Write (ASCII) the histogram table into the output stream.
double binMeanX(int indexX, int indexY) const override
The weighted mean along x of a given bin.
double sumAllBinHeights() const override
Get the sum of all the bins heights (including underflow and overflow bin).
virtual double binRms(int indexX, int indexY) const
The spread (RMS) of this bin.
std::unique_ptr< IMPLEMENTATION > m_rep
Reference to underlying implementation.
int binEntriesY(int indexY) const override
Equivalent to projectionY().binEntries(indexY).
int entries() const override
Get the number or all the entries.
virtual int rIndexY(int index) const
operator methods
double meanY() const override
Returns the mean of the profile, as calculated on filling-time projected on the Y axis...
std::string title() const override
Get the title of the object.
Helper functions to set/get the application return code.
An IAxis represents a binned histogram axis.
bool add(const INTERFACE &h) override
Modifies this profile by adding the contents of profile to it.
const AIDA::IAxis & yAxis() const override
Return the Y axis.
double minBinHeight() const override
Get the minimum height of the in-range bins.
const AIDA::IAxis & xAxis() const override
Return the X axis.
double binHeight(int indexX, int indexY) const override
Total height of the corresponding bin (ie the sum of the weights in this bin).
Common AIDA implementation stuff for histograms and profiles using ROOT implementations.