Reweight configuration option 3

[sarafiorendi]

Description of option 3 is available here.

Ntuples are available at
/eos/cms/store/user/fiorendi//p5prime/2018/skims/newphi/fixBkg/ntuples_with_pt_eta_weights
To train the topological-reweighting BDT (excluding bPt, bEta, and similar) the MC events should be weighted by the variables:
weight*weight_reco_pteta
I also provide the data sample from the sPlot used to get the pt/eta weights
out_distribution_JPSI_2018A_2018D_L1all_preBDT.root and the friend tree containing all the variables
data_charmonium_2018_preBDT_fixBool_asSplot.root

It was not possible to save all the vars in the output-from-splot since there are too many events being at the preselection level.

[sarafiorendi]

Make sure to apply these cuts before using these ntuples:

pass_preselection ==1 (both for data and MC)
trig==1 for MC only

For reference, to get the weights, I applied the selection below

selData = '( ( (bMass*tagB0 + (1-tagB0)*bBarMass)   > {M}-3*{S}   && \
               (bMass*tagB0 + (1-tagB0)*bBarMass)   < {M}+3*{S} ) && \
              ( pass_preselection ==1 ) && \
              (abs(mumuMass - {JPSIM}) < {CUT}*mumuMassE))'\
          .format(M=mc_mass,S=mc_sigma,  JPSIM=JPsiMass_, CUT=nSigma_psiRej)
selMC = selData + ' && (trig==1)'

having defined

mc_sigma = 0.040
mc_mass  = 5.27783 
JPsiMass_ = 3.096916
nSigma_psiRej = 3.

[sarafiorendi]

Just to add, the weight to be used to get the sPlot distribution from data is nsig_sw.
No need to refit the data.

[XuelongQin]

I have summarized the plots of v1 and this option without pteta weights in
https://indico.cern.ch/event/1150965/contributions/4830152/attachments/2427712/4156545/Reweighting_summary.pdf

[sarafiorendi]

thanks! the corresponding weights are these?

[XuelongQin]

This is the reweighting model, you need to use Application.py to generate weight for each MC sample.

[sarafiorendi]

ok, can you save instead directly the weights?
otherwise each time we should re-open the jpsi data/MC and recompute them.
e.g. as done here

[XuelongQin]

Okay, I thought you mean the weights of other q2Bins. i will give the link of root file with weights of Jpsi later.