|
|
|
#include <TestSuite.h>
Public Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
| virtual void | beginJob () |
| virtual void | endJob () |
| TestSuite (const edm::ParameterSet &) | |
| ~TestSuite () | |
Private Attributes | |
| int | bunchcr_ |
| DQMStore * | dbe_ |
| std::string | filename_ |
| int | maxbunch_ |
| int | minbunch_ |
Description: test suite for Mixing Module
Description: global physics test for Mixing Module
Definition at line 37 of file TestSuite.h.
| TestSuite::TestSuite | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 41 of file TestSuite.cc.
References bunchcr_, gather_cfg::cout, and filename_.
: filename_(iConfig.getParameter<std::string>("fileName")), bunchcr_(iConfig.getParameter<int>("BunchNr")), minbunch_(iConfig.getParameter<int>("minBunch")),maxbunch_(iConfig.getParameter<int>("maxBunch")), dbe_(0) { std::cout << "Constructed testSuite , bunchcr "<<bunchcr_<<" filename: "<<filename_<<std::endl; }
| TestSuite::~TestSuite | ( | ) |
Definition at line 48 of file TestSuite.cc.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
}
| void TestSuite::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 71 of file TestSuite.cc.
References DQMStore::book1D(), bunchcr_, dbe_, edm::Event::getByLabel(), edm::Event::getByType(), maxbunch_, and minbunch_.
{
using namespace edm;
// Get input
edm::Handle<CrossingFrame<SimTrack> > cf_track;
edm::Handle<CrossingFrame<SimVertex> > cf_vertex;
edm::Handle<CrossingFrame<PSimHit> > cf_simhit;
edm::Handle<CrossingFrame<PCaloHit> > cf_calohitEcal;
edm::Handle<CrossingFrame<PCaloHit> > cf_calohitHcal;
std::string subdetTracker("g4SimHitsTrackerHitsTECLowTof");
std::string ecalsubdet("g4SimHitsEcalHitsEB");
std::string hcalsubdet("g4SimHitsHcalHits");
iEvent.getByType(cf_track);
iEvent.getByType(cf_vertex);
iEvent.getByLabel("mix",subdetTracker,cf_simhit);
iEvent.getByLabel("mix",ecalsubdet,cf_calohitEcal);
iEvent.getByLabel("mix",hcalsubdet,cf_calohitHcal);
// use MixCollection and its iterator
// Please note that bunch() and getTrigger() are methods of the iterator itself
// while operator-> points to the templated objects!!!!
//track histo
char histotracks[30],sighistotracks[30],histotracksindsig[30],histotracksind[30];
sprintf(histotracks,"Tracks_bcr_%d",bunchcr_);
sprintf(sighistotracks,"SignalTracks_bcr_%d",bunchcr_);
sprintf(histotracksind,"VtxPointers_%d",bunchcr_);
sprintf(histotracksindsig,"VtxPointers_signal_%d",bunchcr_);
MonitorElement * trhist = dbe_->book1D(histotracks,"Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
MonitorElement * trhistsig = dbe_->book1D(sighistotracks,"Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
MonitorElement * trindhist = dbe_->book1D(histotracksind,"Track to Vertex indices",100,0,500);
MonitorElement * trindhistsig = dbe_->book1D(histotracksindsig,"Signal Track to Vertex indices",100,0,500);
std::auto_ptr<MixCollection<SimTrack> > col1(new MixCollection<SimTrack>(cf_track.product()));
MixCollection<SimTrack>::iterator cfi1;
for (cfi1=col1->begin(); cfi1!=col1->end();cfi1++) {
if (cfi1.getTrigger()==0) {
trhist->Fill(cfi1.bunch());
trindhist->Fill(cfi1->vertIndex());
} else {
trindhistsig->Fill(cfi1->vertIndex());
trhistsig->Fill(cfi1.bunch());
}
}
//vertex histo
char histovertices[30], sighistovertices[30],histovertexindices[30],histovertexindicessig[30];
sprintf(histovertices,"Vertices_bcr_%d",bunchcr_);
sprintf(sighistovertices,"SignalVertices_bcr_%d",bunchcr_);
sprintf(histovertexindices,"TrackPointers_%d",bunchcr_);
sprintf(histovertexindicessig,"TrackPointers_signal_%d",bunchcr_);
MonitorElement * vtxhist = dbe_->book1D(histovertices,"Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
MonitorElement * vtxhistsig = dbe_->book1D(sighistovertices,"Bunchcrossings",maxbunch_-minbunch_+1,minbunch_,maxbunch_+1);
MonitorElement * vtxindhist = dbe_->book1D(histovertexindices,"Vertex to Track Indices",100,0,300);
MonitorElement * vtxindhistsig = dbe_->book1D(histovertexindicessig,"Signal Vertex to Track Indices",100,0,300);
std::auto_ptr<MixCollection<SimVertex> > col2(new MixCollection<SimVertex>(cf_vertex.product()));
MixCollection<SimVertex>::iterator cfi2;
for (cfi2=col2->begin(); cfi2!=col2->end();cfi2++) {
if (cfi2.getTrigger()==0) {
vtxhist->Fill(cfi2.bunch());
if (!cfi2->noParent()) vtxindhist->Fill(cfi2->parentIndex());
} else {
vtxhistsig->Fill(cfi2.bunch());
if (!cfi2->noParent()) vtxindhistsig->Fill(cfi2->parentIndex());
}
}
//tracker
int bsp=cf_simhit->getBunchSpace();
char tof[30];
sprintf(tof,"TrackerHit_Tof_bcr_%d",bunchcr_);
MonitorElement * tofhist = dbe_->book1D(tof,"TrackerHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
sprintf(tof,"SignalTrackerHit_Tof_bcr_%d",bunchcr_);
MonitorElement * tofhist_sig = dbe_->book1D(tof,"TrackerHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
std::auto_ptr<MixCollection<PSimHit> > colsh(new MixCollection<PSimHit>(cf_simhit.product()));
MixCollection<PSimHit>::iterator cfish;
for (cfish=colsh->begin(); cfish!=colsh->end();cfish++) {
if (cfish.getTrigger()) {
tofhist_sig->Fill(cfish->timeOfFlight());
}
else {
tofhist->Fill(cfish->timeOfFlight());
}
}
//Ecal
sprintf(tof,"EcalEBHit_Tof_bcr_%d",bunchcr_);
MonitorElement * tofecalhist = dbe_->book1D(tof,"EcalEBHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
sprintf(tof,"SignalEcalEBHit_Tof_bcr_%d",bunchcr_);
MonitorElement * tofecalhist_sig = dbe_->book1D(tof,"EcalEBHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
// std::string ecalsubdet("EcalHitsEB");
std::auto_ptr<MixCollection<PCaloHit> > colecal(new MixCollection<PCaloHit>(cf_calohitEcal.product()));
MixCollection<PCaloHit>::iterator cfiecal;
for (cfiecal=colecal->begin(); cfiecal!=colecal->end();cfiecal++) {
if (cfiecal.getTrigger()) tofecalhist_sig->Fill(cfiecal->time());
else tofecalhist->Fill(cfiecal->time());
}
// Hcal
sprintf(tof,"HcalHit_Tof_bcr_%d",bunchcr_);
MonitorElement * tofhcalhist = dbe_->book1D(tof,"HcalHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
sprintf(tof,"SignalHcalHit_Tof_bcr_%d",bunchcr_);
MonitorElement * tofhcalhist_sig = dbe_->book1D(tof,"HcalHit_ToF",100,float(bsp*minbunch_),float(bsp*maxbunch_)+50.);
// std::string hcalsubdet("HcalHits");
std::auto_ptr<MixCollection<PCaloHit> > colhcal(new MixCollection<PCaloHit>(cf_calohitHcal.product()));
MixCollection<PCaloHit>::iterator cfihcal;
for (cfihcal=colhcal->begin(); cfihcal!=colhcal->end();cfihcal++) {
if (cfihcal.getTrigger()) tofhcalhist_sig->Fill(cfihcal->time());
else tofhcalhist->Fill(cfihcal->time());
}
}
| void TestSuite::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 55 of file TestSuite.cc.
References dbe_, cmsCodeRules::cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), and DQMStore::showDirStructure().
{
// get hold of back-end interface
dbe_ = Service<DQMStore>().operator->();
dbe_->showDirStructure();
dbe_->setCurrentFolder("MixingV/Mixing");
}
| void TestSuite::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 63 of file TestSuite.cc.
References dbe_, filename_, and DQMStore::save().
int TestSuite::bunchcr_ [private] |
Definition at line 49 of file TestSuite.h.
Referenced by analyze(), and TestSuite().
DQMStore* TestSuite::dbe_ [private] |
Definition at line 52 of file TestSuite.h.
Referenced by analyze(), beginJob(), and endJob().
std::string TestSuite::filename_ [private] |
Definition at line 48 of file TestSuite.h.
Referenced by endJob(), and TestSuite().
int TestSuite::maxbunch_ [private] |
Definition at line 51 of file TestSuite.h.
Referenced by analyze().
int TestSuite::minbunch_ [private] |
Definition at line 50 of file TestSuite.h.
Referenced by analyze().
1.7.3