#include <TestMix.h>

Inheritance diagram for edm::TestMix:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

	TestMix (const edm::ParameterSet &)

	~TestMix () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::EDGetTokenT< CrossingFrame< PCaloHit > >	CaloToken1_

edm::EDGetTokenT< CrossingFrame< HepMCProduct > >	HepMCToken_

int	level_

edm::EDGetTokenT< CrossingFrame< SimTrack > >	SimTrackToken_

edm::EDGetTokenT< CrossingFrame< SimVertex > >	SimVertexToken_

std::vector< std::string >	track_containers2_

std::vector< std::string >	track_containers_

edm::EDGetTokenT< CrossingFrame< PSimHit > >	TrackerToken0_

edm::EDGetTokenT< CrossingFrame< PSimHit > >	TrackerToken1_

edm::EDGetTokenT< CrossingFrame< PSimHit > >	TrackerToken2_

edm::EDGetTokenT< CrossingFrame< PSimHit > >	TrackerToken3_

edm::EDGetTokenT< CrossingFrame< PSimHit > >	TrackerToken4_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 43 of file TestMix.h.

Constructor & Destructor Documentation

TestMix::TestMix ( const edm::ParameterSet & iConfig )

explicit

Definition at line 37 of file TestMix.cc.

References CaloToken1_, gather_cfg::cout, HepMCToken_, level_, SimTrackToken_, SimVertexToken_, GlobalPosition_Frontier_DevDB_cff::tag, track_containers2_, track_containers_, TrackerToken0_, TrackerToken1_, TrackerToken2_, TrackerToken3_, and TrackerToken4_.

                                               : 
   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");
 
   edm::InputTag tag = edm::InputTag("mix","g4SimHits");
   
   SimTrackToken_ = consumes<CrossingFrame<SimTrack>>(tag);
   SimVertexToken_ = consumes<CrossingFrame<SimVertex>>(tag);
 
   tag = edm::InputTag("mix","g4SimHitsTrackerHitsTECHighTof");
   TrackerToken0_ = consumes<CrossingFrame<PSimHit>>(tag);
 
   tag = edm::InputTag("mix","g4SimHitsEcalHitsEB");
   CaloToken1_ = consumes<CrossingFrame<PCaloHit>>(tag);
 
   tag = edm::InputTag("mix",track_containers_[0]);
   TrackerToken1_ = consumes<CrossingFrame<PSimHit>>(tag);
 
   tag = edm::InputTag("mix",track_containers_[1]);
   TrackerToken2_ = consumes<CrossingFrame<PSimHit>>(tag);
 
   tag = edm::InputTag("mix",track_containers2_[0]);
   TrackerToken3_ = consumes<CrossingFrame<PSimHit>>(tag);
 
   tag = edm::InputTag("mix",track_containers2_[1]);
   TrackerToken4_ = consumes<CrossingFrame<PSimHit>>(tag);
 
   tag = edm::InputTag("mix","generatorSmeared");
   HepMCToken_ = consumes<CrossingFrame<HepMCProduct>>(tag);
 
 }

TestMix::~TestMix ( )

override

Definition at line 77 of file TestMix.cc.

 {
  
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void TestMix::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 92 of file TestMix.cc.

References MixCollection< T >::MixItr::bunch(), EncodedEventId::bunchCrossing(), CaloToken1_, cuy::col, KineDebug3::count(), gather_cfg::cout, EncodedEventId::event(), CoreSimVertex::eventId(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), MixCollection< T >::MixItr::getSourceType(), MixCollection< T >::MixItr::getTrigger(), HepMCToken_, cuy::ii, SimVertex::parentIndex(), edm::Handle< T >::product(), SimTrackToken_, SimVertexToken_, AlCaHLTBitMon_QueryRunRegistry::string, track_containers2_, track_containers_, TrackerToken0_, TrackerToken1_, TrackerToken2_, TrackerToken3_, and TrackerToken4_.

 {
   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.getByToken(SimTrackToken_,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.getByToken(TrackerToken0_,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::unique_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.getByToken(CaloToken1_,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::unique_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::unique_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.getByToken(SimVertexToken_,cf_simvtx);
   if (!got) std::cout<<" Could not read Simvertices !!!!"<<std::endl;
   else {
     std::cout<<"\n=================== Starting SimVertex access ==================="<<std::endl;
     std::unique_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::unique_ptr<MixCollection<PSimHit> > all_trackhits;
   std::unique_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.getByToken(TrackerToken1_,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.getByToken(TrackerToken2_,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::unique_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.getByToken(TrackerToken3_,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.getByToken(TrackerToken4_,cf_simhit);
     if (got2) {
       cfvec2.push_back(cf_simhit.product());
       all_trackhits2= std::unique_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.getByToken(HepMCToken_,cf_hepmc);
   if (!got) std::cout<<" Could not read HepMCProducts!!!!"<<std::endl;
   else {
     std::unique_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=nullptr;
     col21=new MixCollection<PSimHit>(cf_simhit.product(),std::pair<int,int>(-10,20));
     delete col21;
   }
 }