From b9504e154143121b59b72e04ea3d1a478db5c79b Mon Sep 17 00:00:00 2001 From: Emanuele Nocera Date: Wed, 3 Jun 2020 11:55:07 +0200 Subject: [PATCH 1/5] Added integrability observables --- src/DIS/FKObservables.f | 45 +++++++++++++++++++++++++++++++++++++++ src/DIS/FKSimulator.f | 31 +++++++++++++++++++++++++++ src/DIS/SetFKObservable.f | 3 ++- 3 files changed, 78 insertions(+), 1 deletion(-) diff --git a/src/DIS/FKObservables.f b/src/DIS/FKObservables.f index 901bd2fd..98aad769 100644 --- 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:9).eq."INTEG_XT3")then + FKObservables = xPDF(+2,x) + xpdf(-2,x) + 1 - ( xPDF(+1,x) + xPDF(-1,x) ) +* +**** Non-singlet T8 +* + elseif(obs(1:9).eq."INTEG_XT8")then + FKObservables = xPDF(+2,x) + xpdf(-2,x) + 1 + ( xPDF(+1,x) + xPDF(-1,x) ) + 1 - 2d0 * ( xPDF(+3,x) + xPDF(-3,x) ) +* +**** Non-singlet T15 +* + elseif(obs(1:10).eq."INTEG_XT15")then + FKObservables = xPDF(+2,x) + xpdf(-2,x) + 1 + ( xPDF(+1,x) + xPDF(-1,x) ) + 1 + ( xPDF(+3,x) + xPDF(-3,x) ) + 1 - 3d0 * ( xPDF(+4,x) + xPDF(-4,x) ) +* +* +**** Valence V +* + elseif(obs(1:8).eq."INTEG_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) +* +**** Valence V3 +* + elseif(obs(1:9).eq."INTEG_XV3")then + FKObservables = xPDF(+2,x) - xPDF(-2,x) + 1 - ( xPDF(+1,x) - xPDF(-1,x) ) +* +**** Valence V8 +* + elseif(obs(1:9).eq."INTEG_XV8")then + FKObservables = xPDF(+2,x) - xPDF(-2,x) + 1 + ( xPDF(+1,x) - xPDF(-1,x) ) + 2 - 2d0 * ( xPDF(+3,x) - xPDF(-3,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 index 137de86f..7af300de 100644 --- 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:9).eq."INTEG_XT3")then + FKSimulator = ExternalEvolutionOperator("Ev2Ev",9,i,x,beta) +* +**** Non-singlet T8 +* + elseif(obs(1:9).eq."INTEG_XT8")then + FKSimulator = ExternalEvolutionOperator("Ev2Ev",10,i,x,beta) +* +**** Non-singlet T15 +* + elseif(obs(1:10).eq."INTEG_XT15")then + FKSimulator = ExternalEvolutionOperator("Ev2Ev",11,i,x,beta) +* +* +**** Valence V +* + elseif(obs(1:8).eq."INTEG_XV")then + FKSimulator = ExternalEvolutionOperator("Ev2Ev",3,i,x,beta) +* +**** Valence V3 +* + elseif(obs(1:9).eq."INTEG_XV3")then + FKSimulator = ExternalEvolutionOperator("Ev2Ev",4,i,x,beta) +* +**** Valence V8 +* + elseif(obs(1:9).eq."INTEG_XV8")then + FKSimulator = ExternalEvolutionOperator("Ev2Ev",5,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 index 59677003..c00bebe5 100644 --- a/src/DIS/SetFKObservable.f +++ b/src/DIS/SetFKObservable.f @@ -27,7 +27,8 @@ 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:5).eq."INTEG")then call SetTargetDIS("proton") * **** Light structure function F2light From 1989b083f11493f5058b206651365f9566e09786 Mon Sep 17 00:00:00 2001 From: Emanuele Nocera Date: Wed, 3 Jun 2020 12:38:05 +0200 Subject: [PATCH 2/5] Corrected small bugs --- src/DIS/FKObservables.f | 34 +++++++++++++++++----------------- src/DIS/FKSimulator.f | 12 ++++++------ 2 files changed, 23 insertions(+), 23 deletions(-) diff --git a/src/DIS/FKObservables.f b/src/DIS/FKObservables.f index 98aad769..23e0842e 100644 --- a/src/DIS/FKObservables.f +++ b/src/DIS/FKObservables.f @@ -95,47 +95,47 @@ function FKObservables(x,Q,y) * elseif(obs(1:9).eq."INTEG_XT3")then FKObservables = xPDF(+2,x) + xpdf(-2,x) - 1 - ( xPDF(+1,x) + xPDF(-1,x) ) + 1 - 1d0 * ( xPDF(+1,x) + xPDF(-1,x) ) * **** Non-singlet T8 * elseif(obs(1:9).eq."INTEG_XT8")then FKObservables = xPDF(+2,x) + xpdf(-2,x) - 1 + ( xPDF(+1,x) + xPDF(-1,x) ) + 1 + 1d0 * ( xPDF(+1,x) + xPDF(-1,x) ) 1 - 2d0 * ( xPDF(+3,x) + xPDF(-3,x) ) * **** Non-singlet T15 * elseif(obs(1:10).eq."INTEG_XT15")then FKObservables = xPDF(+2,x) + xpdf(-2,x) - 1 + ( xPDF(+1,x) + xPDF(-1,x) ) - 1 + ( xPDF(+3,x) + xPDF(-3,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 V -* - elseif(obs(1:8).eq."INTEG_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) -* **** Valence V3 * elseif(obs(1:9).eq."INTEG_XV3")then FKObservables = xPDF(+2,x) - xPDF(-2,x) - 1 - ( xPDF(+1,x) - xPDF(-1,x) ) + 1 - 1d0 * ( xPDF(+1,x) - xPDF(-1,x) ) * **** Valence V8 * elseif(obs(1:9).eq."INTEG_XV8")then FKObservables = xPDF(+2,x) - xPDF(-2,x) - 1 + ( xPDF(+1,x) - xPDF(-1,x) ) + 1 + 1d0 * ( xPDF(+1,x) - xPDF(-1,x) ) 2 - 2d0 * ( xPDF(+3,x) - xPDF(-3,x) ) * +* +**** Valence V +* + elseif(obs(1:8).eq."INTEG_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 index 7af300de..7a1b03c9 100644 --- a/src/DIS/FKSimulator.f +++ b/src/DIS/FKSimulator.f @@ -105,12 +105,6 @@ function FKSimulator(x,Q,y,i,beta) elseif(obs(1:10).eq."INTEG_XT15")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",11,i,x,beta) * -* -**** Valence V -* - elseif(obs(1:8).eq."INTEG_XV")then - FKSimulator = ExternalEvolutionOperator("Ev2Ev",3,i,x,beta) -* **** Valence V3 * elseif(obs(1:9).eq."INTEG_XV3")then @@ -121,6 +115,12 @@ function FKSimulator(x,Q,y,i,beta) elseif(obs(1:9).eq."INTEG_XV8")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",5,i,x,beta) * +* +**** Valence V +* + elseif(obs(1:8).eq."INTEG_XV")then + FKSimulator = ExternalEvolutionOperator("Ev2Ev",3,i,x,beta) +* **** Light structure function F2light * elseif(obs(1:7).eq."DIS_F2L")then From 8a4064e145c01f2f438330d64e0971f77475d408 Mon Sep 17 00:00:00 2001 From: Emanuele Nocera Date: Wed, 3 Jun 2020 13:49:44 +0200 Subject: [PATCH 3/5] Final tweaks for integrability observables --- src/DIS/FKObservables.f | 12 ++++++------ src/DIS/FKSimulator.f | 12 ++++++------ 2 files changed, 12 insertions(+), 12 deletions(-) diff --git a/src/DIS/FKObservables.f b/src/DIS/FKObservables.f index 23e0842e..50949038 100644 --- a/src/DIS/FKObservables.f +++ b/src/DIS/FKObservables.f @@ -93,20 +93,20 @@ function FKObservables(x,Q,y) * **** Non-singlet T3 * - elseif(obs(1:9).eq."INTEG_XT3")then + 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:9).eq."INTEG_XT8")then + 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:10).eq."INTEG_XT15")then + 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) ) @@ -114,13 +114,13 @@ function FKObservables(x,Q,y) * **** Valence V3 * - elseif(obs(1:9).eq."INTEG_XV3")then + 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:9).eq."INTEG_XV8")then + 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) ) @@ -128,7 +128,7 @@ function FKObservables(x,Q,y) * **** Valence V * - elseif(obs(1:8).eq."INTEG_XV")then + 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) diff --git a/src/DIS/FKSimulator.f b/src/DIS/FKSimulator.f index 7a1b03c9..fdb4d35d 100644 --- a/src/DIS/FKSimulator.f +++ b/src/DIS/FKSimulator.f @@ -92,33 +92,33 @@ function FKSimulator(x,Q,y,i,beta) * **** Non-singlet T3 * - elseif(obs(1:9).eq."INTEG_XT3")then + elseif(obs(1:7).eq."DIS_XT3")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",9,i,x,beta) * **** Non-singlet T8 * - elseif(obs(1:9).eq."INTEG_XT8")then + elseif(obs(1:7).eq."DIS_XT8")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",10,i,x,beta) * **** Non-singlet T15 * - elseif(obs(1:10).eq."INTEG_XT15")then + elseif(obs(1:8).eq."DIS_XT15")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",11,i,x,beta) * **** Valence V3 * - elseif(obs(1:9).eq."INTEG_XV3")then + elseif(obs(1:6).eq."DIS_XV3")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",4,i,x,beta) * **** Valence V8 * - elseif(obs(1:9).eq."INTEG_XV8")then + elseif(obs(1:6).eq."DIS_XV8")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",5,i,x,beta) * * **** Valence V * - elseif(obs(1:8).eq."INTEG_XV")then + elseif(obs(1:6).eq."DIS_XV")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",3,i,x,beta) * **** Light structure function F2light From b9b9943a6c201875b48193464f97f821abe2a486 Mon Sep 17 00:00:00 2001 From: Emanuele Nocera Date: Wed, 3 Jun 2020 14:45:20 +0200 Subject: [PATCH 4/5] Small tweak --- src/DIS/FKSimulator.f | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/DIS/FKSimulator.f b/src/DIS/FKSimulator.f index fdb4d35d..6bb8f10d 100644 --- a/src/DIS/FKSimulator.f +++ b/src/DIS/FKSimulator.f @@ -107,12 +107,12 @@ function FKSimulator(x,Q,y,i,beta) * **** Valence V3 * - elseif(obs(1:6).eq."DIS_XV3")then + elseif(obs(1:7).eq."DIS_XV3")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",4,i,x,beta) * **** Valence V8 * - elseif(obs(1:6).eq."DIS_XV8")then + elseif(obs(1:7).eq."DIS_XV8")then FKSimulator = ExternalEvolutionOperator("Ev2Ev",5,i,x,beta) * * From c1d751fb3a31c3257f287e07af8127be01fbb5d6 Mon Sep 17 00:00:00 2001 From: Emanuele Nocera Date: Wed, 3 Jun 2020 14:55:24 +0200 Subject: [PATCH 5/5] Ready to go --- src/DIS/SetFKObservable.f | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/src/DIS/SetFKObservable.f b/src/DIS/SetFKObservable.f index c00bebe5..9782ccee 100644 --- a/src/DIS/SetFKObservable.f +++ b/src/DIS/SetFKObservable.f @@ -28,7 +28,12 @@ subroutine SetFKObservable(obs) 5 obs(1:7).eq."DIS_XSB".or. 6 obs(1:7).eq."DIS_XCQ".or. 7 obs(1:7).eq."DIS_XGL".or. - 8 obs(1:5).eq."INTEG")then + 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