RoadMaker Class Reference

#include <RoadMaker.h>

Public Types
enum	GeometryStructure {   FullDetector, FullDetectorII, MTCC, TIF,   TIFTOB, TIFTIB, TIFTIBTOB, TIFTOBTEC,   P5 }
enum	SeedingType { FourRingSeeds, TwoRingSeeds }

Public Member Functions
std::pair< Roads::RoadSeed, Roads::RoadSet >	AddInnerSeedRing (std::pair< Roads::RoadSeed, Roads::RoadSet > input)
std::pair< Roads::RoadSeed, Roads::RoadSet >	AddOuterSeedRing (std::pair< Roads::RoadSeed, Roads::RoadSet > input)
bool	AddRoad (Roads::RoadSeed seed, Roads::RoadSet set)
Roads::RoadSeed	CloneSeed (Roads::RoadSeed seed)
void	collectInnerSeedRings ()
void	collectInnerSeedRings1 ()
void	collectInnerSeedRings2 ()
void	collectInnerTECSeedRings ()
void	collectInnerTECSeedRings1 ()
void	collectInnerTECSeedRings2 ()
void	collectInnerTIBSeedRings ()
void	collectInnerTIBSeedRings1 ()
void	collectInnerTIBSeedRings2 ()
void	collectInnerTIDSeedRings ()
void	collectInnerTIDSeedRings1 ()
void	collectInnerTIDSeedRings2 ()
void	collectInnerTOBSeedRings ()
void	collectInnerTOBSeedRings1 ()
void	collectInnerTOBSeedRings2 ()
void	collectOuterSeedRings ()
void	collectOuterSeedRings1 ()
void	collectOuterTECSeedRings ()
void	collectOuterTECSeedRings1 ()
void	collectOuterTIBSeedRings ()
void	collectOuterTIBSeedRings1 ()
void	collectOuterTOBSeedRings ()
void	collectOuterTOBSeedRings1 ()
bool	CompatibleWithLines (std::vector< std::pair< double, double > > lines, const Ring *ring)
void	constructRoads ()
Roads *	getRoads ()
std::vector< std::pair< double, double > >	LinesThroughRingAndBS (const Ring *ring)
std::vector< std::pair< double, double > >	LinesThroughRings (const Ring *ring1, const Ring *ring2)
bool	RingInBarrel (const Ring *ring)
Roads::RoadSet	RingsCompatibleWithSeed (Roads::RoadSeed seed)
bool	RingsOnSameLayer (const Ring *ring1, const Ring *ring2)
bool	RingsOnSameLayer (std::pair< const Ring *, const Ring * > seed1, std::pair< const Ring *, const Ring * > seed2)
	RoadMaker (const Rings *rings, GeometryStructure structure=FullDetector, SeedingType seedingType=FourRingSeeds)
bool	SameRoadSet (Roads::RoadSet set1, Roads::RoadSet set2)
Roads::RoadSet	SortRingsIntoLayers (std::vector< const Ring * > input)
	~RoadMaker ()

Private Attributes
std::vector< const Ring * >	innerSeedRings1_
std::vector< const Ring * >	innerSeedRings2_
std::vector< const Ring * >	innerSeedRings_
std::vector< const Ring * >	outerSeedRings1_
std::vector< const Ring * >	outerSeedRings_
const Rings *	rings_
Roads *	roads_
SeedingType	seedingType_
GeometryStructure	structure_
float	zBS_

Detailed Description

Definition at line 29 of file RoadMaker.h.

Member Enumeration Documentation

enum RoadMaker::GeometryStructure

Enumerator
FullDetector
FullDetectorII
MTCC
TIF
TIFTOB
TIFTIB
TIFTIBTOB
TIFTOBTEC
P5

Definition at line 33 of file RoadMaker.h.

                         {
    FullDetector,
    FullDetectorII,
    MTCC,
    TIF,
    TIFTOB,
    TIFTIB,
    TIFTIBTOB,
    TIFTOBTEC,
    P5
  };

enum RoadMaker::SeedingType

Enumerator
FourRingSeeds
TwoRingSeeds

Definition at line 45 of file RoadMaker.h.

                   {
    FourRingSeeds,
    TwoRingSeeds
  };

Constructor & Destructor Documentation

RoadMaker::RoadMaker	(	const Rings *	rings,
		RoadMaker::GeometryStructure	structure = FullDetector,
		RoadMaker::SeedingType	seedingType = FourRingSeeds
	)

Definition at line 41 of file RoadMaker.cc.

References constructRoads(), FullDetector, FullDetectorII, MTCC, P5, roads_, structure_, TIF, TIFTIB, TIFTIBTOB, TIFTOB, TIFTOBTEC, and zBS_.

                                               : 
  rings_(rings), 
  structure_(structure),
  seedingType_(seedingType) {

    // define beamspot z extnsion
    zBS_ = 0;
    if(structure_==FullDetector) {
      zBS_ = 3*5.3; // cm
    } else if (structure_==FullDetectorII) {
      zBS_ = 3*5.3; // cm
    } else if (structure_==MTCC) {
      zBS_ = 3*75.; // cm
    } else if (structure_==TIF) {
      zBS_ = 3*75.; // cm
    } else if (structure_==TIFTIB) {
      zBS_ = 5*75.; // cm
    } else if (structure_==TIFTOB) {
      zBS_ = 5*75.; // cm
    } else if (structure_==TIFTIBTOB) {
      zBS_ = 5*75.; // cm
    } else if (structure_==TIFTOBTEC) {
      zBS_ = 5*75.; // cm
    } else if (structure_==P5) {
      zBS_ = 5*75.; // cm
    }

    roads_ = new Roads();

    constructRoads();

  }

RoadMaker::~RoadMaker

(

)

Definition at line 76 of file RoadMaker.cc.

                      { 
  
}

Member Function Documentation

std::pair< Roads::RoadSeed, Roads::RoadSet > RoadMaker::AddInnerSeedRing

(

std::pair< Roads::RoadSeed, Roads::RoadSet >

input

)

Definition at line 2302 of file RoadMaker.cc.

References AddRoad(), CloneSeed(), innerSeedRings_, LaserDQM_cfg::input, query::result, RingsCompatibleWithSeed(), RingsOnSameLayer(), SameRoadSet(), and runtimedef::set().

                                                                                                           {
  //
  // add another inner seed ring
  // check for inner seed ring which is in input RoadSet and on another layer than the first seed ring
  for ( std::vector<const Ring*>::iterator innerSeedRing = innerSeedRings_.begin();
    innerSeedRing != innerSeedRings_.end();
    ++innerSeedRing ) {
    for ( Roads::RoadSet::iterator roadSetVector = input.second.begin();
      roadSetVector != input.second.end();
      ++roadSetVector ) {
      for ( std::vector<const Ring*>::iterator roadSetRing = roadSetVector->begin();
        roadSetRing != roadSetVector->end();
        ++roadSetRing) {
    // check for same ring
    if ( (*innerSeedRing) == (*roadSetRing) ) {
      // check that new ring is not on same layer as previous inner seed rings
      bool onSameLayer = false;
      for ( std::vector<const Ring*>::iterator roadSeedRing = input.first.first.begin();
        roadSeedRing != input.first.first.end();
        ++roadSeedRing ) {
        if ( RingsOnSameLayer((*roadSeedRing),(*innerSeedRing)) ) {
          onSameLayer = true;
        }
      }
      if ( !onSameLayer ) {
          
        Roads::RoadSeed seed = CloneSeed(input.first);
        seed.first.push_back((*innerSeedRing));
          
        Roads::RoadSet set = RingsCompatibleWithSeed(seed);
          
        if ( SameRoadSet(input.second,set) ) {

          AddRoad(seed,set);

          std::pair<Roads::RoadSeed, Roads::RoadSet> result(seed,set);

          return result;
        }
      }
    }
      }
    }
  }

  // if no ring could be added, return input
  return input;

}

std::pair< Roads::RoadSeed, Roads::RoadSet > RoadMaker::AddOuterSeedRing

(

std::pair< Roads::RoadSeed, Roads::RoadSet >

input

)

Definition at line 2352 of file RoadMaker.cc.

References AddRoad(), CloneSeed(), LaserDQM_cfg::input, outerSeedRings_, query::result, RingsCompatibleWithSeed(), RingsOnSameLayer(), and runtimedef::set().

                                                                                                           {
  //
  // add another outer seed ring
  // check for outer seed ring which is in input RoadSet and on another layer than the first seed ring
  for ( std::vector<const Ring*>::iterator outerSeedRing = outerSeedRings_.begin();
    outerSeedRing != outerSeedRings_.end();
    ++outerSeedRing ) {
    for ( Roads::RoadSet::iterator roadSetVector = input.second.begin();
      roadSetVector != input.second.end();
      ++roadSetVector ) {
      for ( std::vector<const Ring*>::iterator roadSetRing = roadSetVector->begin();
        roadSetRing != roadSetVector->end();
        ++roadSetRing) {
    // check for same ring
    if ( (*outerSeedRing) == (*roadSetRing) ) {
      // check that new ring is not on same layer as previous outer seed rings
      bool onSameLayer = false;
      for ( std::vector<const Ring*>::iterator roadSeedRing = input.first.second.begin();
        roadSeedRing != input.first.second.end();
        ++roadSeedRing ) {
        if ( RingsOnSameLayer((*roadSeedRing),(*outerSeedRing)) ) {
          onSameLayer = true;
        }
      }
      if ( !onSameLayer ) {
          
        Roads::RoadSeed seed = CloneSeed(input.first);
        seed.second.push_back((*outerSeedRing));
          
        Roads::RoadSet set = RingsCompatibleWithSeed(seed);
          
        AddRoad(seed,set);

        std::pair<Roads::RoadSeed, Roads::RoadSet> result(seed,set);

        return result;
      }
    }
      }
    }
  }

  // if no ring could be added, return input
  return input;

}

bool RoadMaker::AddRoad	(	Roads::RoadSeed	seed,
		Roads::RoadSet	set
	)

Definition at line 2234 of file RoadMaker.cc.

References Roads::begin(), Roads::end(), Roads::erase(), Roads::insert(), query::result, roads_, and python.multivaluedict::sort().

