GEMTrackMatch Class Reference

#include <GEMTrackMatch.h>

Inheritance diagram for GEMTrackMatch:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &) override=0
void	buildLUT (const int maxChamberId)
void	FillWithTrigger (MonitorElement *me[3], Float_t eta)
void	FillWithTrigger (MonitorElement *me[3][3], Float_t eta, Float_t phi, bool odd[3], bool even[3])
void	FillWithTrigger (MonitorElement *me[4][3], bool array[3][2], Float_t value)
void	FillWithTrigger (MonitorElement *me[4][3][3], bool array[3][2], Float_t eta, Float_t phi, bool odd[3], bool even[3])
	GEMTrackMatch (const edm::ParameterSet &cfg)
std::pair< int, int >	getClosestChambers (const int maxChamberId, int region, float phi)
std::pair< double, double >	getEtaRange (int station, int chamber)
std::pair< double, double >	getEtaRangeForPhi (int station)
bool	isSimTrackGood (const SimTrack &)
void	setGeometry (const GEMGeometry &geom)
	~GEMTrackMatch () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Protected Attributes
edm::ParameterSet	cfg_
float	chamberHeight_
bool	detailPlot_
std::vector< double >	etaRangeForPhi
const GEMGeometry *	gem_geom_
float	maxEta_
float	minEta_
float	minPt_
std::pair< std::vector< float >, std::vector< int > >	negativeLUT_
unsigned int	nstation
std::pair< std::vector< float >, std::vector< int > >	positiveLUT_
float	radiusCenter_
edm::EDGetToken	simHitsToken_
edm::EDGetToken	simTracksToken_
edm::EDGetToken	simVerticesToken_
int	useRoll_

Detailed Description

Definition at line 42 of file GEMTrackMatch.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file GEMTrackMatch.cc.

References edm::ParameterSet::getUntrackedParameter(), maxEta_, minEta_, and minPt_.

{
  minPt_  = ps.getUntrackedParameter<double>("gemMinPt",5.0);
  minEta_ = ps.getUntrackedParameter<double>("gemMinEta",1.55);
  maxEta_ = ps.getUntrackedParameter<double>("gemMaxEta",2.45);
}

GEMTrackMatch::~GEMTrackMatch ( )

override

Definition at line 14 of file GEMTrackMatch.cc.

                              {
}

Member Function Documentation

void GEMTrackMatch::analyze ( const edm::Event &  e, const edm::EventSetup &    )

overridepure virtual

Implemented in GEMRecHitTrackMatch, GEMDigiTrackMatch, and GEMSimTrackMatch.

void GEMTrackMatch::buildLUT

(

const int

maxChamberId

)

Definition at line 140 of file GEMTrackMatch.cc.

References GEMGeometry::chambers(), gem_geom_, negativeLUT_, and positiveLUT_.

Referenced by setGeometry().

{
  edm::LogInfo("GEMTrackMatch")<<"GEMTrackMatch::buildLUT"<<std::endl;
  edm::LogInfo("GEMTrackMatch")<<"max chamber "<<maxChamberId<<"\n";

  std::vector<int> pos_ids, neg_ids;
  std::vector<float> phis_pos;
  std::vector<float> phis_neg;
  LocalPoint lCentre( 0., 0., 0. );

  for(auto it : gem_geom_->chambers()) {
    if(it->id().region()>0) {
      pos_ids.push_back(it->id().rawId());
      edm::LogInfo("GEMTrackMatch")<<"added id = "<<it->id()<<" = "<<it->id().rawId()<<" to pos ids"<<std::endl;
      const BoundPlane& bSurface(it->surface());
      GlobalPoint gCentre(bSurface.toGlobal(lCentre));
      int cphi(gCentre.phi().degrees());
      if (cphi < 0) cphi += 360;
      phis_pos.push_back(cphi);
      edm::LogInfo("GEMTrackMatch")<<"added phi = "<<cphi<<" to phi pos vector"<<std::endl;
    }
    else if(it->id().region()<0) {
      neg_ids.push_back(it->id().rawId());
      edm::LogInfo("GEMTrackMatch")<<"added id = "<<it->id()<<" = "<<it->id().rawId()<<" to neg ids"<<std::endl;
      const BoundPlane& bSurface(it->surface());
      GlobalPoint gCentre(bSurface.toGlobal(lCentre));
      int cphi(gCentre.phi().degrees());
      if (cphi < 0) cphi += 360;
      phis_neg.push_back(cphi);
      edm::LogInfo("GEMTrackMatch")<<"added phi = "<<cphi<<" to phi neg vector"<<std::endl;
    }
  }
  positiveLUT_ = std::make_pair(phis_pos,pos_ids);
  negativeLUT_ = std::make_pair(phis_neg,neg_ids);
}

