#include <EEBeamCaloTask.h>

Inheritance diagram for EEBeamCaloTask:

Public Member Functions
	EEBeamCaloTask (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~EEBeamCaloTask ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...

void	beginJob (void)
	BeginJob. More...

void	beginRun (const edm::Run &r, const edm::EventSetup &c)
	BeginRun. More...

void	cleanup (void)
	Cleanup. More...

void	endJob (void)
	EndJob. More...

void	endRun (const edm::Run &r, const edm::EventSetup &c)
	EndRun. More...

void	reset (void)
	Reset. More...

void	setup (void)
	Setup. More...

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
bool	changed_cry_in_beam_

bool	changed_tb_status_

int	cib_ [12]

int	crystal_step_

MonitorElement *	CrystalInBeam_vs_Event_

MonitorElement *	CrystalsDone_

DQMStore *	dqmStore_

edm::InputTag	EBDigiCollection_

edm::InputTag	EcalRawDataCollection_

edm::InputTag	EcalTBEventHeader_

edm::InputTag	EcalUncalibratedRecHitCollection_

bool	enableCleanup_

int	event_last_reset_

int	evt_after_change_

int	ievt_

bool	init_

int	last_cry_in_beam_

int	lastStableStatus_

MonitorElement *	meBBCaloAllNeededCry_

MonitorElement *	meBBCaloCryOnBeam_

MonitorElement *	meBBCaloCryRead_

MonitorElement *	meBBCaloE3x3_

MonitorElement *	meBBCaloE3x3Moving_

MonitorElement *	meBBCaloEne_ [cryInArray_]

MonitorElement *	meBBCaloGains_ [cryInArray_]

MonitorElement *	meBBCaloMaxEneCry_

MonitorElement *	meBBCaloPulseProf_ [cryInArray_]

MonitorElement *	meBBCaloPulseProfG12_ [cryInArray_]

MonitorElement *	meBBNumCaloCryRead_

MonitorElement *	meEEBCaloBeamCentered_

MonitorElement *	meEEBCaloDesync_

MonitorElement *	meEEBCaloE1MaxCry_

MonitorElement *	meEEBCaloE1vsCry_

MonitorElement *	meEEBCaloE3x3vsCry_

MonitorElement *	meEEBCaloEntriesVsCry_

MonitorElement *	meEEBCaloReadCryErrors_

bool	mergeRuns_

std::string	prefixME_

int	previous_cry_in_beam_

int	previous_ev_num_

int	PreviousCrystalinBeam_ [3]

int	PreviousTableStatus_ [2]

bool	profileArranged_

int	table_step_

MonitorElement *	TableMoving_

bool	wasFakeChange_

Static Private Attributes
static const int	cryInArray_ = 9

static const int	defaultPede_ = 200

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 20 of file EEBeamCaloTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 33 of file EEBeamCaloTask.cc.

References cryInArray_, CrystalInBeam_vs_Event_, CrystalsDone_, dqmStore_, EBDigiCollection_, EcalRawDataCollection_, EcalTBEventHeader_, EcalUncalibratedRecHitCollection_, enableCleanup_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, init_, meBBCaloAllNeededCry_, meBBCaloCryOnBeam_, meBBCaloCryRead_, meBBCaloE3x3_, meBBCaloE3x3Moving_, meBBCaloEne_, meBBCaloGains_, meBBCaloMaxEneCry_, meBBCaloPulseProf_, meBBCaloPulseProfG12_, meBBNumCaloCryRead_, meEEBCaloBeamCentered_, meEEBCaloDesync_, meEEBCaloE1MaxCry_, meEEBCaloE1vsCry_, meEEBCaloE3x3vsCry_, meEEBCaloEntriesVsCry_, meEEBCaloReadCryErrors_, mergeRuns_, cmsCodeRules.cppFunctionSkipper::operator, prefixME_, and TableMoving_.

                                                        {
 
   init_ = false;
 
   dqmStore_ = edm::Service<DQMStore>().operator->();
 
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
 
   mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
 
   EcalTBEventHeader_ = ps.getParameter<edm::InputTag>("EcalTBEventHeader");
   EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
   EBDigiCollection_ = ps.getParameter<edm::InputTag>("EBDigiCollection");
   EcalUncalibratedRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalUncalibratedRecHitCollection");
 
   for (int i = 0; i < cryInArray_ ; i++) {
     meBBCaloPulseProf_[i]=0;
     meBBCaloPulseProfG12_[i]=0;
     meBBCaloGains_[i]=0;
     meBBCaloEne_[i]=0;
 
     //meBBCaloPulseProfMoving_[i]=0;
     //meBBCaloPulseProfG12Moving_[i]=0;
     //meBBCaloGainsMoving_[i]=0;
     //meBBCaloEneMoving_[i]=0;
   }
 
   meBBCaloCryRead_ = 0;
   //meBBCaloCryReadMoving_ = 0;
 
   meBBNumCaloCryRead_ = 0;
   meBBCaloAllNeededCry_ = 0;
 
   meBBCaloE3x3_ = 0;
   meBBCaloE3x3Moving_ = 0;
 
   meBBCaloCryOnBeam_ = 0;
   meBBCaloMaxEneCry_ = 0;
   TableMoving_ = 0;
 
   CrystalsDone_ = 0;
   CrystalInBeam_vs_Event_ = 0;
   meEEBCaloReadCryErrors_ = 0;
   meEEBCaloE1vsCry_ = 0;
   meEEBCaloE3x3vsCry_ = 0;
   meEEBCaloEntriesVsCry_ = 0;
   meEEBCaloBeamCentered_ = 0;
 
   meEEBCaloE1MaxCry_ = 0;
 //   for(int u=0;u<851;u++){
 //     meBBCaloE3x3Cry_[u]=0;
 //     meBBCaloE1Cry_[u]=0;
 //   }
 
   meEEBCaloDesync_ = 0;
 
 }

EEBeamCaloTask::~EEBeamCaloTask ( )

virtual

Destructor.

Definition at line 93 of file EEBeamCaloTask.cc.

                                {
 
 }

Member Function Documentation

void EEBeamCaloTask::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 413 of file EEBeamCaloTask.cc.

References abs, ecalMGPA::adc(), EcalMGPASample::adc(), EcalDCCHeaderBlock::BEAMH2, EcalDCCHeaderBlock::BEAMH4, newFWLiteAna::bin, ExpressReco_HICollisions_FallBack::cerr, changed_cry_in_beam_, changed_tb_status_, cryInArray_, crystal_step_, EcalTBEventHeader::crystalInBeam(), CrystalInBeam_vs_Event_, CrystalsDone_, createTree::dd, defaultPede_, EBDigiCollection_, EcalEndcap, EcalRawDataCollection_, EcalTBEventHeader_, EcalUncalibratedRecHitCollection_, event(), event_last_reset_, EcalTBEventHeader::eventNumber(), MonitorElement::Fill(), EcalMGPASample::gainId(), MonitorElement::getBinContent(), edm::Event::getByLabel(), i, EBDetId::ic(), ievt_, init_, last_cry_in_beam_, lastStableStatus_, LogDebug, meBBCaloAllNeededCry_, meBBCaloCryOnBeam_, meBBCaloCryRead_, meBBCaloE3x3_, meBBCaloE3x3Moving_, meBBCaloEne_, meBBCaloGains_, meBBCaloMaxEneCry_, meBBCaloPulseProf_, meBBCaloPulseProfG12_, meBBNumCaloCryRead_, meEEBCaloBeamCentered_, meEEBCaloDesync_, meEEBCaloE1MaxCry_, meEEBCaloE1vsCry_, meEEBCaloE3x3vsCry_, meEEBCaloEntriesVsCry_, meEEBCaloReadCryErrors_, previous_cry_in_beam_, previous_ev_num_, PreviousCrystalinBeam_, PreviousTableStatus_, edm::Handle< T >::product(), profileArranged_, MonitorElement::Reset(), EcalDataFrame::sample(), MonitorElement::setBinContent(), setup(), Numbers::subDet(), EcalTBEventHeader::syncError(), table_step_, EcalTBEventHeader::tableIsMoving(), TableMoving_, cond::rpcobtemp::temp, and wasFakeChange_.

                                                                      {
 
   bool enable = false;
 
   edm::Handle<EcalRawDataCollection> dcchs;
 
   if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
 
     for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
 
       if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
 
       if ( dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH4 ||
            dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH2 ) enable = true;
     }
 
   } else {
     edm::LogWarning("EEBeamCaloTask") << EcalRawDataCollection_ << " not available";
   }
 
   if ( ! enable ) return;
   if ( ! init_ ) this->setup();
   ievt_++;
 
   edm::Handle<EcalTBEventHeader> pEventHeader;
   const EcalTBEventHeader* evtHeader=0;
 
   if ( e.getByLabel(EcalTBEventHeader_, pEventHeader) ) {
     evtHeader = pEventHeader.product(); // get a ptr to the product
   } else {
     std::cerr << "Error! can't get the product for the event header" << std::endl;
   }
 
   //FIX ME, in the task we should use LV1 instead of ievt_ (prescaling)
   int cry_in_beam = 0;
   bool tb_moving = false;//just for test, to be filled with info from the event
   int event = 0;
 
   if(evtHeader){
     //cry_in_beam = evtHeader->nominalCrystalInBeam();
     cry_in_beam = evtHeader->crystalInBeam();
     tb_moving = evtHeader->tableIsMoving();
     event = evtHeader->eventNumber();
     if( evtHeader->syncError() ) {meEEBCaloDesync_->Fill(crystal_step_);}
   }
   else {
     cry_in_beam =   previous_cry_in_beam_;
     tb_moving = lastStableStatus_;
     event = previous_ev_num_ +10;
   }
 
   previous_cry_in_beam_ = cry_in_beam;
   previous_ev_num_ = event;
 
   //cry_in_beam = 702;//just for test, to be filled with info from the event
 
   bool reset_histos_stable = false;
   bool reset_histos_moving = false;
 
   bool skip_this_event = false;
 
 //   if(ievt_ > 500){tb_moving=true; }
 //   if(ievt_ > 1000){tb_moving=false; cry_in_beam = 703;}
 //   if(ievt_ > 2000){tb_moving=true; }
 //   if(ievt_ > 2500){tb_moving=false; cry_in_beam = 704;}
 //   if(ievt_ == 3000){cry_in_beam = 702;}
 //   if(ievt_ == 3001){cry_in_beam = 703;}
 //   if(ievt_ > 3500){tb_moving=true; }
 
 //   if(ievt_ > 3300){tb_moving=true; }
 //   if(ievt_ > 6100){tb_moving=false; cry_in_beam = 705;}
 //   if(ievt_ == 6201){tb_moving=true; }
 //   if(ievt_ > 9000){tb_moving=true; }
 //   if(ievt_ == 11021){tb_moving=false; }
 //   if(ievt_ > 12100){tb_moving=false; cry_in_beam = 706;}
 //   if(ievt_ > 15320){tb_moving=true; }
 //   if(ievt_ > 15500){tb_moving=false; cry_in_beam = 707;}
 
 
  //  //if(ievt_ > 20){tb_moving=true; }
 //   //  if(ievt_ == 23){tb_moving=true; }
 //   if(ievt_ > 50){tb_moving=false; cry_in_beam = 705;}
 //   //if(ievt_ > 90 ){tb_moving=true; }
 //   if(ievt_ == 65){tb_moving=false; }
 //   if(ievt_ > 110){tb_moving=false; cry_in_beam = 706;}
 //   if(ievt_ == 116){cry_in_beam = 709;}
 //   //if(ievt_ > 115){tb_moving=true; }
 //   if(ievt_ > 150){tb_moving=false; cry_in_beam = 707;}
 
   //   #include "DQM/EcalEndcapMonitorTasks/interface/cry_in_beam_run_ecs73214.h"
 
   if(ievt_ < 3){last_cry_in_beam_ = cry_in_beam;}
 
   if(tb_moving){
 
     TableMoving_->Fill(1);
     if( PreviousTableStatus_[0] == 0 &&  PreviousTableStatus_[1] == 1 && lastStableStatus_ == 0){
       reset_histos_moving=true;
       wasFakeChange_ = false;
       // ! Warning! This works in the assumption that the crystal in beam stay the same
       // while the tb is moving and is set to the new one only when the table
       // reaches the new position
       lastStableStatus_ = 1;
 
     }
     else if( PreviousTableStatus_[1] == 0) {
       skip_this_event=true;
       changed_tb_status_ = true;
       wasFakeChange_ = true;
     }
     // just skip the first event when the table change status
     PreviousTableStatus_[0] = PreviousTableStatus_[1];
     PreviousTableStatus_[1] = 1;
   }//end of if(tb_moving)
 
   else {// table is not moving
 
     TableMoving_->Fill(0);
     if( PreviousTableStatus_[0] == 1 &&  PreviousTableStatus_[1] == 0 && lastStableStatus_ == 1){
       //reset_histos_stable = true;
       wasFakeChange_ = false;
       lastStableStatus_ = 0;
     }
     else if(PreviousTableStatus_[1] == 1) {
       skip_this_event=true;
       changed_tb_status_ = true;
       wasFakeChange_ = true;
     }
     // just skip the first event when the table change status
     PreviousTableStatus_[0]=PreviousTableStatus_[1];
     PreviousTableStatus_[1]=0;
 
     // check also whether cry in beam  has changed
     if(  PreviousCrystalinBeam_[0] == PreviousCrystalinBeam_[1]  &&   PreviousCrystalinBeam_[1] != PreviousCrystalinBeam_[2] && PreviousCrystalinBeam_[2] == cry_in_beam ){
       reset_histos_stable=true;
       wasFakeChange_ = false;
     }
     else if (PreviousCrystalinBeam_[2] != cry_in_beam){
       changed_cry_in_beam_ = true;
       skip_this_event=true;
       wasFakeChange_ = true;
     }
     // }
 
     PreviousCrystalinBeam_[0] = PreviousCrystalinBeam_[1];
     PreviousCrystalinBeam_[1] = PreviousCrystalinBeam_[2];
     PreviousCrystalinBeam_[2] =  cry_in_beam;
   }
   //if (! changed_tb_status_ && ! changed_cry_in_beam_ ){// standard data taking
 
     if( !tb_moving ) {CrystalInBeam_vs_Event_->Fill(event,float(cry_in_beam));}
     else{CrystalInBeam_vs_Event_->Fill(event,-100); }
     if ( !profileArranged_ ){
       float dd=0;
       int mbin =0;
       for( int bin=1; bin < 20001; bin++ ){
     float temp = CrystalInBeam_vs_Event_->getBinContent(bin);
     if(temp>0){ dd= temp+0.01; mbin=bin; break;}
       }
       if(mbin >0) { CrystalInBeam_vs_Event_->Fill(20*mbin-1,dd);}
       profileArranged_ = true;
     }
 
     // }
 //   else{ // here there is either a step or a daq error
 //     // keep 10 events in a buffer waiting to decide for the step or the error
 //     if(tb_moving){cib_[evt_after_change_]=-100;}
 //     else {cib_[evt_after_change_]=cry_in_beam;}
 
 //     if(evt_after_change_ >= 9){
 //       evt_after_change_ =0;
 //       if(wasFakeChange_){// here is an error: put the 10 events in the profile
 //  for(int u=0; u<10; u++){
 //    CrystalInBeam_vs_Event_->Fill(ievt_-9+u , cib_[u]);
 //  }
 //       }
 //       //for a real change just skip the first 10 events after change
 //       changed_tb_status_=false;
 //       changed_cry_in_beam_ = false;
 //     }
 //     else{evt_after_change_ ++;}
 //   }
 
   if(reset_histos_moving){
     edm::LogInfo("EEBeamCaloTask") << "event " << ievt_ << " resetting histos for moving table!! ";
 
     //     meEEBCaloE1vsCry_->setBinContent(crystal_step_ , meBBCaloEne_[4]->getMean() );
     //     meEEBCaloE1vsCry_->setBinError(crystal_step_ , meBBCaloEne_[4]->getRMS() );
     //     meEEBCaloE3x3vsCry_->setBinContent(crystal_step_ , meBBCaloE3x3_->getMean() );
     //     meEEBCaloE3x3vsCry_->setBinError(crystal_step_ , meBBCaloE3x3_->getRMS() );
 
     //     meEEBCaloEntriesVsCry_->setBinContent(crystal_step_ ,  meBBCaloE3x3_->getEntries() );
 
     table_step_++;
 
 
     //here the follwowing histos should be reset
     // for (int u=0;u<cryInArray_;u++){
     //  meBBCaloPulseProfMoving_[u]->Reset();
     //  meBBCaloPulseProfG12Moving_[u]->Reset();
     //  meBBCaloGainsMoving_[u]->Reset();
     //  meBBCaloEneMoving_[u]->Reset();
     // }
     // meBBCaloCryReadMoving_->Reset();
     meBBCaloE3x3Moving_->Reset();
 
   }
 
 
   if(reset_histos_stable){
     if( event - event_last_reset_ > 30){//to be tuned, to avoid a double reset for the change in the table status and
                                         //in the crystal in beam. This works ONLY if the crystal in beam stay the same
                                         // while the table is moving.
                                         //One can also think to remove the reset of the histograms when the table change
                                         // status from moving to stable, and to leave the reset only if the cry_in_beam changes.
 
       edm::LogInfo("EEBeamCaloTask") << "event " << ievt_ << " resetting histos for stable table!! ";
 
       //       meEEBCaloE1vsCry_->setBinContent(crystal_step_ , meBBCaloEne_[4]->getMean() );
       //       meEEBCaloE1vsCry_->setBinError(crystal_step_ , meBBCaloEne_[4]->getRMS() );
       //       meEEBCaloE3x3vsCry_->setBinContent(crystal_step_ , meBBCaloE3x3_->getMean() );
       //       meEEBCaloE3x3vsCry_->setBinError(crystal_step_ , meBBCaloE3x3_->getRMS() );
       //       //
       //       meEEBCaloEntriesVsCry_->setBinContent(crystal_step_ ,  meBBCaloE3x3_->getEntries() );
 
       event_last_reset_ = event;
 
       last_cry_in_beam_ = cry_in_beam;
       crystal_step_++;
 
       //here the follwowing histos should be reset
       for (int u=0;u<cryInArray_;u++){
     meBBCaloPulseProf_[u]->Reset();
     meBBCaloPulseProfG12_[u]->Reset();
     meBBCaloGains_[u]->Reset();
     meBBCaloEne_[u]->Reset();
       }
       meBBCaloCryRead_->Reset();
       meBBCaloE3x3_->Reset();
       meEEBCaloBeamCentered_->Reset();
     }
   }
 
  if(skip_this_event){
    edm::LogInfo("EEBeamCaloTask") << "event " << event <<" analyzed: "<<ievt_ << " : skipping this event!! ";
    return;}
 
  // now CrystalsDone_ contains the crystal on beam at the beginning fo a new step, and not when it has finished !!
  // <5 just to avoid that we skip the event just after the reset and we do not set CrystalsDone_ .
  // if( ievt_ - event_last_reset_ < 5){ CrystalsDone_->setBinContent(cry_in_beam , crystal_step_ );}
  CrystalsDone_->setBinContent(cry_in_beam , crystal_step_ );
   int eta_c = ( cry_in_beam-1)/20 ;
   int phi_c = ( cry_in_beam-1)%20 ;
 
   float xie = eta_c + 0.5;
   float xip = phi_c + 0.5;
   if (!tb_moving) {meBBCaloCryOnBeam_->Fill(xie,xip);}
 
   edm::Handle<EBDigiCollection> digis;
   e.getByLabel(EBDigiCollection_, digis);
   int nebd = digis->size();
 
   meBBNumCaloCryRead_->Fill(nebd);
 
   //matrix 7x7 around cry in beam
   int cry_to_beRead[49]; //0 or -1 for non existing crystals (eg 852)
   for(int u=0;u<49;u++){cry_to_beRead[u]=0;}
   // chech that all the crystals in the 7x7 exist
   for(int de=-3; de<4; de++){
     for(int dp=-3; dp<4; dp++){
       //int cry_num = (phi_c+dp) + 20*(eta_c+de) +1;
       int u = de -7*dp + 24;
       bool existing_cry = (phi_c+dp) >= 0 && (phi_c+dp) <= 19 && (eta_c+de) >=0 && (eta_c+de) <= 84;
       if(!existing_cry){cry_to_beRead[u]=-1;}
     }
   }
 
 
   meEEBCaloEntriesVsCry_->Fill(cry_in_beam);
 
   for ( EBDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr ) {
 
     EBDetId id = digiItr->id();
 
     int ic = id.ic();
     int ie = (ic-1)/20;
     int ip = (ic-1)%20;
 
     int deta_c= ie - eta_c;
     int dphi_c= ip - phi_c;
     if (! tb_moving){meBBCaloCryRead_->Fill(deta_c, dphi_c);}
     //else {meBBCaloCryReadMoving_->Fill(deta_c, dphi_c);}
 
     if(std::abs(deta_c) > 3 || std::abs(dphi_c) > 3){continue;}
     int i_toBeRead = deta_c -7*dphi_c + 24;
     if( i_toBeRead > -1 &&  i_toBeRead <49){
       cry_to_beRead[i_toBeRead]++;
       //if( (ievt_ == 4000 || ievt_ == 13000 || ievt_ == 13002 ) &&   i_toBeRead == 5){ cry_to_beRead[i_toBeRead] -=1;}
     }
 
     if(std::abs(deta_c) > 1 || std::abs(dphi_c) > 1){continue;}
     int i_in_array = deta_c -3*dphi_c + 4;
 
     if( i_in_array < 0 || i_in_array > 8 ){continue;}
 
     EBDataFrame dataframe = (*digiItr);
 
     for (int i = 0; i < 10; i++) {
       int adc = dataframe.sample(i).adc();
       int gainid = dataframe.sample(i).gainId();
       //if( (ievt_ == 15400 || ievt_ == 15600 || ievt_ == 15700 ) &&   i_in_array == 4 && i == 4){ gainid =2;}
       //if( (ievt_ == 15400 || ievt_ == 15600 || ievt_ == 15700 ) &&   i_in_array == 6 && i == 6){ gainid =3;}
 
       if ( gainid == 1 ){// gain 12
     if(! tb_moving){
       meBBCaloPulseProfG12_[i_in_array]->Fill(i,float(adc));
       meBBCaloPulseProf_[i_in_array]->Fill(i,float(adc));
       meBBCaloGains_[i_in_array]->Fill(12);
     }
     //else{
     //  meBBCaloPulseProfG12Moving_[i_in_array]->Fill(i,float(adc));
     //  meBBCaloPulseProfMoving_[i_in_array]->Fill(i,float(adc));
     //  meBBCaloGainsMoving_[i_in_array]->Fill(12);
     //}
       }
       else if ( gainid == 2 ){// gain 6
     float val = (float(adc)-defaultPede_)*2 + defaultPede_;
     if(! tb_moving){
       meBBCaloPulseProf_[i_in_array]->Fill(i,val);
       meBBCaloGains_[i_in_array]->Fill(6);
     }
     //else{
     //  meBBCaloPulseProfMoving_[i_in_array]->Fill(i,val);
     //  meBBCaloGainsMoving_[i_in_array]->Fill(6);
     //}
       }
       else if ( gainid == 3 ){// gain 1
     float val = (float(adc)-defaultPede_)*12 + defaultPede_;
     if(! tb_moving){
     meBBCaloPulseProf_[i_in_array]->Fill(i,val);
     meBBCaloGains_[i_in_array]->Fill(1);
     }
     //else{
     //meBBCaloPulseProfMoving_[i_in_array]->Fill(i,val);
     //meBBCaloGainsMoving_[i_in_array]->Fill(1);
     //}
       }
     }// end of loop over samples
   }// end of loop over digis
 
   //now  if everything was correct cry_to_beRead should be filled with 1 or -1 but not 0
   bool all_cry_readout = true;
 
   // if( ievt_ == 4000 || ievt_ == 13000 || ievt_ == 13002 ) {all_cry_readout = false;}
   if(all_cry_readout){ meBBCaloAllNeededCry_->Fill(1.5);}//bin3
   else {
     meBBCaloAllNeededCry_->Fill(-0.5);//bin1
     if( tb_moving ) {meEEBCaloReadCryErrors_->Fill( crystal_step_+0.5 );}
     else {meEEBCaloReadCryErrors_->Fill( crystal_step_ );}
   }
 
   //the part involving rechits
 
   edm::Handle<EcalUncalibratedRecHitCollection> hits;
   e.getByLabel(EcalUncalibratedRecHitCollection_, hits);
   int neh = hits->size();
   LogDebug("EEBeamCaloTask") << "event " << event <<" analyzed: "<< ievt_ << " hits collection size " << neh;
   float ene3x3=0;
   float maxEne = 0;
   int ieM =-1, ipM = -1;//for the crystal with maximum energy deposition
   float cryInBeamEne =0;
   for ( EcalUncalibratedRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr ) {
 
     EBDetId id = hitItr->id();
 
     int ic = id.ic();
     int ie = (ic-1)/20;
     int ip = (ic-1)%20;
 
     int deta_c= ie - eta_c;
     int dphi_c= ip - phi_c;
 
     int i_in_array = deta_c -3*dphi_c + 4;
 
     float R_ene = hitItr->amplitude();
     if ( R_ene <= 0. ) R_ene = 0.0;
     if(R_ene > maxEne){
       maxEne=R_ene;
       ieM =ie; ipM = ip;
     }
     if(std::abs(deta_c) > 1 || std::abs(dphi_c) > 1){continue;}
     meEEBCaloBeamCentered_->Fill(deta_c,dphi_c,R_ene);
 
     if( i_in_array < 0 || i_in_array > 8 ){continue;}
 
     if(i_in_array == 4){cryInBeamEne = R_ene;}
     if(! tb_moving){meBBCaloEne_[i_in_array]->Fill(R_ene);}
     //else{meBBCaloEneMoving_[i_in_array]->Fill(R_ene);}
     ene3x3 += R_ene;
 
   }//end of loop over rechits
 
   if (!tb_moving){
     meBBCaloE3x3_->Fill(ene3x3);
     meEEBCaloE1vsCry_->Fill(cry_in_beam , cryInBeamEne );
     meEEBCaloE3x3vsCry_->Fill(cry_in_beam, ene3x3 );
     //     if( cry_in_beam > 0 && cry_in_beam < 851){
     //       meBBCaloE3x3Cry_[cry_in_beam]->Fill(ene3x3);
     //       meBBCaloE1Cry_[cry_in_beam]->Fill(cryInBeamEne);
     //     }
     meBBCaloMaxEneCry_->Fill(ieM,ipM);
     meEEBCaloE1MaxCry_->Fill(maxEne);
   }
   else{meBBCaloE3x3Moving_->Fill(ene3x3);}
 }

