EvtRecVeeVertexCnv Class Reference

#include <EvtRecVeeVertexCnv.h>

Inheritance diagram for EvtRecVeeVertexCnv:


Public Types
typedef Ty1	source
typedef Ty2	destination
Public Member Functions
virtual	~EvtRecVeeVertexCnv ()
virtual long	repSvcType () const
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	declareObject (const std::string &fullPath, const CLID &clid, const std::string &treename, const std::string &branchname)
	Store TDS path to link a particular converter to an object on the TDS.
virtual StatusCode	createObj (IOpaqueAddress addr, DataObject &dat)
	Convert the persistent object to transient.
virtual StatusCode	createRep (DataObject pObject, IOpaqueAddress &refpAddress)
	Convert the transient object to the requested representation.
virtual StatusCode	fillRepRefs (IOpaqueAddress pAddress, DataObject pObject)
	Resolve the references of the converted object.
virtual StatusCode	fillObjRefs (IOpaqueAddress pAddress, DataObject pObject)
	Resolve the references of the converted object.
TObject *	getReadObject () const
	get the object to be read
destination *	operator (const source &) const
Static Public Member Functions
static const CLID &	classID ()
static const unsigned char	storageType ()
	Storage type and class ID.
Protected Member Functions
	EvtRecVeeVertexCnv (ISvcLocator *svc)
virtual StatusCode	DataObjectToTObject (DataObject obj, RootAddress addr)
	transformation to root
virtual StatusCode	TObjectToDataObject (DataObject *&obj)
	transformation from root
virtual destination *	convert (const source &) const =0
Protected Attributes
RootCnvSvc *	m_cnvSvc
std::vector< RootCnvSvc::Leaf >	m_leaves
RootInterface *	m_rootInterface
	pointer to the RootInterface
IDataProviderSvc *	m_eds
	pointer to eventdataservice
int	m_branchNr
	the branchNr of this converter for writing
int	m_branchNrDst
int	m_branchNrMc
int	m_branchNrRecon
int	m_branchNrEvtRec
int	m_branchNrEvtHeader
int	m_branchNrEvtNavigator
TObject *	m_objRead
	the object that was read
CLID	CLID_top
	the CLID of the upper converter if any
TArrayS *	m_branchNumbers
	array with number of branches for reading
std::string	m_rootBranchname
	root branchname (may be concatenated of severals)
std::string	m_rootTreename
	each converter knows it's treename
std::string	m_currentFileName
std::vector< void * >	m_adresses
	each converter knows the corresponding adresses
RootEvtSelector *	m_evtsel
Private Attributes
commonData	m_common
	relational maps
TObjArray *	m_evtRecVeeVertexCol
Friends
class	CnvFactory< EvtRecVeeVertexCnv >

Detailed Description

Definition at line 14 of file EvtRecVeeVertexCnv.h.

Member Typedef Documentation

template<class Ty1, class Ty2>

typedef Ty2 Converter< Ty1, Ty2 >::destination [inherited]

Definition at line 19 of file Converter.h.

template<class Ty1, class Ty2>

typedef Ty1 Converter< Ty1, Ty2 >::source [inherited]

Definition at line 18 of file Converter.h.

Constructor & Destructor Documentation

virtual EvtRecVeeVertexCnv::~EvtRecVeeVertexCnv ( ) [inline, virtual]

Definition at line 21 of file EvtRecVeeVertexCnv.h.

00021 {}

EvtRecVeeVertexCnv::EvtRecVeeVertexCnv ( ISvcLocator * svc ) [protected]

Definition at line 20 of file EvtRecVeeVertexCnv.cxx.

References RootEventBaseCnv::m_adresses, m_evtRecVeeVertexCol, and RootEventBaseCnv::m_rootBranchname.

00021    : RootEventBaseCnv(classID(), svc)
00022 {
00023    m_rootBranchname = "m_evtRecVeeVertexCol";
00024    m_adresses.push_back(&m_evtRecVeeVertexCol);
00025    m_evtRecVeeVertexCol = 0;
00026 }

Member Function Documentation

static const CLID& EvtRecVeeVertexCnv::classID ( ) [inline, static]

Definition at line 19 of file EvtRecVeeVertexCnv.h.

References CLID_EvtRecVeeVertexCol.

