Inheritance diagram for SiStripHitEffFromCalibTree:

Public Member Functions
	SiStripHitEffFromCalibTree (const edm::ParameterSet &)

	~SiStripHitEffFromCalibTree ()

Public Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
	ConditionDBWriter (const edm::ParameterSet &iConfig)

virtual	~ConditionDBWriter ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Member Functions
virtual void	algoAnalyze (const edm::Event &e, const edm::EventSetup &c)

virtual void	algoBeginJob ()

virtual void	algoEndJob ()

float	calcPhi (float x, float y)

SiStripBadStrip *	getNewObject ()

void	makeHotColdMaps ()

void	makeSQLite ()

void	makeSummary ()

void	makeTKMap ()

void	SetBadComponents (int i, int component, SiStripQuality::BadComponent &BC, std::stringstream ssV[4][19], int NBadComponent[4][19][4])

void	totalStatistics ()

Private Attributes
unsigned int	_bunchx

float	_ResXSig

int	alllayerfound [35]

int	alllayertotal [35]

map< unsigned int, double >	BadModules

TTree *	CalibTree

TFile *	CalibTreeFile

TString	CalibTreeFilename

unsigned int	doSummary

edm::FileInPath	FileInPath_

edm::Service< TFileService >	fs

int	goodlayerfound [35]

int	goodlayertotal [35]

vector< hit >	hits [23]

vector< TH2F * >	HotColdMaps

int	layerfound [23]

int	layertotal [23]

map< unsigned int, pair < unsigned int, unsigned int > >	modCounter [23]

unsigned int	nModsMin

SiStripQuality *	quality_

SiStripDetInfoFileReader *	reader

float	threshold

TrackerMap *	tkmap

TrackerMap *	tkmapbad

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from ConditionDBWriter< SiStripBadStrip >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...

void	storeOnDbNow ()

cond::Time_t	timeOfLastIOV ()

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 94 of file SiStripHitEffFromCalibTree.cc.

Constructor & Destructor Documentation

SiStripHitEffFromCalibTree::SiStripHitEffFromCalibTree ( const edm::ParameterSet & conf )

explicit

Definition at line 139 of file SiStripHitEffFromCalibTree.cc.

References _bunchx, _ResXSig, CalibTreeFilename, doSummary, FileInPath_, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), nModsMin, quality_, reader, OfflineSiStripConfigDb_cff::SiStripDetInfoFileReader, AlCaHarvesting_cff::SiStripQuality, and threshold.

                                                                                   :
   ConditionDBWriter<SiStripBadStrip>(conf),
   FileInPath_("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat")
 {
   CalibTreeFilename = conf.getParameter<std::string>("CalibTreeFilename"); 
   threshold = conf.getParameter<double>("Threshold");
   nModsMin = conf.getParameter<int>("nModsMin");
   doSummary = conf.getParameter<int>("doSummary");
   _ResXSig = conf.getUntrackedParameter<double>("ResXSig",-1);
   _bunchx = conf.getUntrackedParameter<int>("BunchCrossing",0);
   reader = new SiStripDetInfoFileReader(FileInPath_.fullPath());
   
   quality_ = new SiStripQuality;
 }

SiStripHitEffFromCalibTree::~SiStripHitEffFromCalibTree ( )

Definition at line 154 of file SiStripHitEffFromCalibTree.cc.

154 { }

Member Function Documentation

void SiStripHitEffFromCalibTree::algoAnalyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

privatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 166 of file SiStripHitEffFromCalibTree.cc.

