AlphaTVarProducer Class Reference

#include <AlphaTVarProducer.h>

Inheritance diagram for AlphaTVarProducer:

Public Member Functions
	AlphaTVarProducer (const edm::ParameterSet &)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~AlphaTVarProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Member Functions
double	CalcAlphaT (std::vector< TLorentzVector >)

Static Private Member Functions
static double	alphaT (const double, const double, const double)
static double	CalcHT (const std::vector< TLorentzVector >)
static double	CalcMHT (const std::vector< TLorentzVector >)
static double	deltaHt (const std::vector< double > &)

Private Attributes
edm::InputTag	inputJetTag_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
typedef WorkerT< EDProducer >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 15 of file AlphaTVarProducer.h.

Constructor & Destructor Documentation

AlphaTVarProducer::AlphaTVarProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 26 of file AlphaTVarProducer.cc.

References edm::InputTag::encode(), inputJetTag_, and LogDebug.

                                                                   :
  inputJetTag_    (iConfig.getParameter<edm::InputTag>("inputJetTag")){

  produces<std::vector<double> >();

  LogDebug("") << "Inputs: "
           << inputJetTag_.encode() << " ";
}

AlphaTVarProducer::~AlphaTVarProducer

(

)

Definition at line 35 of file AlphaTVarProducer.cc.

{
}

Member Function Documentation

double AlphaTVarProducer::alphaT ( const double  HT, const double  DHT, const double  MHT  )

staticprivate

Definition at line 103 of file AlphaTVarProducer.cc.

References mathSSE::sqrt().

Referenced by CalcAlphaT(), and produce().

                                                                                    {
    return 0.5 * ( HT - DHT ) / sqrt( HT*HT - MHT*MHT );
}

double AlphaTVarProducer::CalcAlphaT ( std::vector< TLorentzVector >  jets )

private

Definition at line 72 of file AlphaTVarProducer.cc.

References alphaT(), CalcHT(), CalcMHT(), deltaHt(), reco::tau::disc::Eta(), and i.

Referenced by produce().

                                                           {
    std::vector<double> ETs;
    TVector3 MHT = CalcMHT(jets);
    float HT = CalcHT(jets);
    //float HT = 0;
    for(unsigned int i = 0; i < jets.size(); i++){
        if(jets[i].Et() > 50. && fabs(jets[i].Eta()) < 2.5) ETs.push_back(jets[i].Et());
        //HT += jets[i].Et();
    }
    if(ETs.size() < 2.) return 0.0;
    if(ETs.size() > 16.) return 0.0;
    float DHT = deltaHt(ETs);
    
    float AlphaT = alphaT(HT,DHT,MHT.Mag());
    
    return AlphaT;
    
}

double AlphaTVarProducer::CalcHT ( const std::vector< TLorentzVector >  jets )

staticprivate

Definition at line 108 of file AlphaTVarProducer.cc.

References reco::tau::disc::Eta(), and i.

Referenced by CalcAlphaT(), and produce().

                                                       {

    double HT = 0;
    for(unsigned int i = 0; i < jets.size(); i++){
        if(jets[i].Et() > 50. && fabs(jets[i].Eta()) < 2.5) HT += jets[i].Et();
    }
    
    return HT;
    
}

double AlphaTVarProducer::CalcMHT ( const std::vector< TLorentzVector >  jets )

staticprivate

Definition at line 119 of file AlphaTVarProducer.cc.

References reco::tau::disc::Eta(), and i.

Referenced by CalcAlphaT().

                                                               {
    TVector3 MHT;
    for(unsigned int i = 0; i < jets.size(); i++){
        if(jets[i].Et() > 50. && fabs(jets[i].Eta()) < 2.5) MHT -= jets[i].Vect();
    }
    MHT.SetZ(0.0);
    return MHT.Mag();
}

double AlphaTVarProducer::deltaHt ( const std::vector< double > &  ETs )

staticprivate

Definition at line 91 of file AlphaTVarProducer.cc.

References diffTreeTool::diff, i, j, and min.

Referenced by CalcAlphaT().

                                                              {
  if(ETs.size() > 16.) return 9999999;
    std::vector<double> diff( 1<<(ETs.size()-1) , 0. );
    for(unsigned i=0; i < diff.size(); i++)
        for(unsigned j=0; j < ETs.size(); j++)
            diff[i] += ETs[j] * ( 1 - 2 * (int(i>>j)&1) ) ;
    std::vector<double>::const_iterator it;
    double min=9999999;
    for(it = diff.begin(); it !=diff.end(); it++) if(*it<min) min = *it;
    return min;
}

void AlphaTVarProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDProducer.

Definition at line 41 of file AlphaTVarProducer.cc.

References alphaT(), CalcAlphaT(), CalcHT(), edm::Event::getByLabel(), inputJetTag_, j, edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, and query::result.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

{
   using namespace std;
   using namespace edm;
   using namespace reco;


   // get hold of collection of objects
   edm::Handle<reco::CaloJetCollection> calojet_handle; 
   iEvent.getByLabel(inputJetTag_,calojet_handle);

   std::auto_ptr<std::vector<double> > result(new std::vector<double>); 
   // check the the input collections are available
   if (calojet_handle.isValid()){
     std::vector<TLorentzVector> myJets;
     reco::CaloJetCollection::const_iterator jetIt;
     for(jetIt = calojet_handle->begin(); jetIt != calojet_handle->end(); jetIt++){
       TLorentzVector j; j.SetPtEtaPhiE(jetIt->pt(),jetIt->eta(), jetIt->phi(), jetIt->energy());
       myJets.push_back(j);
     }

     double alphaT = CalcAlphaT(myJets);
     double HT     = CalcHT(myJets);

     result->push_back(alphaT);
     result->push_back(HT);
   }
   iEvent.put(result);
}

Member Data Documentation

edm::InputTag AlphaTVarProducer::inputJetTag_

private

Definition at line 29 of file AlphaTVarProducer.h.

Referenced by AlphaTVarProducer(), and produce().