void EEBeamCaloTask::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 97 of file EEBeamCaloTask.cc.

References dqmStore_, ievt_, prefixME_, profileArranged_, DQMStore::rmdir(), and DQMStore::setCurrentFolder().

                                  {
 
   ievt_ = 0;
 
   profileArranged_ = false;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EEBeamCaloTask");
     dqmStore_->rmdir(prefixME_ + "/EEBeamCaloTask");
   }
 
 }

void EEBeamCaloTask::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 110 of file EEBeamCaloTask.cc.

References Numbers::initGeometry(), mergeRuns_, and reset().

                                                                      {
 
   Numbers::initGeometry(c, false);
 
   if ( ! mergeRuns_ ) this->reset();
 
 }

void EEBeamCaloTask::cleanup ( void )

protected

Cleanup.

Definition at line 328 of file EEBeamCaloTask.cc.

References cryInArray_, CrystalInBeam_vs_Event_, CrystalsDone_, dqmStore_, edm::getName(), MonitorElement::getName(), i, init_, meBBCaloAllNeededCry_, meBBCaloCryOnBeam_, meBBCaloCryRead_, meBBCaloE3x3_, meBBCaloE3x3Moving_, meBBCaloEne_, meBBCaloGains_, meBBCaloMaxEneCry_, meBBCaloPulseProf_, meBBCaloPulseProfG12_, meBBNumCaloCryRead_, meEEBCaloBeamCentered_, meEEBCaloDesync_, meEEBCaloE1MaxCry_, meEEBCaloE1vsCry_, meEEBCaloE3x3vsCry_, meEEBCaloEntriesVsCry_, meEEBCaloReadCryErrors_, prefixME_, DQMStore::removeElement(), DQMStore::setCurrentFolder(), and TableMoving_.

