CastorEventDisplay.cc

Go to the documentation of this file.

#include "DQM/CastorMonitor/interface/CastorEventDisplay.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

//***************************************************//
//********** CastorEventDisplay: *******************//
//********** Author: Dmytro Volyanskyy   ************//
//********** Date  : 04.05.2010 (first version) ******// 
//***************************************************//

//==================================================================//
//======================= Constructor ==============================//
//==================================================================//
CastorEventDisplay::CastorEventDisplay() {
  ievt_=0;
}

//==================================================================//
//======================= Destructor ==============================//
//==================================================================//
CastorEventDisplay::~CastorEventDisplay(){
}

void CastorEventDisplay::reset(){
}


//==========================================================//
//========================= setup ==========================//
//==========================================================//

void CastorEventDisplay::setup(const edm::ParameterSet& ps, DQMStore* dbe){
  
  CastorBaseMonitor::setup(ps,dbe);
  baseFolder_ = rootFolder_+"CastorEventDisplay";
  

  offline_             = ps.getUntrackedParameter<bool>("OfflineMode", false); 

  if(fVerbosity>0) std::cout << "CastorEventDisplay::setup (start)" << std::endl;
  
  ievt_=0;  allEnergyEvent=0.; maxEnergyEvent=0.;
  X_pos=0.; Y_pos=0.; Z_pos=0.;
  X_pos_maxE=0.; Y_pos_maxE=0.; Z_pos_maxE=0.;

  if ( m_dbe !=NULL ) {    
  m_dbe->setCurrentFolder(baseFolder_);

  meEVT_ = m_dbe->bookInt("EventDisplay Event Number"); // meEVT_->Fill(ievt_);
  
  meCastor3Dhits = m_dbe->book3D("CASTOR 3D hits- cumulative", "CASTOR 3D hits - cumulative", 30, 1420, 1600, 35, -35, 35, 35, -35, 35);


  if( offline_ ){
  //-- swap z and y axis 
  TH3F* Castor3Dhits = meCastor3Dhits->getTH3F();
  Castor3Dhits->GetXaxis()->SetTitle("Z [cm]"); //-- also swap x and z
  Castor3Dhits->GetYaxis()->SetTitle("X [cm]");
  Castor3Dhits->GetZaxis()->SetTitle("Y [cm]");
  }

  meCastor3DhitsMaxEnergy = m_dbe->book3D("CASTOR 3D hits- event with the largest deposited E", "CASTOR 3D hits- event with the largest deposited E",  30, 1420, 1600, 20, -30, 30, 20, -30, 30); //-- swap z and y axis

  if( offline_ ){
  TH3F* Castor3DhitsMaxEnergy = meCastor3DhitsMaxEnergy->getTH3F();
  Castor3DhitsMaxEnergy->GetXaxis()->SetTitle("Z [cm]"); //-- also swap x and z
  Castor3DhitsMaxEnergy->GetYaxis()->SetTitle("X [cm]");
  Castor3DhitsMaxEnergy->GetZaxis()->SetTitle("Y [cm]");
  Castor3DhitsMaxEnergy->SetDrawOption("LEGO2");
  }

 } 

  else{
  if(fVerbosity>0) std::cout << "CastorEventDisplay::setup - NO DQMStore service" << std::endl; 
  }

  if(fVerbosity>0) std::cout << "CastorEventDisplay::setup (end)" << std::endl;
  
  return;
}

//==========================================================//
//================== processEvent ==========================//
//==========================================================//

void CastorEventDisplay::processEvent(const CastorRecHitCollection& castorHits, const CaloGeometry& caloGeometry ){
  
  if(fVerbosity>0) std::cout << "==>CastorEventDisplay::processEvent !!!"<< std::endl;

   meEVT_->Fill(ievt_);

  if(!m_dbe) { 
    if(fVerbosity>0) std::cout <<"CastorEventDisplay::processEvent => DQMStore is not instantiated !!!"<<std::endl;  
    return; 
  }

   allEnergyEvent=0.;


  CastorRecHitCollection::const_iterator CASTORiter;

     if(castorHits.size()>0)
    {    
       if(fVerbosity>0) std::cout << "==>CastorEventDisplay::processEvent: castorHits.size()>0 !!!" << std::endl; 

       for (CASTORiter=castorHits.begin(); CASTORiter!=castorHits.end(); ++CASTORiter) { 
        
         energy = CASTORiter->energy();    

        HcalCastorDetId CastorID(CASTORiter->detid().rawId());

    //const CaloSubdetectorGeometry* subdetectorGeometry=caloGeometry.getSubdetectorGeometry(CastorID) ; 
        //const CaloCellGeometry* cellGeometry =  subdetectorGeometry->getGeometry(CastorID) ;
        //X_pos = cellGeometry->getPosition().x() ; Y_pos = cellGeometry->getPosition().y() ; 
        //Z_pos = cellGeometry->getPosition().z() ;
       
        X_pos = caloGeometry.getSubdetectorGeometry(CastorID)->getGeometry(CastorID)->getPosition().x() ;
        Y_pos = caloGeometry.getSubdetectorGeometry(CastorID)->getGeometry(CastorID)->getPosition().y() ;
        Z_pos = caloGeometry.getSubdetectorGeometry(CastorID)->getGeometry(CastorID)->getPosition().z() ;
        
        if (energy>0){
         allEnergyEvent= allEnergyEvent+energy ;
         meCastor3Dhits->Fill(std::abs(Z_pos),X_pos,Y_pos, energy);
    }

    if(fVerbosity>0)  std::cout<<"ENERGY="<< energy <<" X_pos="<< X_pos <<" Y_pos="<< Y_pos <<" Z_pos="<< Z_pos << std::endl;   
      } //-- end of for loop


    if(fVerbosity>0)  std::cout<<"TOTAL ENERGY in an event="<< allEnergyEvent << std::endl;



    if (allEnergyEvent > maxEnergyEvent) { 
     maxEnergyEvent = allEnergyEvent;
     meCastor3DhitsMaxEnergy->Reset();
     if(fVerbosity>0) std::cout<<"LARGEST ENERGY in an event="<< maxEnergyEvent << std::endl;
 
     for (CASTORiter=castorHits.begin(); CASTORiter!=castorHits.end(); ++CASTORiter){
       HcalCastorDetId CastorID_maxE(CASTORiter->detid().rawId());
       X_pos_maxE = caloGeometry.getSubdetectorGeometry(CastorID_maxE)->getGeometry(CastorID_maxE)->getPosition().x() ;
       Y_pos_maxE = caloGeometry.getSubdetectorGeometry(CastorID_maxE)->getGeometry(CastorID_maxE)->getPosition().y() ;
       Z_pos_maxE = caloGeometry.getSubdetectorGeometry(CastorID_maxE)->getGeometry(CastorID_maxE)->getPosition().z() ;
       meCastor3DhitsMaxEnergy->Fill(std::abs(Z_pos_maxE),X_pos_maxE,Y_pos_maxE, CASTORiter->energy() );
      }
     meCastor3DhitsMaxEnergy->getTH3F()->SetOption("BOX");
     meCastor3DhitsMaxEnergy->update();
    }
  } 

  else { if(fVerbosity>0) std::cout<<"CastorEventDisplay::processEvent NO Castor RecHits !!!"<<std::endl; }


  if (showTiming) { 
      cpu_timer.stop(); std::cout << " TIMER::CastorEventDisplay -> " << cpu_timer.cpuTime() << std::endl; 
      cpu_timer.reset(); cpu_timer.start();  
    }
    
  ievt_++; 
  
  return;
}