#include <RecMdcTrackCnv.h>
Inheritance diagram for RecMdcTrackCnv:
Public Types | |
typedef Ty2 | destination |
typedef Ty2 | destination |
typedef Ty1 | source |
typedef Ty1 | source |
Public Member Functions | |
virtual StatusCode | createObj (IOpaqueAddress *addr, DataObject *&dat) |
Convert the persistent object to transient. | |
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 | createRep (DataObject *pObject, IOpaqueAddress *&refpAddress) |
Convert the transient object to the requested representation. | |
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. | |
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 | fillObjRefs (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. | |
virtual StatusCode | fillRepRefs (IOpaqueAddress *pAddress, DataObject *pObject) |
Resolve the references of the converted object. | |
virtual StatusCode | fillRepRefs (IOpaqueAddress *pAddress, DataObject *pObject) |
Resolve the references of the converted object. | |
virtual StatusCode | finalize () |
virtual StatusCode | finalize () |
TObject * | getReadObject () const |
get the object to be read | |
TObject * | getReadObject () const |
get the object to be read | |
virtual StatusCode | initialize () |
virtual StatusCode | initialize () |
destination * | operator (const source &) const |
destination * | operator (const source &) const |
virtual long | repSvcType () const |
virtual long | repSvcType () const |
virtual | ~RecMdcTrackCnv () |
virtual | ~RecMdcTrackCnv () |
Static Public Member Functions | |
const CLID & | classID () |
const CLID & | classID () |
const unsigned char | storageType () |
Storage type and class ID. | |
const unsigned char | storageType () |
Storage type and class ID. | |
Protected Member Functions | |
virtual destination * | convert (const source &) const =0 |
virtual destination * | convert (const source &) const =0 |
virtual StatusCode | DataObjectToTObject (DataObject *obj, RootAddress *addr) |
transformation to root | |
virtual StatusCode | DataObjectToTObject (DataObject *obj, RootAddress *addr) |
transformation to root | |
RecMdcTrackCnv (ISvcLocator *svc) | |
RecMdcTrackCnv (ISvcLocator *svc) | |
virtual StatusCode | TObjectToDataObject (DataObject *&obj) |
transformation from root | |
virtual StatusCode | TObjectToDataObject (DataObject *&obj) |
transformation from root | |
Protected Attributes | |
CLID | CLID_top |
the CLID of the upper converter if any | |
std::vector< void * > | m_adresses |
each converter knows the corresponding adresses | |
std::vector< void * > | m_adresses |
each converter knows the corresponding adresses | |
int | m_branchNr |
the branchNr of this converter for writing | |
int | m_branchNrDst |
int | m_branchNrEvtHeader |
int | m_branchNrEvtRec |
int | m_branchNrMc |
int | m_branchNrRecon |
TArrayS * | m_branchNumbers |
array with number of branches for reading | |
TArrayS * | m_branchNumbers |
array with number of branches for reading | |
RootCnvSvc * | m_cnvSvc |
RootCnvSvc * | m_cnvSvc |
std::string | m_currentFileName |
IDataProviderSvc * | m_eds |
pointer to eventdataservice | |
IDataProviderSvc * | m_eds |
pointer to eventdataservice | |
RootEvtSelector * | m_evtsel |
RootEvtSelector * | m_evtsel |
std::vector< RootCnvSvc::Leaf > | m_leaves |
std::vector< RootCnvSvc::Leaf > | m_leaves |
TObject * | m_objRead |
the object that was read | |
TObject * | m_objRead |
the object that was read | |
std::string | m_rootBranchname |
root branchname (may be concatenated of severals) | |
RootInterface * | m_rootInterface |
pointer to the RootInterface | |
RootInterface * | m_rootInterface |
pointer to the RootInterface | |
std::string | m_rootTreename |
each converter knows it's treename | |
Private Attributes | |
commonData | m_common |
relational maps | |
TObjArray * | m_recMdcTrackCol |
root object to be read | |
TObjArray * | m_recMdcTrackCol |
root object to be read | |
Friends | |
class | CnvFactory<RecMdcTrackCnv> |
|
|
|
|
|
|
|
|
|
00025 { };
|
|
00034 : RootEventBaseCnv(classID(), svc) 00035 { 00036 // Here we associate this converter with the /Event path on the TDS. 00037 MsgStream log(msgSvc(), "RecMdcTrackCnv"); 00038 //m_rootTreename ="Rec"; 00039 m_rootBranchname ="m_recMdcTrackCol"; 00040 //declareObject(EventModel::Recon::RecMdcTrackCol, objType(), m_rootTreename, m_rootBranchname); 00041 m_adresses.push_back(&m_recMdcTrackCol); 00042 m_recMdcTrackCol=0; 00043 }
|
|
00025 { };
|
|
|
|
00021 {
00022 return CLID_RecMdcTrackCol;
00023 }
|
|
00021 {
00022 return CLID_RecMdcTrackCol;
00023 }
|
|
|
|
|
|
Convert the persistent object to transient.
|
|
Convert the persistent object to transient.
Reimplemented in EventCnv, and EventCnv. 00157 { 00158 // transform ROOT object to TDS object 00159 MsgStream log(msgSvc(), "RootEventBaseCnv"); 00160 log << MSG::DEBUG << "RootEventBaseCnv::createObj with clid " <<addr->clID()<< endreq; 00161 StatusCode sc; 00162 00163 // add 2005-11-29 00164 // log<<MSG::INFO<<"######### RootEventBaseCnv ::createObj begin of createObj: m_branchNumbers "<<m_branchNumbers->GetSize()<<"###############"<<endreq; 00165 00166 RootAddress *raddr=dynamic_cast<RootAddress *>(addr); 00167 if (!raddr) { 00168 log << MSG::ERROR << "Could not downcast to Root address" << endreq; 00169 return StatusCode::FAILURE; 00170 } 00171 00172 static int temp =0; //control the begin of each files 2005-12-01 00173 static int entryN =0; //control the event number of each files 2005-21-01 00174 static int brN =0; //control munber of branch of the tree; 00175 int lastBrn = brN; 00176 //lastBrn = brN; 00177 static int branchN=0; 00178 static bool isSet=true; 00179 00180 static int entryBefore = 0; 00181 static bool addEntryEachFile = true; 00182 00183 00184 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 00185 00186 if(m_rootInterface->getENDFILE() ) { 00187 entryN = 0; 00188 } 00189 00190 temp++; 00191 00192 delete m_branchNumbers; 00193 m_branchNumbers = new TArrayS(0); 00194 00195 if(temp == branchN) { 00196 temp =0; 00197 } 00198 } 00199 00200 if(m_rootInterface->getENDFILE()) addEntryEachFile = true; 00201 // the 2nd method 00202 if(m_evtsel->getRecId() - entryBefore == 0) { // first event in this file 00203 delete m_branchNumbers; 00204 m_branchNumbers = new TArrayS(0); 00205 } 00206 00207 //new method to initialize the branchNumber 00208 if(m_currentFileName=="") m_currentFileName = m_rootInterface->getCurrentFileName(); 00209 if(!(m_currentFileName == m_rootInterface->getCurrentFileName())){ 00210 m_currentFileName = m_rootInterface->getCurrentFileName(); 00211 delete m_branchNumbers; 00212 m_branchNumbers = new TArrayS(0); 00213 } 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 int eventID = 0; 00251 if(entryBefore > m_evtsel->getRecId()) 00252 eventID = m_evtsel->getRecId() + m_rootInterface->getEntries() - entryBefore; 00253 else if(entryBefore == m_evtsel->getRecId()) eventID = 0; 00254 else log << MSG::ERROR <<"eventId error!!!"<<endreq; 00255 00256 log << MSG::INFO <<" 2nd method set event as : "<<eventID<<endreq; 00257 if(m_evtsel) raddr->setEntryNr(eventID); 00258 // read branch 00259 00260 if (m_branchNumbers->GetSize()>0) { 00261 int nbtot=0,nb; 00262 for (int ib=0;ib<m_branchNumbers->GetSize();ib++) { 00263 //sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),nb); 00264 //change to get branch entry with addr(set address for each entry) liangyt 00265 sc=m_rootInterface->getBranchEntry(m_branchNumbers->At(ib),raddr->getEntryNr(),m_adresses[ib],nb); 00266 if (sc.isFailure()) { 00267 log << MSG::ERROR << "Could not read branch " << raddr->getBranchname(nb) << endreq; 00268 return sc; 00269 } 00270 nbtot+=nb; 00271 } 00272 } 00273 00274 else { // get ROOT object 00275 if (CLID_top) { 00276 IConverter *p=conversionSvc()->converter(CLID_top); 00277 RootEventBaseCnv *cnv=dynamic_cast<RootEventBaseCnv *>(p); 00278 if (!cnv) { 00279 log << MSG::ERROR << "Could not downcast to RootEventBaseCnv " << endreq; 00280 return StatusCode::FAILURE; 00281 } 00282 m_objRead=cnv->getReadObject(); 00283 } 00284 } 00285 00286 //do concrete transformation in derived converter 00287 sc = TObjectToDataObject(refpObject); 00288 if (sc.isFailure()) { 00289 log << MSG::ERROR << "Could not transform object" << endreq; 00290 return sc; 00291 } 00292 00293 // verify if we have to register 00294 IRegistry* ent = addr->registry(); 00295 if ( ent == 0) { 00296 sc=m_eds->registerObject(raddr->getPath(),refpObject); 00297 if (sc.isFailure()) { 00298 log << MSG::ERROR << "Could not register object " << raddr->getPath()<<" status "<<sc.getCode()<<endreq; 00299 } 00300 // } 00301 } 00302 00303 entryN++; 00304 return StatusCode::SUCCESS; 00305 }
|
|
Convert the transient object to the requested representation.
|
|
Convert the transient object to the requested representation.
00078 { 00079 // Purpose and Method: Convert the transient object to ROOT 00080 00081 MsgStream log(msgSvc(), "RootEventBaseCnv"); 00082 00083 StatusCode sc= StatusCode::SUCCESS; 00084 // get the corresponding address 00085 RootAddress *rootaddr; 00086 sc=m_cnvSvc->createAddress(obj,addr); 00087 00088 rootaddr = dynamic_cast<RootAddress *>(addr); 00089 00090 if (sc.isFailure() || !rootaddr ) { 00091 log << MSG::ERROR << "Could not create address for clid " <<obj->clID()<<", objname "<<obj->name()<<endreq; 00092 return StatusCode::FAILURE; 00093 } 00094 00095 // do the real conversion in the derived converter 00096 sc = DataObjectToTObject(obj,rootaddr); 00097 00098 delete addr; 00099 addr = NULL; 00100 00101 if (sc.isFailure()) { 00102 log << MSG::ERROR << "Could not transform object" << endreq; 00103 return sc; 00104 } 00105 00106 return StatusCode::SUCCESS; 00107 }
|
|
transformation to root
Implements RootEventBaseCnv. |
|
transformation to root
Implements RootEventBaseCnv. 00184 { 00185 00186 MsgStream log(msgSvc(), "RecMdcTrackCnv"); 00187 log << MSG::DEBUG << "RecMdcTrackCnv::DataObjectToTObject" << endreq; 00188 00189 RecMdcTrackCol * mdcTrackColTds=dynamic_cast<RecMdcTrackCol *> (obj); 00190 if (!mdcTrackColTds) { 00191 log << MSG::ERROR << "Could not downcast to MdcTrackCol" << endreq; 00192 return StatusCode::FAILURE; 00193 } 00194 DataObject *evt; 00195 StatusCode sc = m_eds->findObject(EventModel::Recon::Event,evt); 00196 if (!sc.isSuccess()) { 00197 log << MSG::ERROR << "Could not get ReconEvent in TDS " << endreq; 00198 return StatusCode::FAILURE; 00199 } 00200 ReconEvent * devtTds=dynamic_cast<ReconEvent *> (evt); 00201 if (!devtTds) { 00202 log << MSG::ERROR << "RecMdcTrackCnv:Could not downcast to TDS ReconEvent" << endreq; 00203 } 00204 00205 IOpaqueAddress *addr; 00206 00207 m_cnvSvc->getRecTrackCnv()->createRep(evt,addr); 00208 TRecTrackEvent *recEvt=m_cnvSvc->getRecTrackCnv()->getWriteObject(); 00209 00210 const TObjArray *m_recMdcTrackCol = recEvt->getRecMdcTrackCol(); 00211 if (!m_recMdcTrackCol) return StatusCode::SUCCESS; 00212 //necessary in case there is I/O at the same time since array is static 00213 recEvt->clearRecMdcTrackCol(); 00214 00215 RecMdcTrackCol::const_iterator recMdcTrack; 00216 for (recMdcTrack = mdcTrackColTds->begin(); recMdcTrack != mdcTrackColTds->end(); recMdcTrack++) { 00217 Int_t trackId =(*recMdcTrack)->trackId(); 00218 Double_t helix[5]; 00219 for(int i=0;i<5;i++){ 00220 log<<MSG::INFO<<" recMdcTrack.helix("<<i<<"): "<<(*recMdcTrack)->helix(i)<<endreq; 00221 helix[i] = (*recMdcTrack)->helix(i); 00222 } 00223 00224 Int_t stat = (*recMdcTrack)->stat(); 00225 Double_t chi = (*recMdcTrack)->chi2(); 00226 Int_t ndof = (*recMdcTrack)->ndof(); 00227 Int_t ns = (*recMdcTrack)->nster(); 00228 00229 Double_t er[15]; 00230 for (Int_t i=0; i<15; i++){ 00231 er[i] = (*recMdcTrack)->err(i); 00232 } 00233 Int_t nh = (*recMdcTrack)->getNhits(); 00234 // vector<Int_t> vecHits = (*recMdcTrack)->vecHits(); 00235 Double_t vx0 = (*recMdcTrack)->getVX0(); 00236 Double_t vy0 = (*recMdcTrack)->getVY0(); 00237 Double_t vz0 = (*recMdcTrack)->getVZ0(); 00238 Double_t fiterm =(*recMdcTrack)->getFiTerm(); 00239 00240 TRecMdcTrack *recMdcTrackRoot = new TRecMdcTrack(); 00241 //m_common.m_recMdcTrackMap[(*recMdcTrack)] = recMdcTrackRoot; 00242 00243 recMdcTrackRoot->setTrackId( trackId ); 00244 recMdcTrackRoot->setHelix(helix); 00245 log<<MSG::INFO<<" test,recMdcTrackRoot.px: "<<recMdcTrackRoot->px() 00246 <<" recMdcTrackRoot.py: "<<recMdcTrackRoot->py() 00247 <<" recMdcTrackRoot.pz: "<<recMdcTrackRoot->pz() 00248 <<endreq; 00249 00250 recMdcTrackRoot->setStat( stat ); 00251 recMdcTrackRoot->setChi2( chi ); 00252 recMdcTrackRoot->setNdof( ndof ); 00253 recMdcTrackRoot->setErr( er ); 00254 recMdcTrackRoot->setNhits( nh ); 00255 recMdcTrackRoot->setNster( ns ); 00256 // recMdcTrackRoot->setVecHits(vecHits); 00257 recMdcTrackRoot->setVX0( vx0 ); 00258 recMdcTrackRoot->setVY0( vy0 ); 00259 recMdcTrackRoot->setVZ0( vz0 ); 00260 recMdcTrackRoot->setFiTerm( fiterm ); 00261 00262 recEvt->addRecMdcTrack(recMdcTrackRoot); 00263 } 00264 00265 00266 return StatusCode::SUCCESS; 00267 }
|
|
Store TDS path to link a particular converter to an object on the TDS.
|
|
Store TDS path to link a particular converter to an object on the TDS.
00150 { 00151 // Purpose and Method: Save the path on the TDS, treename, pathname in the m_leaves vector, 00152 // corresponding to the DataObject that the converter handles. 00153 m_leaves.push_back(RootCnvSvc::Leaf(path, cl,treename,branchname)); 00154 }
|
|
Resolve the references of the converted object.
|
|
Resolve the references of the converted object.
00117 { 00118 // Purpose and Method: Resolve the references of the converted object. 00119 // It is expected that derived classes will override this method. 00120 MsgStream log(msgSvc(), "RootEventBaseCnv"); 00121 return StatusCode::SUCCESS; 00122 }
|
|
Resolve the references of the converted object.
|
|
Resolve the references of the converted object.
00109 { 00110 // Purpose and Method: Resolve the references of the converted object. 00111 // It is expected that derived classes will override this method. 00112 MsgStream log(msgSvc(), "RootEventBaseCnv"); 00113 return StatusCode::SUCCESS; 00114 }
|
|
|
|
00141 { 00142 if ( m_cnvSvc ) { 00143 m_cnvSvc->release(); 00144 m_cnvSvc=0; 00145 } 00146 return Converter::finalize(); 00147 }
|
|
get the object to be read
00124 { return m_objRead;}
|
|
get the object to be read
00124 { return m_objRead;}
|
|
Reimplemented in DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, TrigCnv, DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, and TrigCnv. |
|
Reimplemented in DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, TrigCnv, DigiCnv, DstCnv, EvtHeaderCnv, EvtRecCnv, HltCnv, McCnv, RecTrackCnv, and TrigCnv. 00125 { 00126 00127 StatusCode status = Converter::initialize(); 00128 00129 if ( status.isSuccess() ) { 00130 IService* isvc = 0; 00131 status = serviceLocator()->service("RootCnvSvc", isvc, false); 00132 if ( !status.isSuccess() ) status = serviceLocator()->service("EventCnvSvc", isvc, true); 00133 if ( status.isSuccess() ) { 00134 status = isvc->queryInterface(IID_IRootCnvSvc, (void**)&m_cnvSvc); 00135 } 00136 } 00137 00138 return status; 00139 }
|
|
|
|
|
|
00087 {
00088 return ROOT_StorageType;
00089 }
|
|
00087 {
00088 return ROOT_StorageType;
00089 }
|
|
Storage type and class ID.
00083 {
00084 return ROOT_StorageType;
00085 }
|
|
Storage type and class ID.
00083 {
00084 return ROOT_StorageType;
00085 }
|
|
transformation from root
Implements RootEventBaseCnv. |
|
transformation from root
Implements RootEventBaseCnv. 00045 { 00046 // creation of TDS object from root object 00047 MsgStream log(msgSvc(), "RecMdcTrackCnv"); 00048 log << MSG::DEBUG << "RecMdcTrackCnv::TObjectToDataObject" << endreq; 00049 00050 // create the TDS location for the RecMdcTrack Collection 00051 RecMdcTrackCol* recMdcTrackCol = new RecMdcTrackCol; 00052 refpObject = recMdcTrackCol; 00053 // now convert 00054 if (!m_recMdcTrackCol) return StatusCode::SUCCESS; 00055 TIter mdcTrackIter(m_recMdcTrackCol); 00056 00057 //----------get hit smartRefVector-------- 00058 //HitRefVec *hitRefVec; 00059 //int nTk=0; 00060 //while (mdcTrackIter.Next()){ nTk++; } 00061 //hitRefVec = new HitRefVec[nTk]; 00062 00063 IDataProviderSvc* dataSvc = 0; 00064 StatusCode sc = serviceLocator()->getService ("EventDataSvc", 00065 IID_IDataProviderSvc, (IInterface*&)dataSvc); 00066 if (!sc.isSuccess()) { 00067 log << MSG::FATAL << "Could not get EventDataSvc in RecMdcTrackCnv" << endreq; 00068 return( StatusCode::FAILURE); 00069 } 00070 SmartDataPtr<RecMdcHitCol> recMdcHitCol(dataSvc,"/Event/Recon/RecMdcHitCol"); 00071 if (!recMdcHitCol) { 00072 log << MSG::FATAL << "Could not find RecMdcHitCol" << endreq; 00073 return( StatusCode::FAILURE); 00074 } 00075 //RecMdcHitCol::iterator iter = recMdcHitCol->begin(); 00076 //for (;iter != recMdcHitCol->end(); iter++ ) { 00077 // int tkId = (*iter)->getTrkId(); 00078 // hitRefVec[tkId].push_back(*iter); 00079 //} 00080 00081 //--------------------------------------- 00082 mdcTrackIter = m_recMdcTrackCol; 00083 TRecMdcTrack *recMdcTrackRoot = 0; 00084 while ((recMdcTrackRoot = (TRecMdcTrack*)mdcTrackIter.Next())) { 00085 int trackId = recMdcTrackRoot->trackId(); 00086 int charge = recMdcTrackRoot->charge(); 00087 double helix[5]; 00088 for(int i=0;i<5; i++){ 00089 helix[i] = recMdcTrackRoot->helix(i); 00090 } 00091 double pxy = recMdcTrackRoot->pxy(); 00092 double px = recMdcTrackRoot->px(); 00093 double py = recMdcTrackRoot->py(); 00094 double pz = recMdcTrackRoot->pz(); 00095 double p = recMdcTrackRoot->p(); 00096 double theta = recMdcTrackRoot->theta(); 00097 double phi = recMdcTrackRoot->phi(); 00098 double x = recMdcTrackRoot->x(); 00099 double y = recMdcTrackRoot->y(); 00100 double z = recMdcTrackRoot->z(); 00101 double r = recMdcTrackRoot->r(); 00102 int stat = recMdcTrackRoot->stat(); 00103 double chi = recMdcTrackRoot->chi2(); 00104 int ndof = recMdcTrackRoot->ndof(); 00105 double err[15]; 00106 for (int i=0; i<15; i++){ 00107 err[i] = recMdcTrackRoot->err(i); 00108 } 00109 int nh = recMdcTrackRoot->nhits(); 00110 int ns = recMdcTrackRoot->nster(); 00111 // vector<int> vecHits = recMdcTrackRoot->vecHits(); 00112 double vx0 = recMdcTrackRoot->vx0(); 00113 double vy0 = recMdcTrackRoot->vy0(); 00114 double vz0 = recMdcTrackRoot->vz0(); 00115 double fiterm = recMdcTrackRoot->fiTerm(); 00116 RecMdcTrack *recMdcTrack = new RecMdcTrack(); 00117 m_common.m_rootRecMdcTrackMap[recMdcTrackRoot] = recMdcTrack; 00118 00119 recMdcTrack->setTrackId( trackId ); 00120 recMdcTrack->setCharge( charge ); 00121 recMdcTrack->setHelix( helix ); 00122 recMdcTrack->setPxy( pxy ); 00123 recMdcTrack->setPx( px ); 00124 recMdcTrack->setPy( py ); 00125 recMdcTrack->setPz( pz ); 00126 recMdcTrack->setP( p ); 00127 recMdcTrack->setTheta( theta ); 00128 recMdcTrack->setPhi( phi ); 00129 recMdcTrack->setX( x ); 00130 recMdcTrack->setY( y ); 00131 recMdcTrack->setZ( z ); 00132 recMdcTrack->setR( r ); 00133 recMdcTrack->setStat( stat ); 00134 recMdcTrack->setChi2( chi ); 00135 recMdcTrack->setNdof( ndof ); 00136 recMdcTrack->setError( err ); 00137 recMdcTrack->setNhits( nh ); 00138 recMdcTrack->setNster( ns ); 00139 // recMdcTrack->setVecHits( vecHits ); 00140 //rec 00141 recMdcTrack->setVX0( vx0 ); 00142 recMdcTrack->setVY0( vy0 ); 00143 recMdcTrack->setVZ0( vz0 ); 00144 recMdcTrack->setFiTerm( fiterm ); 00145 00146 HitRefVec theHitRefVec; 00147 00148 RecMdcHitCol::iterator iter = recMdcHitCol->begin(); 00149 for (;iter != recMdcHitCol->end(); iter++){ 00150 00151 //cout<<" (*iter)->getTrkId(): "<<(*iter)->getTrkId()<<endl; 00152 if((*iter)->getTrkId() == trackId){ 00153 SmartRef<RecMdcHit> refhit(*iter); 00154 theHitRefVec.push_back(refhit); 00155 } 00156 } 00157 00158 recMdcTrack->setVecHits(theHitRefVec); 00159 00160 int nhits = recMdcTrack->getVecHits().size(); 00161 00162 //std::cout<<" mdc hits: "<<nhits<< std::endl; 00163 00164 for(int ii=0; ii <nhits ; ii++){ 00165 00166 //cout<<"ddl: "<<(recMdcTrack->getVecHits()[ii])->getDriftDistLeft()<<endl; 00167 //cout<<"erddl: "<<(recMdcTrack->getVecHits()[ii])->getErrDriftDistLeft()<<endl; 00168 Identifier id(recMdcTrack->getVecHits()[ii]->getMdcId()); 00169 int layer = MdcID::layer(id); 00170 int wire = MdcID::wire(id); 00171 //cout<<"layer: "<<layer<<" wire: "<<wire<<endl; 00172 } 00173 recMdcTrackCol->push_back(recMdcTrack); 00174 } 00175 00176 delete m_recMdcTrackCol; 00177 m_recMdcTrackCol = 0; 00178 return StatusCode::SUCCESS; 00179 }
|
|
|
|
the CLID of the upper converter if any
|
|
each converter knows the corresponding adresses
|
|
each converter knows the corresponding adresses
|
|
the branchNr of this converter for writing
|
|
|
|
|
|
|
|
|
|
|
|
array with number of branches for reading
|
|
array with number of branches for reading
|
|
|
|
|
|
relational maps
|
|
|
|
pointer to eventdataservice
|
|
pointer to eventdataservice
|
|
|
|
|
|
|
|
|
|
the object that was read
|
|
the object that was read
|
|
root object to be read
|
|
root object to be read
|
|
root branchname (may be concatenated of severals)
|
|
pointer to the RootInterface
|
|
pointer to the RootInterface
|
|
each converter knows it's treename
|