CMS 3D CMS Logo

EcalPreshowerSimHitsValidation Class Reference

#include <Validation/EcalHits/interface/EcalPreshowerSimHitsValidation.h>

Inheritance diagram for EcalPreshowerSimHitsValidation:

edm::EDAnalyzer

List of all members.

Public Member Functions

 EcalPreshowerSimHitsValidation (const edm::ParameterSet &ps)
 Constructor.
 ~EcalPreshowerSimHitsValidation ()
 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 Types

typedef std::map< uint32_t,
float, std::less< uint32_t > > 		MapType

Private Attributes

DQMStore * dbe_
std::string EEHitsCollection
std::string ESHitsCollection
std::string g4InfoLabel
std::string HepMCLabel
MonitorElement * me2eszmOver1eszm_
MonitorElement * me2eszpOver1eszp_
MonitorElement * meE1alphaE2zm_
MonitorElement * meE1alphaE2zp_
MonitorElement * meEEoverESzm_
MonitorElement * meEEoverESzp_
MonitorElement * meESEnergyHits1zm_
MonitorElement * meESEnergyHits1zp_
MonitorElement * meESEnergyHits2zm_
MonitorElement * meESEnergyHits2zp_
MonitorElement * meEShitLog10Energy_
MonitorElement * meEShitLog10EnergyNorm_
MonitorElement * menESHits1zm_
MonitorElement * menESHits1zp_
MonitorElement * menESHits2zm_
MonitorElement * menESHits2zp_
std::string outputFile_
bool verbose_

Detailed Description

Definition at line 37 of file EcalPreshowerSimHitsValidation.h.

Member Typedef Documentation

typedef std::map<uint32_t,float,std::less<uint32_t> > EcalPreshowerSimHitsValidation::MapType [private]

Definition at line 39 of file EcalPreshowerSimHitsValidation.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 18 of file EcalPreshowerSimHitsValidation.cc.

References dbe_, edm::ParameterSet::getUntrackedParameter(), histo, me2eszmOver1eszm_, me2eszpOver1eszp_, meE1alphaE2zm_, meE1alphaE2zp_, meEEoverESzm_, meEEoverESzp_, meESEnergyHits1zm_, meESEnergyHits1zp_, meESEnergyHits2zm_, meESEnergyHits2zp_, meEShitLog10Energy_, meEShitLog10EnergyNorm_, menESHits1zm_, menESHits1zp_, menESHits2zm_, menESHits2zp_, and verbose_.

00018                                                                                        :
00019 
00020   HepMCLabel(ps.getParameter<std::string>("moduleLabelMC")),
00021   g4InfoLabel(ps.getParameter<std::string>("moduleLabelG4")),
00022   EEHitsCollection(ps.getParameter<std::string>("EEHitsCollection")),
00023   ESHitsCollection(ps.getParameter<std::string>("ESHitsCollection")){
00024 
00025   // verbosity switch
00026   verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
00027  
00028   // get hold of back-end interface
00029   dbe_ = 0;
00030   dbe_ = edm::Service<DQMStore>().operator->();
00031   if ( dbe_ ) {
00032     if ( verbose_ ) { dbe_->setVerbose(1); } 
00033     else            { dbe_->setVerbose(0); }
00034   }
00035 
00036   if ( dbe_ ) {
00037     if ( verbose_ ) dbe_->showDirStructure();
00038   }
00039 
00040 
00041   menESHits1zp_ = 0;     
00042   menESHits2zp_ = 0;     
00043   menESHits1zm_ = 0;     
00044   menESHits2zm_ = 0;     
00045                                     
00046   meESEnergyHits1zp_ = 0;
00047   meESEnergyHits2zp_ = 0;
00048   meESEnergyHits1zm_ = 0;
00049   meESEnergyHits2zm_ = 0;
00050 
00051   meEShitLog10Energy_       = 0;
00052   meEShitLog10EnergyNorm_   = 0;
00053 
00054   meE1alphaE2zp_ = 0;
00055   meE1alphaE2zm_ = 0;
00056   meEEoverESzp_  = 0;
00057   meEEoverESzm_  = 0;
00058 
00059   me2eszpOver1eszp_ = 0;
00060   me2eszmOver1eszm_ = 0;
00061 
00062 
00063   Char_t histo[200];
00064  
00065   if ( dbe_ ) {
00066     dbe_->setCurrentFolder("EcalHitsV/EcalSimHitsValidation");
00067   
00068     sprintf (histo, "ES hits layer 1 multiplicity z+" ) ;
00069     menESHits1zp_ = dbe_->book1D(histo, histo, 50, 0., 50. ) ;
00070 
00071     sprintf (histo, "ES hits layer 2 multiplicity z+" ) ;
00072     menESHits2zp_ = dbe_->book1D(histo, histo, 50, 0., 50. ) ;
00073 
00074     sprintf (histo, "ES hits layer 1 multiplicity z-" ) ;
00075     menESHits1zm_ = dbe_->book1D(histo, histo, 50, 0., 50. ) ;
00076 
00077     sprintf (histo, "ES hits layer 2 multiplicity z-" ) ;
00078     menESHits2zm_ = dbe_->book1D(histo, histo, 50, 0., 50. ) ;
00079 
00080     sprintf (histo, "ES hits energy layer 1 z+" ) ;
00081     meESEnergyHits1zp_ = dbe_->book1D(histo, histo, 100, 0., 0.05 ) ;
00082 
00083     sprintf (histo, "ES hits energy layer 2 z+" ) ;
00084     meESEnergyHits2zp_ = dbe_->book1D(histo, histo, 100, 0., 0.05 ) ;
00085 
00086     sprintf (histo, "ES hits energy layer 1 z-" ) ;
00087     meESEnergyHits1zm_ = dbe_->book1D(histo, histo, 100, 0., 0.05 ) ;
00088 
00089     sprintf (histo, "ES hits energy layer 2 z-" ) ;
00090     meESEnergyHits2zm_ = dbe_->book1D(histo, histo, 100, 0., 0.05 ) ;
00091 
00092     sprintf (histo, "ES hits log10energy spectrum" );
00093     meEShitLog10Energy_ = dbe_->book1D(histo, histo, 140, -10., 4.);
00094 
00095     sprintf (histo, "ES hits log10energy spectrum vs normalized energy" );
00096     meEShitLog10EnergyNorm_ = dbe_->bookProfile(histo, histo, 140, -10., 4., 100, 0., 1.);
00097 
00098     sprintf (histo, "ES E1+07E2 z+" ) ;
00099     meE1alphaE2zp_ = dbe_->book1D(histo, histo, 100, 0., 0.05);
00100 
00101     sprintf (histo, "ES E1+07E2 z-" ) ;
00102     meE1alphaE2zm_ = dbe_->book1D(histo, histo, 100, 0., 0.05);
00103 
00104     sprintf (histo, "EE vs ES z+" ) ;
00105     meEEoverESzp_ = dbe_->bookProfile(histo, histo, 250, 0., 500., 200, 0., 200.);
00106 
00107     sprintf (histo, "EE vs ES z-" ) ;
00108     meEEoverESzm_ = dbe_->bookProfile(histo, histo, 250, 0., 500., 200, 0., 200.);
00109 
00110     sprintf (histo, "ES ene2oEne1 z+" ) ;
00111     me2eszpOver1eszp_ = dbe_->book1D(histo, histo, 50, 0., 10.);
00112 
00113     sprintf (histo, "ES ene2oEne1 z-" ) ;
00114     me2eszmOver1eszm_ = dbe_->book1D(histo, histo, 50, 0., 10.);
00115   }
00116 }

EcalPreshowerSimHitsValidation::~EcalPreshowerSimHitsValidation (  ) 

Destructor.

Definition at line 118 of file EcalPreshowerSimHitsValidation.cc.

00118                                                                {
00119 
00120 }

Member Function Documentation

void EcalPreshowerSimHitsValidation::analyze ( const edm::Event &  e,
const edm::EventSetup &  c 
) [protected, virtual]

Analyze.

Implements edm::EDAnalyzer.

Definition at line 130 of file EcalPreshowerSimHitsValidation.cc.

References EEHitsCollection, ESHitsCollection, edm::EventID::event(), MonitorElement::Fill(), g4InfoLabel, edm::Event::getByLabel(), HepMCLabel, i, edm::Event::id(), int, edm::Handle< T >::isValid(), funct::log(), LogDebug, me2eszmOver1eszm_, me2eszpOver1eszp_, meE1alphaE2zm_, meE1alphaE2zp_, meEEoverESzm_, meEEoverESzp_, meESEnergyHits1zm_, meESEnergyHits1zp_, meESEnergyHits2zm_, meESEnergyHits2zp_, meEShitLog10Energy_, meEShitLog10EnergyNorm_, menESHits1zm_, menESHits1zp_, menESHits2zm_, menESHits2zp_, p, edm::EventID::run(), and funct::tan().

