Migrated observable phi2.

[phys/u1casc.git] / u1casc-ordinary / u1casc-ordinary.cpp

diff --git a/u1casc-ordinary/u1casc-ordinary.cpp b/u1casc-ordinary/u1casc-ordinary.cpp
index f6f05fd13cf91aebcdbdd806c2b92547f4590ad3..4956ee1b304e1fd2d9366eb24a0557570f1a505b 100644 (file)
--- a/u1casc-ordinary/u1casc-ordinary.cpp
+++ b/u1casc-ordinary/u1casc-ordinary.cpp
@@ -10,14 +10,16 @@ sim *Sim;
 
 //const complex<double> _i_ = complex<double>(0.0,1.0);
 
 
 //const complex<double> _i_ = complex<double>(0.0,1.0);
 

-//#include "obs_phi2.hpp"
+#include "obs_phi2.hpp"

 
 o815::comoption specOps[] = {
 
 o815::comoption specOps[] = {

-  { "kappaone", required_argument, NULL, 'k', "set inverse mass kappa_1", "min:max:inc" },
-  { "kappatwo", required_argument, NULL, 'm', "set inverse mass kappa_2", "min:max:inc" },
+  { "kappaone", required_argument, NULL, 'r', "set inverse mass kappa_1", "min:max:inc" },
+  { "kappatwo", required_argument, NULL, 's', "set inverse mass kappa_2", "min:max:inc" },

   { "beta", required_argument, NULL, 'b', "set inverse gauge coupling beta", "min:max:inc" },
   { "lambdaone", required_argument, NULL, 'p', "set quartic coupling lambda_1", "min:max:inc" },
   { "lambdatwo", required_argument, NULL, 'q', "set quartic coupling lambda_2", "min:max:inc" },
   { "beta", required_argument, NULL, 'b', "set inverse gauge coupling beta", "min:max:inc" },
   { "lambdaone", required_argument, NULL, 'p', "set quartic coupling lambda_1", "min:max:inc" },
   { "lambdatwo", required_argument, NULL, 'q', "set quartic coupling lambda_2", "min:max:inc" },

+  { "kappa", required_argument, NULL, 'k', "set kappa_1 = kappa_2 = kappa", "min:max:inc" },
+  { "lambda", required_argument, NULL, 'l', "set lambda_1 = lambda_2 = lambda", "min:max:inc" },

   { "", 0, NULL, 0, "", "" }
 };
 
   { "", 0, NULL, 0, "", "" }
 };
 


@@ -25,8 +27,28 @@ void parseSpecOps()
 {
   for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
     switch(O815->parsedSpecOps[isopt].first) {
 {
   for (int isopt = 0; isopt < O815->parsedSpecOps.size(); isopt++)
     switch(O815->parsedSpecOps[isopt].first) {

-    case 'm':
-      O815->paraQ->addRange("mass", O815->parsedSpecOps[isopt].second);
+    case 'r':
+      O815->paraQ->addRange("kappaone", O815->parsedSpecOps[isopt].second);
+      break;
+    case 's':
+      O815->paraQ->addRange("kappatwo", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'b':
+      O815->paraQ->addRange("beta", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'p':
+      O815->paraQ->addRange("lambdaone", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'q':
+      O815->paraQ->addRange("lambdatwo", O815->parsedSpecOps[isopt].second);
+      break;
+    case 'l':
+      O815->paraQ->addRange("lambdaone", O815->parsedSpecOps[isopt].second);
+      O815->paraQ->linkParas("lambdatwo", "lambdaone");
+      break;
+    case 'k':
+      O815->paraQ->addRange("kappaone", O815->parsedSpecOps[isopt].second);
+      O815->paraQ->linkParas("kappatwo", "kappaone");

       break;
     }
 }
       break;
     }
 }


@@ -39,11 +61,11 @@ void helpHeader()
 void parseLonelyArgs()
 {
   for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {
 void parseLonelyArgs()
 {
   for (vector<char*>::iterator lonit = O815->lonelyArgs.begin(); lonit != O815->lonelyArgs.end(); ++lonit) {

-    /*

     if ( strcmp(*lonit, "phi2") == 0 ) {
       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
       O815->observables.push_back(new obs_phi2(O815));
     }
     if ( strcmp(*lonit, "phi2") == 0 ) {
       *O815->out->log << "MASTER: registered observable: phi2" << endl << flush;
       O815->observables.push_back(new obs_phi2(O815));
     }

+    /*

     else if ( strcmp(*lonit, "phi2_hist") == 0 ) {
       *O815->out->log << "MASTER: registered observable: phi2_hist" << endl << flush;
       O815->observables.push_back(new obs_phi2_hist(O815));
     else if ( strcmp(*lonit, "phi2_hist") == 0 ) {
       *O815->out->log << "MASTER: registered observable: phi2_hist" << endl << flush;
       O815->observables.push_back(new obs_phi2_hist(O815));


@@ -63,7 +85,7 @@ int main (int argc, char *argv[])
   O815->addPara("beta", 1);
   O815->addPara("lambdaone", 1);
   O815->addPara("lambdatwo", 1);
   O815->addPara("beta", 1);
   O815->addPara("lambdaone", 1);
   O815->addPara("lambdatwo", 1);

-  O815->addPara("kappaone", 1);
+  O815->addPara("kappaone", 1);  

   O815->addPara("kappatwo", 1);
 
   parseSpecOps();
   O815->addPara("kappatwo", 1);
 
   parseSpecOps();