Inheritance diagram for ZMuMuUserData:

Public Types
typedef math::XYZVector	Vector

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	ZMuMuUserData (const edm::ParameterSet &)

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
double	alpha_

EDGetTokenT< BeamSpot >	beamSpotToken_

double	beta_

int	counter

string	hltPath_

EDGetTokenT< VertexCollection >	primaryVerticesToken_

EDGetTokenT< std::vector< reco::CompositeCandidate > >	srcToken_

EDGetTokenT< GenParticleMatch >	zGenParticlesMatchToken_

Additional Inherited Members
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::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 29 of file ZMuMuUserData.cc.

Member Typedef Documentation

typedef math::XYZVector ZMuMuUserData::Vector

Definition at line 32 of file ZMuMuUserData.cc.

Constructor & Destructor Documentation

ZMuMuUserData::ZMuMuUserData ( const edm::ParameterSet & cfg )

Definition at line 49 of file ZMuMuUserData.cc.

                                                       :
   srcToken_(consumes<std::vector<reco::CompositeCandidate> > ( cfg.getParameter<InputTag>( "src" ) ) ),
   beamSpotToken_(consumes<BeamSpot> (cfg.getParameter<InputTag>( "beamSpot" ) ) ),
   primaryVerticesToken_(consumes<VertexCollection> (cfg.getParameter<InputTag>( "primaryVertices" ) ) ),
   zGenParticlesMatchToken_(consumes<GenParticleMatch> (cfg.getParameter<InputTag>( "zGenParticlesMatch" ) ) ),
   alpha_(cfg.getParameter<double>("alpha") ),
   beta_(cfg.getParameter<double>("beta") ),
   hltPath_(cfg.getParameter<std::string >("hltPath") ){
   produces<vector<pat::CompositeCandidate> >();
 }

Member Function Documentation

void ZMuMuUserData::produce	(	edm::Event &	evt,
		const edm::EventSetup &
	)

overrideprivate

I fill the dau1 with positive and dau2 with negatove values for the pt, in order to flag the muons used for building zMassSa

I fill the dau1 with negatove and dau2 with positive values for the pt

Definition at line 60 of file ZMuMuUserData.cc.

