GlobalHaloAlgo Class Reference

#include <GlobalHaloAlgo.h>

Public Member Functions
reco::GlobalHaloData	Calculate (const CaloGeometry &TheCaloGeometry, const CSCGeometry &TheCSCGeometry, const reco::CaloMET &TheCaloMET, edm::Handle< edm::View< reco::Candidate > > &TheCaloTowers, edm::Handle< CSCSegmentCollection > &TheCSCSegments, edm::Handle< CSCRecHit2DCollection > &TheCSCRecHits, edm::Handle< reco::MuonCollection > &TheMuons, const reco::CSCHaloData &TheCSCHaloData, const reco::EcalHaloData &TheEcalHaloData, const reco::HcalHaloData &TheHcalHaloData, bool ishlt=false)
	GlobalHaloAlgo ()
void	SetCaloTowerEtThreshold (float EtMin)
void	setDPhicalosegmThresholdforCSCCaloMatchingEB (float x)
void	setDPhicalosegmThresholdforCSCCaloMatchingEE (float x)
void	setDPhicalosegmThresholdforCSCCaloMatchingHB (float x)
void	setDPhicalosegmThresholdforCSCCaloMatchingHE (float x)
void	setDtcalosegmThresholdforCSCCaloMatchingEB (float x)
void	setDtcalosegmThresholdforCSCCaloMatchingEE (float x)
void	setDtcalosegmThresholdforCSCCaloMatchingHB (float x)
void	setDtcalosegmThresholdforCSCCaloMatchingHE (float x)
void	SetEcalMatchingRadius (float min, float max)
void	setEtThresholdforCSCCaloMatchingEB (float x)
void	setEtThresholdforCSCCaloMatchingEE (float x)
void	setEtThresholdforCSCCaloMatchingHB (float x)
void	setEtThresholdforCSCCaloMatchingHE (float x)
void	SetHcalMatchingRadius (float min, float max)
void	SetMaxSegmentTheta (float x)
void	setRcaloMinRsegmHighThresholdforCSCCaloMatchingEB (float x)
void	setRcaloMinRsegmHighThresholdforCSCCaloMatchingEE (float x)
void	setRcaloMinRsegmHighThresholdforCSCCaloMatchingHB (float x)
void	setRcaloMinRsegmHighThresholdforCSCCaloMatchingHE (float x)
void	setRcaloMinRsegmLowThresholdforCSCCaloMatchingEB (float x)
void	setRcaloMinRsegmLowThresholdforCSCCaloMatchingEE (float x)
void	setRcaloMinRsegmLowThresholdforCSCCaloMatchingHB (float x)
void	setRcaloMinRsegmLowThresholdforCSCCaloMatchingHE (float x)
	~GlobalHaloAlgo ()

Private Member Functions
void	AddtoBeamHaloEBEERechits (edm::RefVector< EcalRecHitCollection > &bhtaggedrechits, reco::GlobalHaloData &thehalodata, bool isbarrel)
void	AddtoBeamHaloHBHERechits (edm::RefVector< HBHERecHitCollection > &bhtaggedrechits, reco::GlobalHaloData &thehalodata)
bool	ApplyMatchingCuts (int subdet, bool ishlt, double rhet, double segZ, double rhZ, double segR, double rhR, double segT, double rhT, double segPhi, double rhPhi)
bool	SegmentMatchingEB (reco::GlobalHaloData &thehalodata, const std::vector< reco::HaloClusterCandidateECAL > &haloclustercands, float iZ, float iR, float iT, float iPhi, bool ishlt)
bool	SegmentMatchingEE (reco::GlobalHaloData &thehalodata, const std::vector< reco::HaloClusterCandidateECAL > &haloclustercands, float iZ, float iR, float iT, float iPhi, bool ishlt)
bool	SegmentMatchingHB (reco::GlobalHaloData &thehalodata, const std::vector< reco::HaloClusterCandidateHCAL > &haloclustercands, float iZ, float iR, float iT, float iPhi, bool ishlt)
bool	SegmentMatchingHE (reco::GlobalHaloData &thehalodata, const std::vector< reco::HaloClusterCandidateHCAL > &haloclustercands, float iZ, float iR, float iT, float iPhi, bool ishlt)

Private Attributes
float	dphi_thresh_segvsrh_eb
float	dphi_thresh_segvsrh_ee
float	dphi_thresh_segvsrh_hb
float	dphi_thresh_segvsrh_he
float	dr_highthresh_segvsrh_eb
float	dr_highthresh_segvsrh_ee
float	dr_highthresh_segvsrh_hb
float	dr_highthresh_segvsrh_he
float	dr_lowthresh_segvsrh_eb
float	dr_lowthresh_segvsrh_ee
float	dr_lowthresh_segvsrh_hb
float	dr_lowthresh_segvsrh_he
float	dt_segvsrh_eb
float	dt_segvsrh_ee
float	dt_segvsrh_hb
float	dt_segvsrh_he
float	Ecal_R_Max
float	Ecal_R_Min
float	et_thresh_rh_eb
float	et_thresh_rh_ee
float	et_thresh_rh_hb
float	et_thresh_rh_he
float	Hcal_R_Max
float	Hcal_R_Min
float	max_segment_theta
float	TowerEtThreshold

Detailed Description

Definition at line 60 of file GlobalHaloAlgo.h.

Constructor & Destructor Documentation

GlobalHaloAlgo::GlobalHaloAlgo

(

)

Definition at line 34 of file GlobalHaloAlgo.cc.

{
  // Defaults are "loose"
  Ecal_R_Min = 110.;   // Tight: 200.
  Ecal_R_Max = 330.;   // Tight: 250. 
  Hcal_R_Min = 110.;   // Tight: 220.
  Hcal_R_Max = 490.;   // Tight: 350.
  
}

GlobalHaloAlgo::~GlobalHaloAlgo ( )

inline

Definition at line 65 of file GlobalHaloAlgo.h.

{}

Member Function Documentation

void GlobalHaloAlgo::AddtoBeamHaloEBEERechits ( edm::RefVector< EcalRecHitCollection > &  bhtaggedrechits, reco::GlobalHaloData &  thehalodata, bool  isbarrel  )

private

