CMSSW Reference Manual

EcalPreshowerRecHitsValidation::EcalPreshowerRecHitsValidation ( const edm::ParameterSet & ps )

Constructor.

Definition at line 19 of file EcalPreshowerRecHitsValidation.cc.

References dbe_, EErechitCollection_, EEuncalibrechitCollection_, ESrechitCollection_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), histo, kk, meESEERecHitsEnergy_zm_, meESEERecHitsEnergy_zp_, meESRecHitsEnergy_, meESRecHitsEnergy_zm1st_, meESRecHitsEnergy_zm2nd_, meESRecHitsEnergy_zp1st_, meESRecHitsEnergy_zp2nd_, meESRecHitsMultip_, meESRecHitsMultip_zm1st_, meESRecHitsMultip_zm2nd_, meESRecHitsMultip_zp1st_, meESRecHitsMultip_zp2nd_, meESRecHitsStripOccupancy_zm1st_, meESRecHitsStripOccupancy_zm2nd_, meESRecHitsStripOccupancy_zp1st_, meESRecHitsStripOccupancy_zp2nd_, and verbose_.

00019                                                                                     {
00020 
00021 
00022   // ----------------------
00023   EEuncalibrechitCollection_ = ps.getParameter<edm::InputTag>("EEuncalibrechitCollection");
00024   EErechitCollection_        = ps.getParameter<edm::InputTag>("EErechitCollection");
00025   ESrechitCollection_        = ps.getParameter<edm::InputTag>("ESrechitCollection");
00026 
00027 
00028   // ---------------------- 
00029   // verbosity switch 
00030   verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
00031   
00032   // ----------------------                 
00033   // get hold of back-end interface 
00034   dbe_ = 0;
00035   dbe_ = Service<DQMStore>().operator->();                   
00036   if ( dbe_ ) {
00037     if ( verbose_ ) {
00038       dbe_->setVerbose(1);
00039     } else {
00040       dbe_->setVerbose(0);
00041     }
00042   }                                                                  
00043   if ( dbe_ ) {
00044     if ( verbose_ ) dbe_->showDirStructure();
00045   }
00046 
00047 
00048   // ----------------------  
00049   meESRecHitsEnergy_           = 0;                // total energy
00050   meESRecHitsEnergy_zp1st_     = 0;    
00051   meESRecHitsEnergy_zp2nd_     = 0;    
00052   meESRecHitsEnergy_zm1st_     = 0;    
00053   meESRecHitsEnergy_zm2nd_     = 0;    
00054   meESRecHitsMultip_           = 0;                // total multiplicity
00055   meESRecHitsMultip_zp1st_     = 0;    
00056   meESRecHitsMultip_zp2nd_     = 0;    
00057   meESRecHitsMultip_zm1st_     = 0;    
00058   meESRecHitsMultip_zm2nd_     = 0;    
00059   meESEERecHitsEnergy_zp_      = 0;                // versus EE energy 
00060   meESEERecHitsEnergy_zm_      = 0;      
00061 
00062   for (int kk=0; kk<32; kk++)
00063     { 
00064       meESRecHitsStripOccupancy_zp1st_[kk] = 0;   
00065       meESRecHitsStripOccupancy_zm1st_[kk] = 0;   
00066       meESRecHitsStripOccupancy_zp2nd_[kk] = 0;   
00067       meESRecHitsStripOccupancy_zm2nd_[kk] = 0;   
00068     }
00069 
00070 
00071   // ---------------------- 
00072   Char_t histo[200];
00073   if ( dbe_ ) 
00074     {
00075       dbe_->setCurrentFolder("EcalRecHitsV/EcalPreshowerRecHitsTask");
00076 
00077       sprintf (histo, "ES Energy" );
00078       meESRecHitsEnergy_ = dbe_->book1D(histo, histo, 110, -0.0005, 0.005);
00079       
00080       sprintf (histo, "ES Energy Plane1 Side+" );
00081       meESRecHitsEnergy_zp1st_ = dbe_->book1D(histo, histo, 110, -0.0005, 0.005);
00082 
00083       sprintf (histo, "ES Energy Plane2 Side+");
00084       meESRecHitsEnergy_zp2nd_ = dbe_->book1D(histo, histo, 110, -0.0005, 0.005);
00085      
00086       sprintf (histo, "ES Energy Plane1 Side-");
00087       meESRecHitsEnergy_zm1st_ = dbe_->book1D(histo, histo, 110, -0.0005, 0.005);
00088 
00089       sprintf (histo, "ES Energy Plane2 Side-");
00090       meESRecHitsEnergy_zm2nd_ = dbe_->book1D(histo, histo, 110, -0.0005, 0.005);
00091 
00092       sprintf (histo, "ES Multiplicity" );
00093       meESRecHitsMultip_ = dbe_->book1D(histo, histo, 100, 0., 700.);
00094 
00095       sprintf (histo, "ES Multiplicity Plane1 Side+");
00096       meESRecHitsMultip_zp1st_ = dbe_->book1D(histo, histo, 100, 0., 700.);
00097 
00098       sprintf (histo, "ES Multiplicity Plane2 Side+");
00099       meESRecHitsMultip_zp2nd_ = dbe_->book1D(histo, histo, 100, 0., 700.);
00100 
00101       sprintf (histo, "ES Multiplicity Plane1 Side-");
00102       meESRecHitsMultip_zm1st_ = dbe_->book1D(histo, histo, 100, 0., 700.);
00103 
00104       sprintf (histo, "ES Multiplicity Plane2 Side-");
00105       meESRecHitsMultip_zm2nd_ = dbe_->book1D(histo, histo, 100, 0., 700.);
00106 
00107       sprintf (histo, "Preshower EE vs ES energy Side+");
00108       meESEERecHitsEnergy_zp_ = dbe_->book2D(histo, histo, 100, 0., 0.1, 100, 0., 40.);
00109 
00110       sprintf (histo, "Preshower EE vs ES energy Side-");
00111       meESEERecHitsEnergy_zm_ = dbe_->book2D(histo, histo, 100, 0., 0.1, 100, 0., 40.);
00112 
00113       for (int kk=0; kk<32; kk++)
00114         { 
00115           sprintf(histo, "ES Occupancy Plane1 Side+ Strip%02d", kk+1);    
00116           meESRecHitsStripOccupancy_zp1st_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
00117 
00118           sprintf(histo, "ES Occupancy Plane2 Side+ Strip%02d", kk+1);    
00119           meESRecHitsStripOccupancy_zp2nd_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
00120 
00121           sprintf(histo, "ES Occupancy Plane1 Side- Strip%02d", kk+1);    
00122           meESRecHitsStripOccupancy_zm1st_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
00123 
00124           sprintf(histo, "ES Occupancy Plane2 Side- Strip%02d", kk+1);    
00125           meESRecHitsStripOccupancy_zm2nd_[kk] = dbe_->book2D(histo, histo, 40, 0., 40., 40, 0., 40.);
00126         }
00127     }
00128 }

