TrajectorySeedProducer Class Reference

#include <TrajectorySeedProducer.h>

Inheritance diagram for TrajectorySeedProducer:

Classes
struct	LayerSpec

Public Types
enum	Side { BARREL =0, NEG_ENDCAP =-1, POS_ENDCAP =1 }
enum	SubDet {   NotValid, PXB, PXD, TIB,   TID, TOB, TEC }
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Member Functions
virtual void	beginRun (edm::Run const &run, const edm::EventSetup &es) override
virtual void	produce (edm::Event &e, const edm::EventSetup &es) override
Side	setLayerSpecSide (const std::string &layerSpecSide) const
	TrajectorySeedProducer (const edm::ParameterSet &conf)
virtual	~TrajectorySeedProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
bool	compatibleWithBeamAxis (GlobalPoint &gpos1, GlobalPoint &gpos2, double error, bool forward, unsigned algo) const
void	stateOnDet (const TrajectoryStateOnSurface &ts, unsigned int detid, PTrajectoryStateOnDet &pts) const
	A mere copy (without memory leak) of an existing tracking method. More...

Private Attributes
unsigned int	absMinRecHits
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken
std::vector< unsigned int >	firstHitSubDetectorNumber
std::vector< std::vector < unsigned int > >	firstHitSubDetectors
edm::InputTag	hitProducer
std::vector< double >	maxD0
std::vector< double >	maxZ0
std::vector< unsigned >	minRecHits
bool	newSyntax
std::vector< unsigned int >	numberOfHits
std::vector< double >	originHalfLength
std::vector< double >	originpTMin
std::vector< double >	originRadius
std::vector< edm::InputTag >	primaryVertices
std::vector< double >	pTMin
edm::EDGetTokenT < SiTrackerGSMatchedRecHit2DCollection >	recHitToken
std::vector< edm::EDGetTokenT < reco::VertexCollection > >	recoVertexToken
bool	rejectOverlaps
std::vector< unsigned int >	secondHitSubDetectorNumber
std::vector< std::vector < unsigned int > >	secondHitSubDetectors
bool	seedCleaning
std::vector< std::string >	seedingAlgo
bool	selectMuons
edm::EDGetTokenT < edm::SimTrackContainer >	simTrackToken
edm::EDGetTokenT < edm::SimVertexContainer >	simVertexToken
edm::InputTag	theBeamSpot
const MagneticFieldMap *	theFieldMap
const TrackerGeometry *	theGeometry
std::vector< std::vector < LayerSpec > >	theLayersInSets
const MagneticField *	theMagField
PropagatorWithMaterial *	thePropagator
std::vector< unsigned int >	thirdHitSubDetectorNumber
std::vector< std::vector < unsigned int > >	thirdHitSubDetectors
std::vector< const reco::VertexCollection * >	vertices
double	x0
double	y0
double	z0
std::vector< double >	zVertexConstraint

Additional Inherited Members
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 34 of file TrajectorySeedProducer.h.

Member Enumeration Documentation

enum TrajectorySeedProducer::Side

Enumerator
BARREL
NEG_ENDCAP
POS_ENDCAP

Definition at line 50 of file TrajectorySeedProducer.h.

{ BARREL=0, NEG_ENDCAP=-1, POS_ENDCAP=1};

enum TrajectorySeedProducer::SubDet

Enumerator
NotValid
PXB
PXD
TIB
TID
TOB
TEC

Definition at line 48 of file TrajectorySeedProducer.h.

{ NotValid, PXB, PXD, TIB, TID, TOB, TEC};

Constructor & Destructor Documentation

TrajectorySeedProducer::TrajectorySeedProducer ( const edm::ParameterSet &  conf )

explicit

Definition at line 43 of file TrajectorySeedProducer.cc.

References absMinRecHits, BARREL, beamSpotToken, edm::hlt::Exception, firstHitSubDetectorNumber, firstHitSubDetectors, edm::ParameterSet::getParameter(), hitProducer, i, TrajectorySeedProducer::LayerSpec::idLayer, CombinatorialSeedGeneratorForBeamHalo_cfi::layerList, geometryCSVtoXML::line, maxD0, maxZ0, minRecHits, TrajectorySeedProducer::LayerSpec::name, NEG_ENDCAP, newSyntax, numberOfHits, originHalfLength, originpTMin, originRadius, POS_ENDCAP, primaryVertices, pTMin, PXB, PXD, recHitToken, recoVertexToken, secondHitSubDetectorNumber, secondHitSubDetectors, seedCleaning, seedingAlgo, selectMuons, TrajectorySeedProducer::LayerSpec::side, simTrackToken, simVertexToken, AlCaHLTBitMon_QueryRunRegistry::string, TrajectorySeedProducer::LayerSpec::subDet, TEC, theBeamSpot, theLayersInSets, thirdHitSubDetectorNumber, thirdHitSubDetectors, TIB, TID, TOB, and zVertexConstraint.

                                                                          :thePropagator(0)
{  

  // The input tag for the beam spot
  theBeamSpot = conf.getParameter<edm::InputTag>("beamSpot");

  // The name of the TrajectorySeed Collections
  seedingAlgo = conf.getParameter<std::vector<std::string> >("seedingAlgo");
  for ( unsigned i=0; i<seedingAlgo.size(); ++i )
    produces<TrajectorySeedCollection>(seedingAlgo[i]);

  // The smallest true pT for a track candidate
  pTMin = conf.getParameter<std::vector<double> >("pTMin");
  if ( pTMin.size() != seedingAlgo.size() ) 
    throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
      << " WARNING : pTMin does not have the proper size "
      << std::endl;

  for ( unsigned i=0; i<pTMin.size(); ++i )
    pTMin[i] *= pTMin[i];  // Cut is done of perp2() - CPU saver
  
  // The smallest number of Rec Hits for a track candidate
  minRecHits = conf.getParameter<std::vector<unsigned int> >("minRecHits");
  if ( minRecHits.size() != seedingAlgo.size() ) 
    throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
      << " WARNING : minRecHits does not have the proper size "
      << std::endl;
  // Set the overall number hits to be checked
  absMinRecHits = 0;
  for ( unsigned ialgo=0; ialgo<minRecHits.size(); ++ialgo ) 
    if ( minRecHits[ialgo] > absMinRecHits ) absMinRecHits = minRecHits[ialgo];

  // The smallest true impact parameters (d0 and z0) for a track candidate
  maxD0 = conf.getParameter<std::vector<double> >("maxD0");
  if ( maxD0.size() != seedingAlgo.size() ) 
    throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
      << " WARNING : maxD0 does not have the proper size "
      << std::endl;

  maxZ0 = conf.getParameter<std::vector<double> >("maxZ0");
  if ( maxZ0.size() != seedingAlgo.size() ) 
    throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
      << " WARNING : maxZ0 does not have the proper size "
      << std::endl;

  // The name of the hit producer
  hitProducer = conf.getParameter<edm::InputTag>("HitProducer");

  // The cuts for seed cleaning
  seedCleaning = conf.getParameter<bool>("seedCleaning");

  // Number of hits needed for a seed
  numberOfHits = conf.getParameter<std::vector<unsigned int> >("numberOfHits");
  if ( numberOfHits.size() != seedingAlgo.size() ) 
    throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
      << " WARNING : numberOfHits does not have the proper size "
      << std::endl;

  // Seeding based on muons
  selectMuons = conf.getParameter<bool>("selectMuons");

  // Layers
  newSyntax = conf.getParameter<bool>("newSyntax");
  if (newSyntax) {
    std::vector<std::string> layerList = conf.getParameter<std::vector<std::string> >("layerList");
    //for (unsigned i=0; i<layerList.size();i++) std::cout << "------- Layers = " << layerList[i] << std::endl; 

    for(std::vector<std::string>::const_iterator it=layerList.begin(); it < layerList.end(); ++it) {
      std::vector<LayerSpec> tempResult;
      std::string line = *it;
      std::string::size_type pos=0;
      while (pos != std::string::npos) {
        pos=line.find("+");
        std::string layer = line.substr(0, pos);
        //
        LayerSpec layerSpec;
        layerSpec.name = line.substr(0, pos);;
        // Possible names: BPix(1-3) || FPix(1-2)(pos, neg) || TIB(1-4) || TID(1-3)(pos, neg) || TOB(1-6) || TEC(1-9)(pos, neg)
        // BPix(1-3)
        if (layerSpec.name.substr(0,4)=="BPix" ) {
          layerSpec.subDet=PXB;
          layerSpec.side=BARREL;
          layerSpec.idLayer = std::atoi(layerSpec.name.substr(4,1).c_str());
          if (layerSpec.idLayer>3 || layerSpec.idLayer==0) {
            throw cms::Exception("FastSimulation/Tracking/python")
              << "Bad data naming in IterativeInitialStep_cff.py  iterativeInitialSeeds.layerList" << std::endl
              << "Layers: " << layerSpec.name << ", number needs to be in range 1-3" << std::endl;
          }
        }
        // FPix(1-2)(pos, neg)
        else if (layerSpec.name.substr(0,4)=="FPix" ) {
          layerSpec.subDet=PXD;
          if(layerSpec.name.substr(layerSpec.name.size()-3)=="pos")
            layerSpec.side = POS_ENDCAP;
          else //no validation if it's not neg
            layerSpec.side = NEG_ENDCAP;
          layerSpec.idLayer = std::atoi(layerSpec.name.substr(4,1).c_str());
          if (layerSpec.idLayer>2 || layerSpec.idLayer==0) {
            throw cms::Exception("FastSimulation/Tracking/python")
              << "Bad data naming in IterativeInitialStep_cff.py  iterativeInitialSeeds.layerList" << std::endl
              << "Layers: " << layerSpec.name << ", number needs to be in range 1-2" << std::endl;
          }
        }
        // TIB(1-4)
        else if (layerSpec.name.substr(0,3)=="TIB" || layerSpec.name.substr(0,4)=="MTIB") {
          layerSpec.subDet=TIB;
          layerSpec.side=BARREL;
          layerSpec.idLayer = std::atoi(layerSpec.name.substr(layerSpec.name.find('B')+1,1).c_str());
          if (layerSpec.idLayer>4 || layerSpec.idLayer==0) {
            throw cms::Exception("FastSimulation/Tracking/python")
              << "Bad data naming in IterativeInitialStep_cff.py  iterativeInitialSeeds.layerList" << std::endl
              << "Layers: " << layerSpec.name << ", number needs to be in range 1-4" << std::endl;
          }
        }
        // TID(1-3)(pos, neg)
        else if (layerSpec.name.substr(0,3)=="TID" || layerSpec.name.substr(0,4)=="MTID") {
          layerSpec.subDet=TID;
          if(layerSpec.name.substr(layerSpec.name.size()-3)=="pos")
            layerSpec.side = POS_ENDCAP;
          else
            layerSpec.side = NEG_ENDCAP;
          layerSpec.idLayer = std::atoi(layerSpec.name.substr(layerSpec.name.find('D')+1,1).c_str());
          if (layerSpec.idLayer>3 || layerSpec.idLayer==0) {
            throw cms::Exception("FastSimulation/Tracking/python")
              << "Bad data naming in IterativeInitialStep_cff.py  iterativeInitialSeeds.layerList" << std::endl
              << "Layers: " << layerSpec.name << ", number needs to be in range 1-3" << std::endl;
          }
        }
        // TOB(1-6)
        else if (layerSpec.name.substr(0,3)=="TOB" ) {
          layerSpec.subDet=TOB;
          layerSpec.side=BARREL;
          layerSpec.idLayer = std::atoi(layerSpec.name.substr(3,1).c_str());
          if (layerSpec.idLayer>6 || layerSpec.idLayer==0) {
            throw cms::Exception("FastSimulation/Tracking/python")
              << "Bad data naming in IterativeInitialStep_cff.py  iterativeInitialSeeds.layerList" << std::endl
              << "Layers: " << layerSpec.name << ", number needs to be in range 1-6" << std::endl;
          }
        }
        // TEC(1-9)(pos, neg)
        else if (layerSpec.name.substr(0,3)=="TEC" || layerSpec.name.substr(0,4)=="MTEC") {
          layerSpec.subDet=TEC;
          if(layerSpec.name.substr(layerSpec.name.size()-3)=="pos")
            layerSpec.side = POS_ENDCAP;
          else
            layerSpec.side = NEG_ENDCAP;
          layerSpec.idLayer = std::atoi(layerSpec.name.substr(layerSpec.name.find('C')+1,1).c_str());
          if (layerSpec.idLayer>9 || layerSpec.idLayer==0) {
            throw cms::Exception("FastSimulation/Tracking/python")
              << "Bad data naming in IterativeInitialStep_cff.py  iterativeInitialSeeds.layerList" << std::endl
              << "Layers: " << layerSpec.name << ", number needs to be in range 1-9" << std::endl;
          }
        } 
        else { 
          throw cms::Exception("FastSimulation/Tracking/python")
            << "Bad data naming in IterativeInitialStep_cff.py  iterativeInitialSeeds.layerList" << std::endl
            << "Layer: " << layerSpec.name << ", shouldn't exist" << std::endl
            << "Case sensitive names: BPix FPix TIB MTIB TID MTID TOB TEC MTEC" << std::endl;
        }
        //
        tempResult.push_back(layerSpec);
        line=line.substr(pos+1,std::string::npos); 
      }
      theLayersInSets.push_back(tempResult);
    }

    //prints theLayersInSets
    /* for (std::vector<std::vector<LayerSpec> >::const_iterator it = theLayersInSets.begin(); it != theLayersInSets.end(); ++it) {
      for (std::vector<LayerSpec>::const_iterator is = it->begin(); is != it->end(); ++is) {
        std::cout << is->name << " | " << is->subDet << " | " << is->idLayer << " | " << is->side << std::endl;
      }
      std::cout << "---" << std::endl;
    } */
  } else {
    // TO BE DELETED (AG)
    firstHitSubDetectorNumber = 
      conf.getParameter<std::vector<unsigned int> >("firstHitSubDetectorNumber");
    if ( firstHitSubDetectorNumber.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : firstHitSubDetectorNumber does not have the proper size "
    << std::endl;
    
    std::vector<unsigned int> firstSubDets = 
      conf.getParameter<std::vector<unsigned int> >("firstHitSubDetectors");
    unsigned isub1 = 0;
    unsigned check1 = 0;
    firstHitSubDetectors.resize(seedingAlgo.size());
    for ( unsigned ialgo=0; ialgo<firstHitSubDetectorNumber.size(); ++ialgo ) { 
      check1 += firstHitSubDetectorNumber[ialgo];
      for ( unsigned idet=0; idet<firstHitSubDetectorNumber[ialgo]; ++idet ) { 
    firstHitSubDetectors[ialgo].push_back(firstSubDets[isub1++]);
      }
    }
    if ( firstSubDets.size() != check1 ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : firstHitSubDetectors does not have the proper size (should be " << check1 << ")"
    << std::endl;
    
    
    secondHitSubDetectorNumber = 
      conf.getParameter<std::vector<unsigned int> >("secondHitSubDetectorNumber");
    if ( secondHitSubDetectorNumber.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : secondHitSubDetectorNumber does not have the proper size "
    << std::endl;
    
    std::vector<unsigned int> secondSubDets = 
      conf.getParameter<std::vector<unsigned int> >("secondHitSubDetectors");
    unsigned isub2 = 0;
    unsigned check2 = 0;
    secondHitSubDetectors.resize(seedingAlgo.size());
    for ( unsigned ialgo=0; ialgo<secondHitSubDetectorNumber.size(); ++ialgo ) { 
      check2 += secondHitSubDetectorNumber[ialgo];
      for ( unsigned idet=0; idet<secondHitSubDetectorNumber[ialgo]; ++idet ) { 
    secondHitSubDetectors[ialgo].push_back(secondSubDets[isub2++]);
      }
    }
    if ( secondSubDets.size() != check2 ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : secondHitSubDetectors does not have the proper size (should be " << check2 << ")"
    << std::endl;
    
    thirdHitSubDetectorNumber = 
      conf.getParameter<std::vector<unsigned int> >("thirdHitSubDetectorNumber");
    if ( thirdHitSubDetectorNumber.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : thirdHitSubDetectorNumber does not have the proper size "
    << std::endl;
    
    std::vector<unsigned int> thirdSubDets = 
      conf.getParameter<std::vector<unsigned int> >("thirdHitSubDetectors");
    unsigned isub3 = 0;
    unsigned check3 = 0;
    thirdHitSubDetectors.resize(seedingAlgo.size());
    for ( unsigned ialgo=0; ialgo<thirdHitSubDetectorNumber.size(); ++ialgo ) { 
      check3 += thirdHitSubDetectorNumber[ialgo];
      for ( unsigned idet=0; idet<thirdHitSubDetectorNumber[ialgo]; ++idet ) { 
    thirdHitSubDetectors[ialgo].push_back(thirdSubDets[isub3++]);
      }
    }
    if ( thirdSubDets.size() != check3 ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : thirdHitSubDetectors does not have the proper size (should be " << check3 << ")"
    << std::endl;
    }
    
    originRadius = conf.getParameter<std::vector<double> >("originRadius");
    if ( originRadius.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : originRadius does not have the proper size "
    << std::endl;
    
    originHalfLength = conf.getParameter<std::vector<double> >("originHalfLength");
    if ( originHalfLength.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : originHalfLength does not have the proper size "
    << std::endl;
    
    originpTMin = conf.getParameter<std::vector<double> >("originpTMin");
    if ( originpTMin.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : originpTMin does not have the proper size "
    << std::endl;
    
    primaryVertices = conf.getParameter<std::vector<edm::InputTag> >("primaryVertices");
    if ( primaryVertices.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : primaryVertices does not have the proper size "
    << std::endl;
    
    zVertexConstraint = conf.getParameter<std::vector<double> >("zVertexConstraint");
    if ( zVertexConstraint.size() != seedingAlgo.size() ) 
      throw cms::Exception("FastSimulation/TrajectorySeedProducer ") 
    << " WARNING : zVertexConstraint does not have the proper size "
    << std::endl;
  //removed - }

    // consumes
    beamSpotToken = consumes<reco::BeamSpot>(theBeamSpot);
    edm::InputTag _label("famosSimHits");
    simTrackToken = consumes<edm::SimTrackContainer>(_label);
    simVertexToken = consumes<edm::SimVertexContainer>(_label);
    recHitToken = consumes<SiTrackerGSMatchedRecHit2DCollection>(hitProducer);
    for ( unsigned ialgo=0; ialgo<seedingAlgo.size(); ++ialgo ) {
      _label = edm::InputTag(primaryVertices[ialgo]);
      recoVertexToken.push_back(consumes<reco::VertexCollection>(_label));
    }
}

TrajectorySeedProducer::~TrajectorySeedProducer ( )

virtual

Definition at line 333 of file TrajectorySeedProducer.cc.

References gather_cfg::cout, and thePropagator.

                                                {
  
  if(thePropagator) delete thePropagator;

  // do nothing
#ifdef FAMOS_DEBUG
  std::cout << "TrajectorySeedProducer destructed" << std::endl;
#endif

} 

Member Function Documentation

void TrajectorySeedProducer::beginRun ( edm::Run const &  run, const edm::EventSetup &  es  )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 345 of file TrajectorySeedProducer.cc.

References alongMomentum, g, geometry, edm::EventSetup::get(), theFieldMap, theGeometry, theMagField, and thePropagator.

                                                                        {

  //services
  //  es.get<TrackerRecoGeometryRecord>().get(theGeomSearchTracker);

  edm::ESHandle<MagneticField>          magField;
  edm::ESHandle<TrackerGeometry>        geometry;
  edm::ESHandle<MagneticFieldMap>       magFieldMap;


  es.get<IdealMagneticFieldRecord>().get(magField);
  es.get<TrackerDigiGeometryRecord>().get(geometry);
  es.get<MagneticFieldMapRecord>().get(magFieldMap);

  theMagField = &(*magField);
  theGeometry = &(*geometry);
  theFieldMap = &(*magFieldMap);

  thePropagator = new PropagatorWithMaterial(alongMomentum,0.105,&(*theMagField)); 

  const GlobalPoint g(0.,0.,0.);

}

bool TrajectorySeedProducer::compatibleWithBeamAxis ( GlobalPoint &  gpos1, GlobalPoint &  gpos2, double  error, bool  forward, unsigned  algo  ) const

private

Check that the seed is compatible with a track coming from within a cylinder of radius originRadius, with a decent pT.

Check that the seed is compatible with a track coming from within a cylinder of radius originRadius, with a decent pT, and propagate to the distance of closest approach, for the appropriate charge

Definition at line 828 of file TrajectorySeedProducer.cc.

References gather_cfg::cout, relativeConstraints::error, max(), bookConverter::min, originHalfLength, originpTMin, originRadius, BaseParticlePropagator::propagateToBeamCylinder(), RawParticle::R(), seedCleaning, mathSSE::sqrt(), theFieldMap, vertices, PV3DBase< T, PVType, FrameType >::x(), x0, PV3DBase< T, PVType, FrameType >::y(), y0, PV3DBase< T, PVType, FrameType >::z(), RawParticle::Z(), z0, and zVertexConstraint.

Referenced by produce().

                                                {

  if ( !seedCleaning ) return true;

  // The hits 1 and 2 positions, in HepLorentzVector's
  XYZTLorentzVector thePos1(gpos1.x(),gpos1.y(),gpos1.z(),0.);
  XYZTLorentzVector thePos2(gpos2.x(),gpos2.y(),gpos2.z(),0.);
#ifdef FAMOS_DEBUG
  std::cout << "ThePos1 = " << thePos1 << std::endl;
  std::cout << "ThePos2 = " << thePos2 << std::endl;
#endif


  // Create new particles that pass through the second hit with pT = ptMin 
  // and charge = +/-1
  
  // The momentum direction is by default joining the two hits 
  XYZTLorentzVector theMom2 = (thePos2-thePos1);

  // The corresponding RawParticle, with an (irrelevant) electric charge
  // (The charge is determined in the next step)
  ParticlePropagator myPart(theMom2,thePos2,1.,theFieldMap);

  bool intersect = myPart.propagateToBeamCylinder(thePos1,originRadius[algo]*1.);
  if ( !intersect ) return false;

#ifdef FAMOS_DEBUG
  std::cout << "MyPart R = " << myPart.R() << "\t Z = " << myPart.Z() 
        << "\t pT = " << myPart.Pt() << std::endl;
#endif

  // Check if the constraints are satisfied
  // 1. pT at cylinder with radius originRadius
  if ( myPart.Pt() < originpTMin[algo] ) return false;

  // 2. Z compatible with beam spot size
  if ( fabs(myPart.Z()-z0) > originHalfLength[algo] ) return false;

  // 3. Z compatible with one of the primary vertices (always the case if no primary vertex)
  const reco::VertexCollection* theVertices = vertices[algo];
  if (!theVertices) return true;
  unsigned nVertices = theVertices->size();
  if ( !nVertices || zVertexConstraint[algo] < 0. ) return true;
  // Radii of the two hits with respect to the beam spot position
  double R1 = std::sqrt ( (thePos1.X()-x0)*(thePos1.X()-x0) 
            + (thePos1.Y()-y0)*(thePos1.Y()-y0) );
  double R2 = std::sqrt ( (thePos2.X()-x0)*(thePos2.X()-x0) 
            + (thePos2.Y()-y0)*(thePos2.Y()-y0) );
  // Loop on primary vertices
  for ( unsigned iv=0; iv<nVertices; ++iv ) { 
    // Z position of the primary vertex
    double zV = (*theVertices)[iv].z();
    // Constraints on the inner hit
    double checkRZ1 = forward ?
      (thePos1.Z()-zV+zVertexConstraint[algo]) / (thePos2.Z()-zV+zVertexConstraint[algo]) * R2 : 
      -zVertexConstraint[algo] + R1/R2*(thePos2.Z()-zV+zVertexConstraint[algo]);
    double checkRZ2 = forward ?
      (thePos1.Z()-zV-zVertexConstraint[algo])/(thePos2.Z()-zV-zVertexConstraint[algo]) * R2 :
      +zVertexConstraint[algo] + R1/R2*(thePos2.Z()-zV-zVertexConstraint[algo]);
    double checkRZmin = std::min(checkRZ1,checkRZ2)-3.*error;
    double checkRZmax = std::max(checkRZ1,checkRZ2)+3.*error;
    // Check if the innerhit is within bounds
    bool compat = forward ?
      checkRZmin < R1 && R1 < checkRZmax : 
      checkRZmin < thePos1.Z()-zV && thePos1.Z()-zV < checkRZmax; 
    // If it is, just return ok
    if ( compat ) return compat;
  }
  // Otherwise, return not ok
  return false;

}  

void TrajectorySeedProducer::produce ( edm::Event &  e, const edm::EventSetup &  es  )

overridevirtual

Implements edm::EDProducer.

Definition at line 371 of file TrajectorySeedProducer.cc.

References funct::abs(), absMinRecHits, alongMomentum, beamSpotToken, DeDxDiscriminatorTools::charge(), TrackingRecHit::clone(), compatibleWithBeamAxis(), gather_cfg::cout, relativeConstraints::error, PV3DBase< T, PVType, FrameType >::eta(), firstHitSubDetectors, edm::OwnVector< T, P >::front(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByToken(), TrajectoryStateOnSurface::globalMomentum(), TrackerRecHit::globalPosition(), TrackingRecHit::hit(), TrackerGeometry::idToDet(), TrackerRecHit::isForward(), TrackerRecHit::isOnRequestedDet(), TrackerRecHit::isOnTheSameLayer(), TrajectoryStateOnSurface::isValid(), TrackerRecHit::largerError(), TrackerRecHit::layerNumber(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localMomentum(), LocalTrajectoryError::matrix(), maxD0, maxZ0, minRecHits, CoreSimTrack::momentum(), FreeTrajectoryState::momentum(), newSyntax, numberOfHits, convertSQLitetoXML_cfg::output, AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), CoreSimVertex::position(), position, LocalTrajectoryError::positionError(), edm::ESHandle< class >::product(), pTMin, edm::OwnVector< T, P >::push_back(), edm::Event::put(), DetId::rawId(), recHitToken, recoVertexToken, secondHitSubDetectors, seedingAlgo, selectMuons, RawParticle::setCharge(), simTrackToken, simVertexToken, edm::OwnVector< T, P >::size(), mathSSE::sqrt(), stateOnDet(), TrackerRecHit::subDetId(), GeomDet::surface(), TrajectoryStateOnSurface::surfaceSide(), theGeometry, theLayersInSets, theMagField, thePropagator, thirdHitSubDetectors, CoreSimTrack::type(), vertices, SimTrack::vertIndex(), x, x0, BaseParticlePropagator::xyImpactParameter(), detailsBasic3DVector::y, y0, detailsBasic3DVector::z, z0, and BaseParticlePropagator::zImpactParameter().

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                    {        


  //  if( seedingAlgo[0] ==  "FourthPixelLessPairs") std::cout << "Seed producer in 4th iteration " << std::endl;

#ifdef FAMOS_DEBUG
  std::cout << "################################################################" << std::endl;
  std::cout << " TrajectorySeedProducer produce init " << std::endl;
#endif

  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHand;
  es.get<IdealGeometryRecord>().get(tTopoHand);
  const TrackerTopology *tTopo=tTopoHand.product();


  unsigned nSimTracks = 0;
  unsigned nTracksWithHits = 0;
  unsigned nTracksWithPT = 0;
  unsigned nTracksWithD0Z0 = 0;
  //  unsigned nTrackCandidates = 0;
  PTrajectoryStateOnDet initialState;
  
  // Output
  std::vector<TrajectorySeedCollection*>
    output(seedingAlgo.size(),static_cast<TrajectorySeedCollection*>(0));
  for ( unsigned ialgo=0; ialgo<seedingAlgo.size(); ++ialgo ) { 
    //    std::auto_ptr<TrajectorySeedCollection> p(new TrajectorySeedCollection );
    output[ialgo] = new TrajectorySeedCollection;
  }
  
  // Beam spot
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  e.getByToken(beamSpotToken,recoBeamSpotHandle); 
  math::XYZPoint BSPosition_ = recoBeamSpotHandle->position();

  //not used anymore. take the value from the py

  //double sigmaZ=recoBeamSpotHandle->sigmaZ();
  //double sigmaZ0Error=recoBeamSpotHandle->sigmaZ0Error();
  //double sigmaz0=sqrt(sigmaZ*sigmaZ+sigmaZ0Error*sigmaZ0Error);
  x0 = BSPosition_.X();
  y0 = BSPosition_.Y();
  z0 = BSPosition_.Z();

  // SimTracks and SimVertices
  edm::Handle<edm::SimTrackContainer> theSimTracks;
  e.getByToken(simTrackToken,theSimTracks);
  
  edm::Handle<edm::SimVertexContainer> theSimVtx;
  e.getByToken(simVertexToken,theSimVtx);

#ifdef FAMOS_DEBUG
  std::cout << " Step A: SimTracks found " << theSimTracks->size() << std::endl;
#endif
  
  //  edm::Handle<SiTrackerGSRecHit2DCollection> theGSRecHits;
  edm::Handle<SiTrackerGSMatchedRecHit2DCollection> theGSRecHits;
  e.getByToken(recHitToken, theGSRecHits);
  
  // No tracking attempted if no hits (but put an empty collection in the event)!
#ifdef FAMOS_DEBUG
  std::cout << " Step B: Full GS RecHits found " << theGSRecHits->size() << std::endl;
#endif
  if(theGSRecHits->size() == 0) {
    for ( unsigned ialgo=0; ialgo<seedingAlgo.size(); ++ialgo ) {
      std::auto_ptr<TrajectorySeedCollection> p(output[ialgo]);
      e.put(p,seedingAlgo[ialgo]);
    }
    return;
  }
  
  // Primary vertices
  vertices = std::vector<const reco::VertexCollection*>
    (seedingAlgo.size(),static_cast<const reco::VertexCollection*>(0)); 
  for ( unsigned ialgo=0; ialgo<seedingAlgo.size(); ++ialgo ) { 
    //PAT Attempt!!!! 
   //originHalfLength[ialgo] = 3.*sigmaz0; // Overrides the configuration
    edm::Handle<reco::VertexCollection> aHandle;
    bool isVertexCollection = e.getByToken(recoVertexToken[ialgo],aHandle);
    if (!isVertexCollection ) continue;
    vertices[ialgo] = &(*aHandle);
  }
  
#ifdef FAMOS_DEBUG
  std::cout << " Step C: Loop over the RecHits, track by track " << std::endl;
#endif

  // The vector of simTrack Id's carrying GSRecHits
  const std::vector<unsigned> theSimTrackIds = theGSRecHits->ids();

  // loop over SimTrack Id's
  for ( unsigned tkId=0;  tkId != theSimTrackIds.size(); ++tkId ) {

#ifdef FAMOS_DEBUG
    std::cout << "Track number " << tkId << std::endl;
#endif

    ++nSimTracks;
    unsigned simTrackId = theSimTrackIds[tkId];
    const SimTrack& theSimTrack = (*theSimTracks)[simTrackId]; 
#ifdef FAMOS_DEBUG
    std::cout << "Pt = " << std::sqrt(theSimTrack.momentum().Perp2()) 
          << " eta " << theSimTrack.momentum().Eta()
          << " pdg ID " << theSimTrack.type()
          << std::endl;
#endif

    // Select only muons, if requested
    if (selectMuons && abs(theSimTrack.type()) != 13) continue;
    
    // Check that the sim track comes from the main vertex (loose cut)
    int vertexIndex = theSimTrack.vertIndex();
    const SimVertex& theSimVertex = (*theSimVtx)[vertexIndex]; 
#ifdef FAMOS_DEBUG
    std::cout << " o SimTrack " << theSimTrack << std::endl;
    std::cout << " o SimVertex " << theSimVertex << std::endl;
#endif
    
    BaseParticlePropagator theParticle = 
      BaseParticlePropagator( 
     RawParticle(XYZTLorentzVector(theSimTrack.momentum().px(),
                       theSimTrack.momentum().py(),
                       theSimTrack.momentum().pz(),
                       theSimTrack.momentum().e()),
             XYZTLorentzVector(theSimVertex.position().x(),
                       theSimVertex.position().y(),
                       theSimVertex.position().z(),
                       theSimVertex.position().t())),
                 0.,0.,4.);
    theParticle.setCharge((*theSimTracks)[simTrackId].charge());

    SiTrackerGSMatchedRecHit2DCollection::range theRecHitRange = theGSRecHits->get(simTrackId);
    SiTrackerGSMatchedRecHit2DCollection::const_iterator theRecHitRangeIteratorBegin = theRecHitRange.first;
    SiTrackerGSMatchedRecHit2DCollection::const_iterator theRecHitRangeIteratorEnd   = theRecHitRange.second;
    SiTrackerGSMatchedRecHit2DCollection::const_iterator iterRecHit;
    SiTrackerGSMatchedRecHit2DCollection::const_iterator iterRecHit1;
    SiTrackerGSMatchedRecHit2DCollection::const_iterator iterRecHit2;
    SiTrackerGSMatchedRecHit2DCollection::const_iterator iterRecHit3;

    // Check the number of layers crossed
    unsigned numberOfRecHits = 0;
    TrackerRecHit previousHit, currentHit;
    for ( iterRecHit = theRecHitRangeIteratorBegin; 
      iterRecHit != theRecHitRangeIteratorEnd; 
      ++iterRecHit) { 
      previousHit = currentHit;
      currentHit = TrackerRecHit(&(*iterRecHit),theGeometry,tTopo);
      if ( currentHit.isOnTheSameLayer(previousHit) ) continue;
      ++numberOfRecHits;
      if ( numberOfRecHits == absMinRecHits ) break;
    }

    // Loop on the successive seedings
    for ( unsigned int ialgo = 0; ialgo < seedingAlgo.size(); ++ialgo ) {
#ifdef FAMOS_DEBUG
      std::cout << "Algo " << seedingAlgo[ialgo] << std::endl;
#endif

      // Request a minimum number of RecHits for the track to give a seed.
#ifdef FAMOS_DEBUG
      std::cout << "The number of RecHits = " << numberOfRecHits << std::endl;
#endif
      if ( numberOfRecHits < minRecHits[ialgo] ) continue;
      ++nTracksWithHits;

      // Request a minimum pT for the sim track
      if ( theSimTrack.momentum().Perp2() < pTMin[ialgo] ) continue;
      ++nTracksWithPT;
      
      // Cut on sim track impact parameters
      if ( theParticle.xyImpactParameter(x0,y0) > maxD0[ialgo] ) continue;
      if ( fabs( theParticle.zImpactParameter(x0,y0) - z0 ) > maxZ0[ialgo] ) continue;
      ++nTracksWithD0Z0;
      
      std::vector<TrackerRecHit > 
    theSeedHits(numberOfHits[ialgo],
            static_cast<TrackerRecHit >(TrackerRecHit()));
      TrackerRecHit& theSeedHits0 = theSeedHits[0];
      TrackerRecHit& theSeedHits1 = theSeedHits[1];
      TrackerRecHit& theSeedHits2 = theSeedHits[2];
      
      bool compatible = false;
      
      for ( iterRecHit1 = theRecHitRangeIteratorBegin; iterRecHit1 != theRecHitRangeIteratorEnd; ++iterRecHit1) {
      //std::cout << (*iterRecHit1).localPosition().phi() << " | J - iterRecHit1"<< std::endl; // 
    theSeedHits[0] = TrackerRecHit(&(*iterRecHit1),theGeometry,tTopo);
#ifdef FAMOS_DEBUG
    std::cout << "The first hit position = " << theSeedHits0.globalPosition() << std::endl;
    std::cout << "The first hit subDetId = " << theSeedHits0.subDetId() << std::endl;
    std::cout << "The first hit layer    = " << theSeedHits0.layerNumber() << std::endl;
#endif  
    // Check if inside the requested detectors
    bool isInside = true;
    if (!selectMuons) {
    //(newSyntax) ? std::cout << "TRUE " : std::cout << "FALSE "; //J-
      if (newSyntax)
        isInside = false; // AG placeholder true 
      else
        isInside = theSeedHits0.subDetId() < firstHitSubDetectors[ialgo][0];
      //    bool isInside = theSeedHits0.subDetId() < firstHitSubDetectors[ialgo][0];
      if ( isInside ) continue;
    }

    // Check if on requested detectors
    // bool isOndet =  theSeedHits0.isOnRequestedDet(firstHitSubDetectors[ialgo]);
    bool isOndet = true;
    if (!selectMuons) {
      if (newSyntax)
      isOndet = theSeedHits[0].isOnRequestedDet(theLayersInSets);
      else
        isOndet = theSeedHits0.isOnRequestedDet(firstHitSubDetectors[ialgo], seedingAlgo[ialgo]);
      //std::cout << firstHitSubDetectors[ialgo][0] << " | " << seedingAlgo[ialgo] << " " << std::endl;  //seedingAlgo[iAlgo]: PixelTriplet, LowPtPixelTriplets, PixelPair, DetachedPixelTriplets, MixedTriplets, PixelLessPairs, TobTecLayerPairs......
      //    bool isOndet =  theSeedHits0.isOnRequestedDet(firstHitSubDetectors[ialgo], seedingAlgo[ialgo]);
      //    if ( !isOndet ) break;
      if ( !isOndet ) continue;
    }

#ifdef FAMOS_DEBUG
    std::cout << "Apparently the first hit is on the requested detector! " << std::endl;
#endif
    for ( iterRecHit2 = iterRecHit1+1; iterRecHit2 != theRecHitRangeIteratorEnd; ++iterRecHit2) {
      theSeedHits[1] = TrackerRecHit(&(*iterRecHit2),theGeometry,tTopo);
#ifdef FAMOS_DEBUG
      std::cout << "The second hit position = " << theSeedHits1.globalPosition() << std::endl;
      std::cout << "The second hit subDetId = " << theSeedHits1.subDetId() << std::endl;
      std::cout << "The second hit layer    = " << theSeedHits1.layerNumber() << std::endl;
#endif

      if (!selectMuons) {
        // Check if inside the requested detectors
        if (newSyntax) 
          isInside = false; // AG placeholder true
        else
          isInside = theSeedHits1.subDetId() < secondHitSubDetectors[ialgo][0];
        if ( isInside ) continue;

        // Check if on requested detectors
        if (newSyntax)
        isOndet = theSeedHits[0].isOnRequestedDet(theLayersInSets, theSeedHits[1]);
        else
          isOndet =  theSeedHits1.isOnRequestedDet(secondHitSubDetectors[ialgo], seedingAlgo[ialgo]);
        if ( !isOndet ) break;
      }
      // Check if on the same layer as previous hit
      if ( theSeedHits1.isOnTheSameLayer(theSeedHits0) ) continue;

#ifdef FAMOS_DEBUG
      std::cout << "Apparently the second hit is on the requested detector! " << std::endl;
#endif
      GlobalPoint gpos1 = theSeedHits0.globalPosition();
      GlobalPoint gpos2 = theSeedHits1.globalPosition();
      bool forward = theSeedHits0.isForward();
      double error = std::sqrt(theSeedHits0.largerError()+theSeedHits1.largerError());
      //      compatible = compatibleWithVertex(gpos1,gpos2,ialgo);
      //added out of desperation    
      if(seedingAlgo[0] == "PixelLess" ||  seedingAlgo[0] ==  "TobTecLayerPairs"){
        compatible = true;
        //std::cout << "Algo " << seedingAlgo[0] << "Det/layer = " << theSeedHits0.subDetId() << "/" <<  theSeedHits0.layerNumber() << std::endl;
      } else {
        compatible = compatibleWithBeamAxis(gpos1,gpos2,error,forward,ialgo);
      }
      
#ifdef FAMOS_DEBUG
      std::cout << "Algo" << seedingAlgo[0] << "\t Are the two hits compatible with the PV? " << compatible << std::endl;
#endif
      // Reject non suited pairs
      if ( !compatible ) continue;

#ifdef FAMOS_DEBUG
      std::cout << "Pair kept! " << std::endl;
#endif

      // Leave here if only two hits are required.
      if ( numberOfHits[ialgo] == 2 ) break; 
      
      compatible = false;
      // Check if there is a third satisfying hit otherwise
      for ( iterRecHit3 = iterRecHit2+1; iterRecHit3 != theRecHitRangeIteratorEnd; ++iterRecHit3) {
    
        theSeedHits[2] = TrackerRecHit(&(*iterRecHit3),theGeometry,tTopo);
#ifdef FAMOS_DEBUG
        std::cout << "The third hit position = " << theSeedHits2.globalPosition() << std::endl;
        std::cout << "The third hit subDetId = " << theSeedHits2.subDetId() << std::endl;
        std::cout << "The third hit layer    = " << theSeedHits2.layerNumber() << std::endl;
#endif

        if (!selectMuons) {
          // Check if inside the requested detectors
          if (newSyntax) 
          isInside = false; // AG placeholder
          else 
          isInside = theSeedHits2.subDetId() < thirdHitSubDetectors[ialgo][0];
          if ( isInside ) continue;
        
          // Check if on requested detectors
          if (newSyntax) 
          isOndet = theSeedHits[0].isOnRequestedDet(theLayersInSets, theSeedHits[1], theSeedHits[2]);
          else 
          isOndet =  theSeedHits2.isOnRequestedDet(thirdHitSubDetectors[ialgo], seedingAlgo[ialgo]);
          //        if ( !isOndet ) break;
          if ( !isOndet ) continue;
        }

        // Check if on the same layer as previous hit
        compatible = !(theSeedHits2.isOnTheSameLayer(theSeedHits1));

#ifdef FAMOS_DEBUG
        if ( compatible ) 
          std::cout << "Apparently the third hit is on the requested detector! " << std::endl;
#endif

        if ( compatible ) break;      

      }

      if ( compatible ) break;

    }

    if ( compatible ) break;
  
      }

      // There is no compatible seed for this track with this seeding algorithm 
      // Go to next algo
      if ( !compatible ) continue;
#ifdef FAMOS_DEBUG
      if ( compatible ) 
    std::cout << "@@@ There is at least a compatible seed" << std::endl;
      else
    std::cout << "@@@ There is no compatible seed" << std::endl;
#endif
      

#ifdef FAMOS_DEBUG
      std::cout << "Preparing to create the TrajectorySeed" << std::endl;
#endif
      // The seed is validated -> include in the collection
      // 1) Create the vector of RecHits
      edm::OwnVector<TrackingRecHit> recHits;
      for ( unsigned ih=0; ih<theSeedHits.size(); ++ih ) {
    TrackingRecHit* aTrackingRecHit = theSeedHits[ih].hit()->clone();
    recHits.push_back(aTrackingRecHit);
      }
#ifdef FAMOS_DEBUG
      std::cout << "with " << recHits.size() << " hits." << std::endl;
#endif

      // 2) Create the initial state
      //   a) origin vertex
      GlobalPoint  position((*theSimVtx)[vertexIndex].position().x(),
                (*theSimVtx)[vertexIndex].position().y(),
                (*theSimVtx)[vertexIndex].position().z());
      
      //   b) initial momentum
      GlobalVector momentum( (*theSimTracks)[simTrackId].momentum().x() , 
                 (*theSimTracks)[simTrackId].momentum().y() , 
                 (*theSimTracks)[simTrackId].momentum().z() );
      //   c) electric charge
      float        charge   = (*theSimTracks)[simTrackId].charge();
      //  -> inital parameters
      GlobalTrajectoryParameters initialParams(position,momentum,(int)charge,theMagField);
      //  -> large initial errors
      AlgebraicSymMatrix55 errorMatrix= AlgebraicMatrixID();      
      // errorMatrix = errorMatrix * 10;

      //this line help the fit succeed in the case of pixelless tracks (4th and 5th iteration)
      //for the future: probably the best thing is to use the mini-kalmanFilter
      if(theSeedHits0.subDetId() !=1 || theSeedHits0.subDetId() !=2) errorMatrix = errorMatrix * 0.0000001;



#ifdef FAMOS_DEBUG
      std::cout << "TrajectorySeedProducer: SimTrack parameters " << std::endl;
      std::cout << "\t\t pT  = " << (*theSimTracks)[simTrackId].momentum().Pt() << std::endl;
      std::cout << "\t\t eta = " << (*theSimTracks)[simTrackId].momentum().Eta()  << std::endl;
      std::cout << "\t\t phi = " << (*theSimTracks)[simTrackId].momentum().Phi()  << std::endl;
      std::cout << "TrajectorySeedProducer: AlgebraicSymMatrix " << errorMatrix << std::endl;
#endif
      CurvilinearTrajectoryError initialError(errorMatrix);
      // -> initial state
      FreeTrajectoryState initialFTS(initialParams, initialError);      
#ifdef FAMOS_DEBUG
      std::cout << "TrajectorySeedProducer: FTS momentum " << initialFTS.momentum() << std::endl;
#endif
      // const GeomDetUnit* initialLayer = theGeometry->idToDetUnit( recHits.front().geographicalId() );
      const GeomDet* initialLayer = theGeometry->idToDet( recHits.front().geographicalId() );

      //this is wrong because the FTS is defined at vertex, and it need to be properly propagated.
      //      const TrajectoryStateOnSurface initialTSOS(initialFTS, initialLayer->surface());      

      const TrajectoryStateOnSurface initialTSOS = thePropagator->propagate(initialFTS,initialLayer->surface()) ;
      if (!initialTSOS.isValid()) continue; 

#ifdef FAMOS_DEBUG
      std::cout << "TrajectorySeedProducer: TSOS global momentum "    << initialTSOS.globalMomentum() << std::endl;
      std::cout << "\t\t\tpT = "                                      << initialTSOS.globalMomentum().perp() << std::endl;
      std::cout << "\t\t\teta = "                                     << initialTSOS.globalMomentum().eta() << std::endl;
      std::cout << "\t\t\tphi = "                                     << initialTSOS.globalMomentum().phi() << std::endl;
      std::cout << "TrajectorySeedProducer: TSOS local momentum "     << initialTSOS.localMomentum()  << std::endl;
      std::cout << "TrajectorySeedProducer: TSOS local error "        << initialTSOS.localError().positionError() << std::endl;
      std::cout << "TrajectorySeedProducer: TSOS local error matrix " << initialTSOS.localError().matrix() << std::endl;
      std::cout << "TrajectorySeedProducer: TSOS surface side "       << initialTSOS.surfaceSide()    << std::endl;
#endif
      stateOnDet(initialTSOS, 
         recHits.front().geographicalId().rawId(),
         initialState);
      // Create a new Trajectory Seed    
      output[ialgo]->push_back(TrajectorySeed(initialState, recHits, alongMomentum));
#ifdef FAMOS_DEBUG
      std::cout << "Trajectory seed created ! " << std::endl;
#endif
      break;
      // End of the loop over seeding algorithms
    }
    // End on the loop over simtracks
  }

  for ( unsigned ialgo=0; ialgo<seedingAlgo.size(); ++ialgo ) { 
    std::auto_ptr<TrajectorySeedCollection> p(output[ialgo]);
    e.put(p,seedingAlgo[ialgo]);
  }

}

Side TrajectorySeedProducer::setLayerSpecSide

(

const std::string &

layerSpecSide

)

const

void TrajectorySeedProducer::stateOnDet ( const TrajectoryStateOnSurface &  ts, unsigned int  detid, PTrajectoryStateOnDet &  pts  ) const

private

A mere copy (without memory leak) of an existing tracking method.

should check if corresponds to m

Definition at line 804 of file TrajectorySeedProducer.cc.

References cond::rpcobgas::detid, i, j, gen::k, TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), m, LocalTrajectoryError::matrix(), and TrajectoryStateOnSurface::surfaceSide().

Referenced by produce().

{

  const AlgebraicSymMatrix55& m = ts.localError().matrix();
  
  int dim = 5; 

  float localErrors[15];
  int k = 0;
  for (int i=0; i<dim; ++i) {
    for (int j=0; j<=i; ++j) {
      localErrors[k++] = m(i,j);
    }
  }
  int surfaceSide = static_cast<int>(ts.surfaceSide());

  pts = PTrajectoryStateOnDet( ts.localParameters(),
                   localErrors, detid,
                   surfaceSide);
}

Member Data Documentation

unsigned int TrajectorySeedProducer::absMinRecHits

private

Definition at line 93 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

edm::EDGetTokenT<reco::BeamSpot> TrajectorySeedProducer::beamSpotToken

private

Definition at line 121 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

std::vector<unsigned int> TrajectorySeedProducer::firstHitSubDetectorNumber

private

Definition at line 97 of file TrajectorySeedProducer.h.

Referenced by TrajectorySeedProducer().

std::vector< std::vector<unsigned int> > TrajectorySeedProducer::firstHitSubDetectors

private

Definition at line 100 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

edm::InputTag TrajectorySeedProducer::hitProducer

private

Definition at line 88 of file TrajectorySeedProducer.h.

Referenced by TrajectorySeedProducer().

std::vector<double> TrajectorySeedProducer::maxD0

private

Definition at line 85 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

std::vector<double> TrajectorySeedProducer::maxZ0

private

Definition at line 86 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

std::vector<unsigned> TrajectorySeedProducer::minRecHits

private

Definition at line 87 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

bool TrajectorySeedProducer::newSyntax

private

Definition at line 104 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

std::vector<unsigned int> TrajectorySeedProducer::numberOfHits

private

Definition at line 95 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

std::vector<double> TrajectorySeedProducer::originHalfLength

private

Definition at line 109 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and TrajectorySeedProducer().

std::vector<double> TrajectorySeedProducer::originpTMin

private

Definition at line 110 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and TrajectorySeedProducer().

std::vector<double> TrajectorySeedProducer::originRadius

private

Definition at line 108 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and TrajectorySeedProducer().

std::vector<edm::InputTag> TrajectorySeedProducer::primaryVertices

private

Definition at line 112 of file TrajectorySeedProducer.h.

Referenced by TrajectorySeedProducer().

std::vector<double> TrajectorySeedProducer::pTMin

private

Definition at line 84 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

edm::EDGetTokenT<SiTrackerGSMatchedRecHit2DCollection> TrajectorySeedProducer::recHitToken

private

Definition at line 124 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

std::vector<edm::EDGetTokenT<reco::VertexCollection> > TrajectorySeedProducer::recoVertexToken

private

Definition at line 125 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

bool TrajectorySeedProducer::rejectOverlaps

private

Definition at line 92 of file TrajectorySeedProducer.h.

std::vector<unsigned int> TrajectorySeedProducer::secondHitSubDetectorNumber

private

Definition at line 98 of file TrajectorySeedProducer.h.

Referenced by TrajectorySeedProducer().

std::vector< std::vector<unsigned int> > TrajectorySeedProducer::secondHitSubDetectors

private

Definition at line 101 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

bool TrajectorySeedProducer::seedCleaning

private

Definition at line 91 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and TrajectorySeedProducer().

std::vector<std::string> TrajectorySeedProducer::seedingAlgo

private

Definition at line 94 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

bool TrajectorySeedProducer::selectMuons

private

Definition at line 115 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

edm::EDGetTokenT<edm::SimTrackContainer> TrajectorySeedProducer::simTrackToken

private

Definition at line 122 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

edm::EDGetTokenT<edm::SimVertexContainer> TrajectorySeedProducer::simVertexToken

private

Definition at line 123 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

edm::InputTag TrajectorySeedProducer::theBeamSpot

private

Definition at line 89 of file TrajectorySeedProducer.h.

Referenced by TrajectorySeedProducer().

const MagneticFieldMap* TrajectorySeedProducer::theFieldMap

private

Definition at line 80 of file TrajectorySeedProducer.h.

Referenced by beginRun(), and compatibleWithBeamAxis().

const TrackerGeometry* TrajectorySeedProducer::theGeometry

private

Definition at line 81 of file TrajectorySeedProducer.h.

Referenced by beginRun(), and produce().

std::vector< std::vector<LayerSpec> > TrajectorySeedProducer::theLayersInSets

private

Definition at line 105 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

const MagneticField* TrajectorySeedProducer::theMagField

private

Definition at line 79 of file TrajectorySeedProducer.h.

Referenced by beginRun(), and produce().

PropagatorWithMaterial* TrajectorySeedProducer::thePropagator

private

Definition at line 82 of file TrajectorySeedProducer.h.

Referenced by beginRun(), produce(), and ~TrajectorySeedProducer().

std::vector<unsigned int> TrajectorySeedProducer::thirdHitSubDetectorNumber

private

Definition at line 99 of file TrajectorySeedProducer.h.

Referenced by TrajectorySeedProducer().

std::vector< std::vector<unsigned int> > TrajectorySeedProducer::thirdHitSubDetectors

private

Definition at line 102 of file TrajectorySeedProducer.h.

Referenced by produce(), and TrajectorySeedProducer().

std::vector<const reco::VertexCollection*> TrajectorySeedProducer::vertices

private

Definition at line 117 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and produce().

double TrajectorySeedProducer::x0

private

Definition at line 118 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and produce().

double TrajectorySeedProducer::y0

private

Definition at line 118 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and produce().

double TrajectorySeedProducer::z0

private

Definition at line 118 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and produce().

std::vector<double> TrajectorySeedProducer::zVertexConstraint

private

Definition at line 113 of file TrajectorySeedProducer.h.

Referenced by compatibleWithBeamAxis(), and TrajectorySeedProducer().