Definition at line 531 of file GlobalHaloAlgo.cc.

References reco::GlobalHaloData::GetEBRechits(), reco::GlobalHaloData::GetEERechits(), edm::RefVector< C, T, F >::push_back(), and edm::RefVector< C, T, F >::size().

                                                                                                                                                {
  for(size_t ihit = 0; ihit<bhtaggedrechits.size(); ++ ihit){
    bool alreadyincl = false; 
    edm::Ref<EcalRecHitCollection> rhRef( bhtaggedrechits[ihit] ) ;
    edm::RefVector<EcalRecHitCollection> refrhcoll;
    if(isbarrel) refrhcoll=thehalodata.GetEBRechits();
    else refrhcoll=thehalodata.GetEERechits();
    for(size_t jhit =0; jhit < refrhcoll.size();jhit++){                                                                                                                                                 
      edm::Ref<EcalRecHitCollection> rhitRef( refrhcoll[jhit] ) ;                                                                                                                                        
      if(rhitRef->detid() == rhRef->detid()) alreadyincl=true;                                                                                                                                           
      if(rhitRef->detid() == rhRef->detid()) break;                                                                                                                                                      
    } 
    if(!alreadyincl&&isbarrel)thehalodata.GetEBRechits().push_back(rhRef);
    if(!alreadyincl&&!isbarrel)thehalodata.GetEERechits().push_back(rhRef);
  }  
}

void GlobalHaloAlgo::AddtoBeamHaloHBHERechits ( edm::RefVector< HBHERecHitCollection > &  bhtaggedrechits, reco::GlobalHaloData &  thehalodata  )

private

Definition at line 549 of file GlobalHaloAlgo.cc.

References reco::GlobalHaloData::GetHBHERechits(), and edm::RefVector< C, T, F >::size().

                                                                                                                                 {  
  for(size_t ihit = 0; ihit<bhtaggedrechits.size(); ++ ihit){
    bool alreadyincl = false; 
    edm::Ref<HBHERecHitCollection> rhRef( bhtaggedrechits[ihit] ) ;
    edm::RefVector<HBHERecHitCollection> refrhcoll;
    refrhcoll=thehalodata.GetHBHERechits();
    for(size_t jhit =0; jhit < refrhcoll.size();jhit++){                                                                                                                                                 
      edm::Ref<HBHERecHitCollection> rhitRef( refrhcoll[jhit] ) ;                                                                                                                                        
      if(rhitRef->detid() == rhRef->detid()) alreadyincl=true;                                                                                                                                           
      if(rhitRef->detid() == rhRef->detid()) break;                                                                                                                                                      
    } 
    if(!alreadyincl)thehalodata.GetHBHERechits().push_back(rhRef);
  }
}

bool GlobalHaloAlgo::ApplyMatchingCuts ( int  subdet, bool  ishlt, double  rhet, double  segZ, double  rhZ, double  segR, double  rhR, double  segT, double  rhT, double  segPhi, double  rhPhi  )

private

Definition at line 461 of file GlobalHaloAlgo.cc.

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, EB, EE, HB, HE, and mathSSE::sqrt().

                                                                                                                                                                                 {
  //Std::Absolute time wrt BX
  double tBXrh = rhT+sqrt(rhR*rhR+rhZ*rhZ)/c_cm_per_ns;
  double tBXseg = segT+sqrt(segR*segR+segZ*segZ)/c_cm_per_ns;
  //Time at z=0, under beam halo hypothesis    
  double tcorseg = tBXseg - std::abs(segZ)/c_cm_per_ns;//Outgoing beam halo 
  double tcorsegincbh = tBXseg + std::abs(segZ)/c_cm_per_ns;//Ingoing beam halo
  double truedt[4]={1000,1000,1000,1000};
  //There are four types of segments associated to beam halo, test each hypothesis:
  //IT beam halo, ingoing track
  double twindow_seg = 15;
  if(std::abs(tcorsegincbh) <twindow_seg) truedt[0] =   tBXrh -tBXseg  -std::abs(rhZ-segZ)/c_cm_per_ns; 
  //IT beam halo, outgoing track
  if(std::abs(tcorseg) < twindow_seg) truedt[1] =  tBXseg -tBXrh -std::abs(rhZ-segZ)/c_cm_per_ns;  
  //OT beam halo (from next BX), ingoing track
  if(tcorsegincbh> 25-twindow_seg&& std::abs(tcorsegincbh) <25+twindow_seg) truedt[2] =   tBXrh -tBXseg  -std::abs(rhZ-segZ)/c_cm_per_ns; 
  //OT beam halo (from next BX), outgoing track
  if(tcorseg >25-twindow_seg && tcorseg<25+twindow_seg) truedt[3] =   tBXseg -tBXrh -std::abs(rhZ-segZ)/c_cm_per_ns; 
  
  
  if(subdet==EB){
    if(rhet< et_thresh_rh_eb) return false;
    if(rhet< 20&&ishlt) return false;
    if(std::abs(deltaPhi(rhPhi,segPhi))>dphi_thresh_segvsrh_eb) return false;
    if(rhR-segR< dr_lowthresh_segvsrh_eb)return false;
    if(rhR-segR> dr_highthresh_segvsrh_eb) return false;
    if(std::abs(truedt[0])>dt_segvsrh_eb &&std::abs(truedt[1])>dt_segvsrh_eb &&std::abs(truedt[2])>dt_segvsrh_eb &&std::abs(truedt[3])>dt_segvsrh_eb   )return false;
    return true;
  }


  if(subdet==EE){
    
    if(rhet< et_thresh_rh_ee) return false;
    if(rhet< 20&&ishlt) return false;
    if(std::abs(deltaPhi(rhPhi,segPhi))>dphi_thresh_segvsrh_ee) return false;
    if(rhR-segR< dr_lowthresh_segvsrh_ee)return false;
    if(rhR-segR> dr_highthresh_segvsrh_ee) return false;
    if(std::abs(truedt[0])>dt_segvsrh_ee &&std::abs(truedt[1])>dt_segvsrh_ee &&std::abs(truedt[2])>dt_segvsrh_ee &&std::abs(truedt[3])>dt_segvsrh_ee   )return false;
    return true;
  }

  if(subdet==HB){
    
    if(rhet< et_thresh_rh_hb) return false;
    if(rhet< 20&&ishlt) return false;
    if(std::abs(deltaPhi(rhPhi,segPhi))>dphi_thresh_segvsrh_hb) return false;
    if(rhR-segR< dr_lowthresh_segvsrh_hb)return false;
    if(rhR-segR> dr_highthresh_segvsrh_hb) return false;
    if(std::abs(truedt[0])>dt_segvsrh_hb &&std::abs(truedt[1])>dt_segvsrh_hb &&std::abs(truedt[2])>dt_segvsrh_hb &&std::abs(truedt[3])>dt_segvsrh_hb   )return false;
    return true;
  }

  if(subdet==HE){
    
    if(rhet< et_thresh_rh_he) return false;
    if(rhet< 20&&ishlt) return false;
    if(std::abs(deltaPhi(rhPhi,segPhi))>dphi_thresh_segvsrh_he) return false;
    if(rhR-segR< dr_lowthresh_segvsrh_he)return false;
    if(rhR-segR> dr_highthresh_segvsrh_he) return false;
    if(std::abs(truedt[0])>dt_segvsrh_he &&std::abs(truedt[1])>dt_segvsrh_he &&std::abs(truedt[2])>dt_segvsrh_he &&std::abs(truedt[3])>dt_segvsrh_he   )return false;
    return true;
  }


  return false;
}

