/** \class BlockClasses
 *
 *  Initialization of static fgCompare functions. At the moment a number of 
 *  functions TCompare*** are implemented (see BlockCompare.h). 
 *  Function TCompareXX sorts objects 
 *  by the variable "XX" that MUST be present in the data members of the 
 *  Block TRoot class.
 *  By default most objects are sorted by PT or ET.
 *  Only calorimeter cells and basic cluster are sorted by E
 *  TRootGenParticle particle is not sorted by default
 *  to preserve mother-dautherlinks between particles.
 *
 *  $Date: 2009-01-12 23:36:06 $
 *  $Revision: 1.6 $
 *
 *  
 *  \author S. Ovyn - UCL, Louvain-la-Neuve
 *
 */

#include "Utilities/ExRootAnalysis/interface/BlockClasses.h"
#include "Utilities/ExRootAnalysis/interface/BlockCompare.h"

#include "TLorentzVector.h"

TCompare *TRootGenParticle::fgCompare = 0;
TCompare *TRootJet::fgCompare = TComparePT<TRootJet>::Instance();
TCompare *TRootElectron::fgCompare = TComparePT<TRootElectron>::Instance();
TCompare *TRootMuon::fgCompare = TComparePT<TRootMuon>::Instance();
TCompare *TRootPhoton::fgCompare = TComparePT<TRootPhoton>::Instance();
TCompare *TRootTauJet::fgCompare = TComparePT<TRootTauJet>::Instance();
TCompare *TRootTracks::fgCompare = TComparePT<TRootTracks>::Instance();
TCompare *TRootETmis::fgCompare = 0;
TCompare *TRootCalo::fgCompare = 0;
TCompare *TRootZdcHits::fgCompare = 0;
TCompare *TRootGenEvent:: fgCompare = 0;
TCompare *TRootParticle::fgCompare=0;
TCompare *TRootLHEFParticle::fgCompare = 0;
TCompare *TRootRomanPotHits::fgCompare =0;
TCompare *TRootTrigger::fgCompare =0;

void TRootParticle::Set(const TLorentzVector& momentum) {
	E   = momentum.E();
	Px  = momentum.Px();
	Py  = momentum.Py();
	Pz  = momentum.Pz();
	PT  = momentum.Pt();
	Eta = momentum.Eta();
	Phi = momentum.Phi();
}

void TRootParticle::Set(const float px, const float py, const float pz, const float e) {
        TLorentzVector toFill;
        toFill.SetPxPyPzE(px,py,pz,e);
        E   = e;
        Px  = px;
        Py  = py;
        Pz  = pz;
        PT  = toFill.Pt();
        Eta = toFill.Eta();
        Phi = toFill.Phi();
}