JetPartonCorrector Class Reference

#include <JetPartonCorrector.h>

Inheritance diagram for JetPartonCorrector:

Public Member Functions
virtual double	correction (const LorentzVector &fJet) const
	get correction using Jet information only More...
virtual bool	eventRequired () const
	if correction needs event information More...
	JetPartonCorrector (const edm::ParameterSet &fConfig)
void	setParameters (std::string aCalibrationType, double aJetFinderRadius, int aPartonMixture)
virtual	~JetPartonCorrector ()
Public Member Functions inherited from JetCorrector
virtual double	correction (const reco::Jet &fJet) const =0
	apply correction using Jet information only More...
virtual double	correction (const reco::Jet &fJet, const edm::Event &fEvent, const edm::EventSetup &fSetup) const
	apply correction using all event information More...
virtual double	correction (const reco::Jet &fJet, const edm::RefToBase< reco::Jet > &fJetRef, const edm::Event &fEvent, const edm::EventSetup &fSetup) const
	apply correction using all event information More...
virtual double	correction (const reco::Jet &fJet, const edm::RefToBase< reco::Jet > &fJetRef, const edm::Event &fEvent, const edm::EventSetup &fSetup, LorentzVector &corrected) const
	Apply vectorial correction using all event information. More...
	JetCorrector ()
virtual bool	refRequired () const =0
	if correction needs the jet reference More...
virtual bool	vectorialCorrection () const
	if vectorial correction is provided More...
virtual	~JetCorrector ()

Private Types
typedef std::map< double, JetPartonNamespace::ParametrizationJetParton * >	ParametersMap

Private Attributes
ParametersMap	parametrization
double	theJetFinderRadius
int	thePartonMixture

Additional Inherited Members
Public Types inherited from JetCorrector
typedef reco::Particle::LorentzVector	LorentzVector
Static Public Member Functions inherited from JetCorrector
static const JetCorrector *	getJetCorrector (const std::string &fName, const edm::EventSetup &fSetup)
	retrieve corrector from the event setup. troughs exception if something is missing More...

Detailed Description

Definition at line 21 of file JetPartonCorrector.h.

Member Typedef Documentation

typedef std::map<double,JetPartonNamespace::ParametrizationJetParton *> JetPartonCorrector::ParametersMap

private

Definition at line 36 of file JetPartonCorrector.h.

Constructor & Destructor Documentation

JetPartonCorrector::JetPartonCorrector

(

const edm::ParameterSet &

fConfig

)

Definition at line 146 of file JetPartonCorrector.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
  thePartonMixture = fConfig.getParameter<int>("MixtureType");
  theJetFinderRadius = fConfig.getParameter<double>("Radius");
  setParameters (fConfig.getParameter <std::string> ("tagName"),theJetFinderRadius,thePartonMixture );
}

JetPartonCorrector::~JetPartonCorrector ( )

virtual

Definition at line 153 of file JetPartonCorrector.cc.

References stringResolutionProvider_cfi::parametrization.

{
  for(ParametersMap::iterator ip=parametrization.begin();ip!=parametrization.end();ip++) delete ip->second;  
}

Member Function Documentation

double JetPartonCorrector::correction ( const LorentzVector &  fJet ) const

virtual

get correction using Jet information only

Implements JetCorrector.

Definition at line 192 of file JetPartonCorrector.cc.

References gather_cfg::cout, stringResolutionProvider_cfi::et, stringResolutionProvider_cfi::eta, stringResolutionProvider_cfi::parametrization, and edm::second().

{
  if(parametrization.empty()) { return 1.; }
  
  double et=fJet.Et();
  double eta=fabs(fJet.Eta());
  
  //if(eta<10) { eta=abs(fJet.getY()); }

  double etnew;
  std::map<double,JetPartonNamespace::ParametrizationJetParton*>::const_iterator ip=parametrization.upper_bound(eta);
  if(ip==parametrization.begin()) 
    { 
      etnew=ip->second->value(et,eta); 
    }
  else if(ip==parametrization.end()) 
    {
      etnew=(--ip)->second->value(et,eta);
    }
  else
    {
      double eta2=ip->first;
      double et2=ip->second->value(et,eta);
      ip--;
      double eta1=ip->first;
      double et1=ip->second->value(et,eta);
      
      etnew=(eta2*et1 - eta1*et2 + eta*et2 - eta*et1)/(eta2-eta1); 
    }
   cout<<" JetParton::The new energy found "<<etnew<<" "<<et<<endl;    
  float mScale = 1000.;
   
  if( et > 0.001) mScale = etnew/et;

  return mScale;
  
  
}

virtual bool JetPartonCorrector::eventRequired ( ) const

inlinevirtual

if correction needs event information

Implements JetCorrector.

Definition at line 32 of file JetPartonCorrector.h.

{return false;}

void JetPartonCorrector::setParameters	(	std::string	aCalibrationType,
		double	aJetFinderRadius,
		int	aPartonMixture
	)

Definition at line 158 of file JetPartonCorrector.cc.

References JetPartonNamespace::JetPartonCalibrationParameterSet::eta(), JetPartonNamespace::JetPartonCalibrationParameterSet::neta(), JetPartonNamespace::JetPartonCalibrationParameterSet::parameters(), stringResolutionProvider_cfi::parametrization, muonDTDigis_cfi::pset, edm::second(), JetPartonNamespace::JetPartonCalibrationParameterSet::type(), and JetPartonNamespace::JetPartonCalibrationParameterSet::valid().

{
     
     theJetFinderRadius = aJetFinderRadius;
     thePartonMixture = aPartonMixture;
    
     JetPartonNamespace::JetPartonCalibrationParameterSet pset(aCalibrationType);
      
     if((!pset.valid()) && (aCalibrationType != "no"))
       {
     edm::LogError ("JetPartonCorrector: Jet Corrections not found ") << aCalibrationType << 
       " not found! Cannot apply any correction ... For JetPlusTrack calibration only radii 0.5 and 0.7 are included for JetParton" << endl;
          return;
       }
     if (aCalibrationType=="no") return;
       
       
     map<int,vector<double> > pq;
     map<int,vector<double> > pg;
     map<int,vector<double> > pqcd;
     int iq = 0;
     int ig = 0;
     int iqcd = 0;
    for(int ieta=0; ieta<pset.neta();ieta++)
    {
     if( pset.type(ieta) == 1 ) {pq[iq] = pset.parameters(ieta); iq++;};
     if( pset.type(ieta) == 2 ) {pg[ig] = pset.parameters(ieta);ig++;};
     if( pset.type(ieta) == 3 ) {pqcd[iqcd] = pset.parameters(ieta);iqcd++;};
    }
    
    for(int ieta=0; ieta<iq;ieta++){
      parametrization[pset.eta(ieta)]=new JetPartonNamespace::ParametrizationJetParton(thePartonMixture,(*pq.find(ieta)).second,(*pg.find(ieta)).second,(*pqcd.find(ieta)).second);    
    }
}

Member Data Documentation

ParametersMap JetPartonCorrector::parametrization

private

Definition at line 37 of file JetPartonCorrector.h.

double JetPartonCorrector::theJetFinderRadius

private

Definition at line 39 of file JetPartonCorrector.h.

int JetPartonCorrector::thePartonMixture

private

Definition at line 38 of file JetPartonCorrector.h.