CSCAFEBConnectAnalysis Class Reference

#include <CSCAFEBConnectAnalysis.h>

Public Member Functions
void	analyze (const CSCWireDigiCollection &wirecltn)
	CSCAFEBConnectAnalysis ()
void	done ()
void	setup (const std::string &histoFileName)

Private Member Functions
void	bookForId (int flag, const int &idint, const std::string &ids)
void	hf1ForId (std::map< int, TH1 * > &mp, int flag, const int &id, float &x, float w)
void	hf2ForId (std::map< int, TH2 * > &mp, int flag, const int &id, float &x, float &y, float w)

Private Attributes
const CSCToAFEB	csctoafeb
	Layer, wire to AFEB, channel conversion. More...
TFile *	hist_file
	ROOT hist file. More...
std::map< int, int >	m_csc_list
	Maps. More...
std::map< int, std::vector< std::vector< float > > >	m_res_for_db
std::map< int, std::vector< int > >	m_wire_ev
std::map< int, TH1 * >	mh_Eff
std::map< int, TH2 * >	mh_FirstTime
std::map< int, TH1 * >	mh_LayerNmbPulses
	Histogram maps. More...
std::map< int, TH1 * >	mh_NonPairCrosstalk
std::map< int, TH1 * >	mh_PairCrosstalk
std::map< int, TH1 * >	mh_WireEff
std::map< int, TH1 * >	mh_WireNonPairCrosstalk
std::map< int, TH1 * >	mh_WirePairCrosstalk
int	nmbev
	Statistics. More...
int	nmbev_no_wire
int	nmblayers
std::vector< int >	nmbpulses
int	npulses
int	pulsed_layer

Detailed Description

Author

Definition at line 18 of file CSCAFEBConnectAnalysis.h.

Constructor & Destructor Documentation

CSCAFEBConnectAnalysis()

CSCAFEBConnectAnalysis::CSCAFEBConnectAnalysis

(

)

Definition at line 7 of file CSCAFEBConnectAnalysis.cc.

                                               {
  hist_file = nullptr;  // set to null
 
  nmbev = 0;
  nmbev_no_wire = 0;
  npulses = 0;
  nmblayers = 6;
  nmbpulses.resize(6, 0);
  pulsed_layer = 0;
 
  m_csc_list.clear();
  m_res_for_db.clear();
 
  mh_LayerNmbPulses.clear();
  mh_WireEff.clear();
  mh_WirePairCrosstalk.clear();
  mh_WireNonPairCrosstalk.clear();
  mh_Eff.clear();
  mh_PairCrosstalk.clear();
  mh_NonPairCrosstalk.clear();
 
  mh_FirstTime.clear();
}

References hist_file, m_csc_list, m_res_for_db, mh_Eff, mh_FirstTime, mh_LayerNmbPulses, mh_NonPairCrosstalk, mh_PairCrosstalk, mh_WireEff, mh_WireNonPairCrosstalk, mh_WirePairCrosstalk, nmbev, nmbev_no_wire, nmblayers, nmbpulses, npulses, and pulsed_layer.

Member Function Documentation

analyze()

void CSCAFEBConnectAnalysis::analyze

(

const CSCWireDigiCollection &

wirecltn

)

Store pulses per plane

Plot time bin of the first hit vs wires in layers

Mark given wire as wire having one hit

Fill maps

Fill efficiency map

Fill pair crosstalk map

Fill non-pair crosstalk map

Definition at line 129 of file CSCAFEBConnectAnalysis.cc.

                                                                          {
  std::ostringstream ss;
  std::map<int, std::vector<int> >::iterator viIt;
  std::map<int, std::vector<std::vector<float> > >::iterator vvfIt;
 
  int current_layer;
  float x, y;
  m_wire_ev.clear();
 
 
  nmbev++;
  pulsed_layer++;
  if (pulsed_layer == 7)
    pulsed_layer = 1;
  nmbpulses[pulsed_layer - 1] = nmbpulses[pulsed_layer - 1] + 1;
 
  //Anode wires
 
  CSCWireDigiCollection::DigiRangeIterator wiredetUnitIt;
  if (wirecltn.begin() == wirecltn.end())
    nmbev_no_wire++;
 
  if (wirecltn.begin() != wirecltn.end()) {
    for (wiredetUnitIt = wirecltn.begin(); wiredetUnitIt != wirecltn.end(); ++wiredetUnitIt) {
      const CSCDetId& id = (*wiredetUnitIt).first;
      const int idchamber = id.endcap() * 10000 + id.station() * 1000 + id.ring() * 100 + id.chamber();
      const int idlayer =
          id.endcap() * 100000 + id.station() * 10000 + id.ring() * 1000 + id.chamber() * 10 + id.layer();
 
      const int maxwire = csctoafeb.getMaxWire(id.station(), id.ring());
      std::vector<int> wireplane(maxwire, 0);
 
      const CSCWireDigiCollection::Range& range = (*wiredetUnitIt).second;
      for (CSCWireDigiCollection::const_iterator digiIt = range.first; digiIt != range.second; ++digiIt) {
        int iwire = (*digiIt).getWireGroup();
        if (iwire <= maxwire) {
          if (wireplane[iwire - 1] == 0) {
            wireplane[iwire - 1] = (*digiIt).getBeamCrossingTag() + 1;
 
            x = (id.layer() - 1) * maxwire + iwire;
            y = wireplane[iwire - 1];
            hf2ForId(mh_FirstTime, 1, idchamber, x, y, 1.0);
 
            wireplane[iwire - 1] = 1;
 
          }  // end if wireplane[iwire-1]==0
        }    // end if iwire<=csctoafeb.getMaxWire(id.station(),id.ring()
      }      // end of for digis in layer
 
      if (m_wire_ev.count(idlayer) == 0)
        m_wire_ev[idlayer] = wireplane;
 
    }  // end of cycle on detUnit
 
 
    for (viIt = m_wire_ev.begin(); viIt != m_wire_ev.end(); ++viIt) {
      const int idwirev = (*viIt).first;
      const std::vector<int> wiretemp = (*viIt).second;
      int nsize = 4;
      std::vector<float> zer(nsize, 0);
      vvfIt = m_res_for_db.find(idwirev);
      if (vvfIt == m_res_for_db.end()) {
        for (unsigned int j = 0; j < wiretemp.size(); j++)
          m_res_for_db[idwirev].push_back(zer);
        vvfIt = m_res_for_db.find(idwirev);
      }
      for (unsigned int i = 0; i < (*viIt).second.size(); i++) {
        current_layer = (*viIt).first / 10;
        current_layer = (*viIt).first - current_layer * 10;
        if (pulsed_layer == current_layer) {
          vvfIt->second[i][1] = vvfIt->second[i][1] + (*viIt).second[i];
        }
        if (pulsed_layer == 1 || pulsed_layer == 3 || pulsed_layer == 5)
          if (current_layer == (pulsed_layer + 1))
            vvfIt->second[i][2] = vvfIt->second[i][2] + (*viIt).second[i];
        if (pulsed_layer == 2 || pulsed_layer == 4 || pulsed_layer == 6)
          if (current_layer == (pulsed_layer - 1))
            vvfIt->second[i][2] = vvfIt->second[i][2] + (*viIt).second[i];
        if ((pulsed_layer > 2 && current_layer < 3) ||
            (pulsed_layer != 3 && pulsed_layer != 4 && current_layer > 2 && current_layer < 5) ||
            (pulsed_layer < 5 && current_layer > 4))
          vvfIt->second[i][3] = vvfIt->second[i][3] + (*viIt).second[i];
      }
    }  // end of adding hits to the maps
  }    // end of   if(wirecltn.begin() !=  wirecltn.end())
}  // end of void CSCAFEBConnectAnalysis