References _bunchx, _ResXSig, a, accept(), alllayerfound, alllayertotal, CalibTree, CalibTreeFile, CalibTreeFilename, gather_cfg::cout, SiStripBadStrip::decode(), cond::rpcobgas::detid, end, edm::EventID::event(), SiStripQuality::getBadComponentList(), SiStripBadStrip::getDataVectorBegin(), SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), goodlayerfound, goodlayertotal, hits, i, edm::EventBase::id(), hit::id, errorMatrix2Lands_multiChannel::id, j, gen::k, prof2calltree::l, TOBDetId::layer(), TIBDetId::layer(), edm::EventBase::luminosityBlock(), makeHotColdMaps(), makeSQLite(), makeSummary(), makeTKMap(), modCounter, nevents, nModsMin, quality_, SiStripBadStrip::data::range, reader, edm::EventID::run(), SetBadComponents(), TIDDetId::side(), TECDetId::side(), DetId::subdetId(), SiStripDetId::TEC, threshold, SiStripDetId::TIB, SiStripDetId::TID, edm::EventBase::time(), SiStripDetId::TOB, totalStatistics(), edm::Timestamp::value(), TIDDetId::wheel(), TECDetId::wheel(), x, hit::x, detailsBasic3DVector::y, hit::y, detailsBasic3DVector::z, and hit::z.

                                                                                       {
   //Open the ROOT Calib Tree
   CalibTreeFile = TFile::Open(CalibTreeFilename,"READ");
   CalibTreeFile->cd("anEff"); 
   CalibTree = (TTree*)(gDirectory->Get("traj")) ;
   TLeaf* BadLf = CalibTree->GetLeaf("ModIsBad");
   TLeaf* sistripLf = CalibTree->GetLeaf("SiStripQualBad");
   TLeaf* idLf = CalibTree->GetLeaf("Id");
   TLeaf* acceptLf = CalibTree->GetLeaf("withinAcceptance");
   TLeaf* layerLf = CalibTree->GetLeaf("layer");
   TLeaf* nHitsLf = CalibTree->GetLeaf("nHits");
   TLeaf* xLf = CalibTree->GetLeaf("TrajGlbX");
   TLeaf* yLf = CalibTree->GetLeaf("TrajGlbY");
   TLeaf* zLf = CalibTree->GetLeaf("TrajGlbZ");
   TLeaf* ResXSigLf = CalibTree->GetLeaf("ResXSig");
   TLeaf* BunchLf(0);
   for(int l=0; l < 35; l++) {
     goodlayertotal[l] = 0;
     goodlayerfound[l] = 0;
     alllayertotal[l] = 0;
     alllayerfound[l] = 0;
   }
   if(_bunchx != 0) {
     BunchLf = CalibTree->GetLeaf("bunchx");
   }
   int nevents = CalibTree->GetEntries();
   cout << "Successfully loaded analyze function with " << nevents << " events!\n";
   cout << "A module is bad if efficiency < " << threshold << " and has at least " << nModsMin << " nModsMin." << endl;
 
   //Loop through all of the events
   for(int j =0; j < nevents; j++) {
     CalibTree->GetEvent(j);
     unsigned int isBad = (unsigned int)BadLf->GetValue();
     unsigned int quality = (unsigned int)sistripLf->GetValue();
     unsigned int id = (unsigned int)idLf->GetValue();
     unsigned int accept = (unsigned int)acceptLf->GetValue();
     unsigned int layer = (unsigned int)layerLf->GetValue();
     unsigned int nHits = (unsigned int)nHitsLf->GetValue();
     double x = xLf->GetValue();
     double y = yLf->GetValue();
     double z = zLf->GetValue();
     double resxsig = ResXSigLf->GetValue();
     bool badquality = false;
     if(_bunchx != 0) {
       if(_bunchx != BunchLf->GetValue()) continue;
     }
     //We have two things we want to do, both an XY color plot, and the efficiency measurement
     //First, ignore anything that isn't in acceptance and isn't good quality
     
     //if(quality == 1 || accept != 1 || nHits < 8) continue;
     if(accept != 1 || nHits < 8) continue;
     if(quality == 1) badquality = true;
     
     //Now that we have a good event, we need to look at if we expected it or not, and the location
     //if we didn't
     //Fill the missing hit information first
     bool badflag = false;
     if(_ResXSig < 0) {
       if(isBad == 1) badflag = true;
     }
     else {
       if(isBad == 1 || resxsig > _ResXSig) badflag = true;
     }
     if(badflag && !badquality) {   
       hit temphit;         
       temphit.x = x;
       temphit.y = y;
       temphit.z = z;
       temphit.id = id;
       hits[layer].push_back(temphit);
     } 
     pair<unsigned int, unsigned int> newgoodpair (1,1);
     pair<unsigned int, unsigned int> newbadpair (1,0);
     //First, figure out if the module already exists in the map of maps
     map< unsigned int, pair< unsigned int, unsigned int> >::iterator it = modCounter[layer].find(id);
     if(!badquality) {
       if(it == modCounter[layer].end()) {
         if(badflag) modCounter[layer][id] = newbadpair;
         else modCounter[layer][id] = newgoodpair;
       }
       else {
         ((*it).second.first)++;
         if(!badflag) ((*it).second.second)++;
       }
       //Have to do the decoding for which side to go on (ugh)
       if(layer <= 10) {
         if(!badflag) goodlayerfound[layer]++;
         goodlayertotal[layer]++;
       }
       else if(layer > 10 && layer < 14) {
         if( ((id>>13)&0x3) == 1) {
       if(!badflag) goodlayerfound[layer]++;
           goodlayertotal[layer]++;
     }
     else if( ((id>>13)&0x3) == 2) {
       if(!badflag) goodlayerfound[layer+3]++;
           goodlayertotal[layer+3]++;
     }
       }
       else if(layer > 13 && layer <= 22) {
         if( ((id>>18)&0x3) == 1) {
       if(!badflag) goodlayerfound[layer+3]++;
           goodlayertotal[layer+3]++;
     }
     else if( ((id>>18)&0x3) == 2) {
       if(!badflag) goodlayerfound[layer+12]++;
           goodlayertotal[layer+12]++;
     }
       } 
     }
     //Do the one where we don't exclude bad modules!
     if(layer <= 10) {
       if(!badflag) alllayerfound[layer]++;
       alllayertotal[layer]++;
     }
     else if(layer > 10 && layer < 14) {
       if( ((id>>13)&0x3) == 1) {
     if(!badflag) alllayerfound[layer]++;
         alllayertotal[layer]++;
       }
       else if( ((id>>13)&0x3) == 2) {
         if(!badflag) alllayerfound[layer+3]++;
         alllayertotal[layer+3]++;
       }
     }
     else if(layer > 13 && layer <= 22) {
       if( ((id>>18)&0x3) == 1) {
         if(!badflag) alllayerfound[layer+3]++;
         alllayertotal[layer+3]++;
       }
       else if( ((id>>18)&0x3) == 2) {
         if(!badflag) alllayerfound[layer+12]++;
         alllayertotal[layer+12]++;
       }
     }  
     //At this point, both of our maps are loaded with the correct information
   }
   //CalibTreeFile->Close();
   makeHotColdMaps();
   makeTKMap();
   makeSQLite();
   totalStatistics();
   makeSummary();
   
   //try to write out what's in the quality record
   int NTkBadComponent[4]; //k: 0=BadModule, 1=BadFiber, 2=BadApv, 3=BadStrips
   int NBadComponent[4][19][4];  
   //legend: NBadComponent[i][j][k]= SubSystem i, layer/disk/wheel j, BadModule/Fiber/Apv k
   //     i: 0=TIB, 1=TID, 2=TOB, 3=TEC
   //     k: 0=BadModule, 1=BadFiber, 2=BadApv, 3=BadStrips
   std::stringstream ssV[4][19];
   
   for(int i=0;i<4;++i){
     NTkBadComponent[i]=0;
     for(int j=0;j<19;++j){
       ssV[i][j].str("");
       for(int k=0;k<4;++k)
         NBadComponent[i][j][k]=0;
     }
   }
  
  
   std::vector<SiStripQuality::BadComponent> BC = quality_->getBadComponentList();
    
   for (size_t i=0;i<BC.size();++i){
      
     //&&&&&&&&&&&&&
     //Full Tk
     //&&&&&&&&&&&&&
  
     if (BC[i].BadModule) 
       NTkBadComponent[0]++;
     if (BC[i].BadFibers) 
       NTkBadComponent[1]+= ( (BC[i].BadFibers>>2)&0x1 )+ ( (BC[i].BadFibers>>1)&0x1 ) + ( (BC[i].BadFibers)&0x1 );
     if (BC[i].BadApvs)
       NTkBadComponent[2]+= ( (BC[i].BadApvs>>5)&0x1 )+ ( (BC[i].BadApvs>>4)&0x1 ) + ( (BC[i].BadApvs>>3)&0x1 ) + 
     ( (BC[i].BadApvs>>2)&0x1 )+ ( (BC[i].BadApvs>>1)&0x1 ) + ( (BC[i].BadApvs)&0x1 );
  
     //&&&&&&&&&&&&&&&&&
     //Single SubSystem
     //&&&&&&&&&&&&&&&&&
  
     int component;
     SiStripDetId a(BC[i].detid);
     if ( a.subdetId() == SiStripDetId::TIB ){
       //&&&&&&&&&&&&&&&&&
       //TIB
       //&&&&&&&&&&&&&&&&&
        
       component=TIBDetId(BC[i].detid).layer();
       SetBadComponents(0, component, BC[i], ssV, NBadComponent);          
  
     } else if ( a.subdetId() == SiStripDetId::TID ) {
       //&&&&&&&&&&&&&&&&&
       //TID
       //&&&&&&&&&&&&&&&&&
  
       component=TIDDetId(BC[i].detid).side()==2?TIDDetId(BC[i].detid).wheel():TIDDetId(BC[i].detid).wheel()+3;
       SetBadComponents(1, component, BC[i], ssV, NBadComponent);          
  
     } else if ( a.subdetId() == SiStripDetId::TOB ) {
       //&&&&&&&&&&&&&&&&&
       //TOB
       //&&&&&&&&&&&&&&&&&
  
       component=TOBDetId(BC[i].detid).layer();
       SetBadComponents(2, component, BC[i], ssV, NBadComponent);          
  
     } else if ( a.subdetId() == SiStripDetId::TEC ) {
       //&&&&&&&&&&&&&&&&&
       //TEC
       //&&&&&&&&&&&&&&&&&
  
       component=TECDetId(BC[i].detid).side()==2?TECDetId(BC[i].detid).wheel():TECDetId(BC[i].detid).wheel()+9;
       SetBadComponents(3, component, BC[i], ssV, NBadComponent);          
  
     }    
   }
  
   //&&&&&&&&&&&&&&&&&&
   // Single Strip Info
   //&&&&&&&&&&&&&&&&&&
   float percentage=0;
  
   SiStripQuality::RegistryIterator rbegin = quality_->getRegistryVectorBegin();
   SiStripQuality::RegistryIterator rend   = quality_->getRegistryVectorEnd();
    
   for (SiStripBadStrip::RegistryIterator rp=rbegin; rp != rend; ++rp) {
     unsigned int detid=rp->detid;
  
     int subdet=-999; int component=-999;
     SiStripDetId a(detid);
     if ( a.subdetId() == 3 ){
       subdet=0;
       component=TIBDetId(detid).layer();
     } else if ( a.subdetId() == 4 ) {
       subdet=1;
       component=TIDDetId(detid).side()==2?TIDDetId(detid).wheel():TIDDetId(detid).wheel()+3;
     } else if ( a.subdetId() == 5 ) {
       subdet=2;
       component=TOBDetId(detid).layer();
     } else if ( a.subdetId() == 6 ) {
       subdet=3;
       component=TECDetId(detid).side()==2?TECDetId(detid).wheel():TECDetId(detid).wheel()+9;
     } 
  
     SiStripQuality::Range sqrange = SiStripQuality::Range( quality_->getDataVectorBegin()+rp->ibegin , quality_->getDataVectorBegin()+rp->iend );
      
     percentage=0;
     for(int it=0;it<sqrange.second-sqrange.first;it++){
       unsigned int range=quality_->decode( *(sqrange.first+it) ).range;
       NTkBadComponent[3]+=range;
       NBadComponent[subdet][0][3]+=range;
       NBadComponent[subdet][component][3]+=range;
       percentage+=range;
     }
     if(percentage!=0)
       percentage/=128.*reader->getNumberOfApvsAndStripLength(detid).first;
     if(percentage>1)
       edm::LogError("SiStripQualityStatistics") <<  "PROBLEM detid " << detid << " value " << percentage<< std::endl; 
   }
   //&&&&&&&&&&&&&&&&&&
   // printout
   //&&&&&&&&&&&&&&&&&&
  
   cout << "\n-----------------\nNew IOV starting from run " <<   e.id().run() << " event " << e.id().event() << " lumiBlock " << e.luminosityBlock() << " time " << e.time().value()  << "\n-----------------\n";
   cout << "\n-----------------\nGlobal Info\n-----------------";
   cout << "\nBadComponent \t    Modules \tFibers \tApvs\tStrips\n----------------------------------------------------------------";
   cout << "\nTracker:\t\t"<<NTkBadComponent[0]<<"\t"<<NTkBadComponent[1]<<"\t"<<NTkBadComponent[2]<<"\t"<<NTkBadComponent[3];
   cout << endl;
   cout << "\nTIB:\t\t\t"<<NBadComponent[0][0][0]<<"\t"<<NBadComponent[0][0][1]<<"\t"<<NBadComponent[0][0][2]<<"\t"<<NBadComponent[0][0][3];
   cout << "\nTID:\t\t\t"<<NBadComponent[1][0][0]<<"\t"<<NBadComponent[1][0][1]<<"\t"<<NBadComponent[1][0][2]<<"\t"<<NBadComponent[1][0][3];
   cout << "\nTOB:\t\t\t"<<NBadComponent[2][0][0]<<"\t"<<NBadComponent[2][0][1]<<"\t"<<NBadComponent[2][0][2]<<"\t"<<NBadComponent[2][0][3];
   cout << "\nTEC:\t\t\t"<<NBadComponent[3][0][0]<<"\t"<<NBadComponent[3][0][1]<<"\t"<<NBadComponent[3][0][2]<<"\t"<<NBadComponent[3][0][3];
   cout << "\n";
  
   for (int i=1;i<5;++i)
     cout << "\nTIB Layer " << i   << " :\t\t"<<NBadComponent[0][i][0]<<"\t"<<NBadComponent[0][i][1]<<"\t"<<NBadComponent[0][i][2]<<"\t"<<NBadComponent[0][i][3];
   cout << "\n";
   for (int i=1;i<4;++i)
     cout << "\nTID+ Disk " << i   << " :\t\t"<<NBadComponent[1][i][0]<<"\t"<<NBadComponent[1][i][1]<<"\t"<<NBadComponent[1][i][2]<<"\t"<<NBadComponent[1][i][3];
   for (int i=4;i<7;++i)
     cout << "\nTID- Disk " << i-3 << " :\t\t"<<NBadComponent[1][i][0]<<"\t"<<NBadComponent[1][i][1]<<"\t"<<NBadComponent[1][i][2]<<"\t"<<NBadComponent[1][i][3];
   cout << "\n";
   for (int i=1;i<7;++i)
     cout << "\nTOB Layer " << i   << " :\t\t"<<NBadComponent[2][i][0]<<"\t"<<NBadComponent[2][i][1]<<"\t"<<NBadComponent[2][i][2]<<"\t"<<NBadComponent[2][i][3];
   cout << "\n";
   for (int i=1;i<10;++i)
     cout << "\nTEC+ Disk " << i   << " :\t\t"<<NBadComponent[3][i][0]<<"\t"<<NBadComponent[3][i][1]<<"\t"<<NBadComponent[3][i][2]<<"\t"<<NBadComponent[3][i][3];
   for (int i=10;i<19;++i)
     cout << "\nTEC- Disk " << i-9 << " :\t\t"<<NBadComponent[3][i][0]<<"\t"<<NBadComponent[3][i][1]<<"\t"<<NBadComponent[3][i][2]<<"\t"<<NBadComponent[3][i][3];
   cout << "\n";
  
   cout << "\n----------------------------------------------------------------\n\t\t   Detid  \tModules Fibers Apvs\n----------------------------------------------------------------";
   for (int i=1;i<5;++i)
     cout << "\nTIB Layer " << i << " :" << ssV[0][i].str();
   cout << "\n";
   for (int i=1;i<4;++i)
     cout << "\nTID+ Disk " << i << " :" << ssV[1][i].str();
   for (int i=4;i<7;++i)
     cout << "\nTID- Disk " << i-3 << " :" << ssV[1][i].str();
   cout << "\n";
   for (int i=1;i<7;++i)
     cout << "\nTOB Layer " << i << " :" << ssV[2][i].str();
   cout << "\n";
   for (int i=1;i<10;++i)
     cout << "\nTEC+ Disk " << i << " :" << ssV[3][i].str();
   for (int i=10;i<19;++i)
     cout << "\nTEC- Disk " << i-9 << " :" << ssV[3][i].str();
 
 }

