Skip to content
Snippets Groups Projects
Commit 91951929 authored by Andrew Haas's avatar Andrew Haas Committed by Graeme Stewart
Browse files

'comment out couts' (G4ExternalDecay-00-00-08-01)

parent c838e83a
No related branches found
No related tags found
No related merge requests found
...@@ -94,7 +94,7 @@ void PythiaForDecays::Py1ent(int ip, int kf, double pe, double theta, double phi ...@@ -94,7 +94,7 @@ void PythiaForDecays::Py1ent(int ip, int kf, double pe, double theta, double phi
} }
void PythiaForDecays::DecayRhadrons(int pdgid){ void PythiaForDecays::DecayRhadrons(int pdgid){
std::cout<<"ACH start DecayRhadrons for pdgid "<<pdgid<<std::endl; //std::cout<<"ACH start DecayRhadrons for pdgid "<<pdgid<<std::endl;
int pyl=3;//level to list to output text int pyl=3;//level to list to output text
//only decay daughter particles from the decay that are "near" the hadron //only decay daughter particles from the decay that are "near" the hadron
...@@ -105,7 +105,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){ ...@@ -105,7 +105,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){
//Deal with standard hadrons first... //Deal with standard hadrons first...
if (abs(pdgid)<1000000){ if (abs(pdgid)<1000000){
std::cout<<"ACH decay standard hadron for pdgid "<<pdgid<<std::endl; //std::cout<<"ACH decay standard hadron for pdgid "<<pdgid<<std::endl;
pyjets_.k[0][0]=5;//make it decay right away! pyjets_.k[0][0]=5;//make it decay right away!
//pylist_(&pyl); //pylist_(&pyl);
pyexec_(); pyexec_();
...@@ -118,7 +118,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){ ...@@ -118,7 +118,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){
pyjets_.k[0][0]=6; pyjets_.k[0][0]=6;
//pylist_(&pyl); //pylist_(&pyl);
if ( pdgid/10000==109 || pdgid/1000==1009 || pdgid/100==10009 ){ if ( pdgid/10000==109 || pdgid/1000==1009 || pdgid/100==10009 ){
std::cout<<"ACH decay gluino rhadron .."<<std::endl; //std::cout<<"ACH decay gluino rhadron .."<<std::endl;
Pydat3 m_pydat3; Pydat3 m_pydat3;
int kfgl = 1000021;//gluino int kfgl = 1000021;//gluino
m_pydat3.mdcy(pycomp_(&kfgl),1)=1;// MDCY(PYCOMP(KFGL),1)=1 //unstable gluino m_pydat3.mdcy(pycomp_(&kfgl),1)=1;// MDCY(PYCOMP(KFGL),1)=1 //unstable gluino
...@@ -126,7 +126,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){ ...@@ -126,7 +126,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){
//pylist_(&pyl); //pylist_(&pyl);
} }
else if ( pdgid/1000==1006 || pdgid/100==10006 ){ else if ( pdgid/1000==1006 || pdgid/100==10006 ){
std::cout<<"ACH decay stop rhadron .."<<std::endl; //std::cout<<"ACH decay stop rhadron .."<<std::endl;
Pydat3 m_pydat3; Pydat3 m_pydat3;
int kfgl = 1000006;//stop int kfgl = 1000006;//stop
m_pydat3.mdcy(pycomp_(&kfgl),1)=1;// MDCY(PYCOMP(KFGL),1)=1 //unstable stop m_pydat3.mdcy(pycomp_(&kfgl),1)=1;// MDCY(PYCOMP(KFGL),1)=1 //unstable stop
...@@ -135,7 +135,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){ ...@@ -135,7 +135,7 @@ void PythiaForDecays::DecayRhadrons(int pdgid){
//pylist_(&pyl); //pylist_(&pyl);
} }
else if ( pdgid/1000==1005 || pdgid/100==10005 ){ else if ( pdgid/1000==1005 || pdgid/100==10005 ){
std::cout<<"ACH decay sbottom rhadron .."<<std::endl; //std::cout<<"ACH decay sbottom rhadron .."<<std::endl;
Pydat3 m_pydat3; Pydat3 m_pydat3;
int kfgl = 1000005;//sbottom int kfgl = 1000005;//sbottom
m_pydat3.mdcy(pycomp_(&kfgl),1)=1;// MDCY(PYCOMP(KFGL),1)=1 //unstable sbottom m_pydat3.mdcy(pycomp_(&kfgl),1)=1;// MDCY(PYCOMP(KFGL),1)=1 //unstable sbottom
...@@ -167,6 +167,8 @@ void PythiaForDecays::ImportParticles( std::vector<G4DynamicParticle*> & particl ...@@ -167,6 +167,8 @@ void PythiaForDecays::ImportParticles( std::vector<G4DynamicParticle*> & particl
for (HepMC::GenEvent::particle_const_iterator p=evt->particles_begin(); p!=evt->particles_end(); ++p){ for (HepMC::GenEvent::particle_const_iterator p=evt->particles_begin(); p!=evt->particles_end(); ++p){
if (! (*p) ) continue; if (! (*p) ) continue;
/*
double r=0; double r=0;
if ((*p)->production_vertex()) { if ((*p)->production_vertex()) {
double x=(*p)->production_vertex()->position().x(); double x=(*p)->production_vertex()->position().x();
...@@ -174,6 +176,8 @@ void PythiaForDecays::ImportParticles( std::vector<G4DynamicParticle*> & particl ...@@ -174,6 +176,8 @@ void PythiaForDecays::ImportParticles( std::vector<G4DynamicParticle*> & particl
r=sqrt(x*x+y*y); r=sqrt(x*x+y*y);
} }
std::cout<<"ImportParticles: "<<(*p)->pdg_id()<<", vertex r "<<r<<", status "<<(*p)->status()<<std::endl; std::cout<<"ImportParticles: "<<(*p)->pdg_id()<<", vertex r "<<r<<", status "<<(*p)->status()<<std::endl;
*/
if ( (*p)->status()!=1 ) continue; // stable only if ( (*p)->status()!=1 ) continue; // stable only
//if ( std::abs((*p)->pdg_id())==12 || std::abs((*p)->pdg_id())==14 || std::abs((*p)->pdg_id())==16 ) continue; // skip neutrinos //if ( std::abs((*p)->pdg_id())==12 || std::abs((*p)->pdg_id())==14 || std::abs((*p)->pdg_id())==16 ) continue; // skip neutrinos
G4ThreeVector momentum( (*p)->momentum().x() , (*p)->momentum().y() , (*p)->momentum().z() ); G4ThreeVector momentum( (*p)->momentum().x() , (*p)->momentum().y() , (*p)->momentum().z() );
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment