#include <MySub/LhcTrackAnalyzer/src/LhcTrackAnalyzer.cc>

Inheritance diagram for LhcTrackAnalyzer:

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

	~LhcTrackAnalyzer () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

void	endJob () override

void	SetVarToZero ()

Private Attributes
int	charge_ [nMaxtracks_]

double	chi2_ [nMaxtracks_]

double	chi2ndof_ [nMaxtracks_]

bool	debug_

double	dxy_ [nMaxtracks_]

double	dz_ [nMaxtracks_]

double	eta_ [nMaxtracks_]

int	event_

std::string	filename_

bool	goodbx_

bool	goodvtx_

int	isHighPurity_ [nMaxtracks_]

int	nTracks_

double	phi_ [nMaxtracks_]

double	pt_ [nMaxtracks_]

edm::InputTag	PVtxCollectionTag_

double	qoverp_ [nMaxtracks_]

TFile *	rootFile_

TTree *	rootTree_

int	run_

edm::InputTag	TrackCollectionTag_

int	trkAlgo_ [nMaxtracks_]

int	trkQuality_ [nMaxtracks_]

int	validhits_ [nMaxtracks_][7]

double	xPCA_ [nMaxtracks_]

double	yPCA_ [nMaxtracks_]

double	zPCA_ [nMaxtracks_]

Static Private Attributes
static const int	nMaxtracks_ = 3000

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

Originally written by M.Musich Expanded by A. Bonato

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 54 of file LhcTrackAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 61 of file LhcTrackAnalyzer.cc.

References debug_, filename_, edm::ParameterSet::getParameter(), PVtxCollectionTag_, AlCaHLTBitMon_QueryRunRegistry::string, and TrackCollectionTag_.

 {
   //now do what ever initialization is needed
   debug_ = iConfig.getParameter<bool>("Debug");
   TrackCollectionTag_ = iConfig.getParameter<edm::InputTag>("TrackCollectionTag");
   PVtxCollectionTag_ = iConfig.getParameter<edm::InputTag>("PVtxCollectionTag");
   filename_ = iConfig.getParameter<std::string>("OutputFileName");
 }

LhcTrackAnalyzer::~LhcTrackAnalyzer ( )

override