Referenced by endJob().

                                 {
 
   if ( ! init_ ) return;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EEBeamCaloTask");
     for (int i = 0; i < cryInArray_ ; i++) {
       if ( meBBCaloPulseProf_[i] ) dqmStore_->removeElement( meBBCaloPulseProf_[i]->getName() );
       meBBCaloPulseProf_[i] = 0;
       if ( meBBCaloPulseProfG12_[i] ) dqmStore_->removeElement( meBBCaloPulseProfG12_[i]->getName() );
       meBBCaloPulseProfG12_[i] = 0;
       if ( meBBCaloGains_[i] ) dqmStore_->removeElement( meBBCaloGains_[i]->getName() );
       meBBCaloGains_[i] = 0;
       if ( meBBCaloEne_[i] ) dqmStore_->removeElement( meBBCaloEne_[i]->getName() );
       meBBCaloEne_[i] = 0;
 
 //       if ( meBBCaloPulseProfMoving_[i] ) dqmStore_->removeElement( meBBCaloPulseProfMoving_[i]->getName() );
 //       meBBCaloPulseProfMoving_[i] = 0;
 //       if ( meBBCaloPulseProfG12Moving_[i] ) dqmStore_->removeElement( meBBCaloPulseProfG12Moving_[i]->getName() );
 //       meBBCaloPulseProfG12Moving_[i] = 0;
 //       if ( meBBCaloGainsMoving_[i] ) dqmStore_->removeElement( meBBCaloGainsMoving_[i]->getName() );
 //       meBBCaloGainsMoving_[i] = 0;
 //       if ( meBBCaloEneMoving_[i] ) dqmStore_->removeElement( meBBCaloEneMoving_[i]->getName() );
 //       meBBCaloEneMoving_[i] = 0;
     }
 
 //     dqmStore_->setCurrentFolder(prefixME_ + "/EEBeamCaloTask/EnergyHistos");
 //     for(int u=0; u< 851;u++){
 //       if ( meBBCaloE3x3Cry_[u] ) dqmStore_->removeElement( meBBCaloE3x3Cry_[u]->getName() );
 //       meBBCaloE3x3Cry_[u] = 0;
 //       if ( meBBCaloE1Cry_[u] ) dqmStore_->removeElement( meBBCaloE1Cry_[u]->getName() );
 //       meBBCaloE1Cry_[u] = 0;
 //     }
 
 //     dqmStore_->setCurrentFolder(prefixME_ + "/EEBeamCaloTask");
     if ( meBBCaloCryRead_ ) dqmStore_->removeElement( meBBCaloCryRead_->getName() );
     meBBCaloCryRead_ = 0;
 //    if ( meBBCaloCryReadMoving_ ) dqmStore_->removeElement( meBBCaloCryReadMoving_->getName() );
 //    meBBCaloCryReadMoving_ = 0;
     if ( meBBCaloAllNeededCry_ ) dqmStore_->removeElement( meBBCaloAllNeededCry_->getName() );
     meBBCaloAllNeededCry_ = 0;
     if ( meBBNumCaloCryRead_ ) dqmStore_->removeElement( meBBNumCaloCryRead_->getName() );
     meBBNumCaloCryRead_ = 0;
     if ( meBBCaloE3x3_ ) dqmStore_->removeElement( meBBCaloE3x3_->getName() );
     meBBCaloE3x3_ = 0;
     if ( meBBCaloE3x3Moving_ ) dqmStore_->removeElement( meBBCaloE3x3Moving_->getName() );
     meBBCaloE3x3Moving_ = 0;
     if ( meBBCaloCryOnBeam_ ) dqmStore_->removeElement( meBBCaloCryOnBeam_->getName() );
     meBBCaloCryOnBeam_ = 0;
     if ( meBBCaloMaxEneCry_ ) dqmStore_->removeElement( meBBCaloMaxEneCry_->getName() );
     meBBCaloMaxEneCry_ = 0;
     if ( TableMoving_ ) dqmStore_->removeElement( TableMoving_->getName() );
     TableMoving_ = 0;
     if ( CrystalsDone_ ) dqmStore_->removeElement( CrystalsDone_->getName() );
     CrystalsDone_ = 0;
     if ( CrystalInBeam_vs_Event_ ) dqmStore_->removeElement( CrystalInBeam_vs_Event_->getName() );
     CrystalInBeam_vs_Event_ = 0;
     if( meEEBCaloReadCryErrors_ ) dqmStore_->removeElement( meEEBCaloReadCryErrors_->getName() );
     meEEBCaloReadCryErrors_ = 0;
     if( meEEBCaloE1vsCry_ ) dqmStore_->removeElement( meEEBCaloE1vsCry_->getName() );
     meEEBCaloE1vsCry_ = 0;
     if( meEEBCaloE3x3vsCry_ ) dqmStore_->removeElement( meEEBCaloE3x3vsCry_->getName() );
     meEEBCaloE3x3vsCry_ = 0;
     if( meEEBCaloEntriesVsCry_ )  dqmStore_->removeElement( meEEBCaloEntriesVsCry_->getName() );
     meEEBCaloEntriesVsCry_ = 0;
     if( meEEBCaloBeamCentered_ ) dqmStore_->removeElement( meEEBCaloBeamCentered_->getName() );
     meEEBCaloBeamCentered_ = 0;
     if( meEEBCaloE1MaxCry_ ) dqmStore_->removeElement(meEEBCaloE1MaxCry_->getName() );
     meEEBCaloE1MaxCry_ = 0;
     if( meEEBCaloDesync_ ) dqmStore_->removeElement(meEEBCaloDesync_->getName() );
     meEEBCaloDesync_ = 0;
   }
 
   init_ = false;
 
 }