void SiStripHitEffFromCalibTree::algoBeginJob ( )

privatevirtual

Definition at line 156 of file SiStripHitEffFromCalibTree.cc.

                                               {
   //I have no idea what goes here
   //fs->make<TTree>("HitEffHistos","Tree of the inefficient hit histograms");
 }

void SiStripHitEffFromCalibTree::algoEndJob ( )

privatevirtual

Reimplemented from ConditionDBWriter< SiStripBadStrip >.

Definition at line 161 of file SiStripHitEffFromCalibTree.cc.

                                             {
   //Still have no idea what goes here
 
 }

float SiStripHitEffFromCalibTree::calcPhi	(	float	x,
		float	y
	)

private

Definition at line 838 of file SiStripHitEffFromCalibTree.cc.

References phi, and Pi.

Referenced by makeHotColdMaps().

                                                           {
   float phi = 0;
   float Pi = 3.14159;
   if((x>=0)&&(y>=0)) phi = atan(y/x);
   else if((x>=0)&&(y<=0)) phi = atan(y/x) + 2*Pi;
   else if((x<=0)&&(y>=0)) phi = atan(y/x) + Pi;
   else phi = atan(y/x) + Pi;
   phi = phi*180.0/Pi;
 
   return phi;
 } 

SiStripBadStrip * SiStripHitEffFromCalibTree::getNewObject ( )

