93 std::vector<Provenance const*> theProvenance;
95 for( std::vector<Provenance const*>::const_iterator ip = theProvenance.begin();
96 ip != theProvenance.end(); ip++)
98 cout<<
" Print all module/label names "<<
moduleName(**ip)<<
" "<<(**ip).moduleLabel()<<
99 " "<<(**ip).productInstanceName()<<endl;
120 std::cout<<
" Size of muon collection "<<mucand->size()<<std::endl;
121 for(MuonCollection::const_iterator it = mucand->begin(); it != mucand->end(); it++)
124 std::cout<<
" Pt muon "<<mu->innerMomentum()<<std::endl;
150 cout<<
" We are here "<<endl;
155 std::cout<<
" Tracks size "<<(*tracks).size()<<std::endl;
156 reco::TrackCollection::const_iterator track = tracks->begin ();
158 for (; track != tracks->end (); track++)
160 cout<<
" P track "<<(*track).p()<<
" eta "<<(*track).eta()<<
" phi "<<(*track).phi()<<
" Outer "<<(*track).outerMomentum()<<
" "<<
161 (*track).outerPosition()<<endl;
163 cout<<
" Track extra "<<myextra->outerMomentum()<<
" "<<myextra->outerPosition()<<endl;
175 for (; hite != (ecal.
product())->
end (); hite++)
183 GlobalPoint posE = geo->getPosition((*hite).detid());
185 cout<<
" Energy ECAL "<<(*hite).energy()<<
186 " eta "<<posE.
eta()<<
" phi "<<posE.
phi()<<endl;
188 energyECAL = energyECAL + (*hite).energy();
198 for (; hith != (hbhe.
product())->
end (); hith++)
201 GlobalPoint posH = geo->getPosition((*hith).detid());
203 cout<<
" Energy HCAL "<<(*hith).energy()<<
204 " eta "<<posH.
eta()<<
" phi "<<posH.
phi()<<endl;
206 energyHCAL = energyHCAL + (*hith).energy();
210 cout<<
" Energy ECAL "<< energyECAL<<
" Energy HCAL "<< energyHCAL<<endl;
218 for (; hito != (ho.
product())->
end (); hito++)
229 cout<<
" we are in GammaJetProd area "<<endl;
232 std::cout<<
" Size of ECAL "<<(*ecal).size()<<std::endl;
236 std::cout<<
" Jet size "<<(*jets).size()<<std::endl;
237 reco::CaloJetCollection::const_iterator
jet = jets->begin ();
238 for (; jet != jets->end (); jet++)
240 cout<<
" Et jet "<<(*jet).et()<<
" eta "<<(*jet).eta()<<
" phi "<<(*jet).phi()<<endl;
245 std::cout<<
" Tracks size "<<(*tracks).size()<<std::endl;
251 std::cout<<
" GammaClus size "<<(*eclus).size()<<std::endl;
252 reco::SuperClusterCollection::const_iterator iclus = eclus->begin ();
253 for (; iclus != eclus->end (); iclus++)
255 cout<<
" Et gamma "<<(*iclus).energy()/cosh((*iclus).eta())<<
" eta "<<(*iclus).eta()<<
" phi "<<(*iclus).phi()<<endl;
void getAllProvenance(std::vector< Provenance const * > &provenances) const
double energyECAL(std::vector< DetId > &vdets, edm::Handle< T > &hitsEB, edm::Handle< T > &hitsEE, double ebThr=-100, double eeThr=-100, double tMin=-500, double tMax=500, bool debug=false)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::EDGetTokenT< reco::TrackCollection > tok_tracks_
Geom::Phi< T > phi() const
std::vector< EcalRecHit >::const_iterator const_iterator
edm::EDGetTokenT< reco::MuonCollection > tok_muons_
edm::EDGetTokenT< reco::SuperClusterCollection > tok_gamma_
edm::EDGetTokenT< HORecHitCollection > tok_hoProd_
edm::EDGetTokenT< HORecHitCollection > tok_ho_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheProd_
std::string moduleName(Provenance const &provenance)
edm::EDGetTokenT< HORecHitCollection > tok_horeco_
edm::EDGetTokenT< HOCalibVariableCollection > tok_hovar_
edm::EDGetTokenT< EcalRecHitCollection > tok_ecal_
double energyHCAL(std::vector< DetId > &vdets, edm::Handle< T > &hits, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, bool useRaw=false, bool debug=false)
T const * product() const
std::vector< HOCalibVariables > HOCalibVariableCollection
collection of HOcalibration variabale
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
edm::EDGetTokenT< reco::CaloJetCollection > tok_jets_
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
tuple size
Write out results.