#include <QjetsAdder.h>

Inheritance diagram for QjetsAdder:

Public Member Functions
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup)

	QjetsAdder (const edm::ParameterSet &iConfig)

virtual	~QjetsAdder ()

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 Attributes
double	cutoff_

double	jetRad_

std::string	mJetAlgo_

int	ntrial_

QjetsPlugin	qjetsAlgo_

int	QjetsPreclustering_

edm::Service < edm::RandomNumberGenerator >	rng_

edm::InputTag	src_

edm::EDGetTokenT< edm::View < reco::Jet > >	src_token_

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 13 of file QjetsAdder.h.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 15 of file QjetsAdder.h.

                                                       :
     src_(iConfig.getParameter<edm::InputTag>("src")),
     src_token_(consumes<edm::View<reco::Jet>>(src_)),
     qjetsAlgo_( iConfig.getParameter<double>("zcut"),
         iConfig.getParameter<double>("dcutfctr"),
         iConfig.getParameter<double>("expmin"),
         iConfig.getParameter<double>("expmax"),
         iConfig.getParameter<double>("rigidity")),
     ntrial_(iConfig.getParameter<int>("ntrial")),
     cutoff_(iConfig.getParameter<double>("cutoff")), 
     jetRad_(iConfig.getParameter<double>("jetRad")), 
     mJetAlgo_(iConfig.getParameter<std::string>("jetAlgo")) ,
     QjetsPreclustering_(iConfig.getParameter<int>("preclustering")) 
   {
     produces<edm::ValueMap<float> >("QjetsVolatility");
   }

virtual QjetsAdder::~QjetsAdder ( )

inlinevirtual

Definition at line 32 of file QjetsAdder.h.

32 {}

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 9 of file QjetsAdder.cc.