Referenced by AddInnerSeedRing(), and AddOuterSeedRing().

                              {
  //
  // add road
  // check if seed rings are included in seeds of other roads or if they are new
  // add road if new
  // take road with larger seed and discard older seed
  // assumption: increasing number of seed rings while adding during the program
  //

  // return value
  bool result = true;

  for ( Roads::iterator existingRoad = roads_->begin();
    existingRoad != roads_->end();
    ++existingRoad ) {

    Roads::RoadSeed existingSeed = existingRoad->first;

    // check if existing inner seed is included in new seed
    unsigned int includedInner = 0;
    unsigned int includedOuter = 0;
    for ( std::vector<const Ring*>::iterator existingInnerRing = existingSeed.first.begin();
      existingInnerRing != existingSeed.first.end();
      ++existingInnerRing ) {
      bool ringIncluded = false;
      for ( std::vector<const Ring*>::iterator innerRing = seed.first.begin();
        innerRing != seed.first.end();
        ++innerRing ) {
    if ( (*existingInnerRing) == (*innerRing) ) {
      ringIncluded = true;
    }
      }
      if ( ringIncluded ) {
    ++includedInner;
      }
    }
    for ( std::vector<const Ring*>::iterator existingOuterRing = existingSeed.second.begin();
      existingOuterRing != existingSeed.second.end();
      ++existingOuterRing ) {
      bool ringIncluded = false;
      for ( std::vector<const Ring*>::iterator outerRing = seed.second.begin();
        outerRing != seed.second.end();
        ++outerRing ) {
    if ( (*existingOuterRing) == (*outerRing) ) {
      ringIncluded = true;
    }
      }
      if ( ringIncluded ) {
    ++includedOuter;
      }
    }

    if ( includedInner ==  existingSeed.first.size() &&
     includedOuter ==  existingSeed.second.size() ) {
      // existing road included in new road, remove
      roads_->erase(existingRoad);
    }
  }

  // sort seeds
  std::sort(seed.first.begin(),seed.first.end(),SortRingsByZR());
  std::sort(seed.second.begin(),seed.second.end(),SortRingsByZR());

  roads_->insert(seed,set);
  return result;
}

Roads::RoadSeed RoadMaker::CloneSeed

(

Roads::RoadSeed

seed

)

Definition at line 2212 of file RoadMaker.cc.

References query::result, and relativeConstraints::ring.

Referenced by AddInnerSeedRing(), and AddOuterSeedRing().

                                                     {
  //
  // clone seed
  //

  Roads::RoadSeed result;

  for ( std::vector<const Ring*>::iterator ring = seed.first.begin();
    ring != seed.first.end();
    ++ring ) {
    result.first.push_back((*ring));
  }

  for ( std::vector<const Ring*>::iterator ring = seed.second.begin();
    ring != seed.second.end();
    ++ring ) {
    result.second.push_back((*ring));
  }

  return result;
}

void RoadMaker::collectInnerSeedRings

(

)

Definition at line 258 of file RoadMaker.cc.

References collectInnerTECSeedRings(), collectInnerTIBSeedRings(), collectInnerTIDSeedRings(), collectInnerTOBSeedRings(), FullDetector, FullDetectorII, innerSeedRings_, LogDebug, MTCC, P5, structure_, TIF, TIFTIB, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by constructRoads().

                                      {

  if(structure_==FullDetector || structure_==FullDetectorII) {
    collectInnerTIBSeedRings();
    collectInnerTIDSeedRings();
    collectInnerTECSeedRings(); 
  } else if(structure_==MTCC) {
    collectInnerTIBSeedRings();
  } else if(structure_==TIF) {
    collectInnerTIBSeedRings();
    collectInnerTOBSeedRings();
    collectInnerTIDSeedRings();
    collectInnerTECSeedRings(); 
  } else if(structure_==TIFTOB) {
    collectInnerTOBSeedRings();
  } else if(structure_==TIFTIB) {
    collectInnerTIBSeedRings();
    collectInnerTIDSeedRings();
  } else if(structure_==TIFTIBTOB) {
    collectInnerTIBSeedRings();
    collectInnerTIDSeedRings();
    collectInnerTOBSeedRings();
  } else if(structure_==TIFTOBTEC) {
    collectInnerTOBSeedRings();
    collectInnerTECSeedRings();
  } else if(structure_==P5) {
    collectInnerTOBSeedRings();
    collectInnerTECSeedRings();
  }

  LogDebug("RoadSearch") << "collected " << innerSeedRings_.size() << " inner seed rings"; 
  
}

void RoadMaker::collectInnerSeedRings1

(

)

Definition at line 613 of file RoadMaker.cc.

References collectInnerTECSeedRings1(), collectInnerTIBSeedRings1(), collectInnerTIDSeedRings1(), collectInnerTOBSeedRings1(), FullDetector, FullDetectorII, innerSeedRings1_, LogDebug, MTCC, structure_, TIF, TIFTIB, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by constructRoads().

                                       {

  if(structure_==FullDetector || structure_==FullDetectorII) {
    collectInnerTIBSeedRings1();
    collectInnerTIDSeedRings1();
    collectInnerTECSeedRings1(); 
  } else if(structure_==MTCC) {
    collectInnerTIBSeedRings1();
  } else if(structure_==TIF) {
    collectInnerTIBSeedRings1();
    collectInnerTOBSeedRings1();
    collectInnerTIDSeedRings1();
    collectInnerTECSeedRings1(); 
 } else if(structure_==TIFTOB) {
    collectInnerTOBSeedRings1();
  } else if(structure_==TIFTIB) {
    collectInnerTIBSeedRings1();
    collectInnerTIDSeedRings1();
  } else if(structure_==TIFTIBTOB) {
    collectInnerTIBSeedRings1();
    collectInnerTIDSeedRings1();
    collectInnerTOBSeedRings1();
  } else if(structure_==TIFTOBTEC) {
    collectInnerTOBSeedRings1();
    collectInnerTECSeedRings1();
   }

  LogDebug("RoadSearch") << "collected " << innerSeedRings1_.size() << " inner seed rings"; 
  
}

void RoadMaker::collectInnerSeedRings2

(

)

Definition at line 930 of file RoadMaker.cc.

References collectInnerTECSeedRings2(), collectInnerTIBSeedRings2(), collectInnerTIDSeedRings2(), collectInnerTOBSeedRings2(), FullDetector, FullDetectorII, innerSeedRings2_, LogDebug, MTCC, structure_, TIF, TIFTIB, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by constructRoads().

                                       {

  if(structure_==FullDetector) {
    collectInnerTIBSeedRings2();
    collectInnerTIDSeedRings2();
    collectInnerTECSeedRings2(); 
  } else if(structure_==FullDetectorII) {
    collectInnerTIBSeedRings2();
    collectInnerTIDSeedRings2();
    collectInnerTECSeedRings2(); 
  } else if(structure_==MTCC) {
    collectInnerTIBSeedRings2();
  } else if(structure_==TIF) {
    collectInnerTIBSeedRings2();
    collectInnerTOBSeedRings2();
    collectInnerTIDSeedRings2();
    collectInnerTECSeedRings2(); 
  } else if(structure_==TIFTOB) {
    collectInnerTOBSeedRings2();
  } else if(structure_==TIFTIB) {
    collectInnerTIBSeedRings2();
    collectInnerTIDSeedRings2();
  } else if(structure_==TIFTIBTOB) {
    collectInnerTIBSeedRings2();
    collectInnerTIDSeedRings2();
    collectInnerTOBSeedRings2();
  } else if(structure_==TIFTOBTEC) {
    collectInnerTOBSeedRings2();
    collectInnerTECSeedRings2();
  }


  LogDebug("RoadSearch") << "collected " << innerSeedRings2_.size() << " inner seed rings"; 
  
}

void RoadMaker::collectInnerTECSeedRings

(

)

Definition at line 413 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTECRing(), innerSeedRings_, LogDebug, P5, relativeConstraints::ring, rings_, structure_, TIF, and TIFTOBTEC.

Referenced by collectInnerSeedRings().

                                         {

  // TEC
  unsigned int counter = 0;

  unsigned int fw_bw_min       = 1;
  unsigned int fw_bw_max       = 1;
  unsigned int wheel_min       = 1;
  unsigned int wheel_max       = 9;
  unsigned int ring_min[9];
  unsigned int ring_max[9];
  
  if(structure_==FullDetector || structure_==FullDetectorII) {
    fw_bw_max = 3;
    // TEC WHEEL 1
    ring_min[1] = 1;
    ring_max[1] = 3;
    // TEC WHEEL 2
    ring_min[2] = 1;
    ring_max[2] = 3;
    // TEC WHEEL 3
    ring_min[3] = 1;
    ring_max[3] = 3;
    // TEC WHEEL 4
    ring_min[4] = 2;
    ring_max[4] = 3;
    // TEC WHEEL 5
    ring_min[5] = 2;
    ring_max[5] = 3;
    // TEC WHEEL 6
    ring_min[6] = 2;
    ring_max[6] = 3;   
    // TEC WHEEL 7
    ring_min[7] = 1;
    ring_max[7] = 1;
    // TEC WHEEL 8
    ring_min[8] = 1;
    ring_max[8] = 1;
  } else if (structure_==TIFTOBTEC || structure_==TIF) {
    fw_bw_min = 2;
    fw_bw_max = 3;
    // TEC WHEEL 1
    ring_min[1] = 1;
    ring_max[1] = 4;
    // TEC WHEEL 2
    ring_min[2] = 1;
    ring_max[2] = 4;
    // TEC WHEEL 3
    ring_min[3] = 1;
    ring_max[3] = 4;
    // TEC WHEEL 4
    ring_min[4] = 2;
    ring_max[4] = 4;
    // TEC WHEEL 5
    ring_min[5] = 2;
    ring_max[5] = 4;
    // TEC WHEEL 6
    ring_min[6] = 2;
    ring_max[6] = 4;
    // TEC WHEEL 7
    ring_min[7] = 3;
    ring_max[7] = 4;
    // TEC WHEEL 8
    ring_min[8] = 3;
    ring_max[8] = 4;
  } else if (structure_==P5) {
    fw_bw_min = 1;
    fw_bw_max = 3;
    // TEC WHEEL 1
    ring_min[1] = 4;
    ring_max[1] = 6;
    // TEC WHEEL 2
    ring_min[2] = 4;
    ring_max[2] = 6;
    // TEC WHEEL 3
    ring_min[3] = 4;
    ring_max[3] = 6;
    // TEC WHEEL 4
    ring_min[4] = 4;
    ring_max[4] = 6;
    // TEC WHEEL 5
    ring_min[5] = 4;
    ring_max[5] = 6;
    // TEC WHEEL 6
    ring_min[6] = 4;
    ring_max[6] = 6;
    // TEC WHEEL 7
    ring_min[7] = 4;
    ring_max[7] = 6;
    // TEC WHEEL 8
    ring_min[8] = 4;
    ring_max[8] = 6;
  }

  for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
    for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
      for ( unsigned int ring = ring_min[wheel]; ring < ring_max[wheel]; ++ring ) {
    const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,ring);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TEC inner seed ring with index: " << temp_ring->getindex(); 
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TEC inner seed rings"; 

}

