#include <HitEff.h>

Inheritance diagram for HitEff:

Public Member Functions
bool	check2DPartner (unsigned int iidd, const std::vector< TrajectoryMeasurement > &traj)

double	checkConsistency (const StripClusterParameterEstimator::LocalValues &parameters, double xx, double xerr)

unsigned int	checkLayer (unsigned int iidd, const TrackerTopology *tTopo)

	HitEff (const edm::ParameterSet &conf)

bool	isDoubleSided (unsigned int iidd, const TrackerTopology *tTopo) const

	~HitEff () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	beginJob () override

void	endJob () override

Private Attributes
bool	addCommonMode_

bool	addLumi_

unsigned int	bunchx

float	ClusterLocErrX

float	ClusterLocErrY

float	ClusterLocX

float	ClusterLocY

const edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	clusters_token_

float	ClusterStoN

const edm::EDGetTokenT< reco::TrackCollection >	combinatorialTracks_token_

float	commonMode

const edm::EDGetTokenT< edm::DetSetVector< SiStripRawDigi > >	commonModeToken_

int	compSettings

edm::ParameterSet	conf_

bool	cutOnTracks_

bool	DEBUG

const edm::EDGetTokenT< DetIdCollection >	digis_token_

unsigned int	event

int	events

int	EventTrackCKF

bool	highPurity

unsigned int	Id

float	instLumi

unsigned int	layers

unsigned int	ModIsBad

int	nHits

float	pT

float	PU

float	ResX

float	ResXSig

unsigned int	run

const edm::EDGetTokenT< LumiScalersCollection >	scalerToken_

unsigned int	SiStripQualBad

int	tquality

const edm::EDGetTokenT< MeasurementTrackerEvent >	trackerEvent_token_

unsigned int	trackMultiplicityCut_

TTree *	traj

const edm::EDGetTokenT< std::vector< Trajectory > >	trajectories_token_

float	TrajGlbX

float	TrajGlbY

float	TrajGlbZ

unsigned int	trajHitValid

float	TrajLocAngleX

float	TrajLocAngleY

float	TrajLocErrX

float	TrajLocErrY

float	TrajLocX

float	TrajLocY

const edm::EDGetTokenT< TrajTrackAssociationCollection >	trajTrackAsso_token_

bool	useAllHitsFromTracksWithMissingHits_

bool	useFirstMeas_

bool	useLastMeas_

unsigned int	whatlayer

bool	withinAcceptance

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 48 of file HitEff.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 74 of file HitEff.cc.

References addCommonMode_, addLumi_, compSettings, conf_, cutOnTracks_, DEBUG, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), layers, trackMultiplicityCut_, useAllHitsFromTracksWithMissingHits_, useFirstMeas_, and useLastMeas_.

                                           : 
   scalerToken_( consumes< LumiScalersCollection >(conf.getParameter<edm::InputTag>("lumiScalers")) ),
   commonModeToken_( mayConsume< edm::DetSetVector<SiStripRawDigi> >(conf.getParameter<edm::InputTag>("commonMode")) ),
   combinatorialTracks_token_( consumes< reco::TrackCollection >(conf.getParameter<edm::InputTag>("combinatorialTracks")) ),
   trajectories_token_( consumes< std::vector<Trajectory> >(conf.getParameter<edm::InputTag>("trajectories")) ),
   trajTrackAsso_token_( consumes< TrajTrackAssociationCollection >(conf.getParameter<edm::InputTag>("trajectories")) ),
   clusters_token_( consumes< edmNew::DetSetVector<SiStripCluster> >(conf.getParameter<edm::InputTag>("siStripClusters")) ),
   digis_token_( consumes< DetIdCollection >(conf.getParameter<edm::InputTag>("siStripDigis")) ),
   trackerEvent_token_( consumes< MeasurementTrackerEvent>(conf.getParameter<edm::InputTag>("trackerEvent")) ),
   conf_(conf)
 {
   compSettings=conf_.getUntrackedParameter<int>("CompressionSettings",-1);
   layers =conf_.getParameter<int>("Layer");
   DEBUG = conf_.getParameter<bool>("Debug");
   addLumi_ = conf_.getUntrackedParameter<bool>("addLumi", false);
   addCommonMode_ = conf_.getUntrackedParameter<bool>("addCommonMode", false);
   cutOnTracks_ = conf_.getUntrackedParameter<bool>("cutOnTracks", false);
   trackMultiplicityCut_ = conf.getUntrackedParameter<unsigned int>("trackMultiplicity",100);
   useFirstMeas_ = conf_.getUntrackedParameter<bool>("useFirstMeas", false);
   useLastMeas_ = conf_.getUntrackedParameter<bool>("useLastMeas", false);
   useAllHitsFromTracksWithMissingHits_ = conf_.getUntrackedParameter<bool>("useAllHitsFromTracksWithMissingHits", false);
 }

HitEff::~HitEff ( )

override

Definition at line 98 of file HitEff.cc.

98 { }

Member Function Documentation

void HitEff::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprivate

Definition at line 161 of file HitEff.cc.

