Merge pull request #21 from scarrazza/integrability

Integrability
scarrazza · Jun 17, 2020 · fdfd068 · fdfd068
2 parents a6bd3c1 + c1d751f
commit fdfd068
Show file tree

Hide file tree

Showing 3 changed files with 83 additions and 1 deletion.
diff --git a/src/DIS/FKObservables.f b/src/DIS/FKObservables.f
@@ -91,6 +91,51 @@ function FKObservables(x,Q,y)
       elseif(obs(1:7).eq."DIS_XGL")then
          FKObservables = xPDF(0,x)
 *
+****  Non-singlet T3
+*
+      elseif(obs(1:7).eq."DIS_XT3")then
+         FKObservables = xPDF(+2,x) + xpdf(-2,x)
+     1        - 1d0 * ( xPDF(+1,x) + xPDF(-1,x) )
+*
+****  Non-singlet T8
+*
+      elseif(obs(1:7).eq."DIS_XT8")then
+         FKObservables = xPDF(+2,x) + xpdf(-2,x)
+     1        + 1d0 * ( xPDF(+1,x) + xPDF(-1,x) )
+     1        - 2d0 * ( xPDF(+3,x) + xPDF(-3,x) )
+*
+****  Non-singlet T15
+*
+      elseif(obs(1:8).eq."DIS_XT15")then
+         FKObservables = xPDF(+2,x) + xpdf(-2,x)
+     1        + 1d0 * ( xPDF(+1,x) + xPDF(-1,x) )
+     1        + 1d0 * ( xPDF(+3,x) + xPDF(-3,x) )
+     1        - 3d0 * ( xPDF(+4,x) + xPDF(-4,x) )
+*
+****  Valence V3
+*
+      elseif(obs(1:7).eq."DIS_XV3")then
+         FKObservables = xPDF(+2,x) - xPDF(-2,x) 
+     1        - 1d0 * ( xPDF(+1,x) - xPDF(-1,x) )
+*
+****  Valence V8
+*
+      elseif(obs(1:7).eq."DIS_XV8")then
+         FKObservables = xPDF(+2,x) - xPDF(-2,x) 
+     1        + 1d0 * ( xPDF(+1,x) - xPDF(-1,x) )
+     2        - 2d0 * ( xPDF(+3,x) - xPDF(-3,x) )
+*
+*
+****  Valence V
+*
+      elseif(obs(1:6).eq."DIS_XV")then
+         FKObservables = xPDF(+2,x) - xPDF(-2,x) 
+     1        + xPDF(+1,x) - xPDF(-1,x)
+     1        + xPDF(+3,x) - xPDF(-3,x)
+     1        + xPDF(+4,x) - xPDF(-4,x) 
+     1        + xPDF(+5,x) - xPDF(-5,x)
+     1        + xPDF(+6,x) - xPDF(-6,x) 
+*
 ****  Light structure function F2light
 *
       elseif(obs(1:7).eq."DIS_F2L")then         

diff --git a/src/DIS/FKSimulator.f b/src/DIS/FKSimulator.f
@@ -90,6 +90,37 @@ function FKSimulator(x,Q,y,i,beta)
       elseif(obs(1:7).eq."DIS_XGL")then
          FKSimulator = ExternalEvolutionOperator("Ev2Ph",0,i,x,beta)
 *
+****  Non-singlet T3
+*
+      elseif(obs(1:7).eq."DIS_XT3")then
+         FKSimulator = ExternalEvolutionOperator("Ev2Ev",9,i,x,beta)
+*
+****  Non-singlet T8
+*
+      elseif(obs(1:7).eq."DIS_XT8")then
+         FKSimulator = ExternalEvolutionOperator("Ev2Ev",10,i,x,beta)
+*
+****  Non-singlet T15
+*
+      elseif(obs(1:8).eq."DIS_XT15")then
+         FKSimulator = ExternalEvolutionOperator("Ev2Ev",11,i,x,beta)
+*
+****  Valence V3
+*
+      elseif(obs(1:7).eq."DIS_XV3")then
+         FKSimulator = ExternalEvolutionOperator("Ev2Ev",4,i,x,beta)
+*
+****  Valence V8
+*
+      elseif(obs(1:7).eq."DIS_XV8")then
+         FKSimulator = ExternalEvolutionOperator("Ev2Ev",5,i,x,beta)
+*
+*
+****  Valence V
+*
+      elseif(obs(1:6).eq."DIS_XV")then
+         FKSimulator = ExternalEvolutionOperator("Ev2Ev",3,i,x,beta)
+*
 ****  Light structure function F2light
 *
       elseif(obs(1:7).eq."DIS_F2L")then

diff --git a/src/DIS/SetFKObservable.f b/src/DIS/SetFKObservable.f
@@ -27,7 +27,13 @@ subroutine SetFKObservable(obs)
      4     obs(1:7).eq."DIS_XSQ".or.
      5     obs(1:7).eq."DIS_XSB".or.
      6     obs(1:7).eq."DIS_XCQ".or.
-     7     obs(1:7).eq."DIS_XGL")then
+     7     obs(1:7).eq."DIS_XGL".or.
+     8     obs(1:7).eq."DIS_XT3".or.
+     9     obs(1:7).eq."DIS_XT8".or.
+     1     obs(1:8).eq."DIS_XT15".or.
+     2     obs(1:7).eq."DIS_XV3".or.
+     3     obs(1:7).eq."DIS_XV8".or.
+     4     obs(1:6).eq."DIS_XV")then
          call SetTargetDIS("proton") 
 *
 ****  Light structure function F2light