void RoadMaker::collectInnerTECSeedRings1

(

)

Definition at line 764 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTECRing(), innerSeedRings1_, LogDebug, relativeConstraints::ring, rings_, structure_, TIF, and TIFTOBTEC.

Referenced by collectInnerSeedRings1().

                                          {

  // TEC
  unsigned int counter = 0;

  unsigned int fw_bw_min       = 1;
  unsigned int fw_bw_max       = 1;
  unsigned int wheel_min       = 1;
  unsigned int wheel_max       = 4;
  unsigned int ring_min[9];
  unsigned int ring_max[9];

  if(structure_==FullDetector || structure_==FullDetectorII) {
    fw_bw_min = 1;
    fw_bw_max = 3;
    // TEC WHEEL 1
    ring_min[1] = 1;
    ring_max[1] = 2;
    // TEC WHEEL 2
    ring_min[2] = 1;
    ring_max[2] = 2;
    // TEC WHEEL 3
    ring_min[3] = 1;
    ring_max[3] = 2;
    // TEC WHEEL 4
    ring_min[4] = 1;
    ring_max[4] = 1;
    // TEC WHEEL 5
    ring_min[5] = 1;
    ring_max[5] = 1;
    // TEC WHEEL 6
    ring_min[6] = 1;
    ring_max[6] = 1;
    // TEC WHEEL 7
    ring_min[7] = 1;
    ring_max[7] = 1;
    // TEC WHEEL 8
    ring_min[8] = 1;
    ring_max[8] = 1;
  } else if (structure_==TIFTOBTEC || structure_==TIF) {  
    fw_bw_min = 2;
    fw_bw_max = 3;
    // TEC WHEEL 1
    ring_min[1] = 1;
    ring_max[1] = 2;
    // TEC WHEEL 2
    ring_min[2] = 1;
    ring_max[2] = 2;
    // TEC WHEEL 3
    ring_min[3] = 1;
    ring_max[3] = 2;
    // TEC WHEEL 4
    ring_min[4] = 1;
    ring_max[4] = 1;
    // TEC WHEEL 5
    ring_min[5] = 1;
    ring_max[5] = 1;
    // TEC WHEEL 6
    ring_min[6] = 1;
    ring_max[6] = 1;
    // TEC WHEEL 7
    ring_min[7] = 1;
    ring_max[7] = 1;
    // TEC WHEEL 8
    ring_min[8] = 1;
    ring_max[8] = 1;
  }

  for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
    for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
      for ( unsigned int ring = ring_min[wheel]; ring < ring_max[wheel]; ++ring ) {
    const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,ring);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TEC inner seed ring with index: " << temp_ring->getindex(); 
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TEC inner seed rings"; 

}

void RoadMaker::collectInnerTECSeedRings2

(

)

Definition at line 1095 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTECRing(), innerSeedRings2_, LogDebug, relativeConstraints::ring, rings_, structure_, TIF, and TIFTOBTEC.

Referenced by collectInnerSeedRings2().

                                          {

  // TEC
  unsigned int counter = 0;

  unsigned int fw_bw_min       = 1;
  unsigned int fw_bw_max       = 1;
  unsigned int wheel_min       = 1;
  unsigned int wheel_max       = 9;
  unsigned int ring_min[9];
  unsigned int ring_max[9];

  if(structure_==FullDetector || structure_==FullDetectorII) {
    fw_bw_min = 1;
    fw_bw_max = 3;
    // TEC WHEEL 1
    ring_min[1] = 2;
    ring_max[1] = 3;
    // TEC WHEEL 2
    ring_min[2] = 2;
    ring_max[2] = 3;
    // TEC WHEEL 3
    ring_min[3] = 2;
    ring_max[3] = 3;
    // TEC WHEEL 4
    ring_min[4] = 2;
    ring_max[4] = 3;
    // TEC WHEEL 5
    ring_min[5] = 2;
    ring_max[5] = 3;
    // TEC WHEEL 6
    ring_min[6] = 2;
    ring_max[6] = 3;
    // TEC WHEEL 7
    ring_min[7] = 1;
    ring_max[7] = 1;
    // TEC WHEEL 8
    ring_min[8] = 1;
    ring_max[8] = 1;
  } else if (structure_==TIFTOBTEC || structure_==TIF) {  
    fw_bw_min = 2;
    fw_bw_max = 3;
    // TEC WHEEL 1
    ring_min[1] = 3;
    ring_max[1] = 4;
    // TEC WHEEL 2
    ring_min[2] = 3;
    ring_max[2] = 4;
    // TEC WHEEL 3
    ring_min[3] = 3;
    ring_max[3] = 4;
    // TEC WHEEL 4
    ring_min[4] = 3;
    ring_max[4] = 4;
    // TEC WHEEL 5
    ring_min[5] = 3;
    ring_max[5] = 4;
    // TEC WHEEL 6
    ring_min[6] = 3;
    ring_max[6] = 4; 
    // TEC WHEEL 6
    ring_min[7] = 3;
    ring_max[7] = 4;
    // TEC WHEEL 7
    ring_min[8] = 3;
    ring_max[8] = 4;
  }

  for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
    for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
      for ( unsigned int ring = ring_min[wheel]; ring < ring_max[wheel]; ++ring ) {
    const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,ring);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TEC inner seed ring with index: " << temp_ring->getindex(); 
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TEC inner seed rings"; 

}

void RoadMaker::collectInnerTIBSeedRings

(

)

Definition at line 292 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTIBRing(), innerSeedRings_, LogDebug, MTCC, rings_, structure_, TIF, TIFTIB, and TIFTIBTOB.

Referenced by collectInnerSeedRings().

                                         {

  // TIB
  unsigned int counter = 0, layer_min=0, layer_max=0, fw_bw_min=0, fw_bw_max=0, ext_int_min=0, ext_int_max=0, detector_min = 0, detector_max=0;
  if(structure_==FullDetector || structure_==FullDetectorII) {
    layer_min     = 1;
    layer_max     = 3;
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4;
  } else if(structure_==MTCC) {
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 2;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if (structure_==TIF) { 
    layer_min     = 1;
    layer_max     = 3; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if (structure_==TIFTIB) { 
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if (structure_==TIFTIBTOB) { 
    layer_min     = 1;
    layer_max     = 3; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  }


  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int ext_int = ext_int_min; ext_int < ext_int_max; ++ext_int ) {
    for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
      const Ring* temp_ring = rings_->getTIBRing(layer,fw_bw,ext_int,detector);
      innerSeedRings_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TIB inner seed ring with index: " << temp_ring->getindex(); 
    }    
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TIB inner seed rings"; 

}

void RoadMaker::collectInnerTIBSeedRings1

(

)

Definition at line 644 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTIBRing(), innerSeedRings1_, LogDebug, MTCC, rings_, structure_, TIF, TIFTIB, and TIFTIBTOB.

Referenced by collectInnerSeedRings1().

                                          {

  // TIB
  unsigned int counter = 0, layer_min=0, layer_max=0, fw_bw_min=0, fw_bw_max=0, ext_int_min=0, ext_int_max=0, detector_min=0, detector_max=0;
  if(structure_==FullDetector || structure_==FullDetectorII) {
    layer_min     = 1;
    layer_max     = 2;
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4;
  } else if(structure_==MTCC) {
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 2;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if (structure_==TIF) { 
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if (structure_==TIFTIB) { 
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if (structure_==TIFTIBTOB) { 
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  }

  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int ext_int = ext_int_min; ext_int < ext_int_max; ++ext_int ) {
    for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
      const Ring* temp_ring = rings_->getTIBRing(layer,fw_bw,ext_int,detector);
      innerSeedRings1_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TIB inner seed ring with index: " << temp_ring->getindex(); 
    }    
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TIB inner seed rings"; 

}

void RoadMaker::collectInnerTIBSeedRings2

(

)

Definition at line 966 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTIBRing(), innerSeedRings2_, LogDebug, MTCC, rings_, structure_, TIF, TIFTIB, TIFTIBTOB, and TIFTOBTEC.

Referenced by collectInnerSeedRings2().

                                          {

  // TIB
  unsigned int counter = 0, layer_min = 0, layer_max=0, fw_bw_min=0, fw_bw_max=0, ext_int_min=0, ext_int_max=0, detector_min=0, detector_max=0;
  if(structure_==FullDetector || structure_==FullDetectorII) {
    layer_min     = 2;
    layer_max     = 3;
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4;
  } else if(structure_==MTCC) {
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 2;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if(structure_==TIF) {
    layer_min     = 2;
    layer_max     = 3; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if(structure_==TIFTIB) {
    layer_min     = 1;
    layer_max     = 2; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if(structure_==TIFTIBTOB) {
    layer_min     = 2;
    layer_max     = 3; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  } else if(structure_==TIFTOBTEC) {
    layer_min     = 2;
    layer_max     = 3; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1;
    detector_max  = 4; 
  }

  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int ext_int = ext_int_min; ext_int < ext_int_max; ++ext_int ) {
    for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
      const Ring* temp_ring = rings_->getTIBRing(layer,fw_bw,ext_int,detector);
      innerSeedRings2_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TIB inner seed ring with index: " << temp_ring->getindex(); 
    }    
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TIB inner seed rings"; 

}

void RoadMaker::collectInnerTIDSeedRings

(

)

Definition at line 361 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTIDRing(), innerSeedRings_, LogDebug, relativeConstraints::ring, rings_, structure_, TIF, TIFTIB, and TIFTIBTOB.

Referenced by collectInnerSeedRings().

                                         {

  // TID
  unsigned int counter = 0;

  unsigned int fw_bw_min =0, fw_bw_max = 0,wheel_min = 0, wheel_max = 0, ring_min=0, ring_max = 0;

  if(structure_==FullDetector || structure_==FullDetectorII) {
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 3;
  } else if(structure_==TIF) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 3;
  } else if(structure_==TIFTIBTOB) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 3;
  } else if(structure_==TIFTIB) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 2;
  }

  for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
    for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
      for ( unsigned int ring = ring_min; ring < ring_max; ++ring ) {
    const Ring* temp_ring = rings_->getTIDRing(fw_bw,wheel,ring);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TID inner seed ring with index: " << temp_ring->getindex(); 
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TID inner seed rings"; 

}

void RoadMaker::collectInnerTIDSeedRings1

(

)

Definition at line 712 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTIDRing(), innerSeedRings1_, LogDebug, relativeConstraints::ring, rings_, structure_, TIF, TIFTIB, and TIFTIBTOB.

Referenced by collectInnerSeedRings1().

                                          {

  // TID
  unsigned int counter = 0;

  unsigned int fw_bw_min = 0, fw_bw_max = 0, wheel_min=0, wheel_max = 0, ring_min=0, ring_max = 0;

  if(structure_==FullDetector || structure_==FullDetectorII) {
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 2;
  } else if(structure_==TIF) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 2;
  } else if(structure_==TIFTIBTOB) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 2;
  } else if(structure_==TIFTIB) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 1;
    ring_max        = 2;
  }

  for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
    for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
      for ( unsigned int ring = ring_min; ring < ring_max; ++ring ) {
    const Ring* temp_ring = rings_->getTIDRing(fw_bw,wheel,ring);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TID inner seed ring with index: " << temp_ring->getindex(); 
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TID inner seed rings"; 

}

void RoadMaker::collectInnerTIDSeedRings2

(

)

Definition at line 1043 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTIDRing(), innerSeedRings2_, LogDebug, relativeConstraints::ring, rings_, structure_, TIF, TIFTIB, and TIFTIBTOB.

Referenced by collectInnerSeedRings2().

                                          {

  // TID
  unsigned int counter = 0;

  unsigned int fw_bw_min = 0, fw_bw_max = 0, wheel_min=0, wheel_max = 0, ring_min = 0, ring_max = 0;

  if(structure_==FullDetector) {
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 2;
    ring_max        = 3;
  } else if(structure_==FullDetectorII) {
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 2;
    ring_max        = 4;
  } else if(structure_==TIFTIBTOB || structure_==TIF) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 2;
    ring_max        = 3;
  } else if(structure_==TIFTIB) {
    fw_bw_min       = 2;
    fw_bw_max       = 3;
    wheel_min       = 1;
    wheel_max       = 4;
    ring_min        = 2;
    ring_max        = 2;
  }

  for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
    for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
      for ( unsigned int ring = ring_min; ring < ring_max; ++ring ) {
    const Ring* temp_ring = rings_->getTIDRing(fw_bw,wheel,ring);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TID inner seed ring with index: " << temp_ring->getindex(); 
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TID inner seed rings"; 

}

void RoadMaker::collectInnerTOBSeedRings

(

)

Definition at line 522 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTOBRing(), innerSeedRings_, LogDebug, MTCC, P5, rings_, structure_, TIF, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by collectInnerSeedRings().

                                         {

  // TOB
  unsigned int counter = 0, layer_min=0, layer_max=0, rod_fw_bw_min=0, rod_fw_bw_max=0, detector_min=0, detector_max=0;
  if(structure_==FullDetector || structure_==FullDetectorII) {
    layer_min       = 5;
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==MTCC) { 
    layer_min       = 1; 
    layer_max       = 3;
    rod_fw_bw_min   = 2;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOB || structure_==TIFTIBTOB || structure_==TIF) { 
    layer_min       = 1; 
    layer_max       = 3;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOBTEC) { 
    layer_min       = 1; 
    layer_max       = 3;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==P5) {
    layer_min       = 3; 
    layer_max       = 5;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;  
  }


  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int rod_fw_bw = rod_fw_bw_min; rod_fw_bw < rod_fw_bw_max; ++rod_fw_bw ) {
      for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
    const Ring* temp_ring = rings_->getTOBRing(layer,rod_fw_bw,detector);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
      }    
    }    
  }
  
  if(structure_==FullDetector || structure_==FullDetectorII) { 
    // add most outer rings
    const Ring* temp_ring = rings_->getTOBRing(2,1,6);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(2,2,6);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,1,6);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,2,6);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,1,6);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,2,6);
    innerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 

  }

  LogDebug("RoadSearch") << "collected " << counter << " TOB inner seed rings"; 

}

