EvtDecayTag Class Reference

#include <EvtDecayTag.hh>

Public Member Functions
	EvtDecayTag (EvtParticle *par)
virtual	~EvtDecayTag ()
void	makeTag (EvtParticle *par)
int	TagIndex (EvtId pid)
int	getModeTag ()
int	getMultTag ()
int	getHdrdBase (int index)

Detailed Description

Definition at line 38 of file EvtDecayTag.hh.

Constructor & Destructor Documentation

EvtDecayTag()

EvtDecayTag::EvtDecayTag ( EvtParticle * par )

inline

Definition at line 42 of file EvtDecayTag.hh.

                                  {
    root_par = par;
    _nmode.clear();
    _multiplicity.clear();
    for ( int i = 0; i < 10; i++ ) { _multiplicity.push_back( 0 ); }
    makeTag( root_par );
  }

~EvtDecayTag()

virtual EvtDecayTag::~EvtDecayTag ( )

inlinevirtual

Definition at line 51 of file EvtDecayTag.hh.

{}

Member Function Documentation

getHdrdBase()

int EvtDecayTag::getHdrdBase

(

int

index

)

Definition at line 172 of file EvtDecayTag.cc.

                                        {
  int ten = 10;
  int base;
  if ( index == 0 || index >= 100 && index < 1000 ) { return index; }
  else
  {
    base = index * ten;
    return getHdrdBase( base );
  }
}

Referenced by getHdrdBase().

getModeTag()

int EvtDecayTag::getModeTag

(

)

coding root_parent_mode * 10^9 + daug_0*10^6 + daug_1*10^3

Definition at line 110 of file EvtDecayTag.cc.

                            {
  // get ConExc 74110 mode selected index
  std::vector<int> intFlag = root_par->getIntFlag();
  if ( intFlag.size() > 0 )
  {
    if ( intFlag[0] == 74110 || intFlag[0] == -100 ) return intFlag[1];
  }
 
  /**
   *  coding root_parent_mode * 10^9 + daug_0*10^6 + daug_1*10^3
   **/
  double seg0, seg1, seg2;
 
  int three   = 1000;
  int six     = 1000000;
  int themode = 1000000000;
  //---debugging
  // for(int im=0;im<_nmode.size();im++){std::cout<<"_nmode["<<im<<"]"<<_nmode[im]<<endl;;}
 
  if ( _nmode.size() == 1 )
  {
    themode += _nmode[0] * six;
    return themode;
  }
  else if ( _nmode.size() == 2 )
  {
    seg0 = _nmode[0] * six;
    seg1 = _nmode[1] * three;
    themode += seg0 + seg1;
    return themode;
  }
  else if ( _nmode.size() >= 3 )
  {
    seg0 = _nmode[0] * six;
    seg1 = _nmode[1] * three;
    seg2 = _nmode[2];
    themode += seg0 + seg1 + seg2;
    return themode;
  }
 
  return 0; // pingrg
}

getMultTag()

int EvtDecayTag::getMultTag

(

)

10 digit in the order: N_gamma, N_gamma, N_e, N_mu, N_pi, N_Kaon, N_p, N_Ks, N_Lambda,DecayType

Definition at line 153 of file EvtDecayTag.cc.

                            {
  /**
   * 10 digit in the order: N_gamma, N_gamma, N_e,  N_mu,    N_pi,
   *                        N_Kaon,  N_p,     N_Ks, N_Lambda,DecayType
   **/
  int thetag = 0;
  for ( int i = 1; i < 9; i++ )
  {
    if ( _multiplicity[i] > 9 ) return 0;
  }
  for ( int i = 0; i < 9; i++ )
  {
    int ndx = 8 - i;
    int dig = pow( 10, ndx );
    thetag += dig * _multiplicity[i];
  }
  return thetag;
}

makeTag()

void EvtDecayTag::makeTag

(

EvtParticle *

par

)

Definition at line 28 of file EvtDecayTag.cc.

                                            {
 
  int ndaug = par->getNDaug();
  if ( ndaug < 2 )
  {
    EvtId id     = par->getId();
    int   theTag = TagIndex( id );
    _multiplicity[theTag]++;
  }
  else
  {
    //--
    if ( par == root_par )
    {
      int theMode = root_par->getChannel();
      _nmode.push_back( theMode );
    }
    //--
    int          ND = 0;
    EvtParticle* theD;
    for ( int i = 0; i < ndaug; i++ )
    {
      EvtParticle* theDaug = par->getDaug( i );
      if ( par == root_par && theDaug->getNDaug() != 0 )
      {
        int theMode = theDaug->getChannel();
        _nmode.push_back( theMode );
        ND++;
        theD = theDaug;
      }
 
      int id = EvtPDL::getStdHep( theDaug->getId() );
      if ( id == 310 ) _multiplicity[6]++;  // Ks
      if ( id == 3122 ) _multiplicity[7]++; // Lambda or Lambdabar
      makeTag( theDaug );
    }
    //-- for root particle decays into only one unstable particle, the flag1 further tag the
    // third decay
    if ( ND == 1 )
    {
      int NDD = theD->getNDaug();
      for ( int i = 0; i < NDD; i++ )
      {
        EvtParticle* theDD = theD->getDaug( i );
        if ( theDD->getNDaug() != 0 )
        {
          int ndd = theDD->getChannel();
          _nmode.push_back( ndd );
          break;
        }
      }
    }
    //-------------------
  }
}

Referenced by EvtDecayTag(), and makeTag().

TagIndex()

int EvtDecayTag::TagIndex

(

EvtId

pid

)

gamma : 0 e+/e- : 1 mu+/- : 2 pi+/- : 3 K+/- : 4 p+/- : 5 Ks : 6 Lambda: 7 others: 8

Definition at line 84 of file EvtDecayTag.cc.

                                     {
  /**
   * gamma : 0
   * e+/e- : 1
   * mu+/- : 2
   * pi+/- : 3
   * K+/-  : 4
   * p+/-  : 5
   * Ks    : 6
   * Lambda: 7
   * others: 8
   **/
 
  int id    = EvtPDL::getStdHep( pid );
  int absid = fabs( id );
  if ( absid == 22 ) { return 0; }        // photon
  else if ( absid == 11 ) { return 1; }   // electron
  else if ( absid == 13 ) { return 2; }   // muon
  else if ( absid == 211 ) { return 3; }  // pion
  else if ( absid == 321 ) { return 4; }  // Kaon
  else if ( absid == 2212 ) { return 5; } // pronton / anti-proton
  else if ( absid == 310 ) { return 6; }  // Ks
  else if ( absid == 3122 ) { return 7; } // Lambda / Lambdabar
  else { return 8; }
}

Referenced by makeTag().

---

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

• EvtDecayTag.hh
• EvtDecayTag.cc