#include <TestMixedSource.h>

Inheritance diagram for edm::TestMixedSource:

Public Member Functions
	TestMixedSource (const edm::ParameterSet &)

	~TestMixedSource ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

virtual void	beginJob ()

virtual void	endJob ()

Private Attributes
int	bunchcr_

edm::EDGetTokenT < CrossingFrame< PCaloHit > >	CaloToken1_

std::string	fileName_

edm::EDGetTokenT < CrossingFrame < edm::HepMCProduct > >	HepMCToken_

TFile *	histFile_

TH1I *	histHepMCProduct_bunch_

TH1I *	histPCaloHit_bunch_

TH1I *	histPSimHit_bunchPileups_MuonCSCHits_

TH1I *	histPSimHit_bunchPileups_TrackerHitsTECHighTof_

TH1I *	histPSimHit_bunchSignal_MuonCSCHits_

TH1I *	histPSimHit_bunchSignal_TrackerHitsTECHighTof_

TH1I *	histTrack_bunchPileups_

TH1I *	histTrack_bunchSignal_

TH1I *	histVertex_bunch_

int	maxbunch_

int	minbunch_

edm::EDGetTokenT < CrossingFrame< PSimHit > >	MuonToken_

std::ofstream	outputFile

edm::EDGetTokenT < CrossingFrame< SimTrack > >	SimTrackToken_

edm::EDGetTokenT < CrossingFrame< SimVertex > >	SimVertexToken_

TH1I *	tofhist_

TH1I *	tofhist_sig_

edm::EDGetTokenT < CrossingFrame< PSimHit > >	TrackerToken0_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 45 of file TestMixedSource.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 55 of file TestMixedSource.cc.