Definition at line 72 of file LhcTrackAnalyzer.cc.

                                     {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void LhcTrackAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 82 of file LhcTrackAnalyzer.cc.

References funct::abs(), reco::TrackBase::algoName(), charge_, chi2_, chi2ndof_, gather_cfg::cout, reco::TrackBase::ctf, ztail::d, debug_, reco::TrackBase::detachedQuadStep, reco::TrackBase::detachedTripletStep, dxy_, dz_, eta_, edm::EventID::event(), event_, Exception, edm::Event::getByLabel(), goodbx_, goodvtx_, reco::TrackBase::highPtTripletStep, reco::TrackBase::highPurity, edm::EventBase::id(), reco::TrackBase::initialStep, isHighPurity_, reco::TrackBase::jetCoreRegionalStep, reco::TrackBase::loose, reco::TrackBase::lowPtQuadStep, reco::TrackBase::lowPtTripletStep, reco::TrackBase::mixedTripletStep, nMaxtracks_, nTracks_, phi_, reco::TrackBase::pixelLessStep, reco::TrackBase::pixelPairStep, pt_, PVtxCollectionTag_, qoverp_, rootTree_, edm::EventID::run(), run_, SetVarToZero(), reco::TrackBase::tight, reco::TrackBase::tobTecStep, HLT_2018_cff::track, TrackCollectionTag_, trkAlgo_, trkQuality_, reco::TrackBase::undefAlgorithm, reco::TrackBase::undefQuality, validhits_, badGlobalMuonTaggersAOD_cff::vtx, xPCA_, yPCA_, and zPCA_.

                                                                                   {
   using namespace edm;
   using namespace reco;
   using namespace std;
 
   //=======================================================
   // Initialize Root-tuple variables
   //=======================================================
 
   SetVarToZero();
 
   //=======================================================
   // Retrieve the Track information
   //=======================================================
 
   Handle<TrackCollection> trackCollectionHandle;
   iEvent.getByLabel(TrackCollectionTag_, trackCollectionHandle);
   Handle<VertexCollection> vertexCollectionHandle;
   iEvent.getByLabel(PVtxCollectionTag_, vertexCollectionHandle);
   for (VertexCollection::const_iterator vtx = vertexCollectionHandle->begin(); vtx != vertexCollectionHandle->end();
        ++vtx) {
     if (vtx == vertexCollectionHandle->begin()) {
       if (vtx->isFake())
         goodvtx_ = false;
       else
         goodvtx_ = true;
     } else
       break;
   }
 
   goodbx_ = true;
   // int bx = iEvent.bunchCrossing();
   //if (bx==51 || bx==2724) goodbx_=true;
 
   run_ = iEvent.id().run();
   event_ = iEvent.id().event();
 
   if (debug_)
     cout << "LhcTrackAnalyzer::analyze() looping over " << trackCollectionHandle->size() << "tracks." << endl;
 
   // unsigned int i = 0;
   for (TrackCollection::const_iterator track = trackCollectionHandle->begin(); track != trackCollectionHandle->end();
        ++track) {
     if (nTracks_ >= nMaxtracks_) {
       std::cout << " LhcTrackAnalyzer::analyze() : Warning - Run " << run_ << " Event " << event_
                 << "\tNumber of tracks: " << trackCollectionHandle->size() << " , greater than " << nMaxtracks_
                 << std::endl;
       continue;
     }
     pt_[nTracks_] = track->pt();
     eta_[nTracks_] = track->eta();
     phi_[nTracks_] = track->phi();
     chi2_[nTracks_] = track->chi2();
     chi2ndof_[nTracks_] = track->normalizedChi2();
     charge_[nTracks_] = track->charge();
     qoverp_[nTracks_] = track->qoverp();
     dz_[nTracks_] = track->dz();
     dxy_[nTracks_] = track->dxy();
     xPCA_[nTracks_] = track->vertex().x();
     yPCA_[nTracks_] = track->vertex().y();
     zPCA_[nTracks_] = track->vertex().z();
     validhits_[nTracks_][0] = track->numberOfValidHits();
     validhits_[nTracks_][1] = track->hitPattern().numberOfValidPixelBarrelHits();
     validhits_[nTracks_][2] = track->hitPattern().numberOfValidPixelEndcapHits();
     validhits_[nTracks_][3] = track->hitPattern().numberOfValidStripTIBHits();
     validhits_[nTracks_][4] = track->hitPattern().numberOfValidStripTIDHits();
     validhits_[nTracks_][5] = track->hitPattern().numberOfValidStripTOBHits();
     validhits_[nTracks_][6] = track->hitPattern().numberOfValidStripTECHits();
 
     int myalgo = -88;
     if (track->algo() == reco::TrackBase::undefAlgorithm)
       myalgo = 0;
     if (track->algo() == reco::TrackBase::ctf)
       myalgo = 1;
     if (track->algo() == reco::TrackBase::initialStep)
       myalgo = 4;
     if (track->algo() == reco::TrackBase::lowPtTripletStep)
       myalgo = 5;
     if (track->algo() == reco::TrackBase::pixelPairStep)
       myalgo = 6;
     if (track->algo() == reco::TrackBase::detachedTripletStep)
       myalgo = 7;
     if (track->algo() == reco::TrackBase::mixedTripletStep)
       myalgo = 8;
     if (track->algo() == reco::TrackBase::pixelLessStep)
       myalgo = 9;
     if (track->algo() == reco::TrackBase::tobTecStep)
       myalgo = 10;
     if (track->algo() == reco::TrackBase::jetCoreRegionalStep)
       myalgo = 11;
     // This class is pending the migration to Phase1 tracks
     if (track->algo() == reco::TrackBase::highPtTripletStep || track->algo() == reco::TrackBase::lowPtQuadStep ||
         track->algo() == reco::TrackBase::detachedQuadStep) {
       throw cms::Exception("Not implemented")
           << "LhcTrackAnalyzer does not yet support phase1 tracks, encountered one from algo "
           << reco::TrackBase::algoName(track->algo());
     }
     trkAlgo_[nTracks_] = myalgo;
 
     int myquality = -99;
     if (track->quality(reco::TrackBase::undefQuality))
       myquality = -1;
     if (track->quality(reco::TrackBase::loose))
       myquality = 0;
     if (track->quality(reco::TrackBase::tight))
       myquality = 1;
     if (track->quality(reco::TrackBase::highPurity))
       myquality = 2;
     //if(track->quality(reco::TrackBase::confirmed))myquality=3;
     // if(track->quality(reco::TrackBase::goodIterative))myquality=4;
     // if(track->quality(reco::TrackBase::qualitySize))myquality=5;
     trkQuality_[nTracks_] = myquality;
 
     if (track->quality(reco::TrackBase::highPurity))
       isHighPurity_[nTracks_] = 1;
     else
       isHighPurity_[nTracks_] = 0;
     nTracks_++;
 
   }  //end loop on tracks
 
   for (int d = 0; d < nTracks_; ++d) {
     if (abs(trkQuality_[d]) > 5)
       cout << "MYQUALITY!!! " << trkQuality_[d] << " at track # " << d << "/" << nTracks_ << endl;
   }
 
   rootTree_->Fill();
 }

void LhcTrackAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 212 of file LhcTrackAnalyzer.cc.

References charge_, chi2_, chi2ndof_, dxy_, dz_, eta_, event_, filename_, goodbx_, goodvtx_, isHighPurity_, nTracks_, phi_, pt_, qoverp_, rootFile_, rootTree_, run_, trkAlgo_, trkQuality_, validhits_, xPCA_, yPCA_, and zPCA_.

                                 {
   edm::LogInfo("beginJob") << "Begin Job" << std::endl;
   // Define TTree for output
   rootFile_ = new TFile(filename_.c_str(), "recreate");
   rootTree_ = new TTree("tree", "Lhc Track tree");
 
   // Track Paramters
   rootTree_->Branch("run", &run_, "run/I");
   rootTree_->Branch("event", &event_, "event/I");
   rootTree_->Branch("goodbx", &goodbx_, "goodbx/O");
   rootTree_->Branch("goodvtx", &goodvtx_, "goodvtx/O");
   rootTree_->Branch("nTracks", &nTracks_, "nTracks/I");
   rootTree_->Branch("pt", &pt_, "pt[nTracks]/D");
   rootTree_->Branch("eta", &eta_, "eta[nTracks]/D");
   rootTree_->Branch("phi", &phi_, "phi[nTracks]/D");
   rootTree_->Branch("chi2", &chi2_, "chi2[nTracks]/D");
   rootTree_->Branch("chi2ndof", &chi2ndof_, "chi2ndof[nTracks]/D");
   rootTree_->Branch("charge", &charge_, "charge[nTracks]/I");
   rootTree_->Branch("qoverp", &qoverp_, "qoverp[nTracks]/D");
   rootTree_->Branch("dz", &dz_, "dz[nTracks]/D");
   rootTree_->Branch("dxy", &dxy_, "dxy[nTracks]/D");
   rootTree_->Branch("xPCA", &xPCA_, "xPCA[nTracks]/D");
   rootTree_->Branch("yPCA", &yPCA_, "yPCA[nTracks]/D");
   rootTree_->Branch("zPCA", &zPCA_, "zPCA[nTracks]/D");
   rootTree_->Branch("isHighPurity", &isHighPurity_, "isHighPurity[nTracks]/I");
   rootTree_->Branch("trkQuality", &trkQuality_, "trkQuality[nTracks]/I");
   rootTree_->Branch("trkAlgo", &trkAlgo_, "trkAlgo[nTracks]/I");
   rootTree_->Branch("nValidHits", &validhits_, "nValidHits[nTracks][7]/I");
 }

void LhcTrackAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 243 of file LhcTrackAnalyzer.cc.

References rootFile_.

Referenced by o2olib.O2ORunMgr::executeJob().

                               {
   if (rootFile_) {
     rootFile_->Write();
     rootFile_->Close();
   }
 }

void LhcTrackAnalyzer::SetVarToZero ( )

private

Definition at line 250 of file LhcTrackAnalyzer.cc.

References charge_, chi2_, chi2ndof_, DEFINE_FWK_MODULE, dxy_, dz_, eta_, event_, mps_fire::i, isHighPurity_, dqmiolumiharvest::j, nMaxtracks_, nTracks_, phi_, pt_, qoverp_, run_, trkAlgo_, trkQuality_, validhits_, xPCA_, yPCA_, and zPCA_.

Referenced by analyze().

                                     {
   run_ = -1;
   event_ = -99;
   nTracks_ = 0;
   for (int i = 0; i < nMaxtracks_; ++i) {
     pt_[i] = 0;
     eta_[i] = 0;
     phi_[i] = 0;
     chi2_[i] = 0;
     chi2ndof_[i] = 0;
     charge_[i] = 0;
     qoverp_[i] = 0;
     dz_[i] = 0;
     dxy_[i] = 0;
     xPCA_[i] = 0;
     yPCA_[i] = 0;
     zPCA_[i] = 0;
     trkQuality_[i] = 0;
     trkAlgo_[i] = -1;
     isHighPurity_[i] = -3;
     for (int j = 0; j < 7; j++) {
       validhits_[nTracks_][j] = -1 * j;
     }
   }
 }