References funct::abs(), addCommonMode_, addLumi_, AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, edm::AssociationMap< Tag >::begin(), edmNew::DetSetVector< T >::begin(), Surface::bounds(), edm::EventBase::bunchCrossing(), bunchx, check2DPartner(), checkConsistency(), checkLayer(), vertices_cff::chi2, ClusterLocErrX, ClusterLocErrY, ClusterLocX, ClusterLocY, clusters_token_, ClusterStoN, combinatorialTracks_token_, commonMode, commonModeToken_, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, cutOnTracks_, digis_token_, edmNew::DetSetVector< T >::empty(), edm::AssociationMap< Tag >::end(), edmNew::DetSetVector< T >::end(), edm::EventID::event(), events, EventTrackCKF, edm::EventSetup::get(), SiStripQuality::getBadApvs(), edm::Event::getByLabel(), edm::Event::getByToken(), highPurity, muons_cff::highPurity, mps_fire::i, edm::EventBase::id(), Id, TrackerGeometry::idToDetUnit(), cuy::ii, input, instLumi, ires, isDoubleSided(), edm::HandleBase::isValid(), edm::Ref< C, T, F >::key(), layers, StripClusterParameterEstimator::localParameters(), LogDebug, LayerMeasurements::measurements(), ecalDrivenElectronSeedsParameters_cff::measurementTrackerEvent, ModIsBad, reco::MuonTime::nDof, nHits, StripTopology::nstrips(), AlCaHLTBitMon_ParallelJobs::p, TrapezoidalPlaneBounds::parameters(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, pT, PU, rpcPointValidation_cfi::recHit, TrajectoryMeasurement::recHit(), ResX, ResXSig, edm::EventID::run(), run, scalerToken_, SiStripQualBad, edm::EDCollection< T >::size(), edm::DetSetVector< T >::size(), edmNew::DetSetVector< T >::size(), StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), subdetector, GeomDet::surface(), StripSubdetector::TEC, TrackerTopology::tecSide(), reco::MuonTime::timeAtIpInOut, reco::MuonTime::timeAtIpInOutErr, tmp, trackingTruthProducer_cfi::tracker, trackerEvent_token_, trackMultiplicityCut_, traj, trajectories_token_, TrajGlbX, TrajGlbY, TrajGlbZ, trajHitValid, TrajLocAngleX, TrajLocAngleY, TrajLocErrX, TrajLocErrY, TrajLocX, TrajLocY, trajTrackAsso_token_, useAllHitsFromTracksWithMissingHits_, useFirstMeas_, useLastMeas_, whatlayer, Bounds::width(), and withinAcceptance.

                                                               {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   //  bool DEBUG = false;
 
   LogDebug("SiStripHitEfficiency:HitEff")  << "beginning analyze from HitEff" << endl;
 
   using namespace edm;
   using namespace reco;
   // Step A: Get Inputs 
 
   int run_nr = e.id().run();
   int ev_nr = e.id().event();
   int bunch_nr = e.bunchCrossing();
 
   // Luminosity informations
   edm::Handle<LumiScalersCollection> lumiScalers;
   instLumi=0; PU=0;
   if(addLumi_) {
     e.getByToken(scalerToken_, lumiScalers); 
     if (lumiScalers->begin() != lumiScalers->end()) {
       instLumi = lumiScalers->begin()->instantLumi();
       PU = lumiScalers->begin()->pileup();
     }
   }
   
   // CM
   edm::Handle<edm::DetSetVector<SiStripRawDigi> > commonModeDigis;
   if(addCommonMode_) e.getByToken(commonModeToken_, commonModeDigis);
 
   //CombinatoriaTrack
   edm::Handle<reco::TrackCollection> trackCollectionCKF;
   //edm::InputTag TkTagCKF = conf_.getParameter<edm::InputTag>("combinatorialTracks");
   e.getByToken(combinatorialTracks_token_,trackCollectionCKF);
   
   edm::Handle<std::vector<Trajectory> > TrajectoryCollectionCKF;
   //edm::InputTag TkTrajCKF = conf_.getParameter<edm::InputTag>("trajectories");
   e.getByToken(trajectories_token_,TrajectoryCollectionCKF);
   
   edm::Handle<TrajTrackAssociationCollection> trajTrackAssociationHandle;
   e.getByToken(trajTrackAsso_token_, trajTrackAssociationHandle);
   
   // Clusters
   // get the SiStripClusters from the event
   edm::Handle< edmNew::DetSetVector<SiStripCluster> > theClusters;
   //e.getByLabel("siStripClusters", theClusters);
   e.getByToken(clusters_token_, theClusters);
 
   //get tracker geometry
   edm::ESHandle<TrackerGeometry> tracker;
   es.get<TrackerDigiGeometryRecord>().get(tracker);
   const TrackerGeometry * tkgeom=&(* tracker);
 
   //get Cluster Parameter Estimator
   //std::string cpe = conf_.getParameter<std::string>("StripCPE");
   edm::ESHandle<StripClusterParameterEstimator> parameterestimator;
   es.get<TkStripCPERecord>().get("StripCPEfromTrackAngle", parameterestimator); 
   const StripClusterParameterEstimator &stripcpe(*parameterestimator);
 
   // get the SiStripQuality records
   edm::ESHandle<SiStripQuality> SiStripQuality_;
 //LQ commenting the try/catch that causes problem in 74X calibTree production
 //  try {
 //    es.get<SiStripQualityRcd>().get("forCluster",SiStripQuality_);
 //  }
 //  catch (...) {
     es.get<SiStripQualityRcd>().get(SiStripQuality_);
 //  }
   
   edm::ESHandle<MagneticField> magFieldHandle;
   es.get<IdealMagneticFieldRecord>().get(magFieldHandle);
   const MagneticField* magField_ = magFieldHandle.product();
 
   // get the list of module IDs with FED-detected errors
   edm::Handle< DetIdCollection > fedErrorIds;
   //e.getByLabel("siStripDigis", fedErrorIds );
   e.getByToken(digis_token_, fedErrorIds );
 
   ESHandle<MeasurementTracker> measurementTrackerHandle;
   es.get<CkfComponentsRecord>().get(measurementTrackerHandle);
 
   edm::Handle<MeasurementTrackerEvent> measurementTrackerEvent;
   //e.getByLabel("MeasurementTrackerEvent", measurementTrackerEvent);
   e.getByToken(trackerEvent_token_,measurementTrackerEvent);
 
   edm::ESHandle<Chi2MeasurementEstimatorBase> est;
   es.get<TrackingComponentsRecord>().get("Chi2",est);
 
   edm::ESHandle<Propagator> prop;
   es.get<TrackingComponentsRecord>().get("PropagatorWithMaterial",prop);
   const Propagator* thePropagator = prop.product();
 
   events++;
   
   // *************** SiStripCluster Collection
   const edmNew::DetSetVector<SiStripCluster>& input = *theClusters;
 
   //go through clusters to write out global position of good clusters for the layer understudy for comparison
   // Loop through clusters just to print out locations
   // Commented out to avoid discussion, should really be deleted.
   /*
   for (edmNew::DetSetVector<SiStripCluster>::const_iterator DSViter = input.begin(); DSViter != input.end(); DSViter++) {
     // DSViter is a vector of SiStripClusters located on a single module
     unsigned int ClusterId = DSViter->id();
     DetId ClusterDetId(ClusterId);
     const StripGeomDetUnit * stripdet=(const StripGeomDetUnit*)tkgeom->idToDetUnit(ClusterDetId);
     
     edmNew::DetSet<SiStripCluster>::const_iterator begin=DSViter->begin();
     edmNew::DetSet<SiStripCluster>::const_iterator end  =DSViter->end();
     for(edmNew::DetSet<SiStripCluster>::const_iterator iter=begin;iter!=end;++iter) {
       //iter is a single SiStripCluster
       StripClusterParameterEstimator::LocalValues parameters=stripcpe.localParameters(*iter,*stripdet);
       
       const Surface* surface;
       surface = &(tracker->idToDet(ClusterDetId)->surface());
       LocalPoint lp = parameters.first;
       GlobalPoint gp = surface->toGlobal(lp);
       unsigned int layer = checkLayer(ClusterId, tTopo);
             if(DEBUG) cout << "Found hit in cluster collection layer = " << layer << " with id = " << ClusterId << "   local X position = " << lp.x() << " +- " << sqrt(parameters.second.xx()) << "   matched/stereo/rphi = " << ((ClusterId & 0x3)==0) << "/" << ((ClusterId & 0x3)==1) << "/" << ((ClusterId & 0x3)==2) << endl;
     }
   }
   */
   
   // Tracking 
   const   reco::TrackCollection *tracksCKF=trackCollectionCKF.product();
   LogDebug("SiStripHitEfficiency:HitEff") << "number ckf tracks found = " << tracksCKF->size() << endl;
   //if (tracksCKF->size() == 1 ){
   if (!tracksCKF->empty()) {
     if( cutOnTracks_ && (tracksCKF->size() >= trackMultiplicityCut_) ) return;
     if( cutOnTracks_ ) LogDebug("SiStripHitEfficiency:HitEff") << "starting checking good event with < "<< trackMultiplicityCut_ <<" tracks" << endl;
 
     EventTrackCKF++;  
 
 #ifdef ExtendedCALIBTree
     //get dEdx info if available
     Handle<ValueMap<DeDxData> >          dEdxUncalibHandle;
     if (e.getByLabel("dedxMedianCTF", dEdxUncalibHandle)) {
       const ValueMap<DeDxData> dEdxTrackUncalib = *dEdxUncalibHandle.product();
       
       reco::TrackRef itTrack  = reco::TrackRef( trackCollectionCKF, 0 );
       dedx = dEdxTrackUncalib[itTrack].dEdx();
       dedxNOM  = dEdxTrackUncalib[itTrack].numberOfMeasurements();
     } else {
       dedx = -999.0; dedxNOM = -999;
     }
 
     
 
     //get muon and ecal timing info if available
     Handle<MuonCollection> muH;
     if(e.getByLabel("muonsWitht0Correction",muH)){
       const MuonCollection & muonsT0  =  *muH.product();
       if(!muonsT0.empty()) {
     MuonTime mt0 = muonsT0[0].time();
     timeDT = mt0.timeAtIpInOut; 
     timeDTErr = mt0.timeAtIpInOutErr;
     timeDTDOF = mt0.nDof;
     
     bool hasCaloEnergyInfo = muonsT0[0].isEnergyValid();
     if (hasCaloEnergyInfo) timeECAL = muonsT0[0].calEnergy().ecal_time;
       }
     } else {
       timeDT = -999.0; timeDTErr = -999.0; timeDTDOF = -999; timeECAL = -999.0;
     }
 
 #endif
     // actually should do a loop over all the tracks in the event here
 
     // Looping over traj-track associations to be able to get traj & track informations
     for(TrajTrackAssociationCollection::const_iterator it = trajTrackAssociationHandle->begin();
      it!=trajTrackAssociationHandle->end();
      it++) {
      
       edm::Ref<std::vector<Trajectory> > itraj  = it->key;
       reco::TrackRef itrack  = it->val;
 
       // for each track, fill some variables such as number of hits and momentum
       nHits = itraj->foundHits();
 #ifdef ExtendedCALIBTree      
       nLostHits = itraj->lostHits();
       chi2 = (itraj->chiSquared()/itraj->ndof());
       p = itraj->lastMeasurement().updatedState().globalMomentum().mag();
 #endif
       pT = sqrt( ( itraj->lastMeasurement().updatedState().globalMomentum().x() *
            itraj->lastMeasurement().updatedState().globalMomentum().x()) +
          ( itraj->lastMeasurement().updatedState().globalMomentum().y() *
            itraj->lastMeasurement().updatedState().globalMomentum().y()) );
       
       // track quality
       highPurity = itrack->quality(reco::TrackBase::TrackQuality::highPurity);
 
       
       
       std::vector<TrajectoryMeasurement> TMeas=itraj->measurements();
       vector<TrajectoryMeasurement>::iterator itm;
       double xloc = 0.;
       double yloc = 0.;
       double xErr = 0.;
       double yErr = 0.;
       double angleX = -999.;
       double angleY = -999.;
       double xglob,yglob,zglob;
 
 
       // Check whether the trajectory has some missing hits
       bool hasMissingHits=false;
       for (itm=TMeas.begin();itm!=TMeas.end();itm++){
         auto theHit = (*itm).recHit();
         if(theHit->getType()==TrackingRecHit::Type::missing) hasMissingHits=true;  
       }
     
     
     // Loop on each measurement and take it into consideration
     //--------------------------------------------------------
     
           for (itm=TMeas.begin();itm!=TMeas.end();itm++){
     auto theInHit = (*itm).recHit();
     
     LogDebug("SiStripHitEfficiency:HitEff") << "theInHit is valid = " << theInHit->isValid() << endl;
     
     unsigned int iidd = theInHit->geographicalId().rawId();
     
     unsigned int TKlayers = checkLayer(iidd, tTopo);
     LogDebug("SiStripHitEfficiency:HitEff") << "TKlayer from trajectory: " << TKlayers << "  from module = " << iidd <<  "   matched/stereo/rphi = " << ((iidd & 0x3)==0) << "/" << ((iidd & 0x3)==1) << "/" << ((iidd & 0x3)==2) << endl;
 
     
     // Test first and last points of the trajectory
     // the list of measurements starts from outer layers  !!! This could change -> should add a check
     bool isFirstMeas = (itm==(TMeas.end()-1));
     bool isLastMeas = (itm==(TMeas.begin()));
     
     if(!useFirstMeas_ && isFirstMeas) continue;
     if(!useLastMeas_ && isLastMeas) continue;
     
     // In case of missing hit in the track, check whether to use the other hits or not.
     if(hasMissingHits && theInHit->getType()!=TrackingRecHit::Type::missing && !useAllHitsFromTracksWithMissingHits_) continue;
 
     
         
     // If Trajectory measurement from TOB 6 or TEC 9, skip it because it's always valid they are filled later
     if ( TKlayers == 10 || TKlayers == 22 ) {
       LogDebug("SiStripHitEfficiency:HitEff") << "skipping original TM for TOB 6 or TEC 9" << endl;
       continue;
     }
 
     // Make vector of TrajectoryAtInvalidHits to hold the trajectories
     std::vector<TrajectoryAtInvalidHit> TMs;
     
     // Make AnalyticalPropagator to use in TAVH constructor
     AnalyticalPropagator propagator(magField_,anyDirection); 
     
     // for double sided layers check both sensors--if no hit was found on either sensor surface,
     // the trajectory measurements only have one invalid hit entry on the matched surface
     // so get the TrajectoryAtInvalidHit for both surfaces and include them in the study
     if (isDoubleSided(iidd, tTopo) &&  ((iidd & 0x3)==0) ) {
       // do hit eff check twice--once for each sensor
       //add a TM for each surface
       TMs.push_back(TrajectoryAtInvalidHit(*itm, tTopo, tkgeom, propagator, 1));
       TMs.push_back(TrajectoryAtInvalidHit(*itm, tTopo, tkgeom, propagator, 2));
     } else if ( isDoubleSided(iidd, tTopo) && (!check2DPartner(iidd, TMeas)) ) {
       // if only one hit was found the trajectory measurement is on that sensor surface, and the other surface from
       // the matched layer should be added to the study as well
       TMs.push_back(TrajectoryAtInvalidHit(*itm, tTopo, tkgeom, propagator, 1));
       TMs.push_back(TrajectoryAtInvalidHit(*itm, tTopo, tkgeom, propagator, 2));
       LogDebug("SiStripHitEfficiency:HitEff") << " found a hit with a missing partner" << endl;
     } else {
       //only add one TM for the single surface and the other will be added in the next iteration
       TMs.push_back(TrajectoryAtInvalidHit(*itm, tTopo, tkgeom, propagator));
     }
 
     //Now check for tracks at TOB6 and TEC9
 
     // to make sure we only propagate on the last TOB5 hit check the next entry isn't also in TOB5
     // to avoid bias, make sure the TOB5 hit is valid (an invalid hit on TOB5 could only exist with a valid hit on TOB6)
 
     bool isValid = theInHit->isValid();
     bool isLast = (itm==(TMeas.end()-1));
     bool isLastTOB5 = true;
     if (!isLast) {
       if ( checkLayer((++itm)->recHit()->geographicalId().rawId(), tTopo) == 9 ) isLastTOB5 = false;
       else isLastTOB5 = true;
       --itm;
     }
     
     if ( TKlayers==9 && isValid && isLastTOB5 ) {
       //      if ( TKlayers==9 && itm==TMeas.rbegin()) {
     //    if ( TKlayers==9 && (itm==TMeas.back()) ) {     // to check for only the last entry in the trajectory for propagation
       std::vector< BarrelDetLayer const*> barrelTOBLayers = measurementTrackerHandle->geometricSearchTracker()->tobLayers() ;
       const DetLayer* tob6 = barrelTOBLayers[barrelTOBLayers.size()-1];
       const MeasurementEstimator* estimator = est.product();
       const LayerMeasurements* theLayerMeasurements = new LayerMeasurements(*measurementTrackerHandle, *measurementTrackerEvent);
       const TrajectoryStateOnSurface tsosTOB5 = itm->updatedState();
       vector<TrajectoryMeasurement> tmp = theLayerMeasurements->measurements(*tob6, tsosTOB5, *thePropagator, *estimator);
       
       if ( !tmp.empty()) {
         LogDebug("SiStripHitEfficiency:HitEff") << "size of TM from propagation = " << tmp.size() << endl;
 
         // take the last of the TMs, which is always an invalid hit
         // if no detId is available, ie detId==0, then no compatible layer was crossed
         // otherwise, use that TM for the efficiency measurement
         TrajectoryMeasurement tob6TM(tmp.back());
         const auto& tob6Hit = tob6TM.recHit();
         
         if (tob6Hit->geographicalId().rawId()!=0) {
           LogDebug("SiStripHitEfficiency:HitEff") << "tob6 hit actually being added to TM vector" << endl;
           TMs.push_back(TrajectoryAtInvalidHit(tob6TM, tTopo, tkgeom, propagator));
         }
       }
     }
 
     bool isLastTEC8 = true;
     if (!isLast) {
       if ( checkLayer((++itm)->recHit()->geographicalId().rawId(), tTopo) == 21 ) isLastTEC8 = false;
       else isLastTEC8 = true;
       --itm;
     }
     
     if ( TKlayers==21 && isValid && isLastTEC8 ) {
       
       std::vector< const ForwardDetLayer*> posTecLayers = measurementTrackerHandle->geometricSearchTracker()->posTecLayers() ;
       const DetLayer* tec9pos = posTecLayers[posTecLayers.size()-1];
       std::vector< const ForwardDetLayer*> negTecLayers = measurementTrackerHandle->geometricSearchTracker()->negTecLayers() ;
       const DetLayer* tec9neg = negTecLayers[negTecLayers.size()-1];
       
       const MeasurementEstimator* estimator = est.product();
       const LayerMeasurements* theLayerMeasurements = new LayerMeasurements(*measurementTrackerHandle, *measurementTrackerEvent);
       const TrajectoryStateOnSurface tsosTEC9 = itm->updatedState();
       
       // check if track on positive or negative z
       if (!(iidd ==  StripSubdetector::TEC)) LogDebug("SiStripHitEfficiency:HitEff") << "there is a problem with TEC 9 extrapolation" << endl;
       
       //cout << " tec9 id = " << iidd << " and side = " << tTopo->tecSide(iidd) << endl;
       vector<TrajectoryMeasurement> tmp;
       if ( tTopo->tecSide(iidd) == 1 ) {
         tmp = theLayerMeasurements->measurements(*tec9neg, tsosTEC9, *thePropagator, *estimator);
         //cout << "on negative side" << endl;
       }
       if ( tTopo->tecSide(iidd) == 2 ) {
         tmp = theLayerMeasurements->measurements(*tec9pos, tsosTEC9, *thePropagator, *estimator);
         //cout << "on positive side" << endl;
       }
 
       if ( !tmp.empty()) {
         // take the last of the TMs, which is always an invalid hit
         // if no detId is available, ie detId==0, then no compatible layer was crossed
         // otherwise, use that TM for the efficiency measurement
         TrajectoryMeasurement tec9TM(tmp.back());
         const auto& tec9Hit = tec9TM.recHit();
         
         unsigned int tec9id = tec9Hit->geographicalId().rawId();
         LogDebug("SiStripHitEfficiency:HitEff") << "tec9id = " << tec9id << " is Double sided = " <<  isDoubleSided(tec9id, tTopo) << "  and 0x3 = " << (tec9id & 0x3) << endl;
         
         if (tec9Hit->geographicalId().rawId()!=0) {
           LogDebug("SiStripHitEfficiency:HitEff") << "tec9 hit actually being added to TM vector" << endl;
           // in tec the hit can be single or doubled sided. whenever the invalid hit at the end of vector of TMs is
           // double sided it is always on the matched surface, so we need to split it into the true sensor surfaces
           if (isDoubleSided(tec9id, tTopo)) {
         TMs.push_back(TrajectoryAtInvalidHit(tec9TM, tTopo, tkgeom, propagator, 1));
         TMs.push_back(TrajectoryAtInvalidHit(tec9TM, tTopo, tkgeom, propagator, 2));
           } else 
         TMs.push_back(TrajectoryAtInvalidHit(tec9TM, tTopo, tkgeom, propagator));
         }
       } //else cout << "tec9 tmp empty" << endl;
     }
     
     
     // Modules Constraints
 
     for(std::vector<TrajectoryAtInvalidHit>::const_iterator TM=TMs.begin();TM!=TMs.end();++TM) {
       
       // --> Get trajectory from combinatedState 
       iidd = TM->monodet_id();
       LogDebug("SiStripHitEfficiency:HitEff") << "setting iidd = " << iidd << " before checking efficiency and ";
       
       xloc = TM->localX();
       yloc = TM->localY();
       
       angleX = atan( TM->localDxDz() );
       angleY = atan( TM->localDyDz() );
 
       TrajLocErrX = 0.0; TrajLocErrY = 0.0;
 
       xglob = TM->globalX();
       yglob = TM->globalY();
       zglob = TM->globalZ();
       xErr =  TM->localErrorX();
       yErr =  TM->localErrorY();
 
       TrajGlbX = 0.0; TrajGlbY = 0.0; TrajGlbZ = 0.0;
       withinAcceptance = TM->withinAcceptance();
       
       trajHitValid = TM->validHit();
       int TrajStrip = -1;
 
       // reget layer from iidd here, to account for TOB 6 and TEC 9 TKlayers being off
       TKlayers = checkLayer(iidd, tTopo);
 
       if ( (layers == TKlayers) || (layers == 0) ) {   // Look at the layer not used to reconstruct the track
         whatlayer = TKlayers;
         LogDebug("SiStripHitEfficiency:HitEff") << "Looking at layer under study" << endl;
         ModIsBad = 2; Id = 0; SiStripQualBad = 0; 
         run = 0; event = 0; TrajLocX = 0.0; TrajLocY = 0.0; 
         TrajLocAngleX = -999.0; TrajLocAngleY = -999.0; ResX = 0.0; ResXSig = 0.0;
         ClusterLocX = 0.0; ClusterLocY = 0.0; ClusterLocErrX = 0.0; ClusterLocErrY = 0.0; ClusterStoN = 0.0;
         bunchx = 0; commonMode = -100;
         
         // RPhi RecHit Efficiency 
         
         if (!input.empty() ) {  
           LogDebug("SiStripHitEfficiency:HitEff") << "Checking clusters with size = " << input.size() << endl;
           int nClusters = 0;
           std::vector< std::vector<float> > VCluster_info; //fill with X residual, X residual pull, local X, sig(X), local Y, sig(Y), StoN
           for (edmNew::DetSetVector<SiStripCluster>::const_iterator DSViter = input.begin(); DSViter != input.end(); DSViter++) {
         // DSViter is a vector of SiStripClusters located on a single module
         //if (DEBUG)      cout << "the ID from the DSViter = " << DSViter->id() << endl; 
         unsigned int ClusterId = DSViter->id();
         if (ClusterId == iidd) {
           LogDebug("SiStripHitEfficiency:HitEff") << "found  (ClusterId == iidd) with ClusterId = " << ClusterId << " and iidd = " << iidd << endl;
           DetId ClusterDetId(ClusterId);
           const StripGeomDetUnit * stripdet=(const StripGeomDetUnit*)tkgeom->idToDetUnit(ClusterDetId);
           const StripTopology& Topo  = stripdet->specificTopology();
           
                   float hbedge   = 0.0;
                   float htedge   = 0.0;
                   float hapoth   = 0.0;
                   float uylfac   = 0.0;
                   float uxlden   = 0.0;
                   if(TKlayers>=11) {
                      const BoundPlane& plane = stripdet->surface();
                      const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
              std::array<const float, 4> const & parameterTrap = (*trapezoidalBounds).parameters(); // el bueno aqui
                      hbedge         = parameterTrap[0];
                      htedge         = parameterTrap[1];
                      hapoth         = parameterTrap[3];
                      uylfac   = (htedge-hbedge)/(htedge+hbedge)/hapoth;
                      uxlden   = 1 + yloc*uylfac;
                   }
           
           // Need to know position of trajectory in strip number for selecting the right APV later
           if( TrajStrip==-1 ) {
             int nstrips = Topo.nstrips();
             float pitch = stripdet->surface().bounds().width() / nstrips;
             TrajStrip = xloc/pitch + nstrips/2.0;
             // Need additionnal corrections for endcap
             if(TKlayers>=11) {
               float TrajLocXMid = xloc / (1 + (htedge-hbedge)*yloc/(htedge+hbedge)/hapoth) ; // radialy extrapolated x loc position at middle
               TrajStrip = TrajLocXMid/pitch + nstrips/2.0;
             }
             //cout<<" Layer "<<TKlayers<<" TrajStrip: "<<nstrips<<" "<<pitch<<" "<<TrajStrip<<endl;
           }
 
 
           for(edmNew::DetSet<SiStripCluster>::const_iterator iter=DSViter->begin();iter!=DSViter->end();++iter) {
             //iter is a single SiStripCluster
             StripClusterParameterEstimator::LocalValues parameters=stripcpe.localParameters(*iter,*stripdet);
             float res = (parameters.first.x() - xloc);
             float sigma = checkConsistency(parameters , xloc, xErr);
             // The consistency is probably more accurately measured with the Chi2MeasurementEstimator. To use it
             // you need a TransientTrackingRecHit instead of the cluster
             //theEstimator=       new Chi2MeasurementEstimator(30);
             //const Chi2MeasurementEstimator *theEstimator(100);
             //theEstimator->estimate(TM->tsos(), TransientTrackingRecHit);
 
                     if(TKlayers>=11) {
                        res = parameters.first.x() - xloc/uxlden ; // radialy extrapolated x loc position at middle           
                        sigma = abs(res) / sqrt(parameters.second.xx() + xErr*xErr/uxlden/uxlden + yErr*yErr*xloc*xloc*uylfac*uylfac/uxlden/uxlden/uxlden/uxlden);
                     }
             
             SiStripClusterInfo clusterInfo = SiStripClusterInfo(*iter, es, ClusterId);  
             // signal to noise from SiStripClusterInfo not working in 225. I'll fix this after the interface
             // redesign in 300 -ku
             //float cluster_info[7] = {res, sigma, parameters.first.x(), sqrt(parameters.second.xx()), parameters.first.y(), sqrt(parameters.second.yy()), signal_to_noise};
             std::vector< float > cluster_info;
             cluster_info.push_back(res); 
             cluster_info.push_back(sigma);
             cluster_info.push_back(parameters.first.x()); 
             cluster_info.push_back(sqrt(parameters.second.xx()));
             cluster_info.push_back(parameters.first.y());
             cluster_info.push_back(sqrt(parameters.second.yy()));
             cluster_info.push_back( clusterInfo.signalOverNoise() );
             //cluster_info.push_back( clusterInfo.getSignalOverNoise() );
             VCluster_info.push_back(cluster_info);
             nClusters++;
             LogDebug("SiStripHitEfficiency:HitEff") << "Have ID match. residual = " << VCluster_info.back()[0] << "  res sigma = " << VCluster_info.back()[1] << endl;
             LogDebug("SiStripHitEfficiency:HitEff") << "trajectory measurement compatability estimate = " << (*itm).estimate() << endl;
             LogDebug("SiStripHitEfficiency:HitEff") << "hit position = " << parameters.first.x() << "  hit error = " << sqrt(parameters.second.xx()) << "  trajectory position = " << xloc << "  traj error = " << xErr << endl;
           }
         }
           }
           float FinalResSig = 1000.0;
           float FinalCluster[7]= {1000.0, 1000.0, 0.0, 0.0, 0.0, 0.0, 0.0};
           if (nClusters > 0) {
         LogDebug("SiStripHitEfficiency:HitEff") << "found clusters > 0" << endl;
         if (nClusters > 1) {
           //get the smallest one
           vector< vector<float> >::iterator ires;
           for (ires=VCluster_info.begin(); ires!=VCluster_info.end(); ires++){
             if ( abs((*ires)[1]) < abs(FinalResSig)) {
               FinalResSig = (*ires)[1];
               for (unsigned int i = 0; i<ires->size(); i++) {
             LogDebug("SiStripHitEfficiency:HitEff") << "filling final cluster. i = " << i << " before fill FinalCluster[i]=" << FinalCluster[i] << " and (*ires)[i] =" << (*ires)[i] << endl;
             FinalCluster[i] = (*ires)[i];
             LogDebug("SiStripHitEfficiency:HitEff") << "filling final cluster. i = " << i << " after fill FinalCluster[i]=" << FinalCluster[i] << " and (*ires)[i] =" << (*ires)[i] << endl;
               }
             }
             LogDebug("SiStripHitEfficiency:HitEff") << "iresidual = " << (*ires)[0] << "  isigma = " << (*ires)[1] << "  and FinalRes = " << FinalCluster[0] << endl;
           }
         }
         else {
           FinalResSig = VCluster_info.at(0)[1];
           for (unsigned int i = 0; i<VCluster_info.at(0).size(); i++) {
             FinalCluster[i] = VCluster_info.at(0)[i];
           }
         }
         nClusters=0;
         VCluster_info.clear();
           }
           
           LogDebug("SiStripHitEfficiency:HitEff") << "Final residual in X = " << FinalCluster[0] << "+-" << (FinalCluster[0]/FinalResSig) << endl;
           LogDebug("SiStripHitEfficiency:HitEff") << "Checking location of trajectory: abs(yloc) = " << abs(yloc) << "  abs(xloc) = " << abs(xloc) << endl;
           LogDebug("SiStripHitEfficiency:HitEff") << "Checking location of cluster hit: yloc = " << FinalCluster[4] << "+-" << FinalCluster[5] << "  xloc = " << FinalCluster[2] << "+-" << FinalCluster[3] << endl;
           LogDebug("SiStripHitEfficiency:HitEff") << "Final cluster signal to noise = " << FinalCluster[6] << endl;
           
           float exclusionWidth = 0.4;
           float TOBexclusion = 0.0;
           float TECexRing5 = -0.89;
           float TECexRing6 = -0.56;
           float TECexRing7 = 0.60;
           //Added by Chris Edelmaier to do TEC bonding exclusion
           int subdetector = ((iidd>>25) & 0x7);            
           int ringnumber = ((iidd>>5) & 0x7);
           
           //New TOB and TEC bonding region exclusion zone
           if((TKlayers >= 5 && TKlayers < 11)||((subdetector == 6)&&( (ringnumber >= 5)&&(ringnumber <=7) ))) {
         //There are only 2 cases that we need to exclude for
         float highzone = 0.0;
         float lowzone = 0.0;
         float higherr = yloc + 5.0*yErr;
         float lowerr = yloc - 5.0*yErr;
         if(TKlayers >= 5 && TKlayers < 11) {
           //TOB zone
           highzone = TOBexclusion + exclusionWidth;
           lowzone = TOBexclusion - exclusionWidth;
         } else if (ringnumber == 5) {
           //TEC ring 5
           highzone = TECexRing5 + exclusionWidth;
           lowzone = TECexRing5 - exclusionWidth;
         } else if (ringnumber == 6) {
           //TEC ring 6
           highzone = TECexRing6 + exclusionWidth;
           lowzone = TECexRing6 - exclusionWidth;
         } else if (ringnumber == 7) {
           //TEC ring 7
           highzone = TECexRing7 + exclusionWidth;
           lowzone = TECexRing7 - exclusionWidth;
         }
         //Now that we have our exclusion region, we just have to properly identify it
         if((highzone <= higherr)&&(highzone >= lowerr)) withinAcceptance = false;
         if((lowzone >= lowerr)&&(lowzone <= higherr)) withinAcceptance = false;
         if((higherr <= highzone)&&(higherr >= lowzone)) withinAcceptance = false;
         if((lowerr >= lowzone)&&(lowerr <= highzone)) withinAcceptance = false;
           }
           
           // fill ntuple varibles
           //get global position from module id number iidd
           TrajGlbX = xglob;
           TrajGlbY = yglob;
           TrajGlbZ = zglob; 
   
           TrajLocErrX = xErr;
           TrajLocErrY = yErr;
 
           Id = iidd;
           run = run_nr;
           event = ev_nr;
           bunchx = bunch_nr;
           //if ( SiStripQuality_->IsModuleBad(iidd) ) {
           if ( SiStripQuality_->getBadApvs(iidd)!=0 ) {
         SiStripQualBad = 1; 
         LogDebug("SiStripHitEfficiency:HitEff") << "strip is bad from SiStripQuality" << endl;
           } else {
         SiStripQualBad = 0; 
         LogDebug("SiStripHitEfficiency:HitEff") << "strip is good from SiStripQuality" << endl;
           }
 
           //check for FED-detected errors and include those in SiStripQualBad
           for (unsigned int ii=0;ii< fedErrorIds->size();ii++) {
         if (iidd == (*fedErrorIds)[ii].rawId() )
           SiStripQualBad = 1;
           }
           
           TrajLocX = xloc;
           TrajLocY = yloc;
           TrajLocAngleX = angleX;
           TrajLocAngleY = angleY;
           ResX = FinalCluster[0];
           ResXSig = FinalResSig;
           if (FinalResSig != FinalCluster[1]) LogDebug("SiStripHitEfficiency:HitEff") << "Problem with best cluster selection because FinalResSig = " << FinalResSig << " and FinalCluster[1] = " << FinalCluster[1] << endl;
           ClusterLocX = FinalCluster[2];
           ClusterLocY = FinalCluster[4];
           ClusterLocErrX = FinalCluster[3];
           ClusterLocErrY = FinalCluster[5];
           ClusterStoN = FinalCluster[6];
           
           // CM of APV crossed by traj
           if(addCommonMode_)
           if(commonModeDigis.isValid() && TrajStrip>=0 && TrajStrip<=768) {
             edm::DetSetVector<SiStripRawDigi>::const_iterator digiframe = commonModeDigis->find(iidd);
             if(digiframe != commonModeDigis->end())
               if( (unsigned) TrajStrip/128 < digiframe->data.size())
                 commonMode = digiframe->data.at(TrajStrip/128).adc();
           }
 
           LogDebug("SiStripHitEfficiency:HitEff") << "before check good" << endl;
           
           if ( FinalResSig < 999.0) {  //could make requirement on track/hit consistency, but for
         //now take anything with a hit on the module
         LogDebug("SiStripHitEfficiency:HitEff") << "hit being counted as good " << FinalCluster[0] << " FinalRecHit " << 
           iidd << "   TKlayers  "  <<  TKlayers  << " xloc " <<  xloc << " yloc  " << yloc << " module " << iidd << 
           "   matched/stereo/rphi = " << ((iidd & 0x3)==0) << "/" << ((iidd & 0x3)==1) << "/" << ((iidd & 0x3)==2) << endl;
         ModIsBad = 0;
         traj->Fill();
           }
           else {
         LogDebug("SiStripHitEfficiency:HitEff") << "hit being counted as bad   ######### Invalid RPhi FinalResX " << FinalCluster[0] << " FinalRecHit " << 
           iidd << "   TKlayers  "  <<  TKlayers  << " xloc " <<  xloc << " yloc  " << yloc << " module " << iidd << 
           "   matched/stereo/rphi = " << ((iidd & 0x3)==0) << "/" << ((iidd & 0x3)==1) << "/" << ((iidd & 0x3)==2) << endl;
         ModIsBad = 1;
         traj->Fill();
         
         LogDebug("SiStripHitEfficiency:HitEff") << " RPhi Error " << sqrt(xErr*xErr + yErr*yErr) << " ErrorX " << xErr << " yErr " <<  yErr <<  endl;
           } LogDebug("SiStripHitEfficiency:HitEff") << "after good location check" << endl;
         } LogDebug("SiStripHitEfficiency:HitEff") << "after list of clusters" << endl;
       } LogDebug("SiStripHitEfficiency:HitEff") << "After layers=TKLayers if" << endl;
     } LogDebug("SiStripHitEfficiency:HitEff") << "After looping over TrajAtValidHit list" << endl;
       } LogDebug("SiStripHitEfficiency:HitEff") << "end TMeasurement loop" << endl;
     } LogDebug("SiStripHitEfficiency:HitEff") << "end of trajectories loop" << endl;
   }
 }

