TauJetCorrector Class Reference

#include <TauJetCorrector.h>

Inheritance diagram for TauJetCorrector:

Classes
class	ParametrizationTauJet

Public Member Functions
double	correction (const LorentzVector &fJet) const override
	get correction using Jet information only More...
double	correction (const reco::Jet &) const override
	apply correction using Jet information only More...
bool	eventRequired () const override
	if correction needs event information More...
void	setParameters (std::string, int)
	TauJetCorrector (const edm::ParameterSet &fParameters)
	~TauJetCorrector () override
Public Member Functions inherited from JetCorrector
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, ParametrizationTauJet * >	ParametersMap

Private Attributes
ParametersMap	parametrization
int	type

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

jet energy corrections from Taujet calibration

Definition at line 13 of file TauJetCorrector.h.

Member Typedef Documentation

typedef std::map<double,ParametrizationTauJet *> TauJetCorrector::ParametersMap

private

Definition at line 45 of file TauJetCorrector.h.

Constructor & Destructor Documentation

TauJetCorrector::TauJetCorrector

(

const edm::ParameterSet &

fParameters

)

Definition at line 99 of file TauJetCorrector.cc.

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

{
  type = fConfig.getParameter<int>("TauTriggerType");
  setParameters (fConfig.getParameter <std::string> ("tagName"),type);
}

TauJetCorrector::~TauJetCorrector ( )

override

Definition at line 106 of file TauJetCorrector.cc.

References stringResolutionProvider_cfi::parametrization.

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

Member Function Documentation

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

overridevirtual

get correction using Jet information only

Implements JetCorrector.

Definition at line 178 of file TauJetCorrector.cc.

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

{
  //cout<<" Start Apply Corrections "<<endl;
  if(parametrization.empty()) { return 1.; }
  
    double et=fJet.Et();
    double eta=fabs(fJet.Eta());
    
    //cout<<" Et and eta of jet "<<et<<" "<<eta<<endl;

    double etnew;
    std::map<double,ParametrizationTauJet*>::const_iterator ip=parametrization.upper_bound(eta);
    etnew=(--ip)->second->value(et,eta);

     //cout<<" The new energy found "<<etnew<<" "<<et<<endl;

         float mScale = etnew/et;
            
     return mScale;
}

double TauJetCorrector::correction ( const reco::Jet &  fJet ) const

overridevirtual

apply correction using Jet information only

Implements JetCorrector.

Definition at line 199 of file TauJetCorrector.cc.

References reco::LeafCandidate::p4().

                                                            {
    return correction(fJet.p4());
}

bool TauJetCorrector::eventRequired ( ) const

inlineoverridevirtual

if correction needs event information

Implements JetCorrector.

Definition at line 24 of file TauJetCorrector.h.

{return false;}

void TauJetCorrector::setParameters	(	std::string	aCalibrationType,
		int	itype
	)

Definition at line 111 of file TauJetCorrector.cc.

References JetCalibrationParameterSetTauJet::eta(), createfilelist::int, JetCalibrationParameterSetTauJet::neta(), JetCalibrationParameterSetTauJet::parameters(), stringResolutionProvider_cfi::parametrization, muonDTDigis_cfi::pset, edm::second(), JetCalibrationParameterSetTauJet::type(), and JetCalibrationParameterSetTauJet::valid().

{
     //cout<< " Start to set parameters "<<endl;
     type = itype;
     JetCalibrationParameterSetTauJet pset(aCalibrationType);
     
     if((!pset.valid()) && (aCalibrationType!="no"))
       {
     edm::LogError( "TauJetCorrector: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;
     map<int,double > etaboundx;
     int iq = 0;
     int ig = 0;
     int iqcd = 0;
     int mtype = 0;
    for(int ieta=0; ieta<pset.neta();ieta++)
    {
     if( pset.type(ieta) == 1 ) {pq[iq] = pset.parameters(ieta);    iq++;  mtype=(int)(pset.type(ieta));}
     if( pset.type(ieta) == 2 ) {pg[ig] = pset.parameters(ieta);    ig++;  mtype=(int)(pset.type(ieta));}
     if( pset.type(ieta) == 3 ) {pqcd[iqcd] = pset.parameters(ieta);iqcd++;mtype=(int)(pset.type(ieta));}
     if( pset.type(ieta) == -1 ) {etaboundx[mtype-1] = pset.eta(ieta);}
    }
    
    //cout<<" Number of parameters "<<iq<<" "<<ig<<" "<<iqcd<<endl;
    int mynum = 0;
      for(int ieta=0; ieta<pset.neta();ieta++)
        {
    //cout<<" New parmetrization "<<ieta<<" "<<pset.type(ieta)<<endl;
    
     if ( pset.type(ieta) == -1 ) continue;
     if( ieta < iq+1)
     {
          parametrization[pset.eta(ieta)]=new ParametrizationTauJet(pset.type(ieta),(*pq.find(ieta)).second,
                                                                              (*etaboundx.find(0)).second);
       //cout<<" ALL "<<ieta<<" "<<((*pq.find(ieta)).second)[0]<<" "<<((*pq.find(ieta)).second)[1]<<" "<<
       //((*pq.find(ieta)).second)[2]<<endl;

     }
     if( ieta > iq && ieta < iq + ig + 2 ) 
     {
      mynum = ieta - iq - 1;
          parametrization[pset.eta(ieta)]=new ParametrizationTauJet(pset.type(ieta),(*pg.find(mynum)).second,
                                                                              (*etaboundx.find(1)).second);
       //cout<<" One prong "<<((*pg.find(mynum)).second)[0]<<" "<<((*pg.find(mynum)).second)[1]<<" "<<
       //((*pg.find(mynum)).second)[2]<<endl;
       
     }
     if( ieta > iq + ig + 1) 
     {
      mynum = ieta - iq - ig - 2;
      //cout<<" Mynum "<<mynum<<" "<<ieta<<" "<<pset.type(ieta)<<endl;
          parametrization[pset.eta(ieta)]=new ParametrizationTauJet(pset.type(ieta),(*pqcd.find(mynum)).second,
                                                                              (*etaboundx.find(2)).second);
       //cout<<" Two prongs "<<((*pqcd.find(mynum)).second)[0]<<" "<<((*pqcd.find(mynum)).second)[1]<<" "<<
       //((*pqcd.find(mynum)).second)[2]<<endl;
     }
    }
   //cout<<" Parameters inserted into mAlgorithm "<<endl;                                  
}

Member Data Documentation

ParametersMap TauJetCorrector::parametrization

private

Definition at line 46 of file TauJetCorrector.h.

int TauJetCorrector::type

private

Definition at line 47 of file TauJetCorrector.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputCommands(), core.autovars.NTupleVariable::makeBranch(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and core.AutoHandle.AutoHandle::ReallyLoad().