Phase2OTtiltedBarrelLayer.cc

Go to the documentation of this file.

#include "Phase2OTtiltedBarrelLayer.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DetGroupMerger.h"
#include "CompatibleDetToGroupAdder.h"

#include "TrackingTools/DetLayers/interface/DetLayerException.h"
#include "TrackingTools/DetLayers/interface/MeasurementEstimator.h"
#include "TrackingTools/DetLayers/interface/CylinderBuilderFromDet.h"
#include "Phase2EndcapLayerBuilder.h"

using namespace std;

typedef GeometricSearchDet::DetWithState DetWithState;

Phase2OTtiltedBarrelLayer::Phase2OTtiltedBarrelLayer(std::vector<const Phase2OTBarrelRod*>& innerRods,
                                         std::vector<const Phase2OTBarrelRod*>& outerRods,
                                                     vector<const Phase2EndcapRing*>& negRings,
                                                     vector<const Phase2EndcapRing*>& posRings) :
  Phase2OTBarrelLayer(innerRods,outerRods),
  theNegativeRingsComps(negRings.begin(),negRings.end()),
  thePositiveRingsComps(posRings.begin(),posRings.end())
{

  std::vector<const GeometricSearchDet*> theComps;
  theComps.assign(Phase2OTBarrelLayer::components().begin(),Phase2OTBarrelLayer::components().end());
  theComps.insert(theComps.end(),negRings.begin(),negRings.end());
  theComps.insert(theComps.end(),posRings.begin(),posRings.end());

  for(vector<const GeometricSearchDet*>::const_iterator it=theComps.begin();
      it!=theComps.end();it++){  
    theBasicComps.insert(theBasicComps.end(),   
             (**it).basicComponents().begin(),
             (**it).basicComponents().end());
  }

  BarrelDetLayer::initialize();
  theCylinder = cylinder( theComps );

#ifdef EDM_ML_DEBUG
  //--------- DEBUG INFO --------------
  LogDebug("TkDetLayers") << "==== DEBUG Phase2OTtiltedBarrelLayer =====" ; 
  LogTrace("TkDetLayers") << "Phase2OTtiltedBarrelLayer Cyl r,lenght: "
                          << theCylinder->radius() << " , "
                          << theCylinder->bounds().length();

  for (vector<const GeometricSearchDet*>::const_iterator i=theNegativeRingsComps.begin();
       i != theNegativeRingsComps.end(); i++){
    LogTrace("TkDetLayers") << "negative rings in Phase2OT tilted barrel pos z,perp,eta,phi: " 
                << (**i).position().z()    << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta()  << " , " 
                << (**i).position().phi()  ;
  }
  
  for (vector<const GeometricSearchDet*>::const_iterator i=Phase2OTBarrelLayer::components().begin();
       i != Phase2OTBarrelLayer::components().end(); i++){
    LogTrace("TkDetLayers") << "rods in Phase2OT tilted barrel pos z,perp,eta,phi: " 
                << (**i).position().z()    << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta()  << " , " 
                << (**i).position().phi()  ;
  }

  for (vector<const GeometricSearchDet*>::const_iterator i=thePositiveRingsComps.begin();
       i != thePositiveRingsComps.end(); i++){
    LogTrace("TkDetLayers") << "positive rings in Phase2OT tilted barrel pos z,perp,eta,phi: " 
                << (**i).position().z()    << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta()  << " , " 
                << (**i).position().phi()  ;
  }
  LogTrace("TkDetLayers") << "==== end DEBUG Phase2OTtiltedBarrelLayer =====" ; 
  //----------------------------------- 
#endif
}

Phase2OTtiltedBarrelLayer::~Phase2OTtiltedBarrelLayer(){

  vector<const GeometricSearchDet*>::const_iterator i;
  for (i=theNegativeRingsComps.begin(); i!=theNegativeRingsComps.end(); i++) {
    delete *i;
  }
  for (i=thePositiveRingsComps.begin(); i!=thePositiveRingsComps.end(); i++) {
    delete *i;
  }

} 

void 
Phase2OTtiltedBarrelLayer::groupedCompatibleDetsV( const TrajectoryStateOnSurface& tsos,
                      const Propagator& prop,
                       const MeasurementEstimator& est,
                       std::vector<DetGroup> & result) const {
  vector<DetGroup> closestResultRods;
  vector<DetGroup> closestResultNeg;
  vector<DetGroup> closestResultPos;
  Phase2OTBarrelLayer::groupedCompatibleDetsV(tsos, prop, est, closestResultRods);
  if(tsos.globalPosition().z()<0){
    for(auto& ring : theNegativeRingsComps){
      ring->groupedCompatibleDetsV(tsos, prop, est, closestResultNeg);
    }
  } else {
    for(auto& ring : thePositiveRingsComps){
      ring->groupedCompatibleDetsV(tsos, prop, est, closestResultPos);
    }
  }
  result.assign(closestResultRods.begin(),closestResultRods.end());
  result.insert(result.end(),closestResultPos.begin(),closestResultPos.end());
  result.insert(result.end(),closestResultNeg.begin(),closestResultNeg.end());

#ifdef EDM_ML_DEBUG
  LogDebug("TkDetLayers") << "==== output di Phase2OTtiltedBarrelLayer =====" ; 
  if(closestResultRods.size() != 0){
    for (auto gr : closestResultRods) {
      LogTrace("TkDetLayers") << "New Rod group:";
      for (auto dge : gr) {
        LogTrace("TkDetLayers") << "new det with geom det at r:"<<dge.det()->position().perp()<<" id:"<<dge.det()->geographicalId().rawId()<<" tsos at:" <<dge.trajectoryState().globalPosition();
      }
    }
  }
  if(closestResultNeg.size() != 0){
    for (auto gr : closestResultNeg) {
      LogTrace("TkDetLayers") << "New negative group:";
      for (auto dge : gr) {
        LogTrace("TkDetLayers") << "new det with geom det at r:"<<dge.det()->position().perp()<<" id:"<<dge.det()->geographicalId().rawId()<<" tsos at:" <<dge.trajectoryState().globalPosition();
      }
    }
  }
  if(closestResultPos.size() != 0){
    for (auto gr : closestResultPos) {
      LogTrace("TkDetLayers") << "New positive group:";
      for (auto dge : gr) {
        LogTrace("TkDetLayers") << "new det with geom det at r:"<<dge.det()->position().perp()<<" id:"<<dge.det()->geographicalId().rawId()<<" tsos at:" <<dge.trajectoryState().globalPosition();
      }
    }
  }

  if(result.size() != 0){
    for (auto gr : result) {
      LogTrace("TkDetLayers") << "Total group:";
      for (auto dge : gr) {
        LogTrace("TkDetLayers") << "new det with geom det at r:"<<dge.det()->position().perp()<<" id:"<<dge.det()->geographicalId().rawId()<<" tsos at:" <<dge.trajectoryState().globalPosition();
      }
    }
  } else {
      LogTrace("TkDetLayers") << "result size is zero"; 
  }
#endif
  
}