Save Tree with weight

[wamorkar.t]

Dear Experts,

I have a tree with multiple branches, two of which are “tp_mass” and “genweight”.
If I have to draw “tp_mass” with “genweight” I do:

Tree->Draw("tp_mass","genweight")

After having done this, is it possible to save the tree such that “tp_mass” now has the “genweight” weight applied to it?

Thanks in advance!
Tanvi

[Danilo]

Hi Tanvi,

you can solve this issue with TDataFrame.
In your case:

// assuming these names :)
ROOT::Experimental::TDataFrame f("TreeName", "FileName.root"); 
f.Define("weightedMass", "tp_mass * genweight").Snapshot("myNewTree", "myNewFileName");

If you also want the histogram:

ROOT::Experimental::TDataFrame f("TreeName", "FileName.root"); 
auto h = f.Histo1D(“tp_mass”,“genweight”);
f.Define("weightedMass", "tp_mass * genweight")..Snapshot("myNewTree", "myNewFileName");
h.Draw();

Cheers,
D

[wamorkar.t]

Thanks Danilo! This was extremely helpful.

[wamorkar.t]

Hi Danilo,

Sorry I have another question. I created the reweighted tree. However there is a difference between the two distributions.
In detail:

Case 1 : Draw a branch tp_mass with weight “genTotalweight”
Tree->Draw(“tp_mass”,“genTotalweight”)
gives me Screen Shot 2017-10-11 at 11.28.48 AM

Case 2 : From the reweighted tree I created (using the method you suggested)
On drawing the reweighted branch I get
Screen Shot 2017-10-11 at 11.29.55 AM

Its clear that the two distributions are not the same , whereas I think they should be.
Is there a difference between the two methods of drawing with weights?

Thanks,
Tanvi

[Danilo]

Hi Tanvi,

I misinterpreted your first request: my bad. You cannot save the value and its weight in a single value by multiplying those
The plot is telling us that

Cheers,
D

[Axel]

Right, to give some background: you have values for tp_mass, say 5, 10, 7, 100. These values are stored in the branch tp_mass. Each of those entries has a different weight: 0.1, 0.2, 100, 0.01. No multiplication or addition or whatsoever can combine these, they are completely independent numbers. You will need to keep two branches, and weight each tp_mass entry by the genweight value - as you did with TTree::Draw().

[wamorkar.t]

Yes Thanks for the tip!