EventTagSvc Class Reference

#include <EventTagSvc.h>

Inheritance diagram for EventTagSvc:

List of all members.

Public Member Functions
	EventTagSvc (const std::string &name, ISvcLocator *svc)
	EventTagSvc (const std::string &name, ISvcLocator *svc)
StatusCode	finalize ()
StatusCode	finalize ()
unsigned long long int	getChainCode (Event::McParticle *part)
unsigned long long int	getChainCode (Event::McParticle *part)
unsigned int	getCharmDecayType (Event::McParticle *part)
unsigned int	getCharmDecayType (Event::McParticle *part)
int	getDecayCode (Event::McParticle *part)
int	getDecayCode (Event::McParticle *part)
StatusCode	initialize ()
StatusCode	initialize ()
virtual long int	name2pdg (std::string pname)=0
long int	name2pdg (string pname)
long int	name2pdg (string pname)
string	pdg2name (long int pdg)
string	pdg2name (long int pdg)
virtual StatusCode	queryInterface (const InterfaceID &, void **ppvInterface)
virtual StatusCode	queryInterface (const InterfaceID &, void **ppvInterface)
Static Public Member Functions
const InterfaceID &	interfaceID ()
const InterfaceID &	interfaceID ()
Private Member Functions
int	buildAntiPartTabs ()
int	buildAntiPartTabs ()
int	readDecayTabFile (string fname)
int	readDecayTabFile (string fname)
int	readPdtFile (string fname)
int	readPdtFile (string fname)
Private Attributes
std::set< long int >	chainTrigParticles
std::set< long int >	chainTrigParticles
bool	m_buildAntiTabs
std::map< int, decayModeTab * >	m_decayTabList
std::map< int, decayModeTab * >	m_decayTabList
string	m_decayTabsFile
bool	m_digK0
bool	m_ignoreFSR
bool	m_ignorePhotons
string	m_pdtFile
string	m_userDecayTabsFile
std::map< string, long int >	name2pdgTab
std::map< string, long int >	name2pdgTab
std::map< long int, string >	pdg2nameTab
std::map< long int, string >	pdg2nameTab
std::vector< string >	userChainTrig
std::vector< string >	userChainTrig

---

Constructor & Destructor Documentation

EventTagSvc::EventTagSvc (  const std::string &  name, ISvcLocator *  svc )