Referenced by RootCnvSvc::addConverters().

00019 { return CLID_EvtRecVeeVertexCol; }

template<class Ty1, class Ty2>

virtual destination* Converter< Ty1, Ty2 >::convert ( const source & ) const [protected, pure virtual, inherited]

StatusCode RootEventBaseCnv::createObj	(	IOpaqueAddress *	addr,
		DataObject *&	dat
	)			`[virtual, inherited]`

Convert the persistent object to transient.

Reimplemented in EventCnv.

Definition at line 157 of file RootEventBaseCnv.cxx.

References RootEventBaseCnv::CLID_top, Bes_Common::DEBUG, calibUtil::ERROR, RootInterface::getBranchEntry(), RootAddress::getBranchname(), RootInterface::getCurrentFileName(), RootInterface::getENDFILE(), RootInterface::getEntries(), RootAddress::getEntryNr(), RootAddress::getNrBranches(), RootAddress::getPath(), RootEventBaseCnv::getReadObject(), RootEvtSelector::getRecId(), RootInterface::getSelectFromTag(), RootAddress::getTreename(), Bes_Common::INFO, RootEventBaseCnv::m_adresses, RootEventBaseCnv::m_branchNumbers, RootEventBaseCnv::m_currentFileName, RootEventBaseCnv::m_eds, RootEventBaseCnv::m_evtsel, RootEventBaseCnv::m_objRead, RootEventBaseCnv::m_rootInterface, msgSvc(), RootInterface::setBranchAddress(), RootAddress::setEntryNr(), subSeperate::temp, and RootEventBaseCnv::TObjectToDataObject().