References edm::hlt::Exception, edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), edm::RandomNumberGenerator::getEngine(), reco::Jet::getJetConstituents(), i, cuy::ii, edm::helper::Filler< Map >::insert(), fwrapper::jets, roll_playback::k, m, timingPdfMaker::mean, reco::LeafCandidate::pt(), edm::Event::put(), mathSSE::sqrt(), and edm::Event::streamID().

                                                                         {
   // read input collection
   //edm::Handle<edm::View<pat::Jet> > jets;
   edm::Handle<edm::View<reco::Jet> > jets;
   iEvent.getByToken(src_token_, jets);
 
   // prepare room for output
   std::vector<float> QjetsVolatility;       QjetsVolatility.reserve(jets->size());
 
   for ( typename edm::View<reco::Jet>::const_iterator jetIt = jets->begin() ; jetIt != jets->end() ; ++jetIt ) {
     reco::Jet newCand(*jetIt);
 
     if(newCand.pt()<cutoff_)
     {
       QjetsVolatility.push_back(-1);
       continue;
     }
 
     //refill and recluster
     vector<fastjet::PseudoJet> allconstits;
     //for (unsigned k=0; k < newCand.getPFConstituents().size(); k++){
     for (unsigned k=0; k < newCand.getJetConstituents().size(); k++){
       const edm::Ptr<reco::Candidate> thisParticle = newCand.getJetConstituents().at(k);
       allconstits.push_back( fastjet::PseudoJet( thisParticle->px(), thisParticle->py(), thisParticle->pz(), thisParticle->energy() ) );
     }
 
     fastjet::JetDefinition jetDef(fastjet::cambridge_algorithm, jetRad_);
     if (mJetAlgo_== "AK") jetDef.set_jet_algorithm( fastjet::antikt_algorithm );
     else if (mJetAlgo_ == "CA") jetDef.set_jet_algorithm( fastjet::cambridge_algorithm );
     else throw cms::Exception("GroomedJetFiller") << " unknown jet algorithm " << std::endl;
 
     fastjet::ClusterSequence thisClustering_basic(allconstits, jetDef);
     std::vector<fastjet::PseudoJet> out_jets_basic = sorted_by_pt(thisClustering_basic.inclusive_jets(cutoff_));
     //std::cout << newCand.pt() << " " << out_jets_basic.size() <<std::endl;
     if(out_jets_basic.size()==0){ // jet reclustering failed, most likely due to the higher cutoff. Use a recognizeable default value for this jet
       QjetsVolatility.push_back(-1);
       continue;
     }
     assert( out_jets_basic.size()==1 ); // jet reclustering of one jet should yield exactly one jet at this stage
 
     // setup objects for qjets computation
     fastjet::JetDefinition qjet_def(&qjetsAlgo_);
 
     std::vector<double> qjetmass;
 
     vector<fastjet::PseudoJet> constits;
     unsigned int nqjetconstits = out_jets_basic.at(0).constituents().size(); // there should always be exactly one reclsutered jet => always "at(0)"
     if (nqjetconstits < (unsigned int) QjetsPreclustering_) constits = out_jets_basic.at(0).constituents();
     else constits = out_jets_basic.at(0).associated_cluster_sequence()->exclusive_subjets_up_to(out_jets_basic.at(0),QjetsPreclustering_);
 
     edm::Service<edm::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine* engine = &rng->getEngine(iEvent.streamID());
     qjetsAlgo_.SetRNEngine(engine);
     // create probabilistic recusterings
     for(unsigned int ii = 0 ; ii < (unsigned int) ntrial_ ; ii++){
       //qjetsAlgo_.SetRandSeed(iEvent.id().event()*100 + (jetIt - jets->begin())*ntrial_ + ii );// set random seed for reprudcibility. We need a smarted scheme
       fastjet::ClusterSequence qjet_seq(constits, qjet_def);
       vector<fastjet::PseudoJet> inclusive_jets2 = sorted_by_pt(qjet_seq.inclusive_jets(cutoff_));
       if (inclusive_jets2.size()>0){ // fill the massvalue only if the reclustering was successfull
     qjetmass.push_back(inclusive_jets2[0].m());
       }
 
     }//end loop on trials
 
     if (qjetmass.size()<1) {//protection against dummy case
       QjetsVolatility.push_back(-1);
       continue;
     }
 
     double mean = std::accumulate( qjetmass.begin( ) , qjetmass.end( ) , 0 ) /qjetmass.size() ;
     float totalsquared = 0.;
     for (unsigned int i = 0; i < qjetmass.size(); i++){
       totalsquared += (qjetmass[i] - mean)*(qjetmass[i] - mean) ;
     }
     float variance = sqrt( totalsquared/qjetmass.size() );  
     
     QjetsVolatility.push_back(variance/mean);
   }//end loop on jets
 
   std::auto_ptr<edm::ValueMap<float> > outQJV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerQJV(*outQJV);
   fillerQJV.insert(jets, QjetsVolatility.begin(), QjetsVolatility.end());
   fillerQJV.fill();
 
   iEvent.put(outQJV,"QjetsVolatility");
 }

Member Data Documentation

double QjetsAdder::cutoff_

private

Definition at line 41 of file QjetsAdder.h.

double QjetsAdder::jetRad_

private

Definition at line 42 of file QjetsAdder.h.

std::string QjetsAdder::mJetAlgo_

private

Definition at line 43 of file QjetsAdder.h.

int QjetsAdder::ntrial_

private

Definition at line 40 of file QjetsAdder.h.

QjetsPlugin QjetsAdder::qjetsAlgo_

private

Definition at line 39 of file QjetsAdder.h.

int QjetsAdder::QjetsPreclustering_

private

Definition at line 44 of file QjetsAdder.h.

edm::Service<edm::RandomNumberGenerator> QjetsAdder::rng_

private

Definition at line 45 of file QjetsAdder.h.

edm::InputTag QjetsAdder::src_

private

Definition at line 37 of file QjetsAdder.h.

edm::EDGetTokenT<edm::View<reco::Jet> > QjetsAdder::src_token_

private

Definition at line 38 of file QjetsAdder.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation