[Rivet] Is it possible? Rivet and truth look different

Andy Buckley andy.buckley at cern.ch
Thu Nov 7 17:16:59 GMT 2013


I guess so, but am not sure. And the default ATLAS truth jets probably
apply some restrictions on the clustering inputs such as excluding muons
and neutrinos. I'd be surprised if that made a big difference, though.

I'm also not sure if ATLAS eta/y acceptance cuts are used in the DPD
default output, but that would not be unreasonable. In which case Rivet
would find more jets.

However, I just realised that the curves labelled Sherpa1 and Sherpa2
are not cryptically labelled ATLAS truth jets vs. Rivet, but really
Sherpa1 and Sherpa2... right? So we can't actually see your Rivet output
because you just sent the index.html file, no plots!

Andy


On 07/11/13 13:52, Liron Barak wrote:
> 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
> <mailto: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
> 
> 
> 
> 
> _______________________________________________
> Rivet mailing list
> Rivet at projects.hepforge.org
> http://www.hepforge.org/lists/listinfo/rivet
> 


-- 
Dr Andy Buckley, Royal Society University Research Fellow
Particle Physics Expt Group, University of Glasgow / PH Dept, CERN


More information about the Rivet mailing list