EcalPreshowerRecHitsValidation::~EcalPreshowerRecHitsValidation ( )

Destructor.

Definition at line 130 of file EcalPreshowerRecHitsValidation.cc.

00130                                                                {   
00131 
00132 }

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(), EErechitCollection_, EEuncalibrechitCollection_, edm::SortedCollection< T, SORT >::end(), ESrechitCollection_, MonitorElement::Fill(), edm::SortedCollection< T, SORT >::find(), edm::Event::getByLabel(), edm::Handle< T >::isValid(), meESEERecHitsEnergy_zm_, meESEERecHitsEnergy_zp_, meESRecHitsEnergy_, meESRecHitsEnergy_zm1st_, meESRecHitsEnergy_zm2nd_, meESRecHitsEnergy_zp1st_, meESRecHitsEnergy_zp2nd_, meESRecHitsMultip_, meESRecHitsMultip_zm1st_, meESRecHitsMultip_zm2nd_, meESRecHitsMultip_zp1st_, meESRecHitsMultip_zp2nd_, meESRecHitsStripOccupancy_zm1st_, meESRecHitsStripOccupancy_zm2nd_, meESRecHitsStripOccupancy_zp1st_, meESRecHitsStripOccupancy_zp2nd_, ESDetId::plane(), edm::Handle< T >::product(), ESDetId::six(), ESDetId::siy(), ESDetId::strip(), strip(), ESDetId::zside(), and EEDetId::zside().

00142                                                                                {
00143 
00144   const ESRecHitCollection *ESRecHit = 0;
00145   Handle<ESRecHitCollection> EcalRecHitES;
00146   e.getByLabel( ESrechitCollection_, EcalRecHitES);
00147   if (EcalRecHitES.isValid()) {
00148     ESRecHit = EcalRecHitES.product ();
00149   } else {
00150     return;
00151   }
00152 
00153   bool skipEE = false;
00154   const EERecHitCollection *EERecHit = 0;
00155   Handle<EERecHitCollection> EcalRecHitEE;
00156   e.getByLabel( EErechitCollection_, EcalRecHitEE);
00157   if (EcalRecHitEE.isValid()){   
00158     EERecHit = EcalRecHitEE.product ();  
00159   } else {
00160     skipEE = true;
00161   }
00162 
00163   const EEUncalibratedRecHitCollection *EEUncalibRecHit = 0;
00164   Handle< EEUncalibratedRecHitCollection > EcalUncalibRecHitEE;
00165   e.getByLabel( EEuncalibrechitCollection_, EcalUncalibRecHitEE);
00166   if (EcalUncalibRecHitEE.isValid()) {
00167     EEUncalibRecHit = EcalUncalibRecHitEE.product() ;
00168   } else {
00169     skipEE = true;
00170   }
00171 
00172 
00173 
00174   // ---------------------- 
00175   // loop over RecHits
00176   // multiplicities
00177   int mult_tot       = 0;
00178   int mult_zp1st     = 0;
00179   int mult_zp2nd     = 0;
00180   int mult_zm1st     = 0;
00181   int mult_zm2nd     = 0;
00182 
00183   // energies
00184   float ene_zp1st      = 0.;
00185   float ene_zp2nd      = 0.;
00186   float ene_zm1st      = 0.;
00187   float ene_zm2nd      = 0.;
00188 
00189 
00190   // ES
00191   for (ESRecHitCollection::const_iterator recHit = ESRecHit->begin(); recHit != ESRecHit->end() ; ++recHit)
00192     {
00193       ESDetId ESid = ESDetId(recHit->id());
00194 
00195       int zside = ESid.zside();
00196       int plane = ESid.plane();
00197       int six   = ESid.six();
00198       int siy   = ESid.siy();
00199       int strip = ESid.strip();
00200       
00201       // global
00202       mult_tot++;
00203       if (meESRecHitsEnergy_) meESRecHitsEnergy_ ->Fill(recHit->energy());      
00204       
00205       // side +, plane 1
00206       if ( (zside == +1) && (plane == 1) )
00207         { 
00208           mult_zp1st++;
00209           ene_zp1st += recHit->energy();
00210           if ( meESRecHitsEnergy_zp1st_ )                 { meESRecHitsEnergy_zp1st_                 -> Fill(recHit->energy()); }
00211           if ( meESRecHitsStripOccupancy_zp1st_[strip-1] ){ meESRecHitsStripOccupancy_zp1st_[strip-1]-> Fill( six, siy ); }
00212         }
00213       
00214 
00215       // side +, plane 2
00216       if ( (zside == +1) && (plane == 2) )
00217         { 
00218           mult_zp2nd++; 
00219           ene_zp2nd += recHit->energy();
00220           if ( meESRecHitsEnergy_zp2nd_ )                 { meESRecHitsEnergy_zp2nd_                 -> Fill(recHit->energy()); }
00221           if ( meESRecHitsStripOccupancy_zp2nd_[strip-1] ){ meESRecHitsStripOccupancy_zp2nd_[strip-1]-> Fill( six, siy ); }
00222         }
00223 
00224 
00225       // side -, plane 1
00226       if ( (zside == -1) && (plane == 1) )
00227         { 
00228           mult_zm1st++;
00229           ene_zm1st += recHit->energy(); 
00230           if ( meESRecHitsEnergy_zm1st_ )                 { meESRecHitsEnergy_zm1st_                 -> Fill(recHit->energy()); }
00231           if ( meESRecHitsStripOccupancy_zm1st_[strip-1] ){ meESRecHitsStripOccupancy_zm1st_[strip-1]-> Fill( six, siy ); }
00232         }
00233 
00234 
00235       // side +, plane 2
00236       if ( (zside == -1) && (plane == 2) )
00237         { 
00238           mult_zm2nd ++; 
00239           ene_zm2nd += recHit->energy();
00240           if ( meESRecHitsEnergy_zm2nd_ )                 { meESRecHitsEnergy_zm2nd_                 -> Fill(recHit->energy()); }
00241           if ( meESRecHitsStripOccupancy_zm2nd_[strip-1] ){ meESRecHitsStripOccupancy_zm2nd_[strip-1]-> Fill( six, siy ); }
00242         }
00243 
00244     }  // loop over the ES RecHitCollection
00245 
00246 
00247 
00248   // EE
00249   double zpEE = 0.;
00250   double zmEE = 0.;
00251   if ( ! skipEE ) {
00252     
00253     for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EEUncalibRecHit->begin(); uncalibRecHit != EEUncalibRecHit->end() ; ++uncalibRecHit)
00254       {
00255         EEDetId EEid = EEDetId(uncalibRecHit->id());
00256         int mySide = EEid.zside();
00257         
00258         // Find corresponding recHit
00259         EcalRecHitCollection::const_iterator myRecHit = EERecHit->find(EEid);
00260         
00261         if (myRecHit != EERecHit->end() )
00262           {
00263             if (mySide > 0) { zpEE = zpEE + myRecHit->energy(); }
00264             if (mySide < 0) { zmEE = zmEE + myRecHit->energy(); }
00265           }
00266       } 
00267   }
00268   
00269 
00270   // filling histos
00271   if (meESRecHitsMultip_)            { meESRecHitsMultip_           -> Fill(mult_tot);   }
00272   if (meESRecHitsMultip_zp1st_ )     { meESRecHitsMultip_zp1st_     -> Fill(mult_zp1st); }
00273   if (meESRecHitsMultip_zp2nd_ )     { meESRecHitsMultip_zp2nd_     -> Fill(mult_zp2nd); }
00274   if (meESRecHitsMultip_zm1st_ )     { meESRecHitsMultip_zm1st_     -> Fill(mult_zm1st); }
00275   if (meESRecHitsMultip_zm2nd_ )     { meESRecHitsMultip_zm2nd_     -> Fill(mult_zm2nd); }
00276   if (meESEERecHitsEnergy_zp_)       { meESEERecHitsEnergy_zp_      -> Fill( (ene_zp1st + 0.7*ene_zp2nd)/0.09, zpEE ); }
00277   if (meESEERecHitsEnergy_zm_)       { meESEERecHitsEnergy_zm_      -> Fill( (ene_zm1st + 0.7*ene_zm2nd)/0.09, zmEE ); }
00278 }