void GEMTrackMatch::FillWithTrigger	(	MonitorElement *	me[3],
		Float_t	eta
	)

Definition at line 17 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), mps_fire::i, and nstation.

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), and GEMRecHitTrackMatch::analyze().

{
  for( unsigned int i=0 ; i<nstation ; i++) {
    hist[i]->Fill(eta);
  }
  return;
}

void GEMTrackMatch::FillWithTrigger	(	MonitorElement *	me[3][3],
		Float_t	eta,
		Float_t	phi,
		bool	odd[3],
		bool	even[3]
	)

Definition at line 25 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), plotBeamSpotDB::first, getEtaRange(), mps_fire::i, nstation, edm::second(), and relativeConstraints::station.

{
  for( unsigned int i=0 ; i<nstation ; i++) {
    int station = i+1;
    if ( odd[i] && eta > getEtaRange(station,1).first&& eta < getEtaRange(station,1).second ) {
      hist[i][0]->Fill(phi);
      if( hist[i][1] != nullptr) hist[i][1]->Fill(phi);
    }
    if ( even[i] && eta > getEtaRange(station,2).first&& eta < getEtaRange(station,2).second ) {
      hist[i][0]->Fill(phi);
      if( hist[i][2] != nullptr) hist[i][2]->Fill(phi);
    }
  }
  return;
}

void GEMTrackMatch::FillWithTrigger	(	MonitorElement *	me[4][3],
		bool	array[3][2],
		Float_t	value
	)

Definition at line 41 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), mps_fire::i, and nstation.

{
  for( unsigned int i=0 ; i<nstation ; i++) {
    if ( array[i][0] ) hist[0][i]->Fill(value);
    if ( array[i][1] ) hist[1][i]->Fill(value);
    if ( array[i][0] || array[i][1] ) hist[2][i]->Fill(value);
    if ( array[i][0] && array[i][1] ) hist[3][i]->Fill(value);
  } 
  return;
}

void GEMTrackMatch::FillWithTrigger	(	MonitorElement *	me[4][3][3],
		bool	array[3][2],
		Float_t	eta,
		Float_t	phi,
		bool	odd[3],
		bool	even[3]
	)

Definition at line 52 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), plotBeamSpotDB::first, getEtaRange(), mps_fire::i, nstation, edm::second(), and relativeConstraints::station.