privatevirtual

Implements ConditionDBWriter< SiStripBadStrip >.

Definition at line 821 of file SiStripHitEffFromCalibTree.cc.

References SiStripBadStrip::getDataVectorBegin(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), getGTfromDQMFile::obj, SiStripBadStrip::put(), and quality_.

                                                           {
   //Need this for a Condition DB Writer
   //Initialize a return variable
   SiStripBadStrip* obj=new SiStripBadStrip();
   
   SiStripBadStrip::RegistryIterator rIter=quality_->getRegistryVectorBegin();
   SiStripBadStrip::RegistryIterator rIterEnd=quality_->getRegistryVectorEnd();
   
   for(;rIter!=rIterEnd;++rIter){
     SiStripBadStrip::Range range(quality_->getDataVectorBegin()+rIter->ibegin,quality_->getDataVectorBegin()+rIter->iend);
     if ( ! obj->put(rIter->detid,range) )
     edm::LogError("SiStripHitEffFromCalibTree")<<"[SiStripHitEffFromCalibTree::getNewObject] detid already exists"<<std::endl;
   }
   
   return obj;
 }

void SiStripHitEffFromCalibTree::makeHotColdMaps ( )

private

Definition at line 480 of file SiStripHitEffFromCalibTree.cc.

References begin, calcPhi(), gather_cfg::cout, fs, hits, HotColdMaps, TFileDirectory::make(), and phi.