References pat::PATObject< ObjectType >::addUserFloat(), beamSpotToken_, counter, reco::CompositeCandidate::daughter(), DEFINE_FWK_MODULE, reco::Candidate::eta(), reco::Candidate::get(), edm::Event::getByToken(), mps_fire::i, reco::Muon::isGlobalMuon(), edm::Ref< C, T, F >::isNonnull(), ResonanceBuilder::mass, reco::Candidate::mass(), reco::Candidate::masterClone(), eostools::move(), reco::Candidate::phi(), reco::Candidate::polarP4(), jets_cff::primaryVertices, primaryVerticesToken_, reco::Candidate::pt(), edm::Event::put(), reco::Candidate::rapidity(), srcToken_, z, and zGenParticlesMatchToken_.

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

                                                              {
   Handle<std::vector<reco::CompositeCandidate> > dimuons;
   evt.getByToken(srcToken_,dimuons);
 
   Handle<BeamSpot> beamSpotHandle;
   evt.getByToken(beamSpotToken_, beamSpotHandle);
 
   Handle<VertexCollection> primaryVertices;  // Collection of primary Vertices
   evt.getByToken(primaryVerticesToken_, primaryVertices);
 
 
   bool isMCMatchTrue=false;
 
   Handle<GenParticleMatch> zGenParticlesMatch;
   if(evt.getByToken( zGenParticlesMatchToken_, zGenParticlesMatch )){
     isMCMatchTrue=true;
   }
 
   //cout<<"isMCMatchTrue"<<isMCMatchTrue <<endl;
   unique_ptr<vector<pat::CompositeCandidate> > dimuonColl( new vector<pat::CompositeCandidate> () );
 
 
   for (unsigned int i = 0; i< dimuons->size();++i){
     const CompositeCandidate & z = (*dimuons)[i];
     //CandidateBaseRef zRef = dimuons ->refAt(i);
     edm::Ref<std::vector<reco::CompositeCandidate> > zRef(dimuons, i);
     pat::CompositeCandidate dimuon(z);
 
     float trueMass,truePt,trueEta,truePhi,trueY;
     if (isMCMatchTrue){
     GenParticleRef trueZRef  = (*zGenParticlesMatch)[zRef];
     //CandidateRef trueZRef = trueZIter->val;
     if( trueZRef.isNonnull() ) {
       const Candidate & z = * trueZRef;
       trueMass = z.mass();
       truePt   = z.pt();
       trueEta  = z.eta();
       truePhi  = z.phi();
       trueY    = z.rapidity();
     } else {
       trueMass = -100;
       truePt   = -100;
       trueEta  = -100;
       truePhi  = -100;
       trueY    = -100;
     }
 
     dimuon.addUserFloat("TrueMass",trueMass);
     dimuon.addUserFloat("TruePt",truePt);
     dimuon.addUserFloat("TrueEta",trueEta);
     dimuon.addUserFloat("TruePhi",truePhi);
     dimuon.addUserFloat("TrueY",trueY);
 
     }
     const Candidate * dau1 = z.daughter(0);
     const Candidate * dau2 = z.daughter(1);
     const pat::Muon & mu1 = dynamic_cast<const pat::Muon&>(*dau1->masterClone());
     const pat::Muon & mu2 = dynamic_cast<const pat::Muon&>(*dau2->masterClone());
 
     /*cout<<"mu1 is null? "<<mu1.isMuon()<<endl;
     cout<<"mu2 is null? "<<mu2.isMuon()<<endl;
     cout<<"mu1 is global?"<<mu1.isGlobalMuon()<<endl;
     cout<<"mu2 is global?"<<mu2.isGlobalMuon()<<endl;
     */
 
     if(mu1.isGlobalMuon()==true && mu2.isGlobalMuon()==true){
       TrackRef stAloneTrack1;
       TrackRef stAloneTrack2;
       Vector momentum;
       Candidate::PolarLorentzVector p4_1;
       double mu_mass;
       stAloneTrack1 = dau1->get<TrackRef,reco::StandAloneMuonTag>();
       stAloneTrack2 = dau2->get<TrackRef,reco::StandAloneMuonTag>();
       float zDau1SaEta = stAloneTrack1->eta();
       float zDau2SaEta = stAloneTrack2->eta();
       float zDau1SaPhi = stAloneTrack1->phi();
       float zDau2SaPhi = stAloneTrack2->phi();
       float zDau1SaPt,zDau2SaPt;
       if(counter % 2 == 0) {
     momentum = stAloneTrack1->momentum();
     p4_1 = dau2->polarP4();
     mu_mass = dau1->mass();
     zDau1SaPt = stAloneTrack1->pt();
     zDau2SaPt = - stAloneTrack2->pt();
       }else{
     momentum = stAloneTrack2->momentum();
     p4_1= dau1->polarP4();
     mu_mass = dau2->mass();
     zDau1SaPt = - stAloneTrack1->pt();
     zDau2SaPt =  stAloneTrack2->pt();
       }
 
       Candidate::PolarLorentzVector p4_2(momentum.rho(), momentum.eta(),momentum.phi(), mu_mass);
       double mass = (p4_1+p4_2).mass();
       float zMassSa = mass;
       //cout<<"zMassSa "<<zMassSa;
       dimuon.addUserFloat("MassSa",zMassSa);
       dimuon.addUserFloat("Dau1SaPt",zDau1SaPt);
       dimuon.addUserFloat("Dau2SaPt",zDau2SaPt);
       dimuon.addUserFloat("Dau1SaPhi",zDau1SaPhi);
       dimuon.addUserFloat("Dau2SaPhi",zDau2SaPhi);
       dimuon.addUserFloat("Dau1SaEta",zDau1SaEta);
       dimuon.addUserFloat("Dau2SaEta",zDau2SaEta);
       ++counter;
     }
     dimuonColl->push_back(dimuon);
 
   }
 
 
   evt.put(std::move(dimuonColl));
 }

Member Data Documentation

double ZMuMuUserData::alpha_

private

Definition at line 40 of file ZMuMuUserData.cc.

EDGetTokenT<BeamSpot> ZMuMuUserData::beamSpotToken_

private

Definition at line 37 of file ZMuMuUserData.cc.

Referenced by produce().

double ZMuMuUserData::beta_

private

Definition at line 40 of file ZMuMuUserData.cc.

int ZMuMuUserData::counter

private

Definition at line 42 of file ZMuMuUserData.cc.

Referenced by counter.Counters::__getitem__(), and produce().

string ZMuMuUserData::hltPath_

private

Definition at line 41 of file ZMuMuUserData.cc.

EDGetTokenT<VertexCollection> ZMuMuUserData::primaryVerticesToken_

private

Definition at line 38 of file ZMuMuUserData.cc.

Referenced by produce().

EDGetTokenT<std::vector<reco::CompositeCandidate> > ZMuMuUserData::srcToken_

private

Definition at line 36 of file ZMuMuUserData.cc.

Referenced by produce().

EDGetTokenT<GenParticleMatch> ZMuMuUserData::zGenParticlesMatchToken_

private

Definition at line 39 of file ZMuMuUserData.cc.

Referenced by produce().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation