ShallowGainCalibration Class Reference

#include <ShallowGainCalibration.h>

Inheritance diagram for ShallowGainCalibration:

Public Member Functions
	ShallowGainCalibration (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
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 ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
bool	IsFarFromBorder (TrajectoryStateOnSurface *trajState, const uint32_t detid, const edm::EventSetup *iSetup)
void	produce (edm::Event &, const edm::EventSetup &)
double	thickness (DetId id)

Private Attributes
std::map< DetId, double >	m_thicknessMap
const TrackerGeometry *	m_tracker
std::string	Prefix
std::string	Suffix
edm::InputTag	theTracksLabel

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
typedef WorkerT< EDProducer >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
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::EDProducer
CurrentProcessingContext const *	currentContext () const
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 59 of file ShallowGainCalibration.h.

Constructor & Destructor Documentation

ShallowGainCalibration::ShallowGainCalibration ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 12 of file ShallowGainCalibration.cc.

References Prefix, and Suffix.

  :  theTracksLabel( iConfig.getParameter<edm::InputTag>("Tracks") ),
     Suffix       ( iConfig.getParameter<std::string>("Suffix")    ),
     Prefix       ( iConfig.getParameter<std::string>("Prefix") )
{
  produces <std::vector<int> >            ( Prefix + "trackindex"     + Suffix );
  produces <std::vector<unsigned int> >   ( Prefix + "rawid"          + Suffix );
  produces <std::vector<float> >          ( Prefix + "localdirx"      + Suffix );
  produces <std::vector<float> >          ( Prefix + "localdiry"      + Suffix );
  produces <std::vector<float> >          ( Prefix + "localdirz"      + Suffix );
  produces <std::vector<unsigned short> > ( Prefix + "firststrip"     + Suffix );
  produces <std::vector<unsigned short> > ( Prefix + "nstrips"        + Suffix );
  produces <std::vector<bool> >           ( Prefix + "saturation"     + Suffix );
  produces <std::vector<bool> >           ( Prefix + "overlapping"    + Suffix );
  produces <std::vector<bool> >           ( Prefix + "farfromedge"    + Suffix );
  produces <std::vector<unsigned int> >   ( Prefix + "charge"         + Suffix );
  produces <std::vector<float> >          ( Prefix + "path"           + Suffix );
  produces <std::vector<float> >          ( Prefix + "chargeoverpath" + Suffix );
  produces <std::vector<unsigned char> >  ( Prefix + "amplitude"      + Suffix );
  produces <std::vector<double> >         ( Prefix + "gainused"       + Suffix );
}

Member Function Documentation

bool ShallowGainCalibration::IsFarFromBorder ( TrajectoryStateOnSurface *  trajState, const uint32_t  detid, const edm::EventSetup *  iSetup  )

private

Definition at line 227 of file ShallowGainCalibration.cc.

References Surface::bounds(), gather_cfg::cout, edm::EventSetup::get(), Bounds::length(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), Parameters::parameters, LocalTrajectoryError::positionError(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by produce().

{ 
  edm::ESHandle<TrackerGeometry> tkGeom; iSetup->get<TrackerDigiGeometryRecord>().get( tkGeom );

  LocalPoint  HitLocalPos   = trajState->localPosition();
  LocalError  HitLocalError = trajState->localError().positionError() ;

  const GeomDetUnit* it = tkGeom->idToDetUnit(DetId(detid));
  if (dynamic_cast<const StripGeomDetUnit*>(it)==0 && dynamic_cast<const PixelGeomDetUnit*>(it)==0) {
     std::cout << "this detID doesn't seem to belong to the Tracker" << std::endl;
     return false;
  }

  const BoundPlane plane = it->surface();
  const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
  const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));

  double DistFromBorder = 1.0;    
  double HalfLength     = it->surface().bounds().length() /2.0;

  if(trapezoidalBounds)
  {
      std::array<const float, 4> const & parameters = (*trapezoidalBounds).parameters();
     HalfLength     = parameters[3];
  }else if(rectangularBounds){
     HalfLength     = it->surface().bounds().length() /2.0;
  }else{return false;}

  if (fabs(HitLocalPos.y())+HitLocalError.yy() >= (HalfLength - DistFromBorder) ) return false;

  return true;
}

void ShallowGainCalibration::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Implements edm::EDProducer.

Definition at line 35 of file ShallowGainCalibration.cc.

