RPCSeedrecHitFinder Class Reference

#include <RPCSeedrecHitFinder.h>

Public Member Functions
void	configure (const edm::ParameterSet &iConfig)
void	fillrecHits ()
	RPCSeedrecHitFinder ()
void	setInput (MuonRecHitContainer(&recHits)[12])
void	setLayers (const std::vector< unsigned int > &Layers)
void	setOutput (RPCSeedFinder *Seed)
void	unsetInput ()
	~RPCSeedrecHitFinder ()

Private Types
typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer	ConstMuonRecHitContainer
typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer	ConstMuonRecHitPointer
typedef MuonTransientTrackingRecHit::MuonRecHitContainer	MuonRecHitContainer
typedef MuonTransientTrackingRecHit::MuonRecHitPointer	MuonRecHitPointer

Private Member Functions
void	checkandfill ()
void	complete (unsigned int LayerIndex)
double	getdeltaPhifromrecHits ()

Private Attributes
unsigned int	BxRange
std::vector< int >	ClusterSet
bool	isConfigured
bool	isInputset
bool	isLayerset
bool	isOutputset
std::vector< unsigned int >	LayersinRPC
double	MaxDeltaPhi
MuonRecHitContainer *	recHitsRPC [12]
ConstMuonRecHitContainer	therecHits
RPCSeedFinder *	theSeed

Detailed Description

Definition at line 27 of file RPCSeedrecHitFinder.h.

Member Typedef Documentation

typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer RPCSeedrecHitFinder::ConstMuonRecHitContainer

private

Definition at line 31 of file RPCSeedrecHitFinder.h.

typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer RPCSeedrecHitFinder::ConstMuonRecHitPointer

private

Definition at line 29 of file RPCSeedrecHitFinder.h.

typedef MuonTransientTrackingRecHit::MuonRecHitContainer RPCSeedrecHitFinder::MuonRecHitContainer

private

Definition at line 30 of file RPCSeedrecHitFinder.h.

typedef MuonTransientTrackingRecHit::MuonRecHitPointer RPCSeedrecHitFinder::MuonRecHitPointer

private

Definition at line 28 of file RPCSeedrecHitFinder.h.

Constructor & Destructor Documentation

RPCSeedrecHitFinder::RPCSeedrecHitFinder

(

)

Definition at line 21 of file RPCSeedrecHitFinder.cc.

References RPCSeedFinder::isConfigured, and RPCSeedFinder::isOutputset.

                                         {
  // Initiate the member
  isLayerset = false;
  isConfigured = false;
  isInputset = false;
  isOutputset = false;
  BxRange = 0;
  MaxDeltaPhi = 0;
  ClusterSet.clear();
  LayersinRPC.clear();
  therecHits.clear();
}

RPCSeedrecHitFinder::~RPCSeedrecHitFinder

(

)

Definition at line 34 of file RPCSeedrecHitFinder.cc.

{}

Member Function Documentation

void RPCSeedrecHitFinder::checkandfill ( )

private

Definition at line 220 of file RPCSeedrecHitFinder.cc.

References gather_cfg::cout.

                                       {
  if (therecHits.size() >= 3) {
    theSeed->setrecHits(therecHits);
    theSeed->seed();
  } else
    cout << "Layer less than 3, could not fill a RPCSeedFinder" << endl;
}

void RPCSeedrecHitFinder::complete ( unsigned int  LayerIndex )

private

Definition at line 82 of file RPCSeedrecHitFinder.cc.

