#include <EcalTBValidation.h>

Inheritance diagram for EcalTBValidation:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
	EcalTBValidation (const edm::ParameterSet &)
virtual void	endJob ()
	~EcalTBValidation ()
Private Attributes
int	data_
DQMStore *	dbe_
std::string	digiCollection_
std::string	digiProducer_
std::string	eventHeaderCollection_
std::string	eventHeaderProducer_
std::string	hitCollection_
std::string	hitProducer_
std::string	hodoRecInfoCollection_
std::string	hodoRecInfoProducer_
MonitorElement *	meETBampltdc_
MonitorElement *	meETBe1e25_
MonitorElement *	meETBe1e25vsX_
MonitorElement *	meETBe1e25vsY_
MonitorElement *	meETBe1e9_
MonitorElement *	meETBe1e9vsX_
MonitorElement *	meETBe1e9vsY_
MonitorElement *	meETBe1vsX_
MonitorElement *	meETBe1vsY_
MonitorElement *	meETBe1x1_
MonitorElement *	meETBe1x1_center_
MonitorElement *	meETBe3x3_
MonitorElement *	meETBe3x3_center_
MonitorElement *	meETBe5x5_
MonitorElement *	meETBe5x5_center_
MonitorElement *	meETBe9e25_
MonitorElement *	meETBe9e25vsX_
MonitorElement *	meETBe9e25vsY_
MonitorElement *	meETBhodoX_
MonitorElement *	meETBhodoY_
MonitorElement *	meETBShape_
MonitorElement *	meETBxib_
std::string	tdcRecInfoCollection_
std::string	tdcRecInfoProducer_
bool	verbose_
int	xtalInBeam_

Detailed Description

Definition at line 21 of file EcalTBValidation.h.

Constructor & Destructor Documentation

EcalTBValidation::EcalTBValidation ( const edm::ParameterSet & config ) [explicit]

Definition at line 17 of file EcalTBValidation.cc.

