RecEmcShower.cxx File Reference

#include <iostream>
#include <complex>
#include "EmcRecEventModel/RecEmcShower.h"

Go to the source code of this file.

Functions
ostream &	operator<< (ostream &os, const RecEmcShower &aShower)

Function Documentation

operator<<()

ostream & operator<<	(	ostream &	os,
		const RecEmcShower &	aShower
	)

Definition at line 326 of file RecEmcShower.cxx.

{
  RecEmcFractionMap::const_iterator pFractionMap;
  RecEmcFractionMap::const_iterator ciFractionMap3x3;
  RecEmcFractionMap::const_iterator ciFractionMap5x5;
 
  os<<"------------------RecEmcShower:"<<endl;
  os<<"Track id: "<<aShower.trackId()<<", number of hits: "<<aShower.getSize()
    <<", status: "<<aShower.status()<<", shower id: "<<aShower.getShowerId()
    <<", cluster id: "<<aShower.getClusterId()<<", module: "<<aShower.module()
    <<", time: "<<aShower.time()<<endl;
  os<<"Energy: "<<aShower.energy()<<", de: "<<aShower.dE()
    <<", eseed: "<<aShower.eSeed()
    <<", e3x3: "<<aShower.e3x3()<<", e5x5: "<<aShower.e5x5()
    <<", eall: "<<aShower.getEAll()<<", elepton: "<<aShower.getELepton()<<endl;
  os<<"Position: "<<aShower.position()
    <<", theta: "<<aShower.theta()<<", phi: "<<aShower.phi()
    <<"\ndx: "<<aShower.dx()<<", dy: "<<aShower.dy()<<", dz: "<<aShower.dz()
    <<", dtheta: "<<aShower.dtheta()<<", dphi: "<<aShower.dphi()<<endl;
  os<<"Second moment: "<<aShower.secondMoment()
    <<", lat moment: "<<aShower.latMoment()
    <<", a20 moment: "<<aShower.a20Moment()
    <<", a42 moment: "<<aShower.a42Moment()<<endl;
  os<<"Error matrix: \n"<<aShower.errorMatrix()<<endl;
 
  if(aShower.getSize()>0) {
    os<<"Fraction Map: "<<endl;
    for(pFractionMap=aShower.Begin();
        pFractionMap!=aShower.End();
        pFractionMap++){
      os<< (pFractionMap->second);
    }
  }
 
  RecEmcFractionMap fracMap = aShower.getFractionMap3x3();
  if(fracMap.size()>0) {
    os<<"Fraction Map 3x3: "<<endl; 
    for(ciFractionMap3x3=fracMap.begin();
        ciFractionMap3x3!=fracMap.end();
        ciFractionMap3x3++){
      os<< (ciFractionMap3x3->second);
    }
  }
 
  fracMap = aShower.getFractionMap5x5();
  if(fracMap.size()>0) {
    os<<"Fraction Map 5x5: "<<endl; 
    for(ciFractionMap5x5=fracMap.begin();
        ciFractionMap5x5!=fracMap.end();
        ciFractionMap5x5++){
      os<< (ciFractionMap5x5->second);
    }
  }
 
  os<<endl;
 
  return os;
}