void EEBeamCaloTask::endJob ( void )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 405 of file EEBeamCaloTask.cc.

References cleanup(), enableCleanup_, and ievt_.

                                {
 
   edm::LogInfo("EEBeamCaloTask") << "analyzed " << ievt_ << " events";
 
   if ( enableCleanup_ ) this->cleanup();
 
 }

void EEBeamCaloTask::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

EndRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 118 of file EEBeamCaloTask.cc.

                                                                    {
 
 }

void EEBeamCaloTask::reset ( void )

protected

Reset.

Definition at line 122 of file EEBeamCaloTask.cc.

References cryInArray_, CrystalInBeam_vs_Event_, CrystalsDone_, i, meBBCaloAllNeededCry_, meBBCaloCryOnBeam_, meBBCaloCryRead_, meBBCaloE3x3_, meBBCaloE3x3Moving_, meBBCaloEne_, meBBCaloGains_, meBBCaloMaxEneCry_, meBBCaloPulseProf_, meBBCaloPulseProfG12_, meBBNumCaloCryRead_, meEEBCaloBeamCentered_, meEEBCaloDesync_, meEEBCaloE1MaxCry_, meEEBCaloE1vsCry_, meEEBCaloE3x3vsCry_, meEEBCaloEntriesVsCry_, meEEBCaloReadCryErrors_, MonitorElement::Reset(), and TableMoving_.

