CSCAFEBConnectAnalysis.cc

Go to the documentation of this file.

#include "OnlineDB/CSCCondDB/interface/CSCAFEBConnectAnalysis.h"
#include "OnlineDB/CSCCondDB/interface/CSCToAFEB.h"
#include "OnlineDB/CSCCondDB/interface/CSCOnlineDB.h"
#include "CondFormats/CSCObjects/interface/CSCobject.h"
#include "TMath.h"

CSCAFEBConnectAnalysis::CSCAFEBConnectAnalysis() {
  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();
}

void CSCAFEBConnectAnalysis::setup(const std::string& histoFileName) {
  hist_file = new TFile(histoFileName.c_str(), "RECREATE");
  hist_file->cd();
}

void CSCAFEBConnectAnalysis::bookForId(int flag, const int& idint, const std::string& idstring) {
  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
  }
}

void CSCAFEBConnectAnalysis::hf1ForId(std::map<int, TH1*>& mp, int flag, const int& id, float& x, float w) {
  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);
}

void CSCAFEBConnectAnalysis::hf2ForId(std::map<int, TH2*>& mp, int flag, const int& id, float& x, float& y, float w) {
  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);
}

/* Analyze the hits */
void CSCAFEBConnectAnalysis::analyze(const CSCWireDigiCollection& wirecltn) {
  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

void CSCAFEBConnectAnalysis::done() {
  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;
}