reco::GlobalHaloData GlobalHaloAlgo::Calculate	(	const CaloGeometry &	TheCaloGeometry,
		const CSCGeometry &	TheCSCGeometry,
		const reco::CaloMET &	TheCaloMET,
		edm::Handle< edm::View< reco::Candidate > > &	TheCaloTowers,
		edm::Handle< CSCSegmentCollection > &	TheCSCSegments,
		edm::Handle< CSCRecHit2DCollection > &	TheCSCRecHits,
		edm::Handle< reco::MuonCollection > &	TheMuons,
		const reco::CSCHaloData &	TheCSCHaloData,
		const reco::EcalHaloData &	TheEcalHaloData,
		const reco::HcalHaloData &	TheHcalHaloData,
		bool	ishlt = false
	)

Definition at line 44 of file GlobalHaloAlgo.cc.

References funct::abs(), Abs(), EnergyCorrector::c, CSCGeometry::chamber(), chambers, MuonSubdetId::CSC, CaloTower::emEt(), CaloTower::et(), reco::CSCHaloData::GetCSCTrackImpactPositions(), reco::EcalHaloData::getHaloClusterCandidatesEB(), reco::EcalHaloData::getHaloClusterCandidatesEE(), reco::HcalHaloData::getHaloClusterCandidatesHB(), reco::HcalHaloData::getHaloClusterCandidatesHE(), reco::GlobalHaloData::GetMatchedEcalPhiWedges(), reco::GlobalHaloData::GetMatchedHcalPhiWedges(), reco::EcalHaloData::GetPhiWedges(), reco::HcalHaloData::GetPhiWedges(), CaloTower::hadEt(), i, CSCGeometry::idToDetUnit(), CaloTower::ieta(), CaloTower::iphi(), edm::HandleBase::isValid(), RPCpg::mu, DetId::Muon, PV3DBase< T, PVType, FrameType >::perp2(), PV3DBase< T, PVType, FrameType >::phi(), reco::LeafCandidate::phi(), Phi_To_EcaliPhi(), Phi_To_HcaliPhi(), Pi, reco::LeafCandidate::pt(), alignCSCRings::r, reco::GlobalHaloData::SetHaloPatternFoundEB(), reco::GlobalHaloData::SetHaloPatternFoundEE(), reco::GlobalHaloData::SetHaloPatternFoundHB(), reco::GlobalHaloData::SetHaloPatternFoundHE(), reco::GlobalHaloData::SetMETCorrections(), reco::GlobalHaloData::SetMETOverSumEt(), reco::PhiWedge::SetOverlappingCSCRecHits(), reco::PhiWedge::SetOverlappingCSCSegments(), reco::GlobalHaloData::SetSegmentIsEBCaloMatched(), reco::GlobalHaloData::SetSegmentIsEECaloMatched(), reco::GlobalHaloData::SetSegmentIsHBCaloMatched(), reco::GlobalHaloData::SetSegmentIsHECaloMatched(), mathSSE::sqrt(), reco::MET::sumEt(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::theta(), GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{

  GlobalHaloData TheGlobalHaloData;
  float METOverSumEt = TheCaloMET.sumEt() ? TheCaloMET.pt() / TheCaloMET.sumEt() : 0 ;
  TheGlobalHaloData.SetMETOverSumEt(METOverSumEt);

  //int EcalOverlapping_CSCRecHits[73];
  //int EcalOverlapping_CSCSegments[73];

  int EcalOverlapping_CSCRecHits[361];
  int EcalOverlapping_CSCSegments[361];
  int HcalOverlapping_CSCRecHits[73];
  int HcalOverlapping_CSCSegments[73];

  if( TheCSCSegments.isValid() )
    {
      for(CSCSegmentCollection::const_iterator iSegment = TheCSCSegments->begin(); iSegment != TheCSCSegments->end(); iSegment++) 
    {
      bool EcalOverlap[361];
      bool HcalOverlap[73];
      for( int i = 0 ; i < 361 ; i++ ) 
        {
          EcalOverlap[i] = false;
          if( i < 73 ) HcalOverlap[i] = false;
        }
      
      std::vector<CSCRecHit2D> Hits = iSegment->specificRecHits() ;
      for(std::vector<CSCRecHit2D>::iterator iHit = Hits.begin() ; iHit != Hits.end(); iHit++ )
        {
          DetId TheDetUnitId(iHit->geographicalId());
          if( TheDetUnitId.det() != DetId::Muon ) continue;
          if( TheDetUnitId.subdetId() != MuonSubdetId::CSC ) continue;

          const GeomDetUnit *TheUnit = TheCSCGeometry.idToDetUnit(TheDetUnitId);
          LocalPoint TheLocalPosition = iHit->localPosition();  
          const BoundPlane& TheSurface = TheUnit->surface();
          const GlobalPoint TheGlobalPosition = TheSurface.toGlobal(TheLocalPosition);

          int Hcal_iphi = Phi_To_HcaliPhi( TheGlobalPosition.phi() ) ;
          int Ecal_iphi = Phi_To_EcaliPhi( TheGlobalPosition.phi() ) ;
          float x = TheGlobalPosition.x(); 
          float y = TheGlobalPosition.y();
          
          float r = TMath::Sqrt( x*x + y*y);
          
          if( r < Ecal_R_Max && r > Ecal_R_Min )
        EcalOverlap[Ecal_iphi] = true;
          if( r < Hcal_R_Max && r > Hcal_R_Max ) 
        HcalOverlap[Hcal_iphi] = true;
        }
      for( int i = 0 ; i < 361 ; i++ ) 
        {
          if( EcalOverlap[i] )  EcalOverlapping_CSCSegments[i]++;
          if( i < 73 && HcalOverlap[i] )
        HcalOverlapping_CSCSegments[i]++;
        }
    } 
    }
  if( TheCSCRecHits.isValid() )
    {
      for(CSCRecHit2DCollection::const_iterator iCSCRecHit = TheCSCRecHits->begin();   iCSCRecHit != TheCSCRecHits->end(); iCSCRecHit++ )
    {
      
      DetId TheDetUnitId(iCSCRecHit->geographicalId());
      if( TheDetUnitId.det() != DetId::Muon ) continue;
      if( TheDetUnitId.subdetId() != MuonSubdetId::CSC ) continue;
      
      const GeomDetUnit *TheUnit = TheCSCGeometry.idToDetUnit(TheDetUnitId);
      LocalPoint TheLocalPosition = iCSCRecHit->localPosition();  
      const BoundPlane& TheSurface = TheUnit->surface();
      const GlobalPoint TheGlobalPosition = TheSurface.toGlobal(TheLocalPosition);
      
      int Hcaliphi = Phi_To_HcaliPhi( TheGlobalPosition.phi() ) ;
      int Ecaliphi = Phi_To_EcaliPhi( TheGlobalPosition.phi() ) ;
      float x = TheGlobalPosition.x(); 
      float y = TheGlobalPosition.y();
      
      float r = TMath::Sqrt(x*x + y*y);
      
      if( r < Ecal_R_Max && r > Ecal_R_Min )
        EcalOverlapping_CSCRecHits[Ecaliphi] ++;
      if( r < Hcal_R_Max && r > Hcal_R_Max ) 
        HcalOverlapping_CSCRecHits[Hcaliphi] ++ ;
    }
    }  

  // In development....
  // Get Ecal Wedges
  std::vector<PhiWedge> EcalWedges = TheEcalHaloData.GetPhiWedges();
  
  // Get Hcal Wedges
  std::vector<PhiWedge> HcalWedges = TheHcalHaloData.GetPhiWedges();
  
  //Get Ref to CSC Tracks
  //edm::RefVector<reco::TrackCollection> TheCSCTracks = TheCSCHaloData.GetTracks();
  //for(unsigned int i = 0 ; i < TheCSCTracks.size() ; i++ )
  //edm::Ref<reco::TrackCollection> iTrack( TheCSCTracks, i );
  
  // Get global positions of central most rechit of CSC Halo tracks
  std::vector<GlobalPoint> TheGlobalPositions = TheCSCHaloData.GetCSCTrackImpactPositions();

  // Container to store Ecal/Hcal iPhi values matched to impact point of CSC tracks
  std::vector<int> vEcaliPhi, vHcaliPhi;

  // Keep track of number of calo pointing CSC halo tracks that do not match to Phi wedges
  int N_Unmatched_Tracks = 0;  
  
  for( std::vector<GlobalPoint>::iterator Pos = TheGlobalPositions.begin() ; Pos != TheGlobalPositions.end() ; Pos ++ ) 
    {
      // Calculate global phi coordinate for central most rechit in the track
      float global_phi = Pos->phi();
      float global_r = TMath::Sqrt(Pos->x()*Pos->x() + Pos->y()*Pos->y());
      
      // Convert global phi to iPhi
      int global_EcaliPhi = Phi_To_EcaliPhi( global_phi );
      int global_HcaliPhi = Phi_To_HcaliPhi( global_phi );
      bool MATCHED = false;
      
      //Loop over Ecal Phi Wedges 
      for( std::vector<PhiWedge>::iterator iWedge = EcalWedges.begin() ; iWedge != EcalWedges.end() ; iWedge++ )
    {
      if( (TMath::Abs( global_EcaliPhi - iWedge->iPhi() ) <= 5 ) && (global_r >  Ecal_R_Min && global_r < Ecal_R_Max ) )
        {
          bool StoreWedge = true;
          for( unsigned int i = 0 ; i< vEcaliPhi.size() ; i++ ) if ( vEcaliPhi[i] == iWedge->iPhi() ) StoreWedge = false;
          
          if( StoreWedge ) 
        {
          PhiWedge NewWedge(*iWedge);
          NewWedge.SetOverlappingCSCSegments( EcalOverlapping_CSCSegments[iWedge->iPhi()] );
          NewWedge.SetOverlappingCSCRecHits( EcalOverlapping_CSCRecHits[iWedge->iPhi()] );
          vEcaliPhi.push_back( iWedge->iPhi() );
          TheGlobalHaloData.GetMatchedEcalPhiWedges().push_back( NewWedge );
        }
          MATCHED = true;
        }
    }
      //Loop over Hcal Phi Wedges 
      for( std::vector<PhiWedge>::iterator iWedge = HcalWedges.begin() ; iWedge != HcalWedges.end() ; iWedge++ )
    {
      if(  (TMath::Abs( global_HcaliPhi - iWedge->iPhi() ) <=  2 ) && (global_r > Hcal_R_Min && global_r < Hcal_R_Max ) )
        {
          bool StoreWedge  = true;
          for( unsigned int i = 0 ; i < vHcaliPhi.size() ; i++ ) if(  vHcaliPhi[i] == iWedge->iPhi() ) StoreWedge = false;
          
          if( StoreWedge ) 
        {
          vHcaliPhi.push_back( iWedge->iPhi() ) ;
          PhiWedge NewWedge(*iWedge);
          NewWedge.SetOverlappingCSCSegments( HcalOverlapping_CSCSegments[iWedge->iPhi()] );
          NewWedge.SetOverlappingCSCRecHits(  HcalOverlapping_CSCRecHits[iWedge->iPhi()] );       
          PhiWedge wedge(*iWedge);
          TheGlobalHaloData.GetMatchedHcalPhiWedges().push_back( NewWedge ) ; 
        }
          MATCHED = true;
        }
    }
      if( !MATCHED ) N_Unmatched_Tracks ++;
    }
  
  // Corrections to MEx, MEy
  float dMEx = 0.; 
  float dMEy = 0.;
  // Loop over calotowers and correct the MET for the towers that lie in the trajectory of the CSC Halo Tracks
  for( edm::View<Candidate>::const_iterator iCandidate = TheCaloTowers->begin() ; iCandidate != TheCaloTowers->end() ; iCandidate++ )
    {
      const Candidate* c = &(*iCandidate);
      if ( c )
    {
      const CaloTower* iTower = dynamic_cast<const CaloTower*> (c);
      if( iTower->et() < TowerEtThreshold ) continue;
      if( abs(iTower->ieta()) > 24 )  continue;   // not in barrel/endcap
      int iphi = iTower->iphi();
      for( unsigned int x = 0 ; x < vEcaliPhi.size() ; x++ )
        {
          if( iphi == vEcaliPhi[x] ) 
        {
          dMEx += ( TMath::Cos(iTower->phi())*iTower->emEt() );
          dMEy += ( TMath::Sin(iTower->phi())*iTower->emEt() );
        }
        }
      for( unsigned int x = 0 ; x < vHcaliPhi.size() ; x++ )
        {
          if( iphi == vHcaliPhi[x] ) 
        {
          dMEx += ( TMath::Cos(iTower->phi() )*iTower->hadEt() ) ;
          dMEy += ( TMath::Sin(iTower->phi() )*iTower->hadEt() ) ;
        }
        }
    }
    }
  
  TheGlobalHaloData.SetMETCorrections(dMEx, dMEy);

  

  std::vector<HaloClusterCandidateECAL> hccandEB= TheEcalHaloData.getHaloClusterCandidatesEB();
  std::vector<HaloClusterCandidateECAL> hccandEE= TheEcalHaloData.getHaloClusterCandidatesEE();
  std::vector<HaloClusterCandidateHCAL> hccandHB= TheHcalHaloData.getHaloClusterCandidatesHB();
  std::vector<HaloClusterCandidateHCAL> hccandHE= TheHcalHaloData.getHaloClusterCandidatesHE();

  //CSC-calo matching
  bool ECALBmatched(false), ECALEmatched(false),HCALBmatched(false),HCALEmatched(false);

  if (TheCSCSegments.isValid()) {
    for(CSCSegmentCollection::const_iterator iSegment = TheCSCSegments->begin();
    iSegment != TheCSCSegments->end();
    iSegment++) {
      
      CSCDetId iCscDetID = iSegment->cscDetId();
      bool Segment1IsGood=true;
      
      
      //avoid segments from collision muons
      if( TheMuons.isValid() )
    {
      for(reco::MuonCollection::const_iterator mu = TheMuons->begin(); mu!= TheMuons->end() && (Segment1IsGood)   ; mu++ )
        {

          if( !mu->isTrackerMuon() && !mu->isGlobalMuon() && mu->isStandAloneMuon() ) continue;
          if( !mu->isGlobalMuon() &&  mu->isTrackerMuon() &&  mu->pt()<3) continue; 
          const std::vector<MuonChamberMatch> chambers = mu->matches();
          for(std::vector<MuonChamberMatch>::const_iterator kChamber = chambers.begin();
          kChamber != chambers.end(); kChamber ++ )
        {
          if( kChamber->detector() != MuonSubdetId::CSC ) continue;
          for( std::vector<reco::MuonSegmentMatch>::const_iterator kSegment = kChamber->segmentMatches.begin();
               kSegment != kChamber->segmentMatches.end(); kSegment++ )
            {
              edm::Ref<CSCSegmentCollection> cscSegRef = kSegment->cscSegmentRef;
              CSCDetId kCscDetID = cscSegRef->cscDetId();
              
              if( kCscDetID == iCscDetID ) 
            {
              Segment1IsGood = false;

            }
            }
        }
        }
    }
      if(!Segment1IsGood) continue;
      
      // Get local direction vector; if direction runs parallel to beamline,
      // count this segment as beam halo candidate.
      LocalPoint iLocalPosition = iSegment->localPosition();
      LocalVector iLocalDirection = iSegment->localDirection();

      GlobalPoint iGlobalPosition = TheCSCGeometry.chamber(iCscDetID)->toGlobal(iLocalPosition);
      GlobalVector iGlobalDirection = TheCSCGeometry.chamber(iCscDetID)->toGlobal(iLocalDirection);
      
      float iTheta = iGlobalDirection.theta();
      if (iTheta > max_segment_theta && iTheta < TMath::Pi() - max_segment_theta) continue;
      
      float iPhi = iGlobalPosition.phi();
      float iR =  sqrt(iGlobalPosition.perp2()) ; 
      float iZ = iGlobalPosition.z();
      float iT = iSegment->time();
      
      //CSC-calo matching: 
      //Here, one checks if any halo cluster can be matched to a CSC segment. 
      //The matching uses both geometric (dphi, dR) and timing information (dt).
      //The cut values depend on the subdetector considered (e.g. in HB, Rcalo-Rsegment is allowed to be very negative) 
      
      bool ebmatched =SegmentMatchingEB(TheGlobalHaloData,hccandEB,iZ,iR,iT,iPhi,ishlt);
      bool eematched =SegmentMatchingEE(TheGlobalHaloData,hccandEE,iZ,iR,iT,iPhi,ishlt);
      bool hbmatched =SegmentMatchingHB(TheGlobalHaloData,hccandHB,iZ,iR,iT,iPhi,ishlt);
      bool hematched =SegmentMatchingHE(TheGlobalHaloData,hccandHE,iZ,iR,iT,iPhi,ishlt);
      
      ECALBmatched |= ebmatched;
      ECALEmatched |= eematched;
      HCALBmatched |= hbmatched;
      HCALEmatched |= hematched;
      
    }
  }


  TheGlobalHaloData.SetSegmentIsEBCaloMatched(ECALBmatched);
  TheGlobalHaloData.SetSegmentIsEECaloMatched(ECALEmatched);
  TheGlobalHaloData.SetSegmentIsHBCaloMatched(HCALBmatched);
  TheGlobalHaloData.SetSegmentIsHECaloMatched(HCALEmatched);
  

  //Now checking patterns from EcalHaloData and HcalHaloData: 
  //Simply check whether any cluster has a halo pattern
  //In that case store the rhits in GlobalHaloData 

  bool HaloPatternFoundInEB =false;
  for(auto & hcand : hccandEB){
    if((hcand.getIsHaloFromPattern() &&!ishlt)||(hcand.getIsHaloFromPattern_HLT() &&ishlt)){
      HaloPatternFoundInEB =true;
      edm::RefVector<EcalRecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
      AddtoBeamHaloEBEERechits(bhrhcandidates, TheGlobalHaloData,true);
    } 
  }
  

  bool HaloPatternFoundInEE =false;
  for(auto & hcand : hccandEE){
    if((hcand.getIsHaloFromPattern() &&!ishlt)||(hcand.getIsHaloFromPattern_HLT() &&ishlt)){
      HaloPatternFoundInEE =true;
      edm::RefVector<EcalRecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
      AddtoBeamHaloEBEERechits( bhrhcandidates, TheGlobalHaloData,false);
    } 
  }


  bool HaloPatternFoundInHB =false;
  for(auto & hcand : hccandHB){
    if((hcand.getIsHaloFromPattern() &&!ishlt)||(hcand.getIsHaloFromPattern_HLT() &&ishlt)){
      HaloPatternFoundInHB =true;
      edm::RefVector<HBHERecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
      AddtoBeamHaloHBHERechits(bhrhcandidates, TheGlobalHaloData);
    } 
  }
  

  bool HaloPatternFoundInHE =false;
  for(auto & hcand : hccandHE){
    if((hcand.getIsHaloFromPattern() &&!ishlt)||(hcand.getIsHaloFromPattern_HLT() &&ishlt)){
      HaloPatternFoundInHE =true;
      edm::RefVector<HBHERecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
      AddtoBeamHaloHBHERechits(bhrhcandidates, TheGlobalHaloData);
    } 
  }
  TheGlobalHaloData.SetHaloPatternFoundEB(HaloPatternFoundInEB);
  TheGlobalHaloData.SetHaloPatternFoundEE(HaloPatternFoundInEE);
  TheGlobalHaloData.SetHaloPatternFoundHB(HaloPatternFoundInHB);
  TheGlobalHaloData.SetHaloPatternFoundHE(HaloPatternFoundInHE);
 
  return TheGlobalHaloData;
}

bool GlobalHaloAlgo::SegmentMatchingEB ( reco::GlobalHaloData &  thehalodata, const std::vector< reco::HaloClusterCandidateECAL > &  haloclustercands, float  iZ, float  iR, float  iT, float  iPhi, bool  ishlt  )

private

Definition at line 384 of file GlobalHaloAlgo.cc.

References EB.

                                                                                                                                                                                          {
  bool rhmatchingfound =false;

  for(auto & hcand : haloclustercands){
    
    if(!ApplyMatchingCuts(EB,ishlt, hcand.getSeedEt(), iZ, hcand.getSeedZ(),iR, hcand.getSeedR(), iT,hcand.getSeedTime(), iPhi, hcand.getSeedPhi()))continue;
    
    rhmatchingfound=true;
    
    edm::RefVector<EcalRecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
    
    AddtoBeamHaloEBEERechits(bhrhcandidates, thehalodata,true);
  
  }

  return rhmatchingfound;
}

bool GlobalHaloAlgo::SegmentMatchingEE ( reco::GlobalHaloData &  thehalodata, const std::vector< reco::HaloClusterCandidateECAL > &  haloclustercands, float  iZ, float  iR, float  iT, float  iPhi, bool  ishlt  )

private

Definition at line 403 of file GlobalHaloAlgo.cc.

References EE.

                                                                                                                                                                                          {
  bool rhmatchingfound =false;
  
  for(auto & hcand : haloclustercands){
   
    if(!ApplyMatchingCuts(EE,ishlt, hcand.getSeedEt(), iZ, hcand.getSeedZ(),iR, hcand.getSeedR(), iT,hcand.getSeedTime(), iPhi, hcand.getSeedPhi()))continue;
  
    rhmatchingfound=true;
    
    edm::RefVector<EcalRecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
    
    AddtoBeamHaloEBEERechits(bhrhcandidates, thehalodata,false);
  
  }
  
  return rhmatchingfound;

}

bool GlobalHaloAlgo::SegmentMatchingHB ( reco::GlobalHaloData &  thehalodata, const std::vector< reco::HaloClusterCandidateHCAL > &  haloclustercands, float  iZ, float  iR, float  iT, float  iPhi, bool  ishlt  )

private

Definition at line 422 of file GlobalHaloAlgo.cc.

References HB.

                                                                                                                                                                                          {
  bool rhmatchingfound =false;

  for(auto & hcand : haloclustercands){
   
    if(!ApplyMatchingCuts(HB,ishlt, hcand.getSeedEt(), iZ, hcand.getSeedZ(),iR, hcand.getSeedR(), iT,hcand.getSeedTime(), iPhi, hcand.getSeedPhi()))continue;
  
    rhmatchingfound=true;
    
    edm::RefVector<HBHERecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
    
    AddtoBeamHaloHBHERechits(bhrhcandidates, thehalodata);
  
  }
  
  return rhmatchingfound;

}

bool GlobalHaloAlgo::SegmentMatchingHE ( reco::GlobalHaloData &  thehalodata, const std::vector< reco::HaloClusterCandidateHCAL > &  haloclustercands, float  iZ, float  iR, float  iT, float  iPhi, bool  ishlt  )

private

Definition at line 441 of file GlobalHaloAlgo.cc.

References HE.

                                                                                                                                                                                          {
  bool rhmatchingfound =false;

  for(auto & hcand : haloclustercands){
   
    if(!ApplyMatchingCuts(HE,ishlt, hcand.getSeedEt(),  iZ, hcand.getSeedZ(),iR, hcand.getSeedR(), iT,hcand.getSeedTime(), iPhi, hcand.getSeedPhi()))continue;
  
    rhmatchingfound=true;
    
    edm::RefVector<HBHERecHitCollection> bhrhcandidates = hcand.getBeamHaloRecHitsCandidates();
    
    AddtoBeamHaloHBHERechits(bhrhcandidates, thehalodata);
  
  }
 
  return rhmatchingfound;

}

void GlobalHaloAlgo::SetCaloTowerEtThreshold ( float  EtMin )

inline

Definition at line 75 of file GlobalHaloAlgo.h.

References HLT_25ns10e33_v2_cff::EtMin, and TowerEtThreshold.

{ TowerEtThreshold = EtMin ;}

void GlobalHaloAlgo::setDPhicalosegmThresholdforCSCCaloMatchingEB ( float  x )

inline

Definition at line 85 of file GlobalHaloAlgo.h.

References dphi_thresh_segvsrh_eb, and x.

{dphi_thresh_segvsrh_eb= x;}

void GlobalHaloAlgo::setDPhicalosegmThresholdforCSCCaloMatchingEE ( float  x )

inline

Definition at line 91 of file GlobalHaloAlgo.h.

References dphi_thresh_segvsrh_ee, and x.

{dphi_thresh_segvsrh_ee= x;}

void GlobalHaloAlgo::setDPhicalosegmThresholdforCSCCaloMatchingHB ( float  x )

inline

Definition at line 97 of file GlobalHaloAlgo.h.

References dphi_thresh_segvsrh_hb, and x.

{dphi_thresh_segvsrh_hb= x;}

void GlobalHaloAlgo::setDPhicalosegmThresholdforCSCCaloMatchingHE ( float  x )

inline

Definition at line 103 of file GlobalHaloAlgo.h.

References dphi_thresh_segvsrh_he, and x.

{dphi_thresh_segvsrh_he= x;}

void GlobalHaloAlgo::setDtcalosegmThresholdforCSCCaloMatchingEB ( float  x )

inline

Definition at line 84 of file GlobalHaloAlgo.h.

References dt_segvsrh_eb, and x.

{dt_segvsrh_eb= x;}

void GlobalHaloAlgo::setDtcalosegmThresholdforCSCCaloMatchingEE ( float  x )

inline

Definition at line 90 of file GlobalHaloAlgo.h.

References dt_segvsrh_ee, and x.

{dt_segvsrh_ee= x;}

void GlobalHaloAlgo::setDtcalosegmThresholdforCSCCaloMatchingHB ( float  x )

inline

Definition at line 96 of file GlobalHaloAlgo.h.

References dt_segvsrh_hb, and x.

{dt_segvsrh_hb= x;}

void GlobalHaloAlgo::setDtcalosegmThresholdforCSCCaloMatchingHE ( float  x )

inline

Definition at line 102 of file GlobalHaloAlgo.h.

References dt_segvsrh_he, and x.

{dt_segvsrh_he= x;}

void GlobalHaloAlgo::SetEcalMatchingRadius ( float  min, float  max  )

inline

Definition at line 71 of file GlobalHaloAlgo.h.

References Ecal_R_Max, Ecal_R_Min, bookConverter::max, and min().

{Ecal_R_Min = min ; Ecal_R_Max = max;}

void GlobalHaloAlgo::setEtThresholdforCSCCaloMatchingEB ( float  x )

inline

Definition at line 81 of file GlobalHaloAlgo.h.

References et_thresh_rh_eb, and x.

{et_thresh_rh_eb= x;}

void GlobalHaloAlgo::setEtThresholdforCSCCaloMatchingEE ( float  x )

inline

Definition at line 87 of file GlobalHaloAlgo.h.

References et_thresh_rh_ee, and x.

{et_thresh_rh_ee= x;}

void GlobalHaloAlgo::setEtThresholdforCSCCaloMatchingHB ( float  x )

inline

Definition at line 93 of file GlobalHaloAlgo.h.

References et_thresh_rh_hb, and x.

{et_thresh_rh_hb= x;}

void GlobalHaloAlgo::setEtThresholdforCSCCaloMatchingHE ( float  x )

inline

Definition at line 99 of file GlobalHaloAlgo.h.

References et_thresh_rh_he, and x.

{et_thresh_rh_he= x;}

void GlobalHaloAlgo::SetHcalMatchingRadius ( float  min, float  max  )

inline

Definition at line 73 of file GlobalHaloAlgo.h.

References Hcal_R_Max, Hcal_R_Min, bookConverter::max, and min().

{Hcal_R_Min = min ; Hcal_R_Max = max;}

void GlobalHaloAlgo::SetMaxSegmentTheta ( float  x )

inline

Definition at line 79 of file GlobalHaloAlgo.h.

References max_segment_theta, and x.

{ max_segment_theta = x; }

void GlobalHaloAlgo::setRcaloMinRsegmHighThresholdforCSCCaloMatchingEB ( float  x )

inline

Definition at line 83 of file GlobalHaloAlgo.h.

References dr_highthresh_segvsrh_eb, and x.

{dr_highthresh_segvsrh_eb= x;}

void GlobalHaloAlgo::setRcaloMinRsegmHighThresholdforCSCCaloMatchingEE ( float  x )

inline

Definition at line 89 of file GlobalHaloAlgo.h.

References dr_highthresh_segvsrh_ee, and x.

{dr_highthresh_segvsrh_ee= x;}

void GlobalHaloAlgo::setRcaloMinRsegmHighThresholdforCSCCaloMatchingHB ( float  x )

inline

Definition at line 95 of file GlobalHaloAlgo.h.

References dr_highthresh_segvsrh_hb, and x.

{dr_highthresh_segvsrh_hb= x;}

void GlobalHaloAlgo::setRcaloMinRsegmHighThresholdforCSCCaloMatchingHE ( float  x )

inline

Definition at line 101 of file GlobalHaloAlgo.h.

References dr_highthresh_segvsrh_he, and x.

{dr_highthresh_segvsrh_he= x;}

void GlobalHaloAlgo::setRcaloMinRsegmLowThresholdforCSCCaloMatchingEB ( float  x )

inline

Definition at line 82 of file GlobalHaloAlgo.h.

References dr_lowthresh_segvsrh_eb, and x.

{dr_lowthresh_segvsrh_eb= x;}

void GlobalHaloAlgo::setRcaloMinRsegmLowThresholdforCSCCaloMatchingEE ( float  x )

inline

Definition at line 88 of file GlobalHaloAlgo.h.

References dr_lowthresh_segvsrh_ee, and x.

{dr_lowthresh_segvsrh_ee= x;}

void GlobalHaloAlgo::setRcaloMinRsegmLowThresholdforCSCCaloMatchingHB ( float  x )

inline

Definition at line 94 of file GlobalHaloAlgo.h.

References dr_lowthresh_segvsrh_hb, and x.

{dr_lowthresh_segvsrh_hb= x;}

void GlobalHaloAlgo::setRcaloMinRsegmLowThresholdforCSCCaloMatchingHE ( float  x )

inline

Definition at line 100 of file GlobalHaloAlgo.h.

References dr_lowthresh_segvsrh_he, and x.

{dr_lowthresh_segvsrh_he= x;}

Member Data Documentation

float GlobalHaloAlgo::dphi_thresh_segvsrh_eb

private

Definition at line 118 of file GlobalHaloAlgo.h.

Referenced by setDPhicalosegmThresholdforCSCCaloMatchingEB().

float GlobalHaloAlgo::dphi_thresh_segvsrh_ee

private

Definition at line 124 of file GlobalHaloAlgo.h.

Referenced by setDPhicalosegmThresholdforCSCCaloMatchingEE().

float GlobalHaloAlgo::dphi_thresh_segvsrh_hb

private

Definition at line 130 of file GlobalHaloAlgo.h.

Referenced by setDPhicalosegmThresholdforCSCCaloMatchingHB().

float GlobalHaloAlgo::dphi_thresh_segvsrh_he

private

Definition at line 136 of file GlobalHaloAlgo.h.

Referenced by setDPhicalosegmThresholdforCSCCaloMatchingHE().

float GlobalHaloAlgo::dr_highthresh_segvsrh_eb

private

Definition at line 120 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmHighThresholdforCSCCaloMatchingEB().

float GlobalHaloAlgo::dr_highthresh_segvsrh_ee

private

Definition at line 126 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmHighThresholdforCSCCaloMatchingEE().

float GlobalHaloAlgo::dr_highthresh_segvsrh_hb

private

Definition at line 132 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmHighThresholdforCSCCaloMatchingHB().

float GlobalHaloAlgo::dr_highthresh_segvsrh_he

private

Definition at line 138 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmHighThresholdforCSCCaloMatchingHE().

float GlobalHaloAlgo::dr_lowthresh_segvsrh_eb

private

Definition at line 119 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmLowThresholdforCSCCaloMatchingEB().

float GlobalHaloAlgo::dr_lowthresh_segvsrh_ee

private

Definition at line 125 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmLowThresholdforCSCCaloMatchingEE().

float GlobalHaloAlgo::dr_lowthresh_segvsrh_hb

private

Definition at line 131 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmLowThresholdforCSCCaloMatchingHB().

float GlobalHaloAlgo::dr_lowthresh_segvsrh_he

private

Definition at line 137 of file GlobalHaloAlgo.h.

Referenced by setRcaloMinRsegmLowThresholdforCSCCaloMatchingHE().

float GlobalHaloAlgo::dt_segvsrh_eb

private

Definition at line 121 of file GlobalHaloAlgo.h.

Referenced by setDtcalosegmThresholdforCSCCaloMatchingEB().

float GlobalHaloAlgo::dt_segvsrh_ee

private

Definition at line 127 of file GlobalHaloAlgo.h.

Referenced by setDtcalosegmThresholdforCSCCaloMatchingEE().

float GlobalHaloAlgo::dt_segvsrh_hb

private

Definition at line 133 of file GlobalHaloAlgo.h.

Referenced by setDtcalosegmThresholdforCSCCaloMatchingHB().

float GlobalHaloAlgo::dt_segvsrh_he

private

Definition at line 139 of file GlobalHaloAlgo.h.

Referenced by setDtcalosegmThresholdforCSCCaloMatchingHE().

float GlobalHaloAlgo::Ecal_R_Max

private

Definition at line 109 of file GlobalHaloAlgo.h.

Referenced by SetEcalMatchingRadius().

float GlobalHaloAlgo::Ecal_R_Min

private

Definition at line 108 of file GlobalHaloAlgo.h.

Referenced by SetEcalMatchingRadius().

float GlobalHaloAlgo::et_thresh_rh_eb

private

Definition at line 117 of file GlobalHaloAlgo.h.

Referenced by setEtThresholdforCSCCaloMatchingEB().

float GlobalHaloAlgo::et_thresh_rh_ee

private

Definition at line 123 of file GlobalHaloAlgo.h.

Referenced by setEtThresholdforCSCCaloMatchingEE().

float GlobalHaloAlgo::et_thresh_rh_hb

private

Definition at line 129 of file GlobalHaloAlgo.h.

Referenced by setEtThresholdforCSCCaloMatchingHB().

float GlobalHaloAlgo::et_thresh_rh_he

private

Definition at line 135 of file GlobalHaloAlgo.h.

Referenced by setEtThresholdforCSCCaloMatchingHE().

float GlobalHaloAlgo::Hcal_R_Max

private

Definition at line 111 of file GlobalHaloAlgo.h.

Referenced by SetHcalMatchingRadius().

float GlobalHaloAlgo::Hcal_R_Min

private

Definition at line 110 of file GlobalHaloAlgo.h.

Referenced by SetHcalMatchingRadius().

float GlobalHaloAlgo::max_segment_theta

private

Definition at line 115 of file GlobalHaloAlgo.h.

Referenced by SetMaxSegmentTheta().

float GlobalHaloAlgo::TowerEtThreshold

private

Definition at line 112 of file GlobalHaloAlgo.h.

Referenced by SetCaloTowerEtThreshold().