CTPPSPixelDQMSource Class Reference

Inheritance diagram for CTPPSPixelDQMSource:

Public Member Functions
	CTPPSPixelDQMSource (const edm::ParameterSet &ps)
	~CTPPSPixelDQMSource () override
Public Member Functions inherited from DQMEDAnalyzer
void	beginRun (edm::Run const &, edm::EventSetup const &) final
void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer ()
void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
void	endRun (edm::Run const &run, edm::EventSetup const &eSetup) override
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 Member Functions
int	getDet (int id)
int	getPixPlane (int id)
int	getPlaneIndex (int arm, int station, int rp, int plane)
int	getRPindex (int arm, int station, int rp)
int	prIndex (int rp, int plane)

Private Attributes
int	cluSizeMax
const int	ClusMultMAX = 10
int	ClusMultPlane [RPotsTotalNumber][NplaneMAX]
MonitorElement *	h2CluSize [NArms][NStationMAX]
MonitorElement *	h2ClusMultipl [NArms][NStationMAX]
MonitorElement *	h2HitsMultipl [NArms][NStationMAX]
MonitorElement *	h2HitsMultROC [RPotsTotalNumber]
MonitorElement *	h2PlaneActive [NArms][NStationMAX]
MonitorElement *	h2xyHits [RPotsTotalNumber][NplaneMAX]
MonitorElement *	h2xyROCadc [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	h2xyROCHits [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hBX
MonitorElement *	hBXshort
MonitorElement *	hHitsMult [RPotsTotalNumber][NplaneMAX]
const int	hitMultMAX = 300
int	HitsMultPlane [RPotsTotalNumber][NplaneMAX]
int	HitsMultROC [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hp2HitsMultROC_LS [RPotsTotalNumber]
MonitorElement *	hp2HitsOcc [NArms][NStationMAX]
MonitorElement *	hp2xyADC [RPotsTotalNumber][NplaneMAX]
MonitorElement *	hRPotActivBX [RPotsTotalNumber]
MonitorElement *	hRPotActivBXall [RPotsTotalNumber]
MonitorElement *	hRPotActivBXroc [RPotsTotalNumber]
MonitorElement *	hRPotActivPlanes [RPotsTotalNumber]
MonitorElement *	hRPotActivROCs [RPotsTotalNumber]
const int	IndexNotValid = 0
int	multHitsMax
long int	nEvents = 0
int	pixColMAX = 156
int	pixRowMAX = 160
int	ROCSizeInX = pixRowMAX/2
int	ROCSizeInY = pixColMAX/3
int	RPindexValid [RPotsTotalNumber]
int	RPstatus [NStationMAX][NRPotsMAX]
unsigned int	rpStatusWord = 0x8000
int	StationStatus [NStationMAX]
CTPPSPixelIndices	thePixIndices
edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelCluster > >	tokenCluster
edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelDigi > >	tokenDigi
unsigned int	verbosity

Static Private Attributes
static int	NArms =2
static int	NPlaneBins = NplaneMAX*NRPotBinsInStation
static int	NplaneMAX =6
static int	NROCsMAX = 6
static int	NRPotBinsInStation = RPn_last-RPn_first
static int	NRPotsMAX =6
static int	NStationMAX =3
static int	RPn_first = 3
static int	RPn_last = 4
static int	RPotsTotalNumber =NArms*NStationMAX*NRPotsMAX

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 32 of file CTPPSPixelDQMSource.cc.

Constructor & Destructor Documentation

CTPPSPixelDQMSource::CTPPSPixelDQMSource

(

const edm::ParameterSet &

ps

)

Definition at line 134 of file CTPPSPixelDQMSource.cc.

References edm::ParameterSet::getParameter(), tokenCluster, and tokenDigi.

                                                                  :
  verbosity(ps.getUntrackedParameter<unsigned int>("verbosity", 0)),
  rpStatusWord(ps.getUntrackedParameter<unsigned int>("RPStatusWord",0x8000))
{
 tokenDigi = consumes<DetSetVector<CTPPSPixelDigi> >(ps.getParameter<edm::InputTag>("tagRPixDigi"));
 tokenCluster=consumes<DetSetVector<CTPPSPixelCluster>>(ps.getParameter<edm::InputTag>("tagRPixCluster"));

}

CTPPSPixelDQMSource::~CTPPSPixelDQMSource ( )

override

Definition at line 145 of file CTPPSPixelDQMSource.cc.

{
}

Member Function Documentation

void CTPPSPixelDQMSource::analyze ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

overrideprotectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 356 of file CTPPSPixelDQMSource.cc.

References ecalMGPA::adc(), CTPPSDetId::arm(), edm::DetSet< T >::begin(), edm::EventBase::bunchCrossing(), cluSizeMax, ClusMultPlane, cuy::col, MonitorElement::Fill(), HcalObjRepresent::Fill(), objects.autophobj::float, getDet(), getPixPlane(), getPlaneIndex(), getRPindex(), h2CluSize, h2ClusMultipl, h2HitsMultipl, h2HitsMultROC, h2PlaneActive, h2xyHits, h2xyROCadc, h2xyROCHits, hBX, hBXshort, hHitsMult, HitsMultPlane, HitsMultROC, hp2HitsMultROC_LS, hp2HitsOcc, hp2xyADC, hRPotActivBX, hRPotActivBXall, hRPotActivBXroc, hRPotActivPlanes, hRPotActivROCs, diffTreeTool::index, edm::HandleBase::isValid(), SiStripPI::max, multHitsMax, NArms, nEvents, np, NplaneMAX, NROCsMAX, NRPotsMAX, NStationMAX, AlCaHLTBitMon_ParallelJobs::p, pixColMAX, pixRowMAX, prIndex(), alignCSCRings::r, CTPPSDetId::rp(), RPindexValid, RPstatus, mathSSE::sqrt(), CTPPSDetId::station(), relativeConstraints::station, StationStatus, thePixIndices, tokenCluster, tokenDigi, CTPPSPixelIndices::transformToROC(), verbosity, and DetId::VeryForward.

{
  ++nEvents;
  int lumiId = event.getLuminosityBlock().id().luminosityBlock();
  if(lumiId<0) lumiId=0;

  int RPactivity[NArms][NRPotsMAX], digiSize[NArms][NRPotsMAX];
  for(int arm = 0; arm <2; arm++) { 
    for(int rp=0; rp<NRPotsMAX; rp++) {
      RPactivity[arm][rp] = digiSize[arm][rp] = 0;
    }
  }
  for(int ind=0; ind<2*3*NRPotsMAX; ind++) { 
    for(int p=0; p<NplaneMAX; p++) {
      HitsMultPlane[ind][p] = 0;
      ClusMultPlane[ind][p] = 0;
    }
  }
  for(int ind=0; ind<2*3*NRPotsMAX*NplaneMAX; ind++)  {
    for(int rp=0; rp<NROCsMAX; rp++) {
      HitsMultROC[ind][rp] = 0;
    }
  }
  Handle< DetSetVector<CTPPSPixelDigi> > pixDigi;
  event.getByToken(tokenDigi, pixDigi);

  Handle< DetSetVector<CTPPSPixelCluster> > pixClus;
  event.getByToken(tokenCluster, pixClus);

  hBX->Fill(event.bunchCrossing());
  hBXshort->Fill(event.bunchCrossing());

  bool valid = false;
  valid |= pixDigi.isValid();
//  valid |= Clus.isValid();

  if(!valid && verbosity) LogPrint("CTPPSPixelDQMSource") <<"No valid data in Event "<<nEvents;

  if(pixDigi.isValid()) {
    for(const auto &ds_digi : *pixDigi)
    {
      int idet = getDet(ds_digi.id);
      if(idet != DetId::VeryForward) {
        if(verbosity>1) LogPrint("CTPPSPixelDQMSource") <<"not CTPPS: ds_digi.id"<<ds_digi.id;
        continue;
      }
      //   int subdet = getSubdet(ds_digi.id);

      int plane = getPixPlane(ds_digi.id);
   
      CTPPSDetId theId(ds_digi.id);
      int arm = theId.arm()&0x1;
      int station = theId.station()&0x3;
      int rpot = theId.rp()&0x7;
      RPactivity[arm][rpot] = 1;
      ++digiSize[arm][rpot];

      if(StationStatus[station] && RPstatus[station][rpot]) {

        hp2HitsOcc[arm][station]->Fill(plane,rpot,(int)ds_digi.data.size());
        h2HitsMultipl[arm][station]->Fill(prIndex(rpot,plane),ds_digi.data.size());
        h2PlaneActive[arm][station]->Fill(plane,rpot);

        int index = getRPindex(arm,station,rpot);
        HitsMultPlane[index][plane] += ds_digi.data.size();
        if(RPindexValid[index]) {
          hHitsMult[index][plane]->Fill(ds_digi.data.size());
        }
        int rocHistIndex = getPlaneIndex(arm,station,rpot,plane);

        for(DetSet<CTPPSPixelDigi>::const_iterator dit = ds_digi.begin();
            dit != ds_digi.end(); ++dit) {
          int row = dit->row();
          int col = dit->column();
          int adc = dit->adc();

          if(RPindexValid[index]) {
            h2xyHits[index][plane]->Fill(col,row);
            hp2xyADC[index][plane]->Fill(col,row,adc);
            int colROC, rowROC;
            int trocId;
            if(!thePixIndices.transformToROC(col,row, trocId, colROC, rowROC)) {
              if(trocId>=0 && trocId<NROCsMAX) {
                h2xyROCHits[rocHistIndex][trocId]->Fill(colROC,rowROC);
                h2xyROCadc[rocHistIndex][trocId]->Fill(colROC,rowROC,adc);
                ++HitsMultROC[rocHistIndex][trocId];
              }
            }
          } //end if(RPindexValid[index]) {
        }
        if(int(ds_digi.data.size()) > multHitsMax) multHitsMax = ds_digi.data.size();

      } // end  if(StationStatus[station]) {
    } // end for(const auto &ds_digi : *pixDigi)
  } //if(pixDigi.isValid()) {

  if(pixClus.isValid()) 
    for(const auto &ds : *pixClus)
    {
      int idet = getDet(ds.id);
      if(idet != DetId::VeryForward) {
        if(verbosity>1) LogPrint("CTPPSPixelDQMSource") <<"not CTPPS: cluster.id"<<ds.id;
        continue;
      }
      //   int subdet = getSubdet(ds.id);

      int plane = getPixPlane(ds.id);

      CTPPSDetId theId(ds.id);
      int arm = theId.arm()&0x1;
      int station = theId.station()&0x3;
      int rpot = theId.rp()&0x7;

        int index = getRPindex(arm,station,rpot);
        ++ClusMultPlane[index][plane];

      unsigned int minRow=pixRowMAX, maxRow=0;
      unsigned int minCol=pixColMAX, maxCol=0;
      for (const auto &p : ds) {
        unsigned int max= p.minPixelRow() + p.rowSpan()+1;
        if(minRow > p.minPixelRow()) minRow = p.minPixelRow();
        if(maxRow < max) maxRow = max;
        max= p.minPixelCol() + p.colSpan()+1;
        if(minCol > p.minPixelCol()) minCol = p.minPixelCol();
        if(maxCol < max) maxCol = max;

      }
      int drow = maxRow - minRow;
      int dcol = maxCol - minCol;
      float clusize= sqrt(drow*drow + dcol*dcol);
      if(cluSizeMax < int(clusize)) cluSizeMax = clusize;
      h2CluSize[arm][station]->Fill(prIndex(rpot,plane),clusize);

   } // end if(pixClus.isValid()) for(const auto &ds : *pixClus)

  for(int arm=0; arm<2; arm++) {
    for(int stn=0; stn<NStationMAX; stn++) {
      for(int rp=0; rp<NRPotsMAX; rp++) {
        int index = getRPindex(arm,stn,rp);
        if(RPindexValid[index]==0) continue;
        if(RPactivity[arm][rp]==0) continue;

        for(int p=0; p<NplaneMAX; p++)
          h2ClusMultipl[arm][stn]->Fill(prIndex(rp,p),ClusMultPlane[index][p]);

        int np = 0; 
        for(int p=0; p<NplaneMAX; p++) if(HitsMultPlane[index][p]>0) np++;
    for(int p=0; p<=NplaneMAX; p++) {
      if(p == np) hRPotActivPlanes[index]->Fill(p,1.);
      else hRPotActivPlanes[index]->Fill(p,0.);
    }
        if(np>5) hRPotActivBX[index]->Fill(event.bunchCrossing());
        hRPotActivBXall[index]->Fill(event.bunchCrossing(),float(digiSize[arm][rp]));

        int rocf[NplaneMAX];
        for(int r=0; r<NROCsMAX; r++) rocf[r]=0;
        for(int p=0; p<NplaneMAX; p++) {
          int indp = getPlaneIndex(arm,stn,rp,p);
          for(int r=0; r<NROCsMAX; r++) if(HitsMultROC[indp][r] > 0) ++rocf[r];
          for(int r=0; r<NROCsMAX; r++) { 
            h2HitsMultROC[index]->Fill(p,r,HitsMultROC[indp][r]);
            hp2HitsMultROC_LS[index]->Fill(lumiId,p*NROCsMAX+r,HitsMultROC[indp][r]);
          }
        }
        int max = 0;
        for(int r=0; r<NROCsMAX; r++) 
          if(max < rocf[r]) { max = rocf[r]; }

        for(int r=0; r<NROCsMAX; r++) {
      for(int p=0; p<=NplaneMAX; p++)
        if(p==rocf[r]) hRPotActivROCs[index]->Fill(r,rocf[r],1.);
        else hRPotActivROCs[index]->Fill(r,p,0);
    }

        if(max > 4) hRPotActivBXroc[index]->Fill(event.bunchCrossing());
      } //end for(int rp=0; rp<NRPotsMAX; rp++) {
    }
  } //end for(int arm=0; arm<2; arm++) {

  if((nEvents % 100)) return;
  if(verbosity) LogPrint("CTPPSPixelDQMSource")<<"analyze event "<<nEvents;
}

void CTPPSPixelDQMSource::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 178 of file CTPPSPixelDQMSource.cc.

References CTPPSDetId::armName(), DQMStore::IBooker::book1D(), DQMStore::IBooker::book1DD(), DQMStore::IBooker::book2D(), DQMStore::IBooker::book2DD(), DQMStore::IBooker::bookProfile(), DQMStore::IBooker::bookProfile2D(), DQMStore::IBooker::cd(), ClusMultMAX, MillePedeFileConverter_cfg::e, JetChargeProducer_cfi::exp, getPlaneIndex(), CTPPSDetId::getRPId(), getRPindex(), CTPPSDetId::getStationId(), MonitorElement::getTH2D(), MonitorElement::getTH2F(), MonitorElement::getTProfile2D(), h2CluSize, h2ClusMultipl, h2HitsMultipl, h2HitsMultROC, h2PlaneActive, h2xyHits, h2xyROCadc, h2xyROCHits, hBX, hBXshort, hHitsMult, hitMultMAX, hp2HitsMultROC_LS, hp2HitsOcc, hp2xyADC, hRPotActivBX, hRPotActivBXall, hRPotActivBXroc, hRPotActivPlanes, hRPotActivROCs, mps_fire::i, diffTreeTool::index, cmsBatch::log, n0, pileupCalc::nbins, CTPPSDetId::nFull, CTPPSDetId::nPath, NPlaneBins, NplaneMAX, NROCsMAX, NRPotBinsInStation, CTPPSDetId::nShort, NStationMAX, AlCaHLTBitMon_ParallelJobs::p, pixRowMAX, alignCSCRings::r, ROCSizeInX, ROCSizeInY, RPindexValid, RPn_first, RPn_last, RPstatus, alignCSCRings::s, sd, CTPPSDetId::sdTrackingPixel, DQMStore::IBooker::setCurrentFolder(), CTPPSDetId::setRP(), CTPPSDetId::setStation(), StationStatus, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  ibooker.cd();
  ibooker.setCurrentFolder("CTPPS/TrackingPixel");
  char s[50];
  hBX = ibooker.book1D("events per BX", "ctpps_pixel;Event.BX", 4002, -1.5, 4000. +0.5);
  hBXshort = ibooker.book1D("events per BX(short)", "ctpps_pixel;Event.BX", 102, -1.5, 100. + 0.5);

 for(int arm=0; arm<2; arm++) {
   CTPPSDetId ID(CTPPSDetId::sdTrackingPixel,arm,0);
   string sd, armTitle;
   ID.armName(sd, CTPPSDetId::nPath);
   ID.armName(armTitle, CTPPSDetId::nFull);

   ibooker.setCurrentFolder(sd);

 for(int stn=2; stn<NStationMAX; stn++) {
   if(StationStatus[stn]==0) continue;
   ID.setStation(stn);
   string stnd, stnTitle;

   CTPPSDetId(ID.getStationId()).stationName(stnd, CTPPSDetId::nPath);
   CTPPSDetId(ID.getStationId()).stationName(stnTitle, CTPPSDetId::nFull);

   ibooker.setCurrentFolder(stnd);

   string st = "planes activity";
   string st2 = ": " + stnTitle;

   h2PlaneActive[arm][stn] = ibooker.book2DD(st,st+st2+";Plane #",
        NplaneMAX,0,NplaneMAX, NRPotBinsInStation, RPn_first,RPn_last);
   TH2D *h = h2PlaneActive[arm][stn]->getTH2D();
   h->SetOption("colz");
   TAxis *yah = h->GetYaxis();
   
   st = "hit average multiplicity in planes";
   int PlaneMultCut = 20;
   hp2HitsOcc[arm][stn]= ibooker.bookProfile2D(st,st+st2+";Plane #",
     NplaneMAX, 0, NplaneMAX, NRPotBinsInStation, RPn_first,RPn_last,-1,PlaneMultCut,"");
   TProfile2D *h2 = hp2HitsOcc[arm][stn]->getTProfile2D();
   h2->SetOption("textcolz");
   TAxis *yah2 = h2->GetYaxis();

   st = "hit multiplicity in planes";
   string st3 = ";PlaneIndex(=pixelPot*PlaneMAX + plane)";
   h2HitsMultipl[arm][stn]= ibooker.book2DD(st,st+st2+st3+";multiplicity",
    NPlaneBins,0,NPlaneBins,hitMultMAX,0,hitMultMAX);
   h2HitsMultipl[arm][stn]->getTH2D()->SetOption("colz");

   st = "cluster multiplicity in planes";
   h2ClusMultipl[arm][stn] = ibooker.book2DD(st,st+st2+st3+";multiplicity",
    NPlaneBins,0,NPlaneBins, ClusMultMAX,0,ClusMultMAX);
   h2ClusMultipl[arm][stn]->getTH2D()->SetOption("colz");

   st = "cluster span in planes";
  const int nyClus = 9; //18;
  float xCh[NPlaneBins+1];
  float yClus[nyClus+1];
  for(int i=0; i<=NPlaneBins; i++) xCh[i]=i;
  double n0 = 1; //1./CluSizeMAX;
  double lnA = log(2.);
  yClus[0] = n0; //yClus[1] = n0;
  for(int j=0; j<nyClus; j++) yClus[j+1] = n0*exp(j*lnA);

  h2CluSize[arm][stn] = ibooker.book2D(st,st+st2+st3+";Cluster size",
    NPlaneBins,xCh,nyClus,yClus);
   h2CluSize[arm][stn]->getTH2F()->SetOption("colz");

//--------- RPots ---
   int pixBinW = 4;
     for(int rp=RPn_first; rp<RPn_last; rp++) { // only installed pixel pots
       ID.setRP(rp);
       string rpd, rpTitle;
       CTPPSDetId(ID.getRPId()).rpName(rpTitle, CTPPSDetId::nShort);
    yah->SetBinLabel(rp - RPn_first +1, rpTitle.c_str()); // h
       yah2->SetBinLabel(rp - RPn_first +1, rpTitle.c_str()); //h2

       if(RPstatus[stn][rp]==0) continue;
       int indexP = getRPindex(arm,stn,rp);
       RPindexValid[indexP] = 1;

       CTPPSDetId(ID.getRPId()).rpName(rpTitle, CTPPSDetId::nFull);
       CTPPSDetId(ID.getRPId()).rpName(rpd, CTPPSDetId::nPath);

       ibooker.setCurrentFolder(rpd);

       hRPotActivPlanes[indexP] = 
        ibooker.bookProfile("number of fired planes per event", rpTitle+";nPlanes;Probability",
     NplaneMAX+1, -0.5, NplaneMAX+0.5, -0.5,NplaneMAX+0.5,"");
       hRPotActivBX[indexP] = 
        ibooker.book1D("5 fired planes per BX", rpTitle+";Event.BX", 4002, -1.5, 4000.+0.5);

       h2HitsMultROC[indexP] = ibooker.bookProfile2D("ROCs hits multiplicity per event",
       rpTitle+";plane # ;ROC #", NplaneMAX,-0.5,NplaneMAX-0.5, NROCsMAX,-0.5,NROCsMAX-0.5,
    0.,ROCSizeInX*ROCSizeInY,"");
       h2HitsMultROC[indexP]->getTProfile2D()->SetOption("colztext");
       h2HitsMultROC[indexP]->getTProfile2D()->SetMinimum(1.e-10);

    
       hp2HitsMultROC_LS[indexP]=ibooker.bookProfile2D("ROCs_hits_multiplicity_per_event vs LS",
     rpTitle+";LumiSection;Plane#___ROC#", 1000,0.,1000.,
     NplaneMAX*NROCsMAX,0.,double(NplaneMAX*NROCsMAX),0.,ROCSizeInX*ROCSizeInY,"");
       hp2HitsMultROC_LS[indexP]->getTProfile2D()->SetOption("colz");
       hp2HitsMultROC_LS[indexP]->getTProfile2D()->SetMinimum(1.0e-10);
       hp2HitsMultROC_LS[indexP]->getTProfile2D()->SetCanExtend(TProfile2D::kXaxis);
       TAxis *yahp2 = hp2HitsMultROC_LS[indexP]->getTProfile2D()->GetYaxis();
       for(int p=0; p<NplaneMAX; p++) {
     sprintf(s,"plane%d_0",p);
     yahp2->SetBinLabel(p*NplaneMAX+1,s); 
     for(int r=1; r<NROCsMAX; r++) {
       sprintf(s,"   %d_%d",p,r);
       yahp2->SetBinLabel(p*NplaneMAX+r+1,s);
     }
       }

       hRPotActivROCs[indexP] = ibooker.bookProfile2D("number of fired aligned_ROCs per event",
     rpTitle+";ROC ID;number of fired aligned ROCs", NROCsMAX,-0.5,NROCsMAX-0.5, 
     7,-0.5,6.5, 0,NROCsMAX,"");
       hRPotActivROCs[indexP]->getTProfile2D()->SetOption("colz");

       ibooker.setCurrentFolder(rpd+"/latency");

       hRPotActivBXroc[indexP] = 
        ibooker.book1D("4 fired ROCs per BX", rpTitle+";Event.BX", 4002, -1.5, 4000.+0.5);

       hRPotActivBXall[indexP] =
        ibooker.book1D("hits per BX", rpTitle+";Event.BX", 4002, -1.5, 4000.+0.5);

       int nbins = pixRowMAX/pixBinW;

       for(int p=0; p<NplaneMAX; p++) {
         sprintf(s,"plane_%d",p);
         string pd = rpd+"/"+string(s);
         ibooker.setCurrentFolder(pd);
         string st1 = ": "+rpTitle+"_"+string(s);
         st = "hits position";
         h2xyHits[indexP][p]=ibooker.book2DD(st,st1+";pix col;pix row", 
       nbins,0,pixRowMAX,nbins,0,pixRowMAX);
     h2xyHits[indexP][p]->getTH2D()->SetOption("colz");

     st = "adc average value";
     hp2xyADC[indexP][p]=ibooker.bookProfile2D(st,st1+";pix col;pix row",
       nbins,0,pixRowMAX,nbins,0,pixRowMAX, 0.,512.,"");
     hp2xyADC[indexP][p]->getTProfile2D()->SetOption("colz");

         st = "hits multiplicity";
         hHitsMult[indexP][p]=ibooker.book1DD(st,st1+";number of hits;N / 1 hit",
       hitMultMAX,-0.5,hitMultMAX-0.5);

         ibooker.setCurrentFolder(pd + "/ROCs");
         int index = getPlaneIndex(arm,stn,rp,p);

         for(int roc=0; roc<NROCsMAX; roc++) {
       sprintf(s,"ROC_%d",roc);
       string st2 = st1 + "_" + string(s);
           ibooker.setCurrentFolder(pd + "/ROCs/" + string(s));

           h2xyROCHits[index][roc]=ibooker.book2DD("hits",st2+";pix col;pix row",
        ROCSizeInY,0,ROCSizeInY, ROCSizeInX,0,ROCSizeInX);
           h2xyROCHits[index][roc]->getTH2D()->SetOption("colz");

           string st = "adc average value";
           h2xyROCadc[index][roc]=ibooker.bookProfile2D(st,st2+";pix col;pix row",
        ROCSizeInY,0,ROCSizeInY,ROCSizeInX,0,ROCSizeInX, 0.,512.,"");
           h2xyROCadc[index][roc]->getTProfile2D()->SetOption("colz");
         }
       } // end of for(int p=0; p<NplaneMAX;..

     } // end for(int rp=0; rp<NRPotsMAX;...
   } // end of for(int stn=0; stn<
  } //end of for(int arm=0; arm<2;...

 return;
}

void CTPPSPixelDQMSource::dqmBeginRun ( edm::Run const &  run, edm::EventSetup const &    )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 151 of file CTPPSPixelDQMSource.cc.

References cluSizeMax, CTPPSPixelIndices::getDefaultColDetSize(), CTPPSPixelIndices::getDefaultRowDetSize(), multHitsMax, nEvents, NRPotsMAX, pixColMAX, pixRowMAX, ROCSizeInX, ROCSizeInY, RPindexValid, RPstatus, rpStatusWord, StationStatus, thePixIndices, and verbosity.

{
  if(verbosity) LogPrint("CTPPSPixelDQMSource") <<"RPstatusWord= "<<rpStatusWord;
  nEvents = 0;

  pixRowMAX = thePixIndices.getDefaultRowDetSize();
  pixColMAX = thePixIndices.getDefaultColDetSize();
  ROCSizeInX = pixRowMAX/2;  // ROC row size in pixels = 80
  ROCSizeInY = pixColMAX/3;

 unsigned int rpSts = rpStatusWord<<1;
 for(int stn=0; stn<3; stn++) {
   int stns = 0;
   for(int rp=0; rp<NRPotsMAX; rp++) {
     rpSts = (rpSts >> 1); RPstatus[stn][rp] = rpSts&1;
     if(RPstatus[stn][rp] > 0) stns = 1;
   }
   StationStatus[stn]=stns;
 }
 
 for(int ind=0; ind<2*3*NRPotsMAX; ind++) RPindexValid[ind] = 0;

 multHitsMax = cluSizeMax = -1;
}

void CTPPSPixelDQMSource::endRun ( edm::Run const &  run, edm::EventSetup const &  eSetup  )

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 540 of file CTPPSPixelDQMSource.cc.

References cluSizeMax, DEFINE_FWK_MODULE, multHitsMax, nEvents, edm::RunBase::run(), and verbosity.

{
 if(!verbosity) return; 
 LogPrint("CTPPSPixelDQMSource") 
  <<"end of Run "<<run.run()<<": "<<nEvents<<" events\n"
  <<"multHitsMax= "<<multHitsMax <<"   cluSizeMax= "<<cluSizeMax;
}

int CTPPSPixelDQMSource::getDet ( int  id )

inlineprivate

Definition at line 117 of file CTPPSPixelDQMSource.cc.

References DetId::kDetOffset.

Referenced by analyze().

    { return (id>>DetId::kDetOffset)&0xF; }

int CTPPSPixelDQMSource::getPixPlane ( int  id )

inlineprivate

Definition at line 119 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

    { return ((id>>16)&0x7); }

int CTPPSPixelDQMSource::getPlaneIndex ( int  arm, int  station, int  rp, int  plane  )

inlineprivate

Definition at line 108 of file CTPPSPixelDQMSource.cc.

References getRPindex().

Referenced by analyze(), and bookHistograms().

                                                             {
    if(plane<0 || plane>=NplaneMAX) return(IndexNotValid);
    int rc = getRPindex(arm, station, rp);
    if(rc == IndexNotValid) return(IndexNotValid);
    return(rc*NplaneMAX + plane);
  }

int CTPPSPixelDQMSource::getRPindex ( int  arm, int  station, int  rp  )

inlineprivate

Definition at line 101 of file CTPPSPixelDQMSource.cc.

References relativeConstraints::station.

Referenced by analyze(), bookHistograms(), and getPlaneIndex().

                                               {
    if(arm<0 || station<0 || rp<0) return(IndexNotValid);
    if(arm>1 || station>=NStationMAX || rp>=NRPotsMAX) return(IndexNotValid);
    int rc = (arm*NStationMAX+station)*NRPotsMAX + rp;
    return(rc);
  }

int CTPPSPixelDQMSource::prIndex ( int  rp, int  plane  )

inlineprivate

Definition at line 115 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

    {return((rp - RPn_first)*NplaneMAX + plane);}

Member Data Documentation

int CTPPSPixelDQMSource::cluSizeMax

private

Definition at line 123 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), dqmBeginRun(), and endRun().

const int CTPPSPixelDQMSource::ClusMultMAX = 10

private

Definition at line 57 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::ClusMultPlane[RPotsTotalNumber][NplaneMAX]

private

Definition at line 93 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

MonitorElement* CTPPSPixelDQMSource::h2CluSize[NArms][NStationMAX]

private

Definition at line 74 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2ClusMultipl[NArms][NStationMAX]

private

Definition at line 73 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2HitsMultipl[NArms][NStationMAX]

private

Definition at line 71 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2HitsMultROC[RPotsTotalNumber]

private

Definition at line 82 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2PlaneActive[NArms][NStationMAX]

private

Definition at line 72 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2xyHits[RPotsTotalNumber][NplaneMAX]

private

Definition at line 86 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2xyROCadc[RPotsTotalNumber *NplaneMAX][NROCsMAX]

private

Definition at line 89 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2xyROCHits[RPotsTotalNumber *NplaneMAX][NROCsMAX]

private

Definition at line 88 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hBX

private

Definition at line 68 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement * CTPPSPixelDQMSource::hBXshort

private

Definition at line 68 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hHitsMult[RPotsTotalNumber][NplaneMAX]

private

Definition at line 85 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const int CTPPSPixelDQMSource::hitMultMAX = 300

private

Definition at line 56 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::HitsMultPlane[RPotsTotalNumber][NplaneMAX]

private

Definition at line 92 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

int CTPPSPixelDQMSource::HitsMultROC[RPotsTotalNumber *NplaneMAX][NROCsMAX]

private

Definition at line 91 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

MonitorElement* CTPPSPixelDQMSource::hp2HitsMultROC_LS[RPotsTotalNumber]

private

Definition at line 83 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hp2HitsOcc[NArms][NStationMAX]

private

Definition at line 70 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hp2xyADC[RPotsTotalNumber][NplaneMAX]

private

Definition at line 87 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBX[RPotsTotalNumber]

private

Definition at line 80 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBXall[RPotsTotalNumber]

private

Definition at line 90 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBXroc[RPotsTotalNumber]

private

Definition at line 81 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivPlanes[RPotsTotalNumber]

private

Definition at line 79 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivROCs[RPotsTotalNumber]

private

Definition at line 84 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const int CTPPSPixelDQMSource::IndexNotValid = 0

private

Definition at line 99 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::multHitsMax

private

Definition at line 123 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), dqmBeginRun(), and endRun().

