Add nu Q2 x plots for simulation

Add nu Q2 x plots for simulation Will re-write jpsi kinematics in a simplified way (like nu Q2 x) Will re-write the plotting part to factorize out code

dvmp/analysis/vm_invar.cxx

@@ -106,10 +106,10 @@ int vm_invar(const std::string& config_name) {
     
     // Print canvas to output file
     
-    TCanvas c2{"canvas2", "canvas2", 1800, 600};
-    c2.Divide(3, 1, 0.0001, 0.0001);
+    TCanvas c{"canvas2", "canvas2", 1800, 600};
+    c.Divide(3, 1, 0.0001, 0.0001);
     //pad 1 nu
-    c2.cd(1);
+    c.cd(1);
     //gPad->SetLogx(false);
     //gPad->SetLogy(false);
     auto& hnu = *h_nu_sim;
@@ -135,7 +135,7 @@ int vm_invar(const std::string& config_name) {
     t21->Draw();
     
     //pad 2 Q2
-    c2.cd(2);
+    c.cd(2);
     //gPad->SetLogx(false);
     //gPad->SetLogy(false);
     auto& hQ2 = *h_Q2_sim;
@@ -160,7 +160,7 @@ int vm_invar(const std::string& config_name) {
     t22->Draw();
     
     //pad 1 nu
-    c2.cd(3);
+    c.cd(3);
     //gPad->SetLogx(false);
     //gPad->SetLogy(false);
     auto& hx = *h_x_sim;
@@ -184,7 +184,7 @@ int vm_invar(const std::string& config_name) {
     //tptr1->SetTextColor(plot::kMpOrange);
     t23->Draw();
     
-    c2.Print(fmt::format("{}InvariantQuantities.png", output_prefix).c_str());
+    c.Print(fmt::format("{}InvariantQuantities.png", output_prefix).c_str());
   }
 
   // TODO we're not actually doing an IM fit yet, so for now just return an