|
|
|
#include <RPCSeedrecHitFinder.h>
Definition at line 29 of file RPCSeedrecHitFinder.h.
typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer RPCSeedrecHitFinder::ConstMuonRecHitContainer [private] |
Definition at line 34 of file RPCSeedrecHitFinder.h.
typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer RPCSeedrecHitFinder::ConstMuonRecHitPointer [private] |
Definition at line 32 of file RPCSeedrecHitFinder.h.
typedef MuonTransientTrackingRecHit::MuonRecHitContainer RPCSeedrecHitFinder::MuonRecHitContainer [private] |
Definition at line 33 of file RPCSeedrecHitFinder.h.
typedef MuonTransientTrackingRecHit::MuonRecHitPointer RPCSeedrecHitFinder::MuonRecHitPointer [private] |
Definition at line 31 of file RPCSeedrecHitFinder.h.
| RPCSeedrecHitFinder::RPCSeedrecHitFinder | ( | ) |
Definition at line 23 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 37 of file RPCSeedrecHitFinder.cc.
{
}
| void RPCSeedrecHitFinder::checkandfill | ( | ) | [private] |
Definition at line 233 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 96 of file RPCSeedrecHitFinder.cc.
References abs, begin, RPCRecHit::BunchX(), TrackingRecHit::clone(), RPCRecHit::clusterSize(), gather_cfg::cout, Geom::deltaPhi(), PV3DBase< T, PVType, FrameType >::phi(), colinearityKinematic::Phi, pos, 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.size() == 0)
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 41 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 75 of file RPCSeedrecHitFinder.cc.
References gather_cfg::cout, RPCSeedFinder::isConfigured, RPCSeedFinder::isOutputset, and gen::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 157 of file RPCSeedrecHitFinder.cc.
References gather_cfg::cout, Geom::deltaPhi(), lessPhi(), M_PI, n, and python::multivaluedict::sort().
{
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[RPCLayerNumber] | ) |
Definition at line 52 of file RPCSeedrecHitFinder.cc.
References i, 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 69 of file RPCSeedrecHitFinder.cc.
{
LayersinRPC = Layers;
isLayerset = true;
}
| void RPCSeedrecHitFinder::setOutput | ( | RPCSeedFinder * | Seed | ) |
Definition at line 63 of file RPCSeedrecHitFinder.cc.
References RPCSeedFinder::isOutputset.
{
theSeed = Seed;
isOutputset = true;
}
| void RPCSeedrecHitFinder::unsetInput | ( | ) |
Definition at line 59 of file RPCSeedrecHitFinder.cc.
{
isInputset = false;
}
unsigned int RPCSeedrecHitFinder::BxRange [private] |
Definition at line 53 of file RPCSeedrecHitFinder.h.
std::vector<int> RPCSeedrecHitFinder::ClusterSet [private] |
Definition at line 55 of file RPCSeedrecHitFinder.h.
bool RPCSeedrecHitFinder::isConfigured [private] |
Definition at line 58 of file RPCSeedrecHitFinder.h.
bool RPCSeedrecHitFinder::isInputset [private] |
Definition at line 59 of file RPCSeedrecHitFinder.h.
bool RPCSeedrecHitFinder::isLayerset [private] |
Definition at line 57 of file RPCSeedrecHitFinder.h.
bool RPCSeedrecHitFinder::isOutputset [private] |
Definition at line 60 of file RPCSeedrecHitFinder.h.
std::vector<unsigned int> RPCSeedrecHitFinder::LayersinRPC [private] |
Definition at line 62 of file RPCSeedrecHitFinder.h.
double RPCSeedrecHitFinder::MaxDeltaPhi [private] |
Definition at line 54 of file RPCSeedrecHitFinder.h.
MuonRecHitContainer* RPCSeedrecHitFinder::recHitsRPC[RPCLayerNumber] [private] |
Definition at line 63 of file RPCSeedrecHitFinder.h.
Definition at line 64 of file RPCSeedrecHitFinder.h.
RPCSeedFinder* RPCSeedrecHitFinder::theSeed [private] |
Definition at line 65 of file RPCSeedrecHitFinder.h.
1.7.3