int CTPPSPixelDQMSource::NArms =2

staticprivate

Definition at line 50 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

long int CTPPSPixelDQMSource::nEvents = 0

private

Definition at line 46 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), dqmBeginRun(), endRun(), and looper.Looper::loop().

int CTPPSPixelDQMSource::NPlaneBins = NplaneMAX*NRPotBinsInStation

staticprivate

Definition at line 66 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::NplaneMAX =6

staticprivate

Definition at line 53 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NROCsMAX = 6

staticprivate

Definition at line 54 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NRPotBinsInStation = RPn_last-RPn_first

staticprivate

Definition at line 65 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::NRPotsMAX =6

staticprivate

Definition at line 52 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

int CTPPSPixelDQMSource::NStationMAX =3

staticprivate

Definition at line 51 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::pixColMAX = 156

private

Definition at line 61 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

int CTPPSPixelDQMSource::pixRowMAX = 160

private

Definition at line 60 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::ROCSizeInX = pixRowMAX/2

private

Definition at line 62 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::ROCSizeInY = pixColMAX/3

private

Definition at line 63 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::RPindexValid[RPotsTotalNumber]

private

Definition at line 78 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::RPn_first = 3

staticprivate

Definition at line 55 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::RPn_last = 4

staticprivate

Definition at line 55 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::RPotsTotalNumber =NArms*NStationMAX*NRPotsMAX

staticprivate

Definition at line 76 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::RPstatus[NStationMAX][NRPotsMAX]

private

Definition at line 97 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

unsigned int CTPPSPixelDQMSource::rpStatusWord = 0x8000

private

Definition at line 96 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::StationStatus[NStationMAX]

private

Definition at line 98 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

CTPPSPixelIndices CTPPSPixelDQMSource::thePixIndices

private

Definition at line 59 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

edm::EDGetTokenT< edm::DetSetVector<CTPPSPixelCluster> > CTPPSPixelDQMSource::tokenCluster

private

Definition at line 48 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

edm::EDGetTokenT<edm::DetSetVector<CTPPSPixelDigi> > CTPPSPixelDQMSource::tokenDigi

private

Definition at line 47 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

unsigned int CTPPSPixelDQMSource::verbosity

private

Definition at line 45 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), dqmBeginRun(), and endRun().