Referenced by algoAnalyze().

                                                  {
   cout << "Entering hot cold map generation!\n";
   TStyle* gStyle = new TStyle("gStyle","myStyle");
   gStyle->cd();
   gStyle->SetPalette(1);
   gStyle->SetCanvasColor(kWhite);
   gStyle->SetOptStat(0);
   //Here we make the hot/cold color maps that we love so very much
   //Already have access to the data as a private variable
   //Create all of the histograms in the TFileService 
   TH2F *temph2;
   for(Long_t maplayer = 1; maplayer <=22; maplayer++) {                                           
     //Initialize all of the histograms                                                    
     if(maplayer > 0 && maplayer <= 4) {                                                   
     //We are in the TIB                                                       
     temph2 = fs->make<TH2F>(Form("%s%i","TIB",(int)(maplayer)),"TIB",100,-1,361,100,-100,100);                    
     temph2->GetXaxis()->SetTitle("Phi");
     temph2->GetXaxis()->SetBinLabel(1,TString("360"));
     temph2->GetXaxis()->SetBinLabel(50,TString("180"));
     temph2->GetXaxis()->SetBinLabel(100,TString("0"));                                                
     temph2->GetYaxis()->SetTitle("Global Z");
     temph2->SetOption("colz");                                            
     HotColdMaps.push_back(temph2);                                                    
       }                                                                   
       else if(maplayer > 4 && maplayer <= 10) {                                               
     //We are in the TOB                                                       
     temph2 = fs->make<TH2F>(Form("%s%i","TOB",(int)(maplayer-4)),"TOB",100,-1,361,100,-120,120);                  
     temph2->GetXaxis()->SetTitle("Phi");
     temph2->GetXaxis()->SetBinLabel(1,TString("360"));
     temph2->GetXaxis()->SetBinLabel(50,TString("180"));
     temph2->GetXaxis()->SetBinLabel(100,TString("0"));                                                
     temph2->GetYaxis()->SetTitle("Global Z");
     temph2->SetOption("colz");                                                
     HotColdMaps.push_back(temph2);                                                    
       }                                                                   
       else if(maplayer > 10 && maplayer <= 13) {                                              
     //We are in the TID                                                       
     //Split by +/-                                                            
     temph2 = fs->make<TH2F>(Form("%s%i","TID-",(int)(maplayer-10)),"TID-",100,-100,100,100,-100,100);             
     temph2->GetXaxis()->SetTitle("Global Y");
     temph2->GetXaxis()->SetBinLabel(1,TString("+Y"));
     temph2->GetXaxis()->SetBinLabel(50,TString("0"));
     temph2->GetXaxis()->SetBinLabel(100,TString("-Y"));                                               
     temph2->GetYaxis()->SetTitle("Global X");
     temph2->GetYaxis()->SetBinLabel(1,TString("-X"));
     temph2->GetYaxis()->SetBinLabel(50,TString("0"));
     temph2->GetYaxis()->SetBinLabel(100,TString("+X"));
     temph2->SetOption("colz");                                                
     HotColdMaps.push_back(temph2);                                                    
     temph2 = fs->make<TH2F>(Form("%s%i","TID+",(int)(maplayer-10)),"TID+",100,-100,100,100,-100,100);             
     temph2->GetXaxis()->SetTitle("Global Y");
     temph2->GetXaxis()->SetBinLabel(1,TString("+Y"));
     temph2->GetXaxis()->SetBinLabel(50,TString("0"));
     temph2->GetXaxis()->SetBinLabel(100,TString("-Y"));                                               
     temph2->GetYaxis()->SetTitle("Global X");
     temph2->GetYaxis()->SetBinLabel(1,TString("-X"));
     temph2->GetYaxis()->SetBinLabel(50,TString("0"));
     temph2->GetYaxis()->SetBinLabel(100,TString("+X"));                                               
     temph2->SetOption("colz");
     HotColdMaps.push_back(temph2);                                                    
       }                                                                   
       else if(maplayer > 13) {
         //We are in the TEC
         //Split by +/-
         temph2 = fs->make<TH2F>(Form("%s%i","TEC-",(int)(maplayer-13)),"TEC-",100,-120,120,100,-120,120);
         temph2->GetXaxis()->SetTitle("Global Y");
     temph2->GetXaxis()->SetBinLabel(1,TString("+Y"));
     temph2->GetXaxis()->SetBinLabel(50,TString("0"));
     temph2->GetXaxis()->SetBinLabel(100,TString("-Y"));                                               
     temph2->GetYaxis()->SetTitle("Global X");
     temph2->GetYaxis()->SetBinLabel(1,TString("-X"));
     temph2->GetYaxis()->SetBinLabel(50,TString("0"));
     temph2->GetYaxis()->SetBinLabel(100,TString("+X"));
     temph2->SetOption("colz");
         HotColdMaps.push_back(temph2);
         temph2 = fs->make<TH2F>(Form("%s%i","TEC+",(int)(maplayer-13)),"TEC+",100,-120,120,100,-120,120);
         temph2->GetXaxis()->SetTitle("Global Y");
     temph2->GetXaxis()->SetBinLabel(1,TString("+Y"));
     temph2->GetXaxis()->SetBinLabel(50,TString("0"));
     temph2->GetXaxis()->SetBinLabel(100,TString("-Y"));                                               
     temph2->GetYaxis()->SetTitle("Global X");
     temph2->GetYaxis()->SetBinLabel(1,TString("-X"));
     temph2->GetYaxis()->SetBinLabel(50,TString("0"));
     temph2->GetYaxis()->SetBinLabel(100,TString("+X"));
     temph2->SetOption("colz");
         HotColdMaps.push_back(temph2);
     }
   }
   for(Long_t mylayer = 1; mylayer <= 22; mylayer++) {
     //Determine what kind of plot we want to write out
     //Loop through the entirety of each layer
     //Create an array of the histograms
     vector<hit>::const_iterator iter;
     for(iter = hits[mylayer].begin(); iter != hits[mylayer].end(); iter++) {
       //Looping over the particular layer
       //Fill by 360-x to get the proper location to compare with TKMaps of phi
       //Also global xy is messed up
       if(mylayer > 0 && mylayer <= 4) {
         //We are in the TIB
         float phi = calcPhi(iter->x, iter->y);
         HotColdMaps[mylayer - 1]->Fill(360.-phi,iter->z,1.); 
       }
       else if(mylayer > 4 && mylayer <= 10) {
         //We are in the TOB
         float phi = calcPhi(iter->x,iter->y);
         HotColdMaps[mylayer - 1]->Fill(360.-phi,iter->z,1.);
       }
       else if(mylayer > 10 && mylayer <= 13) {
         //We are in the TID
         //There are 2 different maps here
         int side = (((iter->id)>>13) & 0x3);
     if(side == 1) HotColdMaps[(mylayer - 1) + (mylayer - 11)]->Fill(-iter->y,iter->x,1.); 
         else if(side == 2) HotColdMaps[(mylayer - 1) + (mylayer - 10)]->Fill(-iter->y,iter->x,1.);
         //if(side == 1) HotColdMaps[(mylayer - 1) + (mylayer - 11)]->Fill(iter->x,iter->y,1.); 
         //else if(side == 2) HotColdMaps[(mylayer - 1) + (mylayer - 10)]->Fill(iter->x,iter->y,1.);
       }
       else if(mylayer > 13) {
         //We are in the TEC
         //There are 2 different maps here
         int side = (((iter->id)>>18) & 0x3);
     if(side == 1) HotColdMaps[(mylayer + 2) + (mylayer - 14)]->Fill(-iter->y,iter->x,1.);
         else if(side == 2) HotColdMaps[(mylayer + 2) + (mylayer - 13)]->Fill(-iter->y,iter->x,1.);
         //if(side == 1) HotColdMaps[(mylayer + 2) + (mylayer - 14)]->Fill(iter->x,iter->y,1.);
         //else if(side == 2) HotColdMaps[(mylayer + 2) + (mylayer - 13)]->Fill(iter->x,iter->y,1.);
       }
     }
   }
   cout << "Finished HotCold Map Generation\n";
 }

