TrackerHitAnalyzer Class Reference

#include <TrackerHitAnalyzer.h>

Inheritance diagram for TrackerHitAnalyzer:

Public Member Functions
	TrackerHitAnalyzer (const edm::ParameterSet &ps)
	Constructor. More...
	~TrackerHitAnalyzer ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...
void	beginRun (edm::Run const &, edm::EventSetup const &)
	Begin Run. More...
void	endJob (void)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Attributes
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_pxlBrlHigh_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_pxlBrlLow_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_pxlFwdHigh_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_pxlFwdLow_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTECHigh_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTECLow_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTIBHigh_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTIBLow_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTIDHigh_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTIDLow_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTOBHigh_Token_
edm::EDGetTokenT < edm::PSimHitContainer >	edmPSimHitContainer_siTOBLow_Token_
edm::EDGetTokenT < edm::SimTrackContainer >	edmSimTrackContainerToken_
DQMStore *	fDBE
std::string	fOutputFile
MonitorElement *	h1e [12]
MonitorElement *	h1ex [12]
MonitorElement *	h1ey [12]
MonitorElement *	h1ez [12]
MonitorElement *	h1lx [12]
MonitorElement *	h1ly [12]
MonitorElement *	h2e [12]
MonitorElement *	h2ex [12]
MonitorElement *	h2ey [12]
MonitorElement *	h2ez [12]
MonitorElement *	h2lx [12]
MonitorElement *	h2ly [12]
MonitorElement *	h3e [12]
MonitorElement *	h3ex [12]
MonitorElement *	h3ey [12]
MonitorElement *	h3ez [12]
MonitorElement *	h3lx [12]
MonitorElement *	h3ly [12]
MonitorElement *	h4e [12]
MonitorElement *	h4ex [12]
MonitorElement *	h4ey [12]
MonitorElement *	h4ez [12]
MonitorElement *	h4lx [12]
MonitorElement *	h4ly [12]
MonitorElement *	h5e [12]
MonitorElement *	h5ex [12]
MonitorElement *	h5ey [12]
MonitorElement *	h5ez [12]
MonitorElement *	h5lx [12]
MonitorElement *	h5ly [12]
MonitorElement *	h6e [12]
MonitorElement *	h6ex [12]
MonitorElement *	h6ey [12]
MonitorElement *	h6ez [12]
MonitorElement *	h6lx [12]
MonitorElement *	h6ly [12]
MonitorElement *	htofeta
MonitorElement *	htofeta_profile
MonitorElement *	htofphi
MonitorElement *	htofphi_profile
MonitorElement *	htofr
MonitorElement *	htofr_profile
MonitorElement *	htofz
MonitorElement *	htofz_profile
bool	verbose_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 22 of file TrackerHitAnalyzer.h.

Constructor & Destructor Documentation

TrackerHitAnalyzer::TrackerHitAnalyzer

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 43 of file TrackerHitAnalyzer.cc.

  : verbose_( ps.getUntrackedParameter<bool>( "Verbosity",false ) )
  , edmPSimHitContainer_pxlBrlLow_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "PxlBrlLowSrc" ) ) )
  , edmPSimHitContainer_pxlBrlHigh_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "PxlBrlHighSrc" ) ) )
  , edmPSimHitContainer_pxlFwdLow_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "PxlFwdLowSrc" ) ) )
  , edmPSimHitContainer_pxlFwdHigh_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "PxlFwdHighSrc" ) ) )
  , edmPSimHitContainer_siTIBLow_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTIBLowSrc" ) ) )
  , edmPSimHitContainer_siTIBHigh_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTIBHighSrc" ) ) )
  , edmPSimHitContainer_siTOBLow_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTOBLowSrc" ) ) )
  , edmPSimHitContainer_siTOBHigh_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTOBHighSrc" ) ) )
  , edmPSimHitContainer_siTIDLow_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTIDLowSrc" ) ) )
  , edmPSimHitContainer_siTIDHigh_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTIDHighSrc" ) ) )
  , edmPSimHitContainer_siTECLow_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTECLowSrc" ) ) )
  , edmPSimHitContainer_siTECHigh_Token_( consumes<edm::PSimHitContainer>( ps.getParameter<edm::InputTag>( "SiTECHighSrc" ) ) )
  , edmSimTrackContainerToken_( consumes<edm::SimTrackContainer>( ps.getParameter<edm::InputTag>( "G4TrkSrc" ) ) )
  , fDBE( NULL )
  , fOutputFile( ps.getUntrackedParameter<std::string>( "outputFile", "TrackerHitHisto.root" ) ) {
}

TrackerHitAnalyzer::~TrackerHitAnalyzer

(

)