References csctoafeb, CSCDetId::endcap(), CSCToAFEB::getMaxWire(), hf2ForId(), mps_fire::i, dqmiolumiharvest::j, m_res_for_db, m_wire_ev, mh_FirstTime, nmbev, nmbev_no_wire, nmbpulses, pulsed_layer, FastTimerService_cff::range, relativeConstraints::ring, contentValuesCheck::ss, relativeConstraints::station, x, and y.

Referenced by CSCAFEBAnalyzer::analyze().

bookForId()

void CSCAFEBConnectAnalysis::bookForId ( int  flag, const int &  idint, const std::string &  ids  )

private

Definition at line 37 of file CSCAFEBConnectAnalysis.cc.

                                                                                            {
  hist_file->cd();
 
  std::ostringstream ss;
 
  if (flag == 1) {
    ss << idint << "_Anode_First_Time";
    mh_FirstTime[idint] = new TH2F(ss.str().c_str(), "", 675, 0.0, 675.0, 50, 0.0, 10.0);
    mh_FirstTime[idint]->GetXaxis()->SetTitle("(Layer-1)*Nwires+Wire");
    mh_FirstTime[idint]->GetYaxis()->SetTitle("Anode First Time Bin");
    mh_FirstTime[idint]->SetOption("BOX");
    ss.str("");  // clear
  }
 
  if (flag == 10) {
    ss << "Layer_Nmb_Pulses";
    mh_LayerNmbPulses[idint] = new TH1F(ss.str().c_str(), "", 7, 0.0, 7.0);
    mh_LayerNmbPulses[idint]->GetXaxis()->SetTitle("Layer");
    mh_LayerNmbPulses[idint]->GetYaxis()->SetTitle("Number of pulses");
    ss.str("");  // clear
  }
 
  if (flag == 101) {
    ss << idint << "_Anode_Wire_Eff";
    mh_WireEff[idint] = new TH1F(ss.str().c_str(), "", 675, 0.0, 675.0);
    mh_WireEff[idint]->GetXaxis()->SetTitle("(Layer-1)*Nwires+Wire");
    mh_WireEff[idint]->GetYaxis()->SetTitle("Efficiency");
    ss.str("");  // clear
  }
 
  if (flag == 102) {
    ss << idint << "_Anode_Eff";
    mh_Eff[idint] = new TH1F(ss.str().c_str(), "", 110, -0.05, 1.05);
    mh_Eff[idint]->GetXaxis()->SetTitle("Efficiency");
    mh_Eff[idint]->GetYaxis()->SetTitle("Entries");
    ss.str("");  // clear
  }
 
  if (flag == 201) {
    ss << idint << "_Anode_Wire_Pair_Layer_Crosstalk";
    mh_WirePairCrosstalk[idint] = new TH1F(ss.str().c_str(), "", 675, 0.0, 675.0);
    mh_WirePairCrosstalk[idint]->GetXaxis()->SetTitle("(Layer-1)*Nwires+Wire");
    mh_WirePairCrosstalk[idint]->GetYaxis()->SetTitle("Probability");
    ss.str("");  // clear
  }
 
  if (flag == 202) {
    ss << idint << "_Anode_Pair_Layer_Crosstalk";
    mh_PairCrosstalk[idint] = new TH1F(ss.str().c_str(), "", 70, -0.05, 0.3);
    mh_PairCrosstalk[idint]->GetXaxis()->SetTitle("Probability");
    mh_PairCrosstalk[idint]->GetYaxis()->SetTitle("Entries");
    ss.str("");  // clear
  }
 
  if (flag == 301) {
    ss << idint << "_Anode_Wire_NonPair_Layer_Crosstalk";
    mh_WireNonPairCrosstalk[idint] = new TH1F(ss.str().c_str(), "", 675, 0.0, 675.0);
    mh_WireNonPairCrosstalk[idint]->GetXaxis()->SetTitle("(Layer-1)*Nwires+Wire");
    mh_WireNonPairCrosstalk[idint]->GetYaxis()->SetTitle("Probability");
    ss.str("");  // clear
  }
 
  if (flag == 302) {
    ss << idint << "_Anode_NonPair_Layer_Crosstalk";
    mh_NonPairCrosstalk[idint] = new TH1F(ss.str().c_str(), "", 70, -0.05, 0.3);
    mh_NonPairCrosstalk[idint]->GetXaxis()->SetTitle("Probability");
    mh_NonPairCrosstalk[idint]->GetYaxis()->SetTitle("Entries");
    ss.str("");  // clear
  }
}