void SiStripHitEffFromCalibTree::makeSQLite ( )

private

Definition at line 654 of file SiStripHitEffFromCalibTree.cc.

References BadModules, SiStripQuality::compact(), gather_cfg::cout, SiStripBadStrip::encode(), SiStripQuality::fillBadComponents(), SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), SiStripBadStrip::put(), quality_, reader, and AlCaHarvesting_cff::SiStripQuality.

Referenced by algoAnalyze().

                                             {
   //Generate the SQLite file for use in the Database of the bad modules!
   cout << "Entering SQLite file generation!\n";
   std::vector<unsigned int> BadStripList;
   unsigned short NStrips;
   unsigned int id1;
   SiStripQuality* pQuality = new SiStripQuality;
   //This is the list of the bad strips, use to mask out entire APVs
   //Now simply go through the bad hit list and mask out things that
   //are bad!
   map< unsigned int, double >::const_iterator it;
   for(it = BadModules.begin(); it != BadModules.end(); it++) {
     //We need to figure out how many strips are in this particular module
     //To Mask correctly!
     NStrips=reader->getNumberOfApvsAndStripLength((*it).first).first*128;
     cout << "Number of strips module " << (*it).first << " is " << NStrips << endl;
     BadStripList.push_back(pQuality->encode(0,NStrips,0));
     //Now compact into a single bad module
     id1=(unsigned int)(*it).first;
     cout << "ID1 shoudl match list of modules above " << id1 << endl;
     quality_->compact(id1,BadStripList);
     SiStripQuality::Range range(BadStripList.begin(),BadStripList.end());
     quality_->put(id1,range);
     BadStripList.clear();
   }
   //Fill all the bad components now
   quality_->fillBadComponents();
 }

void SiStripHitEffFromCalibTree::makeSummary ( )

private

Definition at line 697 of file SiStripHitEffFromCalibTree.cc.

References cond::ecalcond::all, alllayerfound, alllayertotal, gather_cfg::cout, newFWLiteAna::found, fs, goodlayerfound, goodlayertotal, i, j, gen::k, diffTwoXMLs::label, create_public_lumi_plots::leg, and TFileDirectory::make().