void RoadMaker::collectInnerTOBSeedRings1

(

)

Definition at line 847 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTOBRing(), innerSeedRings1_, LogDebug, MTCC, rings_, structure_, TIF, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by collectInnerSeedRings1().

                                          {

  // TOB
  unsigned int counter = 0, layer_min=0, layer_max=0, rod_fw_bw_min=0, rod_fw_bw_max=0, detector_min=0, detector_max=0;
  if(structure_==FullDetector || structure_==FullDetectorII) {
    layer_min       = 5;
    layer_max       = 6;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==MTCC) { 
    layer_min       = 1; 
    layer_max       = 3;
    rod_fw_bw_min   = 2;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOB || structure_==TIFTIBTOB || structure_==TIF) { 
    layer_min       = 1; 
    layer_max       = 2;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOBTEC) { 
    layer_min       = 1; 
    layer_max       = 2;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  }

  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int rod_fw_bw = rod_fw_bw_min; rod_fw_bw < rod_fw_bw_max; ++rod_fw_bw ) {
      for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
    const Ring* temp_ring = rings_->getTOBRing(layer,rod_fw_bw,detector);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
      }    
    }    
  }
  
  if(structure_==FullDetector) { 
    // add most outer rings
    const Ring* temp_ring = rings_->getTOBRing(2,1,6);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(2,2,6);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,1,6);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,2,6);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,1,6);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,2,6);
    innerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 

  }

  LogDebug("RoadSearch") << "collected " << counter << " TOB inner seed rings"; 

}

void RoadMaker::collectInnerTOBSeedRings2

(

)

Definition at line 1178 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTOBRing(), innerSeedRings2_, LogDebug, MTCC, rings_, structure_, TIF, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by collectInnerSeedRings2().

                                          {

  // TOB
  unsigned int counter = 0, layer_min=0, layer_max=0, rod_fw_bw_min=0, rod_fw_bw_max=0, detector_min=0, detector_max=0;
  if(structure_==FullDetector || structure_==FullDetectorII) {
    layer_min       = 5;
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_max    = 1;
    detector_max    = 7;
  } else if (structure_==MTCC) { 
    layer_min       = 1; 
    layer_max       = 3;
    rod_fw_bw_min   = 2;
    rod_fw_bw_max   = 3;
    detector_max    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOB || structure_==TIFTIBTOB || structure_==TIFTOBTEC || structure_==TIF) { 
    layer_min       = 2; 
    layer_max       = 3;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_max    = 1;
    detector_max    = 7;
  }

  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int rod_fw_bw = rod_fw_bw_min; rod_fw_bw < rod_fw_bw_max; ++rod_fw_bw ) {
      for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
    const Ring* temp_ring = rings_->getTOBRing(layer,rod_fw_bw,detector);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
      }    
    }    
  }
  
  if(structure_==FullDetector || structure_==FullDetectorII) { 
    // add most outer rings
    const Ring* temp_ring = rings_->getTOBRing(2,1,6);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(2,2,6);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,1,6);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,2,6);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,1,6);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,2,6);
    innerSeedRings2_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 

  }

  LogDebug("RoadSearch") << "collected " << counter << " TOB inner seed rings"; 

}

void RoadMaker::collectOuterSeedRings

(

)

Definition at line 1254 of file RoadMaker.cc.

References collectOuterTECSeedRings(), collectOuterTIBSeedRings(), collectOuterTOBSeedRings(), FullDetector, FullDetectorII, LogDebug, MTCC, outerSeedRings_, P5, structure_, TIF, TIFTIB, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by constructRoads().

                                      {
  
  if(structure_==FullDetector) {
    collectOuterTOBSeedRings();
    collectOuterTECSeedRings();
  } else if(structure_==FullDetectorII) {
    collectOuterTOBSeedRings();
    collectOuterTECSeedRings();
  } else if (structure_==MTCC) {
    collectOuterTOBSeedRings();
  } else if (structure_==TIF) {
    collectOuterTOBSeedRings();
    collectOuterTECSeedRings();
  } else if (structure_==TIFTOB) {
    collectOuterTOBSeedRings();
  } else if (structure_==TIFTIB) {
    collectOuterTIBSeedRings();
  } else if (structure_==TIFTIBTOB) {
    collectOuterTOBSeedRings();
  } else if (structure_==TIFTOBTEC) {
    collectOuterTOBSeedRings();
    collectOuterTECSeedRings();
  } else if (structure_==P5) {
    collectOuterTOBSeedRings();
    collectOuterTECSeedRings();
  }



  LogDebug("RoadSearch") << "collected " << outerSeedRings_.size() << " outer seed rings"; 
}

void RoadMaker::collectOuterSeedRings1

(

)

Definition at line 1620 of file RoadMaker.cc.

References collectOuterTECSeedRings1(), collectOuterTIBSeedRings1(), collectOuterTOBSeedRings1(), FullDetector, FullDetectorII, LogDebug, MTCC, outerSeedRings1_, structure_, TIF, TIFTIB, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by constructRoads().

                                       {
  
  if(structure_==FullDetector) {
    collectOuterTOBSeedRings1();
    collectOuterTECSeedRings1();
  } else if(structure_==FullDetectorII) {
    collectOuterTOBSeedRings1();
    collectOuterTECSeedRings1();
  } else if (structure_==MTCC) {
    collectOuterTOBSeedRings1();
  } else if (structure_==TIF) {
    collectOuterTOBSeedRings1();
    collectOuterTECSeedRings1();
  } else if (structure_==TIFTOB) {
    collectOuterTOBSeedRings1();
  } else if (structure_==TIFTIB) {
    collectOuterTIBSeedRings1();
  } else if (structure_==TIFTIBTOB) {
    collectOuterTOBSeedRings1();
  } else if (structure_==TIFTOBTEC) {
    collectOuterTOBSeedRings1();
    collectOuterTECSeedRings1();

  }

  LogDebug("RoadSearch") << "collected " << outerSeedRings1_.size() << " outer seed rings"; 
}

void RoadMaker::collectOuterTECSeedRings

(

)

Definition at line 1439 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTECRing(), LogDebug, outerSeedRings_, P5, relativeConstraints::ring, rings_, structure_, TIF, and TIFTOBTEC.