References RemoveAddSevLevel::flag, hist_file, mh_Eff, mh_FirstTime, mh_LayerNmbPulses, mh_NonPairCrosstalk, mh_PairCrosstalk, mh_WireEff, mh_WireNonPairCrosstalk, mh_WirePairCrosstalk, and contentValuesCheck::ss.

Referenced by hf1ForId(), and hf2ForId().

done()

void CSCAFEBConnectAnalysis::done

(

)

Fill number of pulses per layer, normalize the non-pair crosstalk, make overal normalization, fill the plots

Fill efficiency plot

Fill pair crosstalk

Fill nonpair crosstalk

Definition at line 223 of file CSCAFEBConnectAnalysis.cc.

                                  {
  float x, y;
 
  //This is for DB transfer
  //  CSCobject *cn = new CSCobject();
  //  condbon *dbon = new condbon();
 
  std::map<int, int>::iterator intIt;
  std::map<int, std::vector<std::vector<float> > >::iterator vvfIt;
  std::cout << "Events analyzed  " << nmbev << std::endl;
  std::cout << "Events no anodes " << nmbev_no_wire << std::endl << std::endl;
 
  std::cout << "Number of pulses per layer" << std::endl;
  for (int i = 0; i < nmblayers; i++)
    std::cout << " " << nmbpulses[i];
  std::cout << "\n" << std::endl;
 
  //  std::vector<float> inputx;
  //  std::vector<float> inputy;
 
 
  for (int i = 0; i < nmblayers; i++) {
    x = i + 1;
    y = nmbpulses[i];
    hf1ForId(mh_LayerNmbPulses, 10, 1, x, y);
  }
  for (vvfIt = m_res_for_db.begin(); vvfIt != m_res_for_db.end(); ++vvfIt) {
    int idlayer = (*vvfIt).first;
    int idchmb = idlayer / 10;
    int layer = idlayer - idchmb * 10;
    for (unsigned int i = 0; i < (*vvfIt).second.size(); i++) {
      (*vvfIt).second[i][0] = nmbpulses[layer - 1];
      (*vvfIt).second[i][3] = (*vvfIt).second[i][3] / 4.0;
      (*vvfIt).second[i][1] = (*vvfIt).second[i][1] / (*vvfIt).second[i][0];
      (*vvfIt).second[i][2] = (*vvfIt).second[i][2] / (*vvfIt).second[i][0];
      (*vvfIt).second[i][3] = (*vvfIt).second[i][3] / (*vvfIt).second[i][0];
 
      x = (layer - 1) * (*vvfIt).second.size() + (i + 1);
 
      y = (*vvfIt).second[i][1];
      hf1ForId(mh_WireEff, 101, idchmb, x, y);
      hf1ForId(mh_Eff, 102, idchmb, y, 1.0);
 
      y = (*vvfIt).second[i][2];
      hf1ForId(mh_WirePairCrosstalk, 201, idchmb, x, y);
      hf1ForId(mh_PairCrosstalk, 202, idchmb, y, 1.0);
 
      y = (*vvfIt).second[i][3];
      hf1ForId(mh_WireNonPairCrosstalk, 301, idchmb, x, y);
      hf1ForId(mh_NonPairCrosstalk, 302, idchmb, y, 1.0);
    }
  }
  std::cout << "Size of map for DB " << m_res_for_db.size() << std::endl;
 
  std::cout << "The following CSCs will go to DB" << std::endl << std::endl;
  for (vvfIt = m_res_for_db.begin(); vvfIt != m_res_for_db.end(); ++vvfIt) {
    int idchmb = (*vvfIt).first / 10;
    if (m_csc_list.count(idchmb) == 0)
      m_csc_list[idchmb] = 0;
    if (m_csc_list.count(idchmb) > 0)
      m_csc_list[idchmb] = m_csc_list[idchmb] + (*vvfIt).second.size();
  }
  int count = 0;
  for (intIt = m_csc_list.begin(); intIt != m_csc_list.end(); ++intIt) {
    count++;
    std::cout << count << " "
              << " CSC " << (*intIt).first << "  " << (*intIt).second << std::endl;
  }
  std::cout << std::endl;
 
  /*
 
  for(vvfIt=m_res_for_db.begin(); vvfIt!=m_res_for_db.end(); 
      ++vvfIt) {
      int idlayer=(*vvfIt).first;
      int size = (*vvfIt).second.size();
      cn->obj[idlayer].resize(size);
      for (unsigned int i=0;i<(*vvfIt).second.size();i++) { 
    std::cout<<idlayer<<" "<<i+1<<"    ";   
        for(int j=0;j<4;j++) std::cout<< (*vvfIt).second[i][j]<<" ";
    std::cout<<std::endl;
 
    cn->obj[idlayer][i].resize(4);
    cn->obj[idlayer][i][0] = (*vvfIt).second[i][0];
    cn->obj[idlayer][i][1] = (*vvfIt).second[i][1];
    cn->obj[idlayer][i][2] = (*vvfIt).second[i][2];
    cn->obj[idlayer][i][3] = (*vvfIt).second[i][3];
      }
  }
 
 
  dbon->cdbon_last_run("afeb_thresholds",&run);
  std::cout<<"Last AFEB thresholds run "<<run<<" for run file "<<myname<<" saved "<<myTime<<std::endl;
  if(debug) dbon->cdbon_write(cn,"afeb_thresholds",run+1,myTime);
*/
 
  if (hist_file != nullptr) {  // if there was a histogram file...
    hist_file->Write();        // write out the histrograms
    delete hist_file;          // close and delete the file
    hist_file = nullptr;       // set to zero to clean up
    std::cout << "Hist. file was closed\n";
  }
  std::cout << " End of CSCAFEBConnectAnalysis" << std::endl;
}