00158                                                            {
00159   // transform ROOT object to TDS object
00160   MsgStream log(msgSvc(), "RootEventBaseCnv");
00161   log << MSG::DEBUG << "RootEventBaseCnv::createObj with clid " <<addr->clID()<< endreq;
00162   StatusCode sc;
00163 
00164   // add 2005-11-29
00165   // log<<MSG::INFO<<"######### RootEventBaseCnv ::createObj begin of createObj:  m_branchNumbers "<<m_branchNumbers->GetSize()<<"###############"<<endreq;
00166 
00167   RootAddress *raddr=dynamic_cast<RootAddress *>(addr);
00168   if (!raddr) {
00169     log << MSG::ERROR << "Could not downcast to Root address" << endreq;
00170     return StatusCode::FAILURE;
00171   }
00172  
00173   static int temp =0;   //control the begin of each files 2005-12-01
00174   static int entryN =0; //control the event number of each files 2005-21-01
00175   static int brN =0; //control munber of branch of the tree;
00176   int lastBrn = brN;
00177   //lastBrn = brN; 
00178   static int branchN=0;
00179   static bool isSet=true;
00180 
00181   static int entryBefore = 0;
00182   static bool addEntryEachFile = true;
00183  
00184  
00185    if(m_rootInterface->getENDFILE() || (temp >0 && temp < branchN)){ // if the file has get the end:y the go to next file to create a new tree
00186     
00187     if(m_rootInterface->getENDFILE() ) {
00188       entryN = 0;
00189     }
00190     
00191     temp++;
00192     
00193     delete m_branchNumbers;
00194     m_branchNumbers = new TArrayS(0);
00195    
00196     if(temp == branchN) {
00197       temp =0;
00198     }
00199    }
00200 
00201    if(m_rootInterface->getENDFILE()) addEntryEachFile = true;
00202    // the 2nd method 
00203    if(m_evtsel->getRecId() - entryBefore == 0) { // first event in this file
00204      delete m_branchNumbers;
00205      m_branchNumbers = new TArrayS(0);
00206    }
00207 
00208    //new method to initialize the branchNumber 
00209    if(m_currentFileName=="") m_currentFileName = m_rootInterface->getCurrentFileName();
00210    if(!(m_currentFileName == m_rootInterface->getCurrentFileName())){ 
00211            m_currentFileName = m_rootInterface->getCurrentFileName();
00212            delete m_branchNumbers;
00213            m_branchNumbers = new TArrayS(0);
00214    }   
00215    //---------------------------------------- 
00216    
00217   if (m_branchNumbers->GetSize()<=0) {
00218     if(isSet) brN++;
00219     int branchNumber;
00220     for (int nb=0;nb<raddr->getNrBranches();nb++) {
00221       sc=m_rootInterface->setBranchAddress(raddr->getTreename().c_str(),raddr->getBranchname(nb).c_str(),m_adresses[nb],branchNumber);
00222       if (!sc.isSuccess()) 
00223       {
00224         if(isSet) brN--;    //liangyt: if fail to retrieve this branch, this will be not a effective branch.
00225         //entryN++;   //liangyt: this is the second method
00226         if(temp>0) temp--;  //temp > 0 means recording effective branch number.
00227         return sc;
00228       }
00229       m_branchNumbers->Set(nb+1);
00230       m_branchNumbers->AddAt(branchNumber,nb);
00231      
00232     }
00233   }
00234 
00236   if(addEntryEachFile&&(m_evtsel->getRecId()>entryBefore)){   // for a new file, add entry for ONLY one time.
00237     entryBefore += m_rootInterface->getEntries();
00238     addEntryEachFile = false;
00239   }
00240 
00241   if(lastBrn == brN && isSet ){
00242     branchN = brN;
00243     isSet=false;    
00244   }
00245 
00246   if(isSet==false) log << MSG::INFO <<" 1st method set event as : "<<int(entryN/branchN)<<endreq;
00247   if(isSet==false) raddr->setEntryNr(int(entryN/branchN));//former method, keep it to be backup.
00248   if(m_evtsel) log << MSG::INFO <<" event id = "<<m_evtsel->getRecId()<<endreq;
00249 
00250 
00251   int eventID = 0;
00252   if(entryBefore > m_evtsel->getRecId())
00253     eventID = m_evtsel->getRecId() + m_rootInterface->getEntries() - entryBefore;
00254   else if(entryBefore == m_evtsel->getRecId()) eventID = 0;
00255   else log << MSG::ERROR <<"eventId error!!!"<<endreq;
00256   
00257   log << MSG::INFO <<" 2nd  method set event as : "<<eventID<<endreq;
00258   
00259   if(m_evtsel) raddr->setEntryNr(eventID);
00260 
00261   //add dengzy
00262   if(m_rootInterface->getSelectFromTag()!=0)
00263   {
00264     if(m_evtsel)
00265       raddr->setEntryNr( m_evtsel->getRecId() );
00266   }//end of add by dengzy
00267 
00268   // read branch
00269   
00270   if (m_branchNumbers->GetSize()>0) {
00271     int nbtot=0,nb;
00272     for (int ib=0;ib<m_branchNumbers->GetSize();ib++) {
00273       //sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),nb);
00274       //change to get branch entry with addr(set address for each entry) liangyt
00275       sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),m_adresses[ib],nb);      
00276       if (sc.isFailure()) {
00277         log << MSG::ERROR << "Could not read branch " << raddr->getBranchname(nb) << endreq;
00278         return sc;
00279       }
00280       nbtot+=nb;
00281     }
00282   }
00283      
00284   else {   // get ROOT object
00285     if (CLID_top) {
00286       IConverter *p=conversionSvc()->converter(CLID_top);
00287       RootEventBaseCnv *cnv=dynamic_cast<RootEventBaseCnv *>(p);
00288       if (!cnv) {
00289         log << MSG::ERROR << "Could not downcast to RootEventBaseCnv "  << endreq;
00290         return StatusCode::FAILURE;
00291       }
00292       m_objRead=cnv->getReadObject();
00293     }
00294   }
00295 
00296   //do concrete transformation in derived converter
00297   sc = TObjectToDataObject(refpObject);
00298   if (sc.isFailure()) {
00299     log << MSG::ERROR << "Could not transform object" << endreq;
00300     return sc;
00301   }
00302 
00303   // verify if we have to register
00304   IRegistry* ent = addr->registry();
00305   if ( ent == 0)   {
00306        sc=m_eds->registerObject(raddr->getPath(),refpObject);
00307        if (sc.isFailure()) {
00308          log << MSG::ERROR << "Could not register object " << raddr->getPath()<<" status "<<sc.getCode()<<endreq;
00309        }
00310        //    }
00311   }
00312 
00313   entryN++;
00314   return StatusCode::SUCCESS;
00315 }