{
  for( unsigned int i=0 ; i<nstation ; i++) {
    int station = i+1;
    if ( odd[i] && eta > getEtaRange(station,1).first&& eta < getEtaRange(station,1).second ) {
      if ( array[i][0] ) { 
        hist[0][i][0]->Fill(phi);
        if ( hist[0][i][1] != nullptr) hist[0][i][1]->Fill(phi); 
      }
      if ( array[i][1] ) { 
        hist[1][i][0]->Fill(phi); 
        if( hist[1][i][1] != nullptr ) hist[1][i][1]->Fill(phi); 
      }
      if ( array[i][0] || array[i][1] ) { 
        hist[2][i][0]->Fill(phi); 
        if ( hist[2][i][1] != nullptr) hist[2][i][1]->Fill(phi); 
      }
      if ( array[i][0] && array[i][1] ) { 
        hist[3][i][0]->Fill(phi); 
        if ( hist[3][i][1] != nullptr) hist[3][i][1]->Fill(phi); 
      }
    }
    if ( even[i] && eta > getEtaRange(station,2).first&& eta < getEtaRange(station,2).second ) {
      if ( array[i][0] ) { 
        hist[0][i][0]->Fill(phi); 
        if ( hist[0][i][2]!=nullptr) hist[0][i][2]->Fill(phi); 
      }
      if ( array[i][1] ) { 
        hist[1][i][0]->Fill(phi); 
        if( hist[1][i][2]!=nullptr ) hist[1][i][2]->Fill(phi); 
      }
      if ( array[i][0] || array[i][1] ) { 
        hist[2][i][0]->Fill(phi); 
        if( hist[2][i][2]!=nullptr) hist[2][i][2]->Fill(phi); 
      }
      if ( array[i][0] && array[i][1] ) { 
        hist[3][i][0]->Fill(phi); 
        if( hist[3][i][2]!=nullptr) hist[3][i][2]->Fill(phi);
      }
    }
  }
  return;
}

std::pair< int, int > GEMTrackMatch::getClosestChambers	(	const int	maxChamberId,
		int	region,
		float	phi
	)

Definition at line 202 of file GEMTrackMatch.cc.

References negativeLUT_, phi, positiveLUT_, and pileupCalc::upper.

{
  auto& phis(positiveLUT_.first);
  auto upper = std::upper_bound(phis.begin(), phis.end(), phi);
  auto& LUT = (region == 1 ? positiveLUT_.second : negativeLUT_.second);
  return std::make_pair(LUT.at(upper - phis.begin()), (LUT.at((upper - phis.begin() + 1)%maxChamberId)));
}

std::pair< double, double > GEMTrackMatch::getEtaRange	(	int	station,
		int	chamber
	)

Definition at line 96 of file GEMTrackMatch.cc.

References GEMGeometry::chambers(), gather_cfg::cout, and gem_geom_.

Referenced by FillWithTrigger().

{
  if( gem_geom_ != nullptr) {
    auto& ch = gem_geom_->chambers().front();
    auto& roll1 = ch->etaPartitions().front(); //.begin();
    auto& roll2 = ch->etaPartitions()[ch->nEtaPartitions()-1];
    const BoundPlane& bSurface1(roll1->surface());
    const BoundPlane& bSurface2(roll2->surface());
    auto& parameters1( roll1->specs()->parameters());
    float height1(parameters1[2]);
    auto& parameters2( roll2->specs()->parameters());
    float height2(parameters2[2]);
    LocalPoint lTop1( 0., height1, 0.);
    GlobalPoint gTop1(bSurface1.toGlobal(lTop1));
    //LocalPoint lBottom1( 0., -height1, 0.);
    //GlobalPoint gBottom1(bSurface1.toGlobal(lBottom1));
    //LocalPoint lTop2( 0., height2, 0.);
    //GlobalPoint gTop2(bSurface2.toGlobal(lTop2));
    LocalPoint lBottom2( 0., -height2, 0.);
    GlobalPoint gBottom2(bSurface2.toGlobal(lBottom2));
    double eta1 = fabs(gTop1.eta()) - 0.01;
    double eta2 = fabs(gBottom2.eta()) + 0.01;
    return std::make_pair(eta1,eta2);
  }
  else { std::cout<<"Failed to get geometry information"<<std::endl;
    return std::make_pair(0,0);
  }
}

std::pair< double, double > GEMTrackMatch::getEtaRangeForPhi

(

int

station

)

Definition at line 210 of file GEMTrackMatch.cc.

References etaRangeForPhi.

{
  std::pair<double, double> range;
  if( station== 0 )      range = std::make_pair( etaRangeForPhi[0],etaRangeForPhi[1]) ; 
  else if( station== 1 ) range = std::make_pair( etaRangeForPhi[2],etaRangeForPhi[3]) ; 
  else if( station== 2 ) range = std::make_pair( etaRangeForPhi[4],etaRangeForPhi[5]) ; 

  return range;
}