References a, SiStripCluster::amplitudes(), DeDxDiscriminatorTools::charge(), shallow::findTrackIndex(), SiStripCluster::firstStrip(), SiStripCluster::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), h, IsFarFromBorder(), edm::ESHandleBase::isValid(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), m_tracker, PV3DBase< T, PVType, FrameType >::mag(), Trajectory::measurements(), NULL, getHLTPrescaleColumns::path, Prefix, edm::Event::put(), Suffix, theTracksLabel, thickness(), testEve_cfg::tracks, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

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

                                                       {
  std::auto_ptr<std::vector<int> >            trackindex    ( new std::vector<int>            );
  std::auto_ptr<std::vector<unsigned int> >   rawid         ( new std::vector<unsigned int>   );
  std::auto_ptr<std::vector<float>  >         localdirx     ( new std::vector<float>          );
  std::auto_ptr<std::vector<float>  >         localdiry     ( new std::vector<float>          );
  std::auto_ptr<std::vector<float>  >         localdirz     ( new std::vector<float>          );
  std::auto_ptr<std::vector<unsigned short> > firststrip    ( new std::vector<unsigned short> );
  std::auto_ptr<std::vector<unsigned short> > nstrips       ( new std::vector<unsigned short> );
  std::auto_ptr<std::vector<bool> >           saturation    ( new std::vector<bool>           );
  std::auto_ptr<std::vector<bool> >           overlapping   ( new std::vector<bool>           );
  std::auto_ptr<std::vector<bool> >           farfromedge   ( new std::vector<bool>           );
  std::auto_ptr<std::vector<unsigned int> >   charge        ( new std::vector<unsigned int>   );
  std::auto_ptr<std::vector<float>  >         path          ( new std::vector<float>          );
  std::auto_ptr<std::vector<float>  >         chargeoverpath( new std::vector<float>          );
  std::auto_ptr<std::vector<unsigned char> >  amplitude     ( new std::vector<unsigned char>  );
  std::auto_ptr<std::vector<double>  >        gainused      ( new std::vector<double>          );

  edm::ESHandle<TrackerGeometry> theTrackerGeometry;         iSetup.get<TrackerDigiGeometryRecord>().get( theTrackerGeometry );  
  m_tracker=&(* theTrackerGeometry );
  edm::ESHandle<SiStripGain> gainHandle;                     iSetup.get<SiStripGainRcd>().get(gainHandle);
  edm::Handle<edm::View<reco::Track> > tracks;               iEvent.getByLabel(theTracksLabel, tracks);   
  edm::Handle<TrajTrackAssociationCollection> associations;  iEvent.getByLabel(theTracksLabel, associations);

  for( TrajTrackAssociationCollection::const_iterator association = associations->begin(); association != associations->end(); association++) {
       const Trajectory*  traj  = association->key.get();
       const reco::Track* track = association->val.get();

       vector<TrajectoryMeasurement> measurements = traj->measurements();
       for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
          TrajectoryStateOnSurface trajState = measurement_it->updatedState();
          if( !trajState.isValid() ) continue;     

          const TrackingRecHit*         hit                 = (*measurement_it->recHit()).hit();
          const SiStripRecHit1D*        sistripsimple1dhit  = dynamic_cast<const SiStripRecHit1D*>(hit);
          const SiStripRecHit2D*        sistripsimplehit    = dynamic_cast<const SiStripRecHit2D*>(hit);
          const SiStripMatchedRecHit2D* sistripmatchedhit   = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);

          const SiStripCluster*   Cluster = NULL;
          uint32_t                DetId = 0;

          for(unsigned int h=0;h<2;h++){
            if(!sistripmatchedhit && h==1){
           continue;
            }else if(sistripmatchedhit  && h==0){
               Cluster = &sistripmatchedhit->monoCluster();
           DetId = sistripmatchedhit->monoId();
            }else if(sistripmatchedhit  && h==1){
               Cluster = &sistripmatchedhit->stereoCluster();;
           DetId = sistripmatchedhit->stereoId();
            }else if(sistripsimplehit){
               Cluster = (sistripsimplehit->cluster()).get();
           DetId = sistripsimplehit->geographicalId().rawId();
            }else if(sistripsimple1dhit){
               Cluster = (sistripsimple1dhit->cluster()).get();
           DetId = sistripsimple1dhit->geographicalId().rawId();
            }else{
               continue;
            }

            LocalVector             trackDirection = trajState.localDirection();
            double                  cosine         = trackDirection.z()/trackDirection.mag();
            const vector<uint8_t>&  Ampls          = Cluster->amplitudes();
        int                     FirstStrip     = Cluster->firstStrip();
            int                     APVId          = FirstStrip/128;
            bool                    Saturation     = false;
            bool                    Overlapping    = false;
            unsigned int            Charge         = 0;
            double                  Path           = (10.0*thickness(DetId))/fabs(cosine);
            double                  PrevGain       = -1;

            if(gainHandle.isValid()){ 
               SiStripApvGain::Range detGainRange = gainHandle->getRange(DetId);
               PrevGain = *(detGainRange.first + APVId);
            }

            for(unsigned int a=0;a<Ampls.size();a++){               
               Charge+=Ampls[a];
               if(Ampls[a] >=254)Saturation =true;
               amplitude->push_back( Ampls[a] );
            }
            double                   ChargeOverPath = (double)Charge / Path ;

            if(FirstStrip==0                                  )Overlapping=true;
            if(FirstStrip==128                                )Overlapping=true;
            if(FirstStrip==256                                )Overlapping=true;
            if(FirstStrip==384                                )Overlapping=true;
            if(FirstStrip==512                                )Overlapping=true;
            if(FirstStrip==640                                )Overlapping=true;

            if(FirstStrip<=127 && FirstStrip+Ampls.size()>127)Overlapping=true;
            if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlapping=true;
            if(FirstStrip<=383 && FirstStrip+Ampls.size()>383)Overlapping=true;
            if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlapping=true;
            if(FirstStrip<=639 && FirstStrip+Ampls.size()>639)Overlapping=true;

            if(FirstStrip+Ampls.size()==127                   )Overlapping=true;
            if(FirstStrip+Ampls.size()==255                   )Overlapping=true;
            if(FirstStrip+Ampls.size()==383                   )Overlapping=true;
            if(FirstStrip+Ampls.size()==511                   )Overlapping=true;
            if(FirstStrip+Ampls.size()==639                   )Overlapping=true;
            if(FirstStrip+Ampls.size()==767                   )Overlapping=true;

            trackindex    ->push_back( shallow::findTrackIndex(tracks, track) ); 
            rawid         ->push_back( DetId );         
            localdirx     ->push_back( trackDirection.x() );
            localdiry     ->push_back( trackDirection.y() );
            localdirz     ->push_back( trackDirection.z() );
            firststrip    ->push_back( FirstStrip );
            nstrips       ->push_back( Ampls.size() );
            saturation    ->push_back( Saturation );
            overlapping   ->push_back( Overlapping );
            farfromedge   ->push_back( IsFarFromBorder(&trajState,DetId, &iSetup) );
            charge        ->push_back( Charge );
            path          ->push_back( Path );
            chargeoverpath->push_back( ChargeOverPath );
            gainused      ->push_back( PrevGain );  
          }
       }
  }

  iEvent.put(trackindex,    Prefix + "trackindex"    + Suffix );
  iEvent.put(rawid     ,    Prefix + "rawid"         + Suffix );
  iEvent.put(localdirx ,    Prefix + "localdirx"     + Suffix );
  iEvent.put(localdiry ,    Prefix + "localdiry"     + Suffix );
  iEvent.put(localdirz ,    Prefix + "localdirz"     + Suffix );
  iEvent.put(firststrip,    Prefix + "firststrip"    + Suffix );
  iEvent.put(nstrips,       Prefix + "nstrips"       + Suffix );
  iEvent.put(saturation,    Prefix + "saturation"    + Suffix );
  iEvent.put(overlapping,   Prefix + "overlapping"   + Suffix );
  iEvent.put(farfromedge,   Prefix + "farfromedge"   + Suffix );
  iEvent.put(charge,        Prefix + "charge"        + Suffix );
  iEvent.put(path,          Prefix + "path"          + Suffix );
  iEvent.put(chargeoverpath,Prefix + "chargeoverpath"+ Suffix );
  iEvent.put(amplitude,     Prefix + "amplitude"     + Suffix );
  iEvent.put(gainused,      Prefix + "gainused"      + Suffix );
}