StatusCode RootEventBaseCnv::createRep	(	DataObject *	pObject,
		IOpaqueAddress *&	refpAddress
	)			`[virtual, inherited]`

Convert the transient object to the requested representation.

Definition at line 78 of file RootEventBaseCnv.cxx.

References RootCnvSvc::createAddress(), RootEventBaseCnv::DataObjectToTObject(), calibUtil::ERROR, RootEventBaseCnv::m_cnvSvc, and msgSvc().

00079                                                            {
00080   // Purpose and Method: Convert the transient object to ROOT
00081 
00082   MsgStream log(msgSvc(), "RootEventBaseCnv");
00083  
00084   StatusCode sc= StatusCode::SUCCESS;
00085   // get the corresponding address 
00086   RootAddress *rootaddr;
00087   sc=m_cnvSvc->createAddress(obj,addr);
00088 
00089   rootaddr = dynamic_cast<RootAddress *>(addr);
00090 
00091   if (sc.isFailure() || !rootaddr  ) { 
00092     log << MSG::ERROR << "Could not create address for clid " <<obj->clID()<<", objname "<<obj->name()<<endreq;
00093     return StatusCode::FAILURE;
00094   }
00095 
00096   //  do the real conversion in the derived converter
00097   sc = DataObjectToTObject(obj,rootaddr);
00098 
00099   delete addr;
00100   addr = NULL;
00101 
00102   if (sc.isFailure()) {
00103     log << MSG::ERROR << "Could not transform object" << endreq;
00104     return sc;
00105   }
00106 
00107   return StatusCode::SUCCESS;
00108 }

StatusCode EvtRecVeeVertexCnv::DataObjectToTObject	(	DataObject *	obj,
		RootAddress *	addr
	)			`[protected, virtual]`

transformation to root

Implements RootEventBaseCnv.

Definition at line 108 of file EvtRecVeeVertexCnv.cxx.

