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 ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
double	CalcAlphaT (const 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_
edm::EDGetTokenT < reco::CaloJetCollection >	inputJetTagToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
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::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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(), edm::ParameterSet::getParameter(), inputJetTag_, inputJetTagToken_, and LogDebug.

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

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

  //set Token(-s)
  inputJetTagToken_ = consumes<reco::CaloJetCollection>(iConfig.getParameter<edm::InputTag>("inputJetTag"));

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

AlphaTVarProducer::~AlphaTVarProducer

(

)

Definition at line 38 of file AlphaTVarProducer.cc.

{
}

Member Function Documentation

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

staticprivate

Definition at line 106 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 ( const std::vector< TLorentzVector > &  jets )

private

Definition at line 75 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),0.0,0.0};
    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 111 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 122 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 94 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 44 of file AlphaTVarProducer.cc.

References alphaT(), CalcAlphaT(), CalcHT(), edm::Event::getByToken(), inputJetTagToken_, 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.getByToken(inputJetTagToken_, 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().

edm::EDGetTokenT<reco::CaloJetCollection> AlphaTVarProducer::inputJetTagToken_

private

Definition at line 32 of file AlphaTVarProducer.h.

Referenced by AlphaTVarProducer(), and produce().