55 std::auto_ptr<BeamHaloSummary> TheBeamHaloSummary(
new BeamHaloSummary() );
69 TheBeamHaloSummary->GetCSCHaloReport()[0] = 1;
78 TheBeamHaloSummary->GetCSCHaloReport()[1] = 1;
82 TheBeamHaloSummary->GetCSCHaloReport()[2] = 1;
88 TheBeamHaloSummary->GetCSCHaloReport()[3] = 1;
97 bool EcalLooseId =
false, EcalTightId =
false;
140 float angle = vm_Angle[cluster];
141 float roundness = vm_Roundness[cluster];
159 TheBeamHaloSummary->GetEcalHaloReport()[0] = 1;
161 TheBeamHaloSummary->GetEcalHaloReport()[1] = 1;
167 const std::vector<PhiWedge> HcalWedges = HcalData.
GetPhiWedges();
168 bool HcalLooseId =
false, HcalTightId =
false;
169 for( std::vector<PhiWedge>::const_iterator iWedge = HcalWedges.begin() ; iWedge != HcalWedges.end() ; iWedge++ )
171 bool HcaliPhi =
false;
186 for(
unsigned int i = 0 ;
i < TheBeamHaloSummary->GetHcaliPhiSuspects().size() ;
i++ )
188 if( iWedge->iPhi() == TheBeamHaloSummary->GetHcaliPhiSuspects()[
i] )
195 TheBeamHaloSummary->GetHcaliPhiSuspects().push_back( iWedge->iPhi() ) ;
199 TheBeamHaloSummary->GetHcalHaloReport()[0] = 1;
201 TheBeamHaloSummary->GetHcalHaloReport()[1] = 1;
206 bool GlobalLooseId =
false;
207 bool GlobalTightId =
false;
213 if( MatchedEcalWedges.size() || MatchedHcalWedges.size() )
214 GlobalLooseId =
true;
217 for( std::vector<PhiWedge>::const_iterator iWedge = MatchedEcalWedges.begin() ; iWedge != MatchedEcalWedges.end(); iWedge ++ )
220 GlobalTightId =
true;
222 GlobalTightId =
true;
225 for( std::vector<PhiWedge>::const_iterator iWedge = MatchedHcalWedges.begin() ; iWedge != MatchedHcalWedges.end(); iWedge ++ )
228 GlobalTightId =
true;
230 GlobalTightId =
true;
234 TheBeamHaloSummary->GetGlobalHaloReport()[0] = 1;
236 TheBeamHaloSummary->GetGlobalHaloReport()[1] = 1;
238 iEvent.
put(TheBeamHaloSummary);
int NumberOfHaloTracks(HaloData::Endcap z=HaloData::both) const
edm::ValueMap< float > & GetShowerShapesRoundness()
float L_HcalPhiWedgeConfidence
float T_EcalPhiWedgeConfidence
edm::InputTag IT_GlobalHaloData
edm::RefVector< reco::SuperClusterCollection > & GetSuperClusters()
edm::ValueMap< float > & GetShowerShapesAngle()
edm::InputTag IT_CSCHaloData
short int NumberOfOutOfTimeTriggers(HaloData::Endcap z=HaloData::both) const
float L_EcalShowerShapesAngle
edm::InputTag IT_HcalHaloData
std::vector< PhiWedge > & GetMatchedEcalPhiWedges()
short int NFlatHaloSegments() const
int L_EcalSuperClusterSize
bool GetSegmentsInBothEndcaps() const
int T_EcalPhiWedgeConstituents
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
short int NTracksSmalldT() const
Abs< T >::type abs(const T &t)
short int NOutOfTimeHits() const
float T_EcalShowerShapesAngle
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
int NumberOfHaloTriggers(HaloData::Endcap z=HaloData::both) const
float L_EcalShowerShapesRoundness
float T_EcalSuperClusterEnergy
float T_HcalPhiWedgeEnergy
edm::InputTag IT_EcalHaloData
int T_EcalSuperClusterSize
T const * product() const
float L_HcalPhiWedgeEnergy
std::vector< PhiWedge > & GetMatchedHcalPhiWedges()
size_type size() const
Size of the RefVector.
float T_HcalPhiWedgeConfidence
const std::vector< PhiWedge > & GetPhiWedges() const
float T_EcalShowerShapesRoundness
float L_EcalSuperClusterEnergy
int T_HcalPhiWedgeConstituents
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
int L_HcalPhiWedgeConstituents