#include <TestMix.h>
Public Member Functions | |
virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
TestMix (const edm::ParameterSet &) | |
~TestMix () | |
Private Attributes | |
int | level_ |
std::vector< std::string > | track_containers2_ |
std::vector< std::string > | track_containers_ |
TestMix::TestMix | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 46 of file TestMix.cc.
References gather_cfg::cout, level_, track_containers2_, and track_containers_.
: level_(iConfig.getUntrackedParameter<int>("PrintLevel")) { std::cout << "Constructed testMix , level "<<level_<<std::endl; track_containers_.push_back("g4SimHitsTrackerHitsTECHighTof"); track_containers_.push_back("g4SimHitsTrackerHitsTECLowTof"); track_containers2_.push_back("g4SimHitsTrackerHitsTECLowTof"); track_containers2_.push_back("g4SimHitsTrackerHitsTECHighTof"); }
TestMix::~TestMix | ( | ) |
Definition at line 60 of file TestMix.cc.
{ // do anything here that needs to be done at desctruction time // (e.g. close files, deallocate resources etc.) }
void TestMix::analyze | ( | const edm::Event & | iEvent, |
const edm::EventSetup & | iSetup | ||
) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 75 of file TestMix.cc.
References MixCollection< T >::begin(), MixCollection< T >::MixItr::bunch(), EncodedEventId::bunchCrossing(), prof2calltree::count, gather_cfg::cout, MixCollection< T >::end(), EncodedEventId::event(), CoreSimVertex::eventId(), edm::Event::getByLabel(), edm::HepMCProduct::GetEvent(), MixCollection< T >::MixItr::getSourceType(), MixCollection< T >::MixItr::getTrigger(), SimVertex::parentIndex(), track_containers2_, and track_containers_.
{ using namespace edm; bool got; int count=0; // test accesses to CrossingFrame // attention: operator-> returns the templated object, but // bunch() and getTrigger() are methods of the iterator itself! // test access to SimTracks directly in CrossingFrame edm::Handle<CrossingFrame<SimTrack> > cf_simtrack; bool gotTracks = iEvent.getByLabel("mix","g4SimHits",cf_simtrack); if (!gotTracks) std::cout<<" Could not read SimTracks!!!!"<<std::endl; // not pointer compatible!!!! // if (gotTracks) { // std::vector<SimTrack>::const_iterator first; // std::vector<SimTrack>::const_iterator last; // cf_simtrack->getPileups(first,last); // unsigned int ic=0; // for (std::vector<SimTrack>::const_iterator it=first;it!=last;it++) { // std::cout<<" pileup SimTrack "<<ic<<" has genpart index "<<(*it).genpartIndex()<<" vertex Index "<<(*it).vertIndex() <<" bunchcrossing "<<cf_simtrack->getBunchCrossing(ic)<<std::endl; // ic++; // } // } const std::string subdet("g4SimHitsTrackerHitsTECHighTof"); edm::Handle<CrossingFrame<PSimHit> > cf_simhit; got = iEvent.getByLabel("mix",subdet,cf_simhit); if (!got) std::cout<<" Could not read SimHits with label "<<subdet<<"!!!!"<<std::endl; else { std::cout<<"\n\n=================== Starting SimHit access, subdet "<<subdet<<" ==================="<<std::endl; std::auto_ptr<MixCollection<PSimHit> > col(new MixCollection<PSimHit>(cf_simhit.product(),std::pair<int,int>(-1,1))); std::cout<<*(col.get())<<std::endl; MixCollection<PSimHit>::iterator cfi; for (cfi=col->begin(); cfi!=col->end();cfi++) { std::cout<<" Hit "<<count<<" has tof "<<cfi->timeOfFlight()<<" trackid "<<cfi->trackId() <<" bunchcr "<<cfi.bunch()<<" trigger "<<cfi.getTrigger()<<", from EncodedEventId: "<<cfi->eventId().bunchCrossing()<<" " <<cfi->eventId().event() <<" bcr from MixCol "<<cfi.bunch()<<std::endl; // std::cout<<" Hit: "<<(*cfi)<<std::endl; count++; } } // test access to CaloHits const std::string subdetcalo("g4SimHitsEcalHitsEB"); edm::Handle<CrossingFrame<PCaloHit> > cf_calo; got = iEvent.getByLabel("mix",subdetcalo,cf_calo); if (!got) std::cout<<" Could not read CaloHits with label "<<subdetcalo<<"!!!!"<<std::endl; else { std::cout<<"\n\n=================== Starting CaloHit access, subdet "<<subdetcalo<<" ==================="<<std::endl; std::auto_ptr<MixCollection<PCaloHit> > colcalo(new MixCollection<PCaloHit>(cf_calo.product(), std::pair<int,int>(-1,1))); std::cout<<*(colcalo.get())<<std::endl; MixCollection<PCaloHit>::iterator cficalo; count=0; for (cficalo=colcalo->begin(); cficalo!=colcalo->end();cficalo++) { std::cout<<" CaloHit "<<count<<" has tof "<<cficalo->time()<<" trackid "<<cficalo->geantTrackId() <<" bunchcr "<<cficalo.bunch()<<" trigger "<<cficalo.getTrigger()<<", from EncodedEventId: "<<cficalo->eventId().bunchCrossing()<<" " <<cficalo->eventId().event() <<std::endl; // std::cout<<" Calo Hit: "<<(*cficalo)<<std::endl; count++; } } // test access to SimTracks if (gotTracks) { std::cout<<"\n=================== Starting SimTrack access ==================="<<std::endl; // edm::Handle<CrossingFrame<SimTrack> > cf_simtrack; // iEvent.getByLabel("mix",cf_simtrack); std::auto_ptr<MixCollection<SimTrack> > col2(new MixCollection<SimTrack>(cf_simtrack.product())); MixCollection<SimTrack>::iterator cfi2; int count2=0; std::cout <<" \nWe got "<<col2->sizeSignal()<<" signal tracks and "<<col2->sizePileup()<<" pileup tracks, total: "<<col2->size()<<std::endl; for (cfi2=col2->begin(); cfi2!=col2->end();cfi2++) { std::cout<<" SimTrack "<<count2<<" has genpart index "<<cfi2->genpartIndex()<<" vertex Index "<<cfi2->vertIndex() <<" bunchcr "<<cfi2.bunch()<<" trigger "<<cfi2.getTrigger()<<", from EncodedEventId: "<<cfi2->eventId().bunchCrossing() <<" "<<cfi2->eventId().event() <<std::endl; count2++; } } // test access to SimVertices edm::Handle<CrossingFrame<SimVertex> > cf_simvtx; got = iEvent.getByLabel("mix","g4SimHits",cf_simvtx); if (!got) std::cout<<" Could not read Simvertices !!!!"<<std::endl; else { std::cout<<"\n=================== Starting SimVertex access ==================="<<std::endl; std::auto_ptr<MixCollection<SimVertex> > col3(new MixCollection<SimVertex>(cf_simvtx.product())); MixCollection<SimVertex>::iterator cfi3; int count3=0; std::cout <<" \nWe got "<<col3->sizeSignal()<<" signal vertices and "<<col3->sizePileup()<<" pileup vertices, total: "<<col3->size()<<std::endl; for (cfi3=col3->begin(); cfi3!=col3->end();cfi3++) { std::cout<<" SimVertex "<<count3<<" has parent index "<<cfi3->parentIndex()<<" bunchcr "<<cfi3.bunch()<<" trigger "<<cfi3.getTrigger()<<", from EncodedEventId: "<<cfi3->eventId().bunchCrossing() <<" "<<cfi3->eventId().event() <<std::endl; SimVertex myvtx=(*cfi3); std::cout<<"Same with op*: "<<count3<<" has parent index "<<myvtx.parentIndex()<<" bunchcr "<<cfi3.bunch()<<" trigger "<<cfi3.getTrigger()<<", from EncodedEventId: "<<myvtx.eventId().bunchCrossing() <<" "<<myvtx.eventId().event() <<std::endl; count3++; } } //test MixCollection constructor with several subdetector names bool got1,got2=false; std::auto_ptr<MixCollection<PSimHit> > all_trackhits; std::auto_ptr<MixCollection<PSimHit> > all_trackhits2; std::cout<<"\n=================== Starting test for coll of several ROU-s ==================="<<std::endl; // edm::Handle<CrossingFrame<PSimHit> > cf_simhit; std::vector<const CrossingFrame<PSimHit> *> cfvec; got1 = iEvent.getByLabel("mix",track_containers_[0],cf_simhit); if (!got1) std::cout<<" Could not read SimHits with label "<<track_containers_[0]<<"!!!!"<<std::endl; else { std::cout<<"\n=================== Starting test for coll of several ROU-s ==================="<<std::endl; cfvec.push_back(cf_simhit.product()); std::cout <<" \nFirst container "<<track_containers_[0]<<" Nr signals "<<cf_simhit->getNrSignals() << ", Nr pileups "<<cf_simhit->getNrPileups() <<std::endl; got2 = iEvent.getByLabel("mix",track_containers_[1],cf_simhit); if (got2) { cfvec.push_back(cf_simhit.product()); std::cout <<" \nSecond container "<<track_containers_[1]<<" Nr signals "<<cf_simhit->getNrSignals() << ", Nr pileups "<<cf_simhit->getNrPileups() <<std::endl; all_trackhits= std::auto_ptr<MixCollection<PSimHit> >(new MixCollection<PSimHit>(cfvec)); std::cout <<" \nFor all containers we got "<<all_trackhits->sizeSignal()<<" signal hits and "<<all_trackhits->sizePileup()<<" pileup hits, total: "<<all_trackhits->size()<<std::endl; MixCollection<PSimHit>::iterator it; int ii=0; for (it=all_trackhits->begin(); it!= all_trackhits->end();it++) { std::cout<<" Hit "<<ii<<" of all hits has tof "<<it->timeOfFlight()<<" trackid "<<it->trackId() <<" bunchcr "<<it.bunch()<<" trigger "<<it.getTrigger()<<", from EncodedEventId: "<<it->eventId().bunchCrossing() <<" "<<it->eventId().event()<<std::endl; ii++; } } } //test the same in different order: should be the same sizes, different order MixCollection<PSimHit>::iterator it2; int ii2=0; std::vector<const CrossingFrame<PSimHit> *> cfvec2; got = iEvent.getByLabel("mix",track_containers2_[0],cf_simhit); if (!got) std::cout<<" Could not read SimHits with label "<<track_containers2_[0]<<"!!!!"<<std::endl; else { cfvec2.push_back(cf_simhit.product()); got2 = iEvent.getByLabel("mix",track_containers2_[1],cf_simhit); if (got2) { cfvec2.push_back(cf_simhit.product()); all_trackhits2= std::auto_ptr<MixCollection<PSimHit> > (new MixCollection<PSimHit>(cfvec2)); std::cout <<" \nSame containers, different order: we got "<<all_trackhits2->sizeSignal()<<" signal hits and "<<all_trackhits2->sizePileup()<<" pileup hits, total: "<<all_trackhits2->size()<<std::endl; for (it2=all_trackhits2->begin(); it2!= all_trackhits2->end();it2++) { std::cout<<" Hit "<<ii2<<" of all hits has tof "<<it2->timeOfFlight()<<" trackid "<<it2->trackId() <<" bunchcr "<<it2.bunch()<<" trigger "<<it2.getTrigger()<<", bcr from Id: "<<it2->eventId().bunchCrossing() <<" evtnr in id "<<it2->eventId().event()<<std::endl; ii2++; } } } //test MixCollection for HepMCProducts //at the same time test getObject method //we should have each line twice //------------------------------------ edm::Handle<CrossingFrame<edm::HepMCProduct> > cf_hepmc; got = iEvent.getByLabel("mix","generator",cf_hepmc); if (!got) std::cout<<" Could not read HepMCProducts!!!!"<<std::endl; else { std::auto_ptr<MixCollection<edm::HepMCProduct> > colhepmc(new MixCollection<edm::HepMCProduct>(cf_hepmc.product())); MixCollection<edm::HepMCProduct>::iterator cfihepmc; int counthepmc=0; std::cout <<" \nWe got "<<colhepmc->sizeSignal()<<" signal hepmc products and "<<colhepmc->sizePileup()<<" pileup hepmcs, total: "<<colhepmc->size()<<std::endl; for (cfihepmc=colhepmc->begin(); cfihepmc!=colhepmc->end();cfihepmc++) { std::cout<<" edm::HepMCProduct "<<counthepmc<<" has event number "<<cfihepmc->GetEvent()->event_number()<<", "<< cfihepmc->GetEvent()->particles_size()<<" particles and "<<cfihepmc->GetEvent()->vertices_size()<<" vertices, bunchcr= "<<cfihepmc.bunch()<<" trigger= "<<cfihepmc.getTrigger() <<" sourcetype= "<<cfihepmc.getSourceType()<<std::endl; HepMCProduct myprod=colhepmc->getObject(counthepmc); std::cout<<"same with getObject:hepmc product "<<counthepmc<<" has event number "<<myprod.GetEvent()->event_number()<<", "<<myprod.GetEvent()->particles_size()<<" particles and "<<myprod.GetEvent()->vertices_size()<<" vertices"<<std::endl; counthepmc++; } } //---------------------------------------------------------------------------- // testing special situations //---------------------------------------------------------------------------- if (got2) { // test reusage of the same iterator int ii3=0; for (it2=all_trackhits2->begin(); it2!= all_trackhits2->end();it2++) ii3++; if (ii3!=ii2) std::cout<<" Problem when re-using iterator!!"<<std::endl; else std::cout<<" \nNo problem when re-using iterator."<<std::endl; } // test access to non-filled collections //cases: 0) ok, collection has members // 1) bunchrange given outside of existent bunchcrossing numbers ==>exc std::cout<<"\n=================== Starting tests for abnormal conditions ==================="<<std::endl; // test case 0 if (got1) { std::cout<<"\n[ Testing abnormal conditions case 0]Should be all ok: registry: "<<all_trackhits->inRegistry()<<" size: "<<all_trackhits->size()<<std::endl; // test case 1 std::cout<<"\n[ Testing abnormal conditions case 1] Should throw an exception " <<std::endl; MixCollection<PSimHit> * col21=0; col21=new MixCollection<PSimHit>(cf_simhit.product(),std::pair<int,int>(-10,20)); } }
int edm::TestMix::level_ [private] |
std::vector<std::string> edm::TestMix::track_containers2_ [private] |
std::vector<std::string> edm::TestMix::track_containers_ [private] |