TauJetCorrector Class Reference

#include <TauJetCorrector.h>

Inheritance diagram for TauJetCorrector:

Classes
class	ParametrizationTauJet

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

virtual

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

virtual

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

virtual

apply correction using Jet information only

Implements JetCorrector.

Definition at line 199 of file TauJetCorrector.cc.

References reco::LeafCandidate::p4().

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

virtual bool TauJetCorrector::eventRequired ( ) const

inlinevirtual

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().