bool GEMTrackMatch::isSimTrackGood

(

const SimTrack &

t

)

Definition at line 126 of file GEMTrackMatch.cc.

References funct::abs(), PVValHelper::eta, maxEta_, minEta_, minPt_, CoreSimTrack::momentum(), SimTrack::noGenpart(), SimTrack::noVertex(), and CoreSimTrack::type().

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), and GEMRecHitTrackMatch::analyze().

{

  // SimTrack selection
  if (t.noVertex())   return false; 
  if (t.noGenpart()) return false;
  if (std::abs(t.type()) != 13) return false; // only interested in direct muon simtracks
  if (t.momentum().pt() < minPt_ ) return false;
  const float eta(std::abs(t.momentum().eta()));
  if (eta > maxEta_ || eta < minEta_ ) return false; // no GEMs could be in such eta
  return true;
}

void GEMTrackMatch::setGeometry

(

const GEMGeometry &

geom

)

Definition at line 177 of file GEMTrackMatch.cc.

References buildLUT(), chamberHeight_, GEMGeometry::chambers(), gem_geom_, relativeConstraints::geom, GEMGeometry::idToDetUnit(), PV3DBase< T, PVType, FrameType >::perp(), radiusCenter_, GEMGeometry::regions(), GEMEtaPartitionSpecs::specificTopology(), GEMEtaPartition::specs(), StripTopology::stripLength(), and useRoll_.

Referenced by GEMSimTrackMatch::bookHistograms(), GEMDigiTrackMatch::bookHistograms(), and GEMRecHitTrackMatch::bookHistograms().

{
  edm::LogInfo("GEMTrackMatch")<<"GEMTrackMatch :: setGeometry"<<std::endl;
  gem_geom_ = &geom;
  GEMDetId chId(gem_geom_->chambers().front()->id());
  useRoll_ = 1;
  const auto top_chamber = static_cast<const GEMEtaPartition*>(geom.idToDetUnit(GEMDetId(chId.region(),chId.ring(),chId.station(),chId.layer(),chId.chamber(),useRoll_)));
  const int nEtaPartitions(gem_geom_->chambers().front()->nEtaPartitions());
  const auto bottom_chamber = static_cast<const GEMEtaPartition*>(geom.idToDetUnit(GEMDetId(chId.region(),chId.ring(),chId.station(),chId.layer(),chId.chamber(),nEtaPartitions)));
  const float top_half_striplength = top_chamber->specs()->specificTopology().stripLength()/2.;
  const float bottom_half_striplength = bottom_chamber->specs()->specificTopology().stripLength()/2.;
  const LocalPoint lp_top(0., top_half_striplength, 0.);
  const LocalPoint lp_bottom(0., -bottom_half_striplength, 0.);
  const GlobalPoint gp_top = top_chamber->toGlobal(lp_top);
  const GlobalPoint gp_bottom = bottom_chamber->toGlobal(lp_bottom);

  radiusCenter_ = (gp_bottom.perp() + gp_top.perp())/2.;
  chamberHeight_ = gp_top.perp() - gp_bottom.perp();

  const int maxChamberId = geom.regions()[0]->stations()[0]->superChambers().size();
  edm::LogInfo("GEMTrackMatch")<<"GEMTrackMatch :: setGeometry --> Calling buildLUT"<<std::endl;
  buildLUT(maxChamberId);
}  

Member Data Documentation

edm::ParameterSet GEMTrackMatch::cfg_

protected

Definition at line 63 of file GEMTrackMatch.h.

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), GEMRecHitTrackMatch::analyze(), GEMDigiTrackMatch::GEMDigiTrackMatch(), GEMRecHitTrackMatch::GEMRecHitTrackMatch(), and GEMSimTrackMatch::GEMSimTrackMatch().