Referenced by algoAnalyze().

                                              {
   //setTDRStyle();
   
   int nLayers = 34;
   
   TH1F *found = fs->make<TH1F>("found","found",nLayers+1,0,nLayers+1);
   TH1F *all = fs->make<TH1F>("all","all",nLayers+1,0,nLayers+1);
   TH1F *found2 = fs->make<TH1F>("found2","found2",nLayers+1,0,nLayers+1);
   TH1F *all2 = fs->make<TH1F>("all2","all2",nLayers+1,0,nLayers+1);
   // first bin only to keep real data off the y axis so set to -1
   found->SetBinContent(0,-1);
   all->SetBinContent(0,1);
   
   TCanvas *c7 =new TCanvas("c7"," test ",10,10,800,600);
   c7->SetFillColor(0);
   c7->SetGrid();
 
   for (Long_t i=1; i< nLayers+1; ++i) {
     if (i==10) i++;
     if (i==25) i++;
     if (i==34) break;
 
     cout << "Fill only good modules layer " << i << ":  S = " << goodlayerfound[i] << "    B = " << goodlayertotal[i] << endl;
     if (goodlayertotal[i] > 5) {
       found->SetBinContent(i,goodlayerfound[i]);
       all->SetBinContent(i,goodlayertotal[i]);
     } else {
       found->SetBinContent(i,0);
       all->SetBinContent(i,10);
     }
     
     cout << "Filling all modules layer " << i << ":  S = " << alllayerfound[i] << "    B = " << alllayertotal[i] << endl;
     if (alllayertotal[i] > 5) {
       found2->SetBinContent(i,alllayerfound[i]);
       all2->SetBinContent(i,alllayertotal[i]);
     } else {
       found2->SetBinContent(i,0);
       all2->SetBinContent(i,10);
     }
 
   }
 
   found->Sumw2();
   all->Sumw2();
 
   found2->Sumw2();
   all2->Sumw2();
 
   TGraphAsymmErrors *gr = new TGraphAsymmErrors(nLayers+1);
   gr->BayesDivide(found,all); 
 
   TGraphAsymmErrors *gr2 = new TGraphAsymmErrors(nLayers+1);
   gr2->BayesDivide(found2,all2);
 
   for(int j = 0; j<nLayers+1; j++){
     gr->SetPointError(j, 0., 0., gr->GetErrorYlow(j),gr->GetErrorYhigh(j) );
     gr2->SetPointError(j, 0., 0., gr2->GetErrorYlow(j),gr2->GetErrorYhigh(j) );
   }
 
   gr->GetXaxis()->SetLimits(0,nLayers);
   gr->SetMarkerColor(2);
   gr->SetMarkerSize(1.2);
   gr->SetLineColor(2);
   gr->SetLineWidth(4);
   gr->SetMarkerStyle(20);
   gr->SetMinimum(0.90);
   gr->SetMaximum(1.001);
   gr->GetYaxis()->SetTitle("Efficiency");
 
   gr2->GetXaxis()->SetLimits(0,nLayers);
   gr2->SetMarkerColor(1);
   gr2->SetMarkerSize(1.2);
   gr2->SetLineColor(1);
   gr2->SetLineWidth(4);
   gr2->SetMarkerStyle(21);
   gr2->SetMinimum(0.90);
   gr2->SetMaximum(1.001);
   gr2->GetYaxis()->SetTitle("Efficiency");
   //cout << "starting labels" << endl;
   //for ( int k=1; k<nLayers+1; k++) {
   for ( Long_t k=1; k<nLayers+1; k++) {
     if (k==10) k++;
     if (k==25) k++;
     if (k==34) break;
     TString label;
     if (k<5) {
       label = TString("TIB  ") + k;
     } else if (k>4&&k<11) {
       label = TString("TOB  ")+(k-4);
     } else if (k>10&&k<14) {
       label = TString("TID- ")+(k-10);
     } else if (k>13&&k<17) {
       label = TString("TID+ ")+(k-13);
     } else if (k>16&&k<26) {
       label = TString("TEC- ")+(k-16);
     } else if (k>25) {
       label = TString("TEC+ ")+(k-25);
     }
     gr->GetXaxis()->SetBinLabel(((k+1)*100)/(nLayers)-2,label);
     gr2->GetXaxis()->SetBinLabel(((k+1)*100)/(nLayers)-2,label);
   }
   
   gr->Draw("AP");
   gr->GetXaxis()->SetNdivisions(36);
 
   c7->cd();
   TPad *overlay = new TPad("overlay","",0,0,1,1);
   overlay->SetFillStyle(4000);
   overlay->SetFillColor(0);
   overlay->SetFrameFillStyle(4000);
   overlay->Draw("same");
   overlay->cd();
   gr2->Draw("AP");
 
   TLegend *leg = new TLegend(0.70,0.20,0.92,0.39);
   leg->AddEntry(gr,"Good Modules","p");
   leg->AddEntry(gr2,"All Modules","p");
   leg->SetTextSize(0.020);
   leg->SetFillColor(0);
   leg->Draw("same");
   
   c7->SaveAs("Summary.png");
 }

void SiStripHitEffFromCalibTree::makeTKMap ( )

private

Definition at line 610 of file SiStripHitEffFromCalibTree.cc.

References BadModules, begin, gather_cfg::cout, TrackerMap::fill(), TrackerMap::fillc(), first, i, layerfound, layertotal, modCounter, nModsMin, TrackerMap::save(), edm::second(), threshold, tkmap, and tkmapbad.