Referenced by beginRun().

                                {
 
     for (int i = 0; i < cryInArray_ ; i++) {
       if ( meBBCaloPulseProf_[i] ) meBBCaloPulseProf_[i]->Reset();
       if ( meBBCaloPulseProfG12_[i] ) meBBCaloPulseProfG12_[i]->Reset();
       if ( meBBCaloGains_[i] ) meBBCaloGains_[i]->Reset();
       if ( meBBCaloEne_[i] ) meBBCaloEne_[i]->Reset();
 
 //       if ( meBBCaloPulseProfMoving_[i] ) meBBCaloPulseProfMoving_[i]->Reset();
 //       if ( meBBCaloPulseProfG12Moving_[i] ) meBBCaloPulseProfG12Moving_[i]->Reset();
 //       if ( meBBCaloGainsMoving_[i] ) meBBCaloGainsMoving_[i]->Reset();
 //       if ( meBBCaloEneMoving_[i] ) meBBCaloEneMoving_[i]->Reset();
     }
 
 //     for(int u=0; u< 1701;u++){
 //       if ( meBBCaloE3x3Cry_[u] ) meBBCaloE3x3Cry_[u]->Reset();
 //       if ( meBBCaloE1Cry_[u] ) meBBCaloE1Cry_[u]->Reset();
 //     }
 
     if ( meBBCaloCryRead_ ) meBBCaloCryRead_->Reset();
 //    if ( meBBCaloCryReadMoving_ ) meBBCaloCryReadMoving_->Reset();
     if ( meBBCaloAllNeededCry_ ) meBBCaloAllNeededCry_->Reset();
     if ( meBBNumCaloCryRead_ ) meBBNumCaloCryRead_->Reset();
     if ( meBBCaloE3x3_ ) meBBCaloE3x3_->Reset();
     if ( meBBCaloE3x3Moving_ ) meBBCaloE3x3Moving_->Reset();
     if ( meBBCaloCryOnBeam_ ) meBBCaloCryOnBeam_->Reset();
     if ( meBBCaloMaxEneCry_ ) meBBCaloMaxEneCry_->Reset();
     if ( TableMoving_ ) TableMoving_->Reset();
     if ( CrystalsDone_ ) CrystalsDone_->Reset();
     if ( CrystalInBeam_vs_Event_ ) CrystalInBeam_vs_Event_->Reset();
     if( meEEBCaloReadCryErrors_ ) meEEBCaloReadCryErrors_->Reset();
     if( meEEBCaloE1vsCry_ ) meEEBCaloE1vsCry_->Reset();
     if( meEEBCaloE3x3vsCry_ ) meEEBCaloE3x3vsCry_->Reset();
     if( meEEBCaloEntriesVsCry_ )  meEEBCaloEntriesVsCry_->Reset();
     if( meEEBCaloBeamCentered_ ) meEEBCaloBeamCentered_->Reset();
     if( meEEBCaloE1MaxCry_ ) meEEBCaloE1MaxCry_->Reset();
     if( meEEBCaloDesync_ ) meEEBCaloDesync_->Reset();
 
 }

