|
[Rivet] Is it possible? Rivet and truth look differentLiron Barak lironbarak83 at gmail.comThu Nov 7 12:52:10 GMT 2013
Hi Andy, The "ROOT" analysis.... is only making the truthD3PDs from evgens with: Reco_trf.py preExec='from D3PDMakerConfig.D3PDMakerFlags import D3PDMakerFlags;D3PDMakerFlags.TruthWriteHadrons.set_Value_and_Lock(True)' inputEVNTFile=1111.pool.root outputNTUP_TRUTHFile=11111.truth.1.root autoConfiguration=everything Then I simply do stupid draw command: sherpa1.Draw('jet_AntiKt4TruthJets_n','','hist&&norm') sherpa2.Draw('jet_AntiKt4TruthJets_n','','hist&&norm&&same') So I didn't apply any isolation here..... I would guess it is also all the jets, no? Thanks Liron On Thu, Nov 7, 2013 at 1:05 PM, Andy Buckley <andy.buckley at cern.ch> wrote: > On 07/11/13 09:48, Liron Barak wrote: > > Hi Andy, > > > > Now that Rivet is working for me, thanks! > > I see something weird.... > > I took the official sherpa samples (made by sherpa 1) and compared to my > > private sherpa2 samples... I used MC_TTBAR in river to compare (with > > evgens) and simple Draw of the truth jet branch with the truthD3PDs I > > made... > > Please see the attached.... the Jet_N for both samples behave opposite > > in the 2 plots, no? > > Hi Liron, > > Can you check if there is lepton/jet isolation in the ROOT analysis? In > Rivet we make jets from *everything* provided to the FastJets, so hard > leptons will count as jets unless they are identified and removed. > > This is something which would be nice to improve, as is the definition > of jet inputs, but unlike the ATLAS truth jets we need to be generic > rather than expt-specific. I hope we'll provide an AtlasJets projection > to make this easier at some point, but even ATLAS itself doesn't have a > single truth-jet constituent definition yet! > > Andy > > -- > Dr Andy Buckley, Royal Society University Research Fellow > Particle Physics Expt Group, University of Glasgow / PH Dept, CERN > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://www.hepforge.org/lists-archive/rivet/attachments/20131107/853912ba/attachment.html>
More information about the Rivet mailing list |