Destructor.

Definition at line 314 of file TrackerHitAnalyzer.cc.

{
   // don't try to delete any pointers - they're handled by DQM machinery
}

Member Function Documentation

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

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 348 of file TrackerHitAnalyzer.cc.

References Geom::Phi< T >::degrees(), cond::rpcobgas::detid, edmPSimHitContainer_pxlBrlHigh_Token_, edmPSimHitContainer_pxlBrlLow_Token_, edmPSimHitContainer_pxlFwdHigh_Token_, edmPSimHitContainer_pxlFwdLow_Token_, edmPSimHitContainer_siTECHigh_Token_, edmPSimHitContainer_siTECLow_Token_, edmPSimHitContainer_siTIBHigh_Token_, edmPSimHitContainer_siTIBLow_Token_, edmPSimHitContainer_siTIDHigh_Token_, edmPSimHitContainer_siTIDLow_Token_, edmPSimHitContainer_siTOBHigh_Token_, edmPSimHitContainer_siTOBLow_Token_, edmSimTrackContainerToken_, PV3DBase< T, PVType, FrameType >::eta(), eta(), edm::EventID::event(), MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), h1e, h1ex, h1ey, h1ez, h1lx, h1ly, h2e, h2ex, h2ey, h2ez, h2lx, h2ly, h3e, h3ex, h3ey, h3ez, h3lx, h3ly, h4e, h4ex, h4ey, h4ez, h4lx, h4ly, h5e, h5ex, h5ey, h5ez, h5lx, h5ly, h6e, h6ex, h6ey, h6ez, h6lx, h6ly, htofeta, htofphi, htofr, htofz, edm::EventBase::id(), edm::HandleBase::isValid(), create_public_lumi_plots::log, PV3DBase< T, PVType, FrameType >::mag(), siStripFEDMonitor_P5_cff::Max, siStripFEDMonitor_P5_cff::Min, AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::phi(), reco::return(), edm::EventID::run(), mathSSE::sqrt(), funct::tan(), theta(), patCandidatesForDimuonsSequences_cff::tracker, and PV3DBase< T, PVType, FrameType >::z().

{

  edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
   
  // iterator to access containers
  edm::PSimHitContainer::const_iterator itHit;
  // get Pixel Barrel information
  // extract low container
  edm::Handle<edm::PSimHitContainer> PxlBrlLowContainer;
  e.getByToken( edmPSimHitContainer_pxlBrlLow_Token_, PxlBrlLowContainer );
  if (!PxlBrlLowContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze")
      << "Unable to find TrackerHitsPixelBarrelLowTof in event!";
    return;
  }  
  // extract high container
  edm::Handle<edm::PSimHitContainer> PxlBrlHighContainer;
  e.getByToken( edmPSimHitContainer_pxlBrlHigh_Token_, PxlBrlHighContainer );
  if (!PxlBrlHighContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze")
      << "Unable to find TrackerHitsPixelBarrelHighTof in event!";
    return;
  }
  // get Pixel Forward information
  // extract low container
  edm::Handle<edm::PSimHitContainer> PxlFwdLowContainer;
  e.getByToken( edmPSimHitContainer_pxlFwdLow_Token_, PxlFwdLowContainer );
  if (!PxlFwdLowContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze")
      << "Unable to find TrackerHitsPixelEndcapLowTof in event!";
    return;
  }
  // extract high container
  edm::Handle<edm::PSimHitContainer> PxlFwdHighContainer;
  e.getByToken( edmPSimHitContainer_pxlFwdHigh_Token_, PxlFwdHighContainer );
  if (!PxlFwdHighContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze")
      << "Unable to find TrackerHitsPixelEndcapHighTof in event!";
    return;
  }
  
  // get Silicon TIB information
  // extract TIB low container
  edm::Handle<edm::PSimHitContainer> SiTIBLowContainer;
  e.getByToken( edmPSimHitContainer_siTIBLow_Token_, SiTIBLowContainer );
  if (!SiTIBLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTIBLowTof in event!";
    return;
  }
  // extract TIB high container
  edm::Handle<edm::PSimHitContainer> SiTIBHighContainer;
  e.getByToken( edmPSimHitContainer_siTIBHigh_Token_, SiTIBHighContainer );
  if (!SiTIBHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTIBHighTof in event!";
    return;
  }
  // get Silicon TOB information
  // extract TOB low container
  edm::Handle<edm::PSimHitContainer> SiTOBLowContainer;
  e.getByToken( edmPSimHitContainer_siTOBLow_Token_, SiTOBLowContainer );
  if (!SiTOBLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTOBLowTof in event!";
    return;
  }
  // extract TOB high container
  edm::Handle<edm::PSimHitContainer> SiTOBHighContainer;
  e.getByToken( edmPSimHitContainer_siTOBHigh_Token_, SiTOBHighContainer );
  if (!SiTOBHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTOBHighTof in event!";
    return;
  }
  
  // get Silicon TID information
  // extract TID low container
  edm::Handle<edm::PSimHitContainer> SiTIDLowContainer;
  e.getByToken( edmPSimHitContainer_siTIDLow_Token_, SiTIDLowContainer );
  if (!SiTIDLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTIDLowTof in event!";
    return;
  }
  // extract TID high container
  edm::Handle<edm::PSimHitContainer> SiTIDHighContainer;
  e.getByToken( edmPSimHitContainer_siTIDHigh_Token_, SiTIDHighContainer );
  if (!SiTIDHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTIDHighTof in event!";
    return;
  }
  // get Silicon TEC information
  // extract TEC low container
  edm::Handle<edm::PSimHitContainer> SiTECLowContainer;
  e.getByToken( edmPSimHitContainer_siTECLow_Token_, SiTECLowContainer );
  if (!SiTECLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTECLowTof in event!";
    return;
  }
  // extract TEC high container
  edm::Handle<edm::PSimHitContainer> SiTECHighContainer;
  e.getByToken( edmPSimHitContainer_siTECHigh_Token_, SiTECHighContainer );
  if (!SiTECHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze")
      << "Unable to find TrackerHitsTECHighTof in event!";
    return;
  }
  
  // get G4Track information
  
  edm::Handle<edm::SimTrackContainer> G4TrkContainer;
  e.getByToken( edmSimTrackContainerToken_, G4TrkContainer );
  if (!G4TrkContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze")
      << "Unable to find SimTrack in event!";
    return;
  }

  // Get geometry information

  edm::ESHandle<TrackerGeometry> tracker;
  c.get<TrackerDigiGeometryRecord>().get( tracker );


  
  int ir = -100;
  edm::SimTrackContainer::const_iterator itTrk;
  for (itTrk = G4TrkContainer->begin(); itTrk != G4TrkContainer->end(); 
       ++itTrk) {

//    std::cout << "itTrk = "<< itTrk << std::endl;
    double eta =0, p =0;
    const CLHEP::HepLorentzVector& G4Trk = CLHEP::HepLorentzVector(itTrk->momentum().x(),
                                                     itTrk->momentum().y(),
                                                     itTrk->momentum().z(),
                                                     itTrk->momentum().e());  
    p =sqrt(G4Trk[0]*G4Trk[0]+G4Trk[1]*G4Trk[1]+G4Trk[2]*G4Trk[2]);
      if ( p == 0) 
          edm::LogError("TrackerHitAnalyzer::analyze") 
          << "TrackerTest::INFO: Primary has p = 0 ";
      else {
          double costheta  = G4Trk[2]/p;
          double theta = acos(TMath::Min(TMath::Max(costheta, -1.),1.));
          eta = -log(tan(theta/2));          
          
          if (eta>0.0 && eta<=0.5) ir = 0;
          if (eta>0.5 && eta<=1.0) ir = 1;
          if (eta>1.0 && eta<=1.5) ir = 2;
          if (eta>1.5 && eta<=2.0) ir = 3;
          if (eta>2.0 && eta<=2.5) ir = 4;
          if (eta>2.5) ir = 5;
      
          if (eta>-0.5 && eta<= 0.0) ir = 6;
          if (eta>-1.0 && eta<=-0.5) ir = 7;
          if (eta>-1.5 && eta<=-1.0) ir = 8;
          if (eta>-2.0 && eta<=-1.5) ir = 9;
          if (eta>-2.5 && eta<=-2.0) ir = 10;
          if (eta<=-2.5) ir = 11;
//          edm::LogInfo("EventInfo") << " eta = " << eta << " ir = " << ir;
//    std::cout << " " << std::endl;
//          std::cout << "eta " << eta << " ir = " << ir << std::endl;                  
//    std::cout << " " << std::endl;
      }
  }   
  // get Pixel information
  for (itHit = PxlBrlLowContainer->begin(); itHit != PxlBrlLowContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h5e[ir]->Fill(itHit->energyLoss());
    h5ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
    h5ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
    h5ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
    h5lx[ir]->Fill(itHit->localPosition().x());
    h5ly[ir]->Fill(itHit->localPosition().y());
  }
  for (itHit = PxlBrlHighContainer->begin(); itHit != PxlBrlHighContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());


   h5e[ir]->Fill(itHit->energyLoss());
   h5ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h5ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h5ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h5lx[ir]->Fill(itHit->localPosition().x());
   h5ly[ir]->Fill(itHit->localPosition().y());
  }
  for (itHit = PxlFwdLowContainer->begin(); itHit != PxlFwdLowContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());


   h6e[ir]->Fill(itHit->energyLoss());
   h6ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h6ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h6ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h6lx[ir]->Fill(itHit->localPosition().x());
   h6ly[ir]->Fill(itHit->localPosition().y());
  }  
  for (itHit = PxlFwdHighContainer->begin(); itHit != PxlFwdHighContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());


   h6e[ir]->Fill(itHit->energyLoss());
   h6ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h6ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h6ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h6lx[ir]->Fill(itHit->localPosition().x());
   h6ly[ir]->Fill(itHit->localPosition().y());
  }
  // get TIB information
  for (itHit = SiTIBLowContainer->begin(); itHit != SiTIBLowContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());


   h1e[ir]->Fill(itHit->energyLoss());
   h1ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h1ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h1ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h1lx[ir]->Fill(itHit->localPosition().x());
   h1ly[ir]->Fill(itHit->localPosition().y());
  }
  for (itHit = SiTIBHighContainer->begin(); itHit != SiTIBHighContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());


   h1e[ir]->Fill(itHit->energyLoss());
   h1ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h1ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h1ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h1lx[ir]->Fill(itHit->localPosition().x());
   h1ly[ir]->Fill(itHit->localPosition().y());
  }
  // get TOB information
  for (itHit = SiTOBLowContainer->begin(); itHit != SiTOBLowContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());


   h2e[ir]->Fill(itHit->energyLoss());
   h2ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h2ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h2ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h2lx[ir]->Fill(itHit->localPosition().x());
   h2ly[ir]->Fill(itHit->localPosition().y());
  }  
  for (itHit = SiTOBHighContainer->begin(); itHit != SiTOBHighContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

   h2e[ir]->Fill(itHit->energyLoss());
   h2ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h2ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h2ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h2lx[ir]->Fill(itHit->localPosition().x());
   h2ly[ir]->Fill(itHit->localPosition().y());
  }
  // get TID information
  for (itHit = SiTIDLowContainer->begin(); itHit != SiTIDLowContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

   h3e[ir]->Fill(itHit->energyLoss());
   h3ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h3ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h3ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h3lx[ir]->Fill(itHit->localPosition().x());
   h3ly[ir]->Fill(itHit->localPosition().y());
  }  
  for (itHit = SiTIDHighContainer->begin(); itHit != SiTIDHighContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

   h3e[ir]->Fill(itHit->energyLoss());
   h3ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h3ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h3ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h3lx[ir]->Fill(itHit->localPosition().x());
   h3ly[ir]->Fill(itHit->localPosition().y());
  }
  // get TEC information
  for (itHit = SiTECLowContainer->begin(); itHit != SiTECLowContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

   h4e[ir]->Fill(itHit->energyLoss());
   h4ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h4ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h4ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h4lx[ir]->Fill(itHit->localPosition().x());
   h4ly[ir]->Fill(itHit->localPosition().y());
  }  
  for (itHit = SiTECHighContainer->begin(); itHit != SiTECHighContainer->end(); ++itHit) {
    DetId detid=DetId(itHit->detUnitId());
    const GeomDetUnit * det=(const GeomDetUnit*)tracker->idToDetUnit( detid );
    GlobalPoint gpos=det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

   h4e[ir]->Fill(itHit->energyLoss());
   h4ex[ir]->Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
   h4ey[ir]->Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
   h4ez[ir]->Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
   h4lx[ir]->Fill(itHit->localPosition().x());
   h4ly[ir]->Fill(itHit->localPosition().y());
  }
  
      
   return ;
   
}