void HitEff::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 100 of file HitEff.cc.

References addCommonMode_, addLumi_, bunchx, vertices_cff::chi2, ClusterLocErrX, ClusterLocErrY, ClusterLocX, ClusterLocY, ClusterStoN, commonMode, compSettings, events, EventTrackCKF, TFileService::file(), highPurity, Id, instLumi, TFileService::make(), ModIsBad, nHits, AlCaHLTBitMon_ParallelJobs::p, pT, PU, ResX, ResXSig, run, SiStripQualBad, tquality, traj, TrajGlbX, TrajGlbY, TrajGlbZ, trajHitValid, TrajLocAngleX, TrajLocAngleY, TrajLocErrX, TrajLocErrY, TrajLocX, TrajLocY, whatlayer, and withinAcceptance.

                      {
 
   edm::Service<TFileService> fs;
   if(compSettings>0){
     edm::LogInfo("SiStripHitEfficiency:HitEff")<<"the compressions settings are:"<< compSettings << std::endl;
     fs->file().SetCompressionSettings(compSettings);
   }
 
   traj = fs->make<TTree>("traj","tree of trajectory positions");
   #ifdef ExtendedCALIBTree
   traj->Branch("timeDT",&timeDT,"timeDT/F");
   traj->Branch("timeDTErr",&timeDTErr,"timeDTErr/F");
   traj->Branch("timeDTDOF",&timeDTDOF,"timeDTDOF/I");
   traj->Branch("timeECAL",&timeECAL,"timeECAL/F");
   traj->Branch("dedx",&dedx,"dedx/F");
   traj->Branch("dedxNOM",&dedxNOM,"dedxNOM/I"); 
   traj->Branch("nLostHits",&nLostHits,"nLostHits/I");
   traj->Branch("chi2",&chi2,"chi2/F");
   traj->Branch("p",&p,"p/F");
   #endif
   traj->Branch("TrajGlbX",&TrajGlbX,"TrajGlbX/F");
   traj->Branch("TrajGlbY",&TrajGlbY,"TrajGlbY/F");
   traj->Branch("TrajGlbZ",&TrajGlbZ,"TrajGlbZ/F");
   traj->Branch("TrajLocX",&TrajLocX,"TrajLocX/F");
   traj->Branch("TrajLocY",&TrajLocY,"TrajLocY/F");
   traj->Branch("TrajLocAngleX",&TrajLocAngleX,"TrajLocAngleX/F");
   traj->Branch("TrajLocAngleY",&TrajLocAngleY,"TrajLocAngleY/F");
   traj->Branch("TrajLocErrX",&TrajLocErrX,"TrajLocErrX/F");
   traj->Branch("TrajLocErrY",&TrajLocErrY,"TrajLocErrY/F");
   traj->Branch("ClusterLocX",&ClusterLocX,"ClusterLocX/F");
   traj->Branch("ClusterLocY",&ClusterLocY,"ClusterLocY/F");
   traj->Branch("ClusterLocErrX",&ClusterLocErrX,"ClusterLocErrX/F");
   traj->Branch("ClusterLocErrY",&ClusterLocErrY,"ClusterLocErrY/F");
   traj->Branch("ClusterStoN",&ClusterStoN,"ClusterStoN/F");
   traj->Branch("ResX",&ResX,"ResX/F");
   traj->Branch("ResXSig",&ResXSig,"ResXSig/F");
   traj->Branch("ModIsBad",&ModIsBad,"ModIsBad/i");
   traj->Branch("SiStripQualBad",&SiStripQualBad,"SiStripQualBad/i");
   traj->Branch("withinAcceptance",&withinAcceptance,"withinAcceptance/O");
   traj->Branch("nHits",&nHits,"nHits/I");
   traj->Branch("pT",&pT,"pT/F");
   traj->Branch("highPurity",&highPurity,"highPurity/O");
   traj->Branch("trajHitValid", &trajHitValid, "trajHitValid/i");
   traj->Branch("Id",&Id,"Id/i");
   traj->Branch("run",&run,"run/i");
   traj->Branch("event",&event,"event/i");
   traj->Branch("layer",&whatlayer,"layer/i");
   traj->Branch("tquality",&tquality,"tquality/I");
   traj->Branch("bunchx",&bunchx,"bunchx/I");
   if(addLumi_) {
     traj->Branch("instLumi",&instLumi,"instLumi/F");
     traj->Branch("PU",&PU,"PU/F");
   }
   if(addCommonMode_) traj->Branch("commonMode",&commonMode,"commonMode/F");
 
   events = 0;
   EventTrackCKF = 0;
   
 }

bool HitEff::check2DPartner	(	unsigned int	iidd,
		const std::vector< TrajectoryMeasurement > &	traj
	)

Definition at line 853 of file HitEff.cc.

Referenced by analyze().

                                                                                            {
   unsigned int partner_iidd = 0;
   bool found2DPartner = false;
   // first get the id of the other detector
   if ((iidd & 0x3)==1) partner_iidd = iidd+1;
   if ((iidd & 0x3)==2) partner_iidd = iidd-1;
   // next look in the trajectory measurements for a measurement from that detector
   // loop through trajectory measurements to find the partner_iidd
   for (std::vector<TrajectoryMeasurement>::const_iterator iTM=traj.begin(); iTM!=traj.end(); ++iTM) {
     if (iTM->recHit()->geographicalId().rawId()==partner_iidd) {
       found2DPartner = true;
     }
   }
   return found2DPartner;
 }

double HitEff::checkConsistency	(	const StripClusterParameterEstimator::LocalValues &	parameters,
		double	xx,
		double	xerr
	)

Definition at line 814 of file HitEff.cc.

References funct::abs(), relativeConstraints::error, mathSSE::sqrt(), and geometryCSVtoXML::xx.

Referenced by analyze().

                                                                                                                    {
   double error = sqrt(parameters.second.xx() + xerr*xerr);
   double separation = abs(parameters.first.x() - xx);
   double consistency = separation/error;
   return consistency;
 }

unsigned int HitEff::checkLayer	(	unsigned int	iidd,
		const TrackerTopology *	tTopo
	)

Definition at line 869 of file HitEff.cc.

References DEFINE_FWK_MODULE, digitizers_cfi::strip, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecWheel(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidWheel(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

Referenced by analyze().

                                                                                 {
   StripSubdetector strip=StripSubdetector(iidd);
   unsigned int subid=strip.subdetId();
   if (subid ==  StripSubdetector::TIB) { 
     
     return tTopo->tibLayer(iidd);
   }
   if (subid ==  StripSubdetector::TOB) { 
     
     return tTopo->tobLayer(iidd) + 4 ; 
   }
   if (subid ==  StripSubdetector::TID) { 
     
     return tTopo->tidWheel(iidd) + 10;
   }
   if (subid ==  StripSubdetector::TEC) { 
     
     return tTopo->tecWheel(iidd) + 13 ; 
   }
   return 0;
 }

void HitEff::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 806 of file HitEff.cc.

References events, EventTrackCKF, LogDebug, and traj.

Referenced by o2olib.O2ORunMgr::executeJob().

                    {
   traj->GetDirectory()->cd();
   traj->Write();  
   
   LogDebug("SiStripHitEfficiency:HitEff") << " Events Analysed             " <<  events          << endl;
   LogDebug("SiStripHitEfficiency:HitEff") << " Number Of Tracked events    " <<  EventTrackCKF   << endl;
 }

bool HitEff::isDoubleSided	(	unsigned int	iidd,
		const TrackerTopology *	tTopo
	)		const

Definition at line 821 of file HitEff.cc.

References digitizers_cfi::strip, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecRing(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

Referenced by analyze().

                                                                                 {
   StripSubdetector strip=StripSubdetector(iidd);
   unsigned int subid=strip.subdetId();
   unsigned int layer = 0;
   if (subid ==  StripSubdetector::TIB) { 
     
     layer = tTopo->tibLayer(iidd);
     if (layer == 1 || layer == 2) return true;
     else return false;
   }
   else if (subid ==  StripSubdetector::TOB) { 
     
     layer = tTopo->tobLayer(iidd) + 4 ; 
     if (layer == 5 || layer == 6) return true;
     else return false;
   }
   else if (subid ==  StripSubdetector::TID) { 
     
     layer = tTopo->tidRing(iidd) + 10;
     if (layer == 11 || layer == 12) return true;
     else return false;
   }
   else if (subid ==  StripSubdetector::TEC) { 
     
     layer = tTopo->tecRing(iidd) + 13 ; 
     if (layer == 14 || layer == 15 || layer == 18) return true;
     else return false;
   }
   else
     return false;
 }

Member Data Documentation

bool HitEff::addCommonMode_

private

Definition at line 68 of file HitEff.h.

Referenced by analyze(), beginJob(), and HitEff().

bool HitEff::addLumi_

private

Definition at line 67 of file HitEff.h.

Referenced by analyze(), beginJob(), and HitEff().

unsigned int HitEff::bunchx

private

Definition at line 111 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::ClusterLocErrX

private

Definition at line 105 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::ClusterLocErrY

private

Definition at line 105 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::ClusterLocX

private

Definition at line 105 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::ClusterLocY

private

Definition at line 105 of file HitEff.h.

Referenced by analyze(), and beginJob().

const edm::EDGetTokenT< edmNew::DetSetVector<SiStripCluster> > HitEff::clusters_token_

private

Definition at line 78 of file HitEff.h.

Referenced by analyze().

float HitEff::ClusterStoN

private

Definition at line 105 of file HitEff.h.

Referenced by analyze(), and beginJob().

const edm::EDGetTokenT< reco::TrackCollection > HitEff::combinatorialTracks_token_

private

Definition at line 75 of file HitEff.h.

Referenced by analyze().

float HitEff::commonMode

private

Definition at line 114 of file HitEff.h.

Referenced by analyze(), and beginJob().

const edm::EDGetTokenT<edm::DetSetVector<SiStripRawDigi> > HitEff::commonModeToken_

private

Definition at line 65 of file HitEff.h.

Referenced by analyze().

int HitEff::compSettings

private

Definition at line 87 of file HitEff.h.

Referenced by beginJob(), and HitEff().

edm::ParameterSet HitEff::conf_

private

Definition at line 82 of file HitEff.h.

Referenced by HitEff().

bool HitEff::cutOnTracks_

private

Definition at line 69 of file HitEff.h.

Referenced by analyze(), and HitEff().

bool HitEff::DEBUG

private

Definition at line 89 of file HitEff.h.

Referenced by HitEff().

const edm::EDGetTokenT<DetIdCollection> HitEff::digis_token_

private

Definition at line 79 of file HitEff.h.

Referenced by analyze().

unsigned int HitEff::event

private

Definition at line 111 of file HitEff.h.

Referenced by Types.EventID::cppID(), looper.Looper::process(), and core.AutoHandle.AutoHandle::product().

int HitEff::events

private

Definition at line 85 of file HitEff.h.

Referenced by eventsfwlite.Events::__getattr__(), analyze(), beginJob(), endJob(), looper.Looper::loop(), and looper.Looper::process().

int HitEff::EventTrackCKF

private

Definition at line 85 of file HitEff.h.

Referenced by analyze(), beginJob(), and endJob().

bool HitEff::highPurity

private

Definition at line 108 of file HitEff.h.

Referenced by analyze(), and beginJob().

unsigned int HitEff::Id

private

Definition at line 107 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::instLumi

private

Definition at line 113 of file HitEff.h.

Referenced by analyze(), and beginJob().

unsigned int HitEff::layers

private

Definition at line 88 of file HitEff.h.

Referenced by analyze(), and HitEff().

unsigned int HitEff::ModIsBad

private

Definition at line 107 of file HitEff.h.

Referenced by analyze(), and beginJob().

int HitEff::nHits

private

Definition at line 109 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::pT

private

Definition at line 110 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::PU

private

Definition at line 113 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::ResX

private

Definition at line 106 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::ResXSig

private

Definition at line 106 of file HitEff.h.

Referenced by analyze(), and beginJob().

unsigned int HitEff::run

private

Definition at line 111 of file HitEff.h.

Referenced by generateEDF.LumiInfo::__str__(), MatrixUtil.InputInfo::__str__(), DTWorkflow.DTWorkflow::all(), analyze(), beginJob(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), MatrixUtil.InputInfo::das(), personalPlayback.Playback::do_create_lumi(), personalPlayback.Playback::do_exec(), o2olib.O2OTool::execute(), generateEDF.LumiInfo::fixXingInfo(), MatrixUtil.InputInfo::lumiRanges(), MatrixUtil.InputInfo::queries(), and dqm_interfaces.DirWalkerFile::walk().

const edm::EDGetTokenT<LumiScalersCollection> HitEff::scalerToken_

private

Definition at line 64 of file HitEff.h.

Referenced by analyze().

unsigned int HitEff::SiStripQualBad

private

Definition at line 107 of file HitEff.h.

Referenced by analyze(), and beginJob().

int HitEff::tquality

private

Definition at line 112 of file HitEff.h.

Referenced by beginJob().

const edm::EDGetTokenT<MeasurementTrackerEvent> HitEff::trackerEvent_token_

private

Definition at line 80 of file HitEff.h.

Referenced by analyze().

unsigned int HitEff::trackMultiplicityCut_

private

Definition at line 70 of file HitEff.h.

Referenced by analyze(), and HitEff().

TTree* HitEff::traj

private

Definition at line 84 of file HitEff.h.

Referenced by analyze(), beginJob(), and endJob().

const edm::EDGetTokenT< std::vector<Trajectory> > HitEff::trajectories_token_

private

Definition at line 76 of file HitEff.h.

Referenced by analyze().

float HitEff::TrajGlbX

private

Definition at line 102 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajGlbY

private

Definition at line 102 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajGlbZ

private

Definition at line 102 of file HitEff.h.

Referenced by analyze(), and beginJob().

unsigned int HitEff::trajHitValid

private

Definition at line 111 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajLocAngleX

private

Definition at line 103 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajLocAngleY

private

Definition at line 103 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajLocErrX

private

Definition at line 104 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajLocErrY

private

Definition at line 104 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajLocX

private

Definition at line 103 of file HitEff.h.

Referenced by analyze(), and beginJob().

float HitEff::TrajLocY

private

Definition at line 103 of file HitEff.h.

Referenced by analyze(), and beginJob().

const edm::EDGetTokenT< TrajTrackAssociationCollection > HitEff::trajTrackAsso_token_

private

Definition at line 77 of file HitEff.h.

Referenced by analyze().

bool HitEff::useAllHitsFromTracksWithMissingHits_

private

Definition at line 73 of file HitEff.h.

Referenced by analyze(), and HitEff().

bool HitEff::useFirstMeas_

private

Definition at line 71 of file HitEff.h.

Referenced by analyze(), and HitEff().

bool HitEff::useLastMeas_

private

Definition at line 72 of file HitEff.h.

Referenced by analyze(), and HitEff().

unsigned int HitEff::whatlayer

private

Definition at line 90 of file HitEff.h.

Referenced by analyze(), and beginJob().

bool HitEff::withinAcceptance

private

Definition at line 107 of file HitEff.h.

Referenced by analyze(), and beginJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation