94 std::vector<StableProvenance const*> theProvenance;
96 for(
auto const& provenance : theProvenance)
98 cout<<
" Print all module/label names "<<provenance->moduleName()<<
" "<<provenance->moduleLabel()<<
99 " "<<provenance->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++) {
184 cout<<
" Energy ECAL "<<(*hite).energy()<<
185 " eta "<<posE.
eta()<<
" phi "<<posE.
phi()<<endl;
187 energyECAL = energyECAL + (*hite).energy();
197 for (; hith != (hbhe.
product())->
end (); hith++) {
201 cout<<
" Energy HCAL "<<(*hith).energy()<<
202 " eta "<<posH.
eta()<<
" phi "<<posH.
phi()<<endl;
204 energyHCAL = energyHCAL + (*hith).energy();
208 cout<<
" Energy ECAL "<< energyECAL<<
" Energy HCAL "<< energyHCAL<<endl;
216 for (; hito != (ho.
product())->
end (); hito++)
227 cout<<
" we are in GammaJetProd area "<<endl;
230 std::cout<<
" Size of ECAL "<<(*ecal).size()<<std::endl;
234 std::cout<<
" Jet size "<<(*jets).size()<<std::endl;
235 reco::CaloJetCollection::const_iterator
jet = jets->begin ();
236 for (; jet != jets->end (); jet++)
238 cout<<
" Et jet "<<(*jet).et()<<
" eta "<<(*jet).eta()<<
" phi "<<(*jet).phi()<<endl;
243 std::cout<<
" Tracks size "<<(*tracks).size()<<std::endl;
249 std::cout<<
" GammaClus size "<<(*eclus).size()<<std::endl;
250 reco::SuperClusterCollection::const_iterator iclus = eclus->begin ();
251 for (; iclus != eclus->end (); iclus++)
253 cout<<
" Et gamma "<<(*iclus).energy()/cosh((*iclus).eta())<<
" eta "<<(*iclus).eta()<<
" phi "<<(*iclus).phi()<<endl;
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
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_
edm::EDGetTokenT< HORecHitCollection > tok_horeco_
edm::EDGetTokenT< HOCalibVariableCollection > tok_hovar_
edm::EDGetTokenT< EcalRecHitCollection > tok_ecal_
const GlobalPoint & getPosition(const DetId &id) const
Get the position of a given detector id.
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_
T const * product() const
void getAllStableProvenance(std::vector< StableProvenance const * > &provenances) const