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

(

)

Definition at line 7 of file CSCAFEBConnectAnalysis.cc.

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.

                                               {

hist_file=0; // 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();

}

Member Function Documentation

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 138 of file CSCAFEBConnectAnalysis.cc.

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

Referenced by CSCAFEBAnalyzer::analyze().

                                                                          {

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::bookForId ( int  flag, const int &  idint, const std::string &  ids  )

private

Definition at line 39 of file CSCAFEBConnectAnalysis.cc.

References hist_file, mh_Eff, mh_FirstTime, mh_LayerNmbPulses, mh_NonPairCrosstalk, mh_PairCrosstalk, mh_WireEff, mh_WireNonPairCrosstalk, and mh_WirePairCrosstalk.

Referenced by hf1ForId(), and hf2ForId().

                                                                         {
  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::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 243 of file CSCAFEBConnectAnalysis.cc.

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().

                                  {

  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!=0) { // if there was a histogram file...
    hist_file->Write(); // write out the histrograms
    delete hist_file; // close and delete the file
    hist_file=0; // set to zero to clean up
    std::cout << "Hist. file was closed\n";
  }
  std::cout<<" End of CSCAFEBConnectAnalysis"<<std::endl;  
}

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

private

Definition at line 112 of file CSCAFEBConnectAnalysis.cc.

References bookForId(), and h.

Referenced by done().

                                  {

  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  )

private

Definition at line 124 of file CSCAFEBConnectAnalysis.cc.

References bookForId(), and h.

Referenced by analyze().

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

void CSCAFEBConnectAnalysis::setup

(

const std::string &

histoFileName

)

open the histogram file

Definition at line 33 of file CSCAFEBConnectAnalysis.cc.

References hist_file.

Referenced by CSCAFEBAnalyzer::CSCAFEBAnalyzer().

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

Member Data Documentation

const CSCToAFEB CSCAFEBConnectAnalysis::csctoafeb

private

Layer, wire to AFEB, channel conversion.

Definition at line 46 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze().

TFile* CSCAFEBConnectAnalysis::hist_file

private

ROOT hist file.

Definition at line 49 of file CSCAFEBConnectAnalysis.h.

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

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

private

Maps.

Definition at line 41 of file CSCAFEBConnectAnalysis.h.

Referenced by CSCAFEBConnectAnalysis(), and done().

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

private

Definition at line 43 of file CSCAFEBConnectAnalysis.h.

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

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

private

Definition at line 42 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze().

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

private

Definition at line 54 of file CSCAFEBConnectAnalysis.h.

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

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

private

Definition at line 59 of file CSCAFEBConnectAnalysis.h.

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

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

private

Histogram maps.

Definition at line 52 of file CSCAFEBConnectAnalysis.h.

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

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

private

Definition at line 58 of file CSCAFEBConnectAnalysis.h.

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

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

private

Definition at line 56 of file CSCAFEBConnectAnalysis.h.

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

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

private

Definition at line 53 of file CSCAFEBConnectAnalysis.h.

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

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

private

Definition at line 57 of file CSCAFEBConnectAnalysis.h.

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

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

private

Definition at line 55 of file CSCAFEBConnectAnalysis.h.

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

int CSCAFEBConnectAnalysis::nmbev

private

Statistics.

Definition at line 32 of file CSCAFEBConnectAnalysis.h.

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

int CSCAFEBConnectAnalysis::nmbev_no_wire

private

Definition at line 33 of file CSCAFEBConnectAnalysis.h.

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

int CSCAFEBConnectAnalysis::nmblayers

private

Definition at line 35 of file CSCAFEBConnectAnalysis.h.

Referenced by CSCAFEBConnectAnalysis(), and done().

std::vector<int> CSCAFEBConnectAnalysis::nmbpulses

private

Definition at line 37 of file CSCAFEBConnectAnalysis.h.

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

int CSCAFEBConnectAnalysis::npulses

private

Definition at line 34 of file CSCAFEBConnectAnalysis.h.

Referenced by CSCAFEBConnectAnalysis().

int CSCAFEBConnectAnalysis::pulsed_layer

private

Definition at line 36 of file CSCAFEBConnectAnalysis.h.

Referenced by analyze(), and CSCAFEBConnectAnalysis().