References CaloToken1_, HepMCToken_, histHepMCProduct_bunch_, histPCaloHit_bunch_, histPSimHit_bunchPileups_MuonCSCHits_, histPSimHit_bunchPileups_TrackerHitsTECHighTof_, histPSimHit_bunchSignal_MuonCSCHits_, histPSimHit_bunchSignal_TrackerHitsTECHighTof_, histTrack_bunchPileups_, histTrack_bunchSignal_, histVertex_bunch_, HLT_25ns14e33_v1_cff::InputTag, maxbunch_, minbunch_, MuonToken_, SimTrackToken_, SimVertexToken_, GlobalPosition_Frontier_DevDB_cff::tag, tofhist_, tofhist_sig_, and TrackerToken0_.

 : fileName_(iConfig.getParameter<std::string>("fileName")), minbunch_(iConfig.getParameter<int>("minBunch")),maxbunch_(iConfig.getParameter<int>("maxBunch"))
 { 
   
   histTrack_bunchSignal_ = new TH1I("histoTrackSignal","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
   histTrack_bunchPileups_ = new TH1I("histoTrackPileups","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
 
   // Vertex
   histVertex_bunch_ = new TH1I("histoVertex","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
   
   // PCaloHit
   histPCaloHit_bunch_ = new TH1I("histoPCaloHit","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
   
   // PSimHit
   histPSimHit_bunchSignal_TrackerHitsTECHighTof_ = new TH1I("histoPSimHitTrackerHitsTECHighTofSignal","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
   histPSimHit_bunchPileups_TrackerHitsTECHighTof_ = new TH1I("histoPSimHitTrackerHitsTECHighTofPileups","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
   histPSimHit_bunchSignal_MuonCSCHits_ = new TH1I("histoPSimHitMuonCSCHitsSignal","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
   histPSimHit_bunchPileups_MuonCSCHits_ = new TH1I("histoPSimHitMuonCSCHitsPileups","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
 
   int bsp = 25;//bunchspace
   tofhist_ = new TH1I ("TrackerHit_Tof_bcr","TrackerHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
   tofhist_sig_ = new TH1I("SignalTrackerHit_Tof_bcr","TrackerHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
 
 
   // HepMCProduct
   histHepMCProduct_bunch_ = new TH1I("histoHepMCProduct","Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
 
   // Tokens
 
   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","g4SimHitsMuonCSCHits");
   MuonToken_ = consumes<CrossingFrame<PSimHit>>(tag);
 
   tag = edm::InputTag("mix","generator");
   HepMCToken_ = consumes<CrossingFrame<edm::HepMCProduct>>(tag);
 
 
 }

TestMixedSource::~TestMixedSource ( )

Definition at line 105 of file TestMixedSource.cc.

References gather_cfg::cout.

 { std::cout << " Destructor TestMixed"  << std::endl;
 
 }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 117 of file TestMixedSource.cc.

References a, b, MixCollection< T >::MixItr::bunch(), EnergyCorrector::c, CaloToken1_, cuy::col, gather_cfg::cout, ztail::d, alignCSCRings::e, edm::Event::getByToken(), MixCollection< T >::MixItr::getTrigger(), HepMCToken_, histHepMCProduct_bunch_, histPCaloHit_bunch_, histPSimHit_bunchPileups_MuonCSCHits_, histPSimHit_bunchPileups_TrackerHitsTECHighTof_, histPSimHit_bunchSignal_MuonCSCHits_, histPSimHit_bunchSignal_TrackerHitsTECHighTof_, histTrack_bunchPileups_, histTrack_bunchSignal_, histVertex_bunch_, MuonToken_, outputFile, edm::Handle< T >::product(), SimTrackToken_, SimVertexToken_, AlCaHLTBitMon_QueryRunRegistry::string, tofhist_, tofhist_sig_, and TrackerToken0_.

 { 
   // test SimTracks
   //----------------------
   edm::Handle<CrossingFrame<SimTrack> > cf_simtrack;
   bool gotTracks = iEvent.getByToken(SimTrackToken_,cf_simtrack);
   if (!gotTracks)  outputFile<<" Could not read SimTracks!!!!" 
            << " Please, check if the object SimTracks has been declared in the"
            << " MixingModule configuration file."<<std::endl;
   
   if (gotTracks) {
     outputFile<<"\n=================== Starting SimTrack access ==================="<<std::endl;
 
     std::auto_ptr<MixCollection<SimTrack> > col1(new MixCollection<SimTrack>(cf_simtrack.product()));
     MixCollection<SimTrack>::iterator cfi1;
     int count1=0;
     std::cout <<" \nWe got "<<col1->sizeSignal()<<" signal tracks and "<<col1->sizePileup()<<" pileup tracks, total: "<<col1->size()<<std::endl;
     for (cfi1=col1->begin(); cfi1!=col1->end();cfi1++) {
       //std::cout << " BUNCH cfi1.bunch() = " << cfi1.bunch() << std::endl;
       if (cfi1.getTrigger()==0){
       histTrack_bunchPileups_->Fill(cfi1.bunch());
       }
       
       if (cfi1.getTrigger()==1){
       histTrack_bunchSignal_->Fill(cfi1.bunch());
       }
 
       
       int a = count1%4;
       if (a==3){
       outputFile<<" SimTrack "<<count1<<" has genpart index  "<<cfi1->genpartIndex()<<" vertex Index "<<cfi1->vertIndex() <<" bunchcr "<<cfi1.bunch()<<" trigger "<<cfi1.getTrigger()<<", from EncodedEventId: "<<cfi1->eventId().bunchCrossing() <<" "<<cfi1->eventId().event() <<std::endl;
       }
       count1++; 
     }
   }
 
   
   // test SimVertices
   //---------------------
   edm::Handle<CrossingFrame<SimVertex> > cf_simvtx;
   bool gotSimVertex = iEvent.getByToken(SimVertexToken_,cf_simvtx);
   if (!gotSimVertex) outputFile<<" Could not read Simvertices !!!!"<<std::endl;
   else {
     outputFile<<"\n=================== Starting SimVertex access ==================="<<std::endl;
     std::auto_ptr<MixCollection<SimVertex> > col2(new MixCollection<SimVertex>(cf_simvtx.product()));
     MixCollection<SimVertex>::iterator cfi2;
     int count2=0;
     outputFile <<" \nWe got "<<col2->sizeSignal()<<" signal vertices and "<<col2->sizePileup()<<" pileup vertices, total: "<<col2->size()<<std::endl;
     for (cfi2=col2->begin(); cfi2!=col2->end();cfi2++) {
       histVertex_bunch_->Fill(cfi2.bunch());
       int b = count2%4;
       if (count2 == 0 || b==3){
       //outputFile<<" SimVertex "<<count2<<" has parent index  "<<cfi2->parentIndex()<<" bunchcr "<<cfi2.bunch()<<" trigger "<<cfi2.getTrigger()<<", from EncodedEventId: "<<cfi2->eventId().bunchCrossing() <<" "<<cfi2->eventId().event() <<std::endl;
       }
       //SimVertex myvtx=(*cfi2);
       //outputFile<<"Same with op*: "<<count2<<" has parent index  "<<myvtx.parentIndex()<<" bunchcr "<<cfi2.bunch()<<" trigger "<<cfi2.getTrigger()<<", from EncodedEventId: "<<myvtx.eventId().bunchCrossing() <<" "<<myvtx.eventId().event() <<std::endl;
       count2++; 
     }
   }
   
   
   //test HepMCProducts
   //------------------------------------
   edm::Handle<CrossingFrame<edm::HepMCProduct> > cf_hepmc;
   bool gotHepMCP = iEvent.getByToken(HepMCToken_,cf_hepmc);
   if (!gotHepMCP) std::cout<<" Could not read HepMCProducts!!!!"<<std::endl;
   else {
     outputFile<<"\n=================== Starting HepMCProduct access ==================="<<std::endl;
     std::auto_ptr<MixCollection<edm::HepMCProduct> > colhepmc(new MixCollection<edm::HepMCProduct>(cf_hepmc.product()));
     MixCollection<edm::HepMCProduct>::iterator cfihepmc;
     
     int count3=0;
     outputFile <<" \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++) {
       histHepMCProduct_bunch_->Fill(cfihepmc.bunch());
       int c = count3%4;
       if (count3==0 || c==3){
       //outputFile<<" edm::HepMCProduct "<<count3<<" 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(count3);
       //outputFile<<"same with getObject:hepmc product   "<<count3<<" has event number "<<myprod.GetEvent()->event_number()<<", "<<myprod.GetEvent()->particles_size()<<" particles and "<<myprod.GetEvent()->vertices_size()<<" vertices"<<std::endl;
       count3++;
     }
   }
 
 
   // test CaloHits
   //--------------------------  
   const std::string subdetcalo("g4SimHitsEcalHitsEB");
   edm::Handle<CrossingFrame<PCaloHit> > cf_calo;
   bool gotPCaloHit = iEvent.getByToken(CaloToken1_,cf_calo);
   if (!gotPCaloHit) outputFile<<" Could not read CaloHits with label "<<subdetcalo<<"!!!!"<<std::endl;
   else {
     outputFile<<"\n\n=================== Starting CaloHit access, subdet "<<subdetcalo<<"  ==================="<<std::endl;
     std::auto_ptr<MixCollection<PCaloHit> > colcalo(new MixCollection<PCaloHit>(cf_calo.product()));
     //outputFile<<*(colcalo.get())<<std::endl;
     MixCollection<PCaloHit>::iterator cficalo;
     int count4=0;
     for (cficalo=colcalo->begin(); cficalo!=colcalo->end();cficalo++) {
       histPCaloHit_bunch_->Fill(cficalo.bunch());
       int d = count4%4;
       if (count4==0 || d==3){
       //outputFile<<" CaloHit "<<count4<<" has tof "<<cficalo->time()<<" trackid "<<cficalo->geantTrackId() <<" bunchcr "<<cficalo.bunch()<<" trigger "<<cficalo.getTrigger()<<", from EncodedEventId: "<<cficalo->eventId().bunchCrossing()<<" " <<cficalo->eventId().event() <<std::endl;
       }
       count4++;
     }
   }
 
   
   // test PSimHit for one particular subdet
   //---------------------------------------
   
   
   const std::string subdet("g4SimHitsTrackerHitsTECHighTof");
   edm::Handle<CrossingFrame<PSimHit> > cf_simhit;
   bool gotPSimHit = iEvent.getByToken(TrackerToken0_,cf_simhit);
   if (!gotPSimHit) outputFile<<" Could not read SimHits with label "<<subdet<<"!!!!"<<std::endl;
   else {
     outputFile<<"\n\n=================== Starting SimHit access, subdet "<<subdet<<"  ==================="<<std::endl;
 
     std::auto_ptr<MixCollection<PSimHit> > col(new MixCollection<PSimHit>(cf_simhit.product()));
     //outputFile<<*(col.get())<<std::endl;
     MixCollection<PSimHit>::iterator cfi;
     int count5=0;
     for (cfi=col->begin(); cfi!=col->end();cfi++) {
       
       // Signal
       if (cfi.getTrigger()==1){ 
         histPSimHit_bunchSignal_TrackerHitsTECHighTof_->Fill(cfi.bunch());
         tofhist_sig_->Fill(cfi->timeOfFlight());
       }
       
       // Pileups
       if (cfi.getTrigger()==0){
        histPSimHit_bunchPileups_TrackerHitsTECHighTof_->Fill(cfi.bunch());
        std::cout << " cfi->timeOfFlight() = " << cfi->timeOfFlight() << std::endl;
        tofhist_->Fill(cfi->timeOfFlight()); 
       }
       
       int e = count5%4;
       if (e==3){
       outputFile<<" Hit "<<count5<<" 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;
       }
       count5++;
     }
   }
   
   const std::string subdet1("g4SimHitsMuonCSCHits");
   edm::Handle<CrossingFrame<PSimHit> > cf_simhit1;
   bool gotPSimHit1 = iEvent.getByToken(MuonToken_,cf_simhit1);
   if (!gotPSimHit1) outputFile<<" Could not read SimHits with label "<<subdet1<<"!!!!"<<std::endl;
   else {
     outputFile<<"\n\n=================== Starting SimHit access, subdet "<<subdet1<<"  ==================="<<std::endl;
 
     std::auto_ptr<MixCollection<PSimHit> > col(new MixCollection<PSimHit>(cf_simhit1.product()));
     //outputFile<<*(col.get())<<std::endl;
     MixCollection<PSimHit>::iterator cfi;
     int count5=0;
     for (cfi=col->begin(); cfi!=col->end();cfi++) {
       
       if (cfi.getTrigger()==1) histPSimHit_bunchSignal_MuonCSCHits_->Fill(cfi.bunch());
       
       if (cfi.getTrigger()==0) histPSimHit_bunchPileups_MuonCSCHits_->Fill(cfi.bunch());
       
       int e = count5%4;
       if (e==3){
       outputFile<<" Hit "<<count5<<" 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;
       }
       count5++;
     }
   }
   
 
 
 }

void TestMixedSource::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 296 of file TestMixedSource.cc.

References gather_cfg::cout, and outputFile.

 { 
    outputFile.open("test.log");
    if (!outputFile.is_open())
    {
      std::cout << "Unable to open file!" << std::endl;
    } 
    
 }

void TestMixedSource::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 308 of file TestMixedSource.cc.

References gather_cfg::cout, fileName_, histFile_, histHepMCProduct_bunch_, histPCaloHit_bunch_, histPSimHit_bunchPileups_MuonCSCHits_, histPSimHit_bunchPileups_TrackerHitsTECHighTof_, histPSimHit_bunchSignal_MuonCSCHits_, histPSimHit_bunchSignal_TrackerHitsTECHighTof_, histTrack_bunchPileups_, histTrack_bunchSignal_, histVertex_bunch_, outputFile, tree::t, tofhist_, and tofhist_sig_.

                         {
   
   outputFile.close();
   
   char t[30];
   sprintf(t,"%s",fileName_.c_str());
   std::cout << " fileName = " << t << std::endl; 
   histFile_=new TFile(t,"RECREATE");
   histTrack_bunchPileups_->Write();
   histTrack_bunchSignal_->Write();
   histVertex_bunch_->Write();
   histPCaloHit_bunch_->Write();
   histPSimHit_bunchPileups_TrackerHitsTECHighTof_->Write();
   histPSimHit_bunchSignal_TrackerHitsTECHighTof_->Write();
   tofhist_sig_->Write();
   tofhist_->Write();
   histPSimHit_bunchSignal_MuonCSCHits_->Write();
   histPSimHit_bunchPileups_MuonCSCHits_->Write();
   histHepMCProduct_bunch_->Write();
   histFile_->Write();
   histFile_->Close();
 
 }