References data_, dbe_, digiCollection_, digiProducer_, eventHeaderCollection_, eventHeaderProducer_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hitCollection_, hitProducer_, hodoRecInfoCollection_, hodoRecInfoProducer_, meETBampltdc_, meETBe1e25_, meETBe1e25vsX_, meETBe1e25vsY_, meETBe1e9_, meETBe1e9vsX_, meETBe1e9vsY_, meETBe1vsX_, meETBe1vsY_, meETBe1x1_, meETBe1x1_center_, meETBe3x3_, meETBe3x3_center_, meETBe5x5_, meETBe5x5_center_, meETBe9e25_, meETBe9e25vsX_, meETBe9e25vsY_, meETBhodoX_, meETBhodoY_, meETBShape_, meETBxib_, cmsCodeRules::cppFunctionSkipper::operator, tdcRecInfoCollection_, tdcRecInfoProducer_, verbose_, and xtalInBeam_.

                                                                  {  

  data_                  = config.getUntrackedParameter<int>("data",-1000);
  xtalInBeam_            = config.getUntrackedParameter<int>("xtalInBeam",-1000);

  digiCollection_        = config.getParameter<std::string>("digiCollection");
  digiProducer_          = config.getParameter<std::string>("digiProducer");
  hitCollection_         = config.getParameter<std::string>("hitCollection");
  hitProducer_           = config.getParameter<std::string>("hitProducer");
  hodoRecInfoCollection_ = config.getParameter<std::string>("hodoRecInfoCollection");
  hodoRecInfoProducer_   = config.getParameter<std::string>("hodoRecInfoProducer");
  tdcRecInfoCollection_  = config.getParameter<std::string>("tdcRecInfoCollection");
  tdcRecInfoProducer_    = config.getParameter<std::string>("tdcRecInfoProducer");
  eventHeaderCollection_ = config.getParameter<std::string>("eventHeaderCollection");
  eventHeaderProducer_   = config.getParameter<std::string>("eventHeaderProducer");

  //rootfile_              = config.getUntrackedParameter<std::string>("rootfile","EcalTBValidation.root");

  // verbosity...
  verbose_ = config.getUntrackedParameter<bool>("verbose", false);

  dbe_ = edm::Service<DQMStore>().operator->();
  if( dbe_ ) {
    if( verbose_ ) {
      dbe_->setVerbose(1);
      dbe_->showDirStructure();
    }
    else {
      dbe_->setVerbose(0);
    }
  }
    
  meETBxib_ = 0;
  meETBampltdc_ = 0;
  meETBShape_ = 0;
  meETBhodoX_ = 0;
  meETBhodoY_ = 0;
  meETBe1x1_ = 0;
  meETBe3x3_ = 0;
  meETBe5x5_ = 0;
  meETBe1e9_ = 0;
  meETBe1e25_ = 0;
  meETBe9e25_ = 0;
  meETBe1x1_center_ = 0;
  meETBe3x3_center_ = 0;
  meETBe5x5_center_ = 0;
  meETBe1vsX_ = 0;
  meETBe1vsY_ = 0;
  meETBe1e9vsX_ = 0;
  meETBe1e9vsY_ = 0;
  meETBe1e25vsX_ = 0;
  meETBe1e25vsY_ = 0;
  meETBe9e25vsX_ = 0;
  meETBe9e25vsY_ = 0;

  if( dbe_ ) {

    std::string hname;   
    dbe_->setCurrentFolder( "EcalRecHitsV/EcalTBValidationTask" );

    hname = "xtal in beam position";
    meETBxib_          = dbe_->book2D( hname, hname, 85, 0., 85., 20,0., 20. );
    hname = "Max Amplitude vs TDC offset";
    meETBampltdc_      = dbe_->book2D( hname, hname, 100, 0., 1., 1000, 0., 4000. );
    hname = "Beam Profile X";
    meETBhodoX_        = dbe_->book1D( hname, hname, 100, -20., 20. );
    hname = "Beam Profile Y";
    meETBhodoY_        = dbe_->book1D( hname, hname, 100, -20., 20. );
    hname = "E1x1 energy";
    meETBe1x1_         = dbe_->book1D( hname, hname, 1000, 0., 4000. );
    hname = "E3x3 energy";
    meETBe3x3_         = dbe_->book1D( hname, hname, 1000, 0., 4000. );
    hname = "E5x5 energy";
    meETBe5x5_         = dbe_->book1D( hname, hname, 1000, 0., 4000. );
    hname = "E1x1 energy center";
    meETBe1x1_center_  = dbe_->book1D( hname, hname, 1000, 0., 4000. );
    hname = "E3x3 energy center";
    meETBe3x3_center_  = dbe_->book1D( hname, hname, 1000, 0., 4000. );
    hname = "E5x5 energy center";
    meETBe5x5_center_  = dbe_->book1D( hname, hname, 1000, 0., 4000. );
    hname = "E1 over E9 ratio";
    meETBe1e9_         = dbe_->book1D( hname, hname, 600, 0., 1.2 );
    hname = "E1 over E25 ratio";
    meETBe1e25_        = dbe_->book1D( hname, hname, 600, 0., 1.2 );
    hname = "E9 over E25 ratio";
    meETBe9e25_        = dbe_->book1D( hname, hname, 600, 0., 1.2 );
    hname = "E1 vs X";
    meETBe1vsX_        = dbe_->book2D( hname, hname, 80, -20, 20, 1000, 0., 4000. );
    hname = "E1 vs Y";
    meETBe1vsY_        = dbe_->book2D( hname, hname, 80, -20, 20, 1000, 0., 4000. );  
    hname = "E1 over E9 vs X";
    meETBe1e9vsX_      = dbe_->book2D( hname, hname, 80, -20, 20, 600, 0., 1.2 );
    hname = "E1 over E9 vs Y";
    meETBe1e9vsY_      = dbe_->book2D( hname, hname, 80, -20, 20, 600, 0., 1.2 );
    hname = "E1 over E25 vs X";
    meETBe1e25vsX_     = dbe_->book2D( hname, hname, 80, -20, 20, 600, 0., 1.2 );
    hname = "E1 over E25 vs Y";
    meETBe1e25vsY_     = dbe_->book2D( hname, hname, 80, -20, 20, 600, 0., 1.2 );
    hname = "E9 over E25 vs X";
    meETBe9e25vsX_     = dbe_->book2D( hname, hname, 80, -20, 20, 600, 0., 1.2 );
    hname = "E9 over E25 vs Y";
    meETBe9e25vsY_     = dbe_->book2D( hname, hname, 80, -20, 20, 600, 0., 1.2 );
    hname = "Xtal in Beam Shape";
    meETBShape_        = dbe_->book2D( hname, hname, 250, 0, 10, 350, 0, 3500 );
  }

}