void TrackerHitAnalyzer::beginRun ( edm::Run const &  , edm::EventSetup const &    )

protectedvirtual

Begin Run.

Reimplemented from edm::EDAnalyzer.

Definition at line 62 of file TrackerHitAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), fDBE, h1e, h1ex, h1ey, h1ez, h1lx, h1ly, h2e, h2ex, h2ey, h2ez, h2lx, h2ly, h3e, h3ex, h3ey, h3ez, h3lx, h3ly, h4e, h4ex, h4ey, h4ez, h4lx, h4ly, h5e, h5ex, h5ey, h5ez, h5lx, h5ly, h6e, h6ex, h6ey, h6ez, h6lx, h6ly, htofeta, htofphi, htofr, htofz, i, NULL, cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), DQMStore::showDirStructure(), and verbose_.

                                                                       {
  fDBE = edm::Service<DQMStore>().operator->();
    
  Char_t  hname1[50], htitle1[80];
  Char_t  hname2[50], htitle2[80];
  Char_t  hname3[50], htitle3[80];
  Char_t  hname4[50], htitle4[80];
  Char_t  hname5[50], htitle5[80];
  Char_t  hname6[50], htitle6[80];
   
  if ( fDBE ) {
     if ( verbose_ ) {
       fDBE->setVerbose(1);
     } else {
       fDBE->setVerbose(0);
     }
  }
                                                                        
  if ( fDBE) {
    if ( verbose_ ) fDBE->showDirStructure();
  } 

  if ( fDBE != NULL ) {
//   fDBE->setCurrentFolder("TrackerHitsV/TrackerHitTask");
     
     // is there any way to record CPU Info ???
     // if so, it can be done once - via beginJob() 
 int nbin = 5000;   

    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/");
    htofeta  = fDBE->book2D ("tof_eta", "Time of flight vs eta", nbin , -3.0 , 3.0,200,-100,100);
    htofphi  = fDBE->book2D("tof_phi", "Time of flight vs phi", nbin,-180,180,200,-100,100);
    htofr  = fDBE->book2D("tof_r", "Time of flight vs r", nbin , 0 , 300, 200, -100,100);
    htofz  = fDBE->book2D("tof_z", "Time of flight vs z", nbin , -280 , 280, 200, -100,100);


 const float E2NEL = 1.; 
 
 const char *Region[] = {"005","051","115","152","225","253",
                       "-050","-105","-151","-215","-252","-325"};  
  nbin = 10000;   
      
// Energy loss histograms
   for(int i=0; i<12; i++) {
  
    sprintf (htitle1,"Energy loss in TIB %s", Region[i]);
    sprintf (htitle2,"Energy loss in TOB %s", Region[i]);
    sprintf (htitle3,"Energy loss in TID %s", Region[i]);
    sprintf (htitle4,"Energy loss in TEC %s", Region[i]);
    sprintf (htitle5,"Energy loss in BPIX %s", Region[i]);
    sprintf (htitle6,"Energy loss in FPIX %s", Region[i]);
    
    sprintf (hname1,"Eloss_TIB_%i",i+1);
    sprintf (hname2,"Eloss_TOB_%i",i+1);
    sprintf (hname3,"Eloss_TID_%i",i+1);
    sprintf (hname4,"Eloss_TEC_%i",i+1);
    sprintf (hname5,"Eloss_BPIX_%i",i+1);
    sprintf (hname6,"Eloss_FPIX_%i",i+1);
   
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
    h1e[i]  = fDBE->book1D (hname1, htitle1, nbin , 0.0 , 0.001*E2NEL);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
    h2e[i]  = fDBE->book1D (hname2, htitle2, nbin , 0.0 , 0.001*E2NEL);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
    h3e[i]  = fDBE->book1D (hname3, htitle3, nbin , 0.0 , 0.001*E2NEL);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
    h4e[i]  = fDBE->book1D (hname4, htitle4, nbin , 0.0 , 0.001*E2NEL);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
    h5e[i]  = fDBE->book1D (hname5, htitle5, nbin , 0.0 , 0.001*E2NEL);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
    h6e[i]  = fDBE->book1D (hname6, htitle6, nbin , 0.0 , 0.001*E2NEL);
   
   }

// limits
const float high[] = {0.03, 0.03, 0.02, 0.03, 0.03, 0.03};
const float low[] = {-0.03, -0.03, -0.02, -0.03, -0.03, -0.03};
   
   for(int i=0; i<12; i++) {
  
    sprintf (htitle1,"Entryx-Exitx in TIB %s", Region[i]);
    sprintf (htitle2,"Entryx-Exitx in TOB %s", Region[i]);
    sprintf (htitle3,"Entryx-Exitx in TID %s", Region[i]);
    sprintf (htitle4,"Entryx-Exitx in TEC %s", Region[i]);
    sprintf (htitle5,"Entryx-Exitx in BPIX %s", Region[i]);
    sprintf (htitle6,"Entryx-Exitx in FPIX %s", Region[i]);
    
    sprintf (hname1,"Entryx-Exitx_TIB_%i",i+1);
    sprintf (hname2,"Entryx-Exitx_TOB_%i",i+1);
    sprintf (hname3,"Entryx-Exitx_TID_%i",i+1);
    sprintf (hname4,"Entryx-Exitx_TEC_%i",i+1);
    sprintf (hname5,"Entryx-Exitx_BPIX_%i",i+1);
    sprintf (hname6,"Entryx-Exitx_FPIX_%i",i+1);
   
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
    h1ex[i]  = fDBE->book1D (hname1, htitle1, nbin , low[0] , high[0]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
    h2ex[i]  = fDBE->book1D (hname2, htitle2, nbin , low[1] , high[1]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
    h3ex[i]  = fDBE->book1D (hname3, htitle3, nbin , low[2] , high[2]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
    h4ex[i]  = fDBE->book1D (hname4, htitle4, nbin , low[3] , high[3]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
    h5ex[i]  = fDBE->book1D (hname5, htitle5, nbin , low[4] , high[4]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
    h6ex[i]  = fDBE->book1D (hname6, htitle6, nbin , low[5] , high[5]);
   
   }

const float high0[] = {0.05, 0.06, 0.03, 0.03, 0.03, 0.03};
const float low0[] = {-0.05, -0.06, -0.03, -0.03, -0.03, -0.03};

   for(int i=0; i<12; i++) {
  
    sprintf (htitle1,"Entryy-Exity in TIB %s", Region[i]);
    sprintf (htitle2,"Entryy-Exity in TOB %s", Region[i]);
    sprintf (htitle3,"Entryy-Exity in TID %s", Region[i]);
    sprintf (htitle4,"Entryy-Exity in TEC %s", Region[i]);
    sprintf (htitle5,"Entryy-Exity in BPIX %s", Region[i]);
    sprintf (htitle6,"Entryy-Exity in FPIX %s", Region[i]);
    
    sprintf (hname1,"Entryy-Exity_TIB_%i",i+1);
    sprintf (hname2,"Entryy-Exity_TOB_%i",i+1);
    sprintf (hname3,"Entryy-Exity_TID_%i",i+1);
    sprintf (hname4,"Entryy-Exity_TEC_%i",i+1);
    sprintf (hname5,"Entryy-Exity_BPIX_%i",i+1);
    sprintf (hname6,"Entryy-Exity_FPIX_%i",i+1);
   
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
    h1ey[i]  = fDBE->book1D (hname1, htitle1, nbin , low0[0] , high0[0]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
    h2ey[i]  = fDBE->book1D (hname2, htitle2, nbin , low0[1] , high0[1]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
    h3ey[i]  = fDBE->book1D (hname3, htitle3, nbin , low0[2] , high0[2]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
    h4ey[i]  = fDBE->book1D (hname4, htitle4, nbin , low0[3] , high0[3]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
    h5ey[i]  = fDBE->book1D (hname5, htitle5, nbin , low0[4] , high0[4]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
    h6ey[i]  = fDBE->book1D (hname6, htitle6, nbin , low0[5] , high0[5]);
   
   }

const float high1[] = {0.05, 0.06, 0.05, 0.06, 0.05, 0.05};
const float low1[]  = {0.,0.,0.,0.,0.,0.};

  for(int i=0; i<12; i++) {
  
    sprintf (htitle1,"abs(Entryz-Exitz) in TIB %s", Region[i]);
    sprintf (htitle2,"abs(Entryz-Exitz) in TOB %s", Region[i]);
    sprintf (htitle3,"abs(Entryz-Exitz) in TID %s", Region[i]);
    sprintf (htitle4,"abs(Entryz-Exitz) in TEC %s", Region[i]);
    sprintf (htitle5,"abs(Entryz-Exitz) in BPIX %s", Region[i]);
    sprintf (htitle6,"abs(Entryz-Exitz) in FPIX %s", Region[i]);
    
    sprintf (hname1,"Entryz-Exitz_TIB_%i",i+1);
    sprintf (hname2,"Entryz-Exitz_TOB_%i",i+1);
    sprintf (hname3,"Entryz-Exitz_TID_%i",i+1);
    sprintf (hname4,"Entryz-Exitz_TEC_%i",i+1);
    sprintf (hname5,"Entryz-Exitz_BPIX_%i",i+1);
    sprintf (hname6,"Entryz-Exitz_FPIX_%i",i+1);
   
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
    h1ez[i]  = fDBE->book1D (hname1, htitle1, nbin , low1[0] , high1[0]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
    h2ez[i]  = fDBE->book1D (hname2, htitle2, nbin , low1[1] , high1[1]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
    h3ez[i]  = fDBE->book1D (hname3, htitle3, nbin , low1[2] , high1[2]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
    h4ez[i]  = fDBE->book1D (hname4, htitle4, nbin , low1[3] , high1[3]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
    h5ez[i]  = fDBE->book1D (hname5, htitle5, nbin , low1[4] , high1[4]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
    h6ez[i]  = fDBE->book1D (hname6, htitle6, nbin , low1[5] , high1[5]);
   
   }


const float high2[] = {3.2, 5.0, 5.5, 6.2, 0.85, 0.5};   
const float low2[] = {-3.2, -5.0, -5.5, -6.2, -0.85, -0.5};   

   for(int i=0; i<12; i++) {
  
    sprintf (htitle1,"Localx in TIB %s", Region[i]);
    sprintf (htitle2,"Localx in TOB %s", Region[i]);
    sprintf (htitle3,"Localx in TID %s", Region[i]);
    sprintf (htitle4,"Localx in TEC %s", Region[i]);
    sprintf (htitle5,"Localx in BPIX %s", Region[i]);
    sprintf (htitle6,"Localx in FPIX %s", Region[i]);
    
    sprintf (hname1,"Localx_TIB_%i",i+1);
    sprintf (hname2,"Localx_TOB_%i",i+1);
    sprintf (hname3,"Localx_TID_%i",i+1);
    sprintf (hname4,"Localx_TEC_%i",i+1);
    sprintf (hname5,"Localx_BPIX_%i",i+1);
    sprintf (hname6,"Localx_FPIX_%i",i+1);
   
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
    h1lx[i]  = fDBE->book1D (hname1, htitle1, nbin , low2[0] , high2[0]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
    h2lx[i]  = fDBE->book1D (hname2, htitle2, nbin , low2[1] , high2[1]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
    h3lx[i]  = fDBE->book1D (hname3, htitle3, nbin , low2[2] , high2[2]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
    h4lx[i]  = fDBE->book1D (hname4, htitle4, nbin , low2[3] , high2[3]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
    h5lx[i]  = fDBE->book1D (hname5, htitle5, nbin , low2[4] , high2[4]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
    h6lx[i]  = fDBE->book1D (hname6, htitle6, nbin , low2[5] , high2[5]);
   
   }


const float high3[] = {6.0, 10., 5.6, 10.5, 3.4, 0.52};
const float low3[] = {-6.0, -10., -5.6, -10.5, -3.4, -0.52};

   for(int i=0; i<12; i++) {
  
    sprintf (htitle1,"Localy in TIB %s", Region[i]);
    sprintf (htitle2,"Localy in TOB %s", Region[i]);
    sprintf (htitle3,"Localy in TID %s", Region[i]);
    sprintf (htitle4,"Localy in TEC %s", Region[i]);
    sprintf (htitle5,"Localy in BPIX %s", Region[i]);
    sprintf (htitle6,"Localy in FPIX %s", Region[i]);
    
    sprintf (hname1,"Localy_TIB_%i",i+1);
    sprintf (hname2,"Localy_TOB_%i",i+1);
    sprintf (hname3,"Localy_TID_%i",i+1);
    sprintf (hname4,"Localy_TEC_%i",i+1);
    sprintf (hname5,"Localy_BPIX_%i",i+1);
    sprintf (hname6,"Localy_FPIX_%i",i+1);
   
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
    h1ly[i]  = fDBE->book1D (hname1, htitle1, nbin , low3[0] , high3[0]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
    h2ly[i]  = fDBE->book1D (hname2, htitle2, nbin , low3[1] , high3[1]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
    h3ly[i]  = fDBE->book1D (hname3, htitle3, nbin , low3[2] , high3[2]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
    h4ly[i]  = fDBE->book1D (hname4, htitle4, nbin , low3[3] , high3[3]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
    h5ly[i]  = fDBE->book1D (hname5, htitle5, nbin , low3[4] , high3[4]);
    fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
    h6ly[i]  = fDBE->book1D (hname6, htitle6, nbin , low3[5] , high3[5]);
   
   }
   
   
  }
}

void TrackerHitAnalyzer::endJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 321 of file TrackerHitAnalyzer.cc.

References DQMStore::bookProfile(), fDBE, fOutputFile, DQMStore::get(), MonitorElement::getTH2F(), htofeta, htofeta_profile, htofphi, htofphi_profile, htofr, htofr_profile, htofz, htofz_profile, reco::return(), DQMStore::save(), and DQMStore::setCurrentFolder().

{
  //before check that histos are there....

  // check if ME still there (and not killed by MEtoEDM for memory saving)
  if( fDBE )
    {
      // check existence of first histo in the list
      if (! fDBE->get("TrackerHitsV/TrackerHit/tof_eta")) return;
    }
  else
    return;

  fDBE->setCurrentFolder("TrackerHitsV/TrackerHit/");
  htofeta_profile  = fDBE->bookProfile ("tof_eta_profile",htofeta->getTH2F()->ProfileX());
  htofphi_profile  = fDBE->bookProfile("tof_phi_profile", htofphi->getTH2F()->ProfileX());
  htofr_profile  = fDBE->bookProfile("tof_r_profile",htofr->getTH2F()->ProfileX());
  htofz_profile  = fDBE->bookProfile("tof_z_profile", htofz->getTH2F()->ProfileX());


  if ( fOutputFile.size() != 0 && fDBE ) fDBE->save(fOutputFile);

  return ;

}

Member Data Documentation

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlBrlHigh_Token_

private

Definition at line 51 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlBrlLow_Token_

private

Definition at line 51 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlFwdHigh_Token_

private

Definition at line 52 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlFwdLow_Token_

private

Definition at line 52 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTECHigh_Token_

private

Definition at line 56 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTECLow_Token_

private

Definition at line 56 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIBHigh_Token_

private

Definition at line 53 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIBLow_Token_

private

Definition at line 53 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIDHigh_Token_

private

Definition at line 55 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIDLow_Token_

private

Definition at line 55 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTOBHigh_Token_

private

Definition at line 54 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTOBLow_Token_

private

Definition at line 54 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edm::EDGetTokenT<edm::SimTrackContainer> TrackerHitAnalyzer::edmSimTrackContainerToken_

private

Definition at line 57 of file TrackerHitAnalyzer.h.

Referenced by analyze().

DQMStore* TrackerHitAnalyzer::fDBE

private

Definition at line 59 of file TrackerHitAnalyzer.h.

Referenced by beginRun(), and endJob().

std::string TrackerHitAnalyzer::fOutputFile

private

Definition at line 111 of file TrackerHitAnalyzer.h.

Referenced by endJob().

MonitorElement* TrackerHitAnalyzer::h1e[12]

private

Definition at line 69 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h1ex[12]

private

Definition at line 76 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h1ey[12]

private

Definition at line 83 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h1ez[12]

private

Definition at line 90 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h1lx[12]

private

Definition at line 97 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h1ly[12]

private

Definition at line 104 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h2e[12]

private

Definition at line 70 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h2ex[12]

private

Definition at line 77 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h2ey[12]

private

Definition at line 84 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h2ez[12]

private

Definition at line 91 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h2lx[12]

private

Definition at line 98 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h2ly[12]

private

Definition at line 105 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h3e[12]

private

Definition at line 71 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h3ex[12]

private

Definition at line 78 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h3ey[12]

private

Definition at line 85 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h3ez[12]

private

Definition at line 92 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h3lx[12]

private

Definition at line 99 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h3ly[12]

private

Definition at line 106 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h4e[12]

private

Definition at line 72 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h4ex[12]

private

Definition at line 79 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h4ey[12]

private

Definition at line 86 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h4ez[12]

private

Definition at line 93 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h4lx[12]

private

Definition at line 100 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h4ly[12]

private

Definition at line 107 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h5e[12]

private

Definition at line 73 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h5ex[12]

private

Definition at line 80 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h5ey[12]

private

Definition at line 87 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h5ez[12]

private

Definition at line 94 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h5lx[12]

private

Definition at line 101 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h5ly[12]

private

Definition at line 108 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h6e[12]

private

Definition at line 74 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h6ex[12]

private

Definition at line 81 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h6ey[12]

private

Definition at line 88 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h6ez[12]

private

Definition at line 95 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h6lx[12]

private

Definition at line 102 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::h6ly[12]

private

Definition at line 109 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackerHitAnalyzer::htofeta

private

Definition at line 61 of file TrackerHitAnalyzer.h.

Referenced by analyze(), beginRun(), and endJob().

MonitorElement* TrackerHitAnalyzer::htofeta_profile

private

Definition at line 65 of file TrackerHitAnalyzer.h.

Referenced by endJob().

MonitorElement* TrackerHitAnalyzer::htofphi

private

Definition at line 62 of file TrackerHitAnalyzer.h.

Referenced by analyze(), beginRun(), and endJob().

MonitorElement* TrackerHitAnalyzer::htofphi_profile

private

Definition at line 66 of file TrackerHitAnalyzer.h.

Referenced by endJob().

MonitorElement* TrackerHitAnalyzer::htofr

private

Definition at line 63 of file TrackerHitAnalyzer.h.

Referenced by analyze(), beginRun(), and endJob().

MonitorElement* TrackerHitAnalyzer::htofr_profile

private

Definition at line 67 of file TrackerHitAnalyzer.h.

Referenced by endJob().

MonitorElement* TrackerHitAnalyzer::htofz

private

Definition at line 64 of file TrackerHitAnalyzer.h.

Referenced by analyze(), beginRun(), and endJob().

MonitorElement* TrackerHitAnalyzer::htofz_profile

private

Definition at line 68 of file TrackerHitAnalyzer.h.

Referenced by endJob().

bool TrackerHitAnalyzer::verbose_

private

Definition at line 49 of file TrackerHitAnalyzer.h.

Referenced by beginRun().