---

DeDxTools Namespace Reference

Classes
struct	RawHits
Functions
Measurement1D	beta (const reco::TrackDeDxHits &trackWithHits)
	Default implementation simply returns the beta obtained inverting bethe bloch for the energy obtained with dedx().
float	betaFromInvertedBetheBloch (float dedx)
double	genericAverage (const reco::DeDxHitCollection &, float expo=1.)
float	massFromBeta (float beta, const reco::Track &)
void	trajectoryRawHits (const edm::Ref< std::vector< Trajectory > > &trajectory, vector< RawHits > &hits, bool usePixel, bool useStrip)
void	trajectoryRawHits (const edm::Ref< std::vector< Trajectory > > &trajectory, std::vector< RawHits > &hits, bool usePixel, bool useStrip)

---

Function Documentation

Measurement1D DeDxTools::beta

(

const reco::TrackDeDxHits &

trackWithHits

)

Default implementation simply returns the beta obtained inverting bethe bloch for the energy obtained with dedx().

The error is computed from bethe-bloch first derivative and ...? (Landau width)/sqrt(n) ??????

Referenced by SiStripRandomLorentzAngle::algoBeginJob(), PxCPEdbReader::beginJob(), DreamSD::cherenkovDeposit_(), MultipleScatteringUpdator::compute(), VolumeEnergyLossEstimator::computeBetheBloch(), EnergyLossUpdator::computeBetheBloch(), SiStripGainFromData::ComputeChargeOverPath(), HouseholderDecomposition::decompose(), GenericHouseholder::decompose(), DTSurveyConvert::endJob(), VolumeMultipleScatteringEstimator::estimate(), DDEcalBarrelAlgo::execute(), DTSurvey::FillWheelInfo(), PairProductionSimulator::gbteth(), BremsstrahlungSimulator::gbteth(), CastorSD::getEnergyDeposit(), ZdcSD::getEnergyDeposit(), NtupleManager::GetGlobalAngles(), HCalSD::getHitPMT(), TangentCircle::getPosition(), TopologyWorker::ludbrb(), main(), HcalRecHitsMaker::noiseInfCfromDB(), GflashEMShowerProfile::parameterization(), SiStripLAFakeESSource::produce(), SiPixelGaussianSmearingRecHitConverterAlgorithm::smearHit(), DDI::Polyhedra::volume(), DDEcalBarrelAlgo::web(), LMFLaserBluePrimDat::writeArrayDB(), and LMFLaserPrimDat::writeArrayDB().

float DeDxTools::betaFromInvertedBetheBloch

(

float

dedx

)

double DeDxTools::genericAverage	(	const reco::DeDxHitCollection &	hits,
		float	expo = 1.
	)

Definition at line 100 of file DeDxTools.cc.

References i, n, funct::pow(), and HLT_VtxMuL3::result.

Referenced by GenericAverageDeDxEstimator::dedx().

{
 double result=0;
 size_t n = hits.size();
 for(size_t i = 0; i< n; i ++)
 {
    result+=pow(hits[i].charge(),expo); 
 }
 return (n>0)?pow(result/n,1./expo):0.;
}

float DeDxTools::massFromBeta	(	float	beta,
		const reco::Track &
	)

void DeDxTools::trajectoryRawHits	(	const edm::Ref< std::vector< Trajectory > > &	trajectory,
		vector< RawHits > &	hits,
		bool	usePixel,
		bool	useStrip
	)

Definition at line 16 of file DeDxTools.cc.

References DeDxTools::RawHits::angleCosine, DeDxTools::RawHits::charge, SiStripRecHit2D::cluster(), DeDxTools::RawHits::detId, TrackingRecHit::geographicalId(), TrajectoryStateOnSurface::isValid(), it, TrajectoryStateOnSurface::localDirection(), PV3DBase< T, PVType, FrameType >::mag(), DeDxTools::RawHits::trajectoryMeasurement, and PV3DBase< T, PVType, FrameType >::z().

  {

    //    vector<RawHits> hits;

    const vector<TrajectoryMeasurement> & measurements = trajectory->measurements();
    for(vector<TrajectoryMeasurement>::const_iterator it = measurements.begin(); it!=measurements.end(); it++){

      //FIXME: check that "updated" State is the best one (wrt state in the middle) 
      TrajectoryStateOnSurface trajState=it->updatedState();
      if( !trajState.isValid()) continue;
     
      const TrackingRecHit * recHit=(*it->recHit()).hit();

       LocalVector trackDirection = trajState.localDirection();
       double cosine = trackDirection.z()/trackDirection.mag();
              
       if(const SiStripMatchedRecHit2D* matchedHit=dynamic_cast<const SiStripMatchedRecHit2D*>(recHit)){
           if(!useStrip) continue;

           RawHits mono,stereo; 
           mono.trajectoryMeasurement = &(*it);
           stereo.trajectoryMeasurement = &(*it);
           mono.angleCosine = cosine; 
           stereo.angleCosine = cosine;
           const std::vector<uint8_t> &  amplitudes = matchedHit->monoHit()->cluster()->amplitudes(); 
           mono.charge = accumulate(amplitudes.begin(), amplitudes.end(), 0);
       
           const std::vector<uint8_t> & amplitudesSt = matchedHit->stereoHit()->cluster()->amplitudes();
           stereo.charge = accumulate(amplitudesSt.begin(), amplitudesSt.end(), 0);
           
           mono.detId= matchedHit->monoHit()->geographicalId();
           stereo.detId= matchedHit->stereoHit()->geographicalId();

           hits.push_back(mono);
           hits.push_back(stereo);

        }else if(const ProjectedSiStripRecHit2D* projectedHit=dynamic_cast<const ProjectedSiStripRecHit2D*>(recHit)) {
           if(!useStrip) continue;

           const SiStripRecHit2D* singleHit=&(projectedHit->originalHit());
           RawHits mono;

           mono.trajectoryMeasurement = &(*it);

           mono.angleCosine = cosine; 
           const std::vector<uint8_t> & amplitudes = singleHit->cluster()->amplitudes();
           mono.charge = accumulate(amplitudes.begin(), amplitudes.end(), 0);
           mono.detId= singleHit->geographicalId();
           hits.push_back(mono);
      
        }else if(const SiStripRecHit2D* singleHit=dynamic_cast<const SiStripRecHit2D*>(recHit)){
           if(!useStrip) continue;

           RawHits mono;
               
           mono.trajectoryMeasurement = &(*it);

           mono.angleCosine = cosine; 
           const std::vector<uint8_t> & amplitudes = singleHit->cluster()->amplitudes();
           mono.charge = accumulate(amplitudes.begin(), amplitudes.end(), 0);
           mono.detId= singleHit->geographicalId();
           hits.push_back(mono);
      
        }else if(const SiPixelRecHit* pixelHit=dynamic_cast<const SiPixelRecHit*>(recHit)){
           if(!usePixel) continue;

           RawHits pixel;

           pixel.trajectoryMeasurement = &(*it);

           pixel.angleCosine = cosine; 
           pixel.charge = pixelHit->cluster()->charge();;
           pixel.detId= pixelHit->geographicalId();
           hits.push_back(pixel);
       }
       
    }
    // return hits;
  }

void DeDxTools::trajectoryRawHits	(	const edm::Ref< std::vector< Trajectory > > &	trajectory,
		std::vector< RawHits > &	hits,
		bool	usePixel,
		bool	useStrip
	)

Referenced by DeDxEstimatorProducer::produce().

---