References KineDebug3::count(), gather_cfg::cout, hf1ForId(), hist_file, mps_fire::i, m_csc_list, m_res_for_db, mh_Eff, mh_LayerNmbPulses, mh_NonPairCrosstalk, mh_PairCrosstalk, mh_WireEff, mh_WireNonPairCrosstalk, mh_WirePairCrosstalk, nmbev, nmbev_no_wire, nmblayers, nmbpulses, x, and y.

Referenced by CSCAFEBAnalyzer::endJob().

hf1ForId()

void CSCAFEBConnectAnalysis::hf1ForId ( std::map< int, TH1 * > &  mp, int  flag, const int &  id, float &  x, float  w  )

private

Definition at line 108 of file CSCAFEBConnectAnalysis.cc.

                                                                                                       {
  std::map<int, TH1*>::iterator h;
  h = mp.find(id);
  if (h == mp.end()) {
    bookForId(flag, id, "");
    h = mp.find(id);
  }
  h->second->Fill(x, w);
}

References bookForId(), RemoveAddSevLevel::flag, h, w, and x.

Referenced by done().

hf2ForId()

void CSCAFEBConnectAnalysis::hf2ForId ( std::map< int, TH2 * > &  mp, int  flag, const int &  id, float &  x, float &  y, float  w  )

private

Definition at line 118 of file CSCAFEBConnectAnalysis.cc.

                                                                                                                 {
  std::map<int, TH2*>::iterator h;
  h = mp.find(id);
  if (h == mp.end()) {
    bookForId(flag, id, "");
    h = mp.find(id);
  }
  h->second->Fill(x, y, w);
}

References bookForId(), RemoveAddSevLevel::flag, h, w, x, and y.

Referenced by analyze().

setup()

void CSCAFEBConnectAnalysis::setup

(

const std::string &

histoFileName

)

open the histogram file

