EmcMcHitCnv Class Reference

#include <EmcMcHitCnv.h>

Inheritance diagram for EmcMcHitCnv:


Public Types
typedef Ty1	source
typedef Ty2	destination
Public Member Functions
virtual	~EmcMcHitCnv ()
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
	EmcMcHitCnv (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_emcMcHitCol
	root object to be read
Friends
class	CnvFactory< EmcMcHitCnv >

Detailed Description

Definition at line 19 of file EmcMcHitCnv.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 EmcMcHitCnv::~EmcMcHitCnv ( ) [inline, virtual]

Definition at line 29 of file EmcMcHitCnv.h.

00029 { };

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

Definition at line 28 of file EmcMcHitCnv.cxx.

References RootEventBaseCnv::m_adresses, m_emcMcHitCol, RootEventBaseCnv::m_rootBranchname, and msgSvc().

00029 : RootEventBaseCnv(classID(), svc)
00030 {
00031     // Here we associate this converter with the /Event path on the TDS.
00032     MsgStream log(msgSvc(), "EmcMcHitCnv");
00033     // log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
00034     m_rootBranchname ="m_emcMcHitCol";
00035     //declareObject(EventModel::MC::EmcMcHitCol, objType(), m_rootTreename, m_rootBranchname);
00036     m_adresses.push_back(&m_emcMcHitCol);
00037     m_emcMcHitCol=0;
00038 }

Member Function Documentation

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

Definition at line 24 of file EmcMcHitCnv.h.

References CLID_EmcMcHitCol.

Referenced by RootCnvSvc::addConverters().

00025   {
00026        return CLID_EmcMcHitCol; 
00027   }

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 EmcMcHitCnv::DataObjectToTObject	(	DataObject *	obj,
		RootAddress *	addr
	)			`[protected, virtual]`

transformation to root

Implements RootEventBaseCnv.

Definition at line 132 of file EmcMcHitCnv.cxx.

References TMcEvent::addEmcMc(), TMcEvent::clearEmcMcHitCol(), RootEventBaseCnv::createRep(), Bes_Common::DEBUG, EventModel::MC::EmcMcHitCol, calibUtil::ERROR, EventModel::MC::Event, TMcEvent::getEmcMcHitCol(), RootCnvSvc::getMcCnv(), McCnv::getWriteObject(), RootEventBaseCnv::m_cnvSvc, RootEventBaseCnv::m_eds, m_emcMcHitCol, msgSvc(), TEmcMc::setDepositEnergy(), TEmcMc::setHitEmc(), TEmcMc::setHitMap(), TEmcMc::setId(), TEmcMc::setPDGCharge(), TEmcMc::setPDGCode(), TEmcMc::setPositionX(), TEmcMc::setPositionY(), TEmcMc::setPositionZ(), TEmcMc::setPx(), TEmcMc::setPy(), TEmcMc::setPz(), TEmcMc::setTime(), and TEmcMc::setTrackIndex().

00132                                                                                  {
00133 
00134     MsgStream log(msgSvc(), "EmcMcHitCnv");
00135     log << MSG::DEBUG << "EmcMcHitCnv::DataObjectToTObject" << endreq;
00136     StatusCode sc=StatusCode::SUCCESS;
00137 
00138     EmcMcHitCol * emcMcHitCnvTds=dynamic_cast<EmcMcHitCol *> (obj);
00139     if (!emcMcHitCnvTds) {
00140         log << MSG::ERROR << "Could not downcast to EmcMcHitCol" << endreq;
00141         return StatusCode::FAILURE;
00142     }
00143 
00144     DataObject *evt;
00145     m_eds->findObject(EventModel::MC::Event,evt);
00146     if (evt==NULL) {
00147         log << MSG::ERROR << "Could not get McEvent in TDS "  << endreq;
00148         return StatusCode::FAILURE;
00149     }
00150     McEvent * devtTds=dynamic_cast<McEvent *> (evt);
00151     if (!devtTds) {
00152         log << MSG::ERROR << "EmcMcHitCnv:Could not downcast to TDS McEvent" << endreq;
00153     }
00154     IOpaqueAddress *addr;
00155 
00156     m_cnvSvc->getMcCnv()->createRep(evt,addr); 
00157     TMcEvent *McEvt=m_cnvSvc->getMcCnv()->getWriteObject();
00158 
00159     const TObjArray *m_emcMcHitCol = McEvt->getEmcMcHitCol();
00160     if (!m_emcMcHitCol) return sc;
00161     McEvt->clearEmcMcHitCol(); //necessary in case there is I/O at the same time since array is static
00162     EmcMcHitCol::const_iterator emcMcTds;
00163 
00164     for (emcMcTds = emcMcHitCnvTds->begin(); emcMcTds != emcMcHitCnvTds->end(); emcMcTds++) {
00165         UInt_t id  = ((*emcMcTds)->identify()).get_value() ;
00166         UInt_t trackIndex = (*emcMcTds) ->getTrackIndex();
00167 
00168         int   hitEmc = (*emcMcTds) ->getHitEmc();
00169         int   PDGCode = (*emcMcTds) ->getPDGCode();
00170         double PDGCharge = (*emcMcTds) ->getPDGCharge();
00171         double time = (*emcMcTds) ->getTime();
00172         std::map<Identifier, Double_t>  TdshitMap = (*emcMcTds)->getHitMap();
00173 
00174         map<Int_t, Double_t>  HitMap;
00175         map<Identifier, Double_t>::iterator iHitMap;
00176         for(iHitMap=TdshitMap.begin();
00177                 iHitMap!=TdshitMap.end();
00178                 iHitMap++) {
00179 
00180             Int_t id(iHitMap->first.get_value());
00181             pair<Int_t, Double_t> Hit(id, iHitMap->second);
00182             HitMap.insert(Hit);
00183 
00184             //RecEmcHitCol::iterator iHit;
00185             //for(iHit=emcRecHitCol->begin();
00186             //        iHit!=emcRecHitCol->end();
00187             //        iHit++) {
00188 
00189             //    RecEmcID idHit((*iHit)->getCellId());
00190 
00191             //    if(id==idHit) {
00192             //        RecEmcFraction frac(*(*iHit));
00193             //        frac.Fraction(iHitMap->second);
00194             //        emcShowerTds->Insert(frac);
00195             //        break;
00196             //    }
00197             //} //RecEmcHitCol
00198         } //CellIdMap
00199 
00200         Double_t xPosition =   (*emcMcTds) ->getPositionX() ;
00201         Double_t yPosition =   (*emcMcTds) ->getPositionY() ;
00202         Double_t zPosition =   (*emcMcTds) ->getPositionZ() ;
00203         Double_t px=   (*emcMcTds) ->getPx() ;
00204         Double_t py =   (*emcMcTds) ->getPy() ;
00205         Double_t pz =   (*emcMcTds) ->getPz() ;
00206         Double_t depositEnergy = (*emcMcTds) ->getDepositEnergy() ;
00207 
00208         TEmcMc *emcMcRoot = new TEmcMc();
00209         //m_common.m_emcMcHitMap[(*emcMcTds)] = emcMcRoot;
00210 
00211         emcMcRoot->setId(id);
00212         emcMcRoot->setHitMap(HitMap);
00213 
00214         emcMcRoot->setHitEmc(hitEmc);
00215         emcMcRoot->setPDGCode(PDGCode);
00216         emcMcRoot->setPDGCharge(PDGCharge);
00217         emcMcRoot->setTime(time);
00218         //emcMcRoot->setHitMap(hitMap);
00219 
00220         emcMcRoot->setTrackIndex(trackIndex);
00221         emcMcRoot->setPositionX(xPosition);
00222         emcMcRoot->setPositionY(yPosition);
00223         emcMcRoot->setPositionZ(zPosition);
00224         emcMcRoot->setPx(px);
00225         emcMcRoot->setPy(py);
00226         emcMcRoot->setPz(pz);
00227         emcMcRoot->setDepositEnergy(depositEnergy);
00228 
00229         McEvt->addEmcMc(emcMcRoot);
00230     }
00231 
00232     return StatusCode::SUCCESS;
00233 }

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 EmcMcHitCnv::TObjectToDataObject ( DataObject *& obj ) [protected, virtual]

transformation from root

Implements RootEventBaseCnv.

Definition at line 40 of file EmcMcHitCnv.cxx.

00040                                                                    {
00041   // creation of TDS object from root object
00042 
00043     MsgStream log(msgSvc(), "EmcMcHitCnv");
00044     log << MSG::DEBUG << "EmcMcHitCnv::TObjectToDataObject" << endreq;
00045     StatusCode sc=StatusCode::SUCCESS;
00046 
00047     // create the TDS location for the MdcMc Collection
00048     EmcMcHitCol* emcMcTdsCol = new EmcMcHitCol;
00049     refpObject=emcMcTdsCol;
00050 
00051 
00052     // now convert
00053     if (!m_emcMcHitCol) return sc;
00054     TIter emcMcIter(m_emcMcHitCol);
00055     TEmcMc *emcMcRoot = 0;
00056     while ((emcMcRoot = (TEmcMc*)emcMcIter.Next())) {
00057       // unsigned int id = emcMcRoot ->getId();
00058       Identifier id ;
00059          unsigned int id_root = emcMcRoot ->getId();
00060          id = id_root;
00061                 unsigned int trackIndex = emcMcRoot ->getTrackIndex();
00062 
00063     int   hitEmc = emcMcRoot ->getHitEmc();
00064     int   PDGCode = emcMcRoot ->getPDGCode();
00065     double PDGCharge = emcMcRoot ->getPDGCharge();
00066     double time = emcMcRoot ->getTime();
00067 
00068                 double xPosition =   emcMcRoot ->getPositionX() ;
00069                 double yPosition =   emcMcRoot ->getPositionY() ;
00070                 double zPosition =   emcMcRoot ->getPositionZ() ;
00071                 double px=   emcMcRoot ->getPx() ;
00072     double py =   emcMcRoot ->getPy() ;
00073     double pz =   emcMcRoot ->getPz() ;
00074     double depositEnergy = emcMcRoot ->getDepositEnergy() ;
00075 
00076     EmcMcHit *emcMcTds = new EmcMcHit();
00077     m_common.m_rootEmcMcHitMap[emcMcRoot] = emcMcTds;
00078 
00079     emcMcTds->setIdentifier(id);
00080 
00081     map<Int_t, Double_t>  hitMap = emcMcRoot ->getHitMap();
00082     map<Identifier, Double_t>  TdshitMap;
00083     map<Int_t, Double_t> HitMap = emcMcRoot->getHitMap();
00084     map<Int_t, Double_t>::iterator iHitMap;
00085     for(iHitMap=HitMap.begin();
00086             iHitMap!=HitMap.end();
00087             iHitMap++) {
00088 
00089         Identifier id(iHitMap->first);
00090         pair<Identifier, Double_t> TdsHit(id, iHitMap->second);
00091         TdshitMap.insert(TdsHit);
00092 
00093         //RecEmcHitCol::iterator iHit;
00094         //for(iHit=emcRecHitCol->begin();
00095         //        iHit!=emcRecHitCol->end();
00096         //        iHit++) {
00097 
00098         //    RecEmcID idHit((*iHit)->getCellId());
00099 
00100         //    if(id==idHit) {
00101         //        RecEmcFraction frac(*(*iHit));
00102         //        frac.Fraction(iHitMap->second);
00103         //        emcShowerTds->Insert(frac);
00104         //        break;
00105         //    }
00106         //} //RecEmcHitCol
00107     } //CellIdMap
00108     emcMcTds->setHitMap(TdshitMap);
00109 
00110     emcMcTds->setHitEmc(hitEmc);
00111     emcMcTds->setPDGCode(PDGCode);
00112     emcMcTds->setPDGCharge(PDGCharge);
00113     emcMcTds->setTime(time);
00114 
00115     emcMcTds->setTrackIndex(trackIndex);
00116     emcMcTds->setPositionX(xPosition);
00117     emcMcTds->setPositionY(yPosition);
00118     emcMcTds->setPositionZ(zPosition);
00119     emcMcTds->setPx(px);
00120     emcMcTds->setPy(py);
00121     emcMcTds->setPz(pz);
00122     emcMcTds->setDepositEnergy(depositEnergy);
00123 
00124     emcMcTdsCol->push_back(emcMcTds);
00125     }
00126     //m_emcMcHitCol->Delete();  // wensp add 2005/12/30
00127     delete m_emcMcHitCol;
00128     m_emcMcHitCol = 0;
00129     return StatusCode::SUCCESS;
00130 }