|
|
|
#include <EcalPreshowerRecHitsValidation.h>
Definition at line 35 of file EcalPreshowerRecHitsValidation.h.
| EcalPreshowerRecHitsValidation::EcalPreshowerRecHitsValidation | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 19 of file EcalPreshowerRecHitsValidation.cc.
References dbe_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), timingPdfMaker::histo, and cmsCodeRules::cppFunctionSkipper::operator.
{
// ----------------------
EEuncalibrechitCollection_ = ps.getParameter<edm::InputTag>("EEuncalibrechitCollection");
EErechitCollection_ = ps.getParameter<edm::InputTag>("EErechitCollection");
ESrechitCollection_ = ps.getParameter<edm::InputTag>("ESrechitCollection");
// ----------------------
// verbosity switch
verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
// ----------------------
// get hold of back-end interface
dbe_ = 0;
dbe_ = Service<DQMStore>().operator->();
if ( dbe_ ) {
if ( verbose_ ) {
dbe_->setVerbose(1);
} else {
dbe_->setVerbose(0);
}
}
if ( dbe_ ) {
if ( verbose_ ) dbe_->showDirStructure();
}
// ----------------------
meESRecHitsEnergy_ = 0; // total energy
meESRecHitsEnergy_zp1st_ = 0;
meESRecHitsEnergy_zp2nd_ = 0;
meESRecHitsEnergy_zm1st_ = 0;
meESRecHitsEnergy_zm2nd_ = 0;
meESRecHitsMultip_ = 0; // total multiplicity
meESRecHitsMultip_zp1st_ = 0;
meESRecHitsMultip_zp2nd_ = 0;
meESRecHitsMultip_zm1st_ = 0;
meESRecHitsMultip_zm2nd_ = 0;
meESEERecHitsEnergy_zp_ = 0; // versus EE energy
meESEERecHitsEnergy_zm_ = 0;
for (int kk=0; kk<32; kk++)
{
meESRecHitsStripOccupancy_zp1st_[kk] = 0;
meESRecHitsStripOccupancy_zm1st_[kk] = 0;
meESRecHitsStripOccupancy_zp2nd_[kk] = 0;
meESRecHitsStripOccupancy_zm2nd_[kk] = 0;
}
// ----------------------
Char_t histo[200];
if ( dbe_ )
{
dbe_->setCurrentFolder("EcalRecHitsV/EcalPreshowerRecHitsTask");
sprintf (histo, "ES Energy" );
meESRecHitsEnergy_ = dbe_->book1D(histo, histo, 210, -0.0005, 0.01);
sprintf (histo, "ES Energy Plane1 Side+" );
meESRecHitsEnergy_zp1st_ = dbe_->book1D(histo, histo, 210, -0.0005, 0.01);
sprintf (histo, "ES Energy Plane2 Side+");
meESRecHitsEnergy_zp2nd_ = dbe_->book1D(histo, histo, 210, -0.0005, 0.01);
sprintf (histo, "ES Energy Plane1 Side-");
meESRecHitsEnergy_zm1st_ = dbe_->book1D(histo, histo, 210, -0.0005, 0.01);
sprintf (histo, "ES Energy Plane2 Side-");
meESRecHitsEnergy_zm2nd_ = dbe_->book1D(histo, histo, 210, -0.0005, 0.01);
sprintf (histo, "ES Multiplicity" );
meESRecHitsMultip_ = dbe_->book1D(histo, histo, 100, 0., 700.);
sprintf (histo, "ES Multiplicity Plane1 Side+");
meESRecHitsMultip_zp1st_ = dbe_->book1D(histo, histo, 100, 0., 700.);
sprintf (histo, "ES Multiplicity Plane2 Side+");
meESRecHitsMultip_zp2nd_ = dbe_->book1D(histo, histo, 100, 0., 700.);
sprintf (histo, "ES Multiplicity Plane1 Side-");
meESRecHitsMultip_zm1st_ = dbe_->book1D(histo, histo, 100, 0., 700.);
sprintf (histo, "ES Multiplicity Plane2 Side-");
meESRecHitsMultip_zm2nd_ = dbe_->book1D(histo, histo, 100, 0., 700.);
sprintf (histo, "Preshower EE vs ES energy Side+");
meESEERecHitsEnergy_zp_ = dbe_->book2D(histo, histo, 100, 0., 0.2, 100, 0., 150.);
sprintf (histo, "Preshower EE vs ES energy Side-");
meESEERecHitsEnergy_zm_ = dbe_->book2D(histo, histo, 100, 0., 0.2, 100, 0., 150.);
for (int kk=0; kk<32; kk++)
{
sprintf(histo, "ES Occupancy Plane1 Side+ Strip%02d", kk+1);
meESRecHitsStripOccupancy_zp1st_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
sprintf(histo, "ES Occupancy Plane2 Side+ Strip%02d", kk+1);
meESRecHitsStripOccupancy_zp2nd_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
sprintf(histo, "ES Occupancy Plane1 Side- Strip%02d", kk+1);
meESRecHitsStripOccupancy_zm1st_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
sprintf(histo, "ES Occupancy Plane2 Side- Strip%02d", kk+1);
meESRecHitsStripOccupancy_zm2nd_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
}
}
}
| EcalPreshowerRecHitsValidation::~EcalPreshowerRecHitsValidation | ( | ) |
| void EcalPreshowerRecHitsValidation::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Analyze.
Implements edm::EDAnalyzer.
Definition at line 142 of file EcalPreshowerRecHitsValidation.cc.
References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), HcalObjRepresent::Fill(), edm::SortedCollection< T, SORT >::find(), edm::Event::getByLabel(), edm::HandleBase::isValid(), ESDetId::plane(), edm::Handle< T >::product(), ESDetId::six(), ESDetId::siy(), ESDetId::strip(), strip(), ESDetId::zside(), and EEDetId::zside().
{
const ESRecHitCollection *ESRecHit = 0;
Handle<ESRecHitCollection> EcalRecHitES;
e.getByLabel( ESrechitCollection_, EcalRecHitES);
if (EcalRecHitES.isValid()) {
ESRecHit = EcalRecHitES.product ();
} else {
return;
}
bool skipEE = false;
const EERecHitCollection *EERecHit = 0;
Handle<EERecHitCollection> EcalRecHitEE;
e.getByLabel( EErechitCollection_, EcalRecHitEE);
if (EcalRecHitEE.isValid()){
EERecHit = EcalRecHitEE.product ();
} else {
skipEE = true;
}
const EEUncalibratedRecHitCollection *EEUncalibRecHit = 0;
Handle< EEUncalibratedRecHitCollection > EcalUncalibRecHitEE;
e.getByLabel( EEuncalibrechitCollection_, EcalUncalibRecHitEE);
if (EcalUncalibRecHitEE.isValid()) {
EEUncalibRecHit = EcalUncalibRecHitEE.product() ;
} else {
skipEE = true;
}
// ----------------------
// loop over RecHits
// multiplicities
int mult_tot = 0;
int mult_zp1st = 0;
int mult_zp2nd = 0;
int mult_zm1st = 0;
int mult_zm2nd = 0;
// energies
float ene_zp1st = 0.;
float ene_zp2nd = 0.;
float ene_zm1st = 0.;
float ene_zm2nd = 0.;
// ES
for (ESRecHitCollection::const_iterator recHit = ESRecHit->begin(); recHit != ESRecHit->end() ; ++recHit)
{
ESDetId ESid = ESDetId(recHit->id());
int zside = ESid.zside();
int plane = ESid.plane();
int six = ESid.six();
int siy = ESid.siy();
int strip = ESid.strip();
// global
mult_tot++;
if (meESRecHitsEnergy_) meESRecHitsEnergy_ ->Fill(recHit->energy());
// side +, plane 1
if ( (zside == +1) && (plane == 1) )
{
mult_zp1st++;
ene_zp1st += recHit->energy();
if ( meESRecHitsEnergy_zp1st_ ) { meESRecHitsEnergy_zp1st_ -> Fill(recHit->energy()); }
if ( meESRecHitsStripOccupancy_zp1st_[strip-1] ){ meESRecHitsStripOccupancy_zp1st_[strip-1]-> Fill( six, siy ); }
}
// side +, plane 2
if ( (zside == +1) && (plane == 2) )
{
mult_zp2nd++;
ene_zp2nd += recHit->energy();
if ( meESRecHitsEnergy_zp2nd_ ) { meESRecHitsEnergy_zp2nd_ -> Fill(recHit->energy()); }
if ( meESRecHitsStripOccupancy_zp2nd_[strip-1] ){ meESRecHitsStripOccupancy_zp2nd_[strip-1]-> Fill( six, siy ); }
}
// side -, plane 1
if ( (zside == -1) && (plane == 1) )
{
mult_zm1st++;
ene_zm1st += recHit->energy();
if ( meESRecHitsEnergy_zm1st_ ) { meESRecHitsEnergy_zm1st_ -> Fill(recHit->energy()); }
if ( meESRecHitsStripOccupancy_zm1st_[strip-1] ){ meESRecHitsStripOccupancy_zm1st_[strip-1]-> Fill( six, siy ); }
}
// side +, plane 2
if ( (zside == -1) && (plane == 2) )
{
mult_zm2nd ++;
ene_zm2nd += recHit->energy();
if ( meESRecHitsEnergy_zm2nd_ ) { meESRecHitsEnergy_zm2nd_ -> Fill(recHit->energy()); }
if ( meESRecHitsStripOccupancy_zm2nd_[strip-1] ){ meESRecHitsStripOccupancy_zm2nd_[strip-1]-> Fill( six, siy ); }
}
} // loop over the ES RecHitCollection
// EE
double zpEE = 0.;
double zmEE = 0.;
if ( ! skipEE ) {
for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EEUncalibRecHit->begin(); uncalibRecHit != EEUncalibRecHit->end() ; ++uncalibRecHit)
{
EEDetId EEid = EEDetId(uncalibRecHit->id());
int mySide = EEid.zside();
// Find corresponding recHit
EcalRecHitCollection::const_iterator myRecHit = EERecHit->find(EEid);
if (myRecHit != EERecHit->end() )
{
if (mySide > 0) { zpEE = zpEE + myRecHit->energy(); }
if (mySide < 0) { zmEE = zmEE + myRecHit->energy(); }
}
}
}
// filling histos
if (meESRecHitsMultip_) { meESRecHitsMultip_ -> Fill(mult_tot); }
if (meESRecHitsMultip_zp1st_ ) { meESRecHitsMultip_zp1st_ -> Fill(mult_zp1st); }
if (meESRecHitsMultip_zp2nd_ ) { meESRecHitsMultip_zp2nd_ -> Fill(mult_zp2nd); }
if (meESRecHitsMultip_zm1st_ ) { meESRecHitsMultip_zm1st_ -> Fill(mult_zm1st); }
if (meESRecHitsMultip_zm2nd_ ) { meESRecHitsMultip_zm2nd_ -> Fill(mult_zm2nd); }
if (meESEERecHitsEnergy_zp_) { meESEERecHitsEnergy_zp_ -> Fill( (ene_zp1st + 0.7*ene_zp2nd)/0.09, zpEE ); }
if (meESEERecHitsEnergy_zm_) { meESEERecHitsEnergy_zm_ -> Fill( (ene_zm1st + 0.7*ene_zm2nd)/0.09, zmEE ); }
}
| void EcalPreshowerRecHitsValidation::beginJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 134 of file EcalPreshowerRecHitsValidation.cc.
{
}
| void EcalPreshowerRecHitsValidation::endJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 138 of file EcalPreshowerRecHitsValidation.cc.
{
}
DQMStore* EcalPreshowerRecHitsValidation::dbe_ [private] |
Definition at line 60 of file EcalPreshowerRecHitsValidation.h.
Definition at line 63 of file EcalPreshowerRecHitsValidation.h.
Definition at line 62 of file EcalPreshowerRecHitsValidation.h.
Definition at line 64 of file EcalPreshowerRecHitsValidation.h.
Definition at line 77 of file EcalPreshowerRecHitsValidation.h.
Definition at line 76 of file EcalPreshowerRecHitsValidation.h.
Definition at line 66 of file EcalPreshowerRecHitsValidation.h.
Definition at line 69 of file EcalPreshowerRecHitsValidation.h.
Definition at line 70 of file EcalPreshowerRecHitsValidation.h.
Definition at line 67 of file EcalPreshowerRecHitsValidation.h.
Definition at line 68 of file EcalPreshowerRecHitsValidation.h.
Definition at line 71 of file EcalPreshowerRecHitsValidation.h.
Definition at line 74 of file EcalPreshowerRecHitsValidation.h.
Definition at line 75 of file EcalPreshowerRecHitsValidation.h.
Definition at line 72 of file EcalPreshowerRecHitsValidation.h.
Definition at line 73 of file EcalPreshowerRecHitsValidation.h.
Definition at line 79 of file EcalPreshowerRecHitsValidation.h.
Definition at line 81 of file EcalPreshowerRecHitsValidation.h.
Definition at line 78 of file EcalPreshowerRecHitsValidation.h.
Definition at line 80 of file EcalPreshowerRecHitsValidation.h.
bool EcalPreshowerRecHitsValidation::verbose_ [private] |
Definition at line 58 of file EcalPreshowerRecHitsValidation.h.
1.7.3