rtang/PtolemyGUI
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Compare commits

base: rtang:b92e795919c10e83335ebdf5f115b9cf8a3524ad
Branches Tags
rtang:master
..
compare: rtang:939b63a79f354795f6202c1eed0619f48504d579
Branches Tags
rtang:master

No commits in common. "b92e795919c10e83335ebdf5f115b9cf8a3524ad" and "939b63a79f354795f6202c1eed0619f48504d579" have entirely different histories.
b92e795919 ... 939b63a79f

1 changed files with 16 additions and 69 deletions
Show Stats Expand all files Collapse all files

85
Cleopatra/InFileCreator.h
View File

@ -47,19 +47,6 @@ int GetLValue(string spdf){
return -1; return -1;
} }
float GetSpinValue(string spin){
size_t slashPos = spin.find('/');
if (slashPos != string::npos) {
int num = atoi(spin.substr(0, slashPos).c_str());
int den = atoi(spin.substr(slashPos + 1).c_str());
if (den != 0) return static_cast<float>(num) / static_cast<float>(den);
} else {
return atof(spin.c_str());
}
return -1.0;
}
int InFileCreator(string readFile, string infile, double angMin, double angMax, double angStep) { int InFileCreator(string readFile, string infile, double angMin, double angMax, double angStep) {
@ -132,20 +119,12 @@ int InFileCreator(string readFile, string infile, double angMin, double angMax,
if( iso_a.A == 6 && iso_a.Z == 3 && iso_b.A == 2 && iso_b.Z == 1 ) isReactionSupported = true; if( iso_a.A == 6 && iso_a.Z == 3 && iso_b.A == 2 && iso_b.Z == 1 ) isReactionSupported = true;
if( isReactionSupported == false ){ if( isReactionSupported == false ){
printf("\033[31m ===> Ignored. Reaction type not supported. \033[0m\n"); printf(" ===> Ignored. Reaction type not supported. \n");
continue; continue;
} }
// Continues to decode the input string // Continues to decode the input string
string gsSpinA = str0[1]; string gsSpinA = str0[1];
string gsSpinA = gsSpinparityA.substr(0, gsSpinparityA.length()-1);
string gsParityA = gsSpinparityA.substr(gsSpinparityA.length()-1);
// determine numeric ground-state spin value: half-integer if contains '/', otherwise integer
double gsSpinAValue = GetSpinValue(gsSpinA);
int gsParityAValue = (gsParityA == "+" ? 1 : -1);
string orbital = str0[2]; string orbital = str0[2];
string spinParity = str0[3]; string spinParity = str0[3];
@ -153,9 +132,6 @@ int InFileCreator(string readFile, string infile, double angMin, double angMax,
string spin = spinParity.substr(0, lenSpinParity-1); string spin = spinParity.substr(0, lenSpinParity-1);
string parity = spinParity.substr(lenSpinParity-1); string parity = spinParity.substr(lenSpinParity-1);
double spinValue = GetSpinValue(spin);
int parityValue = (parity == "+" ? 1 : -1);
string Ex = str0[4]; string Ex = str0[4];
string reactionEnergy = str0[5]; string reactionEnergy = str0[5];
string potential = str0[6]; string potential = str0[6];
@ -169,26 +145,26 @@ int InFileCreator(string readFile, string infile, double angMin, double angMax,
/// check is iso_A or iso_B exist in the mass table /// check is iso_A or iso_B exist in the mass table
if( iso_A.Mass == -404 || iso_B.Mass == -404 ){ if( iso_A.Mass == -404 || iso_B.Mass == -404 ){
printf("\033[31m ===> Error! mass does not found. \033[0m\n"); printf(" ===> Error! mass does not found. \n");
continue; continue;
} }
/// check reaction valid by balancing the A and Z number; /// check reaction valid by balancing the A and Z number;
if( iso_A.A + iso_a.A != iso_B.A + iso_b.A || iso_A.Z + iso_a.Z != iso_B.Z + iso_b.Z ) { if( iso_A.A + iso_a.A != iso_B.A + iso_b.A || iso_A.Z + iso_a.Z != iso_B.Z + iso_b.Z ) {
printf("\033[31m====> ERROR! A-number or Z-number not balanced. \033[0m\n"); printf("====> ERROR! A-number or Z-number not balanced. \n");
Isotope isotopeK(iso_A.A + iso_a.A - iso_b.A, iso_A.Z + iso_a.Z - iso_b.Z); Isotope isotopeK(iso_A.A + iso_a.A - iso_b.A, iso_A.Z + iso_a.Z - iso_b.Z);
printf(" try : %s(%s,%s)%s ??\n", iso_A.Name.c_str(), iso_a.Name.c_str(), iso_b.Name.c_str(), isotopeK.Name.c_str()); printf(" try : %s(%s,%s)%s ??\n", iso_A.Name.c_str(), iso_a.Name.c_str(), iso_b.Name.c_str(), isotopeK.Name.c_str());
continue; continue;
} }
if( isTransferReaction && potential.length() != 2 ){ if( isTransferReaction && potential.length() != 2 ){
printf("\033[31m====> ERROR! Potential input should be 2 charaters! skipped. \033[0m\n"); printf("====> ERROR! Potential input should be 2 charaters! skipped. \n");
continue; continue;
} }
string node ; string node ;
string jStr ; string jValue ;
string lStr ; string lValue ;
int spdf = -1; int spdf = -1;
if( isTransferReaction ) { if( isTransferReaction ) {
@ -197,46 +173,17 @@ int InFileCreator(string readFile, string infile, double angMin, double angMax,
// single nucleon transfer // single nucleon transfer
if( abs(iso_a.A - iso_b.A) == 1 ){ if( abs(iso_a.A - iso_b.A) == 1 ){
lStr = orbital.substr(1,1); lValue = orbital.substr(1,1);
jStr = orbital.substr(2); jValue = orbital.substr(2);
///printf(" l : %s, j : %s \n", lStr.c_str(), jStr.c_str()); ///printf(" l : %s, j : %s \n", lValue.c_str(), jValue.c_str());
spdf = GetLValue(lStr); spdf = GetLValue(lValue);
//Check parity, gsParityAValue * parityValue should be equal to (-1)^l
int expectedParity = gsParityAValue * parityValue;
int calculatedParity = ( (spdf %2) ==0 ) ? 1 : -1;
if( expectedParity != calculatedParity ){
printf("\033[31m ===> skipped. Parity mismatch. (expected %d but got %d) \033[0m\n", expectedParity, calculatedParity);
continue;
}
//Check converation of angular momentum, spinValue should be in |gsSpinAValue - jValue| to |gsSpinAValue + jValue|
float jValue = GetSpinValue(jStr);
//generate all possible spinValue
std::vector<float> possibleSpinValues;
for (float s = fabs(gsSpinAValue - jValue); s <= (gsSpinAValue + jValue); s += 1.0f) {
possibleSpinValues.push_back(s);
}
bool spinMatch = false;
for (float possibleSpin : possibleSpinValues) {
if (fabs(possibleSpin - spinValue) < 1e-6) {
spinMatch = true;
break;
}
}
if (!spinMatch) {
printf("\033[31m ===> skipped. Angular momentum mismatch. \033[0m\n");
continue;
}
} }
// two-nucleons transfer, v vc // two-nucleons transfer, v vc
if( abs(iso_a.A - iso_b.A) == 2 || abs(iso_a.A - iso_b.A) == 4 ){ if( abs(iso_a.A - iso_b.A) == 2 || abs(iso_a.A - iso_b.A) == 4 ){
size_t posEq = orbital.find('='); size_t posEq = orbital.find('=');
lStr = orbital.substr(posEq+1,1); lValue = orbital.substr(posEq+1,1);
spdf=atoi(lStr.c_str()); spdf=atoi(lValue.c_str());
} }
// (6Li,d) // (6Li,d)
@ -245,11 +192,11 @@ int InFileCreator(string readFile, string infile, double angMin, double angMax,
} }
if( abs(iso_a.A - iso_b.A) == 0 ){ if( abs(iso_a.A - iso_b.A) == 0 ){
printf("\033[31m ===> skipped. p-n exchange reaction does not support. \033[0m\n"); printf(" ===? skipped. p-n exchange reaction does not support. \n");
} }
if( spdf == -1 ){ if( spdf == -1 ){
printf("\033[31m ===> skipped. Not reconginzed orbital-label. (user input : l=%s | %s) \033[0m\n", lStr.c_str(), orbital.c_str()); printf(" ===> skipped. Not reconginzed orbital-label. (user input : l=%s | %s) \n", lValue.c_str(), orbital.c_str());
continue; continue;
} }
} }
@ -406,9 +353,9 @@ int InFileCreator(string readFile, string infile, double angMin, double angMax,
if( abs(iso_a.A-iso_b.A) == 1 ){ if( abs(iso_a.A-iso_b.A) == 1 ){
fprintf(file_out, "JBIGA=%s\n", gsSpinA.c_str()); fprintf(file_out, "JBIGA=%s\n", gsSpinA.c_str());
if( isAboveThreshold ) { if( isAboveThreshold ) {
fprintf(file_out, "nodes=%s l=%d jp=%s E=-.2 $node is n-1, set binding 200 keV \n", node.c_str(), spdf, jStr.c_str()); fprintf(file_out, "nodes=%s l=%d jp=%s E=-.2 $node is n-1, set binding 200 keV \n", node.c_str(), spdf, jValue.c_str());
}else{ }else{
fprintf(file_out, "nodes=%s l=%d jp=%s $node is n-1 \n", node.c_str(), spdf, jStr.c_str()); fprintf(file_out, "nodes=%s l=%d jp=%s $node is n-1 \n", node.c_str(), spdf, jValue.c_str());
} }
fprintf(file_out, "r0=1.25 a=.65 \n"); fprintf(file_out, "r0=1.25 a=.65 \n");
fprintf(file_out, "vso=6 rso0=1.10 aso=.65 \n"); fprintf(file_out, "vso=6 rso0=1.10 aso=.65 \n");