Referenced by collectOuterSeedRings().

                                         {

  // TEC
  unsigned int counter = 0;

  if(structure_==FullDetector) {
    // outer and second outer ring in all wheels except those treated in the following
    unsigned int fw_bw_min       = 1;
    unsigned int fw_bw_max       = 3;
    unsigned int wheel_min       = 1;
    unsigned int wheel_max       = 8;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
    const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    temp_ring = rings_->getTECRing(fw_bw,wheel,7);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
      }
    }
    
    // add last two wheels
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    unsigned int wheel_start     = 8;
    wheel_max       = 10;
    unsigned int second_ring_min[10];
    unsigned int second_ring_max[10];
    
    // TEC WHEEL 1
    second_ring_min[1] = 1;
    second_ring_max[1] = 8;
    // TEC WHEEL 2
    second_ring_min[2] = 1;
    second_ring_max[2] = 8;
    // TEC WHEEL 3
    second_ring_min[3] = 1;
    second_ring_max[3] = 8;
    // TEC WHEEL 4
    second_ring_min[4] = 2;
    second_ring_max[4] = 8;
    // TEC WHEEL 5
    second_ring_min[5] = 2;
    second_ring_max[5] = 8;
    // TEC WHEEL 6
    second_ring_min[6] = 2;
    second_ring_max[6] = 8;
    // TEC WHEEL 7
    second_ring_min[7] = 3;
    second_ring_max[7] = 8;
    // TEC WHEEL 8
    second_ring_min[8] = 3;
    second_ring_max[8] = 8;
    // TEC WHEEL 9
    second_ring_min[9] = 4;
    second_ring_max[9] = 8;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_start; wheel < wheel_max; ++wheel ) {
    for ( unsigned int second_ring = second_ring_min[wheel]; second_ring < second_ring_max[wheel]; ++second_ring ) {
      const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,second_ring);
      outerSeedRings_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
    }    
      }
    }
  } else if(structure_==FullDetectorII) {
    // outer and second outer ring in all wheels except those treated in the following
    unsigned int fw_bw_min       = 1;
    unsigned int fw_bw_max       = 3;
    unsigned int wheel_min       = 1;
    unsigned int wheel_max       = 8;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
    const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,5);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    temp_ring = rings_->getTECRing(fw_bw,wheel,4);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
      }
    }
    
    // add last two wheels
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    unsigned int wheel_start     = 8;
    wheel_max       = 10;
    unsigned int second_ring_min[10];
    unsigned int second_ring_max[10];
    
    // TEC WHEEL 1
    second_ring_min[1] = 1;
    second_ring_max[1] = 8;
    // TEC WHEEL 2
    second_ring_min[2] = 1;
    second_ring_max[2] = 8;
    // TEC WHEEL 3
    second_ring_min[3] = 1;
    second_ring_max[3] = 8;
    // TEC WHEEL 4
    second_ring_min[4] = 2;
    second_ring_max[4] = 8;
    // TEC WHEEL 5
    second_ring_min[5] = 2;
    second_ring_max[5] = 8;
    // TEC WHEEL 6
    second_ring_min[6] = 2;
    second_ring_max[6] = 8;
    // TEC WHEEL 7
    second_ring_min[7] = 3;
    second_ring_max[7] = 8;
    // TEC WHEEL 8
    second_ring_min[8] = 3;
    second_ring_max[8] = 6;
    // TEC WHEEL 9
    second_ring_min[9] = 4;
    second_ring_max[9] = 6;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_start; wheel < wheel_max; ++wheel ) {
    for ( unsigned int second_ring = second_ring_min[wheel]; second_ring < second_ring_max[wheel]; ++second_ring ) {
      const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,second_ring);
      outerSeedRings_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
    }    
      }
    }
  } else if (structure_==TIFTOBTEC || structure_==TIF) {
    unsigned int fw_bw_min = 2;
    unsigned int fw_bw_max = 3;
    unsigned int wheel_min = 1;
    unsigned int wheel_max = 10;
    unsigned int ring_min  = 6;
    unsigned int ring_max  = 8;

    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
    for ( unsigned int ring = ring_min; ring < ring_max; ++ring ) {
      const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,ring);
      outerSeedRings_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
    }
      }
    }
    
  } else if (structure_==P5) {
    unsigned int fw_bw_min = 1;
    unsigned int fw_bw_max = 3;
    unsigned int wheel_min = 1;
    unsigned int wheel_max = 10;
    unsigned int ring_min  = 6;
    unsigned int ring_max  = 8;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
    for ( unsigned int ring = ring_min; ring < ring_max; ++ring ) {
      const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,ring);
      outerSeedRings_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
    }
      }
    }
    
  }

  LogDebug("RoadSearch") << "collected " << counter << " TEC outer seed rings"; 
  
}

void RoadMaker::collectOuterTECSeedRings1

(

)

Definition at line 1805 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTECRing(), LogDebug, outerSeedRings1_, relativeConstraints::ring, rings_, structure_, TIF, and TIFTOBTEC.

Referenced by collectOuterSeedRings1().

                                          {

  // TEC
  unsigned int counter = 0;

  if(structure_==FullDetector) {
    // outer ring in all wheels except those treated in the following
    unsigned int fw_bw_min       = 1;
    unsigned int fw_bw_max       = 3;
    unsigned int wheel_min       = 1;
    unsigned int wheel_max       = 9;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
    const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,7);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
      }
    }
    
    // add last two wheels
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    unsigned int wheel_start     = 9;
    wheel_max       = 10;
    unsigned int second_ring_min[10];
    unsigned int second_ring_max[10];
    
    // TEC WHEEL 1
    second_ring_min[1] = 1;
    second_ring_max[1] = 8;
    // TEC WHEEL 2
    second_ring_min[2] = 1;
    second_ring_max[2] = 8;
    // TEC WHEEL 3
    second_ring_min[3] = 1;
    second_ring_max[3] = 8;
    // TEC WHEEL 4
    second_ring_min[4] = 2;
    second_ring_max[4] = 8;
    // TEC WHEEL 5
    second_ring_min[5] = 2;
    second_ring_max[5] = 8;
    // TEC WHEEL 6
    second_ring_min[6] = 2;
    second_ring_max[6] = 8;
    // TEC WHEEL 7
    second_ring_min[7] = 3;
    second_ring_max[7] = 8;
    // TEC WHEEL 8
    second_ring_min[8] = 3;
    second_ring_max[8] = 8;
    // TEC WHEEL 9
    second_ring_min[9] = 4;
    second_ring_max[9] = 8;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_start; wheel < wheel_max; ++wheel ) {
    for ( unsigned int second_ring = second_ring_min[wheel]; second_ring < second_ring_max[wheel]; ++second_ring ) {
      const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,second_ring);
      outerSeedRings1_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
    }    
      }
    }
  } else if(structure_==FullDetectorII) {
    // outer ring in all wheels except those treated in the following
    unsigned int fw_bw_min       = 1;
    unsigned int fw_bw_max       = 3;
    unsigned int wheel_min       = 1;
    unsigned int wheel_max       = 9;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
    const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,5);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
      }
    }
    
    // add last two wheels
    fw_bw_min       = 1;
    fw_bw_max       = 3;
    wheel_min       = 1;
    unsigned int wheel_start     = 9;
    wheel_max       = 10;
    unsigned int second_ring_min[10];
    unsigned int second_ring_max[10];
    
    // TEC WHEEL 1
    second_ring_min[1] = 1;
    second_ring_max[1] = 8;
    // TEC WHEEL 2
    second_ring_min[2] = 1;
    second_ring_max[2] = 8;
    // TEC WHEEL 3
    second_ring_min[3] = 1;
    second_ring_max[3] = 8;
    // TEC WHEEL 4
    second_ring_min[4] = 2;
    second_ring_max[4] = 8;
    // TEC WHEEL 5
    second_ring_min[5] = 2;
    second_ring_max[5] = 8;
    // TEC WHEEL 6
    second_ring_min[6] = 2;
    second_ring_max[6] = 8;
    // TEC WHEEL 7
    second_ring_min[7] = 3;
    second_ring_max[7] = 8;
    // TEC WHEEL 8
    second_ring_min[8] = 3;
    second_ring_max[8] = 8;
    // TEC WHEEL 9
    second_ring_min[9] = 4;
    second_ring_max[9] = 6;
    
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_start; wheel < wheel_max; ++wheel ) {
    for ( unsigned int second_ring = second_ring_min[wheel]; second_ring < second_ring_max[wheel]; ++second_ring ) {
      const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,second_ring);
      outerSeedRings1_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
    }    
      }
    }
  } else if (structure_==TIFTOBTEC || structure_==TIF) {
    unsigned int fw_bw_min = 2;
    unsigned int fw_bw_max = 3;
    unsigned int wheel_min = 1;
    unsigned int wheel_max = 10;
    unsigned int ring_min  = 7;
    unsigned int ring_max  = 8;

    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int wheel = wheel_min; wheel < wheel_max; ++wheel ) {
    for ( unsigned int ring = ring_min; ring < ring_max; ++ring ) {
      const Ring* temp_ring = rings_->getTECRing(fw_bw,wheel,ring);
      outerSeedRings1_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TEC outer seed ring with index: " << temp_ring->getindex(); 
    }
      }
    }
    
  }
  LogDebug("RoadSearch") << "collected " << counter << " TEC outer seed rings"; 

}

void RoadMaker::collectOuterTIBSeedRings

(

)

Definition at line 1286 of file RoadMaker.cc.

References Ring::getindex(), Rings::getTIBRing(), LogDebug, outerSeedRings_, rings_, structure_, and TIFTIB.

Referenced by collectOuterSeedRings().

                                         {

  // TIB
  unsigned int counter = 0, layer_min=0, layer_max=0, fw_bw_min=0, fw_bw_max=0, ext_int_min=0, ext_int_max=0, detector_min=0, detector_max=0;
  if(structure_==TIFTIB) {
    layer_min     = 3; 
    layer_max     = 5; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1; 
    detector_max  = 4; 
  } 

  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int ext_int = ext_int_min; ext_int < ext_int_max; ++ext_int ) {
    for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
      const Ring* temp_ring = rings_->getTIBRing(layer,fw_bw,ext_int,detector);
      outerSeedRings_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TIB outer seed ring with index: " << temp_ring->getindex(); 
    }    
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TIB outer seed rings"; 

}

void RoadMaker::collectOuterTIBSeedRings1

(

)

Definition at line 1648 of file RoadMaker.cc.

References Ring::getindex(), Rings::getTIBRing(), LogDebug, outerSeedRings1_, rings_, structure_, and TIFTIB.

Referenced by collectOuterSeedRings1().

                                          {

  // TIB
  unsigned int counter = 0, layer_min=0, layer_max=0, fw_bw_min=0, fw_bw_max=0, ext_int_min=0, ext_int_max=0, detector_min=0, detector_max=0;
  if(structure_==TIFTIB) {
    layer_min     = 4; 
    layer_max     = 5; 
    fw_bw_min     = 1;
    fw_bw_max     = 3;
    ext_int_min   = 1;
    ext_int_max   = 3;
    detector_min  = 1; 
    detector_max  = 4; 
  } 

  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int fw_bw = fw_bw_min; fw_bw < fw_bw_max; ++fw_bw ) {
      for ( unsigned int ext_int = ext_int_min; ext_int < ext_int_max; ++ext_int ) {
    for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
      const Ring* temp_ring = rings_->getTIBRing(layer,fw_bw,ext_int,detector);
      outerSeedRings1_.push_back(temp_ring);
      ++counter;
      LogDebug("RoadSearch") << "collected TIB outer seed ring with index: " << temp_ring->getindex(); 
    }    
      }    
    }
  }

  LogDebug("RoadSearch") << "collected " << counter << " TIB outer seed rings"; 

}