EcalTBValidation::~EcalTBValidation ( )

Definition at line 126 of file EcalTBValidation.cc.

{}

Member Function Documentation

void EcalTBValidation::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)		`[virtual]`

Implements edm::EDAnalyzer.

Definition at line 132 of file EcalTBValidation.cc.

References EcalMGPASample::adc(), benchmark_cfg::cerr, data_, digiCollection_, digiProducer_, jptDQMConfig_cff::eMax, EBDetId::ETAPHIMODE, eventHeaderProducer_, MonitorElement::Fill(), HcalObjRepresent::Fill(), hitCollection_, hitProducer_, hodoRecInfoCollection_, hodoRecInfoProducer_, EBDetId::ieta(), EBDetId::iphi(), EBDetId::ism(), meETBampltdc_, meETBe1e25_, meETBe1e25vsX_, meETBe1e25vsY_, meETBe1e9_, meETBe1e9vsX_, meETBe1e9vsY_, meETBe1vsX_, meETBe1vsY_, meETBe1x1_, meETBe1x1_center_, meETBe3x3_, meETBe3x3_center_, meETBe5x5_, meETBe5x5_center_, meETBe9e25_, meETBe9e25vsX_, meETBe9e25vsY_, meETBhodoX_, meETBhodoY_, meETBShape_, meETBxib_, edm::Handle< T >::product(), EcalDataFrame::sample(), EcalDataFrame::size(), EBDetId::SMCRYSTALMODE, tdcRecInfoCollection_, tdcRecInfoProducer_, EBDetId::validDetId(), and xtalInBeam_.

                                                                                  {

  using namespace edm;
  using namespace cms;

  // digis
  const EBDigiCollection* theDigis=0;
  Handle<EBDigiCollection> pdigis;
  event.getByLabel(digiProducer_, digiCollection_, pdigis);
  if(pdigis.isValid()){
    theDigis = pdigis.product(); 
  } 
  else {    
    std::cerr << "Error! can't get the product " << digiCollection_.c_str() << std::endl;
    return;
  }

  // rechits
  const EBUncalibratedRecHitCollection* theHits=0;  
  Handle<EBUncalibratedRecHitCollection> phits;
  event.getByLabel(hitProducer_, hitCollection_, phits);
  if(phits.isValid()){
    theHits = phits.product(); 
  } 
  else {
    std::cerr << "Error! can't get the product " << hitCollection_.c_str() << std::endl;
    return;
  }

  // hodoscopes
  const EcalTBHodoscopeRecInfo* theHodo=0;  
  Handle<EcalTBHodoscopeRecInfo> pHodo;
  event.getByLabel(hodoRecInfoProducer_, hodoRecInfoCollection_, pHodo);
  if(pHodo.isValid()){ 
    theHodo = pHodo.product(); 
  }
  else{ 
    std::cerr << "Error! can't get the product " << hodoRecInfoCollection_.c_str() << std::endl;
    return;
  }
  
  // tdc
  const EcalTBTDCRecInfo* theTDC=0;
  Handle<EcalTBTDCRecInfo> pTDC;
  event.getByLabel(tdcRecInfoProducer_, tdcRecInfoCollection_, pTDC);
  if(pTDC.isValid()){
    theTDC = pTDC.product(); 
  }
  else{ 
    std::cerr << "Error! can't get the product " << tdcRecInfoCollection_.c_str() << std::endl;
    return;
  }

  // event header
  const EcalTBEventHeader* evtHeader=0;
  Handle<EcalTBEventHeader> pEventHeader;
  event.getByLabel(eventHeaderProducer_ , pEventHeader);
  if(pEventHeader.isValid()){
    evtHeader = pEventHeader.product(); 
  }
  else{ 
    std::cerr << "Error! can't get the product " << eventHeaderProducer_.c_str() << std::endl;
    return;
  }
  
 
  // -----------------------------------------------------------------------
  // xtal-in-beam
  EBDetId xtalInBeamId(1,xtalInBeam_, EBDetId::SMCRYSTALMODE);
  if (xtalInBeamId==EBDetId(0)){ return; }
  int xibEta = xtalInBeamId.ieta();
  int xibPhi = xtalInBeamId.iphi();

  // skipping events with moving table (in data)
  if (data_ && (evtHeader->tableIsMoving())) return;
  
  // amplitudes
  EBDetId Xtals5x5[25];
  for (unsigned int icry=0;icry<25;icry++){
    unsigned int row    = icry/5;
    unsigned int column = icry%5;
    int ieta = xtalInBeamId.ieta()+column-2;
    int iphi = xtalInBeamId.iphi()+row-2;
    if(EBDetId::validDetId(ieta, iphi)){ 
      EBDetId tempId(ieta, iphi,EBDetId::ETAPHIMODE);
      if (tempId.ism()==1) 
        Xtals5x5[icry] = tempId;
      else
        Xtals5x5[icry] = EBDetId(0);
    } else {
      Xtals5x5[icry] = EBDetId(0);   
    }
  } 
 
  // matrices
  double ampl1x1 = 0.;
  double ampl3x3 = 0.;
  double ampl5x5 = 0.;
  for (unsigned int icry=0;icry<25;icry++) {
    if (!Xtals5x5[icry].null()){
      double theAmpl = (theHits->find(Xtals5x5[icry]))->amplitude();
      ampl5x5 += theAmpl;
      if (icry==12){ampl1x1 = theAmpl;}
      if (icry==6 || icry==7 || icry==8 || icry==11 || icry==12 || icry==13 || icry==16 || icry==17 || icry==18){ampl3x3 += theAmpl;}
    }}


  // pulse shape
  double sampleSave[10];
  for(int ii=0; ii < 10; ++ii){ sampleSave[ii] = 0.0; }
  EBDigiCollection::const_iterator thisDigi = theDigis->find(xtalInBeamId);
  int sMax = -1;
  double eMax = 0.;
  if (thisDigi != theDigis->end()){
    EBDataFrame myDigi = (*thisDigi);
    for (int sample=0; sample < myDigi.size(); ++sample){
      double analogSample = myDigi.sample(sample).adc();
      sampleSave[sample]  = analogSample;
      if ( eMax < analogSample ) {
        eMax = analogSample;
        sMax = sample;
      }
    }
  }
  
  // beam profile
  double xBeam = theHodo->posX();
  double yBeam = theHodo->posY();

  // filling histos
  meETBxib_      -> Fill(xibEta, xibPhi);
  meETBhodoX_    -> Fill(xBeam);
  meETBhodoY_    -> Fill(yBeam);
  meETBampltdc_  -> Fill(theTDC->offset(),ampl1x1);
  meETBe1x1_     -> Fill(ampl1x1);
  meETBe3x3_     -> Fill(ampl3x3);
  meETBe5x5_     -> Fill(ampl5x5);
  meETBe1e9_     -> Fill(ampl1x1/ampl3x3);
  meETBe1e25_    -> Fill(ampl1x1/ampl5x5);
  meETBe9e25_    -> Fill(ampl3x3/ampl5x5);
  meETBe1vsX_    -> Fill(xBeam,ampl1x1);
  meETBe1vsY_    -> Fill(yBeam,ampl1x1);
  meETBe1e9vsX_  -> Fill(xBeam,ampl1x1/ampl3x3);
  meETBe1e9vsY_  -> Fill(yBeam,ampl1x1/ampl3x3);
  meETBe1e25vsX_ -> Fill(xBeam,ampl1x1/ampl5x5);
  meETBe1e25vsY_ -> Fill(yBeam,ampl1x1/ampl5x5);
  meETBe9e25vsX_ -> Fill(xBeam,ampl3x3/ampl5x5);
  meETBe9e25vsY_ -> Fill(yBeam,ampl3x3/ampl5x5);

  for(int ii=0; ii < 10; ++ii){ meETBShape_->Fill(double(ii)+theTDC->offset(),sampleSave[ii]); }

  if ( (fabs(xBeam)<2.5) && (fabs(yBeam)<2.5) ){ 
    meETBe1x1_center_  -> Fill(ampl1x1);
    meETBe3x3_center_  -> Fill(ampl3x3);
    meETBe5x5_center_  -> Fill(ampl5x5);
  }
}