---

Atlfast::JetMaker Class Reference

#include <JetMaker.h>

Collaboration diagram for Atlfast::JetMaker:

[legend]List of all members.

Public Methods
	JetMaker (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~JetMaker ()
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Jet *	create (Cluster *)
bool	isAcceptable (Jet *)
void	addMuons (Jet *, MuonCollection &)
MsgStream &	dumpParameters (MsgStream &) const
Private Types
typedef std::vector< Cluster * >	localClusterCollection
typedef localClusterCollection::iterator	localClusterIterator
typedef std::vector< ReconstructedParticle * >	MuonCollection
typedef MuonCollection::iterator	MuonIterator
Private Methods
int	getJetTags (MCparticleCollection &mcParticles, const PartitionCondition::BelowThresholdDeltaR &tagSelector) const
int	getTauTags (MCparticleCollection &mcParticles, const PartitionCondition::BelowThresholdDeltaR &tagSelector, double jetPt) const
void	tag (Jet *jet, MCparticleCollection &mcParticles_b, MCparticleCollection &mcParticles_c, MCparticleCollection &mcParticles_tau) const
HepLorentzVector	addUnusedCells (MsgStream &) const
Private Attributes
double	m_minPT
double	m_maxEta
bool	m_doSmearing
double	m_barrelForwardEta
double	m_rconeb
double	m_rconef
double	m_bPtMin
double	m_bMaxDeltaR
double	m_cPtMin
double	m_cMaxDeltaR
double	m_tauPtMin
double	m_tauMaxDeltaR
double	m_etaTagMax
double	m_tauJetPtRatio
ISmearer *	m_smearer
ISmearer *	m_cellSmearer
TesIO *	m_tesIO
std::string	m_inputLocation
std::string	m_outputLocation
std::string	m_muonLocation
std::string	m_mcTruthLocation
std::string	m_unusedCellLocation
std::string	m_missingMomentumLocation
KinematicHelper	m_kinehelp

---

Member Typedef Documentation

typedef std::vector<Cluster*> Atlfast::JetMaker::localClusterCollection [private]

Definition at line 98 of file JetMaker.h. Referenced by execute().

typedef localClusterCollection::iterator Atlfast::JetMaker::localClusterIterator [private]

Definition at line 99 of file JetMaker.h. Referenced by execute().

typedef std::vector<ReconstructedParticle*> Atlfast::JetMaker::MuonCollection [private]

Definition at line 100 of file JetMaker.h. Referenced by execute().

typedef MuonCollection::iterator Atlfast::JetMaker::MuonIterator [private]

Definition at line 101 of file JetMaker.h. Referenced by addMuons().

---

Constructor & Destructor Documentation

Atlfast::JetMaker::JetMaker (  const std::string &    name, ISvcLocator *    pSvcLocator )

Definition at line 83 of file JetMaker.cxx. : Algorithm( name, pSvcLocator ){ // Setting the parameter defaults. m_minPT = 10; // Min jet trans momentum m_maxEta = 5; //Max jet eta m_doSmearing = true; //smearing is done m_rconeb = 0.400; m_rconef = 0.400; m_bPtMin = 5.0; //min pt for bquark to tag jet m_bMaxDeltaR = 0.2; //max dR for bquark to tag jet m_cPtMin = 5.0; //min pt for cquark to tag jet m_cMaxDeltaR = 0.2; m_tauPtMin = 10.0; m_tauMaxDeltaR = 0.3; m_etaTagMax = 2.5; m_tauJetPtRatio = 0.9; // Default paths for entities in the TES m_inputLocation = "/Event/AtlfastClusters"; m_outputLocation = "/Event/AtlfastJets"; m_muonLocation = "/Event/AtlfastNonIsolatedMuons"; m_unusedCellLocation = "/Event/AtlfastUnusedCells"; m_missingMomentumLocation = "/Event/AtlfastMissingMomentum"; m_mcTruthLocation = "/Event/McEventCollection"; // This is how you declare the paramemters to Gaudi so that // they can be over-written via the job options file declareProperty( "MinimumPT", m_minPT ) ; declareProperty( "MaximumEta", m_maxEta ) ; declareProperty( "DoSmearing", m_doSmearing ); declareProperty( "RconeB", m_rconeb ); declareProperty( "RconeF", m_rconef ); declareProperty( "bPtMin", m_bPtMin ); declareProperty( "bMaxDeltaR", m_bMaxDeltaR ); declareProperty( "cPtMin", m_cPtMin ); declareProperty( "cMaxDeltaR", m_cMaxDeltaR ); declareProperty( "tauPtMin", m_tauPtMin ); declareProperty( "tauMaxDeltaR", m_tauMaxDeltaR ); declareProperty( "etaTagMax", m_etaTagMax ); declareProperty( "tauJetPtRatio", m_tauJetPtRatio); declareProperty( "InputLocation", m_inputLocation ) ; declareProperty( "OutputLocation", m_outputLocation ) ; declareProperty( "UnusedCellLocation", m_unusedCellLocation ) ; declareProperty( "MuonLocation", m_muonLocation ); declareProperty( "MissingMomentumLocation", m_missingMomentumLocation ); declareProperty( "McTruthLocation", m_mcTruthLocation ); }

Atlfast::JetMaker::~JetMaker (    )  [virtual]

Definition at line 139 of file JetMaker.cxx. References m_cellSmearer, m_smearer, and m_tesIO. { MsgStream log( messageService(), name() ) ; log << MSG::INFO << "Destructor Called" << endreq; if (m_smearer) { log << MSG::INFO << "Deleting jet Smearer" << endreq; delete m_smearer; } if (m_cellSmearer) { log << MSG::INFO << "Deleting Cell Smearer" << endreq; delete m_cellSmearer; } if (m_tesIO) { delete m_tesIO; } }

---

Member Function Documentation

int Atlfast::JetMaker::getJetTags (  MCparticleCollection &    mcParticles, const PartitionCondition::BelowThresholdDeltaR &    tagSelector )  const [private]

Definition at line 546 of file JetMaker.cxx. { MsgStream log( messageService(), name() ) ; MCparticleCollectionCIter src ; for( src = mcParticles.begin(); src != mcParticles.end(); ++src ){ if(tagSelector(SimpleKinematic(*src))) return (*src)->pdg_id(); } return 0; }

int Atlfast::JetMaker::getTauTags (  MCparticleCollection &    mcParticles, const PartitionCondition::BelowThresholdDeltaR &    tagSelector, double    jetPt )  const [private]

Definition at line 560 of file JetMaker.cxx. { MsgStream log( messageService(), name() ) ; MCparticleCollectionCIter src ; for( src = mcParticles.begin() ; src != mcParticles.end() ; ++src ) { double tauPt = (*src)->momentum().perp(); double neuPt=0.; // Find Tau children HepMC::GenVertex* endVertex = (*src)->end_vertex(); if(!endVertex) return 0; HepMC::GenVertex::particle_iterator firstChild = endVertex->particles_begin(HepMC::children); HepMC::GenVertex::particle_iterator endChild = endVertex->particles_end(HepMC::children); HepMC::GenVertex::particle_iterator thisChild = firstChild; // Find Pt of Taus tau-neutrino for(; thisChild!=endChild; ++thisChild){ if( abs((*thisChild)->pdg_id()) == 16){ neuPt = (*thisChild)->momentum().perp(); } } if(tagSelector(SimpleKinematic(*src)) && ((tauPt-neuPt)/jetPt) >= m_tauJetPtRatio) return (*src)->pdg_id(); } return 0; }

void Atlfast::JetMaker::tag (  Jet *    jet, MCparticleCollection &    mcParticles_b, MCparticleCollection &    mcParticles_c, MCparticleCollection &    mcParticles_tau )  const [private]

Definition at line 512 of file JetMaker.cxx. Referenced by execute(). { //See if jet can be b-tagged //Construct an object to select IKinematics within DeltaR of the current jet PartitionCondition::BelowThresholdDeltaR bTagSelector(jet, m_bMaxDeltaR); int btag = getJetTags(mcParticles_b,bTagSelector); if(btag != 0 ) { jet->setBTag(); return; } //else see if jet can be c-tagged //Construct an object to select IKinematics within DeltaR of the current jet PartitionCondition::BelowThresholdDeltaR cTagSelector(jet, m_cMaxDeltaR); int ctag = getJetTags( mcParticles_c, cTagSelector); if(ctag != 0 ) { jet->setCTag(); return; } //finally see if jet can be tau-tagged //Construct an object to select IKinematics within DeltaR of the current jet PartitionCondition::BelowThresholdDeltaR tauTagSelector(jet, m_tauMaxDeltaR); int tauTag=getTauTags( mcParticles_tau, tauTagSelector, jet->pT()); if(tauTag !=0 ) jet->setTauTag(tauTag); return; }

HepLorentzVector Atlfast::JetMaker::addUnusedCells (  MsgStream &    )  const [private]

Definition at line 598 of file JetMaker.cxx. References m_cellSmearer, m_tesIO, and m_unusedCellLocation. Referenced by execute(). { HepLorentzVector temp(0., 0., 0., 0.); DataHandle<ITwoCptCellCollection> unusedCells; if( ! m_tesIO->getDH( unusedCells, m_unusedCellLocation) ){ log << MSG::INFO << "No unused cells in TES " << endreq; return temp; } log << MSG::DEBUG << "No unused cells in TES " << endreq; return std::accumulate(unusedCells->begin(), unusedCells->end(), temp, SmearCell(m_cellSmearer) ); }

StatusCode Atlfast::JetMaker::initialize (    )

Definition at line 161 of file JetMaker.cxx. References Atlfast::GlobalEventData::barrelForwardEta(), Atlfast::GlobalEventData::lumi(), m_barrelForwardEta, m_bMaxDeltaR, m_bPtMin, m_cellSmearer, m_cMaxDeltaR, m_cPtMin, m_doSmearing, m_etaTagMax, m_inputLocation, m_maxEta, m_mcTruthLocation, m_minPT, m_missingMomentumLocation, m_muonLocation, m_outputLocation, m_rconeb, m_rconef, m_smearer, m_tauJetPtRatio, m_tauMaxDeltaR, m_tauPtMin, m_tesIO, and Atlfast::GlobalEventData::randSeed(). { MsgStream log( messageService(), name() ) ; log << MSG::DEBUG << "instantiating a JetSmearer" << endreq; // m_tesIO = new TesIO(eventDataService()); m_tesIO = new TesIO(); //get the Global Event Data using singleton pattern GlobalEventData* ged = GlobalEventData::Instance(); int lumi = ged->lumi(); int randSeed = ged->randSeed() ; m_barrelForwardEta = ged->barrelForwardEta(); if(m_doSmearing){ m_smearer = new JetSmearer(randSeed, lumi, m_rconeb, m_rconef, m_barrelForwardEta); m_cellSmearer = new JetSmearer(randSeed, lumi, 0., 0., m_barrelForwardEta); } else { m_smearer = NULL; m_cellSmearer = NULL; } HeaderPrinter hp("Atlfast JetMaker:", log); hp.add("TES Locations: "); hp.add(" Muons from ", m_muonLocation); hp.add(" Clusters from ", m_inputLocation); hp.add(" GeneratorParticles from ", m_mcTruthLocation); hp.add(" Jets to ", m_outputLocation); hp.add(" unused cell,cluster sum to ", m_missingMomentumLocation); hp.add("Smearing on ", m_doSmearing); hp.add("Cluster min Pt ", m_minPT); hp.add("Cluster max Eta ", m_maxEta); hp.add("Luminosity ", lumi); hp.add("Initial random seed ", randSeed); hp.add("End cap Cone size ", m_rconef); hp.add("Barrel Cone size ", m_rconeb); hp.add("Eta of start of end cap ", m_barrelForwardEta); hp.add("Parameters for tags ", " "); hp.add(" Max eta ", m_etaTagMax); hp.add(" b-quark min pt ", m_bPtMin); hp.add(" b-quark max DeltaR ", m_bMaxDeltaR); hp.add(" c-quark min pt ", m_cPtMin); hp.add(" c-quark max DeltaR ", m_cMaxDeltaR); hp.add(" tau lep min pt ", m_tauPtMin); hp.add(" tau lep max DeltaR ", m_tauMaxDeltaR); hp.add(" tau/Jet Pt Ratio ", m_tauJetPtRatio); hp.print(); return StatusCode::SUCCESS ; }

StatusCode Atlfast::JetMaker::execute (    )

Definition at line 252 of file JetMaker.cxx. References addMuons(), addUnusedCells(), create(), isAcceptable(), localClusterCollection, localClusterIterator, m_bPtMin, m_cPtMin, m_etaTagMax, m_inputLocation, m_missingMomentumLocation, m_muonLocation, m_outputLocation, m_tauPtMin, m_tesIO, MuonCollection, and tag(). { MsgStream log( messageService(), name() ) ; std::string mess; // make a message logging stream //............................. // Make some locals stores which be used to keep pointers to // all of the entities from the event store which are needed by // this algorithm. These are all local and are typedefed in this class localClusterCollection my_Clusters ; MuonCollection my_Muons ; //................................ // if( ! m_tesIO->copy( my_Clusters, m_inputLocation ) ) { if( ! m_tesIO->copy<ClusterCollection> ( my_Clusters, m_inputLocation ) ) { log << MSG::INFO << "No Clusters in TES " << endreq; } else { log << MSG::DEBUG << my_Clusters.size()<<" Clusters from TES " << endreq; } // get non-isolated muons for jet construction if( ! m_tesIO->copy<ReconstructedParticleCollection> ( my_Muons, m_muonLocation ) ) { log << MSG::DEBUG << "No muons in TES " << endreq; } // ...................... // Make a container object which will store pointers to all isolated // Jetss which are successfully created. Since it is going to go // into the event store, it has to be a special Gaudi type of container. // This is typedefined in the entity include file // As it is going into the TES it needs to be created on the heap with new. JetCollection* myJets = new JetCollection ; //......................................................... // From each Cluster create a reconstructed Jet. // If all requirements are satisfied add to the collection Jet* candidate ; localClusterIterator src ; HepLorentzVector missingMomentum(0.,0.,0.,0.); const int BQUARK(ParticleCodes::BQUARK); const int CQUARK(ParticleCodes::CQUARK); //Construct an object to select b, c and tau quarks with kine conditions const HepMC_helper::SelectJetTag bSelector(BQUARK, m_bPtMin, m_etaTagMax); const HepMC_helper::SelectJetTag cSelector(CQUARK, m_cPtMin, m_etaTagMax); const HepMC_helper::SelectTauTag tauSelector(m_tauPtMin, m_etaTagMax); MCparticleCollection mcParticles_b ; MCparticleCollection mcParticles_c ; MCparticleCollection mcParticles_tau ; //get quarks used in tagging from TES using tesio m_tesIO->getMC( mcParticles_b, &bSelector ) ; m_tesIO->getMC( mcParticles_c, &cSelector ) ; m_tesIO->getMC( mcParticles_tau, &tauSelector ) ; for( src = my_Clusters.begin() ; src != my_Clusters.end() ; ++src ) { IAssociationManager* iAssocMan = *src; ReconstructedParticle rp; if (iAssocMan->associations(rp).size() == 0) { // Make a new jet candidate = this->create( *src ); //Make a temp HepLorentzVector from candidates 4-mom //before muon is added in HepLorentzVector temp4Vec = candidate->momentum(); // add isolated muons this->addMuons(candidate, my_Muons); // add to missing momentum vector if not used in a jet. // remember, these are smeared momenta if smearing is enabled if( this->isAcceptable( candidate ) ) { // if the jet is acceptable flavor tag the jet this->tag(candidate,mcParticles_b,mcParticles_c,mcParticles_tau); log << MSG::DEBUG << "Pushing back a Jet" << endreq; myJets->push_back( candidate ) ; log<<MSG::DEBUG<<"jet "<<*candidate<<endreq; }else { log << MSG::DEBUG << "Jet fails cuts" << endreq; //add in temp4Vec to avoid muon double counting missingMomentum += temp4Vec; delete candidate ; } } } log<<MSG::DEBUG<<"Missing momentum: clusters "<<missingMomentum<<endreq; missingMomentum += addUnusedCells(log); log<<MSG::DEBUG<<"Missing momentum: cells+clstrs " <<missingMomentum<<endreq; //...................................... // Register any Jets which have been successfilly created in the // transient event store. TesIoStat stat; log<<MSG::DEBUG<<"Storing "<< myJets->size() <<" Jets"<<endreq; stat = m_tesIO->store( myJets, m_outputLocation ) ; if(!stat){ log<<MSG::ERROR<<"Could not store jets"<<endreq; return stat; } MissingMomentum* mm = new MissingMomentum(missingMomentum); stat = m_tesIO->store(mm,m_missingMomentumLocation ) ; std::string message = stat? "Stored missing mom":"Failed to missing mom"; log<<MSG::DEBUG<<message<<endreq; return StatusCode::SUCCESS; }

StatusCode Atlfast::JetMaker::finalize (    )

Definition at line 227 of file JetMaker.cxx. { MsgStream log( messageService(), name() ) ; log << MSG::INFO << "finalizing" << endreq; return StatusCode::SUCCESS ; }

Jet * Atlfast::JetMaker::create (  Cluster *    )

Definition at line 389 of file JetMaker.cxx. References m_smearer, and Atlfast::ISmearer::smear(). Referenced by execute(). { Jet* jet; MsgStream log( messageService(), name() ) ; HepLorentzVector vec = src->momentum(); if (m_doSmearing) { // Ask our Smearer make a smeared 4-vector vec = m_smearer->smear(vec); } jet = new Jet(vec, *src); log << MSG::DEBUG << "Unsmeared Cluster: " << *src << endreq; log << MSG::DEBUG << "Smeared Jet : " << jet << endreq ; // return Jet return jet; }

bool Atlfast::JetMaker::isAcceptable (  Jet *    )

Definition at line 484 of file JetMaker.cxx. References m_maxEta, and m_minPT. Referenced by execute(). { MsgStream log( messageService(), name() ) ; // Apply kinematic criteria to the candidate jet if( candidate->pT() < m_minPT ) { log << MSG::DEBUG << "Candidate failed pt cut: pt " << candidate->pT() << " min was " << m_minPT << endreq; return false ; } if( abs(candidate->eta()) > m_maxEta ) { log << MSG::DEBUG << "Candidate failed max eta cut: eta " << candidate->eta() << " max was " << m_maxEta << endreq; return false ; } return true ; }

void Atlfast::JetMaker::addMuons (  Jet *   , MuonCollection &    )

Definition at line 416 of file JetMaker.cxx. References m_barrelForwardEta, m_kinehelp, m_rconeb, m_rconef, and MuonIterator. Referenced by execute(). { if (muons.size() == 0) return; MsgStream log( messageService(), name() ) ; MuonIterator first = muons.begin() ; MuonIterator last = muons.end() ; // Sort by distance from the jet sort( first, last, SortAttribute::DeltaR(jet) ); // check that closest muon is within cone // if it is, add muon to jet float muonJetRdist = m_kinehelp.deltaR( jet, *first ) ; log << MSG::DEBUG << "distance in R of muon to Jet " << muonJetRdist << endreq; if ( ( abs((*first)->eta()) <= m_barrelForwardEta && muonJetRdist < m_rconeb ) || ( abs((*first)->eta()) > m_barrelForwardEta && muonJetRdist < m_rconef ) ) { log << MSG::DEBUG << "Adding muon to jet with momentum " << (*first)->momentum() << endreq; log << MSG::DEBUG << "Old mom: " << *jet << endreq; jet->setMomentum(jet->momentum() + (*first)->momentum()); log << MSG::DEBUG << "New mom: " << *jet << endreq ; // create association IAssociationManager* ia = jet; ia->associate(*first); // remove muon from list muons.erase(first); } else { log << MSG::DEBUG << "Muon not in R cone" << endreq; } }

MsgStream& Atlfast::JetMaker::dumpParameters (  MsgStream &    )  const

---

Member Data Documentation

double Atlfast::JetMaker::m_minPT [private]

Definition at line 109 of file JetMaker.h. Referenced by initialize(), and isAcceptable().

double Atlfast::JetMaker::m_maxEta [private]

Definition at line 110 of file JetMaker.h. Referenced by initialize(), and isAcceptable().

bool Atlfast::JetMaker::m_doSmearing [private]

Definition at line 111 of file JetMaker.h. Referenced by initialize().

double Atlfast::JetMaker::m_barrelForwardEta [private]

Definition at line 113 of file JetMaker.h. Referenced by addMuons(), and initialize().

double Atlfast::JetMaker::m_rconeb [private]

Definition at line 114 of file JetMaker.h. Referenced by addMuons(), and initialize().

double Atlfast::JetMaker::m_rconef [private]

Definition at line 115 of file JetMaker.h. Referenced by addMuons(), and initialize().

double Atlfast::JetMaker::m_bPtMin [private]

Definition at line 117 of file JetMaker.h. Referenced by execute(), and initialize().

double Atlfast::JetMaker::m_bMaxDeltaR [private]

Definition at line 118 of file JetMaker.h. Referenced by initialize().

double Atlfast::JetMaker::m_cPtMin [private]

Definition at line 119 of file JetMaker.h. Referenced by execute(), and initialize().

double Atlfast::JetMaker::m_cMaxDeltaR [private]

Definition at line 120 of file JetMaker.h. Referenced by initialize().

double Atlfast::JetMaker::m_tauPtMin [private]

Definition at line 121 of file JetMaker.h. Referenced by execute(), and initialize().

double Atlfast::JetMaker::m_tauMaxDeltaR [private]

Definition at line 122 of file JetMaker.h. Referenced by initialize().

double Atlfast::JetMaker::m_etaTagMax [private]

Definition at line 123 of file JetMaker.h. Referenced by execute(), and initialize().

double Atlfast::JetMaker::m_tauJetPtRatio [private]

Definition at line 124 of file JetMaker.h. Referenced by initialize().

ISmearer* Atlfast::JetMaker::m_smearer [private]

Definition at line 132 of file JetMaker.h. Referenced by create(), initialize(), and ~JetMaker().

ISmearer* Atlfast::JetMaker::m_cellSmearer [private]

Definition at line 133 of file JetMaker.h. Referenced by addUnusedCells(), initialize(), and ~JetMaker().

TesIO* Atlfast::JetMaker::m_tesIO [private]

Definition at line 136 of file JetMaker.h. Referenced by addUnusedCells(), execute(), initialize(), and ~JetMaker().

std::string Atlfast::JetMaker::m_inputLocation [private]

Definition at line 144 of file JetMaker.h. Referenced by execute(), and initialize().

std::string Atlfast::JetMaker::m_outputLocation [private]

Definition at line 145 of file JetMaker.h. Referenced by execute(), and initialize().

std::string Atlfast::JetMaker::m_muonLocation [private]

Definition at line 146 of file JetMaker.h. Referenced by execute(), and initialize().

std::string Atlfast::JetMaker::m_mcTruthLocation [private]

Definition at line 147 of file JetMaker.h. Referenced by initialize().

std::string Atlfast::JetMaker::m_unusedCellLocation [private]

Definition at line 148 of file JetMaker.h. Referenced by addUnusedCells().

std::string Atlfast::JetMaker::m_missingMomentumLocation [private]

Definition at line 149 of file JetMaker.h. Referenced by execute(), and initialize().

KinematicHelper Atlfast::JetMaker::m_kinehelp [private]

Definition at line 155 of file JetMaker.h. Referenced by addMuons().

---

The documentation for this class was generated from the following files:

• JetMaker.h
• JetMaker.cxx

---