00108                                                                                          {
00109   MsgStream log(msgSvc(), "EvtRecVeeVertexCnv");
00110   log << MSG::DEBUG << "EvtRecVeeVertexCnv::DataObjectToTObject" << endreq;
00111 
00112   EvtRecVeeVertexCol* evtRecVeeVertexCol = dynamic_cast<EvtRecVeeVertexCol*> (obj);
00113   if ( ! evtRecVeeVertexCol ) {
00114     log << MSG::ERROR << "Could not downcast to EvtRecVeeVertexCol" << endreq;
00115     return StatusCode::FAILURE;
00116   }
00117 
00118   DataObject* evt;
00119   m_eds->findObject(EventModel::EvtRec::Event, evt);
00120   if ( evt == NULL ) {
00121      log << MSG::ERROR << "Could not get EvtRecObject in TDS" << endreq;
00122      return StatusCode::FAILURE;
00123   }
00124   EvtRecObject* devtTds = dynamic_cast<EvtRecObject*> (evt);
00125   if ( ! devtTds ) {
00126     log << MSG::ERROR << "EvtRecVeeVertexCnv: Could not downcast to TDS EvtRecObject" << endreq;
00127   }
00128 
00129   IOpaqueAddress *addr;
00130   m_cnvSvc->getEvtRecCnv()->createRep(evt, addr);
00131   TEvtRecObject* recEvt = m_cnvSvc->getEvtRecCnv()->getWriteObject();
00132 
00133   const TObjArray* m_evtRecVeeVertexCol = recEvt->getEvtRecVeeVertexCol();
00134   if ( ! m_evtRecVeeVertexCol ) return StatusCode::SUCCESS;
00135 
00136   IDataProviderSvc* dataSvc = 0;
00137   StatusCode sc = serviceLocator()->getService("EventDataSvc",
00138         IDataProviderSvc::interfaceID(), (IInterface*&)dataSvc);
00139   if ( sc.isFailure() ) {
00140      log << MSG::FATAL << "Could not get EventDataSvc in EvtRecVeeVertexCnv" << endreq;
00141      return sc;
00142   }
00143 
00144   EvtRecTrackCol::iterator    evtRecTrackColbegin, evtRecTrackColend;
00145 
00146   SmartDataPtr<EvtRecTrackCol> evtRecTrackCol(dataSvc, EventModel::EvtRec::EvtRecTrackCol);
00147   if ( ! evtRecTrackCol) {
00148      log << MSG::INFO << "Could not find EvtRecTrackCol" << endreq;
00149   }
00150   else {
00151      evtRecTrackColbegin = evtRecTrackCol->begin();
00152      evtRecTrackColend = evtRecTrackCol->end();
00153   }
00154 
00155   // convert
00156   recEvt->clearEvtRecVeeVertexCol();
00157   EvtRecVeeVertexCol::const_iterator evtRecVeeVertex = evtRecVeeVertexCol->begin();
00158 
00159   for ( ; evtRecVeeVertex != evtRecVeeVertexCol->end(); evtRecVeeVertex++) {
00160      EvtRecVeeVertex* ptr = *evtRecVeeVertex;
00161      TEvtRecVeeVertex* evtRecVeeVertexRoot = new TEvtRecVeeVertex();
00162      // set ...
00163      evtRecVeeVertexRoot->setVertexId( ptr->vertexId() );
00164      evtRecVeeVertexRoot->setVertexType( ptr->vertexType() );
00165      evtRecVeeVertexRoot->setChi2( ptr->chi2() );
00166      evtRecVeeVertexRoot->setNdof( ptr->ndof() );
00167      evtRecVeeVertexRoot->setMass( ptr->mass() );
00168 
00169      double wTrackParameter[7];
00170      for (int i = 0; i < 7; i++) {
00171        wTrackParameter[i] = ptr->w()[i];
00172      }
00173      evtRecVeeVertexRoot->setW( wTrackParameter );
00174      double EwTrackParameter[28];
00175      int i = 0;
00176      for (int j = 0; j < 7; j++) {
00177       for (int k = j; k < 7; k++) {
00178         EwTrackParameter[i] = ptr->Ew()[j][k];
00179         i++;
00180        }
00181      }
00182      assert(i == 28);
00183      evtRecVeeVertexRoot->setEw( EwTrackParameter );
00184      int pairId[2];
00185      pairId[0] = ptr->pair(0);
00186      pairId[1] = ptr->pair(1);
00187      evtRecVeeVertexRoot->setPair( pairId );
00188      evtRecVeeVertexRoot->setNCharge( ptr->nCharge() );
00189      evtRecVeeVertexRoot->setNTracks( ptr->nTracks() );
00190      
00191      int daughterId[2];
00192      if ( ptr->daughter(0).target() != NULL && ptr->daughter(1).target() != NULL ) {
00193        assert(evtRecTrackColbegin != evtRecTrackColend);
00194        EvtRecTrackCol::iterator it0 = find(evtRecTrackColbegin, 
00195                                          evtRecTrackColend, ptr->daughter(0).target());
00196        EvtRecTrackCol::iterator it1 = find(evtRecTrackColbegin, 
00197                                          evtRecTrackColend, ptr->daughter(1).target());
00198        assert(it0 != evtRecTrackColend);
00199        assert(it1 != evtRecTrackColend);
00200        daughterId[0] = it0 - evtRecTrackColbegin;
00201        daughterId[1] = it1 - evtRecTrackColbegin;
00202        evtRecVeeVertexRoot->setDaughter( daughterId );
00203      }
00204      recEvt->addEvtRecVeeVertex(evtRecVeeVertexRoot);
00205   }
00206 
00207   return StatusCode::SUCCESS;
00208 }

void RootEventBaseCnv::declareObject	(	const std::string &	fullPath,
		const CLID &	clid,
		const std::string &	treename,
		const std::string &	branchname
	)			`[inherited]`

Store TDS path to link a particular converter to an object on the TDS.

Definition at line 150 of file RootEventBaseCnv.cxx.

References RootEventBaseCnv::m_leaves.

Referenced by RecMucRecHitCnv::RecMucRecHitCnv().

00151                                                                                       {
00152     // Purpose and Method:  Save the path on the TDS, treename, pathname in the m_leaves vector, 
00153     //   corresponding to the DataObject that the converter handles.
00154    m_leaves.push_back(RootCnvSvc::Leaf(path, cl,treename,branchname));
00155 }

StatusCode RootEventBaseCnv::fillObjRefs	(	IOpaqueAddress *	pAddress,
		DataObject *	pObject
	)			`[virtual, inherited]`

Resolve the references of the converted object.

Definition at line 117 of file RootEventBaseCnv.cxx.

References msgSvc().

00118                                                  {
00119     // Purpose and Method:  Resolve the references of the converted object.
00120     //     It is expected that derived classes will override this method.
00121     MsgStream log(msgSvc(), "RootEventBaseCnv");
00122     return StatusCode::SUCCESS;
00123 }

StatusCode RootEventBaseCnv::fillRepRefs	(	IOpaqueAddress *	pAddress,
		DataObject *	pObject
	)			`[virtual, inherited]`

Resolve the references of the converted object.

Definition at line 109 of file RootEventBaseCnv.cxx.

References msgSvc().

00110                                                  {
00111     // Purpose and Method:  Resolve the references of the converted object.
00112     //     It is expected that derived classes will override this method.
00113     MsgStream log(msgSvc(), "RootEventBaseCnv");
00114     return StatusCode::SUCCESS;
00115 }

StatusCode RootEventBaseCnv::finalize ( ) [virtual, inherited]

Definition at line 142 of file RootEventBaseCnv.cxx.

References RootEventBaseCnv::m_cnvSvc.

00142                                         {
00143     if ( m_cnvSvc )     {
00144         m_cnvSvc->release();
00145         m_cnvSvc=0;
00146     }
00147     return Converter::finalize();
00148 }

TObject* RootEventBaseCnv::getReadObject ( ) const [inline, inherited]

get the object to be read

Definition at line 125 of file RootEventBaseCnv.h.

References RootEventBaseCnv::m_objRead.

Referenced by RootEventBaseCnv::createObj().

00125 { return m_objRead;}

StatusCode RootEventBaseCnv::initialize ( ) [virtual, inherited]

Reimplemented in DigiCnv, DstCnv, EvtHeaderCnv, EvtNavigatorCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, and TrigCnv.

Definition at line 126 of file RootEventBaseCnv.cxx.

References IID_IRootCnvSvc, and RootEventBaseCnv::m_cnvSvc.

Referenced by TrigCnv::initialize(), RecTrackCnv::initialize(), McCnv::initialize(), HltCnv::initialize(), EvtRecCnv::initialize(), EvtNavigatorCnv::initialize(), EvtHeaderCnv::initialize(), DstCnv::initialize(), and DigiCnv::initialize().

00126                                           {
00127  
00128   StatusCode status = Converter::initialize();
00129 
00130   if ( status.isSuccess() )   {
00131     IService* isvc = 0;
00132     status = serviceLocator()->service("RootCnvSvc", isvc, false);
00133     if ( !status.isSuccess() )   status = serviceLocator()->service("EventCnvSvc", isvc, true);
00134     if ( status.isSuccess() )   {
00135       status = isvc->queryInterface(IID_IRootCnvSvc, (void**)&m_cnvSvc);
00136     }
00137   }
00138   
00139   return status;
00140 }

template<class Ty1, class Ty2>

destination* Converter< Ty1, Ty2 >::operator ( const source & ) const [inline, inherited]

virtual long RootEventBaseCnv::repSvcType ( ) const [inline, virtual, inherited]

Definition at line 88 of file RootEventBaseCnv.h.

References ROOT_StorageType.

00088                                     { 
00089       return ROOT_StorageType; 
00090     }

static const unsigned char RootEventBaseCnv::storageType ( ) [inline, static, inherited]

Storage type and class ID.

Definition at line 84 of file RootEventBaseCnv.h.

References ROOT_StorageType.

00084                                                   {
00085       return ROOT_StorageType;
00086     }

StatusCode EvtRecVeeVertexCnv::TObjectToDataObject ( DataObject *& obj ) [protected, virtual]

transformation from root

Implements RootEventBaseCnv.

Definition at line 28 of file EvtRecVeeVertexCnv.cxx.

