library "FebrileSeizure" version '0.0.0'
using FHIR version '4.0.0'
include "FHIRHelpers" version '4.0.0' called FHIRHelpers
codesystem "ICD-10-CM": 'http://hl7.org/fhir/sid/icd-10-cm'
codesystem "CVX": 'http://hl7.org/fhir/sid/cvx'
codesystem "RxNorm": 'http://www.nlm.nih.gov/research/umls/rxnorm'
codesystem "NDC": 'http://hl7.org/fhir/sid/ndc'
codesystem "CPT": 'http://www.ama-assn.org/go/cpt'
codesystem "SNOMED": 'http://snomed.info/sct'

define "Covid Vaccine Drug Exposure": Concept {
Code '207' from "CVX",
Code '208' from "CVX",
Code '210' from "CVX",
Code '212' from "CVX",
Code '213' from "CVX",
Code '213' from "CVX",
Code '2470234' from "RxNorm",
Code '2468235' from "RxNorm",
Code '59267100003' from "NDC",
Code '59267100002' from "NDC",
Code '80777027310' from "NDC",
Code '80777027399' from "NDC",
Code '59676058005' from "NDC",
Code '00310122210' from "NDC",
Code '59676058015' from "NDC",
Code '59267100001' from "NDC"
}

define "Covid Vaccine Procedure": Concept {
Code '0002A' from "CPT",
Code '0001A' from "CPT",
Code '91300' from "CPT",
Code '0012A' from "CPT",
Code '0011A' from "CPT",
Code '90301' from "CPT",
Code '0022A' from "CPT",
Code '0021A' from "CPT",
Code '91302' from "CPT"
}

define "FebrileSeizure": Concept{
  Code 'R56.00' from "ICD-10-CM",
  Code 'R56.01' from "ICD-10-CM",
  Code 'G40.89' from "ICD-10-CM",
  Code 'R50.81' from "ICD-10-CM",
  Code 'R50.82' from "ICD-10-CM",
  Code 'R50.83' from "ICD-10-CM",
  Code 'R50.84' from "ICD-10-CM",
  Code 'R50.9' from "ICD-10-CM",
  Code 'R56.9' from "ICD-10-CM",
  Code '432354000' from "SNOMED",
  Code '433083002' from "SNOMED",
  Code '91175000' from "SNOMED",
  Code '386661006' from "SNOMED",
  Code '130091000119103' from "SNOMED",
  Code '123471000119103' from "SNOMED",
  Code '73301000' from "SNOMED",
  Code '386661006' from "SNOMED",
  Code '91175000' from "SNOMED"
  }

context Patient

define "FebrileSeizureCondition":
  ConditionLookBack([Condition: "FebrileSeizure"], 365 days)

define "CovidImmunization": [Immunization: "Covid Vaccine Drug Exposure"]

define "CovidProcedure": [Procedure: "Covid Vaccine Procedure"]

define "MostRecentCovidI": MostRecentImmunizationDT("CovidImmunization")

define "MostRecentCovidP": MostRecentProcedureDT("CovidProcedure")

define "LastCovidInstanceDT":
  case
        when not exists("MostRecentCovidI") then "MostRecentCovidP"
        when not exists("MostRecentCovidP") then "MostRecentCovidI"
        when "MostRecentCovidI" after "MostRecentCovidP" then "MostRecentCovidI"
        else "MostRecentCovidP"
  end

define "RiskPeriod":
  Interval[ Last("LastCovidInstanceDT"), Last("LastCovidInstanceDT") + 1 day]

define "CleanPeriod":
  Interval[ Last("LastCovidInstanceDT") - 42 days, Last("LastCovidInstanceDT")]

define "FebrileSeizureinRiskPeriod": "FebrileSeizureCondition" cond
  where (cond.onset as FHIR.dateTime).value in "RiskPeriod"
  or PeriodToInterval(cond.onset as FHIR.Period) overlaps "RiskPeriod"

define "FebrileSeizureinCleanPeriod": "FebrileSeizureCondition" cond
    where (cond.onset as FHIR.dateTime).value in "CleanPeriod"
    or PeriodToInterval(cond.onset as FHIR.Period) overlaps "CleanPeriod"

define "InPopulation":
  exists(FebrileSeizureinRiskPeriod) and not exists(FebrileSeizureinCleanPeriod)

define function MostRecentImmunizationDT(ImmunizationList List<FHIR.Immunization>):
  from ImmunizationList I
  let result: Last(ImmunizationList I sort by (occurrence as FHIR.dateTime).value)
    return (result.occurrence as FHIR.dateTime).value

define function MostRecentProcedureDT(ProcedureList List<FHIR.Procedure>):
  from ProcedureList P
  let result:Last(ProcedureList P sort by (performed as FHIR.dateTime).value)
    return (result.performed as FHIR.dateTime).value

define function ConditionLookBack(CondList List<Condition>, LookBack System.Quantity):
  CondList C
    let LookBackInterval: Interval[Now() - LookBack, Now()]
    where (C.onset as FHIR.dateTime).value in LookBackInterval
      or PeriodToInterval(C.onset as FHIR.Period) overlaps LookBackInterval
      and Last(C.verificationStatus.coding.code.value) = 'confirmed'

define function PeriodToInterval(period FHIR.Period):
  if period is null then
    null
  else
    Interval[period."start".value, period."end".value]