double ShallowGainCalibration::thickness ( DetId  id )

private

Definition at line 261 of file ShallowGainCalibration.cc.

References TrackerGeometry::idToDetUnit(), m_thicknessMap, and m_tracker.

Referenced by produce().

{
 map<DetId,double>::iterator th=m_thicknessMap.find(id);
 if(th!=m_thicknessMap.end())
   return (*th).second;
 else {
   double detThickness=1.;
   //compute thickness normalization
   const GeomDetUnit* it = m_tracker->idToDetUnit(DetId(id));
   bool isPixel = dynamic_cast<const PixelGeomDetUnit*>(it)!=0;
   bool isStrip = dynamic_cast<const StripGeomDetUnit*>(it)!=0;
   if (!isPixel && ! isStrip) {
   //FIXME throw exception
      edm::LogWarning("DeDxHitsProducer") << "\t\t this detID doesn't seem to belong to the Tracker";
      detThickness = 1.;
  }else{
      detThickness = it->surface().bounds().thickness();
  }

   m_thicknessMap[id]=detThickness;//computed value
   return detThickness;
 }

}

Member Data Documentation

std::map<DetId,double> ShallowGainCalibration::m_thicknessMap

private

Definition at line 74 of file ShallowGainCalibration.h.

Referenced by thickness().

const TrackerGeometry* ShallowGainCalibration::m_tracker

private

Definition at line 73 of file ShallowGainCalibration.h.

Referenced by produce(), and thickness().

std::string ShallowGainCalibration::Prefix

private

Definition at line 65 of file ShallowGainCalibration.h.

Referenced by produce(), and ShallowGainCalibration().

std::string ShallowGainCalibration::Suffix

private

Definition at line 64 of file ShallowGainCalibration.h.

Referenced by produce(), and ShallowGainCalibration().

edm::InputTag ShallowGainCalibration::theTracksLabel

private

Definition at line 63 of file ShallowGainCalibration.h.

Referenced by produce().