References funct::abs(), SplitLinear::begin, RPCRecHit::BunchX(), L1TStage2uGTEmulatorClient_cff::BX, TrackingRecHit::clone(), RPCRecHit::clusterSize(), PixelTestBeamValidation_cfi::ClusterSize, gather_cfg::cout, srCondWrite_cfg::deltaPhi, TrackingRecHit::hit(), colinearityKinematic::Phi, PV3DBase< T, PVType, FrameType >::phi(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                          {
  for (MuonRecHitContainer::const_iterator it = recHitsRPC[LayersinRPC[LayerIndex]]->begin();
       it != recHitsRPC[LayersinRPC[LayerIndex]]->end();
       it++) {
    cout << "Completing layer[" << LayersinRPC[LayerIndex] << "]." << endl;

    // Check validation
    if (!(*it)->isValid())
      continue;

    // Check BX range, be sure there is only RPCRecHit in the MuonRecHitContainer when use the dynamic_cast
    TrackingRecHit* thisTrackingRecHit = (*it)->hit()->clone();
    // Should also delete the RPCRecHit object cast by dynamic_cast<> ?
    RPCRecHit* thisRPCRecHit = dynamic_cast<RPCRecHit*>(thisTrackingRecHit);
    int BX = thisRPCRecHit->BunchX();
    int ClusterSize = thisRPCRecHit->clusterSize();
    delete thisTrackingRecHit;
    // Check BX
    if ((unsigned int)abs(BX) > BxRange)
      continue;
    // Check cluster size
    bool Clustercheck = false;
    if (ClusterSet.empty())
      Clustercheck = true;
    for (std::vector<int>::const_iterator CluIter = ClusterSet.begin(); CluIter != ClusterSet.end(); CluIter++)
      if (ClusterSize == (*CluIter))
        Clustercheck = true;
    if (Clustercheck != true)
      continue;
    // Check the recHits Phi range
    GlobalPoint pos = (*it)->globalPosition();
    double Phi = pos.phi();
    cout << "Phi: " << Phi << endl;
    // The recHits should locate in some phi range
    therecHits.push_back(*it);
    double deltaPhi = getdeltaPhifromrecHits();
    cout << "Delta phi: " << deltaPhi << endl;
    therecHits.pop_back();
    if (deltaPhi > MaxDeltaPhi)
      continue;

    // If pass all, add to the seed
    therecHits.push_back(*it);
    cout << "RecHit's global position: " << pos.x() << ", " << pos.y() << ", " << pos.z() << endl;

    // Check if this recHit is the last one in the seed
    // If it is the last one, calculate the seed
    if (LayerIndex == (LayersinRPC.size() - 1)) {
      cout << "Check and fill one seed." << endl;
      checkandfill();
    }
    // If it is not the last one, continue to fill the seed from other layers
    else
      complete(LayerIndex + 1);

    // Remember to pop the recHit before add another one from the same layer!
    therecHits.pop_back();
  }
}

void RPCSeedrecHitFinder::configure

(

const edm::ParameterSet &

iConfig

)

Definition at line 36 of file RPCSeedrecHitFinder.cc.

References edm::ParameterSet::getParameter(), and RPCSeedFinder::isConfigured.

                                                                  {
  // Set the configuration
  BxRange = iConfig.getParameter<unsigned int>("BxRange");
  MaxDeltaPhi = iConfig.getParameter<double>("MaxDeltaPhi");
  ClusterSet = iConfig.getParameter<std::vector<int> >("ClusterSet");

  // Set the signal open
  isConfigured = true;
}

void RPCSeedrecHitFinder::fillrecHits

(

)

Definition at line 63 of file RPCSeedrecHitFinder.cc.

References gather_cfg::cout, RPCSeedFinder::isConfigured, RPCSeedFinder::isOutputset, and isotrackApplyRegressor::k.

                                      {
  if (isLayerset == false || isConfigured == false || isOutputset == false || isInputset == false) {
    cout << "Not set the IO or not configured yet" << endl;
    return;
  }
  cout << "Now fill recHits from Layers: ";
  for (unsigned int k = 0; k < LayersinRPC.size(); k++)
    cout << LayersinRPC[k] << " ";
  cout << endl;
  unsigned int LayerIndex = 0;
  therecHits.clear();
  complete(LayerIndex);

  // Unset the signal
  LayersinRPC.clear();
  therecHits.clear();
  isLayerset = false;
}

double RPCSeedrecHitFinder::getdeltaPhifromrecHits ( )

private

Definition at line 142 of file RPCSeedrecHitFinder.cc.

References gather_cfg::cout, srCondWrite_cfg::deltaPhi, lessPhi(), M_PI, and dqmiodumpmetadata::n.

                                                   {
  ConstMuonRecHitContainer sortRecHits = therecHits;
  sort(sortRecHits.begin(), sortRecHits.end(), lessPhi);
  cout << "Sorted recHit's Phi: ";
  for (ConstMuonRecHitContainer::const_iterator iter = sortRecHits.begin(); iter != sortRecHits.end(); iter++)
    cout << (*iter)->globalPosition().phi() << ", ";
  cout << endl;
  // Calculate the deltaPhi, take care Geom::Phi always in range [-pi,pi)
  // In case of some deltaPhi larger then Pi, use value() in Geom::Phi to get the true value in radians of Phi, then do the calculation
  double deltaPhi = 0;
  if (sortRecHits.size() <= 1)
    return deltaPhi;
  if (sortRecHits.size() == 2) {
    ConstMuonRecHitContainer::const_iterator iter1 = sortRecHits.begin();
    ConstMuonRecHitContainer::const_iterator iter2 = sortRecHits.begin();
    iter2++;
    deltaPhi = (((*iter2)->globalPosition().phi().value() - (*iter1)->globalPosition().phi().value()) > M_PI)
                   ? (2 * M_PI - ((*iter2)->globalPosition().phi().value() - (*iter1)->globalPosition().phi().value()))
                   : ((*iter2)->globalPosition().phi().value() - (*iter1)->globalPosition().phi().value());
    return deltaPhi;
  } else {
    deltaPhi = 2 * M_PI;
    int n = 0;
    for (ConstMuonRecHitContainer::const_iterator iter = sortRecHits.begin(); iter != sortRecHits.end(); iter++) {
      cout << "Before this loop deltaPhi is " << deltaPhi << endl;
      n++;
      double deltaPhi_more = 0;
      double deltaPhi_less = 0;
      if (iter == sortRecHits.begin()) {
        cout << "Calculateing frist loop..." << endl;
        ConstMuonRecHitContainer::const_iterator iter_more = ++iter;
        --iter;
        ConstMuonRecHitContainer::const_iterator iter_less = sortRecHits.end();
        --iter_less;
        cout << "more_Phi: " << (*iter_more)->globalPosition().phi()
             << ", less_Phi: " << (*iter_less)->globalPosition().phi()
             << ", iter_Phi: " << (*iter)->globalPosition().phi() << endl;
        deltaPhi_more =
            (2 * M_PI) - ((*iter_more)->globalPosition().phi().value() - (*iter)->globalPosition().phi().value());
        deltaPhi_less = (*iter_less)->globalPosition().phi().value() - (*iter)->globalPosition().phi().value();
      } else if (iter == (--sortRecHits.end())) {
        cout << "Calculateing last loop..." << endl;
        ConstMuonRecHitContainer::const_iterator iter_less = --iter;
        ++iter;
        ConstMuonRecHitContainer::const_iterator iter_more = sortRecHits.begin();
        cout << "more_Phi: " << (*iter_more)->globalPosition().phi()
             << ", less_Phi: " << (*iter_less)->globalPosition().phi()
             << ", iter_Phi: " << (*iter)->globalPosition().phi() << endl;
        deltaPhi_less =
            (2 * M_PI) - ((*iter)->globalPosition().phi().value() - (*iter_less)->globalPosition().phi().value());
        deltaPhi_more = (*iter)->globalPosition().phi().value() - (*iter_more)->globalPosition().phi().value();
      } else {
        cout << "Calculateing " << n << "st loop..." << endl;
        ConstMuonRecHitContainer::const_iterator iter_less = --iter;
        ++iter;
        ConstMuonRecHitContainer::const_iterator iter_more = ++iter;
        --iter;
        cout << "more_Phi: " << (*iter_more)->globalPosition().phi()
             << ", less_Phi: " << (*iter_less)->globalPosition().phi()
             << ", iter_Phi: " << (*iter)->globalPosition().phi() << endl;
        deltaPhi_less =
            (2 * M_PI) - ((*iter)->globalPosition().phi().value() - (*iter_less)->globalPosition().phi().value());
        deltaPhi_more =
            (2 * M_PI) - ((*iter_more)->globalPosition().phi().value() - (*iter)->globalPosition().phi().value());
      }
      if (deltaPhi > deltaPhi_more)
        deltaPhi = deltaPhi_more;
      if (deltaPhi > deltaPhi_less)
        deltaPhi = deltaPhi_less;

      cout << "For this loop deltaPhi_more is " << deltaPhi_more << endl;
      cout << "For this loop deltaPhi_less is " << deltaPhi_less << endl;
      cout << "For this loop deltaPhi is " << deltaPhi << endl;
    }
    return deltaPhi;
  }
}

void RPCSeedrecHitFinder::setInput

(

MuonRecHitContainer(&)

recHits[12]

)

Definition at line 46 of file RPCSeedrecHitFinder.cc.

References mps_fire::i, HLT_FULL_cff::recHits, and RPCLayerNumber.

                                                                                 {
  for (unsigned int i = 0; i < RPCLayerNumber; i++)
    recHitsRPC[i] = &recHits[i];
  isInputset = true;
}

void RPCSeedrecHitFinder::setLayers

(

const std::vector< unsigned int > &

Layers

)

Definition at line 58 of file RPCSeedrecHitFinder.cc.

                                                                         {
  LayersinRPC = Layers;
  isLayerset = true;
}

void RPCSeedrecHitFinder::setOutput

(

RPCSeedFinder *

Seed

)

Definition at line 53 of file RPCSeedrecHitFinder.cc.

References RPCSeedFinder::isOutputset.

                                                       {
  theSeed = Seed;
  isOutputset = true;
}

void RPCSeedrecHitFinder::unsetInput

(

)

Definition at line 52 of file RPCSeedrecHitFinder.cc.

{ isInputset = false; }

Member Data Documentation

unsigned int RPCSeedrecHitFinder::BxRange

private

Definition at line 51 of file RPCSeedrecHitFinder.h.

std::vector<int> RPCSeedrecHitFinder::ClusterSet

private

Definition at line 53 of file RPCSeedrecHitFinder.h.

bool RPCSeedrecHitFinder::isConfigured

private

Definition at line 56 of file RPCSeedrecHitFinder.h.

bool RPCSeedrecHitFinder::isInputset

private

Definition at line 57 of file RPCSeedrecHitFinder.h.

bool RPCSeedrecHitFinder::isLayerset

private

Definition at line 55 of file RPCSeedrecHitFinder.h.

bool RPCSeedrecHitFinder::isOutputset

private

Definition at line 58 of file RPCSeedrecHitFinder.h.

std::vector<unsigned int> RPCSeedrecHitFinder::LayersinRPC

private

Definition at line 60 of file RPCSeedrecHitFinder.h.

double RPCSeedrecHitFinder::MaxDeltaPhi

private

Definition at line 52 of file RPCSeedrecHitFinder.h.

MuonRecHitContainer* RPCSeedrecHitFinder::recHitsRPC[12]

private

Definition at line 61 of file RPCSeedrecHitFinder.h.

ConstMuonRecHitContainer RPCSeedrecHitFinder::therecHits

private

Definition at line 62 of file RPCSeedrecHitFinder.h.

RPCSeedFinder* RPCSeedrecHitFinder::theSeed

private

Definition at line 63 of file RPCSeedrecHitFinder.h.