Fork me on GitHub

Changeset d244bc9 in git


Ignore:
Timestamp:
Dec 4, 2014, 3:25:06 PM (10 years ago)
Author:
Michele <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
51f7063
Parents:
6fb1a5d
Message:

added exception in FastJetFinder when winner-take-all recombination scheme is used that prevents jets from having mass = 0

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/FastJetFinder.cc

    r6fb1a5d rd244bc9  
    223223  Candidate *candidate, *constituent;
    224224  TLorentzVector momentum;
     225 
     226  TLorentzVector constmomentum;
     227 
    225228  Double_t deta, dphi, detaMax, dphiMax;
    226229  Double_t time, weightTime, avTime;
     
    286289  for(itOutputList = outputList.begin(); itOutputList != outputList.end(); ++itOutputList)
    287290  {
    288     momentum.SetPxPyPzE(itOutputList->px(), itOutputList->py(), itOutputList->pz(), itOutputList->E());
     291    jet = *itOutputList;
     292    if(fJetAlgorithm == 7) jet = join(jet.constituents());
     293   
     294    momentum.SetPxPyPzE(jet.px(), jet.py(), jet.pz(), jet.E());
     295   
    289296    area.reset(0.0, 0.0, 0.0, 0.0);
    290297    if(fAreaDefinition) area = itOutputList->area_4vector();
     
    297304    inputList.clear();
    298305    inputList = sequence->constituents(*itOutputList);
     306
     307    constmomentum.SetPxPyPzE(0.0,0.0,0.0,0.0);;
    299308
    300309    for(itInputList = inputList.begin(); itInputList != inputList.end(); ++itInputList)
    301310    {
    302311      constituent = static_cast<Candidate*>(fInputArray->At(itInputList->user_index()));
     312
     313      constmomentum += constituent->Momentum;
    303314
    304315      deta = TMath::Abs(momentum.Eta() - constituent->Momentum.Eta());
Note: See TracChangeset for help on using the changeset viewer.