61 std::auto_ptr<BeamHaloSummary> TheBeamHaloSummary(
new BeamHaloSummary() );
76 TheBeamHaloSummary->GetCSCHaloReport()[0] = 1;
85 TheBeamHaloSummary->GetCSCHaloReport()[1] = 1;
89 TheBeamHaloSummary->GetCSCHaloReport()[2] = 1;
95 TheBeamHaloSummary->GetCSCHaloReport()[3] = 1;
105 bool EcalLooseId =
false, EcalTightId =
false;
148 float angle = vm_Angle[cluster];
149 float roundness = vm_Roundness[cluster];
167 TheBeamHaloSummary->GetEcalHaloReport()[0] = 1;
169 TheBeamHaloSummary->GetEcalHaloReport()[1] = 1;
177 const std::vector<PhiWedge> HcalWedges = HcalData.
GetPhiWedges();
178 bool HcalLooseId =
false, HcalTightId =
false;
179 for( std::vector<PhiWedge>::const_iterator iWedge = HcalWedges.begin() ; iWedge != HcalWedges.end() ; iWedge++ )
181 bool HcaliPhi =
false;
196 for(
unsigned int i = 0 ;
i < TheBeamHaloSummary->GetHcaliPhiSuspects().size() ;
i++ )
198 if( iWedge->iPhi() == TheBeamHaloSummary->GetHcaliPhiSuspects()[
i] )
205 TheBeamHaloSummary->GetHcaliPhiSuspects().push_back( iWedge->iPhi() ) ;
209 TheBeamHaloSummary->GetHcalHaloReport()[0] = 1;
211 TheBeamHaloSummary->GetHcalHaloReport()[1] = 1;
218 bool GlobalLooseId =
false;
219 bool GlobalTightId =
false;
225 if( MatchedEcalWedges.size() || MatchedHcalWedges.size() )
226 GlobalLooseId =
true;
229 for( std::vector<PhiWedge>::const_iterator iWedge = MatchedEcalWedges.begin() ; iWedge != MatchedEcalWedges.end(); iWedge ++ )
232 GlobalTightId =
true;
234 GlobalTightId =
true;
237 for( std::vector<PhiWedge>::const_iterator iWedge = MatchedHcalWedges.begin() ; iWedge != MatchedHcalWedges.end(); iWedge ++ )
240 GlobalTightId =
true;
242 GlobalTightId =
true;
246 TheBeamHaloSummary->GetGlobalHaloReport()[0] = 1;
248 TheBeamHaloSummary->GetGlobalHaloReport()[1] = 1;
250 iEvent.
put(TheBeamHaloSummary);
int NumberOfHaloTracks(HaloData::Endcap z=HaloData::both) const
edm::ValueMap< float > & GetShowerShapesRoundness()
edm::EDGetTokenT< HcalHaloData > hcalhalodata_token_
float L_HcalPhiWedgeConfidence
float T_EcalPhiWedgeConfidence
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::RefVector< reco::SuperClusterCollection > & GetSuperClusters()
edm::ValueMap< float > & GetShowerShapesAngle()
short int NumberOfOutOfTimeTriggers(HaloData::Endcap z=HaloData::both) const
float L_EcalShowerShapesAngle
std::vector< PhiWedge > & GetMatchedEcalPhiWedges()
short int NFlatHaloSegments() const
int L_EcalSuperClusterSize
bool GetSegmentsInBothEndcaps() const
edm::EDGetTokenT< EcalHaloData > ecalhalodata_token_
edm::EDGetTokenT< GlobalHaloData > globalhalodata_token_
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
int NumberOfHaloTriggers(HaloData::Endcap z=HaloData::both) const
float L_EcalShowerShapesRoundness
float T_EcalSuperClusterEnergy
float T_HcalPhiWedgeEnergy
int T_EcalSuperClusterSize
T const * product() const
float L_HcalPhiWedgeEnergy
std::vector< PhiWedge > & GetMatchedHcalPhiWedges()
edm::EDGetTokenT< CSCHaloData > cschalodata_token_
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