00028                                                                           {
00029   // creation of TDS object from root object
00030   MsgStream log(msgSvc(), "EvtRecVeeVertexCnv");
00031   log << MSG::DEBUG << "EvtRecVeeVertexCnv::TObjectToDataObject" << endreq;
00032 
00033   // create the TDS location for the EvtRecVeeVertex Collection
00034   EvtRecVeeVertexCol* evtRecVeeVertexCol = new EvtRecVeeVertexCol;
00035   refpObject = evtRecVeeVertexCol;
00036 
00037   if ( ! m_evtRecVeeVertexCol ) return StatusCode::SUCCESS;
00038 
00039   IDataProviderSvc* dataSvc = 0;
00040   StatusCode sc = serviceLocator()->getService("EventDataSvc",
00041         IDataProviderSvc::interfaceID(), (IInterface*&)dataSvc);
00042   if ( sc.isFailure() ) {
00043      log << MSG::FATAL << "Could not get EventDataSvc in EvtRecVeeVertexCnv" << endreq;
00044      return sc;
00045   }
00046 
00047   SmartDataPtr<EvtRecTrackCol> evtRecTrackCol(dataSvc, EventModel::EvtRec::EvtRecTrackCol);
00048   if ( ! evtRecTrackCol) {
00049      log << MSG::INFO << "Could not find EvtRecTrackCol" << endreq;
00050   }
00051 
00052   // convert
00053   TIter evtRecVeeVertexIter(m_evtRecVeeVertexCol);
00054   TEvtRecVeeVertex* evtRecVeeVertexRoot = 0;
00055   while ( (evtRecVeeVertexRoot = (TEvtRecVeeVertex*)evtRecVeeVertexIter.Next() ) ) {
00056      EvtRecVeeVertex* evtRecVeeVertex = new EvtRecVeeVertex();
00057      m_common.m_rootEvtRecVeeVertexMap[evtRecVeeVertexRoot] = evtRecVeeVertex;
00058      // set ....
00059      evtRecVeeVertex->setVertexId( evtRecVeeVertexRoot->vertexId() );
00060      evtRecVeeVertex->setVertexType( evtRecVeeVertexRoot->vertexType() );
00061      evtRecVeeVertex->setChi2( evtRecVeeVertexRoot->chi2() );
00062      evtRecVeeVertex->setNdof( evtRecVeeVertexRoot->ndof() );
00063      evtRecVeeVertex->setMass( evtRecVeeVertexRoot->mass() );
00064 
00065      HepVector wTrackParameter(7, 0);
00066      for (int i = 0; i < 7; i++) {
00067        wTrackParameter[i] = evtRecVeeVertexRoot->w(i);
00068      }
00069      evtRecVeeVertex->setW( wTrackParameter );
00070 
00071      HepSymMatrix EwTrackParameter(7, 0);
00072      int i = 0;
00073      for (int j = 0; j < 7; j++) {
00074       for (int k = j; k < 7; k++) {
00075         EwTrackParameter[j][k] = EwTrackParameter[k][j] = evtRecVeeVertexRoot->Ew(i);
00076         i++;
00077        }
00078      }
00079      evtRecVeeVertex->setEw( EwTrackParameter );
00080 
00081      std::pair<int, int> pairId;
00082      pairId.first = evtRecVeeVertexRoot->pair(0);
00083      pairId.second = evtRecVeeVertexRoot->pair(1);
00084      evtRecVeeVertex->setPair( pairId );
00085      evtRecVeeVertex->setNCharge( evtRecVeeVertexRoot->nCharge() );
00086      evtRecVeeVertex->setNTracks( evtRecVeeVertexRoot->nTracks() );
00087      int evtRecTrackId0 = evtRecVeeVertexRoot->daughter(0);
00088      if ( evtRecTrackId0 >= 0) {
00089         evtRecVeeVertex->addDaughter(
00090             dynamic_cast<EvtRecTrack*>(evtRecTrackCol->containedObject(evtRecTrackId0)), 0
00091               );
00092      }
00093      int evtRecTrackId1 = evtRecVeeVertexRoot->daughter(1);
00094      if ( evtRecTrackId1 >= 0) {
00095         evtRecVeeVertex->addDaughter(
00096             dynamic_cast<EvtRecTrack*>(evtRecTrackCol->containedObject(evtRecTrackId1)), 1
00097               );
00098      }
00099 
00100      evtRecVeeVertexCol->push_back(evtRecVeeVertex);
00101   }
00102 
00103   delete m_evtRecVeeVertexCol;
00104   m_evtRecVeeVertexCol = 0;
00105   return StatusCode::SUCCESS;
00106 }