EvtPHOTOS Class Reference

#include <EvtPHOTOS.hh>

Inheritance diagram for EvtPHOTOS:

Public Member Functions
void	doRadCorr (EvtParticle *p)

Detailed Description

Definition at line 30 of file EvtPHOTOS.hh.

Member Function Documentation

doRadCorr()

void EvtPHOTOS::doRadCorr ( EvtParticle * p )

virtual

Implements EvtAbsRadCorr.

Definition at line 44 of file EvtPHOTOS.cc.

                                          {
 
  static int first = 1;
 
  // added by Lange Jan4,2000
  static EvtId GAMM = EvtPDL::getId( "gammaFSR" ); // pingrg, 2009/02/19
 
  if ( first )
  {
 
    first = 0;
    phoini_();
  }
 
  int entry, eventnum, numparticle, istat, partnum, mother;
  int daugfirst, dauglast;
 
  int numparticlephotos;
 
  double px, py, pz, e, m, x, y, z, t;
 
  static EvtId    dq  = EvtPDL::getId( "d" );
  static EvtId    adq = EvtPDL::getId( "anti-d" );
  static EvtId    uq  = EvtPDL::getId( "u" );
  static EvtId    auq = EvtPDL::getId( "anti-u" );
  static EvtId    sq  = EvtPDL::getId( "s" );
  static EvtId    asq = EvtPDL::getId( "anti-s" );
  static EvtId    cq  = EvtPDL::getId( "c" );
  static EvtId    acq = EvtPDL::getId( "anti-c" );
  static EvtId    bq  = EvtPDL::getId( "b" );
  static EvtId    abq = EvtPDL::getId( "anti-b" );
  static EvtId    tq  = EvtPDL::getId( "t" );
  static EvtId    atq = EvtPDL::getId( "anti-t" );
  static EvtIdSet quarks( dq, adq, uq, auq, sq, asq, cq, acq, bq, abq, tq, atq );
 
  px = 0.0;
  py = 0.0;
  pz = 0.0;
  e  = p->mass();
  m  = p->mass();
  x  = 0.0;
  y  = 0.0;
  z  = 0.0;
  t  = 0.0;
 
  entry       = 1;
  eventnum    = 1;
  numparticle = 1;
  istat       = 2;
  partnum     = EvtPDL::getStdHep( p->getId() );
  mother      = 0;
  daugfirst   = 2;
  dauglast    = 1 + p->getNDaug();
 
  begevtgenstorex_( &entry, &eventnum, &numparticle, &istat, &partnum, &mother, &daugfirst,
                    &dauglast, &px, &py, &pz, &e, &m, &x, &y, &z, &t );
 
  int i;
  //  std::cout << EvtPDL::name(p->getId()) <<" " ;
  for ( i = 0; i < p->getNDaug(); i++ )
  {
 
    // No quarks to photos
    if ( quarks.contains( p->getDaug( i )->getId() ) == 1 ) continue;
 
    px = p->getDaug( i )->getP4().get( 1 );
    py = p->getDaug( i )->getP4().get( 2 );
    pz = p->getDaug( i )->getP4().get( 3 );
    e  = p->getDaug( i )->getP4().get( 0 );
    m  = p->getDaug( i )->mass();
    x  = 0.0;
    y  = 0.0;
    z  = 0.0;
    t  = 0.0;
 
    //    std::cout << EvtPDL::name(p->getDaug(i)->getId()) << " " ;
    entry += 1;
    eventnum = 1;
    numparticle += 1;
    istat     = 1;
    partnum   = EvtPDL::getStdHep( p->getDaug( i )->getId() );
    mother    = 1;
    daugfirst = 0;
    dauglast  = 0;
 
    begevtgenstorex_( &entry, &eventnum, &numparticle, &istat, &partnum, &mother, &daugfirst,
                      &dauglast, &px, &py, &pz, &e, &m, &x, &y, &z, &t );
  }
  //  std::cout << std::endl;
  // can't use heplst since the common block used by the BaBar
  // implementation of PHOTOS  is renamed due to real*4 vs real*8
  // problems.
 
  // int mlst=1;
 
  // heplst_(&mlst);
 
  entry = 1;
 
  //  report(INFO,"EvtGen") << "Doing photos " << EvtPDL::name(p->getId()) << endl;
  photos_( &entry );
  //  report(INFO,"EvtGen") << "done\n";
  begevtgengetx_( &entry, &eventnum, &numparticlephotos, &istat, &partnum, &mother, &daugfirst,
                  &dauglast, &px, &py, &pz, &e, &m, &x, &y, &z, &t );
 
  // report(INFO,"EvtGen") << "numparticlephotos:"<<numparticlephotos<<endl;
 
  if ( numparticle == numparticlephotos ) return;
 
  EvtVector4R new4mom;
 
  int np;
 
  for ( i = 0; i < p->getNDaug(); i++ )
  {
 
    entry = i + 2;
 
    begevtgengetx_( &entry, &eventnum, &np, &istat, &partnum, &mother, &daugfirst, &dauglast,
                    &px, &py, &pz, &e, &m, &x, &y, &z, &t );
 
    // this is needed to ensure that photos does not
    // change the masses. But it will violate energy conservation!
    double mp = p->getDaug( i )->mass();
    e         = sqrt( mp * mp + px * px + py * py + pz * pz );
 
    new4mom.set( e, px, py, pz );
 
    p->getDaug( i )->setP4( new4mom );
  }
 
  for ( entry = numparticle + 1; entry <= numparticlephotos; entry++ )
  {
 
    begevtgengetx_( &entry, &eventnum, &np, &istat, &partnum, &mother, &daugfirst, &dauglast,
                    &px, &py, &pz, &e, &m, &x, &y, &z, &t );
 
    new4mom.set( e, px, py, pz );
 
    // new4mom.dump();
 
    EvtPhotonParticle* gamma;
    gamma = new EvtPhotonParticle;
    gamma->init( GAMM, new4mom );
    //    report(INFO,"EvtGen") << gamma << " " << p << " "<< px << " " << py << " " << pz << "
    //    " << p->getNDaug() << " " << EvtPDL::name(p->getId())<<" " << entry << " "
    //    <<numparticlephotos<<endl;
    gamma->addDaug( p );
 
    //    p->getDaug(i)->set_type(EvtSpinType::PHOTON);
  }
  return;
}

---

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

• EvtPHOTOS.hh
• EvtPHOTOS.cc