Definition at line 31 of file CSCAFEBConnectAnalysis.cc.

                                                                 {
  hist_file = new TFile(histoFileName.c_str(), "RECREATE");
  hist_file->cd();
}

References hist_file, and mix_2012_lumiLevel_15_20_50ns_PoissonOOTPU_cfi::histoFileName.

Referenced by CSCAFEBAnalyzer::CSCAFEBAnalyzer().

Member Data Documentation

csctoafeb

const CSCToAFEB CSCAFEBConnectAnalysis::csctoafeb

private

Layer, wire to AFEB, channel conversion.

Definition at line 44 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze().

hist_file

TFile* CSCAFEBConnectAnalysis::hist_file

private

ROOT hist file.

Definition at line 47 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), done(), and setup().

m_csc_list

std::map<int, int> CSCAFEBConnectAnalysis::m_csc_list

private

Maps.

Definition at line 39 of file CSCAFEBConnectAnalysis.h.

Referenced by CSCAFEBConnectAnalysis(), and done().

m_res_for_db

std::map<int, std::vector<std::vector<float> > > CSCAFEBConnectAnalysis::m_res_for_db

private

Definition at line 41 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze(), CSCAFEBConnectAnalysis(), and done().

m_wire_ev

std::map<int, std::vector<int> > CSCAFEBConnectAnalysis::m_wire_ev

private

Definition at line 40 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze().

mh_Eff

std::map<int, TH1*> CSCAFEBConnectAnalysis::mh_Eff

private

Definition at line 52 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), and done().

mh_FirstTime

std::map<int, TH2*> CSCAFEBConnectAnalysis::mh_FirstTime

private

Definition at line 57 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze(), bookForId(), and CSCAFEBConnectAnalysis().

mh_LayerNmbPulses

std::map<int, TH1*> CSCAFEBConnectAnalysis::mh_LayerNmbPulses

private

Histogram maps.

Definition at line 50 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), and done().

mh_NonPairCrosstalk

std::map<int, TH1*> CSCAFEBConnectAnalysis::mh_NonPairCrosstalk

private

Definition at line 56 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), and done().

mh_PairCrosstalk

std::map<int, TH1*> CSCAFEBConnectAnalysis::mh_PairCrosstalk

private

Definition at line 54 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), and done().

mh_WireEff

std::map<int, TH1*> CSCAFEBConnectAnalysis::mh_WireEff

private

Definition at line 51 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), and done().

mh_WireNonPairCrosstalk

std::map<int, TH1*> CSCAFEBConnectAnalysis::mh_WireNonPairCrosstalk

private

Definition at line 55 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), and done().

mh_WirePairCrosstalk

std::map<int, TH1*> CSCAFEBConnectAnalysis::mh_WirePairCrosstalk

private

Definition at line 53 of file CSCAFEBConnectAnalysis.h.

Referenced by bookForId(), CSCAFEBConnectAnalysis(), and done().

nmbev

int CSCAFEBConnectAnalysis::nmbev

private

Statistics.

Definition at line 30 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze(), CSCAFEBConnectAnalysis(), and done().

nmbev_no_wire

int CSCAFEBConnectAnalysis::nmbev_no_wire

private

Definition at line 31 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze(), CSCAFEBConnectAnalysis(), and done().

nmblayers

int CSCAFEBConnectAnalysis::nmblayers

private

Definition at line 33 of file CSCAFEBConnectAnalysis.h.

Referenced by CSCAFEBConnectAnalysis(), and done().

nmbpulses

std::vector<int> CSCAFEBConnectAnalysis::nmbpulses

private

Definition at line 35 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze(), CSCAFEBConnectAnalysis(), and done().

npulses

int CSCAFEBConnectAnalysis::npulses

private

Definition at line 32 of file CSCAFEBConnectAnalysis.h.

Referenced by CSCAFEBConnectAnalysis().

pulsed_layer

int CSCAFEBConnectAnalysis::pulsed_layer

private

Definition at line 34 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze(), and CSCAFEBConnectAnalysis().