189 std::vector<Provenance const*> theProvenance;
201 (*myout_part)<<
"Event "<<iEvent.
id().
run()<<
" "<<iEvent.
id().
event()<<endl;
202 (*myout_jet)<<
"Event "<<iEvent.
id().
run()<<
" "<<iEvent.
id().
event()<<endl;
203 (*myout_hcal)<<
"Event "<<iEvent.
id().
run()<<
" "<<iEvent.
id().
event()<<endl;
204 (*myout_ecal)<<
"Event "<<iEvent.
id().
run()<<
" "<<iEvent.
id().
event()<<endl;
205 (*myout_photon)<<
"Event "<<iEvent.
id().
run()<<
" "<<iEvent.
id().
event()<<endl;
208 std::vector<edm::InputTag>::const_iterator ic;
212 double etlost = -100.1;
220 reco::CaloJetCollection::const_iterator
jet = jets->begin ();
221 cout<<
" Size of Calo jets "<<jets->size()<<endl;
224 if(jets->size() > 0 )
227 for (; jet != jets->end (); jet++)
229 cout<<
" Jet et "<<(*jet).et()<<
" "<<(*jet).eta()<<
" "<<(*jet).phi()<<endl;
231 if(ij<4) (*myout_jet)<<jettype<<
" "<<
reco<<
" "<<ij<<
" "<<(*jet).et()<<
" "<<(*jet).eta()<<
" "<<(*jet).phi()
232 <<
" "<<iEvent.
id().
event()<<endl;
246 cout<<
" Calojets are missed "<<endl;
251 cout<<
" We filled CaloJet part "<<jetexist<<endl;
253 if( jetexist < 0 ) (*myout_jet)<<jetexist<<
" "<<
reco<<
" "<<etlost
254 <<
" "<<etlost<<
" "<<etlost
255 <<
" "<<iEvent.
id().
event()<<endl;
258 std::vector<edm::InputTag>::const_iterator
i;
259 vector<std::pair<DetId, double> > theRecHits;
267 recHit != (*ec).end(); ++recHit)
272 theRecHits.push_back(std::pair<DetId, double>(recHit->detid(), recHit->energy()));
274 if( (*recHit).energy()>
ecut[recHit->detid().subdetId()-1][0] )
275 (*
myout_ecal)<<recHit->detid().subdetId()<<
" "<<(*recHit).energy()<<
" "<<pos.
phi()<<
" "<<pos.
eta()
276 <<
" "<<iEvent.
id().
event()<<endl;
282 cout<<
" Ecal collection is missed "<<endl;
287 cout<<
" Fill EcalRecHits "<<endl;
297 hbheItr != (*hbhe).end(); ++hbheItr)
301 (*myout_hcal)<<
id.subdetId()<<
" "<<(*hbheItr).energy()<<
" "<<pos.
phi()<<
302 " "<<pos.
eta()<<
" "<<iEvent.
id().
event()<<endl;
303 theRecHits.push_back(std::pair<DetId, double>(hbheItr->detid(), hbheItr->energy()));
308 cout<<
" HBHE collection is missed "<<endl;
314 for(
int i = 0;
i<9;
i++)
331 for(reco::SuperClusterCollection::const_iterator aClus = correctedSuperClusters->begin();
332 aClus != correctedSuperClusters->end(); aClus++) {
333 double vet = aClus->energy()/cosh(aClus->eta());
334 cout<<
" Supercluster " << ij<<
" Et "<< vet <<
" energy "<<aClus->energy()<<
" eta "<<aClus->eta()<<
" Cut "<<
CutOnEgammaEnergy_<<endl;
338 float gammaiso_ecal[9] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
339 for(vector<std::pair<DetId, double> >::const_iterator it = theRecHits.begin(); it != theRecHits.end(); it++)
344 double deta = fabs(eta-aClus->eta());
345 double dphi = fabs(phi-aClus->phi());
346 if(dphi>4.*atan(1.)) dphi = 8.*atan(1.)-dphi;
347 double dr =
sqrt(deta*deta+dphi*dphi);
350 if( fabs(aClus->eta()) > 1.47 ) rmin = 0.07*(fabs(aClus->eta())-.47)*1.2;
351 if( fabs(aClus->eta()) > 2.2 ) rmin = 0.07*(fabs(aClus->eta())-.47)*1.4;
355 for (
int i = 0;
i<3;
i++)
357 for (
int j = 0;
j<3;
j++)
362 if(it->first.subdetId() == 1) ecutn =
ecut[0][
j];
363 if(it->first.subdetId() == 2) ecutn =
ecut[1][
j];
364 if( dr>rmin && dr<
risol[
i])
366 if((*it).second > ecutn) gammaiso_ecal[itype_ecal] = gammaiso_ecal[itype_ecal]+(*it).second/cosh(eta);
373 if( dr>rmin && dr<
risol[i])
375 if((*it).first > ecutn)
377 gammaiso_ecal[itype_ecal] = gammaiso_ecal[itype_ecal]+(*it).second/cosh(eta);
392 for (
int ii = 0;
ii<9 ;
ii++)
402 (*myout_photon)<<ij<<
" "<<
barrel<<
" "<<vet<<
" "<<aClus->eta()<<
" "<<aClus->phi()<<
" "<<iEvent.
id().
event()<<endl;
403 (*myout_photon)<<ij<<
" "<<gammaiso_ecal[0]<<
" "<<gammaiso_ecal[1] <<
" "<<gammaiso_ecal[2]<<
" "<<gammaiso_ecal[3]
404 <<
" "<<gammaiso_ecal[4]<<
" "<<gammaiso_ecal[5]<<
" "<<gammaiso_ecal[6]<<
" "<<gammaiso_ecal[7]<<
" "<<gammaiso_ecal[8]<<endl;
411 cout<<
" Ecal barrel clusters are missed "<<endl;
416 cout<<
" After iso cuts "<<jetexist<<endl;
418 double ecluslost = -100.1;
419 if(jetexist<0) (*myout_photon)<<jetexist<<
" "<<
barrel<<
" "<<ecluslost<<
" "<<ecluslost
420 <<
" "<<ecluslost<<
" "<<iEvent.
id().
event()<<endl;
422 cout<<
" Event is ready "<<endl;
EventNumber_t event() const
void getAllProvenance(std::vector< Provenance const * > &provenances) const
bool getByToken(EDGetToken token, Handle< PROD > &result) const
double CutOnEgammaEnergy_
Geom::Phi< T > phi() const
std::vector< EcalRecHit >::const_iterator const_iterator
edm::EDGetTokenT< EcalRecHitCollection > tok_ecal_
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
std::ofstream * myout_ecal
const GlobalPoint & getPosition(const DetId &id) const
Get the position of a given detector id.
T const * product() const
T const * product() const
edm::EDGetTokenT< reco::CaloJetCollection > tok_jets_
edm::EDGetTokenT< reco::SuperClusterCollection > tok_egamma_
float GammaIsoEcal[9][20]
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_