void EcalPreshowerRecHitsValidation::beginJob ( const edm::EventSetup & c ) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 134 of file EcalPreshowerRecHitsValidation.cc.

00134                                                                 {  
00135 
00136 }

void EcalPreshowerRecHitsValidation::endJob ( void ) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 138 of file EcalPreshowerRecHitsValidation.cc.

00138                                            {
00139 
00140 }

DQMStore* EcalPreshowerRecHitsValidation::dbe_ [private]

Definition at line 60 of file EcalPreshowerRecHitsValidation.h.

Referenced by EcalPreshowerRecHitsValidation().

edm::InputTag EcalPreshowerRecHitsValidation::EErechitCollection_ [private]

Definition at line 63 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

edm::InputTag EcalPreshowerRecHitsValidation::EEuncalibrechitCollection_ [private]

Definition at line 62 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

edm::InputTag EcalPreshowerRecHitsValidation::ESrechitCollection_ [private]

Definition at line 64 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESEERecHitsEnergy_zm_ [private]

Definition at line 77 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESEERecHitsEnergy_zp_ [private]

Definition at line 76 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsEnergy_ [private]

Definition at line 66 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsEnergy_zm1st_ [private]

Definition at line 69 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsEnergy_zm2nd_ [private]

Definition at line 70 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsEnergy_zp1st_ [private]

