{ #include #include #include #include "TString.h" // Header file for the classes stored in the TTree if any. #include "TLorentzVector.h" #include "vector" #include "vector" #include "vector" #include #include #include TTree::SetMaxTreeSize(200000000000LL); TFile *tfile = new TFile("hist-tree_mce.361293.root", "RECREATE"); TTree* tree1; TFile *f = TFile::Open("hist-mc16_13TeV.361293.MGaMcAtNloPy8EG_NNPDF30LO_A14NNPDF23LO_WZ_lvll_FxFx.deriv.DAOD_STDM5.e7958_s3126_r10724_p4097.root"); TDirectory * dir = (TDirectory*)f->Get("hist-mc16_13TeV.361293.MGaMcAtNloPy8EG_NNPDF30LO_A14NNPDF23LO_WZ_lvll_FxFx.deriv.DAOD_STDM5.e7958_s3126_r10724_p4097.root:/Nominal"); dir->GetObject("TruthInclusiveRegions_tree_DBGN_OBJMINDELTAR(kTruthWLepton,kTruthZLepton2,0.3)",tree1); tree1->GetBranch("genWeightdyn=-1muR=010000E+01muF=010000E+01")->SetTitle("genWeightdyn1muR_010000E01muF_010000E01"); tree1->GetBranch("genWeightdyn=-1muR=010000E+01muF=010000E+01")->SetName("genWeightdyn1muR_010000E01muF_010000E01"); tree1->GetLeaf("genWeightdyn=-1muR=010000E+01muF=010000E+01")->SetTitle("genWeightdyn1muR_010000E01muF_010000E01"); tree1->GetLeaf("genWeightdyn=-1muR=010000E+01muF=010000E+01")->SetName("genWeightdyn1muR_010000E01muF_010000E01"); tree1->GetBranch("genWeightdyn=-1muR=020000E+01muF=010000E+01")->SetTitle("genWeightdyn1muR_020000E01muF_010000E01"); tree1->GetBranch("genWeightdyn=-1muR=020000E+01muF=010000E+01")->SetName("genWeightdyn1muR_020000E01muF_010000E01"); tree1->GetLeaf("genWeightdyn=-1muR=020000E+01muF=010000E+01")->SetTitle("genWeightdyn1muR_020000E01muF_010000E01"); tree1->GetLeaf("genWeightdyn=-1muR=020000E+01muF=010000E+01")->SetName("genWeightdyn1muR_020000E01muF_010000E01"); tree1->GetBranch("genWeightdyn=-1muR=050000E+00muF=010000E+01")->SetTitle("genWeightdyn1muR_050000E00muF_010000E01"); tree1->GetBranch("genWeightdyn=-1muR=050000E+00muF=010000E+01")->SetName("genWeightdyn1muR_050000E00muF_010000E01"); tree1->GetLeaf("genWeightdyn=-1muR=050000E+00muF=010000E+01")->SetTitle("genWeightdyn1muR_050000E00muF_010000E01"); tree1->GetLeaf("genWeightdyn=-1muR=050000E+00muF=010000E+01")->SetName("genWeightdyn1muR_050000E00muF_010000E01"); tree1->GetBranch("genWeightdyn=-1muR=010000E+01muF=020000E+01")->SetTitle("genWeightdyn1muR_010000E01muF_020000E01"); tree1->GetBranch("genWeightdyn=-1muR=010000E+01muF=020000E+01")->SetName("genWeightdyn1muR_010000E01muF_020000E01"); tree1->GetLeaf("genWeightdyn=-1muR=010000E+01muF=020000E+01")->SetTitle("genWeightdyn1muR_010000E01muF_020000E01"); tree1->GetLeaf("genWeightdyn=-1muR=010000E+01muF=020000E+01")->SetName("genWeightdyn1muR_010000E01muF_020000E01"); tree1->GetBranch("genWeightdyn=-1muR=020000E+01muF=020000E+01")->SetTitle("genWeightdyn1muR_020000E01muF_020000E01"); tree1->GetBranch("genWeightdyn=-1muR=020000E+01muF=020000E+01")->SetName("genWeightdyn1muR_020000E01muF_020000E01"); tree1->GetLeaf("genWeightdyn=-1muR=020000E+01muF=020000E+01")->SetTitle("genWeightdyn1muR_020000E01muF_020000E01"); tree1->GetLeaf("genWeightdyn=-1muR=020000E+01muF=020000E+01")->SetName("genWeightdyn1muR_020000E01muF_020000E01"); tree1->GetBranch("genWeightdyn=-1muR=050000E+00muF=020000E+01")->SetTitle("genWeightdyn1muR_050000E00muF_020000E01"); tree1->GetBranch("genWeightdyn=-1muR=050000E+00muF=020000E+01")->SetName("genWeightdyn1muR_050000E00muF_020000E01"); tree1->GetLeaf("genWeightdyn=-1muR=050000E+00muF=020000E+01")->SetTitle("genWeightdyn1muR_050000E00muF_020000E01"); tree1->GetLeaf("genWeightdyn=-1muR=050000E+00muF=020000E+01")->SetName("genWeightdyn1muR_050000E00muF_020000E01"); tree1->GetBranch("genWeightdyn=-1muR=010000E+01muF=050000E+00")->SetTitle("genWeightdyn1muR_010000E01muF_050000E00"); tree1->GetBranch("genWeightdyn=-1muR=010000E+01muF=050000E+00")->SetName("genWeightdyn1muR_010000E01muF_050000E00"); tree1->GetLeaf("genWeightdyn=-1muR=010000E+01muF=050000E+00")->SetTitle("genWeightdyn1muR_010000E01muF_050000E00"); tree1->GetLeaf("genWeightdyn=-1muR=010000E+01muF=050000E+00")->SetName("genWeightdyn1muR_010000E01muF_050000E00"); tree1->GetBranch("genWeightdyn=-1muR=020000E+01muF=050000E+00")->SetTitle("genWeightdyn1muR_020000E01muF_050000E00"); tree1->GetBranch("genWeightdyn=-1muR=020000E+01muF=050000E+00")->SetName("genWeightdyn1muR_020000E01muF_050000E00"); tree1->GetLeaf("genWeightdyn=-1muR=020000E+01muF=050000E+00")->SetTitle("genWeightdyn1muR_020000E01muF_050000E00"); tree1->GetLeaf("genWeightdyn=-1muR=020000E+01muF=050000E+00")->SetName("genWeightdyn1muR_020000E01muF_050000E00"); tree1->GetBranch("genWeightdyn=-1muR=050000E+00muF=050000E+00")->SetTitle("genWeightdyn1muR_050000E00muF_050000E00"); tree1->GetBranch("genWeightdyn=-1muR=050000E+00muF=050000E+00")->SetName("genWeightdyn1muR_050000E00muF_050000E00"); tree1->GetLeaf("genWeightdyn=-1muR=050000E+00muF=050000E+00")->SetTitle("genWeightdyn1muR_050000E00muF_050000E00"); tree1->GetLeaf("genWeightdyn=-1muR=050000E+00muF=050000E+00")->SetName("genWeightdyn1muR_050000E00muF_050000E00"); string name[] = { "genWeightPDF=260000NNPDF30_nlo_as_0118","genWeightPDF=260001NNPDF30_nlo_as_0118", "genWeightPDF=260002NNPDF30_nlo_as_0118","genWeightPDF=260003NNPDF30_nlo_as_0118", "genWeightPDF=260004NNPDF30_nlo_as_0118","genWeightPDF=260005NNPDF30_nlo_as_0118", "genWeightPDF=260006NNPDF30_nlo_as_0118","genWeightPDF=260007NNPDF30_nlo_as_0118", "genWeightPDF=260008NNPDF30_nlo_as_0118","genWeightPDF=260009NNPDF30_nlo_as_0118", "genWeightPDF=260010NNPDF30_nlo_as_0118","genWeightPDF=260011NNPDF30_nlo_as_0118", "genWeightPDF=260012NNPDF30_nlo_as_0118","genWeightPDF=260013NNPDF30_nlo_as_0118", "genWeightPDF=260014NNPDF30_nlo_as_0118","genWeightPDF=260015NNPDF30_nlo_as_0118", "genWeightPDF=260016NNPDF30_nlo_as_0118","genWeightPDF=260017NNPDF30_nlo_as_0118", "genWeightPDF=260018NNPDF30_nlo_as_0118","genWeightPDF=260019NNPDF30_nlo_as_0118", "genWeightPDF=260020NNPDF30_nlo_as_0118","genWeightPDF=260021NNPDF30_nlo_as_0118", "genWeightPDF=260022NNPDF30_nlo_as_0118","genWeightPDF=260023NNPDF30_nlo_as_0118", "genWeightPDF=260024NNPDF30_nlo_as_0118","genWeightPDF=260025NNPDF30_nlo_as_0118", "genWeightPDF=260026NNPDF30_nlo_as_0118","genWeightPDF=260027NNPDF30_nlo_as_0118", "genWeightPDF=260028NNPDF30_nlo_as_0118","genWeightPDF=260029NNPDF30_nlo_as_0118", "genWeightPDF=260030NNPDF30_nlo_as_0118","genWeightPDF=260031NNPDF30_nlo_as_0118", "genWeightPDF=260032NNPDF30_nlo_as_0118","genWeightPDF=260033NNPDF30_nlo_as_0118", "genWeightPDF=260034NNPDF30_nlo_as_0118","genWeightPDF=260035NNPDF30_nlo_as_0118", "genWeightPDF=260036NNPDF30_nlo_as_0118","genWeightPDF=260037NNPDF30_nlo_as_0118", "genWeightPDF=260038NNPDF30_nlo_as_0118","genWeightPDF=260039NNPDF30_nlo_as_0118", "genWeightPDF=260040NNPDF30_nlo_as_0118","genWeightPDF=260041NNPDF30_nlo_as_0118", "genWeightPDF=260042NNPDF30_nlo_as_0118","genWeightPDF=260043NNPDF30_nlo_as_0118", "genWeightPDF=260044NNPDF30_nlo_as_0118","genWeightPDF=260045NNPDF30_nlo_as_0118", "genWeightPDF=260046NNPDF30_nlo_as_0118","genWeightPDF=260047NNPDF30_nlo_as_0118", "genWeightPDF=260048NNPDF30_nlo_as_0118","genWeightPDF=260049NNPDF30_nlo_as_0118", "genWeightPDF=260050NNPDF30_nlo_as_0118","genWeightPDF=260051NNPDF30_nlo_as_0118", "genWeightPDF=260052NNPDF30_nlo_as_0118","genWeightPDF=260053NNPDF30_nlo_as_0118", "genWeightPDF=260054NNPDF30_nlo_as_0118","genWeightPDF=260055NNPDF30_nlo_as_0118", "genWeightPDF=260056NNPDF30_nlo_as_0118","genWeightPDF=260057NNPDF30_nlo_as_0118", "genWeightPDF=260058NNPDF30_nlo_as_0118","genWeightPDF=260059NNPDF30_nlo_as_0118", "genWeightPDF=260060NNPDF30_nlo_as_0118","genWeightPDF=260061NNPDF30_nlo_as_0118", "genWeightPDF=260062NNPDF30_nlo_as_0118","genWeightPDF=260063NNPDF30_nlo_as_0118", "genWeightPDF=260064NNPDF30_nlo_as_0118","genWeightPDF=260065NNPDF30_nlo_as_0118", "genWeightPDF=260066NNPDF30_nlo_as_0118","genWeightPDF=260067NNPDF30_nlo_as_0118", "genWeightPDF=260068NNPDF30_nlo_as_0118","genWeightPDF=260069NNPDF30_nlo_as_0118", "genWeightPDF=260070NNPDF30_nlo_as_0118","genWeightPDF=260071NNPDF30_nlo_as_0118", "genWeightPDF=260072NNPDF30_nlo_as_0118","genWeightPDF=260073NNPDF30_nlo_as_0118", "genWeightPDF=260074NNPDF30_nlo_as_0118","genWeightPDF=260075NNPDF30_nlo_as_0118", "genWeightPDF=260076NNPDF30_nlo_as_0118","genWeightPDF=260077NNPDF30_nlo_as_0118", "genWeightPDF=260078NNPDF30_nlo_as_0118","genWeightPDF=260079NNPDF30_nlo_as_0118", "genWeightPDF=260080NNPDF30_nlo_as_0118","genWeightPDF=260081NNPDF30_nlo_as_0118", "genWeightPDF=260082NNPDF30_nlo_as_0118","genWeightPDF=260083NNPDF30_nlo_as_0118", "genWeightPDF=260084NNPDF30_nlo_as_0118","genWeightPDF=260085NNPDF30_nlo_as_0118", "genWeightPDF=260086NNPDF30_nlo_as_0118","genWeightPDF=260087NNPDF30_nlo_as_0118", "genWeightPDF=260088NNPDF30_nlo_as_0118","genWeightPDF=260089NNPDF30_nlo_as_0118", "genWeightPDF=260090NNPDF30_nlo_as_0118","genWeightPDF=260091NNPDF30_nlo_as_0118", "genWeightPDF=260092NNPDF30_nlo_as_0118","genWeightPDF=260093NNPDF30_nlo_as_0118", "genWeightPDF=260094NNPDF30_nlo_as_0118","genWeightPDF=260095NNPDF30_nlo_as_0118", "genWeightPDF=260096NNPDF30_nlo_as_0118","genWeightPDF=260097NNPDF30_nlo_as_0118", "genWeightPDF=260098NNPDF30_nlo_as_0118","genWeightPDF=260099NNPDF30_nlo_as_0118", "genWeightPDF=260100NNPDF30_nlo_as_0118"}; int nNames = sizeof(name)/sizeof(name[0]); for (int in =0 ; in< nNames; in++) { char buffer [33]; string S= std::to_string(in); std::string genName="genPDF_"+S; const char *nname=name[in].c_str(); const char *title=genName.c_str(); tree1->GetBranch(nname)->SetTitle(title); tree1->GetBranch(nname)->SetName(title); tree1->GetLeaf(nname)->SetTitle(title); tree1->GetLeaf(nname)->SetTitle(title); } tfile->cd(); tree1->CloneTree()->Write(); tree1->Print(); delete(f); delete(tfile); }