00130                                                                                      {
00131 
00132   edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
00133   
00134   edm::Handle<edm::HepMCProduct> MCEvt;
00135   e.getByLabel(HepMCLabel, MCEvt);
00136 
00137   edm::Handle<edm::PCaloHitContainer> EcalHitsEE;
00138   e.getByLabel(g4InfoLabel,EEHitsCollection,EcalHitsEE);
00139 
00140   edm::Handle<edm::PCaloHitContainer> EcalHitsES;
00141   e.getByLabel(g4InfoLabel,ESHitsCollection,EcalHitsES);
00142 
00143   std::vector<PCaloHit> theEECaloHits;  
00144   if( EcalHitsEE.isValid() ) {
00145     theEECaloHits.insert(theEECaloHits.end(), EcalHitsEE->begin(), EcalHitsEE->end());
00146   }
00147 
00148   std::vector<PCaloHit> theESCaloHits;
00149   if( EcalHitsES.isValid() ) {
00150     theESCaloHits.insert(theESCaloHits.end(), EcalHitsES->begin(), EcalHitsES->end());
00151   }
00152 
00153   double ESEnergy_ = 0.;
00154   //std::map<unsigned int, std::vector<PCaloHit>,std::less<unsigned int> > CaloHitMap;
00155 
00156   // endcap
00157   double EEetzp_ = 0.;
00158   double EEetzm_ = 0.;
00159   for (std::vector<PCaloHit>::iterator isim = theEECaloHits.begin(); isim != theEECaloHits.end(); ++isim){
00160     EEDetId eeid (isim->id()) ;
00161     if (eeid.zside() > 0 ) EEetzp_ += isim->energy();
00162     if (eeid.zside() < 0 ) EEetzm_ += isim->energy();
00163   }
00164   
00165   
00166   uint32_t nESHits1zp = 0;
00167   uint32_t nESHits1zm = 0;
00168   uint32_t nESHits2zp = 0;
00169   uint32_t nESHits2zm = 0;
00170   double ESet1zp_ = 0.;
00171   double ESet2zp_ = 0.;
00172   double ESet1zm_ = 0.;
00173   double ESet2zm_ = 0.;
00174   std::vector<double> econtr(140, 0. );
00175   
00176   for (std::vector<PCaloHit>::iterator isim = theESCaloHits.begin();
00177        isim != theESCaloHits.end(); ++isim){
00178     //CaloHitMap[(*isim).id()].push_back((*isim));
00179     
00180     ESDetId esid (isim->id()) ;
00181     
00182     LogDebug("HitInfo")
00183       << " CaloHit " << isim->getName() << "\n" 
00184       << " DetID = "<<isim->id()<< " ESDetId: z side " << esid.zside() << "  plane " << esid.plane() << esid.six() << ',' << esid.siy() << ':' << esid.strip() << "\n"
00185       << " Time = " << isim->time() << "\n"
00186       << " Track Id = " << isim->geantTrackId() << "\n"
00187       << " Energy = " << isim->energy();
00188     
00189     ESEnergy_ += isim->energy();
00190     meEShitLog10Energy_->Fill(log10(isim->energy()));
00191   
00192     int log10i = int( ( log10(isim->energy()) + 10. ) * 10. );
00193     if( log10i >=0 && log10i < 140 ) econtr[log10i] += isim->energy();
00194   
00195     if (esid.plane() == 1 ) { 
00196       if (esid.zside() > 0 ) {
00197         nESHits1zp++ ;
00198         ESet1zp_ += isim->energy();
00199         if (meESEnergyHits1zp_) meESEnergyHits1zp_->Fill(isim->energy()) ; 
00200       }
00201       else if (esid.zside() < 0 ) {
00202         nESHits1zm++ ; 
00203         ESet1zm_ += isim->energy();
00204         if (meESEnergyHits1zm_) meESEnergyHits1zm_->Fill(isim->energy()) ; 
00205       }
00206     }
00207     else if (esid.plane() == 2 ) {
00208       if (esid.zside() > 0 ) {
00209         nESHits2zp++ ; 
00210         ESet2zp_ += isim->energy();
00211         if (meESEnergyHits2zp_) meESEnergyHits2zp_->Fill(isim->energy()) ; 
00212       }
00213       else if (esid.zside() < 0 ) {
00214         nESHits2zm++ ; 
00215         ESet2zm_ += isim->energy();
00216         if (meESEnergyHits2zm_) meESEnergyHits2zm_->Fill(isim->energy()) ; 
00217       }
00218     }
00219     
00220   }
00221   
00222   if (menESHits1zp_) menESHits1zp_->Fill(nESHits1zp);
00223   if (menESHits1zm_) menESHits1zm_->Fill(nESHits1zm);
00224   
00225   if (menESHits2zp_) menESHits2zp_->Fill(nESHits2zp);
00226   if (menESHits2zm_) menESHits2zm_->Fill(nESHits2zm);
00227   
00228   if( meEShitLog10EnergyNorm_ && ESEnergy_ != 0 ) {
00229     for( int i=0; i<140; i++ ) {
00230       meEShitLog10EnergyNorm_->Fill( -10.+(float(i)+0.5)/10., econtr[i]/ESEnergy_ );
00231     }
00232   }
00233 
00234 
00235   for ( HepMC::GenEvent::particle_const_iterator p = MCEvt->GetEvent()->particles_begin();
00236         p != MCEvt->GetEvent()->particles_end(); ++p ) {
00237     
00238     double htheta = (*p)->momentum().theta();
00239     double heta = -log(tan(htheta * 0.5));
00240     
00241     if ( heta > 1.653 && heta < 2.6 ) {
00242 
00243       if (meE1alphaE2zp_) meE1alphaE2zp_->Fill(ESet1zp_+0.7*ESet2zp_);
00244       if (meEEoverESzp_)  meEEoverESzp_ ->Fill((ESet1zp_+0.7*ESet2zp_)/0.00009, EEetzp_);
00245       if ((me2eszpOver1eszp_) && (ESet1zp_ != 0.)) me2eszpOver1eszp_->Fill(ESet2zp_/ESet1zp_);
00246     }
00247     if ( heta < -1.653 && heta > -2.6 ) {
00248       if (meE1alphaE2zm_) meE1alphaE2zm_->Fill(ESet1zm_+0.7*ESet2zm_);
00249       if (meEEoverESzm_)  meEEoverESzm_ ->Fill((ESet1zm_+0.7*ESet2zm_)/0.00009, EEetzm_);
00250       if ((me2eszmOver1eszm_) && (ESet1zm_ != 0.)) me2eszmOver1eszm_->Fill(ESet2zm_/ESet1zm_);
00251     }
00252   }
00253 }

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

