#include <GlobalTest.h>

Inheritance diagram for GlobalTest:

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

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	GlobalTest (const edm::ParameterSet &)

	~GlobalTest () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Attributes
MonitorElement *	caloEnergyEBH_ [nMaxH]

MonitorElement *	caloEnergyEEH_ [nMaxH]

edm::EDGetTokenT< CrossingFrame< SimTrack > >	cfTrackToken_

edm::EDGetTokenT< CrossingFrame< SimTrack > >	cfVertexToken_

std::string	filename_

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	g4SimHits_EB_Token_

edm::EDGetTokenT< CrossingFrame< PCaloHit > >	g4SimHits_EE_Token_

TFile *	histfile_

char *	labels [nrHistos]

int	maxbunch_

int	minbunch_

MonitorElement *	nrPileupsH_ [nMaxH]

MonitorElement *	nrTracksH_ [nMaxH]

MonitorElement *	nrVerticesH_ [nMaxH]

MonitorElement *	trackPartIdH_ [nMaxH]

Static Private Attributes
static const int	nMaxH = 10

static const int	nrHistos = 6

Detailed Description

Description: test suite for Mixing Module

Definition at line 45 of file GlobalTest.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 39 of file GlobalTest.cc.

References gather_cfg::cout, filename_, g4SimHits_EB_Token_, g4SimHits_EE_Token_, maxbunch_, minbunch_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : filename_(iConfig.getParameter<std::string>("fileName")),
       minbunch_(iConfig.getParameter<int>("minBunch")),
       maxbunch_(iConfig.getParameter<int>("maxBunch")),
       cfTrackToken_(consumes<CrossingFrame<SimTrack>>(iConfig.getParameter<edm::InputTag>("cfTrackTag"))),
       cfVertexToken_(consumes<CrossingFrame<SimTrack>>(iConfig.getParameter<edm::InputTag>("cfVertexTag"))) {
   std::string ecalsubdetb("");
   std::string ecalsubdete("g4SimHitsEcalHitsEE");
   g4SimHits_EB_Token_ = consumes<CrossingFrame<PCaloHit>>(edm::InputTag("mix", "g4SimHitsEcalHitsEB"));
   g4SimHits_EE_Token_ = consumes<CrossingFrame<PCaloHit>>(edm::InputTag("mix", "g4SimHitsEcalHitsEE"));
 
   std::cout << "Constructed GlobalTest, filename: " << filename_ << " minbunch: " << minbunch_
             << ", maxbunch: " << maxbunch_ << std::endl;
 }

GlobalTest::~GlobalTest ( )

override

Definition at line 54 of file GlobalTest.cc.

References mps_fire::i, and labels.

                         {
   for (int i = 0; i < 6; i++)
     delete[] labels[i];
 }

Member Function Documentation

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

override

Definition at line 114 of file GlobalTest.cc.

References MixCollection< T >::MixItr::bunch(), caloEnergyEBH_, caloEnergyEEH_, cfTrackToken_, cfVertexToken_, MonitorElement::Fill(), g4SimHits_EB_Token_, g4SimHits_EE_Token_, edm::Event::getByToken(), mps_fire::i, maxbunch_, minbunch_, nrTracksH_, nrVerticesH_, edm::Handle< T >::product(), AlCaHLTBitMon_QueryRunRegistry::string, and trackPartIdH_.

                                                                             {
   using namespace edm;
   using namespace std;
 
   // Get input
   edm::Handle<CrossingFrame<SimTrack>> cf_track;
   edm::Handle<CrossingFrame<SimTrack>> cf_vertex;
   edm::Handle<CrossingFrame<PCaloHit>> cf_calohitE;
   edm::Handle<CrossingFrame<PCaloHit>> cf_calohitB;
   std::string ecalsubdetb("g4SimHitsEcalHitsEB");
   std::string ecalsubdete("g4SimHitsEcalHitsEE");
   iEvent.getByToken(cfTrackToken_, cf_track);
   iEvent.getByToken(cfVertexToken_, cf_vertex);
   iEvent.getByToken(g4SimHits_EB_Token_, cf_calohitB);
   iEvent.getByToken(g4SimHits_EE_Token_, cf_calohitE);
 
   // number of events/bcr ??
 
   // number of tracks
   for (int i = minbunch_; i <= maxbunch_; ++i) {
     nrTracksH_[i - minbunch_]->Fill(cf_track->getNrPileups(i));
   }
 
   // number of vertices
   for (int i = minbunch_; i <= maxbunch_; ++i) {
     nrVerticesH_[i - minbunch_]->Fill(cf_vertex->getNrPileups(i));
   }
 
   // part id for each track
   std::unique_ptr<MixCollection<SimTrack>> coltr(new MixCollection<SimTrack>(cf_track.product()));
   MixCollection<SimTrack>::iterator cfitr;
   for (cfitr = coltr->begin(); cfitr != coltr->end(); cfitr++) {
     trackPartIdH_[cfitr.bunch() - minbunch_]->Fill(cfitr->type());
   }
 
   // energy sum
   double sumE[10] = {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.};
   std::unique_ptr<MixCollection<PCaloHit>> colecalb(new MixCollection<PCaloHit>(cf_calohitB.product()));
   MixCollection<PCaloHit>::iterator cfiecalb;
   for (cfiecalb = colecalb->begin(); cfiecalb != colecalb->end(); cfiecalb++) {
     sumE[cfiecalb.bunch() - minbunch_] += cfiecalb->energy();
     //      if (cfiecal.getTrigger())    tofecalhist_sig->Fill(cfiecal->time());
     //      else    tofecalhist->Fill(cfiecal->time());
   }
   for (int i = minbunch_; i <= maxbunch_; ++i) {
     caloEnergyEBH_[i - minbunch_]->Fill(sumE[i - minbunch_]);
   }
   double sumEE[10] = {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.};
   std::unique_ptr<MixCollection<PCaloHit>> colecale(new MixCollection<PCaloHit>(cf_calohitE.product()));
   MixCollection<PCaloHit>::iterator cfiecale;
   for (cfiecale = colecale->begin(); cfiecale != colecale->end(); cfiecale++) {
     sumEE[cfiecale.bunch() - minbunch_] += cfiecale->energy();
     //      if (cfiecal.getTrigger())    tofecalhist_sig->Fill(cfiecal->time());
     //      else    tofecalhist->Fill(cfiecal->time());
   }
   for (int i = minbunch_; i <= maxbunch_; ++i) {
     caloEnergyEEH_[i - minbunch_]->Fill(sumEE[i - minbunch_]);
   }
 }

