void AddTiStates(double E);

void Kine(){
  
  
 NPL::Reaction K("47K(d,p)48K@362");
  
 // K states
 auto g = K.GetKinematicLine3();
 g->SetLineColor(kAzure+7);
 g->SetLineWidth(3);
 g->Draw("ac");
 K.SetExcitationHeavy(4);
 g = K.GetKinematicLine3();
 g->SetLineColor(kAzure+7);
 g->SetLineWidth(2);
 g->SetLineStyle(1);
 g->Draw("c");

 AddTiStates(0); 

 AddTiStates(0); 
 AddTiStates(0.969); 
 AddTiStates(2.2292); 
 AddTiStates(2.419); 
 AddTiStates(3.223); 
 AddTiStates(3.332); 
 AddTiStates(3.622); 
 AddTiStates(4.388); 
 AddTiStates(4.458); 
 AddTiStates(4.719); 
 AddTiStates(4.852); 
 AddTiStates(5.151); 
   
  }

void AddTiStates(double E){
  
 NPL::Reaction Ti("47Ti(d,p)48Ti@362");
   // Ti states
 Ti.SetExcitationHeavy(E);
 auto g = Ti.GetKinematicLine3();
 g->SetLineWidth(1);
 g->SetLineStyle(2);
 g->Draw("c");


  
  }