#include <CosmicHitPairGeneratorFromLayerPair.h>

Public Member Functions
	CosmicHitPairGeneratorFromLayerPair (const LayerWithHits inner, const LayerWithHits outer, const edm::EventSetup &iSetup)

void	hitPairs (const TrackingRegion &ar, OrderedHitPairs &ap, const edm::EventSetup &iSetup)

const LayerWithHits *	innerLayer () const

const LayerWithHits *	outerLayer () const

	~CosmicHitPairGeneratorFromLayerPair ()

Private Attributes
const DetLayer *	innerlay

const DetLayer *	outerlay

const LayerWithHits *	theInnerLayer

const LayerWithHits *	theOuterLayer

const TrackerGeometry *	trackerGeometry

const TransientTrackingRecHitBuilder *	TTRHbuilder

Detailed Description

Definition at line 40 of file CosmicHitPairGeneratorFromLayerPair.h.

Constructor & Destructor Documentation

CosmicHitPairGeneratorFromLayerPair::CosmicHitPairGeneratorFromLayerPair	(	const LayerWithHits *	inner,
		const LayerWithHits *	outer,
		const edm::EventSetup &	iSetup
	)

Definition at line 15 of file CosmicHitPairGeneratorFromLayerPair.cc.

References edm::EventSetup::get(), edm::ESHandle< T >::product(), trackingTruthProducer_cfi::tracker, and trackerGeometry.

   : TTRHbuilder(nullptr),trackerGeometry(nullptr),
     //theLayerCache(*layerCache), 
     theOuterLayer(outer), theInnerLayer(inner)
 {
 
   edm::ESHandle<TrackerGeometry> tracker;
   iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
   trackerGeometry = tracker.product();
 }

CosmicHitPairGeneratorFromLayerPair::~CosmicHitPairGeneratorFromLayerPair ( )

Definition at line 28 of file CosmicHitPairGeneratorFromLayerPair.cc.

28 {}

Member Function Documentation

void CosmicHitPairGeneratorFromLayerPair::hitPairs	(	const TrackingRegion &	ar,
		OrderedHitPairs &	ap,
		const edm::EventSetup &	iSetup
	)

Definition at line 29 of file CosmicHitPairGeneratorFromLayerPair.cc.

References funct::abs(), edm::EventSetup::get(), BaseTrackerRecHit::globalPosition(), LayerWithHits::layer(), LayerWithHits::recHits(), AlCaHLTBitMon_QueryRunRegistry::string, theInnerLayer, theOuterLayer, and PV3DBase< T, PVType, FrameType >::z().

 {
 //  static int NSee = 0; static int Ntry = 0; static int Nacc = 0;
 
 
   typedef OrderedHitPair::InnerRecHit InnerHit;
   typedef OrderedHitPair::OuterRecHit OuterHit;
 
 
   if (theInnerLayer->recHits().empty()) return;
 
   if (theOuterLayer->recHits().empty()) return;
   //  const DetLayer* innerlay=theOuterLayer->layer();
   // const BarrelDetLayer *pippo=dynamic_cast<const BarrelDetLayer*>(theOuterLayer->layer());
 
 
   // ************ Daniele
 
   const DetLayer* blay1;
   const DetLayer* blay2;
   blay1 = dynamic_cast<const BarrelDetLayer*>(theInnerLayer->layer());
   blay2 = dynamic_cast<const BarrelDetLayer*>(theOuterLayer->layer());
 
 
   bool seedfromoverlaps= false;
   bool InTheBarrel = false;
   bool InTheForward = false;
   if (blay1 && blay2) {
     InTheBarrel = true;
   }
   else  InTheForward = true;
 
   if (InTheBarrel){
     float radius1 =dynamic_cast<const BarrelDetLayer*>(theInnerLayer->layer())->specificSurface().radius();
     float radius2 =dynamic_cast<const BarrelDetLayer*>(theOuterLayer->layer())->specificSurface().radius();
      seedfromoverlaps=(abs(radius1-radius2)<0.1) ? true : false;
   }
 
  
   vector<OrderedHitPair> allthepairs;
   std::string builderName = "WithTrackAngle";
   edm::ESHandle<TransientTrackingRecHitBuilder> builder;
   iSetup.get<TransientRecHitRecord>().get(builderName, builder);
 
   
 
   for( auto ohh=theOuterLayer->recHits().begin();ohh!=theOuterLayer->recHits().end();ohh++){
     for(auto ihh=theInnerLayer->recHits().begin();ihh!=theInnerLayer->recHits().end();ihh++){
      auto oh = static_cast<BaseTrackerRecHit const * const>(*ohh);
      auto ih = static_cast<BaseTrackerRecHit const * const>(*ihh);
       
       float z_diff =ih->globalPosition().z()-oh->globalPosition().z();
       float inny=ih->globalPosition().y();
       float outy=oh->globalPosition().y();
       float innx=ih->globalPosition().x();
       float outx=oh->globalPosition().x();;
       float dxdy=abs((outx-innx)/(outy-inny));
       float DeltaR=oh->globalPosition().perp()-ih->globalPosition().perp();;
       
       if( InTheBarrel && (abs(z_diff)<30) // && (outy > 0.) && (inny > 0.)
       //&&((abs(inny-outy))<30) 
       &&(dxdy<2)
       &&(inny*outy>0)
       && (abs(DeltaR)>0)) {
 
     //  cout << " ******** sono dentro inthebarrel *********** " << endl;
     if (seedfromoverlaps){
       //this part of code works for MTCC
       // for the other geometries must be verified
       //Overlaps in the difference in z is decreased and the difference in phi is
       //less than 0.05
       if ((DeltaR<0)&&(abs(z_diff)<18)&&(abs(ih->globalPosition().phi()-oh->globalPosition().phi())<0.05)&&(dxdy<2)) result.push_back( OrderedHitPair(ih, oh));
     }
     else  result.push_back( OrderedHitPair( ih, oh));
       
     
 
 
       }
       if( InTheForward &&  (abs(z_diff) > 1.)) {
     //  cout << " ******** sono dentro intheforward *********** " << endl;
     result.push_back( OrderedHitPair(ih,oh));
       }
     }
   }
 
 //   stable_sort(allthepairs.begin(),allthepairs.end(),CompareHitPairsY(iSetup));
 //   //Seed from overlaps are saved only if 
 //   //no others have been saved
 
 //   if (allthepairs.size()>0) {
 //     if (seedfromoverlaps) {
 //       if (result.size()==0) result.push_back(allthepairs[0]);
 //     }
 //     else result.push_back(allthepairs[0]);
 //   }
 
 
 }