Definition at line 67 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsEnergy_zp2nd_ [private]

Definition at line 68 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsMultip_ [private]

Definition at line 71 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsMultip_zm1st_ [private]

Definition at line 74 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsMultip_zm2nd_ [private]

Definition at line 75 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsMultip_zp1st_ [private]

Definition at line 72 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsMultip_zp2nd_ [private]

Definition at line 73 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsStripOccupancy_zm1st_[36] [private]

Definition at line 79 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsStripOccupancy_zm2nd_[36] [private]

Definition at line 81 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsStripOccupancy_zp1st_[36] [private]

Definition at line 78 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

MonitorElement* EcalPreshowerRecHitsValidation::meESRecHitsStripOccupancy_zp2nd_[36] [private]

Definition at line 80 of file EcalPreshowerRecHitsValidation.h.

Referenced by analyze(), and EcalPreshowerRecHitsValidation().

bool EcalPreshowerRecHitsValidation::verbose_ [private]

Definition at line 58 of file EcalPreshowerRecHitsValidation.h.

Referenced by EcalPreshowerRecHitsValidation().


Public Member Functions
	EcalPreshowerRecHitsValidation (const edm::ParameterSet &ps)
	Constructor.
	~EcalPreshowerRecHitsValidation ()
	Destructor.
Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze.
void	beginJob (const edm::EventSetup &c)
void	endJob (void)
Private Attributes
DQMStore *	dbe_
edm::InputTag	EErechitCollection_
edm::InputTag	EEuncalibrechitCollection_
edm::InputTag	ESrechitCollection_
MonitorElement *	meESEERecHitsEnergy_zm_
MonitorElement *	meESEERecHitsEnergy_zp_
MonitorElement *	meESRecHitsEnergy_
MonitorElement *	meESRecHitsEnergy_zm1st_
MonitorElement *	meESRecHitsEnergy_zm2nd_
MonitorElement *	meESRecHitsEnergy_zp1st_
MonitorElement *	meESRecHitsEnergy_zp2nd_
MonitorElement *	meESRecHitsMultip_
MonitorElement *	meESRecHitsMultip_zm1st_
MonitorElement *	meESRecHitsMultip_zm2nd_
MonitorElement *	meESRecHitsMultip_zp1st_
MonitorElement *	meESRecHitsMultip_zp2nd_
MonitorElement *	meESRecHitsStripOccupancy_zm1st_ [36]
MonitorElement *	meESRecHitsStripOccupancy_zm2nd_ [36]
MonitorElement *	meESRecHitsStripOccupancy_zp1st_ [36]
MonitorElement *	meESRecHitsStripOccupancy_zp2nd_ [36]
bool	verbose_

EcalPreshowerRecHitsValidation Class Reference

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation