EventTagSvc Class Reference

#include <EventTagSvc.h>

Inheritance diagram for EventTagSvc:

Public Member Functions
	EventTagSvc (const std::string &name, ISvcLocator *svc)
StatusCode	initialize ()
std::string	pdg2name (long int pdg)
long int	name2pdg (std::string pname)
StatusCode	finalize ()
int	getDecayCode (Event::McParticle *part)
unsigned long long int	getChainCode (Event::McParticle *part)
unsigned int	getCharmDecayType (Event::McParticle *part)

Detailed Description

Definition at line 26 of file EventTagSvc.h.

Constructor & Destructor Documentation

EventTagSvc()

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

Definition at line 24 of file EventTagSvc.cxx.

                                                                 :
  base_class(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);
}

Member Function Documentation

finalize()

StatusCode EventTagSvc::finalize

(

)

Definition at line 186 of file EventTagSvc.cxx.

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

getChainCode()

unsigned long long int EventTagSvc::getChainCode

(

Event::McParticle *

part

)

Definition at line 809 of file EventTagSvc.cxx.

                                                                     {
#else
unsigned long long int EventTagSvc::getChainCode (TMcParticle * part){
   if (!m_TMcEvent) {
      cerr << "(ERROR) EventTag: you must set TMcEvent first!" << endl;
      return -1;
   }
#endif
  unsigned long long int code=0;
  int shift=0;
#ifndef BEAN
  vector<Event::McParticle*> chainVect;
#else
  vector < TMcParticle * >chainVect;
#endif
  while(shift<57){
    code |=  (getDecayCode(part)&0xFF)<<shift;
#ifndef BEAN
    SmartRefVector<Event::McParticle> dref=part->daughterList();
#else
    vector < int >dref = part->getDaughters ();
#endif
    chainVect.clear();
#ifndef BEAN
    for(SmartRefVector<Event::McParticle>::iterator it=dref.begin();
     it!=dref.end();it++){
      int abspdg=abs((*it)->particleProperty());
#else
    for (unsigned int i = 0; i < dref.size (); i++) {
      const TMcParticle *daughter = m_TMcEvent->getMcParticle (dref[i]);
 
      int abspdg = abs (daughter->getParticleID ());
#endif
      if(((((abspdg%1000)/10)<45) && (((abspdg%1000)/10)>40)) // select particles with *4?? pdg code
         ||(chainTrigParticles.find(abspdg)!=chainTrigParticles.end())) {
#ifndef BEAN
        chainVect.push_back(*it);
#else
        chainVect.push_back (const_cast < TMcParticle * >(daughter));
#endif
      }
    }
    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;
}

Referenced by getChainCode().

getCharmDecayType()

unsigned int EventTagSvc::getCharmDecayType

(

Event::McParticle *

part

)

Definition at line 723 of file EventTagSvc.cxx.

                                                                 {
#else
unsigned int EventTagSvc::getCharmDecayType (TMcParticle * part){
   if (!m_TMcEvent) {
      cerr << "(ERROR) EventTag: you must set TMcEvent first!" << endl;
      return -1;
   }
#endif
  unsigned int a=0;
  //  int type=0;
  if(! part) return 0xF;
#ifndef BEAN
  SmartRefVector<Event::McParticle> dref=part->daughterList();
#else
  vector < int >dref = part->getDaughters ();
#endif
  if(dref.size()==0) return 0;
  int nhad=0;
  long int pdg1=0;
  bool dd=false;
#ifndef BEAN
  for(SmartRefVector<Event::McParticle>::iterator it=dref.begin();
   it!=dref.end();it++){
    int pdg=abs((*it)->particleProperty());
#else
  for (unsigned int i = 0; i < dref.size (); i++) {
    const TMcParticle *daughter = m_TMcEvent->getMcParticle (dref[i]);
 
    int pdg = abs (daughter->getParticleID ());
#endif
    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
}

Referenced by getChainCode(), and getCharmDecayType().

getDecayCode()

int EventTagSvc::getDecayCode

(

Event::McParticle *

part

)

Definition at line 590 of file EventTagSvc.cxx.

                                                  {
  MsgStream log(msgSvc(), name());
#else
int EventTagSvc::getDecayCode (TMcParticle * part){
   if (!m_TMcEvent) {
      cerr << "(ERROR) EventTag: you must set TMcEvent first!" << endl;
      return -1;
   }
#endif
  keyVector key;
 
#ifndef BEAN
  long int mother_pdg=part->particleProperty();
#else
  long int mother_pdg = part->getParticleID ();
#endif
 
  map<int,decayModeTab*>::iterator itlist=m_decayTabList.find(mother_pdg);
  if(itlist==m_decayTabList.end())return 0;
  decayModeTab* dTab=itlist->second;
 
#ifndef BEAN
  SmartRefVector<Event::McParticle> dref=part->daughterList();
#else
  vector < int >dref = part->getDaughters ();
#endif
  if(dref.size()==0)return 0;
  key.clear();
 
#ifndef BEAN
  for(SmartRefVector<Event::McParticle>::iterator it=dref.begin();
   it!=dref.end();it++){
    int pdg=(*it)->particleProperty();
#else
   for (unsigned int i = 0; i < dref.size (); i++) {
      const TMcParticle *daughter = m_TMcEvent->getMcParticle (dref[i]);
 
      int pdg = daughter->getParticleID ();
#endif
    if(m_digK0&&(abs(pdg)==311)){
#ifndef BEAN
      SmartRefVector<Event::McParticle> dref1=(*it)->daughterList();
#else
      vector < int >dref1 = daughter->getDaughters ();
#endif
      if(dref1.size()!=1){
#ifndef BEAN
        log<< MSG::WARNING<<"unknown decay of K0, aborting"<<endreq;
#else
        cout << "unknown decay of K0, aborting" << endl;
#endif
        return 0;
      }
#ifndef BEAN
      pdg=dref1[0]->particleProperty();
#else
      pdg = m_TMcEvent->getMcParticle (dref1[0])->getParticleID ();
#endif
    }
    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;
}

Referenced by getChainCode(), and getDecayCode().

initialize()

StatusCode EventTagSvc::initialize

(

)

Definition at line 60 of file EventTagSvc.cxx.

                                  {
#ifndef BEAN
  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;
  }
#else
   if( verbose ) {
      cout << "EventTagSvc initialize()" << endl;
      cout << "         pdgFile=" << m_pdtFile << endl;
      cout << "         decayCodesFile=" << m_decayTabsFile << endl;
      cout << "         userDecayTabsFile=" << m_userDecayTabsFile << endl;
      cout << "         buildAntiPartTable=" << m_buildAntiTabs << endl;
      cout << "         Ignore Photons=" << m_ignorePhotons << endl;
      cout << " Chain trigger particles:" << endl;
 
      for (std::vector < string >::iterator it = userChainTrig.begin ();
        it != userChainTrig.end (); it++) {
         cout << "                         " << (*it) << endl;
      }
   }
#endif
 
  //  static const bool CREATEIFNOTTHERE(true);
  if(readPdtFile(m_pdtFile.c_str())){
#ifndef BEAN
    log << MSG::ERROR<<"Can not read pdgFile "<< m_pdtFile<<endreq;
    return StatusCode::StatusCode::FAILURE;
#else
    cerr << "Can not read pdgFile " << m_pdtFile << endl;
    return false;
#endif
  }
 
  if(readDecayTabFile(m_decayTabsFile)){
#ifndef BEAN
    log << MSG::ERROR<<"Can not read decayTabsFile "<< m_decayTabsFile<<endreq;
    return StatusCode::StatusCode::FAILURE;
#else
    cerr << "Can not read decayTabsFile " << m_decayTabsFile << endl;
    return false;
#endif
  }
 
  if(m_userDecayTabsFile.size()){
    if(readDecayTabFile(m_userDecayTabsFile)){
#ifndef BEAN
      log << MSG::ERROR<<"Can not read userDecayTabsFile "<<m_userDecayTabsFile
          <<endreq;
      return StatusCode::StatusCode::FAILURE;
#else
      cerr << "Can not read userDecayTabsFile " << m_userDecayTabsFile << endl;
      return false;
#endif
    }
  }
 
  if(m_buildAntiTabs){
    if(buildAntiPartTabs()){
#ifndef BEAN
      log << MSG::ERROR<<"fail to build decay code table for anti-particles"
          <<endreq;
      return StatusCode::StatusCode::FAILURE;
#else
      cerr << "fail to build decay code table for anti-particles" << endl;
      return false;
#endif
    }
  }
 
  // fill the set with trigger chain particles
  for(std::vector<string>::iterator it=userChainTrig.begin();
   it!=userChainTrig.end();it++){
    long int pdg=labs(name2pdg(*it));
    if(pdg==0){
#ifndef BEAN
      log << MSG::ERROR<<"Unknown ChainTrigger particle "<< (*it)
          <<" ..ignoring"<<endreq;
#else
      cerr << "Unknown ChainTrigger particle " << (*it)
           <<  " ..ignoring" << endl;
#endif
      continue;
    }
    chainTrigParticles.insert(pdg);
  }
 
#ifndef BEAN
  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;
  }
#else
   if (verbose) {
      cout << "sorted chain trigger particles:" << endl;
      for (set < long int >::iterator it = chainTrigParticles.begin ();
        it != chainTrigParticles.end (); it++) {
         cout << "                         " << (*it) << "  "
              << pdg2name(*it) << endl;
      }
   }
#endif
 
 
#ifndef BEAN
  return StatusCode::SUCCESS;
#else
  return true;
#endif
}

name2pdg()

long int EventTagSvc::name2pdg

(

std::string

pname

)

Definition at line 214 of file EventTagSvc.cxx.

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

Referenced by initialize().

pdg2name()

string EventTagSvc::pdg2name

(

long int

pdg

)

Definition at line 208 of file EventTagSvc.cxx.

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

Referenced by getDecayCode(), and initialize().

---

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

• EventTagSvc.h
• EventTagSvc.cxx