void writevectordouble(vector<double> vec, TTree* tree, int entry, std::string leaf) {

  tree->GetEntry(entry);
  if (entry == 0) tree->Branch(leaf.c_str(), &vec);
  tree->Fill();

}

void writevectorint(vector<int> vec, TTree* tree, int entry, std::string leaf) {

  tree->GetEntry(entry);
  if (entry == 0) tree->Branch(leaf.c_str(), &vec);
  tree->Fill();

}

void writeint(int integer, TTree* tree, int entry, std::string leaf) {

  tree->GetEntry(entry);
  if (entry == 0) tree->Branch(leaf.c_str(), &integer);
  tree->Fill();

}