LhcTrackAnalyzer Class Reference

#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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (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::LhcTrackAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 61 of file LhcTrackAnalyzer.cc.

{
  //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");
}

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

~LhcTrackAnalyzer()

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

analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 82 of file LhcTrackAnalyzer.cc.

                                                                                  {
  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();
}

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_, event_, Exception, goodbx_, goodvtx_, reco::TrackBase::highPtTripletStep, reco::TrackBase::highPurity, iEvent, 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_, 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_.

beginJob()

void LhcTrackAnalyzer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 212 of file LhcTrackAnalyzer.cc.

                                {
  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");
}

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_.

endJob()

void LhcTrackAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 243 of file LhcTrackAnalyzer.cc.

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

References rootFile_.

Referenced by o2olib.O2ORunMgr::executeJob().

SetVarToZero()

void LhcTrackAnalyzer::SetVarToZero ( )

private

Definition at line 250 of file LhcTrackAnalyzer.cc.

                                    {
  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;
    }
  }
}

References charge_, chi2_, chi2ndof_, 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().

Member Data Documentation

charge_

int LhcTrackAnalyzer::charge_[nMaxtracks_]

private

Definition at line 87 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

chi2_

double LhcTrackAnalyzer::chi2_[nMaxtracks_]

private

Definition at line 85 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

chi2ndof_

double LhcTrackAnalyzer::chi2ndof_[nMaxtracks_]

private

Definition at line 86 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

debug_

bool LhcTrackAnalyzer::debug_

private

Definition at line 67 of file LhcTrackAnalyzer.h.

Referenced by analyze(), and LhcTrackAnalyzer().

dxy_

double LhcTrackAnalyzer::dxy_[nMaxtracks_]

private

Definition at line 90 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

dz_

double LhcTrackAnalyzer::dz_[nMaxtracks_]

private

Definition at line 89 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

eta_

double LhcTrackAnalyzer::eta_[nMaxtracks_]

private

Definition at line 83 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

event_

int LhcTrackAnalyzer::event_

private

Definition at line 81 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

filename_

std::string LhcTrackAnalyzer::filename_

private

Definition at line 70 of file LhcTrackAnalyzer.h.

Referenced by beginJob(), and LhcTrackAnalyzer().

goodbx_

bool LhcTrackAnalyzer::goodbx_

private

Definition at line 98 of file LhcTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

goodvtx_

bool LhcTrackAnalyzer::goodvtx_

private

Definition at line 99 of file LhcTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

isHighPurity_

int LhcTrackAnalyzer::isHighPurity_[nMaxtracks_]

private

Definition at line 96 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

nMaxtracks_

const int LhcTrackAnalyzer::nMaxtracks_ = 3000

staticprivate

Definition at line 78 of file LhcTrackAnalyzer.h.

Referenced by analyze(), and SetVarToZero().

nTracks_

int LhcTrackAnalyzer::nTracks_

private

Definition at line 79 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

phi_

double LhcTrackAnalyzer::phi_[nMaxtracks_]

private

Definition at line 84 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

pt_

double LhcTrackAnalyzer::pt_[nMaxtracks_]

private

Definition at line 82 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

PVtxCollectionTag_

edm::InputTag LhcTrackAnalyzer::PVtxCollectionTag_

private

Definition at line 66 of file LhcTrackAnalyzer.h.

Referenced by analyze(), and LhcTrackAnalyzer().

qoverp_

double LhcTrackAnalyzer::qoverp_[nMaxtracks_]

private

Definition at line 88 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

rootFile_

TFile* LhcTrackAnalyzer::rootFile_

private

Definition at line 71 of file LhcTrackAnalyzer.h.

Referenced by beginJob(), and endJob().

rootTree_

TTree* LhcTrackAnalyzer::rootTree_

private

Definition at line 72 of file LhcTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

run_

int LhcTrackAnalyzer::run_

private

Definition at line 80 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

TrackCollectionTag_

edm::InputTag LhcTrackAnalyzer::TrackCollectionTag_

private

Definition at line 65 of file LhcTrackAnalyzer.h.

Referenced by analyze(), and LhcTrackAnalyzer().

trkAlgo_

int LhcTrackAnalyzer::trkAlgo_[nMaxtracks_]

private

Definition at line 94 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

trkQuality_

int LhcTrackAnalyzer::trkQuality_[nMaxtracks_]

private

Definition at line 95 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

validhits_

int LhcTrackAnalyzer::validhits_[nMaxtracks_][7]

private

Definition at line 97 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

xPCA_

double LhcTrackAnalyzer::xPCA_[nMaxtracks_]

private

Definition at line 91 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

yPCA_

double LhcTrackAnalyzer::yPCA_[nMaxtracks_]

private

Definition at line 92 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

zPCA_

double LhcTrackAnalyzer::zPCA_[nMaxtracks_]

private

Definition at line 93 of file LhcTrackAnalyzer.h.

Referenced by analyze(), beginJob(), and SetVarToZero().