Referenced by algoAnalyze().

                                            {
   cout << "Entering TKMap generation!\n";
   tkmap = new TrackerMap("  Detector Inefficiency  ");
   tkmapbad = new TrackerMap("  Inefficient Modules  ");
   for(Long_t i = 1; i <= 22; i++) {
     layertotal[i] = 0;
     layerfound[i] = 0;
     //Loop over every layer, extracting the information from
     //the map of the efficiencies
     map<unsigned int, pair<unsigned int, unsigned int> >::const_iterator ih;
     for( ih = modCounter[i].begin(); ih != modCounter[i].end(); ih++) {
       //We should be in the layer in question, and looping over all of the modules in said layer
       //Generate the list for the TKmap, and the bad module list
       double myeff = (double)(((*ih).second).second)/(((*ih).second).first);
       if ( ((((*ih).second).first) >= nModsMin) && (myeff < threshold) ) {
         //We have a bad module, put it in the list!
     BadModules[(*ih).first] = myeff;
     tkmapbad->fillc((*ih).first,255,0,0);
     cout << "Layer " << i << " module " << (*ih).first << " efficiency " << myeff << " " << (((*ih).second).second) << "/" << (((*ih).second).first) << endl;
       }
       else {
         //Fill the bad list with empty results for every module
         tkmapbad->fillc((*ih).first,255,255,255);
       }
       if((((*ih).second).first) < 100 ) {
         cout << "Module " << (*ih).first << " layer " << i << " is under occupancy at " << (((*ih).second).first) << endl;
       }
       //Put any module into the TKMap
       //Should call module ID, and then 1- efficiency for that module
       //if((*ih).first == 369137820) {
       //  cout << "Module 369137820 has 1-eff of " << 1.-myeff << endl;
     //cout << "Which is " << ((*ih).second).second << "/" << ((*ih).second).first << endl;
       //}
       tkmap->fill((*ih).first,1.-myeff);
       //Find the total number of hits in the module
       layertotal[i] += int(((*ih).second).first);
       layerfound[i] += int(((*ih).second).second);
     }
   }
   tkmap->save(true, 0, 0, "SiStripHitEffTKMap.png");
   tkmapbad->save(true, 0, 0, "SiStripHitEffTKMapBad.png");
   cout << "Finished TKMap Generation\n";
 }

void SiStripHitEffFromCalibTree::SetBadComponents	(	int	i,
		int	component,
		SiStripQuality::BadComponent &	BC,
		std::stringstream	ssV[4][19],
		int	NBadComponent[4][19][4]
	)

private

Definition at line 850 of file SiStripHitEffFromCalibTree.cc.

References SiStripQuality::BadComponent::BadApvs, SiStripQuality::BadComponent::BadFibers, SiStripQuality::BadComponent::BadModule, SiStripQuality::BadComponent::detid, SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), i, and reader.

Referenced by algoAnalyze().

                                                                                                                                                              {
  
   int napv=reader->getNumberOfApvsAndStripLength(BC.detid).first;
  
   ssV[i][component] << "\n\t\t " 
             << BC.detid 
             << " \t " << BC.BadModule << " \t " 
             << ( (BC.BadFibers)&0x1 ) << " ";
   if (napv==4)
     ssV[i][component] << "x " <<( (BC.BadFibers>>1)&0x1 );
    
   if (napv==6)
     ssV[i][component] << ( (BC.BadFibers>>1)&0x1 ) << " "
               << ( (BC.BadFibers>>2)&0x1 );
     ssV[i][component] << " \t " 
               << ( (BC.BadApvs)&0x1 ) << " " 
               << ( (BC.BadApvs>>1)&0x1 ) << " ";
   if (napv==4) 
     ssV[i][component] << "x x " << ( (BC.BadApvs>>2)&0x1 ) << " " 
               << ( (BC.BadApvs>>3)&0x1 );
   if (napv==6) 
     ssV[i][component] << ( (BC.BadApvs>>2)&0x1 ) << " " 
               << ( (BC.BadApvs>>3)&0x1 ) << " " 
               << ( (BC.BadApvs>>4)&0x1 ) << " " 
               << ( (BC.BadApvs>>5)&0x1 ) << " "; 
  
   if (BC.BadApvs){
     NBadComponent[i][0][2]+= ( (BC.BadApvs>>5)&0x1 )+ ( (BC.BadApvs>>4)&0x1 ) + ( (BC.BadApvs>>3)&0x1 ) + 
       ( (BC.BadApvs>>2)&0x1 )+ ( (BC.BadApvs>>1)&0x1 ) + ( (BC.BadApvs)&0x1 );
     NBadComponent[i][component][2]+= ( (BC.BadApvs>>5)&0x1 )+ ( (BC.BadApvs>>4)&0x1 ) + ( (BC.BadApvs>>3)&0x1 ) + 
       ( (BC.BadApvs>>2)&0x1 )+ ( (BC.BadApvs>>1)&0x1 ) + ( (BC.BadApvs)&0x1 );
   }
   if (BC.BadFibers){ 
     NBadComponent[i][0][1]+= ( (BC.BadFibers>>2)&0x1 )+ ( (BC.BadFibers>>1)&0x1 ) + ( (BC.BadFibers)&0x1 );
     NBadComponent[i][component][1]+= ( (BC.BadFibers>>2)&0x1 )+ ( (BC.BadFibers>>1)&0x1 ) + ( (BC.BadFibers)&0x1 );
   }   
   if (BC.BadModule){
     NBadComponent[i][0][0]++;
     NBadComponent[i][component][0]++;
   }
 }

void SiStripHitEffFromCalibTree::totalStatistics ( )

private

Definition at line 683 of file SiStripHitEffFromCalibTree.cc.

References gather_cfg::cout, i, layerfound, and layertotal.

Referenced by algoAnalyze().

                                                  {
   //Calculate the statistics by layer
   int totalfound = 0;
   int totaltotal = 0;
   double layereff;
   for(Long_t i=1; i<=22; i++) {
     layereff = double(layerfound[i])/double(layertotal[i]);
     cout << "Layer " << i << " has total efficiency " << layereff << " " << layerfound[i] << "/" << layertotal[i] << endl;
     totalfound += layerfound[i];
     totaltotal += layertotal[i];
   }
   cout << "The total efficiency is " << double(totalfound)/double(totaltotal) << endl;
 }