|
|
|
Definition at line 115 of file SiStripGainFromCalibTree.cc.
| SiStripGainFromCalibTree::SiStripGainFromCalibTree | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 190 of file SiStripGainFromCalibTree.cc.
References AllowSaturation, CalibrationLevel, FirstSetOfConstants, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_calibrationPath, MaxChi2OverNDF, MaxMPVError, MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackMomentum, MinNrEntries, MinTrackEta, MinTrackHits, MinTrackMomentum, OldGainRemoving, OutputGains, useCalibration, Validation, and VInputFiles.
: ConditionDBWriter<SiStripApvGain>(iConfig) { OutputGains = iConfig.getParameter<std::string>("OutputGains"); MinNrEntries = iConfig.getUntrackedParameter<double> ("minNrEntries" , 20); MaxMPVError = iConfig.getUntrackedParameter<double> ("maxMPVError" , 500.0); MaxChi2OverNDF = iConfig.getUntrackedParameter<double> ("maxChi2OverNDF" , 5.0); MinTrackMomentum = iConfig.getUntrackedParameter<double> ("minTrackMomentum" , 3.0); MaxTrackMomentum = iConfig.getUntrackedParameter<double> ("maxTrackMomentum" , 99999.0); MinTrackEta = iConfig.getUntrackedParameter<double> ("minTrackEta" , -5.0); MaxTrackEta = iConfig.getUntrackedParameter<double> ("maxTrackEta" , 5.0); MaxNrStrips = iConfig.getUntrackedParameter<unsigned>("maxNrStrips" , 2); MinTrackHits = iConfig.getUntrackedParameter<unsigned>("MinTrackHits" , 8); MaxTrackChiOverNdf = iConfig.getUntrackedParameter<double> ("MaxTrackChiOverNdf" , 3); AllowSaturation = iConfig.getUntrackedParameter<bool> ("AllowSaturation" , false); FirstSetOfConstants = iConfig.getUntrackedParameter<bool> ("FirstSetOfConstants", true); Validation = iConfig.getUntrackedParameter<bool> ("Validation" , false); OldGainRemoving = iConfig.getUntrackedParameter<bool> ("OldGainRemoving" , false); CalibrationLevel = iConfig.getUntrackedParameter<int> ("CalibrationLevel" , 0); VInputFiles = iConfig.getParameter<vector<string> > ("InputFiles"); useCalibration = iConfig.getUntrackedParameter<bool>("UseCalibration", false); m_calibrationPath = iConfig.getUntrackedParameter<string>("calibrationPath"); }
| SiStripGainFromCalibTree::~SiStripGainFromCalibTree | ( | ) |
Definition at line 699 of file SiStripGainFromCalibTree.cc.
{
}
| void SiStripGainFromCalibTree::algoAnalyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Reimplemented from ConditionDBWriter< SiStripApvGain >.
Definition at line 321 of file SiStripGainFromCalibTree.cc.
{
}
| void SiStripGainFromCalibTree::algoAnalyzeTheTree | ( | ) | [private] |
Definition at line 363 of file SiStripGainFromCalibTree.cc.
References AllowSaturation, sistrip::APV, APVsColl, stAPVGain::CalibGain, DeDxDiscriminatorTools::charge(), Charge_Vs_Index, Charge_Vs_Index_Absolute, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, ERun, stAPVGain::Eta, funct::false, FirstSetOfConstants, i, stAPVGain::Index, MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackMomentum, MinTrackEta, MinTrackHits, MinTrackMomentum, NCluster, NEvent, NTrack, NULL, OldGainRemoving, scaleCards::path, stAPVGain::PreviousGain, cond::runnumber, alignCSCRings::s, SRun, stAPVGain::SubDet, StripSubdetector::TEC, stAPVGain::Thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, diffTreeTool::tree, funct::true, useCalibration, Validation, and VInputFiles.
Referenced by algoBeginJob().
{
for(unsigned int i=0;i<VInputFiles.size();i++){
printf("Openning file %3i/%3i --> %s\n",i+1, (int)VInputFiles.size(), (char*)(VInputFiles[i].c_str())); fflush(stdout);
TChain* tree = new TChain("gainCalibrationTree/tree");
tree->Add(VInputFiles[i].c_str());
TString EventPrefix("");
TString EventSuffix("");
TString TrackPrefix("track");
TString TrackSuffix("");
TString CalibPrefix("GainCalibration");
TString CalibSuffix("");
unsigned int eventnumber = 0; tree->SetBranchAddress(EventPrefix + "event" + EventSuffix, &eventnumber , NULL);
unsigned int runnumber = 0; tree->SetBranchAddress(EventPrefix + "run" + EventSuffix, &runnumber , NULL);
std::vector<bool>* TrigTech = 0; tree->SetBranchAddress(EventPrefix + "TrigTech" + EventSuffix, &TrigTech , NULL);
std::vector<double>* trackchi2ndof = 0; tree->SetBranchAddress(TrackPrefix + "chi2ndof" + TrackSuffix, &trackchi2ndof , NULL);
std::vector<float>* trackp = 0; tree->SetBranchAddress(TrackPrefix + "momentum" + TrackSuffix, &trackp , NULL);
std::vector<float>* trackpt = 0; tree->SetBranchAddress(TrackPrefix + "pt" + TrackSuffix, &trackpt , NULL);
std::vector<double>* tracketa = 0; tree->SetBranchAddress(TrackPrefix + "eta" + TrackSuffix, &tracketa , NULL);
std::vector<double>* trackphi = 0; tree->SetBranchAddress(TrackPrefix + "phi" + TrackSuffix, &trackphi , NULL);
std::vector<unsigned int>* trackhitsvalid = 0; tree->SetBranchAddress(TrackPrefix + "hitsvalid" + TrackSuffix, &trackhitsvalid, NULL);
std::vector<int>* trackindex = 0; tree->SetBranchAddress(CalibPrefix + "trackindex" + CalibSuffix, &trackindex , NULL);
std::vector<unsigned int>* rawid = 0; tree->SetBranchAddress(CalibPrefix + "rawid" + CalibSuffix, &rawid , NULL);
std::vector<float>* localdirx = 0; tree->SetBranchAddress(CalibPrefix + "localdirx" + CalibSuffix, &localdirx , NULL);
std::vector<float>* localdiry = 0; tree->SetBranchAddress(CalibPrefix + "localdiry" + CalibSuffix, &localdiry , NULL);
std::vector<float>* localdirz = 0; tree->SetBranchAddress(CalibPrefix + "localdirz" + CalibSuffix, &localdirz , NULL);
std::vector<unsigned short>* firststrip = 0; tree->SetBranchAddress(CalibPrefix + "firststrip" + CalibSuffix, &firststrip , NULL);
std::vector<unsigned short>* nstrips = 0; tree->SetBranchAddress(CalibPrefix + "nstrips" + CalibSuffix, &nstrips , NULL);
std::vector<bool>* saturation = 0; tree->SetBranchAddress(CalibPrefix + "saturation" + CalibSuffix, &saturation , NULL);
std::vector<bool>* overlapping = 0; tree->SetBranchAddress(CalibPrefix + "overlapping" + CalibSuffix, &overlapping , NULL);
std::vector<bool>* farfromedge = 0; tree->SetBranchAddress(CalibPrefix + "farfromedge" + CalibSuffix, &farfromedge , NULL);
std::vector<unsigned int>* charge = 0; tree->SetBranchAddress(CalibPrefix + "charge" + CalibSuffix, &charge , NULL);
std::vector<float>* path = 0; tree->SetBranchAddress(CalibPrefix + "path" + CalibSuffix, &path , NULL);
std::vector<float>* chargeoverpath = 0; tree->SetBranchAddress(CalibPrefix + "chargeoverpath" + CalibSuffix, &chargeoverpath, NULL);
std::vector<unsigned char>* amplitude = 0; tree->SetBranchAddress(CalibPrefix + "amplitude" + CalibSuffix, &litude , NULL);
std::vector<double>* gainused = 0; tree->SetBranchAddress(CalibPrefix + "gainused" + CalibSuffix, &gainused , NULL);
printf("Number of Events = %i + %i = %i\n",NEvent,(unsigned int)tree->GetEntries(),(unsigned int)(NEvent+tree->GetEntries()));
printf("Progressing Bar :0%% 20%% 40%% 60%% 80%% 100%%\n");
printf("Looping on the Tree :");
int TreeStep = tree->GetEntries()/50;if(TreeStep<=1)TreeStep=1;
for (unsigned int ientry = 0; ientry < tree->GetEntries(); ientry++) {
// for (unsigned int ientry = 0; ientry < tree->GetEntries() && ientry<50000; ientry++) {
if(ientry%TreeStep==0){printf(".");fflush(stdout);}
tree->GetEntry(ientry);
if(runnumber<SRun)SRun=runnumber;
if(runnumber>ERun)ERun=runnumber;
NEvent++;
NTrack+=(*trackp).size();
unsigned int FirstAmplitude=0;
for(unsigned int i=0;i<(*chargeoverpath).size();i++){
FirstAmplitude+=(*nstrips)[i];
int TI = (*trackindex)[i];
if((*tracketa )[TI] < MinTrackEta )continue;
if((*tracketa )[TI] > MaxTrackEta )continue;
if((*trackp )[TI] < MinTrackMomentum )continue;
if((*trackp )[TI] > MaxTrackMomentum )continue;
if((*trackhitsvalid)[TI] < MinTrackHits )continue;
if((*trackchi2ndof )[TI] > MaxTrackChiOverNdf )continue;
if((*farfromedge)[i] == false )continue;
if((*overlapping)[i] == true )continue;
if((*saturation )[i] && !AllowSaturation)continue;
if((*nstrips )[i] > MaxNrStrips )continue;
NCluster++;
stAPVGain* APV = APVsColl[((*rawid)[i]<<3) | ((*firststrip)[i]/128)];
// printf("detId=%7i run=%7i event=%9i charge=%5i cs=%3i\n",(*rawid)[i],runnumber,eventnumber,(*charge)[i],(*nstrips)[i]);
//double trans = atan2((*localdiry)[i],(*localdirx)[i])*(180/3.14159265);
//double alpha = acos ((*localdirx)[i] / sqrt( pow((*localdirx)[i],2) + pow((*localdirz)[i],2) ) ) * (180/3.14159265);
//double beta = acos ((*localdiry)[i] / sqrt( pow((*localdirx)[i],2) + pow((*localdirz)[i],2) ) ) * (180/3.14159265);
//printf("NStrip = %i : Charge = %i --> Path = %f --> ChargeOverPath=%f\n",(*nstrips)[i],(*charge)[i],(*path)[i],(*chargeoverpath)[i]);
//printf("Amplitudes: ");
//for(unsigned int a=0;a<(*nstrips)[i];a++){printf("%i ",(*amplitude)[FirstAmplitude+a]);}
//printf("\n");
int Charge = 0;
if(useCalibration || !FirstSetOfConstants){
bool Saturation = false;
for(unsigned int s=0;s<(*nstrips)[i];s++){
int StripCharge = (*amplitude)[FirstAmplitude-(*nstrips)[i]+s];
if(useCalibration && !FirstSetOfConstants){ StripCharge=(int)(StripCharge*(APV->PreviousGain/APV->CalibGain));
}else if(useCalibration){ StripCharge=(int)(StripCharge/APV->CalibGain);
}else if(!FirstSetOfConstants){ StripCharge=(int)(StripCharge*APV->PreviousGain);}
if(StripCharge>1024){
StripCharge = 255;
Saturation = true;
}else if(StripCharge>254){
StripCharge = 254;
Saturation = true;
}
Charge += StripCharge;
}
if(Saturation && !AllowSaturation)continue;
}else{
Charge = (*charge)[i];
}
//printf("ChargeDifference = %i Vs %i with Gain = %f\n",(*charge)[i],Charge,APV->CalibGain);
double ClusterChargeOverPath = ( (double) Charge )/(*path)[i] ;
if(Validation) {ClusterChargeOverPath/=(*gainused)[i];}
if(OldGainRemoving){ClusterChargeOverPath*=(*gainused)[i];}
Charge_Vs_Index_Absolute->Fill(APV->Index,Charge);
Charge_Vs_Index ->Fill(APV->Index,ClusterChargeOverPath);
if(APV->SubDet==StripSubdetector::TIB){ Charge_Vs_PathlengthTIB ->Fill((*path)[i],Charge);
}else if(APV->SubDet==StripSubdetector::TOB){ Charge_Vs_PathlengthTOB ->Fill((*path)[i],Charge);
}else if(APV->SubDet==StripSubdetector::TID){
if(APV->Eta<0){ Charge_Vs_PathlengthTIDM ->Fill((*path)[i],Charge);
}else if(APV->Eta>0){ Charge_Vs_PathlengthTIDP ->Fill((*path)[i],Charge);
}
}else if(APV->SubDet==StripSubdetector::TEC){
if(APV->Eta<0){
if(APV->Thickness<0.04){ Charge_Vs_PathlengthTECM1->Fill((*path)[i],Charge);
}else if(APV->Thickness>0.04){ Charge_Vs_PathlengthTECM2->Fill((*path)[i],Charge);
}
}else if(APV->Eta>0){
if(APV->Thickness<0.04){ Charge_Vs_PathlengthTECP1->Fill((*path)[i],Charge);
}else if(APV->Thickness>0.04){ Charge_Vs_PathlengthTECP2->Fill((*path)[i],Charge);
}
}
}
}// END OF ON-CLUSTER LOOP
}printf("\n");// END OF EVENT LOOP
}
}
| void SiStripGainFromCalibTree::algoBeginJob | ( | const edm::EventSetup & | iSetup | ) | [private, virtual] |
Reimplemented from ConditionDBWriter< SiStripApvGain >.
Definition at line 223 of file SiStripGainFromCalibTree.cc.
References algoAnalyzeTheTree(), algoComputeMPVandGain(), sistrip::APV, stAPVGain::APVId, APVsColl, APVsCollOrdered, BAD, PV3DBase< T, PVType, FrameType >::basicVector(), stAPVGain::Bin, BoundSurface::bounds(), Charge_Vs_Index, Charge_Vs_Index_Absolute, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, stAPVGain::DetId, ERun, stAPVGain::Eta, cmsRelvalreport::exit, FirstSetOfConstants, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, edm::EventSetup::get(), GOOD, i, stAPVGain::Index, stAPVGain::isMasked, edm::ESHandleBase::isValid(), j, MakeCalibrationMap(), NCluster, stAPVGain::NEntries, NEvent, StripTopology::nstrips(), NTrack, stAPVGain::Phi, GeomDet::position(), stAPVGain::PreviousGain, stAPVGain::R, DetId::rawId(), StripGeomDetUnit::specificTopology(), SRun, stAPVGain::SubDet, DetId::subdetId(), GeomDet::surface(), StripSubdetector::TEC, tfs, stAPVGain::Thickness, Bounds::thickness(), StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, stAPVGain::x, stAPVGain::y, and stAPVGain::z.
{
Charge_Vs_Index = tfs->make<TH2F>("Charge_Vs_Index" , "Charge_Vs_Index" , 72785, 0 , 72784,1000,0,2000);
Charge_Vs_Index_Absolute = tfs->make<TH2F>("Charge_Vs_Index_Absolute" , "Charge_Vs_Index_Absolute" , 72785, 0 , 72784, 500,0,2000);
Charge_Vs_PathlengthTIB = tfs->make<TH2F>("Charge_Vs_PathlengthTIB" , "Charge_Vs_PathlengthTIB" , 20 , 0.3 , 1.3 , 250,0,2000);
Charge_Vs_PathlengthTOB = tfs->make<TH2F>("Charge_Vs_PathlengthTOB" , "Charge_Vs_PathlengthTOB" , 20 , 0.3 , 1.3 , 250,0,2000);
Charge_Vs_PathlengthTIDP = tfs->make<TH2F>("Charge_Vs_PathlengthTIDP" , "Charge_Vs_PathlengthTIDP" , 20 , 0.3 , 1.3 , 250,0,2000);
Charge_Vs_PathlengthTIDM = tfs->make<TH2F>("Charge_Vs_PathlengthTIDM" , "Charge_Vs_PathlengthTIDM" , 20 , 0.3 , 1.3 , 250,0,2000);
Charge_Vs_PathlengthTECP1 = tfs->make<TH2F>("Charge_Vs_PathlengthTECP1", "Charge_Vs_PathlengthTECP1", 20 , 0.3 , 1.3 , 250,0,2000);
Charge_Vs_PathlengthTECP2 = tfs->make<TH2F>("Charge_Vs_PathlengthTECP2", "Charge_Vs_PathlengthTECP2", 20 , 0.3 , 1.3 , 250,0,2000);
Charge_Vs_PathlengthTECM1 = tfs->make<TH2F>("Charge_Vs_PathlengthTECM1", "Charge_Vs_PathlengthTECM1", 20 , 0.3 , 1.3 , 250,0,2000);
Charge_Vs_PathlengthTECM2 = tfs->make<TH2F>("Charge_Vs_PathlengthTECM2", "Charge_Vs_PathlengthTECM2", 20 , 0.3 , 1.3 , 250,0,2000);
edm::ESHandle<TrackerGeometry> tkGeom;
iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
vector<GeomDet*> Det = tkGeom->dets();
edm::ESHandle<SiStripGain> gainHandle;
iSetup.get<SiStripGainRcd>().get(gainHandle);
if(!gainHandle.isValid()){printf("\n#####################\n\nERROR --> gainHandle is not valid\n\n#####################\n\n");exit(0);}
// size_t numberOfTag = gainHandle->getNumberOfTags();
for(unsigned int i=0;i<gainHandle->getNumberOfTags();i++){
printf("Reccord %i --> Rcd Name = %s Label Name = %s\n",i,gainHandle->getRcdName(i).c_str(), gainHandle->getLabelName(i).c_str());
}
edm::ESHandle<SiStripQuality> SiStripQuality_;
iSetup.get<SiStripQualityRcd>().get(SiStripQuality_);
unsigned int Index=0;
for(unsigned int i=0;i<Det.size();i++){
DetId Detid = Det[i]->geographicalId();
int SubDet = Detid.subdetId();
if( SubDet == StripSubdetector::TIB || SubDet == StripSubdetector::TID ||
SubDet == StripSubdetector::TOB || SubDet == StripSubdetector::TEC ){
StripGeomDetUnit* DetUnit = dynamic_cast<StripGeomDetUnit*> (Det[i]);
if(!DetUnit)continue;
const StripTopology& Topo = DetUnit->specificTopology();
unsigned int NAPV = Topo.nstrips()/128;
for(unsigned int j=0;j<NAPV;j++){
stAPVGain* APV = new stAPVGain;
APV->Index = Index;
APV->Bin = Charge_Vs_Index->GetXaxis()->FindBin(APV->Index);
APV->DetId = Detid.rawId();
APV->APVId = j;
APV->SubDet = SubDet;
APV->FitMPV = -1;
APV->FitMPVErr = -1;
APV->FitWidth = -1;
APV->FitWidthErr = -1;
APV->FitChi2 = -1;
APV->Gain = -1;
APV->PreviousGain = 1;
APV->x = DetUnit->position().basicVector().x();
APV->y = DetUnit->position().basicVector().y();
APV->z = DetUnit->position().basicVector().z();
APV->Eta = DetUnit->position().basicVector().eta();
APV->Phi = DetUnit->position().basicVector().phi();
APV->R = DetUnit->position().basicVector().transverse();
APV->Thickness = DetUnit->surface().bounds().thickness();
APV->NEntries = 0;
APV->isMasked = SiStripQuality_->IsApvBad(Detid.rawId(),j);
if(!FirstSetOfConstants){
if(gainHandle->getNumberOfTags()!=2){printf("ERROR: NUMBER OF GAIN TAG IS EXPECTED TO BE 2\n");fflush(stdout);exit(0);};
APV->PreviousGain = gainHandle->getApvGain(APV->APVId,gainHandle->getRange(APV->DetId, 1),1);
printf("DETID = %7i APVID=%1i Previous Gain=%8.4f\n",APV->DetId,APV->APVId,APV->PreviousGain);
}
APVsCollOrdered.push_back(APV);
APVsColl[(APV->DetId<<3) | APV->APVId] = APV;
Index++;
}
}
}
MakeCalibrationMap();
NEvent = 0;
NTrack = 0;
NCluster = 0;
SRun = 1<<31;
ERun = 0;
GOOD = 0;
BAD = 0;
algoAnalyzeTheTree();
algoComputeMPVandGain();
}
| void SiStripGainFromCalibTree::algoComputeMPVandGain | ( | ) | [private] |
Definition at line 511 of file SiStripGainFromCalibTree.cc.
References sistrip::APV, stAPVGain::APVId, APVsColl, BAD, stAPVGain::Bin, CalibrationLevel, Charge_Vs_Index, stAPVGain::DetId, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, getPeakOfLandau(), GOOD, Exhume::I, i, stAPVGain::NEntries, NULL, and storeOnTree().
Referenced by algoBeginJob().
{
unsigned int I=0;
TH1D* Proj = NULL;
double FitResults[5];
double MPVmean = 300;
printf("Progressing Bar :0%% 20%% 40%% 60%% 80%% 100%%\n");
printf("Fitting Charge Distribution :");
int TreeStep = APVsColl.size()/50;
for(__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual >::iterator it = APVsColl.begin();it!=APVsColl.end();it++,I++){
if(I%TreeStep==0){printf(".");fflush(stdout);}
//if(I>1000)break;
stAPVGain* APV = it->second;
Proj = (TH1D*)(Charge_Vs_Index->ProjectionY("",APV->Bin,APV->Bin,"e"));
if(!Proj)continue;
if(CalibrationLevel==0){
}else if(CalibrationLevel==1){
int SecondAPVId = APV->APVId;
if(SecondAPVId%2==0){ SecondAPVId = SecondAPVId+1; }else{ SecondAPVId = SecondAPVId-1; }
stAPVGain* APV2 = APVsColl[(APV->DetId<<3) | SecondAPVId];
TH1D* Proj2 = (TH1D*)(Charge_Vs_Index->ProjectionY("",APV2->Bin,APV2->Bin,"e"));
if(Proj2){Proj->Add(Proj2,1);delete Proj2;}
}else if(CalibrationLevel==2){
for(unsigned int i=0;i<6;i++){
__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual >::iterator tmpit;
tmpit = APVsColl.find((APV->DetId<<3) | i);
if(tmpit==APVsColl.end())continue;
stAPVGain* APV2 = tmpit->second;
if(APV2->DetId != APV->DetId || APV2->APVId == APV->APVId)continue;
TH1D* Proj2 = (TH1D*)(Charge_Vs_Index->ProjectionY("",APV2->Bin,APV2->Bin,"e"));
// if(Proj2 && APV->DetId==369171124)printf("B) DetId %6i APVId %1i --> NEntries = %f\n",APV2->DetId, APV2->APVId, Proj2->GetEntries());
if(Proj2){Proj->Add(Proj2,1);delete Proj2;}
}
}else{
CalibrationLevel = 0;
printf("Unknown Calibration Level, will assume %i\n",CalibrationLevel);
}
getPeakOfLandau(Proj,FitResults);
APV->FitMPV = FitResults[0];
APV->FitMPVErr = FitResults[1];
APV->FitWidth = FitResults[2];
APV->FitWidthErr = FitResults[3];
APV->FitChi2 = FitResults[4];
APV->NEntries = Proj->GetEntries();
if(APV->FitMPV>0){
APV->Gain = APV->FitMPV / MPVmean;
GOOD++;
}else{
APV->Gain = 1;
BAD++;
}
if(APV->Gain<=0) APV->Gain = 1;
// if(!FirstSetOfConstants) APV->Gain *= APV->PreviousGain;
//printf("%5i/%5i: %6i - %1i %5E Entries --> MPV = %f +- %f\n",I,APVsColl.size(),APV->DetId, APV->APVId, Proj->GetEntries(), FitResults[0], FitResults[1]);fflush(stdout);
delete Proj;
}printf("\n");
storeOnTree();
}
| void SiStripGainFromCalibTree::algoEndJob | ( | ) | [private, virtual] |
Reimplemented from ConditionDBWriter< SiStripApvGain >.
Definition at line 326 of file SiStripGainFromCalibTree.cc.
{
}
| SiStripApvGain * SiStripGainFromCalibTree::getNewObject | ( | ) | [private, virtual] |
Implements ConditionDBWriter< SiStripApvGain >.
Definition at line 671 of file SiStripGainFromCalibTree.cc.
References a, sistrip::APV, APVsCollOrdered, stAPVGain::DetId, stAPVGain::Gain, NULL, VarParsing::obj, and SiStripApvGain::put().
{
SiStripApvGain* obj = new SiStripApvGain();
std::vector<float>* theSiStripVector = NULL;
unsigned int PreviousDetId = 0;
for(unsigned int a=0;a<APVsCollOrdered.size();a++)
{
stAPVGain* APV = APVsCollOrdered[a];
if(APV==NULL){ printf("Bug\n"); continue; }
if(APV->DetId != PreviousDetId){
if(theSiStripVector!=NULL){
SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
if ( !obj->put(PreviousDetId,range) ) printf("Bug to put detId = %i\n",PreviousDetId);
}
theSiStripVector = new std::vector<float>;
PreviousDetId = APV->DetId;
}
theSiStripVector->push_back(APV->Gain);
}
if(theSiStripVector!=NULL){
SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
if ( !obj->put(PreviousDetId,range) ) printf("Bug to put detId = %i\n",PreviousDetId);
}
return obj;
}
| void SiStripGainFromCalibTree::getPeakOfLandau | ( | TH1 * | InputHisto, |
| double * | FitResults, | ||
| double | LowRange = 50, |
||
| double | HighRange = 5400 |
||
| ) | [private] |
Definition at line 330 of file SiStripGainFromCalibTree.cc.
References MinNrEntries.
Referenced by algoComputeMPVandGain().
{
FitResults[0] = -0.5; //MPV
FitResults[1] = 0; //MPV error
FitResults[2] = -0.5; //Width
FitResults[3] = 0; //Width error
FitResults[4] = -0.5; //Fit Chi2/NDF
if( InputHisto->GetEntries() < MinNrEntries)return;
// perform fit with standard landau
TF1* MyLandau = new TF1("MyLandau","landau",LowRange, HighRange);
MyLandau->SetParameter(1,300);
InputHisto->Fit("MyLandau","0QR WW");
// MPV is parameter 1 (0=constant, 1=MPV, 2=Sigma)
FitResults[0] = MyLandau->GetParameter(1); //MPV
FitResults[1] = MyLandau->GetParError(1); //MPV error
FitResults[2] = MyLandau->GetParameter(2); //Width
FitResults[3] = MyLandau->GetParError(2); //Width error
FitResults[4] = MyLandau->GetChisquare() / MyLandau->GetNDF(); //Fit Chi2/NDF
delete MyLandau;
}
| bool SiStripGainFromCalibTree::IsGoodLandauFit | ( | double * | FitResults | ) | [private] |
Definition at line 355 of file SiStripGainFromCalibTree.cc.
References MaxChi2OverNDF, and MaxMPVError.
{
if(FitResults[0] < 2 )return false;
if(FitResults[1] > MaxMPVError )return false;
if(FitResults[4] > MaxChi2OverNDF)return false;
return true;
}
| void SiStripGainFromCalibTree::MakeCalibrationMap | ( | ) | [private] |
Definition at line 703 of file SiStripGainFromCalibTree.cc.
References sistrip::APV, APVsColl, stAPVGain::CalibGain, m_calibrationPath, and useCalibration.
Referenced by algoBeginJob().
{
if(!useCalibration)return;
TChain* t1 = new TChain("SiStripCalib/APVGain");
t1->Add(m_calibrationPath.c_str());
unsigned int tree_DetId;
unsigned char tree_APVId;
double tree_Gain;
t1->SetBranchAddress("DetId" ,&tree_DetId );
t1->SetBranchAddress("APVId" ,&tree_APVId );
t1->SetBranchAddress("Gain" ,&tree_Gain );
for (unsigned int ientry = 0; ientry < t1->GetEntries(); ientry++) {
t1->GetEntry(ientry);
stAPVGain* APV = APVsColl[(tree_DetId<<3) | (unsigned int)tree_APVId];
APV->CalibGain = tree_Gain;
}
}
| void SiStripGainFromCalibTree::storeOnTree | ( | ) | [private] |
Definition at line 578 of file SiStripGainFromCalibTree.cc.
References a, sistrip::APV, stAPVGain::APVId, APVsColl, APVsCollOrdered, BAD, stAPVGain::Bin, stAPVGain::DetId, stAPVGain::Eta, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, GOOD, stAPVGain::Index, stAPVGain::isMasked, NCluster, stAPVGain::NEntries, NEvent, NTrack, NULL, OutputGains, stAPVGain::Phi, stAPVGain::PreviousGain, stAPVGain::R, stAPVGain::SubDet, tfs, stAPVGain::Thickness, stAPVGain::x, stAPVGain::y, and stAPVGain::z.
Referenced by algoComputeMPVandGain().
{
unsigned int tree_Index;
unsigned int tree_Bin;
unsigned int tree_DetId;
unsigned char tree_APVId;
unsigned char tree_SubDet;
float tree_x;
float tree_y;
float tree_z;
float tree_Eta;
float tree_R;
float tree_Phi;
float tree_Thickness;
float tree_FitMPV;
float tree_FitMPVErr;
float tree_FitWidth;
float tree_FitWidthErr;
float tree_FitChi2NDF;
double tree_Gain;
double tree_PrevGain;
double tree_NEntries;
bool tree_isMasked;
TTree* MyTree;
MyTree = tfs->make<TTree> ("APVGain","APVGain");
MyTree->Branch("Index" ,&tree_Index ,"Index/i");
MyTree->Branch("Bin" ,&tree_Bin ,"Bin/i");
MyTree->Branch("DetId" ,&tree_DetId ,"DetId/i");
MyTree->Branch("APVId" ,&tree_APVId ,"APVId/b");
MyTree->Branch("SubDet" ,&tree_SubDet ,"SubDet/b");
MyTree->Branch("x" ,&tree_x ,"x/F");
MyTree->Branch("y" ,&tree_y ,"y/F");
MyTree->Branch("z" ,&tree_z ,"z/F");
MyTree->Branch("Eta" ,&tree_Eta ,"Eta/F");
MyTree->Branch("R" ,&tree_R ,"R/F");
MyTree->Branch("Phi" ,&tree_Phi ,"Phi/F");
MyTree->Branch("Thickness" ,&tree_Thickness ,"Thickness/F");
MyTree->Branch("FitMPV" ,&tree_FitMPV ,"FitMPV/F");
MyTree->Branch("FitMPVErr" ,&tree_FitMPVErr ,"FitMPVErr/F");
MyTree->Branch("FitWidth" ,&tree_FitWidth ,"FitWidth/F");
MyTree->Branch("FitWidthErr" ,&tree_FitWidthErr,"FitWidthErr/F");
MyTree->Branch("FitChi2NDF" ,&tree_FitChi2NDF ,"FitChi2NDF/F");
MyTree->Branch("Gain" ,&tree_Gain ,"Gain/D");
MyTree->Branch("PrevGain" ,&tree_PrevGain ,"PrevGain/D");
MyTree->Branch("NEntries" ,&tree_NEntries ,"NEntries/D");
MyTree->Branch("isMasked" ,&tree_isMasked ,"isMasked/O");
FILE* Gains = fopen(OutputGains.c_str(),"w");
fprintf(Gains,"NEvents = %i\n",NEvent);
fprintf(Gains,"NTracks = %i\n",NTrack);
fprintf(Gains,"NClusters = %i\n",NCluster);
fprintf(Gains,"Number of APVs = %lu\n",static_cast<unsigned long>(APVsColl.size()));
fprintf(Gains,"GoodFits = %i BadFits = %i ratio = %f\n",GOOD,BAD,(100.0*GOOD)/(GOOD+BAD));
for(unsigned int a=0;a<APVsCollOrdered.size();a++){
stAPVGain* APV = APVsCollOrdered[a];
if(APV==NULL)continue;
// printf( "%i | %i | PreviousGain = %7.5f NewGain = %7.5f\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain);
fprintf(Gains,"%i | %i | PreviousGain = %7.5f NewGain = %7.5f\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain);
tree_Index = APV->Index;
tree_Bin = APV->Bin;
tree_DetId = APV->DetId;
tree_APVId = APV->APVId;
tree_SubDet = APV->SubDet;
tree_x = APV->x;
tree_y = APV->y;
tree_z = APV->z;
tree_Eta = APV->Eta;
tree_R = APV->R;
tree_Phi = APV->Phi;
tree_Thickness = APV->Thickness;
tree_FitMPV = APV->FitMPV;
tree_FitMPVErr = APV->FitMPVErr;
tree_FitWidth = APV->FitWidth;
tree_FitWidthErr= APV->FitWidthErr;
tree_FitChi2NDF = APV->FitChi2;
tree_Gain = APV->Gain;
tree_PrevGain = APV->PreviousGain;
tree_NEntries = APV->NEntries;
tree_isMasked = APV->isMasked;
MyTree->Fill();
}
fclose(Gains);
}
bool SiStripGainFromCalibTree::AllowSaturation [private] |
Definition at line 149 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual > SiStripGainFromCalibTree::APVsColl [private] |
Definition at line 187 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), algoBeginJob(), algoComputeMPVandGain(), MakeCalibrationMap(), and storeOnTree().
std::vector<stAPVGain*> SiStripGainFromCalibTree::APVsCollOrdered [private] |
Definition at line 186 of file SiStripGainFromCalibTree.cc.
Referenced by algoBeginJob(), getNewObject(), and storeOnTree().
unsigned int SiStripGainFromCalibTree::BAD [private] |
Definition at line 178 of file SiStripGainFromCalibTree.cc.
Referenced by algoBeginJob(), algoComputeMPVandGain(), and storeOnTree().
int SiStripGainFromCalibTree::CalibrationLevel [private] |
Definition at line 153 of file SiStripGainFromCalibTree.cc.
Referenced by algoComputeMPVandGain(), and SiStripGainFromCalibTree().
TH2F* SiStripGainFromCalibTree::Charge_Vs_Index [private] |
Definition at line 161 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), algoBeginJob(), and algoComputeMPVandGain().
TH2F* SiStripGainFromCalibTree::Charge_Vs_Index_Absolute [private] |
Definition at line 162 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTECM1 [private] |
Definition at line 169 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTECM2 [private] |
Definition at line 170 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTECP1 [private] |
Definition at line 167 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTECP2 [private] |
Definition at line 168 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTIB [private] |
Definition at line 163 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTIDM [private] |
Definition at line 166 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTIDP [private] |
Definition at line 165 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
TH2F* SiStripGainFromCalibTree::Charge_Vs_PathlengthTOB [private] |
Definition at line 164 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
unsigned int SiStripGainFromCalibTree::ERun [private] |
Definition at line 176 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
bool SiStripGainFromCalibTree::FirstSetOfConstants [private] |
Definition at line 150 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), algoBeginJob(), and SiStripGainFromCalibTree().
unsigned int SiStripGainFromCalibTree::GOOD [private] |
Definition at line 177 of file SiStripGainFromCalibTree.cc.
Referenced by algoBeginJob(), algoComputeMPVandGain(), and storeOnTree().
string SiStripGainFromCalibTree::m_calibrationPath [private] |
Definition at line 156 of file SiStripGainFromCalibTree.cc.
Referenced by MakeCalibrationMap(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MaxChi2OverNDF [private] |
Definition at line 141 of file SiStripGainFromCalibTree.cc.
Referenced by IsGoodLandauFit(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MaxMPVError [private] |
Definition at line 140 of file SiStripGainFromCalibTree.cc.
Referenced by IsGoodLandauFit(), and SiStripGainFromCalibTree().
unsigned int SiStripGainFromCalibTree::MaxNrStrips [private] |
Definition at line 146 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MaxTrackChiOverNdf [private] |
Definition at line 148 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MaxTrackEta [private] |
Definition at line 145 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MaxTrackMomentum [private] |
Definition at line 143 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MinNrEntries [private] |
Definition at line 139 of file SiStripGainFromCalibTree.cc.
Referenced by getPeakOfLandau(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MinTrackEta [private] |
Definition at line 144 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
unsigned int SiStripGainFromCalibTree::MinTrackHits [private] |
Definition at line 147 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
double SiStripGainFromCalibTree::MinTrackMomentum [private] |
Definition at line 142 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
unsigned int SiStripGainFromCalibTree::NCluster [private] |
Definition at line 174 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), algoBeginJob(), and storeOnTree().
unsigned int SiStripGainFromCalibTree::NEvent [private] |
Definition at line 172 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), algoBeginJob(), and storeOnTree().
unsigned int SiStripGainFromCalibTree::NTrack [private] |
Definition at line 173 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), algoBeginJob(), and storeOnTree().
bool SiStripGainFromCalibTree::OldGainRemoving [private] |
Definition at line 152 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
std::string SiStripGainFromCalibTree::OutputGains [private] |
Definition at line 158 of file SiStripGainFromCalibTree.cc.
Referenced by SiStripGainFromCalibTree(), and storeOnTree().
unsigned int SiStripGainFromCalibTree::SRun [private] |
Definition at line 175 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and algoBeginJob().
Definition at line 137 of file SiStripGainFromCalibTree.cc.
Referenced by algoBeginJob(), and storeOnTree().
bool SiStripGainFromCalibTree::useCalibration [private] |
Definition at line 155 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), MakeCalibrationMap(), and SiStripGainFromCalibTree().
bool SiStripGainFromCalibTree::Validation [private] |
Definition at line 151 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
vector<string> SiStripGainFromCalibTree::VInputFiles [private] |
Definition at line 159 of file SiStripGainFromCalibTree.cc.
Referenced by algoAnalyzeTheTree(), and SiStripGainFromCalibTree().
1.7.3