EvtPHOTOS Class Reference

#include <EvtPHOTOS.hh>

Inheritance diagram for EvtPHOTOS:

Public Member Functions
void	doRadCorr (EvtParticle *p)
virtual void	doRadCorr (EvtParticle *p)=0

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 50 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