---

ExtKtJetAlgorithmWrapper Class Reference

#include <RecoJets/JetAlgorithms/interface/ExtKtJetAlgorithmWrapper.h>

List of all members.

Public Member Functions
	ExtKtJetAlgorithmWrapper (const edm::ParameterSet &ps)
	ExtKtJetAlgorithmWrapper ()
void	run (const JetReco::InputCollection &fInput, JetReco::OutputCollection *fOutput) const
	~ExtKtJetAlgorithmWrapper ()
Private Attributes
int	theAngle
int	theColType
double	theDcut
int	theMode
int	theNjets
double	thePtMin
int	theRecom
double	theRparam

---

Detailed Description

Definition at line 25 of file ExtKtJetAlgorithmWrapper.h.

---

Constructor & Destructor Documentation

ExtKtJetAlgorithmWrapper::ExtKtJetAlgorithmWrapper

(

)

Definition at line 29 of file ExtKtJetAlgorithmWrapper.cc.

                                                  {
  
}

ExtKtJetAlgorithmWrapper::ExtKtJetAlgorithmWrapper

(

const edm::ParameterSet &

ps

)

Definition at line 33 of file ExtKtJetAlgorithmWrapper.cc.

References edm::ParameterSet::getParameter(), theAngle, theColType, theDcut, theMode, theNjets, thePtMin, theRecom, and theRparam.

{
  
  //Getting information from the ParameterSet (reading config files):
  theRparam=ps.getParameter<double>("ExtKtRParam");
  thePtMin=ps.getParameter<double>("PtMin");
  theDcut=ps.getParameter<double>("dcut");
  theNjets=ps.getParameter<int>("njets");
  theAngle=ps.getParameter<int>("KtAngle");
  theRecom=ps.getParameter<int>("KtRecom");
  
  
  theColType=4; //set to pp-collision

  //configuring algorithm 
  theMode=-1;
  
  if ((theDcut==-1)&&(theNjets==-1)) theMode=0;
  else if ((theDcut!=-1)&&(theNjets==-1)&&(thePtMin==-1)) theMode=1;
  else if ((theDcut==-1)&&(theNjets!=-1)&&(thePtMin==-1)) theMode=2;
  else LogWarning("ExtKtJetDefinition")<<"[ExtKtJetWrapper] Wrong Configuration! - will not produce KtJets!";
}

ExtKtJetAlgorithmWrapper::~ExtKtJetAlgorithmWrapper

(

)

[inline]

Definition at line 30 of file ExtKtJetAlgorithmWrapper.h.

{};

---

Member Function Documentation

void ExtKtJetAlgorithmWrapper::run	(	const JetReco::InputCollection &	fInput,
		JetReco::OutputCollection *	fOutput
	)			const

Definition at line 56 of file ExtKtJetAlgorithmWrapper.cc.

References p, p4, theAngle, theColType, theDcut, theJets, theMode, theNjets, thePtMin, theRecom, theRparam, and StDecayID::tID.

Referenced by cms::ExtKtPilupSubtractionJetProducer::runAlgorithm().

                                                                           {
  if (theMode==-1) return;
   int index_=0;
   vector<KtJet::KtLorentzVector> theInput;
   //map<KtJet::KtLorentzVector,const std::vector<FJCand>> theMap;
   map<const unsigned int,JetReco::InputItem> theIDMap;
   for (JetReco::InputCollection::const_iterator inputCand=fInput.begin();
        inputCand!=fInput.end();inputCand++){
      double px=(*inputCand)->px();
      double py=(*inputCand)->py();
      double pz=(*inputCand)->pz();
      double E=(*inputCand)->energy();
      KtJet::KtLorentzVector p(px,py,pz,E);
      //p.set_user_index(index_);
      const unsigned int lvID=p.getID();
      theIDMap[lvID]=(*inputCand);
      theInput.push_back(p);
      index_++;
   }
   
   //Construct KtEvent
   KtJet::KtEvent *ktev;
   if (theMode==0) ktev=new KtJet::KtEvent(theInput,theColType,theAngle,theRecom,theRparam);
   else ktev=new KtJet::KtEvent(theInput,theColType,theAngle,theRecom);

   if (theMode==1) ktev->findJetsD(theDcut);
   else if (theMode==2) ktev->findJetsN(theNjets);
   
   //getJets
   vector<KtJet::KtLorentzVector>* theJets=new vector<KtJet::KtLorentzVector>(ktev->getJetsPt());
   


   // run the jet clustering with the above jet definition
   
  //make CMSSW-Objects:
   JetReco::InputCollection jetConst;
   for (std::vector<KtJet::KtLorentzVector>::const_iterator itJet=theJets->begin();
        itJet!=theJets->end();itJet++){
     const std::vector<const KtJet::KtLorentzVector*>  jet_constituents = (*itJet).getConstituents();
     for (std::vector<const KtJet::KtLorentzVector*>::const_iterator itConst=jet_constituents.begin();
          itConst!=jet_constituents.end();++itConst){
       const unsigned int tID=(*itConst)->getID();
       jetConst.push_back(theIDMap[tID]);
     }
     double px=(*itJet).px();
     double py=(*itJet).py();
     double pz=(*itJet).pz();
     double E=(*itJet).e();
     math::XYZTLorentzVector p4(px,py,pz,E);
     if (p4.Pt()>=thePtMin) fOutput->push_back(ProtoJet(p4,jetConst));
     jetConst.clear();
   }
   delete theJets;
   delete ktev;
   
}

---

Member Data Documentation

int ExtKtJetAlgorithmWrapper::theAngle [private]

Definition at line 37 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

int ExtKtJetAlgorithmWrapper::theColType [private]

Definition at line 41 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

double ExtKtJetAlgorithmWrapper::theDcut [private]

Definition at line 35 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

int ExtKtJetAlgorithmWrapper::theMode [private]

Definition at line 33 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

int ExtKtJetAlgorithmWrapper::theNjets [private]

Definition at line 36 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

double ExtKtJetAlgorithmWrapper::thePtMin [private]

Definition at line 34 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

int ExtKtJetAlgorithmWrapper::theRecom [private]

Definition at line 38 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

double ExtKtJetAlgorithmWrapper::theRparam [private]

Definition at line 39 of file ExtKtJetAlgorithmWrapper.h.

Referenced by ExtKtJetAlgorithmWrapper(), and run().

---

The documentation for this class was generated from the following files:

• RecoJets/JetAlgorithms/interface/ExtKtJetAlgorithmWrapper.h
• RecoJets/JetAlgorithms/src/ExtKtJetAlgorithmWrapper.cc

---