In This Package:

Ana15 Class Reference

Ana15 -- an example to show how to use Fifteen package. More...

#include <Ana15.h>

Inheritance diagram for Ana15:

[legend]Collaboration diagram for Ana15:[legend]List of all members.

Public Member Functions
	Ana15 (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~Ana15 ()
virtual StatusCode	initialize ()
virtual StatusCode	execute ()
virtual StatusCode	finalize ()
virtual StatusCode	sysExecute ()
void	put (IDataProviderSvc *svc, DataObject *object, const std::string &address, const bool useRootInTES=true) const
void	put (DataObject *object, const std::string &address, const bool useRootInTES=true) const
Gaudi::Utils::GetData< TYPE >::return_type	get (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
Gaudi::Utils::GetData< TYPE >::return_type	get (const std::string &location, const bool useRootInTES=true) const
TYPE *	getDet (IDataProviderSvc *svc, const std::string &location) const
TYPE *	getDet (const std::string &location) const
bool	exist (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
bool	exist (const std::string &location, const bool useRootInTES=true) const
bool	existDet (IDataProviderSvc *svc, const std::string &location) const
bool	existDet (const std::string &location) const
TYPE *	getOrCreate (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
TYPE *	getOrCreate (const std::string &location, const bool useRootInTES=true) const
bool	registerContext () const
INTupleSvc *	evtColSvc () const
IAlgContextSvc *	contextSvc () const
TOOL *	tool (const std::string &type, const std::string &name, const IInterface *parent=0, bool create=true) const
TOOL *	tool (const std::string &type, const IInterface *parent=0, bool create=true) const
SERVICE *	svc (const std::string &name, const bool create=true) const
IUpdateManagerSvc *	updMgrSvc () const
IDataProviderSvc *	fastContainersSvc () const
StatusCode	Error (const std::string &msg, const StatusCode st=StatusCode::FAILURE, const size_t mx=10) const
StatusCode	Warning (const std::string &msg, const StatusCode st=StatusCode::FAILURE, const size_t mx=10) const
StatusCode	Print (const std::string &msg, const StatusCode st=StatusCode::SUCCESS, const MSG::Level lev=MSG::INFO) const
StatusCode	Assert (const bool ok, const std::string &message="", const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Assert (const bool ok, const char *message, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Exception (const std::string &msg, const GaudiException &exc, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Exception (const std::string &msg, const std::exception &exc, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
StatusCode	Exception (const std::string &msg="no message", const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
MsgStream &	msgStream (const MSG::Level level) const
MsgStream &	always () const
MsgStream &	fatal () const
MsgStream &	err () const
MsgStream &	error () const
MsgStream &	warning () const
MsgStream &	info () const
MsgStream &	debug () const
MsgStream &	verbose () const
MsgStream &	msg () const
const Statistics &	counters () const
StatEntity &	counter (const std::string &tag) const
MSG::Level	msgLevel () const
bool	msgLevel (const MSG::Level level) const
void	resetMsgStream () const
bool	typePrint () const
bool	propsPrint () const
bool	statPrint () const
bool	errorsPrint () const
long	printStat (const MSG::Level level=MSG::ALWAYS) const
long	printErrors (const MSG::Level level=MSG::ALWAYS) const
long	printProps (const MSG::Level level=MSG::ALWAYS) const
void	registerCondition (const std::string &condition, StatusCode(CallerClass::*mf)()=NULL)
void	registerCondition (const std::string &condition, CondType *&condPtrDest, StatusCode(CallerClass::*mf)()=NULL)
void	registerCondition (char *condition, StatusCode(CallerClass::*mf)()=NULL)
void	registerCondition (TargetClass *condition, StatusCode(CallerClass::*mf)()=NULL)
StatusCode	runUpdate ()
TransientFastContainer< T > *	getFastContainer (const std::string &location, typename TransientFastContainer< T >::size_type initial=0)
StatusCode	release (const IInterface *interface) const
unsigned long	release ()
const std::string &	context () const
const std::string &	rootInTES () const
double	globalTimeOffset () const
virtual StatusCode	sysStart ()
virtual StatusCode	sysInitialize ()
virtual StatusCode	sysReinitialize ()
virtual StatusCode	sysRestart ()
virtual StatusCode	sysStop ()
virtual StatusCode	sysFinalize ()
virtual StatusCode	sysBeginRun ()
virtual StatusCode	sysEndRun ()
virtual const std::string &	name () const
virtual const std::string &	version () const
virtual StatusCode	configure ()
virtual StatusCode	terminate ()
virtual StatusCode	start ()
virtual StatusCode	stop ()
virtual StatusCode	reinitialize ()
virtual StatusCode	restart ()
virtual bool	isExecuted () const
virtual void	setExecuted (bool state)
virtual void	resetExecuted ()
virtual StatusCode	beginRun ()
virtual StatusCode	endRun ()
virtual Gaudi::StateMachine::State	FSMState () const
virtual Gaudi::StateMachine::State	targetFSMState () const
virtual bool	isEnabled () const
virtual bool	filterPassed () const
virtual void	setFilterPassed (bool state)
StatusCode	service (const std::string &name, T *&psvc, bool createIf=true) const
StatusCode	service (const std::string &svcType, const std::string &svcName, T *&psvc) const
void	setOutputLevel (int level)
IAuditorSvc *	auditorSvc () const
IChronoStatSvc *	chronoSvc () const
IChronoStatSvc *	chronoStatService () const
IDataProviderSvc *	detSvc () const
IDataProviderSvc *	detDataService () const
IConversionSvc *	detCnvSvc () const
IConversionSvc *	detDataCnvService () const
IDataProviderSvc *	eventSvc () const
IDataProviderSvc *	evtSvc () const
IDataProviderSvc *	eventDataService () const
IConversionSvc *	eventCnvSvc () const
IConversionSvc *	eventDataCnvService () const
IHistogramSvc *	histoSvc () const
IHistogramSvc *	histogramDataService () const
IMessageSvc *	msgSvc () const
IMessageSvc *	messageService () const
INTupleSvc *	ntupleSvc () const
INTupleSvc *	ntupleService () const
IRndmGenSvc *	randSvc () const
IToolSvc *	toolSvc () const
IExceptionSvc *	exceptionSvc () const
ISvcLocator *	serviceLocator () const
ISvcLocator *	svcLoc () const
StatusCode	createSubAlgorithm (const std::string &type, const std::string &name, Algorithm *&pSubAlg)
std::vector< Algorithm * > *	subAlgorithms () const
virtual StatusCode	setProperty (const Property &p)
virtual StatusCode	setProperty (const std::string &s)
virtual StatusCode	setProperty (const std::string &n, const std::string &v)
StatusCode	setProperty (const std::string &name, const TYPE &value)
virtual StatusCode	getProperty (Property *p) const
virtual const Property &	getProperty (const std::string &name) const
virtual StatusCode	getProperty (const std::string &n, std::string &v) const
virtual const std::vector< Property * > &	getProperties () const
StatusCode	setProperties ()
Property *	declareProperty (const std::string &name, T &property, const std::string &doc="none") const
Property *	declareRemoteProperty (const std::string &name, IProperty *rsvc, const std::string &rname="") const
IMonitorSvc *	monitorSvc () const
void	declareInfo (const std::string &name, const T &var, const std::string &desc) const
void	declareInfo (const std::string &name, const std::string &format, const void *var, int size, const std::string &desc) const
unsigned long	addRef ()
StatusCode	queryInterface (const InterfaceID &riid, void **)
Static Public Member Functions
static const InterfaceID &	interfaceID ()
static const InterfaceID &	interfaceID ()
static const InterfaceID &	interfaceID ()
Public Attributes
	SUCCESS
	NO_INTERFACE
	VERSMISMATCH
	LAST_ERROR
Protected Types
typedef std::map< std::string, StatEntity >	Statistics
typedef std::map< std::string, unsigned int >	Counter
typedef std::vector< IAlgTool * >	AlgTools
typedef std::pair< IInterface *, std::string >	ServiceEntry
typedef std::vector< ServiceEntry >	Services
Protected Member Functions
StatusCode	releaseTool (const IAlgTool *tool) const
StatusCode	releaseSvc (const IInterface *svc) const
bool	isInitialized () const
bool	isFinalized () const
int	outputLevel () const
IntegerProperty &	outputLevelProperty ()
void	initOutputLevel (Property &prop)
Static Protected Attributes
static const bool	IgnoreRootInTES
static const bool	UseRootInTES
Private Attributes
std::string	m_TopStage
std::string	m_path
ValidationTree *	m_valiTree
	create a root tree

---

Detailed Description

Ana15 -- an example to show how to use Fifteen package.

Mar. 30, 2009 created by Zhe Wang

Definition at line 19 of file Ana15.h.

---

Constructor & Destructor Documentation

Ana15::Ana15	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 61 of file Ana15.cc.

                                                           :
  GaudiAlgorithm(name, pSvcLocator)
{
  declareProperty("TopStage",m_TopStage="SingleLoader","Name of top stage");  
}

Ana15::~Ana15

(

)

[virtual]

Definition at line 69 of file Ana15.cc.

{
}

---

Member Function Documentation

StatusCode Ana15::initialize

(

)

[virtual]

Reimplemented from GaudiAlgorithm.

Definition at line 73 of file Ana15.cc.

{
  this->GaudiAlgorithm::initialize();

  // init path
  if(m_TopStage=="SingleLoader") m_path = ReadoutHeader::defaultLocation();
  if(m_TopStage=="TrigRead")     m_path = SimReadoutHeader::defaultLocation();
  if(m_TopStage=="Electronic")   m_path = ElecHeader::defaultLocation();
  if(m_TopStage=="Detector")     m_path = SimHeader::defaultLocation();
  if(m_TopStage=="Kinematic")    m_path = GenHeader::defaultLocation();
  
  // init root tree
  m_valiTree=new ValidationTree;
  m_valiTree->create();

  StatusCode status;

  return StatusCode::SUCCESS;
}

StatusCode Ana15::execute

(

)

[virtual]

Reimplemented from GaudiAlgorithm.

Definition at line 93 of file Ana15.cc.

{
  double m_e=5.109989e-1;  // MeV
  double m_alpha=3727.379; // MeV

  m_valiTree->reset();

  StatusCode   sc;
  DataObject   *pObject;
  sc=eventSvc()->retrieveObject(m_path,pObject);     // get a new data
  if(sc.isFailure()) {
    error()<<"Failed to read data"<<endreq;
    return sc;
  }

  const HeaderObject     *pHeader=0;
  pHeader=dynamic_cast<HeaderObject*>(pObject);

  info()<<"to grep: new data arrived at time: "<<pHeader->earliest()<<endreq;

  const GenHeader        *pGenHeader=0;
  const SimHeader        *pSimHeader=0;
  const ElecHeader       *pElecHeader=0;
  const SimTrigHeader    *pSimTrigHeader=0;
  const SimReadoutHeader *pSimReadoutHeader=0;
  const ReadoutHeader    *pReadoutHeader=0;

  const IHeader* firstHeader=0;

  // get the header object pointer of the top stage
  if(m_TopStage=="SingleLoader") {
    // ReadoutHeader
    pReadoutHeader=dynamic_cast<ReadoutHeader*>(pObject);
  } else if(m_TopStage=="TrigRead") { 
    // SimReadoutHeader
    pSimReadoutHeader=dynamic_cast<SimReadoutHeader*>(pObject);
  } else if(m_TopStage=="Electronic") {
    // ElecHeader
    pElecHeader=dynamic_cast<ElecHeader*>(pObject);
  } else if (m_TopStage=="Detector") {
    // SimHeader
    pSimHeader=dynamic_cast<SimHeader*>(pObject);
  } else if (m_TopStage=="Kinematic") {
    // GenHeader
    pGenHeader=dynamic_cast<GenHeader*>(pObject);
  }

  // get the header object pointer of it lower stages
  // only the first one, for low event rate, it is enough
  if(pReadoutHeader!=0) {
    // SimReadoutHeader
    if(pReadoutHeader->inputHeaders().size() != 0) {
      firstHeader=*((pReadoutHeader->inputHeaders()).begin());
      pSimReadoutHeader= dynamic_cast<const SimReadoutHeader*>(firstHeader);
    }
  }
  if(pSimReadoutHeader!=0) {
    debug() << "Number of input SimTrigHeaders:" << pSimReadoutHeader->inputHeaders().size() << endreq;
    // SimTrigHeader
    if(pSimReadoutHeader->inputHeaders().size() != 0) {
      firstHeader=*((pSimReadoutHeader->inputHeaders()).begin());
      pSimTrigHeader= dynamic_cast<const SimTrigHeader*>(firstHeader);
    }
  }
  if(pSimTrigHeader!=0) {
    // ElecHeader
    if(pSimTrigHeader->inputHeaders().size() != 0) {
      firstHeader=*((pSimTrigHeader->inputHeaders()).begin());
      pElecHeader= dynamic_cast<const ElecHeader*>(firstHeader);
    }
  }
  if(pElecHeader!=0) {
    // SimHeader
    if(pElecHeader->inputHeaders().size() != 0) {
      firstHeader=*((pElecHeader->inputHeaders()).begin());
      pSimHeader= dynamic_cast<const SimHeader*>(firstHeader);
    }
  }
  if(pSimHeader!=0) {
    // GenHeader
    if(pSimHeader->inputHeaders().size() != 0) {
      firstHeader=*((pSimHeader->inputHeaders()).begin());
      pGenHeader= dynamic_cast<const GenHeader*>(firstHeader);
    }
  }


  // event information
  m_valiTree->time=pHeader->earliest().GetSeconds();

  // check ReadoutHeader information
  if(pReadoutHeader!=0) {
    m_valiTree->execSing=pReadoutHeader->execNumber();
    
    // adc tdc
    const Readout* pReadout = pReadoutHeader->readout();
    Detector det = pReadout->detector();
    
    if(det.detectorId() == DetectorId::kAD1 ||  // AD information
       det.detectorId() == DetectorId::kAD2 ||
       det.detectorId() == DetectorId::kAD3 ||
       det.detectorId() == DetectorId::kAD4) {
      
      // convert it to PMT Crate Readout
      const ReadoutPmtCrate* pmtcrate = 0;
      pmtcrate = dynamic_cast<const ReadoutPmtCrate*>(pReadout);
      
      // channnel map
      const ReadoutPmtCrate::PmtChannelReadouts& channelMap
          = pmtcrate->channelReadout();
      
      ReadoutPmtCrate::PmtChannelReadouts::const_iterator cmci;

      for(cmci=channelMap.begin();cmci!=channelMap.end();++cmci) {
        const FeeChannelId&      channelId = cmci->first;
        const ReadoutPmtChannel& aChannel  = cmci->second;

        const vector<int>&          tdc=aChannel.tdc();
        const vector<int>&          adc=aChannel.adc();

        int raw_tdc=tdc[0];
        int raw_adc=adc[0];
        
        m_valiTree->adc_sum+=raw_adc;

        verbose()<<channelId<<" tdc "<<raw_tdc<<" adc "<<raw_adc<<endreq;
      }
    }  // end of AD crate
  } // end of pReadoutHeader

  // check SimReadoutHeader information
  if(pSimReadoutHeader!=0) {
    m_valiTree->execTrig=pSimReadoutHeader->execNumber();

    const SimReadoutHeader::SimReadoutContainer 
      simReadout=pSimReadoutHeader->readouts();
    debug()<<"No. of SimReadout: "<<simReadout.size()<<endreq;
    debug()<<pSimReadoutHeader->execNumber()<<endreq;
  } // end of pSimReadoutHeader

  // check electrontic simulation information if available
  if(pElecHeader!=0) {
    m_valiTree->execElec=pElecHeader->execNumber();
  } // end of pElecHeader

  // check trigger header information if available
  if(pSimTrigHeader!=0) {
    debug()<<pSimTrigHeader->commandHeader()->collections().size();
  }

  // check detector simulation information if available
  if(pSimHeader!=0) {
    m_valiTree->execDets=pSimHeader->execNumber();
    const SimHitHeader* simhitheader = pSimHeader->hits();
    const SimHitHeader::hc_map& hcmap = simhitheader->hitCollection();
    SimHitHeader::hc_map::const_iterator it;

    m_valiTree->SimHitsEntries=0;
    int index(0);

    debug()<<"size of hit collection "<< hcmap.size()<<endreq;
    for (it=hcmap.begin(); it != hcmap.end(); ++it) {

      Detector det(it->first);
      debug() << "Got hit collection from " << det.detName()<<" id= "
              << det.siteDetPackedData()<<endreq;

      const std::vector<SimHit*>& hitvec=it->second->collection();
      std::vector<SimHit*>::const_iterator it_hvec;
      
      debug() <<"size of hit vector "<< hitvec.size()<<endreq;
      for (it_hvec=hitvec.begin(); it_hvec!=hitvec.end(); ++it_hvec) {
        
        index=m_valiTree->SimHitsEntries;
        if(index<MaxSimHitEntries){

          m_valiTree->hitTime[index]  = (*it_hvec)->hitTime();
          m_valiTree->hitx[index]     = (*it_hvec)->localPos().x();
          m_valiTree->hity[index]     = (*it_hvec)->localPos().y();
          m_valiTree->hitz[index]     = (*it_hvec)->localPos().z();
          m_valiTree->sensID[index]   = (*it_hvec)->sensDetId();
          m_valiTree->weight[index]   = (*it_hvec)->weight();
          
          // optical photon's ancestor's pdg
          if((*it_hvec)->ancestor().track()) {
            m_valiTree->ancestorPdg[index] = (*it_hvec)->ancestor().track()->particle();
          }
          else {
            m_valiTree->ancestorPdg[index] = 0;
          }
          //optical photon's wavelength
          const SimPmtHit* pmthit = static_cast<const SimPmtHit*>(*it_hvec);
          m_valiTree->wavelength[index]   = (*pmthit).wavelength();
          m_valiTree->polx[index]   = (*pmthit).pol().x();
          m_valiTree->poly[index]   = (*pmthit).pol().y();
          m_valiTree->polz[index]   = (*pmthit).pol().z();
          m_valiTree->px[index]   = (*pmthit).dir().x();
          m_valiTree->py[index]   = (*pmthit).dir().y();
          m_valiTree->pz[index]   = (*pmthit).dir().z();
          
          m_valiTree->SimHitsEntries = index+1;
          
        }
        else {
          debug()<< " Reached the max hit limit!!!!!" <<endl;
        }
      }
    }
  } // end of pSimHeader

  //  if(0==1) {
  if(pGenHeader!=0) {
    m_valiTree->execKine=pGenHeader->execNumber();
    debug()<<"Generator name: "<<pGenHeader->generatorName()<<endreq;
    debug()<<"execNum: "<<pGenHeader->execNumber()<<endreq;

    const HepMC::GenEvent* event=pGenHeader->event();
    //  event->print();

    HepMC::GenEvent::particle_const_iterator p;
    int count=0;
    for ( p = event->particles_begin(); p!=event->particles_end(); ++p ) {      
      (*p)->print();
      if((*p)->production_vertex()) {   // only outgoing particle has a production vertex
        count++;

        if(count==1) {
          
          m_valiTree->trk1_pdgId=(*p)->pdg_id();
          
          m_valiTree->trk1_px=(*p)->momentum().px();
          m_valiTree->trk1_py=(*p)->momentum().py();
          m_valiTree->trk1_pz=(*p)->momentum().pz();
          m_valiTree->trk1_e =(*p)->momentum().e();
          
          m_valiTree->trk1_x=(*p)->production_vertex()->position().x();
          m_valiTree->trk1_y=(*p)->production_vertex()->position().y();
          m_valiTree->trk1_z=(*p)->production_vertex()->position().z();
          m_valiTree->trk1_t=(*p)->production_vertex()->position().t();

          debug()<<"trk1_pdgId: "<<m_valiTree->trk1_pdgId<<endreq;
          
          if(m_valiTree->trk1_pdgId==22) { // gamma
            m_valiTree->trk1_ke=m_valiTree->trk1_e;
            m_valiTree->trk1_ve=m_valiTree->trk1_e;
          }else if (m_valiTree->trk1_pdgId==11) { // e-
            m_valiTree->trk1_ke=m_valiTree->trk1_e-m_e;
            m_valiTree->trk1_ve=m_valiTree->trk1_e-m_e;
          }else if (m_valiTree->trk1_pdgId==-11) { // e+
            m_valiTree->trk1_ke=m_valiTree->trk1_e-m_e;
            m_valiTree->trk1_ve=m_valiTree->trk1_e+m_e;
          }else if (m_valiTree->trk1_pdgId==1000020040) { // alpha
            m_valiTree->trk1_ke=m_valiTree->trk1_e-m_alpha;
            m_valiTree->trk1_ve=m_valiTree->trk1_e-m_alpha;
          }
        }  // outgoing particle 1
        if(count==2) {
          m_valiTree->trk2_pdgId=(*p)->pdg_id();
          
          m_valiTree->trk2_px=(*p)->momentum().px();
          m_valiTree->trk2_py=(*p)->momentum().py();
          m_valiTree->trk2_pz=(*p)->momentum().pz();
          m_valiTree->trk2_e =(*p)->momentum().e();
          
          m_valiTree->trk2_x=(*p)->production_vertex()->position().x();
          m_valiTree->trk2_y=(*p)->production_vertex()->position().y();
          m_valiTree->trk2_z=(*p)->production_vertex()->position().z();
          m_valiTree->trk2_t=(*p)->production_vertex()->position().t();
          if(m_valiTree->trk2_pdgId==22) { // gamma
            m_valiTree->trk2_ke=m_valiTree->trk2_e;
            m_valiTree->trk2_ve=m_valiTree->trk2_e;
          }else if (m_valiTree->trk2_pdgId==11) { // e-
            m_valiTree->trk2_ke=m_valiTree->trk2_e-m_e;
            m_valiTree->trk2_ve=m_valiTree->trk2_e-m_e;
          }else if (m_valiTree->trk2_pdgId==-11) { // e+
            m_valiTree->trk2_ke=m_valiTree->trk2_e-m_e;
            m_valiTree->trk2_ve=m_valiTree->trk2_e+m_e;
          }else if (m_valiTree->trk2_pdgId==1000020040) { // alpha
            m_valiTree->trk2_ke=m_valiTree->trk2_e-m_alpha;
            m_valiTree->trk2_ve=m_valiTree->trk2_e-m_alpha;
          }
        }  // outgoing particle 2
      } // end of outgoing particles
    }
  } // end of pGenHeader

  m_valiTree->fill();
  return StatusCode::SUCCESS;
}

StatusCode Ana15::finalize

(

)

[virtual]

Reimplemented from GaudiAlgorithm.

Definition at line 392 of file Ana15.cc.

{
  m_valiTree->close();
  delete m_valiTree;

  return this->GaudiAlgorithm::finalize();
}

---

Member Data Documentation

std::string Ana15::m_TopStage [private]

Definition at line 32 of file Ana15.h.

std::string Ana15::m_path [private]

Definition at line 33 of file Ana15.h.

ValidationTree* Ana15::m_valiTree [private]

create a root tree

Definition at line 36 of file Ana15.h.

---

The documentation for this class was generated from the following files:

• /home/dayabay/installs/trunk_2011_04_11_opt/NuWa-trunk/dybgaudi/Tutorial/Sim15/src/Ana15.h
• /home/dayabay/installs/trunk_2011_04_11_opt/NuWa-trunk/dybgaudi/Tutorial/Sim15/src/Ana15.cc

---