void GlobalTest::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

override

Definition at line 59 of file GlobalTest.cc.

References DQMStore::IBooker::book1D(), caloEnergyEBH_, caloEnergyEEH_, mps_fire::i, cuy::ii, label, labels, maxbunch_, minbunch_, nrPileupsH_, nrTracksH_, nrVerticesH_, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and trackPartIdH_.

                                                                                                {
   using namespace std;
 
   ibooker.setCurrentFolder("MixingV/Mixing");
   // book histos
   std::string NrPileupEvts = "NrPileupEvts";
   size_t NrPileupEvtsSize = NrPileupEvts.size() + 1;
   std::string NrVertices = "NrVertices";
   size_t NrVerticesSize = NrVertices.size() + 1;
   std::string NrTracks = "NrTracks";
   size_t NrTracksSize = NrTracks.size() + 1;
   std::string TrackPartId = "TrackPartId";
   size_t TrackPartIdSize = TrackPartId.size() + 1;
   std::string CaloEnergyEB = "CaloEnergyEB";
   size_t CaloEnergyEBSize = CaloEnergyEB.size() + 1;
   std::string CaloEnergyEE = "CaloEnergyEE";
   size_t CaloEnergyEESize = CaloEnergyEE.size() + 1;
 
   labels[0] = new char[NrPileupEvtsSize];
   strncpy(labels[0], NrPileupEvts.c_str(), NrPileupEvtsSize);
   labels[1] = new char[NrVerticesSize];
   strncpy(labels[1], NrVertices.c_str(), NrVerticesSize);
   labels[2] = new char[NrTracksSize];
   strncpy(labels[2], NrTracks.c_str(), NrTracksSize);
   labels[3] = new char[TrackPartIdSize];
   strncpy(labels[3], TrackPartId.c_str(), TrackPartIdSize);
   labels[4] = new char[CaloEnergyEBSize];
   strncpy(labels[4], CaloEnergyEB.c_str(), CaloEnergyEBSize);
   labels[5] = new char[CaloEnergyEESize];
   strncpy(labels[5], CaloEnergyEE.c_str(), CaloEnergyEESize);
 
   // FIXME: test for max nr of histos
   for (int i = minbunch_; i <= maxbunch_; ++i) {
     int ii = i - minbunch_;
     char label[50];
     sprintf(label, "%s_%d", labels[0], i);
     nrPileupsH_[ii] = ibooker.book1D(label, label, 100, 0, 100);
     sprintf(label, "%s_%d", labels[1], i);
     nrVerticesH_[ii] = ibooker.book1D(label, label, 100, 0, 5000);
     sprintf(label, "%s_%d", labels[2], i);
     nrTracksH_[ii] = ibooker.book1D(label, label, 100, 0, 10000);
     sprintf(label, "%s_%d", labels[3], i);
     trackPartIdH_[ii] = ibooker.book1D(label, label, 100, 0, 100);
     sprintf(label, "%s_%d", labels[4], i);
     caloEnergyEBH_[ii] = ibooker.book1D(label, label, 100, 0., 1000.);
     sprintf(label, "%s_%d", labels[5], i);
     caloEnergyEEH_[ii] = ibooker.book1D(label, label, 100, 0., 1000.);
   }
 }