: Service(name,svc) { declareProperty("pdgFile", m_pdtFile="pdt.table"); declareProperty("decayCodesFile", m_decayTabsFile="decay.codes"); // declareProperty("userDecayCodesFile", m_userDecayTabsFile="userdecay.codes"); declareProperty("userDecayCodesFile", m_userDecayTabsFile=""); declareProperty("buildAntiPartTable", m_buildAntiTabs=true); declareProperty("ignorePhotons", m_ignorePhotons=true); declareProperty("ignoreFSR", m_ignoreFSR=true); declareProperty("chainTrigParticles", userChainTrig); declareProperty("digK0", m_digK0=true); }

EventTagSvc::EventTagSvc (  const std::string &  name, ISvcLocator *  svc )

---

Member Function Documentation

int EventTagSvc::buildAntiPartTabs (   )  [private]

int EventTagSvc::buildAntiPartTabs (   )  [private]

{ MsgStream log(msgSvc(), name()); decayModeTab* dTab; decayModeTab* adTab; keyVector key,akey; int mother,amother,code; for(std::map<int,decayModeTab*>::iterator itl=m_decayTabList.begin(); itl!=m_decayTabList.end();itl++){ mother=itl->first; amother=-mother; if(m_decayTabList.find(amother)!=m_decayTabList.end())continue ;// aparticle tab already exist if(pdg2name(mother).size()==0){ log<<MSG::ERROR<<" buildAntiPartTabs :: Unknown pdg code for mother particle "<<endreq; return 1; } if(pdg2name(amother).size()==0)continue; //Aparticle==particle?? log<<MSG::INFO<<" buildAntiPartTabs :: create new table for particle " << amother<<"("<<pdg2name(amother)<<")"<<endreq; dTab=itl->second; adTab=new decayModeTab; m_decayTabList[amother]=adTab; key.clear(); for(decayModeTab::iterator it=dTab->begin();it!=dTab->end();it++){ key=it->first; code=it->second; for(keyVector::iterator it=key.begin();it!=key.end();it++){ if(pdg2name(*it).size()==0){ log<<MSG::ERROR<<" buildAntiPartTabs :: Unknown pdg code" << "for particle in keyVector"<<endreq; return 1; } if(pdg2name(-(*it)).size()==0)continue ; // true neutral ?? (*it) *= -1; } sort(key.begin(),key.end()); (*adTab)[key]=code; } } return 0; }

StatusCode EventTagSvc::finalize (   )

StatusCode EventTagSvc::finalize (   )

{ // Part 1: Get the messaging service, print where you are MsgStream log(msgSvc(), name()); log << MSG::INFO << "EventTagSvc finalize()" << endreq; return StatusCode::SUCCESS; }

unsigned long long int EventTagSvc::getChainCode (  Event::McParticle *  part  )  [virtual]

Implements IEventTagSvc.

unsigned long long int EventTagSvc::getChainCode (  Event::McParticle *  part  )  [virtual]

Implements IEventTagSvc. { unsigned long long int code=0; int shift=0; vector<Event::McParticle*> chainVect; while(shift<57){ code |= (getDecayCode(part)&0xFF)<<shift; SmartRefVector<Event::McParticle> dref=part->daughterList(); chainVect.clear(); for(SmartRefVector<Event::McParticle>::iterator it=dref.begin();it!=dref.end();it++){ int abspdg=abs((*it)->particleProperty()); if(((((abspdg%1000)/10)<45) && (((abspdg%1000)/10)>40)) // select particles with *4?? pdg code ||(chainTrigParticles.find(abspdg)!=chainTrigParticles.end())) { chainVect.push_back(*it); } } if(chainVect.size()==0 ||chainVect.size()>2)break; if(chainVect.size()==1){part=chainVect[0];} else{ sort(chainVect.begin(),chainVect.end(),mcPartPrtLess); code |= ((getDecayCode(chainVect[0])&0xFF)<<(shift+8)); code |= ((getDecayCode(chainVect[1])&0xFF)<<(shift+16)); break; } shift+=8; } return code; }

unsigned int EventTagSvc::getCharmDecayType (  Event::McParticle *  part  )  [virtual]

Implements IEventTagSvc.

unsigned int EventTagSvc::getCharmDecayType (  Event::McParticle *  part  )  [virtual]

Implements IEventTagSvc. { unsigned int a=0; // int type=0; if(! part) return 0xF; SmartRefVector<Event::McParticle> dref=part->daughterList(); if(dref.size()==0) return 0; int nhad=0; long int pdg1=0; bool dd=false; for(SmartRefVector<Event::McParticle>::iterator it=dref.begin(); it!=dref.end();it++){ int pdg=abs((*it)->particleProperty()); if((pdg==22)&&m_ignorePhotons)continue; if(pdg==22){a|=1;} else if( (pdg>10)&&(pdg<19) ){a |=(1<<1);} else if( pdg>110) { if( ((pdg%1000)/100==4)&&((pdg%100)/10==4)){ a |= 1<<3; } else{ a |= 1<<2; nhad+=1; if(nhad==1)pdg1=pdg; dd=(pdg1==pdg); } } else a |=1<<4; } if(a>0xF)return 0xf; //unknown particle if(a==4){ // hadronic if((nhad!=2)||(!dd))return 6; // long int pdg1=dref[0]->particleProperty(); // long int pdg2=dref[1]->particleProperty(); // if(abs(pdg1)!=abs(pdg2)) return 6; switch (abs(pdg1)){ case 411: return 7; break; case 421: return 8; break; case 413: return 9; break; case 423: return 10; break; case 431: return 11; break; case 433: return 12; break; default: return 6; break; } } if((a&0xE)==2)return 1; // leptonic if((a&0x6)==0x6)return 2; // semileptonic if((a&0xE)==0xC)return 4; // hadronic charmonium transition if((a&0xE)==0x8)return 3; // radiative charmonium transition !!use 0XE mask for correct work with "ignore photons" options; if((a&1)==1)return 5; // radiative return 0xf; //unknown }

int EventTagSvc::getDecayCode (  Event::McParticle *  part  )  [virtual]

Implements IEventTagSvc.

int EventTagSvc::getDecayCode (  Event::McParticle *  part  )  [virtual]

Implements IEventTagSvc. { MsgStream log(msgSvc(), name()); keyVector key; long int mother_pdg=part->particleProperty(); map<int,decayModeTab*>::iterator itlist=m_decayTabList.find(mother_pdg); if(itlist==m_decayTabList.end())return 0; decayModeTab* dTab=itlist->second; SmartRefVector<Event::McParticle> dref=part->daughterList(); if(dref.size()==0)return 0; key.clear(); for(SmartRefVector<Event::McParticle>::iterator it=dref.begin();it!=dref.end();it++){ int pdg=(*it)->particleProperty(); if(m_digK0&&(abs(pdg)==311)){ SmartRefVector<Event::McParticle> dref1=(*it)->daughterList(); if(dref1.size()!=1){ log<< MSG::WARNING<<"unknown decay of K0, aborting"<<endreq; return 0; } pdg=dref1[0]->particleProperty(); } if((abs(pdg)==22)&&m_ignorePhotons)continue; if((pdg==-22)&&m_ignoreFSR)continue; key.push_back(pdg); } sort(key.begin(),key.end()); if(m_ignorePhotons&&(key.size()==0))key.push_back(22); decayModeTab::iterator ittab=dTab->find(key); if(ittab==dTab->end())return 0; return ittab->second; }

StatusCode EventTagSvc::initialize (   )

StatusCode EventTagSvc::initialize (   )

{ Service::initialize(); MsgStream log(msgSvc(), name()); setProperties(); log << MSG::INFO << "EventTagSvc initialize()" << endreq; log << MSG::INFO << " pdgFile="<<m_pdtFile<<endreq; log << MSG::INFO << " decayCodesFile="<<m_decayTabsFile<<endreq; log << MSG::INFO << " userDecayTabsFile="<<m_userDecayTabsFile<<endreq; log << MSG::INFO << " buildAntiPartTable="<<m_buildAntiTabs<<endreq; log << MSG::INFO << " Ignore Photons="<<m_ignorePhotons<<endreq; log << MSG::INFO << " Chain trigger particles:"<<endreq; for(std::vector<string>::iterator it=userChainTrig.begin();it!=userChainTrig.end();it++){ log << MSG::INFO<<" " <<(*it)<<endreq; } // static const bool CREATEIFNOTTHERE(true); if(readPdtFile(m_pdtFile.c_str())){ log << MSG::ERROR<<"Can not read pdgFile "<< m_pdtFile<<endreq; return StatusCode::StatusCode::FAILURE; } if(readDecayTabFile(m_decayTabsFile)){ log << MSG::ERROR<<"Can not read decayTabsFile "<< m_decayTabsFile<<endreq; return StatusCode::StatusCode::FAILURE; } if(m_userDecayTabsFile.size()){ if(readDecayTabFile(m_userDecayTabsFile)){ log << MSG::ERROR<<"Can not read userDecayTabsFile "<<m_userDecayTabsFile<<endreq; return StatusCode::StatusCode::FAILURE; } } if(m_buildAntiTabs){ if(buildAntiPartTabs()){ log << MSG::ERROR<<"fail to build decay code table for anti-particles"<<endreq; return StatusCode::StatusCode::FAILURE; } } // fill the set with trigger chain particles for(std::vector<string>::iterator it=userChainTrig.begin();it!=userChainTrig.end();it++){ long int pdg=abs(name2pdg(*it)); if(pdg==0){ log << MSG::ERROR<<"Unknown ChainTrigger particle "<< (*it)<<" ..ignoring"<<endreq; continue; } chainTrigParticles.insert(pdg); } log << MSG::DEBUG << "sorted chain trigger particles:"<<endreq; for(set<long int>::iterator it=chainTrigParticles.begin();it!=chainTrigParticles.end();it++){ log << MSG::DEBUG<<" "<<(*it)<<" " <<pdg2name(*it)<<endreq; } return StatusCode::SUCCESS; }

const InterfaceID& IEventTagSvc::interfaceID (   )  [static, inherited]

const InterfaceID & IEventTagSvc::interfaceID (   )  [inline, static, inherited]

{ static const InterfaceID _IID("IEventtagsvc", 1, 0); // static const InterfaceID _IID(7091, 1, 0); return _IID; }

virtual long int IEventTagSvc::name2pdg (  std::string  pname  )  [pure virtual, inherited]

long int EventTagSvc::name2pdg (  string  pname  )

long int EventTagSvc::name2pdg (  string  pname  )

{ map<string,long int>::iterator it=name2pdgTab.find(pname); if(it!=name2pdgTab.end())return it->second; return 0; }

string EventTagSvc::pdg2name (  long int  pdg  )  [virtual]

Implements IEventTagSvc.

string EventTagSvc::pdg2name (  long int  pdg  )  [virtual]

Implements IEventTagSvc. { map<long int,string>::iterator it=pdg2nameTab.find(pdg); if(it!=pdg2nameTab.end())return it->second; return string(""); }

virtual StatusCode EventTagSvc::queryInterface (  const InterfaceID &  , void **  ppvInterface )  [virtual]

StatusCode EventTagSvc::queryInterface (  const InterfaceID &  , void **  ppvInterface )  [virtual]

{ if ( IEventTagSvc::interfaceID().versionMatch(riid) ) { *ppvInterface = (IEventTagSvc*)this; } else { // Interface is not directly available: try out a base class return Service::queryInterface(riid, ppvInterface); } addRef(); return StatusCode::SUCCESS; }

int EventTagSvc::readDecayTabFile (  string  fname  )  [private]

int EventTagSvc::readDecayTabFile (  string  fname  )  [private]

{ MsgStream log(msgSvc(), name()); ifstream infile; infile.open(fname.c_str()); if(!infile.is_open()){ log<<MSG::ERROR<< "Can not open decay table file :" <<fname<<endreq; infile.close(); return 1; } log<<MSG::INFO<<"Reading decay codes from file " << fname<<endreq; string ln; long int pdg,mo_pdg; int dcode; int nstr=0; char* pch; char str[256],mother[256]; char whitespaces[]=" \t\f\v\n\r"; while (! infile.eof() ){ getline (infile,ln); nstr++; if(ln.empty())continue; if(ln[0]=='#'||ln[0]=='*')continue; if(ln.size()>255){ log<<MSG::ERROR<<"line "<< nstr <<"in " <<fname<<" is too big" <<endreq; infile.close(); return 1; } strcpy(str,ln.c_str()); pch = strtok (str,whitespaces); if(pch==NULL)continue; if(strcmp(pch,"Decay")!=0){ log<<MSG::ERROR<<"No open decay, unrecognized token in line " <<nstr <<" of " <<fname<<endreq; infile.close(); return 1; } pch = strtok (NULL, whitespaces); if(pch==NULL) { log<<MSG::ERROR<< "Can not find particles in line " << nstr <<" of " <<fname<<endreq; infile.close(); return 1; } mo_pdg=name2pdg(string(pch)); strcpy(mother,pch); if(mo_pdg==0){ log<<MSG::ERROR<< "Unrecognized particle "<<pch<<" in line " << nstr <<" of " <<fname<<endreq; infile.close(); return 1; } log<<MSG::INFO<<"START DECAY TABLE FOR " << mother<<", pdg_code:"<<mo_pdg<<endl; decayModeTab* dTab; map<int,decayModeTab*>::iterator itlist=m_decayTabList.find(mo_pdg); if(itlist!=m_decayTabList.end()){ log<<MSG::WARNING<< "add new info into existing decay tab for "<<mother<<endreq; dTab=itlist->second; } else{ log<<MSG::INFO<< "Create new decay code table for "<<mother<<endreq; dTab=new decayModeTab; m_decayTabList[mo_pdg]=dTab; } while(true){ getline (infile,ln); if( infile.eof() ){ log<<MSG::ERROR << "Enexpected EOF for dacay code file " <<fname <<endreq; infile.close(); return 1; } nstr++; if(ln.empty())continue; if(ln[0]=='#')continue; if(ln.size()>255){ log<<MSG::ERROR<<"line "<< nstr <<"in decay table file " <<fname<<" is too big" <<endreq; infile.close(); return 1; } strcpy(str,ln.c_str()); pch = strtok (str,whitespaces); if(pch==NULL)continue; if(strcmp(pch,"Enddecay")==0){ log<<MSG::INFO <<" End of "<< mother <<" decay table" <<endreq; break; } log<<MSG::DEBUG<<"parsing line "<<ln<<endreq; dcode=atoi(pch); if(dcode<=0||dcode>255){ log<<MSG::ERROR<< " wrong decay code in line "<<nstr <<", file "<<fname<< endreq; infile.close(); return 1; } // decayModeTab::iterator itTab=dTab->find(code); //if( itTab != dTab->end()){ log<<MSG::DEBUG<<" dcode:"<<dcode<<", particles: "; pch = strtok (NULL, whitespaces); if(pch==NULL){ // if( log<<MSG::ERROR<< " no decay found in line "<<nstr <<", file "<<fname<< endreq; infile.close(); return 1; } // keyVector* pKey=new keyVector; keyVector key; key.clear(); // keyVector* pKey=&key; while (pch != NULL){ log<<MSG::DEBUG <<pch<<" "; pdg=name2pdg(string(pch)); if(pdg==0){ log<<MSG::ERROR<<endl<< "Unrecognized particle " <<pch <<" in line " << nstr <<" of " <<fname<<endreq; infile.close(); return 1; } // pKey->push_back(pdg); key.push_back(pdg); pch = strtok (NULL, whitespaces); } // sort(pKey->begin(),pKey->end()); sort(key.begin(),key.end()); log<<MSG::DEBUG<<" sorted pdg" ; // for(keyVector::iterator it=pKey->begin(); it!=pKey->end(); it++) for(keyVector::iterator it=key.begin(); it!=key.end(); it++) log<<MSG::DEBUG <<" "<<(*it); log<<MSG::DEBUG<<endreq; // decayModeTab::iterator itdec=dTab->find(*pKey); decayModeTab::iterator itdec=dTab->find(key); if(itdec!=dTab->end()){ if(itdec->second !=dcode){ log<<MSG::ERROR << "line "<<nstr<<" of file " <<fname<<"has an error : "<<endreq ; log<<MSG::ERROR<<"decay of "<<mother<< " to"; // for(keyVector::iterator it=pKey->begin(); it!=pKey->end(); it++) for(keyVector::iterator it=key.begin(); it!=key.end(); it++) log<<MSG::ERROR<<" "<<pdg2name(*it); log<<MSG::ERROR<<" is already in table. Aborting."<<endreq; log<<MSG::ERROR<<"check line "<<nstr<<" of " <<fname<<endreq; infile.close(); return 1; } } // (*dTab)[*pKey]=dcode; (*dTab)[key]=dcode; } } // end if file read loop infile.close(); return 0; }

int EventTagSvc::readPdtFile (  string  fname  )  [private]

int EventTagSvc::readPdtFile (  string  fname  )  [private]

{ MsgStream log(msgSvc(), name()); ifstream infile; infile.open(fname.c_str()); if(!infile.is_open()){ log<<MSG::ERROR<< "Can not open pdg-code file :" <<fname<<endreq; return 1; } log<<MSG::INFO<<"reading pdg codes table from "<<fname<<endreq; string ln; int nstr=0; long int pcode; string cmd,dummy,pname,pcode_str; stringstream* lns; while (! infile.eof() ){ getline (infile,ln); nstr++; log<<MSG::DEBUG<<"parsing line:"<<ln<<endreq; if(ln.empty())continue; if(ln[0]=='*')continue; lns=new stringstream; (*lns)<<ln; cmd=""; (*lns)>>cmd; if(!strcmp(cmd.c_str(),"add")){ pcode=0; pname=""; pcode_str=""; (*lns)>>dummy>>dummy>>pname>>pcode_str; pcode=atol(pcode_str.c_str()); if(pcode==0||pcode==LONG_MAX||pcode==LONG_MIN){ log<<MSG::ERROR<<"Unrecognized particle code in line " <<nstr<<" of file" <<fname<<endreq; infile.close(); return 1; } log<<MSG::DEBUG<<" find particle: " <<pname<<" pgd:"<<pcode <<endreq; if(pdg2nameTab.find(pcode)!=pdg2nameTab.end()){ log<<MSG::ERROR<< "pdg code "<<pcode <<" is repeated in line "<<nstr<<". Aborting"<<endreq; infile.close(); return 1; } if(name2pdgTab.find(pname)!=name2pdgTab.end()){ log<<MSG::ERROR<< "particle name "<<pcode <<" is repeated in line "<<nstr<<". Aborting"<<endreq; infile.close(); return 1; } pdg2nameTab[pcode]=pname; name2pdgTab[pname]=pcode; } delete lns; } infile.close(); return 0; }

---

Member Data Documentation

std::set<long int> EventTagSvc::chainTrigParticles [private]

std::set<long int> EventTagSvc::chainTrigParticles [private]

bool EventTagSvc::m_buildAntiTabs [private]

std::map<int,decayModeTab*> EventTagSvc::m_decayTabList [private]

std::map<int,decayModeTab*> EventTagSvc::m_decayTabList [private]

string EventTagSvc::m_decayTabsFile [private]

bool EventTagSvc::m_digK0 [private]

bool EventTagSvc::m_ignoreFSR [private]

bool EventTagSvc::m_ignorePhotons [private]

string EventTagSvc::m_pdtFile [private]

string EventTagSvc::m_userDecayTabsFile [private]

std::map<string,long int> EventTagSvc::name2pdgTab [private]

std::map<string,long int> EventTagSvc::name2pdgTab [private]

std::map<long int,string> EventTagSvc::pdg2nameTab [private]

std::map<long int,string> EventTagSvc::pdg2nameTab [private]

std::vector<string> EventTagSvc::userChainTrig [private]

std::vector<string> EventTagSvc::userChainTrig [private]

---

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

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Event/EventTag/EventTag-00-01-05/EventTag/EventTagSvc.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/InstallArea/include/EventTag/EventTag/EventTagSvc.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Event/EventTag/EventTag-00-01-05/src/EventTagSvc.cxx

---