# Inserting data

pappsodata <- data.frame(annee=rep(2001:2022,2),
                     
                     dd=c(0.1, 0.3, 0.4, 0.7),
                     
                     pappso=c(0.1, 0.1, 0.1, 0.1)
                     )

espace_dd <- c(c(0.1, 0.3, 0.4, 0.7),c(0.8,0.9,1,1.2,1.8,2), c(2.5,3,3.5,4,4.5,5), c(7,8,9,10,12,14))

espace_pappso <- c(c(0.1, 0.1, 0.1, 0.1),c(0.3,0.6,1,1.5,2,2.3), c(2.6,3.5,5,6.5,8,9), c(12,15,20,30,40,50))

pappsodata$type <- c(rep("disques durs", 22), rep("PAPPSO", 22))

pappsodata$espace <- c(espace_dd,espace_pappso)

# Bar Plot
library(ggplot2)


plt <-ggplot(data=pappsodata[which(pappsodata$type == "PAPPSO"),], aes(x=annee, y=espace, colour = type, shape = type)) +
  
  geom_line() +
  geom_point() +
  labs(x = "Années", y = "To",color = "Stockage", shape = "Stockage") +
  geom_label(aes(x = 2005, y = 5, label = "LTQ"), 
             hjust = +1.3, 
             vjust = 0.2, 
             colour = "#FAAB18", 
             label.size = NA, 
             family="Helvetica", 
             size = 6) +
  geom_curve(aes(x = 2004, y = 5, xend = 2005, yend = 0.31), 
             colour = "#FAAB18", 
             size=0.5, 
             curvature = -0.2, 
             arrow = arrow(length = unit(0.03, "npc"))) +
  geom_label(aes(x = 2007, y = 7, label = "Orbitrap"), 
             hjust = +0.5, 
             vjust = 0.2, 
             colour = "#FAAB18", 
             label.size = NA, 
             family="Helvetica", 
             size = 6) +
  geom_curve(aes(x = 2007, y = 7, xend = 2007, yend = 1), 
             colour = "#FAAB18", 
             size=0.5, 
             curvature = -0.2, 
             arrow = arrow(length = unit(0.03, "npc"))) +
  geom_label(aes(x = 2008, y = 10, label = "LTQ-XL"), 
             hjust = +0.5, 
             vjust = 0.2, 
             colour = "#FAAB18", 
             label.size = NA, 
             family="Helvetica", 
             size = 6) +
  geom_curve(aes(x = 2008, y = 10, xend = 2008, yend = 1.5), 
             colour = "#FAAB18", 
             size=0.5, 
             curvature = -0.2, 
             arrow = arrow(length = unit(0.03, "npc"))) +
  geom_label(aes(x = 2011, y = 13, label = "Q-Exactive"), 
             hjust = +0.5, 
             vjust = 0.2, 
             colour = "#FAAB18", 
             label.size = NA, 
             family="Helvetica", 
             size = 6) +
  geom_curve(aes(x = 2011, y = 13, xend = 2011, yend = 2.4), 
             colour = "#FAAB18", 
             size=0.5, 
             curvature = -0.2, 
             arrow = arrow(length = unit(0.03, "npc"))) +
  geom_label(aes(x = 2015, y = 18, label = "Lumos"), 
             hjust = +0.5, 
             vjust = -0.1, 
             colour = "#FAAB18", 
             label.size = NA, 
             family="Helvetica", 
             size = 6) +
  geom_curve(aes(x = 2015, y = 18, xend = 2015, yend = 8), 
             colour = "#FAAB18", 
             size=0.5, 
             curvature = -0.2, 
             arrow = arrow(length = unit(0.03, "npc"))) +
  geom_label(aes(x = 2020, y = 30, label = "timsTOF pro"), 
             hjust = +1.2, 
             vjust = -0.1, 
             colour = "#FAAB18", 
             label.size = NA, 
             family="Helvetica", 
             size = 6) +
  geom_curve(aes(x = 2018, y = 30, xend = 2019, yend = 20), 
             colour = "#FAAB18", 
             size=0.5, 
             curvature = -0.2, 
             arrow = arrow(length = unit(0.03, "npc"))) +
  geom_label(aes(x = 2008, y = 45, label = "Instruments"), 
             hjust = +1.2, 
             vjust = -0.1, 
             colour = "#FAAB18", 
             label.size = NA, 
             family="Helvetica", 
             size = 6) +
  theme_classic()

plt

ggsave("deluge.png", width = 8, height = 5)