float GEMTrackMatch::chamberHeight_

protected

Definition at line 77 of file GEMTrackMatch.h.

Referenced by setGeometry().

bool GEMTrackMatch::detailPlot_

protected

Definition at line 81 of file GEMTrackMatch.h.

Referenced by GEMDigiTrackMatch::analyze(), GEMDigiTrackMatch::bookHistograms(), GEMDigiTrackMatch::GEMDigiTrackMatch(), and GEMRecHitTrackMatch::GEMRecHitTrackMatch().

std::vector< double > GEMTrackMatch::etaRangeForPhi

protected

Definition at line 71 of file GEMTrackMatch.h.

Referenced by getEtaRangeForPhi().

const GEMGeometry* GEMTrackMatch::gem_geom_

protected

Definition at line 79 of file GEMTrackMatch.h.

Referenced by buildLUT(), getEtaRange(), and setGeometry().

float GEMTrackMatch::maxEta_

protected

Definition at line 76 of file GEMTrackMatch.h.

Referenced by GEMSimTrackMatch::bookHistograms(), GEMDigiTrackMatch::bookHistograms(), GEMRecHitTrackMatch::bookHistograms(), GEMTrackMatch(), and isSimTrackGood().

float GEMTrackMatch::minEta_

protected

Definition at line 75 of file GEMTrackMatch.h.

Referenced by GEMSimTrackMatch::bookHistograms(), GEMDigiTrackMatch::bookHistograms(), GEMRecHitTrackMatch::bookHistograms(), GEMTrackMatch(), and isSimTrackGood().

float GEMTrackMatch::minPt_

protected

Definition at line 74 of file GEMTrackMatch.h.

Referenced by GEMTrackMatch(), and isSimTrackGood().

std::pair<std::vector<float>,std::vector<int> > GEMTrackMatch::negativeLUT_

protected

Definition at line 69 of file GEMTrackMatch.h.

Referenced by buildLUT(), and getClosestChambers().

unsigned int GEMTrackMatch::nstation

protected

Definition at line 80 of file GEMTrackMatch.h.

Referenced by GEMSimTrackMatch::bookHistograms(), GEMDigiTrackMatch::bookHistograms(), GEMRecHitTrackMatch::bookHistograms(), and FillWithTrigger().

std::pair<std::vector<float>,std::vector<int> > GEMTrackMatch::positiveLUT_

protected

Definition at line 68 of file GEMTrackMatch.h.

Referenced by buildLUT(), and getClosestChambers().

float GEMTrackMatch::radiusCenter_

protected

Definition at line 77 of file GEMTrackMatch.h.

Referenced by setGeometry().

edm::EDGetToken GEMTrackMatch::simHitsToken_

protected

Definition at line 64 of file GEMTrackMatch.h.

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), GEMRecHitTrackMatch::analyze(), GEMDigiTrackMatch::GEMDigiTrackMatch(), GEMRecHitTrackMatch::GEMRecHitTrackMatch(), and GEMSimTrackMatch::GEMSimTrackMatch().

edm::EDGetToken GEMTrackMatch::simTracksToken_

protected

Definition at line 65 of file GEMTrackMatch.h.

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), GEMRecHitTrackMatch::analyze(), GEMDigiTrackMatch::GEMDigiTrackMatch(), GEMRecHitTrackMatch::GEMRecHitTrackMatch(), and GEMSimTrackMatch::GEMSimTrackMatch().

edm::EDGetToken GEMTrackMatch::simVerticesToken_

protected

Definition at line 66 of file GEMTrackMatch.h.

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), GEMRecHitTrackMatch::analyze(), GEMDigiTrackMatch::GEMDigiTrackMatch(), GEMRecHitTrackMatch::GEMRecHitTrackMatch(), and GEMSimTrackMatch::GEMSimTrackMatch().

int GEMTrackMatch::useRoll_

protected

Definition at line 78 of file GEMTrackMatch.h.

Referenced by setGeometry().