void RoadMaker::collectOuterTOBSeedRings

(

)

Definition at line 1318 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTOBRing(), LogDebug, MTCC, outerSeedRings_, P5, rings_, structure_, TIF, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by collectOuterSeedRings().

                                         {

  // TOB
  unsigned int counter = 0, layer_min=0, layer_max=0, rod_fw_bw_min=0, rod_fw_bw_max=0, detector_min=0, detector_max=0;
  if(structure_==FullDetector) {
    layer_min       = 5;
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if(structure_==FullDetectorII) {
    layer_min       = 1;
    layer_max       = 2;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==MTCC) { 
    layer_min       = 1; 
    layer_max       = 3;
    rod_fw_bw_min   = 2;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIF) { 
    layer_min       = 5; 
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOB) { 
    layer_min       = 5; 
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTIBTOB) { 
    layer_min       = 5; 
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOBTEC) { 
    layer_min       = 5; 
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==P5) { 
    layer_min       = 5; 
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  }

  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int rod_fw_bw = rod_fw_bw_min; rod_fw_bw < rod_fw_bw_max; ++rod_fw_bw ) {
      for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
    const Ring* temp_ring = rings_->getTOBRing(layer,rod_fw_bw,detector);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
      }    
    }    
  }
  
  if(structure_==FullDetector) { 
    // add most outer rings
    const Ring* temp_ring = rings_->getTOBRing(1,1,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(1,2,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(2,1,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(2,2,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,1,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,2,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,1,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,2,6);
    outerSeedRings_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 

  }

  LogDebug("RoadSearch") << "collected " << counter << " TOB outer seed rings"; 

}

void RoadMaker::collectOuterTOBSeedRings1

(

)

Definition at line 1680 of file RoadMaker.cc.

References FullDetector, FullDetectorII, Ring::getindex(), Rings::getTOBRing(), LogDebug, MTCC, outerSeedRings1_, rings_, structure_, TIF, TIFTIBTOB, TIFTOB, and TIFTOBTEC.

Referenced by collectOuterSeedRings1().

                                          {

  // TOB
  unsigned int counter = 0, layer_min=0, layer_max=0, rod_fw_bw_min=0, rod_fw_bw_max=0, detector_min=0, detector_max=0;
  if(structure_==FullDetector) {
    layer_min       = 6;
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if(structure_==FullDetectorII) {
    layer_min       = 2;
    layer_max       = 3;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==MTCC) { 
    layer_min       = 1; 
    layer_max       = 3;
    rod_fw_bw_min   = 2;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIF) { 
    layer_min       = 6;
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
   } else if (structure_==TIFTOB) { 
    layer_min       = 5;
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTIBTOB) { 
    layer_min       = 6;
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  } else if (structure_==TIFTOBTEC) { 
    layer_min       = 6; 
    layer_max       = 7;
    rod_fw_bw_min   = 1;
    rod_fw_bw_max   = 3;
    detector_min    = 1;
    detector_max    = 7;
  }


  for ( unsigned int layer = layer_min; layer < layer_max; ++layer ) {
    for ( unsigned int rod_fw_bw = rod_fw_bw_min; rod_fw_bw < rod_fw_bw_max; ++rod_fw_bw ) {
      for ( unsigned int detector = detector_min; detector < detector_max; ++detector ) {
    const Ring* temp_ring = rings_->getTOBRing(layer,rod_fw_bw,detector);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
      }    
    }    
  }
  
  if(structure_==FullDetector) { 
    // add most outer rings
    const Ring* temp_ring = rings_->getTOBRing(1,1,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(1,2,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(2,1,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(2,2,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,1,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(3,2,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,1,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
    
    temp_ring = rings_->getTOBRing(4,2,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 

    temp_ring = rings_->getTOBRing(5,1,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 
      
    temp_ring = rings_->getTOBRing(5,2,6);
    outerSeedRings1_.push_back(temp_ring);
    ++counter;
    LogDebug("RoadSearch") << "collected TOB outer seed ring with index: " << temp_ring->getindex(); 

  }

  LogDebug("RoadSearch") << "collected " << counter << " TOB outer seed rings"; 

}

bool RoadMaker::CompatibleWithLines	(	std::vector< std::pair< double, double > >	lines,
		const Ring *	ring
	)

Definition at line 2129 of file RoadMaker.cc.

References Ring::getrmax(), Ring::getrmin(), Ring::getzmax(), Ring::getzmin(), geometryCSVtoXML::line, groupFilesInBlocks::lines, m, alignCSCRings::r, query::result, and detailsBasic3DVector::z.

Referenced by constructRoads(), and RingsCompatibleWithSeed().

                                      {
  //
  // check compatibility of ring z extension with extrapolation of lines to ring radius
  //

  // return value
  bool result = true;

  // calculate zmin, zmax at the radius of the ring
  double zmin = 999;
  double zmax = -zmin;
  for (int m=0; m<2; ++m) {
    double r = m==0 ? ring->getrmin() : ring->getrmax();
    for (std::vector<std::pair<double,double> >::iterator line = lines.begin();
     line != lines.end();
     ++line ) {
      double z = line->first*r + line->second;
      if (zmin>z) zmin = z;
      if (zmax<z) zmax = z;
    }
  }
  if (ring->getzmax()<zmin || ring->getzmin()>zmax) {
    result = false;
  }

  return result;
}

void RoadMaker::constructRoads

(

)

Definition at line 80 of file RoadMaker.cc.

References abs, collectInnerSeedRings(), collectInnerSeedRings1(), collectInnerSeedRings2(), collectOuterSeedRings(), collectOuterSeedRings1(), CompatibleWithLines(), diffTreeTool::diff, FourRingSeeds, innerSeedRings1_, innerSeedRings2_, innerSeedRings_, Roads::insert(), LinesThroughRingAndBS(), LinesThroughRings(), outerSeedRings1_, outerSeedRings_, alignCSCRings::r, relativeConstraints::ring, RingInBarrel(), RingsCompatibleWithSeed(), RingsOnSameLayer(), roads_, seedingType_, runtimedef::set(), Roads::size(), findQualityFiles::size, slope, python.multivaluedict::sort(), TwoRingSeeds, and detailsBasic3DVector::z.

Referenced by RoadMaker().

                               {

  // fill vector of inner Rings
  collectInnerSeedRings();
  collectInnerSeedRings1();
  collectInnerSeedRings2();

  // fill vector of outer Rings
  collectOuterSeedRings();
  collectOuterSeedRings1();

  if ( seedingType_ == TwoRingSeeds ) {
    // loop over inner-outer ring pairs
    for ( std::vector<const Ring*>::iterator innerRing = innerSeedRings_.begin(); 
      innerRing != innerSeedRings_.end(); 
      ++innerRing ) {
      // draw lines (z = a*r + b) through all corners and beam spot
      std::vector<std::pair<double,double> > linesInnerSeedRingAndBS = LinesThroughRingAndBS((*innerRing));
      for ( std::vector<const Ring*>::iterator outerRing = outerSeedRings_.begin(); 
        outerRing != outerSeedRings_.end(); 
        ++outerRing ) {
    // check if ring is compatible with extrapolation of lines
    if ( CompatibleWithLines(linesInnerSeedRingAndBS, (*outerRing) ) ) {
      // construct road
      Roads::RoadSeed seed;

      // fill first inner seed ring
      seed.first.push_back((*innerRing));

      // fill first outer seed ring
      seed.second.push_back((*outerRing));

      Roads::RoadSet set = RingsCompatibleWithSeed(seed);

      // sort seeds
      std::sort(seed.first.begin(),seed.first.end(),SortRingsByZR());
      std::sort(seed.second.begin(),seed.second.end(),SortRingsByZR());
        
      roads_->insert(seed,set);
    }
      }
    }
    
  } else if (seedingType_ == FourRingSeeds ) {

    // loop over inner-outer ring pairs
    for ( std::vector<const Ring*>::iterator innerRing1 = innerSeedRings1_.begin(); 
      innerRing1 != innerSeedRings1_.end(); 
      ++innerRing1 ) {
      
      // draw lines (z = a*r + b) through all corners and beam spot
      std::vector<std::pair<double,double> > linesInnerSeedRing1AndBS = LinesThroughRingAndBS((*innerRing1));
      
      for ( std::vector<const Ring*>::iterator innerRing2 = innerSeedRings2_.begin(); 
        innerRing2 != innerSeedRings2_.end(); 
        ++innerRing2 ) {
    
    // check if ring is compatible with extrapolation of lines
    if ( CompatibleWithLines(linesInnerSeedRing1AndBS, (*innerRing2) ) ) {

      // draw lines (z = a*r + b) through all corners and beam spot
      std::vector<std::pair<double,double> > linesInnerSeedRing2AndBS = LinesThroughRingAndBS((*innerRing2));

      // draw lines (z = a*r + b) through all corners of both inner seed rings
      std::vector<std::pair<double,double> > linesInnerSeedRing1AndInnerSeedRing2 = LinesThroughRings(*innerRing1,*innerRing2);

      for ( std::vector<const Ring*>::iterator outerRing1 = outerSeedRings1_.begin(); 
        outerRing1 != outerSeedRings1_.end(); 
        ++outerRing1 ) {

        // check if ring is compatible with extrapolation of lines
        if ( CompatibleWithLines(linesInnerSeedRing1AndBS, (*outerRing1) ) &&
         CompatibleWithLines(linesInnerSeedRing2AndBS, (*outerRing1) ) &&
         CompatibleWithLines(linesInnerSeedRing1AndInnerSeedRing2, (*outerRing1) )  ) {

          std::vector<std::pair<double,double> > linesOuterSeedRing1AndBS = LinesThroughRingAndBS((*outerRing1));
          std::vector<std::pair<double,double> > linesInnerSeedRing1AndOuterSeedRing1 = LinesThroughRings(*innerRing1,*outerRing1);
          std::vector<std::pair<double,double> > linesInnerSeedRing2AndOuterSeedRing1 = LinesThroughRings(*innerRing2,*outerRing1);

          // construct road
          Roads::RoadSeed seed;

          // fill first inner seed ring
          seed.first.push_back((*innerRing1));
          seed.first.push_back((*innerRing2));

          // fill first outer seed ring
          seed.second.push_back((*outerRing1));


          std::vector<const Ring*> outerRing2Candidates;
          // try to add second outer seed ring from all outer seed rings if compatible with other three seeds
          for ( std::vector<const Ring*>::iterator outerRing = outerSeedRings_.begin(); 
            outerRing != outerSeedRings_.end(); 
            ++outerRing ) {
        if ( CompatibleWithLines(linesInnerSeedRing1AndBS, (*outerRing) ) &&
             CompatibleWithLines(linesInnerSeedRing2AndBS, (*outerRing) ) &&
             CompatibleWithLines(linesOuterSeedRing1AndBS, (*outerRing) ) &&
             CompatibleWithLines(linesInnerSeedRing1AndInnerSeedRing2, (*outerRing) ) &&
             CompatibleWithLines(linesInnerSeedRing1AndOuterSeedRing1, (*outerRing) ) &&
             CompatibleWithLines(linesInnerSeedRing2AndOuterSeedRing1, (*outerRing) ) ) {
          if ( (*outerRing)->getindex() < (*outerRing1)->getindex() ) {
            if ( !RingsOnSameLayer((*outerRing),(*outerRing1)) ) {
              outerRing2Candidates.push_back((*outerRing));
            }
          }
        }
          }

          unsigned int size = outerRing2Candidates.size();
          if ( size > 0 ) {
        if ( size == 1 ) {
          seed.second.push_back(outerRing2Candidates[0]);
        } else {
          // extrapolate center of outerRing1 throigh 0,0 to candidates
          // candidate with center clostest to extrapolation willl be taken
          const Ring *selectedRing = 0;
          double z_2 = ((*outerRing1)->getzmax()+(*outerRing1)->getzmin())/2;
          double r_2 = ((*outerRing1)->getrmax()+(*outerRing1)->getrmin())/2;
          double z_1 = ((*innerRing1)->getzmax()+(*innerRing1)->getzmin())/2;
          double r_1 = ((*innerRing1)->getrmax()+(*innerRing1)->getrmin())/2;
          if ( RingInBarrel(*outerRing1) ) {
            double slope = (z_2-z_1) /(r_2-r_1);
            double intercept = z_1 - slope * r_1;
            double zDifference = 999;
            for ( std::vector<const Ring*>::iterator ring = outerRing2Candidates.begin();
              ring != outerRing2Candidates.end();
              ++ring ) {
              double z = slope * (((*ring)->getrmax()+(*ring)->getrmin())/2) + intercept;
              double diff = std::abs(z-(((*ring)->getzmax()+(*ring)->getzmin())/2));
              if ( diff < zDifference ) {
            selectedRing = *ring;
            zDifference = diff;
              }
            }
          } else {
            double slope = (r_2-r_1) /(z_2-z_1);
            double intercept = r_1 - slope * z_1;
            double rDifference = 999;
            for ( std::vector<const Ring*>::iterator ring = outerRing2Candidates.begin();
              ring != outerRing2Candidates.end();
              ++ring ) {
              double r = slope * (((*ring)->getzmax()+(*ring)->getzmin())/2) + intercept;
              double diff = std::abs(r-(((*ring)->getrmax()+(*ring)->getrmin())/2));
              if ( diff < rDifference ) {
            selectedRing = *ring;
            rDifference = diff;
              }
            }
          }
          seed.second.push_back(selectedRing);
        }
          }

          
          // collect all rings compatible with temporary seed from inner seed ring 1 and outer seed ring 1
          Roads::RoadSeed tempSeed;
          tempSeed.first.push_back((*innerRing1));
          tempSeed.second.push_back((*outerRing1));
          Roads::RoadSet set = RingsCompatibleWithSeed(tempSeed);

          // sort seeds
          std::sort(seed.first.begin(),seed.first.end(),SortRingsByZR());
          std::sort(seed.second.begin(),seed.second.end(),SortRingsByZR());
        
          roads_->insert(seed,set);

        }
      }
    }
      }
    }
  }
  
  edm::LogInfo("RoadSearch") << "Constructed " << roads_->size() << " roads.";

}

Roads* RoadMaker::getRoads ( )

inline

Definition at line 86 of file RoadMaker.h.

References roads_.

Referenced by RoadMapMakerESProducer::produce().

{ return roads_; }

std::vector< std::pair< double, double > > RoadMaker::LinesThroughRingAndBS

(

const Ring *

ring

)

Definition at line 2067 of file RoadMaker.cc.

References Ring::getrmax(), Ring::getrmin(), Ring::getzmax(), Ring::getzmin(), alignCSCRings::r, query::result, detailsBasic3DVector::z, and zBS_.

Referenced by constructRoads(), and RingsCompatibleWithSeed().

                                                                                     {
  //
  // calculate lines through all 4 corners of the ring and the BS
  //
   
  // return value
  std::vector<std::pair<double,double> > result;
   
  double z = 0;
  double r = 0;;
  
  for (int ic1 = 0; ic1<4; ++ic1) {
    switch (ic1) {
    case 0: z = ring->getzmin(); r = ring->getrmin(); break;
    case 1: z = ring->getzmin(); r = ring->getrmax(); break;
    case 2: z = ring->getzmax(); r = ring->getrmax(); break;
    case 3: z = ring->getzmax(); r = ring->getrmin(); break;
    }
    for (int ib = 0; ib<2; ++ib) {
      double zb = zBS_*(2*ib-1);
      result.push_back(std::pair<double,double>((z-zb)/r,zb));
    }
  }

  return result;
}

std::vector< std::pair< double, double > > RoadMaker::LinesThroughRings	(	const Ring *	ring1,
		const Ring *	ring2
	)

Definition at line 2094 of file RoadMaker.cc.

References Ring::getrmax(), Ring::getrmin(), Ring::getzmax(), Ring::getzmin(), diffTwoXMLs::r1, diffTwoXMLs::r2, and query::result.

Referenced by constructRoads(), and RingsCompatibleWithSeed().

                                                       {
  //
  // calculate lines through all 4 corners of the rings
  //

  // return value
  std::vector<std::pair<double,double> > result;

  double z1 = 0;
  double r1 = 0;;
  double z2 = 0;
  double r2 = 0;

  for (int ic1 = 0; ic1<4; ++ic1) {
    switch (ic1) {
    case 0: z1 = ring1->getzmin(); r1 = ring1->getrmin(); break;
    case 1: z1 = ring1->getzmin(); r1 = ring1->getrmax(); break;
    case 2: z1 = ring1->getzmax(); r1 = ring1->getrmax(); break;
    case 3: z1 = ring1->getzmax(); r1 = ring1->getrmin(); break;
    }
    for (int ic2 = 0; ic2<4; ++ic2) {
      switch (ic2) {
      case 0: z2 = ring2->getzmin(); r2 = ring2->getrmin(); break;
      case 1: z2 = ring2->getzmin(); r2 = ring2->getrmax(); break;
      case 2: z2 = ring2->getzmax(); r2 = ring2->getrmax(); break;
      case 3: z2 = ring2->getzmax(); r2 = ring2->getrmin(); break;
      }
      result.push_back(std::pair<double,double>((z2 - z1)/(r2 - r1),(r2*z1 - r1*z2)/(r2 - r1)));
    }
  }
  
  return result;
}

bool RoadMaker::RingInBarrel

(

const Ring *

ring

)

Definition at line 2032 of file RoadMaker.cc.

References Ring::getFirst(), query::result, DetId::subdetId(), StripSubdetector::TIB, and StripSubdetector::TOB.

Referenced by constructRoads().

                                             {
  //
  // check if the ring is a TIB or TOB ring
  //

  // return value
  bool result = false;
  
  // get first DetId of ring
  const DetId ringDetId = ring->getFirst();

  // check if both rings belong to same subdetector
  if ( (unsigned int)ringDetId.subdetId() == StripSubdetector::TIB ) {
    result = true;
  } else if ( (unsigned int)ringDetId.subdetId() == StripSubdetector::TOB ) {
    result = true;
  }
  
  return result;
}

Roads::RoadSet RoadMaker::RingsCompatibleWithSeed

(

Roads::RoadSeed

seed

)

Definition at line 2158 of file RoadMaker.cc.

References Rings::begin(), CompatibleWithLines(), Rings::end(), geometryCSVtoXML::line, groupFilesInBlocks::lines, LinesThroughRingAndBS(), LinesThroughRings(), relativeConstraints::ring, rings_, and SortRingsIntoLayers().

Referenced by AddInnerSeedRing(), AddOuterSeedRing(), and constructRoads().

                                                                  {
  //
  // collect all rings which are compatible with the seed
  //
  
  // return value
  std::vector<const Ring*> tempRings;

  // calculate lines
  std::vector<std::vector<std::pair<double,double> > > lines;

  for ( std::vector<const Ring*>::iterator innerRing = seed.first.begin();
    innerRing != seed.first.end();
    ++innerRing) {
    // calculate lines between inner seed rings and BS
    lines.push_back(LinesThroughRingAndBS(*innerRing));

    for ( std::vector<const Ring*>::iterator outerRing = seed.second.begin();
      outerRing != seed.second.end();
      ++outerRing) {
      // calculate lines between inner and outer seed rings
      lines.push_back(LinesThroughRings((*innerRing),(*outerRing)));
    }
  }

  for ( std::vector<const Ring*>::iterator outerRing = seed.second.begin();
    outerRing != seed.second.end();
    ++outerRing) {
    // calculate lines between outer seed rings and BS
    lines.push_back(LinesThroughRingAndBS(*outerRing));
  }

  for ( Rings::const_iterator ring = rings_->begin();
    ring != rings_->end();
    ++ring ) {
    bool compatible = true;
    for ( std::vector<std::vector<std::pair<double,double> > >::iterator line = lines.begin();
      line != lines.end();
      ++line ) {
      if ( !CompatibleWithLines(*line, &(ring->second))) {
    compatible = false;
      }
    }
    if ( compatible ) {
      tempRings.push_back(&(ring->second));
    }
  }

  

  return SortRingsIntoLayers(tempRings);
}

bool RoadMaker::RingsOnSameLayer	(	const Ring *	ring1,
		const Ring *	ring2
	)

Definition at line 1960 of file RoadMaker.cc.

References Ring::getFirst(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), query::result, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, and StripSubdetector::TOB.

Referenced by AddInnerSeedRing(), AddOuterSeedRing(), constructRoads(), RingsOnSameLayer(), and SortRingsIntoLayers().

                                                                     {
  //
  // check whether two input rings are on the same layer
  //

  // return value
  bool result = false;
  
  // get first DetId of ring
  const DetId ring1DetId = ring1->getFirst();
  const DetId ring2DetId = ring2->getFirst();

  // check if both rings belong to same subdetector
  if ( (unsigned int)ring1DetId.subdetId() == StripSubdetector::TIB && 
       (unsigned int)ring2DetId.subdetId() == StripSubdetector::TIB ) {
    // make TIBDetId instance
    TIBDetId ring1DetIdTIB(ring1DetId.rawId());
    TIBDetId ring2DetIdTIB(ring2DetId.rawId());
    // check whether both rings are on the same TIB layer
    if ( ring1DetIdTIB.layer() == ring2DetIdTIB.layer() ) {
      result = true;
    }
  } else if ( (unsigned int)ring1DetId.subdetId() == StripSubdetector::TOB &&
          (unsigned int)ring2DetId.subdetId() == StripSubdetector::TOB ) {
    // make TOBDetId instance
    TOBDetId ring1DetIdTOB(ring1DetId.rawId());
    TOBDetId ring2DetIdTOB(ring2DetId.rawId());
    // check whether both rings are on the same TOB layer
    if ( ring1DetIdTOB.layer() == ring2DetIdTOB.layer() ) {
      result = true;
    }
  } else if ( (unsigned int)ring1DetId.subdetId() == StripSubdetector::TID && 
          (unsigned int)ring2DetId.subdetId() == StripSubdetector::TID) {
    // make TIDDetId instance
    TIDDetId ring1DetIdTID(ring1DetId.rawId());
    TIDDetId ring2DetIdTID(ring2DetId.rawId());
    // check whether both rings are on the same TID wheel
    if ( ring1DetIdTID.wheel() == ring2DetIdTID.wheel() ) {
      result = true;
    }
  } else if ( (unsigned int)ring1DetId.subdetId() == StripSubdetector::TEC &&
          (unsigned int)ring2DetId.subdetId() == StripSubdetector::TEC ) {
    // make TECDetId instance
    TECDetId ring1DetIdTEC(ring1DetId.rawId());
    TECDetId ring2DetIdTEC(ring2DetId.rawId());
    // check whether both rings are on the same TEC wheel
    if ( ring1DetIdTEC.wheel() == ring2DetIdTEC.wheel() ) {
      result = true;
    }
  } else if ( (unsigned int)ring1DetId.subdetId() == PixelSubdetector::PixelBarrel && 
          (unsigned int)ring2DetId.subdetId() == PixelSubdetector::PixelBarrel) {
    // make PXBDetId instance
    PXBDetId ring1DetIdPXB(ring1DetId.rawId());
    PXBDetId ring2DetIdPXB(ring2DetId.rawId());
    // check whether both rings are on the same PXB layer
    if ( ring1DetIdPXB.layer() == ring2DetIdPXB.layer() ) {
      result = true;
    }
  } else if ( (unsigned int)ring1DetId.subdetId() == PixelSubdetector::PixelEndcap &&
          (unsigned int)ring2DetId.subdetId() == PixelSubdetector::PixelEndcap) {
    // make PXFDetId instance
    PXFDetId ring1DetIdPXF(ring1DetId.rawId());
    PXFDetId ring2DetIdPXF(ring2DetId.rawId());
    // check whether both rings are on the same PXF disk
    if ( ring1DetIdPXF.disk() == ring2DetIdPXF.disk() ) {
      result = true;
    }
  }
  
  return result;
}

bool RoadMaker::RingsOnSameLayer	(	std::pair< const Ring *, const Ring * >	seed1,
		std::pair< const Ring *, const Ring * >	seed2
	)

Definition at line 2053 of file RoadMaker.cc.

References query::result, and RingsOnSameLayer().

                                                                                                               {
  //
  // check whether two input seeds (two ring seeds) are on the same layer
  //
  
  // return value
  bool result = false;
  
  result = RingsOnSameLayer(seed1.first,seed2.first) &&
    RingsOnSameLayer(seed1.second,seed2.second);
  
  return result;
}

bool RoadMaker::SameRoadSet	(	Roads::RoadSet	set1,
		Roads::RoadSet	set2
	)

Definition at line 2399 of file RoadMaker.cc.

Referenced by AddInnerSeedRing().

                                                                 {
  //
  // check if input roadsets contains exactly the same rings
  //

  // calculate how many rings are in both input sets
  unsigned int nRingsSet1 = 0;
  unsigned int nRingsSet2 = 0;
  for ( Roads::RoadSet::iterator vector1 = set1.begin();
    vector1 != set1.end();
    ++vector1 ) {
    nRingsSet1 += vector1->size();
  }

  for ( Roads::RoadSet::iterator vector2 = set2.begin();
    vector2 != set2.end();
    ++vector2 ) {
    nRingsSet2 += vector2->size();
  }

  // if one of the input sets has more rings than the other, they cannot be the same
  if ( nRingsSet1 != nRingsSet2 ) {
    return false;
  }

  bool different = false;
  for ( Roads::RoadSet::iterator vector1 = set1.begin();
    vector1 != set1.end();
    ++vector1 ) {
    for ( std::vector<const Ring*>::iterator ring1 = vector1->begin();
      ring1 != vector1->end();
      ++ring1 ) {
      bool included = false;
      for ( Roads::RoadSet::iterator vector2 = set2.begin();
        vector2 != set2.end();
        ++vector2 ) {
    for ( std::vector<const Ring*>::iterator ring2 = vector2->begin();
          ring2 != vector2->end();
          ++ring2 ) {
      if ( (*ring1) == (*ring2) ) {
        included = true;
      }
    }
      }
      if ( !included ) {
    different = true;
      }
    }
  }

  return !different;
}

Roads::RoadSet RoadMaker::SortRingsIntoLayers

(

std::vector< const Ring * >

input

)

Definition at line 2452 of file RoadMaker.cc.

References RecoTauValidation_cfi::reference, relativeConstraints::ring, RingsOnSameLayer(), runtimedef::set(), python.multivaluedict::sort(), and tmp.

Referenced by RingsCompatibleWithSeed().

                                                                        {
  //
  // sort rings from input into layer structure of RoadSet
  //

  // return value
  Roads::RoadSet set;

  // sort rings in input by their center in rz, do it twice
  std::sort(input.begin(),input.end(),SortRingsByZR());

  const Ring *reference = (*(input.begin()));
  std::vector<const Ring*> tmp;
  tmp.push_back(reference);
  for (std::vector<const Ring*>::iterator ring = ++input.begin();
       ring != input.end();
       ++ring ) {
    if ( RingsOnSameLayer(reference,(*ring)) ) {
      reference = (*ring);
      tmp.push_back(reference);
    } else {
      set.push_back(tmp);
      tmp.clear();
      reference = (*ring);
      tmp.push_back(reference);
    }
  }

  if ( tmp.size() > 0 ) {
    set.push_back(tmp);
  }

  // order layers in set
  std::sort(set.begin(),set.end(),SortLayersByZR());

//   set.push_back(input);
  return set;
}

Member Data Documentation

std::vector<const Ring*> RoadMaker::innerSeedRings1_

private

Definition at line 115 of file RoadMaker.h.

Referenced by collectInnerSeedRings1(), collectInnerTECSeedRings1(), collectInnerTIBSeedRings1(), collectInnerTIDSeedRings1(), collectInnerTOBSeedRings1(), and constructRoads().

std::vector<const Ring*> RoadMaker::innerSeedRings2_

private

Definition at line 116 of file RoadMaker.h.

Referenced by collectInnerSeedRings2(), collectInnerTECSeedRings2(), collectInnerTIBSeedRings2(), collectInnerTIDSeedRings2(), collectInnerTOBSeedRings2(), and constructRoads().

std::vector<const Ring*> RoadMaker::innerSeedRings_

private

Definition at line 114 of file RoadMaker.h.

Referenced by AddInnerSeedRing(), collectInnerSeedRings(), collectInnerTECSeedRings(), collectInnerTIBSeedRings(), collectInnerTIDSeedRings(), collectInnerTOBSeedRings(), and constructRoads().

std::vector<const Ring*> RoadMaker::outerSeedRings1_

private

Definition at line 118 of file RoadMaker.h.

Referenced by collectOuterSeedRings1(), collectOuterTECSeedRings1(), collectOuterTIBSeedRings1(), collectOuterTOBSeedRings1(), and constructRoads().

std::vector<const Ring*> RoadMaker::outerSeedRings_

private

Definition at line 117 of file RoadMaker.h.

Referenced by AddOuterSeedRing(), collectOuterSeedRings(), collectOuterTECSeedRings(), collectOuterTIBSeedRings(), collectOuterTOBSeedRings(), and constructRoads().

const Rings* RoadMaker::rings_

private

Definition at line 108 of file RoadMaker.h.

Referenced by collectInnerTECSeedRings(), collectInnerTECSeedRings1(), collectInnerTECSeedRings2(), collectInnerTIBSeedRings(), collectInnerTIBSeedRings1(), collectInnerTIBSeedRings2(), collectInnerTIDSeedRings(), collectInnerTIDSeedRings1(), collectInnerTIDSeedRings2(), collectInnerTOBSeedRings(), collectInnerTOBSeedRings1(), collectInnerTOBSeedRings2(), collectOuterTECSeedRings(), collectOuterTECSeedRings1(), collectOuterTIBSeedRings(), collectOuterTIBSeedRings1(), collectOuterTOBSeedRings(), collectOuterTOBSeedRings1(), and RingsCompatibleWithSeed().

Roads* RoadMaker::roads_

private

Definition at line 110 of file RoadMaker.h.

Referenced by AddRoad(), constructRoads(), getRoads(), and RoadMaker().

SeedingType RoadMaker::seedingType_

private

Definition at line 112 of file RoadMaker.h.

Referenced by constructRoads().

GeometryStructure RoadMaker::structure_

private

Definition at line 111 of file RoadMaker.h.

Referenced by collectInnerSeedRings(), collectInnerSeedRings1(), collectInnerSeedRings2(), collectInnerTECSeedRings(), collectInnerTECSeedRings1(), collectInnerTECSeedRings2(), collectInnerTIBSeedRings(), collectInnerTIBSeedRings1(), collectInnerTIBSeedRings2(), collectInnerTIDSeedRings(), collectInnerTIDSeedRings1(), collectInnerTIDSeedRings2(), collectInnerTOBSeedRings(), collectInnerTOBSeedRings1(), collectInnerTOBSeedRings2(), collectOuterSeedRings(), collectOuterSeedRings1(), collectOuterTECSeedRings(), collectOuterTECSeedRings1(), collectOuterTIBSeedRings(), collectOuterTIBSeedRings1(), collectOuterTOBSeedRings(), collectOuterTOBSeedRings1(), and RoadMaker().

float RoadMaker::zBS_

private

Definition at line 120 of file RoadMaker.h.

Referenced by LinesThroughRingAndBS(), and RoadMaker().