void EEBeamCaloTask::setup ( void )

protected

Setup.

Definition at line 162 of file EEBeamCaloTask.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), changed_cry_in_beam_, changed_tb_status_, cib_, cryInArray_, crystal_step_, CrystalInBeam_vs_Event_, CrystalsDone_, dqmStore_, event_last_reset_, evt_after_change_, trackerHits::histo, i, init_, last_cry_in_beam_, lastStableStatus_, meBBCaloAllNeededCry_, meBBCaloCryOnBeam_, meBBCaloCryRead_, meBBCaloE3x3_, meBBCaloE3x3Moving_, meBBCaloEne_, meBBCaloGains_, meBBCaloMaxEneCry_, meBBCaloPulseProf_, meBBCaloPulseProfG12_, meBBNumCaloCryRead_, meEEBCaloBeamCentered_, meEEBCaloDesync_, meEEBCaloE1MaxCry_, meEEBCaloE1vsCry_, meEEBCaloE3x3vsCry_, meEEBCaloEntriesVsCry_, meEEBCaloReadCryErrors_, prefixME_, previous_cry_in_beam_, PreviousCrystalinBeam_, PreviousTableStatus_, profileArranged_, MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), table_step_, TableMoving_, and wasFakeChange_.

Referenced by analyze().

                               {
 
   init_ = true;
   profileArranged_= false;
   char histo[200];
 
   PreviousTableStatus_[0]=0;//let's start with stable...
   PreviousTableStatus_[1]=0;//let's start with stable...
 
   PreviousCrystalinBeam_[0] = 0;
   PreviousCrystalinBeam_[1] = 0;
   PreviousCrystalinBeam_[2] = -1;
   // PreviousCrystalinBeam_[2] = -1 is needed to have a correct step vs cry matching
   lastStableStatus_=0;
   for(int u=0;u<10;u++){cib_[u]=0;}
   changed_tb_status_= false;
   changed_cry_in_beam_ = false;
   evt_after_change_ =0;
   wasFakeChange_= false;
   table_step_=1;
   crystal_step_=1;
   event_last_reset_ = 0;
   last_cry_in_beam_ = 0;
   previous_cry_in_beam_ = 1;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EEBeamCaloTask");
 
     for (int i = 0; i < cryInArray_ ; i++) {
       sprintf(histo, "EEBCT pulse profile cry %01d", i+1);
       //considering the gain the range is 4096*12 ~ 50000
       meBBCaloPulseProf_[i] = dqmStore_->bookProfile(histo, histo, 10,0.,10.,50000,0.,50000.,"s");
 
       sprintf(histo, "EEBCT pulse profile in G12 cry %01d", i+1);
       meBBCaloPulseProfG12_[i] = dqmStore_->bookProfile(histo, histo, 10,0.,10.,4096,0.,4096.,"s");
       meBBCaloPulseProfG12_[i]->setAxisTitle("#sample", 1);
       meBBCaloPulseProfG12_[i]->setAxisTitle("ADC", 2);
 
       sprintf(histo, "EEBCT found gains cry %01d", i+1);
       meBBCaloGains_[i] =  dqmStore_->book1D(histo,histo,14,0.,14.);
       meBBCaloGains_[i]->setAxisTitle("gain", 1);
       // g1-> bin 2, g6-> bin 7, g12-> bin 13
 
       sprintf(histo, "EEBCT rec energy cry %01d", i+1);
       meBBCaloEne_[i] =  dqmStore_->book1D(histo,histo,500,0.,9000.);
       meBBCaloEne_[i]->setAxisTitle("rec ene (ADC)", 1);
       //9000 ADC in G12 equivalent is about 330 GeV
 
 
 //       sprintf(histo, "EEBCT pulse profile moving table cry %01d", i+1);
 //       //considering the gain the range is 4096*12 ~ 50000
 //       meBBCaloPulseProfMoving_[i] = dqmStore_->bookProfile(histo, histo, 10,0.,10.,50000,0.,50000.,"s");
 
 //       sprintf(histo, "EEBCT pulse profile in G12 moving table cry %01d", i+1);
 //       meBBCaloPulseProfG12Moving_[i] = dqmStore_->bookProfile(histo, histo, 10,0.,10.,4096,0.,4096.,"s");
 
 //       sprintf(histo, "EEBCT found gains moving table cry %01d", i+1);
 //       meBBCaloGainsMoving_[i] =  dqmStore_->book1D(histo,histo,14,0.,14.);
 //       // g1-> bin 2, g6-> bin 7, g12-> bin 13
 
 //       sprintf(histo, "EEBCT rec energy moving table cry %01d", i+1);
 //       meBBCaloEneMoving_[i] =  dqmStore_->book1D(histo,histo,2000,0.,9000.);
 //       //9000 ADC in G12 equivalent is about 330 GeV
 
     }
 
 //     dqmStore_->setCurrentFolder(prefixME_ + "/EEBeamCaloTask/EnergyHistos");
 //     for(int u=0; u< 851;u++){
 //       sprintf(histo, "EEBCT rec Ene sum 3x3 cry: %04d",u);
 //       meBBCaloE3x3Cry_[u] = dqmStore_->book1D(histo,histo,1000,0.,4500.);
 
 //       sprintf(histo, "EEBCT rec Energy1 cry: %04d",u);
 //       meBBCaloE1Cry_[u] = dqmStore_->book1D(histo,histo,1000,0.,4500.);
 //     }
 
 //     dqmStore_->setCurrentFolder(prefixME_ + "/EEBeamCaloTask");
     sprintf(histo, "EEBCT readout crystals");
     meBBCaloCryRead_ =  dqmStore_->book2D(histo,histo,9,-4.,5.,9,-4.,5.);
     //matrix of readout crystal around cry in beam
 
     //sprintf(histo, "EEBCT readout crystals table moving");
     //meBBCaloCryReadMoving_ =  dqmStore_->book2D(histo,histo,9,-4.,5.,9,-4.,5.);
     //matrix of readout crystal around cry in beam
 
     sprintf(histo, "EEBCT all needed crystals readout");
     meBBCaloAllNeededCry_ = dqmStore_->book1D(histo,histo,3,-1.,2.);
     // not all needed cry are readout-> bin 1, all needed cry are readout-> bin 3
 
     sprintf(histo, "EEBCT readout crystals number");
     meBBNumCaloCryRead_ = dqmStore_->book1D(histo,histo,851,0.,851.);
     meBBNumCaloCryRead_->setAxisTitle("number of read crystals", 1);
 
     sprintf(histo, "EEBCT rec Ene sum 3x3");
     meBBCaloE3x3_ = dqmStore_->book1D(histo,histo,500,0.,9000.);
     meBBCaloE3x3_->setAxisTitle("rec ene (ADC)", 1);
     //9000 ADC in G12 equivalent is about 330 GeV
 
     sprintf(histo, "EEBCT rec Ene sum 3x3 table moving");
     meBBCaloE3x3Moving_ = dqmStore_->book1D(histo,histo,500,0.,9000.);
     //9000 ADC in G12 equivalent is about 330 GeV
 
     sprintf(histo, "EEBCT crystal on beam");
     meBBCaloCryOnBeam_ = dqmStore_->book2D(histo, histo, 85, 0., 85., 20, 0., 20.);
 
     sprintf(histo, "EEBCT crystal with maximum rec energy");
     meBBCaloMaxEneCry_ = dqmStore_->book2D(histo, histo, 85, 0., 85., 20, 0., 20.);
 
     sprintf(histo, "EEBCT table is moving");
     TableMoving_ = dqmStore_->book1D(histo,histo,2,0.,1.1);
     TableMoving_->setAxisTitle("table status (0=stable, 1=moving)", 1);
     //table is moving-> bin 2, table is not moving-> bin 1
 
     sprintf(histo, "EEBCT crystals done");
     CrystalsDone_ = dqmStore_->book1D(histo,histo,850,1.,851.);
     CrystalsDone_->setAxisTitle("crystal", 1);
     CrystalsDone_->setAxisTitle("step in the scan", 2);
     //for a crystal done the corresponing bin is filled with the step in the
     //autoscan pertainig to the given crystales
 
     sprintf(histo, "EEBCT crystal in beam vs event");
     CrystalInBeam_vs_Event_ = dqmStore_->bookProfile(histo, histo, 20000,0.,400000.,1802,-101.,851.,"s");
     CrystalInBeam_vs_Event_->setAxisTitle("event", 1);
     CrystalInBeam_vs_Event_->setAxisTitle("crystal in beam", 2);
     // 1 bin each 20 events
     // when table is moving for just one events fill with -100
 
     sprintf(histo, "EEBCT readout crystals errors");
     meEEBCaloReadCryErrors_ = dqmStore_->book1D(histo, histo, 425,1.,86.);
     meEEBCaloReadCryErrors_->setAxisTitle("step in the scan", 1);
 
     sprintf(histo, "EEBCT average rec energy in the single crystal");
     //meEEBCaloE1vsCry_ = dqmStore_->book1D(histo, histo, 85,1.,86.);
     meEEBCaloE1vsCry_ = dqmStore_->bookProfile(histo, histo, 850,1.,851.,500,0.,9000.,"s");
     meEEBCaloE1vsCry_->setAxisTitle("crystal", 1);
     meEEBCaloE1vsCry_->setAxisTitle("rec energy (ADC)", 2);
 
     sprintf(histo, "EEBCT average rec energy in the 3x3 array");
     //meEEBCaloE3x3vsCry_= dqmStore_->book1D(histo, histo,85,1.,86.);
     meEEBCaloE3x3vsCry_ = dqmStore_->bookProfile(histo, histo, 850,1.,851.,500,0.,9000.,"s");
     meEEBCaloE3x3vsCry_->setAxisTitle("crystal", 1);
     meEEBCaloE3x3vsCry_->setAxisTitle("rec energy (ADC)", 2);
 
     sprintf(histo, "EEBCT number of entries");
     meEEBCaloEntriesVsCry_ = dqmStore_->book1D(histo, histo,850,1.,851.);
     meEEBCaloEntriesVsCry_->setAxisTitle("crystal", 1);
     meEEBCaloEntriesVsCry_->setAxisTitle("number of events (prescaled)", 2);
 
     sprintf(histo, "EEBCT energy deposition in the 3x3");
     meEEBCaloBeamCentered_ = dqmStore_->book2D(histo, histo,3,-1.5,1.5,3,-1.5,1.5);
     meEEBCaloBeamCentered_->setAxisTitle("\\Delta \\eta", 1);
     meEEBCaloBeamCentered_->setAxisTitle("\\Delta \\phi", 2);
 
     sprintf(histo, "EEBCT E1 in the max cry");
     meEEBCaloE1MaxCry_= dqmStore_->book1D(histo,histo,500,0.,9000.);
     meEEBCaloE1MaxCry_->setAxisTitle("rec Ene (ADC)", 1);
 
     sprintf(histo, "EEBCT Desynchronization vs step");
     meEEBCaloDesync_= dqmStore_->book1D(histo, histo, 85 ,1.,86.);
     meEEBCaloDesync_->setAxisTitle("step", 1);
     meEEBCaloDesync_->setAxisTitle("Desynchronized events", 2);
 
   }
 
 }