Page MenuHomeHEPForge
Files F19244206

EvtPDL.cpp
No OneTemporary
Actions

Size
9 KB
Referenced Files
None
Subscribers
None

EvtPDL.cpp
View Options

//--------------------------------------------------------------------------
//
// Environment:
// This software is part of the EvtGen package developed jointly
// for the BaBar and CLEO collaborations. If you use all or part
// of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
// Copyright (C) 1998 Caltech, UCSB
//
// Module: EvtPDL.cc
//
// Description: routines to store particle properties in EvtPDL structure.
//
// Modification history:
//
// DJL/RYD September 25, 1996 Module created
//
//------------------------------------------------------------------------
//
#include "EvtGenBase/EvtPatches.hh"
#include <iostream>
#include <fstream>
#include <stdlib.h>
#include <ctype.h>
#include <string.h>
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtPartProp.hh"
#include "EvtGenBase/EvtId.hh"
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtReport.hh"
using std::endl;
using std::fstream;
using std::ifstream;
static int first=1;
unsigned int EvtPDL::_firstAlias;
int EvtPDL::_nentries;
std::map<std::string,int> EvtPDL::_particleNameLookup;
EvtPDL::EvtPDL() {
if (first!=0) {
first=0;
_nentries=0;
_firstAlias=999999;
}
}
EvtPDL::~EvtPDL(){
}
void EvtPDL::read(const char* fname)
{
readPDT(fname);
}
void EvtPDL::readPDT(const std::string fname){
ifstream indec;
indec.open(fname.c_str());
char cmnd[100];
char xxxx[100];
char pname[100];
int stdhepid;
double mass;
double pwidth;
double pmaxwidth;
int chg3;
int spin2;
double ctau;
int lundkc;
EvtId i;
if (!indec) {
EvtGenReport(EVTGEN_ERROR,"EvtGen") << "Could not open:"<<fname.c_str()<<"EvtPDL"<<endl;
return;
}
do{
char ch,ch1;
do{
indec.get(ch);
if (ch=='\n') indec.get(ch);
if (ch!='*') {
indec.putback(ch);
}
else{
while (indec.get(ch1),ch1!='\n');
}
} while(ch=='*');
indec >> cmnd;
if (strcmp(cmnd,"end")){
if (!strcmp(cmnd,"add")){
indec >> xxxx;
indec >> xxxx;
indec >> pname;
indec >> stdhepid;
indec >> mass;
indec >> pwidth;
indec >> pmaxwidth;
indec >> chg3;
indec >> spin2;
indec >> ctau;
indec >> lundkc;
i=EvtId(_nentries,_nentries);
EvtPartProp tmp;
tmp.setSpinType(EvtSpinType::SCALAR);
if (spin2==0) tmp.setSpinType(EvtSpinType::SCALAR);
if (spin2==1) tmp.setSpinType(EvtSpinType::DIRAC);
if (spin2==2) tmp.setSpinType(EvtSpinType::VECTOR);
if (spin2==3) tmp.setSpinType(EvtSpinType::RARITASCHWINGER);
if (spin2==4) tmp.setSpinType(EvtSpinType::TENSOR);
if (spin2==5) tmp.setSpinType(EvtSpinType::SPIN5HALF);
if (spin2==6) tmp.setSpinType(EvtSpinType::SPIN3);
if (spin2==7) tmp.setSpinType(EvtSpinType::SPIN7HALF);
if (spin2==8) tmp.setSpinType(EvtSpinType::SPIN4);
if (spin2==2 && mass < 0.0001 ) tmp.setSpinType(EvtSpinType::PHOTON);
if (spin2==1 && mass < 0.0001 ) tmp.setSpinType(EvtSpinType::NEUTRINO);
if (!strcmp(pname,"string")){
tmp.setSpinType(EvtSpinType::STRING);
}
if (!strcmp(pname,"vpho")){
tmp.setSpinType(EvtSpinType::VECTOR);
}
tmp.setId(i);
tmp.setIdChgConj(EvtId(-1,-1));
tmp.setStdHep(stdhepid);
tmp.setLundKC(lundkc);
tmp.setName(pname);
if (_particleNameLookup.find(std::string(pname))!=
_particleNameLookup.end()) {
EvtGenReport(EVTGEN_ERROR,"EvtGen")<<"The particle name:"<<pname<<" is already defined."<<endl;
EvtGenReport(EVTGEN_ERROR,"EvtGen") << "Will terminate execution.";
::abort();
}
_particleNameLookup[std::string(pname)]=_nentries;
tmp.setctau(ctau);
tmp.setChg3(chg3);
tmp.initLineShape(mass,pwidth,pmaxwidth);
partlist().push_back(tmp);
_nentries++;
}
// if find a set read information and discard it
if (!strcmp(cmnd,"set")){
indec >> xxxx;
indec >> xxxx;
indec >> xxxx;
indec >> xxxx;
}
}
}while(strcmp(cmnd,"end"));
setUpConstsPdt();
}
void EvtPDL::aliasChgConj(EvtId a,EvtId abar){
if (EvtPDL::chargeConj(EvtId(a.getId(),a.getId()))!=
EvtId(abar.getId(),abar.getId())) {
EvtGenReport(EVTGEN_ERROR,"EvtGen")<<"Can't charge conjugate the two aliases:"
<<EvtPDL::name(a).c_str()<<" and "<<EvtPDL::name(abar).c_str()<<endl;
::abort();
}
partlist()[a.getAlias()].setIdChgConj(abar);
partlist()[abar.getAlias()].setIdChgConj(a);
}
EvtId EvtPDL::chargeConj(EvtId id){
// EvtId idchg=partlist()[id.getAlias()].getIdChgConj();
int index=id.getAlias();
EvtId idchg;
if (index>-1) {
idchg=partlist()[id.getAlias()].getIdChgConj();
}
if (idchg!=EvtId(-1,-1)) return idchg;
if (id.getId()!=id.getAlias()){
if (chargeConj(EvtId(id.getId(),id.getId()))==EvtId(id.getId(),id.getId())){
partlist()[id.getAlias()].setIdChgConj(id);
return id;
}
}
if (id.getAlias()!=id.getId()) {
EvtGenReport(EVTGEN_ERROR,"EvtGen")<<"Trying to charge conjugate alias particle:"
<<name(id).c_str()<<" without defining the alias!"<<endl;
::abort();
}
for (size_t i=0;i<partlist().size();i++){
if (partlist()[i].getStdHep()==-partlist()[id.getId()].getStdHep()){
partlist()[id.getId()].setIdChgConj(partlist()[i].getId());
return partlist()[i].getId();
}
}
partlist()[id.getId()].setIdChgConj(id);
return id;
}
EvtId EvtPDL::evtIdFromStdHep(int stdhep){
for (size_t i=0;i<partlist().size();i++){
if (partlist()[i].getStdHep()==stdhep)
return partlist()[i].getId();
}
return EvtId(-1,-1);
}
void EvtPDL::alias(EvtId num,const std::string& newname){
if ( _firstAlias < partlist().size() ) {
for(size_t i=_firstAlias;i<partlist().size();i--){
if (newname==partlist()[i].getName()){
EvtGenReport(EVTGEN_WARNING,"EvtGen")<<"Redefining alias:"<<newname.c_str()<<" will be ignored!"<<endl;
return;
}
}
}
else{
_firstAlias=partlist().size();
}
partlist().push_back(partlist()[num.getId()]);
int entry=partlist().size()-1;
partlist()[entry].setName(newname);
if (_particleNameLookup.find(std::string(newname))!=
_particleNameLookup.end()){
EvtGenReport(EVTGEN_ERROR,"EvtGen")<<"The particle name:"<<newname<<" is already defined."<<endl;
EvtGenReport(EVTGEN_ERROR,"EvtGen") << "Will terminate execution.";
::abort();
}
_particleNameLookup[std::string(newname)]=entry;
partlist()[entry].setId(EvtId(num.getId(),entry));
//Lange - Dec7, 2003. Unset the charge conjugate.
partlist()[entry].setIdChgConj(EvtId(-1,-1));
}
EvtId EvtPDL::getId(const std::string& name ){
std::map<std::string,int>::iterator it=_particleNameLookup.find(std::string(name));
if (it==_particleNameLookup.end()) return EvtId(-1,-1);
return partlist()[it->second].getId();
}
void EvtPDL::setUpConstsPdt(){
}
// Function to get EvtId from LundKC ( == Pythia Hep Code , KF )
EvtId EvtPDL::evtIdFromLundKC(int pythiaId){
unsigned int i;
for (i=0;i<partlist().size();i++){
if (partlist()[i].getLundKC()==pythiaId)
return partlist()[i].getId();
}
return EvtId(-1,-1);
}
double EvtPDL::getMeanMass(EvtId i ) {
return partlist()[i.getId()].getMass();
}
double EvtPDL::getMass(EvtId i ) {
return partlist()[i.getId()].rollMass();
}
double EvtPDL::getRandMass(EvtId i, EvtId *parId, int nDaug, EvtId *dauId,
EvtId *othDaugId,double maxMass,
double *dauMasses ) {
return partlist()[i.getId()].getRandMass(parId,nDaug,dauId,
othDaugId,maxMass,dauMasses);
}
double EvtPDL::getMassProb(EvtId i, double mass, double massPar, int nDaug,
double *massDau) {
return partlist()[i.getId()].getMassProb(mass,massPar,nDaug,massDau);
}
double EvtPDL::getMaxMass(EvtId i ) {
return partlist()[i.getId()].getMassMax();
}
double EvtPDL::getMinMass(EvtId i ) {
return partlist()[i.getId()].getMassMin();
}
double EvtPDL::getMaxRange(EvtId i ) {
return partlist()[i.getId()].getMaxRange();
}
double EvtPDL::getWidth(EvtId i ) {
return partlist()[i.getId()].getWidth();
}
double EvtPDL::getctau(EvtId i ) {
return partlist()[i.getId()].getctau();
}
int EvtPDL::getStdHep(EvtId id ) {
return partlist()[id.getId()].getStdHep();
}
int EvtPDL::getLundKC(EvtId id ) {
return partlist()[id.getId()].getLundKC();
}
int EvtPDL::chg3(EvtId i ) {
return partlist()[i.getId()].getChg3();
}
EvtSpinType::spintype EvtPDL::getSpinType(EvtId i ) {
return partlist()[i.getId()].getSpinType();
}
std::string EvtPDL::name(EvtId i) {
return partlist()[i.getAlias()].getName();
}
size_t EvtPDL::entries() {
return partlist().size();
}
EvtId EvtPDL::getEntry(int i) {
return partlist()[i].getId();
}
void EvtPDL::reSetMass(EvtId i, double mass) {
partlist()[i.getId()].reSetMass(mass);
}
void EvtPDL::reSetWidth(EvtId i, double width) {
partlist()[i.getId()].reSetWidth(width);
}
void EvtPDL::reSetMassMin(EvtId i, double mass) {
partlist()[i.getId()].reSetMassMin(mass);
}
void EvtPDL::reSetMassMax(EvtId i,double mass) {
partlist()[i.getId()].reSetMassMax(mass);
}
void EvtPDL::reSetBlatt(EvtId i,double blatt) {
partlist()[i.getId()].reSetBlatt(blatt);
}
void EvtPDL::reSetBlattBirth(EvtId i,double blatt) {
partlist()[i.getId()].reSetBlattBirth(blatt);
}
void EvtPDL::includeBirthFactor(EvtId i,bool yesno) {
partlist()[i.getId()].includeBirthFactor(yesno);
}
void EvtPDL::includeDecayFactor(EvtId i,bool yesno) {
partlist()[i.getId()].includeDecayFactor(yesno);
}
void EvtPDL::changeLS(EvtId i, std::string &newLS ) {
partlist()[i.getId()].newLineShape(newLS);
}
void EvtPDL::setPWForDecay(EvtId i, int spin, EvtId d1, EvtId d2) {
partlist()[i.getId()].setPWForDecay(spin,d1,d2);
}
void EvtPDL::setPWForBirthL(EvtId i, int spin, EvtId par, EvtId othD) {
partlist()[i.getId()].setPWForBirthL(spin,par,othD);
}

File Metadata

Mime Type
text/x-c
Expires
Tue, Sep 30, 4:40 AM (15 h, 32 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6528986
Default Alt Text
EvtPDL.cpp (9 KB)

Event Timeline

HEPForge.org · .