Reimplemented from edm::EDAnalyzer.

Definition at line 122 of file EcalPreshowerSimHitsValidation.cc.

00122                                                                    {
00123 
00124 }

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

Reimplemented from edm::EDAnalyzer.

Definition at line 126 of file EcalPreshowerSimHitsValidation.cc.

00126                                            {
00127 
00128 }

Member Data Documentation

DQMStore* EcalPreshowerSimHitsValidation::dbe_ [private]

Definition at line 69 of file EcalPreshowerSimHitsValidation.h.

Referenced by EcalPreshowerSimHitsValidation().

std::string EcalPreshowerSimHitsValidation::EEHitsCollection [private]

Definition at line 64 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze().

std::string EcalPreshowerSimHitsValidation::ESHitsCollection [private]

Definition at line 65 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze().

std::string EcalPreshowerSimHitsValidation::g4InfoLabel [private]

Definition at line 63 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze().

std::string EcalPreshowerSimHitsValidation::HepMCLabel [private]

Definition at line 62 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze().

MonitorElement* EcalPreshowerSimHitsValidation::me2eszmOver1eszm_ [private]

Definition at line 95 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::me2eszpOver1eszp_ [private]

Definition at line 94 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meE1alphaE2zm_ [private]

Definition at line 89 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meE1alphaE2zp_ [private]

Definition at line 88 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meEEoverESzm_ [private]

Definition at line 92 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meEEoverESzp_ [private]

Definition at line 91 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meESEnergyHits1zm_ [private]

Definition at line 85 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meESEnergyHits1zp_ [private]

Definition at line 82 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meESEnergyHits2zm_ [private]

Definition at line 86 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meESEnergyHits2zp_ [private]

Definition at line 83 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meEShitLog10Energy_ [private]

Definition at line 79 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::meEShitLog10EnergyNorm_ [private]

Definition at line 80 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::menESHits1zm_ [private]

Definition at line 76 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::menESHits1zp_ [private]

Definition at line 73 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::menESHits2zm_ [private]

Definition at line 77 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

MonitorElement* EcalPreshowerSimHitsValidation::menESHits2zp_ [private]

Definition at line 74 of file EcalPreshowerSimHitsValidation.h.

Referenced by analyze(), and EcalPreshowerSimHitsValidation().

std::string EcalPreshowerSimHitsValidation::outputFile_ [private]

Definition at line 71 of file EcalPreshowerSimHitsValidation.h.

bool EcalPreshowerSimHitsValidation::verbose_ [private]

Definition at line 67 of file EcalPreshowerSimHitsValidation.h.

Referenced by EcalPreshowerSimHitsValidation().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:19:48 2009 for CMSSW by  doxygen 1.5.4