|
|
|
Definition at line 85 of file SiStripGainFromData.cc.
| SiStripGainFromData::SiStripGainFromData | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 292 of file SiStripGainFromData.cc.
References AlgoMode, AllowSaturation, CalibrationLevel, CheckIfFileExist, CheckLocalAngle, dqmStore_, dqmStore_infile, FirstSetOfConstants, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), MaxChi2OverNDF, MaxMPVError, MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackMomentum, MinNrEntries, MinTrackEta, MinTrackHits, MinTrackMomentum, cppFunctionSkipper::operator, OutputGains, OutputHistos, TrajToTrackLabel, TrajToTrackProducer, Validation, and VInputFiles.
: ConditionDBWriter<SiStripApvGain>(iConfig) { AlgoMode = iConfig.getParameter<std::string>("AlgoMode"); OutputGains = iConfig.getParameter<std::string>("OutputGains"); OutputHistos = iConfig.getParameter<std::string>("OutputHistos"); TrajToTrackProducer = iConfig.getParameter<std::string>("TrajToTrackProducer"); TrajToTrackLabel = iConfig.getParameter<std::string>("TrajToTrackLabel"); CheckLocalAngle = iConfig.getUntrackedParameter<bool> ("checkLocalAngle" , false); MinNrEntries = iConfig.getUntrackedParameter<unsigned>("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); CheckIfFileExist = iConfig.getUntrackedParameter<bool> ("CheckIfFileExist" , false); CalibrationLevel = iConfig.getUntrackedParameter<int> ("CalibrationLevel" , 0); if( strcmp(AlgoMode.c_str(),"WriteOnDB")==0 ) VInputFiles = iConfig.getParameter<vector<string> >("VInputFiles"); dqmStore_ = edm::Service<DQMStore>().operator->(); dqmStore_infile = edm::Service<DQMStore>().operator->(); //if( OutputHistos!="" ) // dqmStore_->open(OutputHistos.c_str(), true); }
| SiStripGainFromData::~SiStripGainFromData | ( | ) |
Definition at line 332 of file SiStripGainFromData.cc.
{
}
| void SiStripGainFromData::algoAnalyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Reimplemented from ConditionDBWriter< SiStripApvGain >.
Definition at line 1000 of file SiStripGainFromData.cc.
References AlgoMode, edm::AssociationMap< Tag >::begin(), Trajectory::chiSquared(), SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), ComputeChargeOverPath(), EEvent, edm::AssociationMap< Tag >::end(), ERun, reco::TrackBase::eta(), ETimestamp, edm::EventID::event(), Trajectory::foundHits(), edm::Event::getByLabel(), HTrackChi2OverNDF, HTrackHits, edm::EventBase::id(), iEvent, iEvent_, TrajectoryStateOnSurface::isValid(), MaxTrackChiOverNdf, MaxTrackEta, MaxTrackMomentum, Trajectory::measurements(), MinTrackEta, MinTrackHits, MinTrackMomentum, SiStripMatchedRecHit2D::monoCluster(), NEvent, reco::TrackBase::p(), edm::Handle< T >::product(), reco::TrackBase::pt(), Trajectory::recHits(), edm::EventID::run(), SEvent, SRun, SiStripMatchedRecHit2D::stereoCluster(), STimestamp, edm::EventBase::time(), Tracks_P_Vs_Eta, Tracks_Pt_Vs_Eta, TrajToTrackLabel, TrajToTrackProducer, and edm::Timestamp::value().
{
if( strcmp(AlgoMode.c_str(),"WriteOnDB")==0 ) return;
if(NEvent==0){
SRun = iEvent.id().run();
SEvent = iEvent.id().event();
STimestamp = iEvent.time().value();
}
ERun = iEvent.id().run();
EEvent = iEvent.id().event();
ETimestamp = iEvent.time().value();
NEvent++;
iEvent_ = &iEvent;
Handle<TrajTrackAssociationCollection> trajTrackAssociationHandle;
iEvent.getByLabel(TrajToTrackProducer, TrajToTrackLabel, trajTrackAssociationHandle);
const TrajTrackAssociationCollection TrajToTrackMap = *trajTrackAssociationHandle.product();
for(TrajTrackAssociationCollection::const_iterator it = TrajToTrackMap.begin(); it!=TrajToTrackMap.end(); ++it) {
const Track track = *it->val;
const Trajectory traj = *it->key;
if(track.p()<MinTrackMomentum || track.p()>MaxTrackMomentum || track.eta()<MinTrackEta || track.eta()>MaxTrackEta) continue;
Tracks_Pt_Vs_Eta->Fill(fabs(track.eta()),track.pt());
Tracks_P_Vs_Eta ->Fill(fabs(track.eta()),track.p());
//BEGIN TO COMPUTE NDOF FOR TRACKS NO IMPLEMENTED BEFORE 200pre3
int ndof =0;
const Trajectory::RecHitContainer transRecHits = traj.recHits();
for(Trajectory::RecHitContainer::const_iterator rechit = transRecHits.begin(); rechit != transRecHits.end(); ++rechit)
if ((*rechit)->isValid()) ndof += (*rechit)->dimension();
ndof -= 5;
//END TO COMPUTE NDOF FOR TRACKS NO IMPLEMENTED BEFORE 200pre3
HTrackChi2OverNDF->Fill(traj.chiSquared()/ndof);
if(traj.chiSquared()/ndof>MaxTrackChiOverNdf)continue;
vector<TrajectoryMeasurement> measurements = traj.measurements();
HTrackHits->Fill(traj.foundHits());
if(traj.foundHits()<(int)MinTrackHits)continue;
/*
//BEGIN TO COMPUTE #MATCHEDRECHIT IN THE TRACK
int NMatchedHit = 0;
for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
TrajectoryStateOnSurface trajState = measurement_it->updatedState();
if( !trajState.isValid() ) continue;
const TrackingRecHit* hit = (*measurement_it->recHit()).hit();
const SiStripMatchedRecHit2D* sistripmatchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
if(sistripmatchedhit)NMatchedHit++;
// NMatchedHit++;
}
//END TO COMPUTE #MATCHEDRECHIT IN THE TRACK
if(NMatchedHit<2){
printf("NOT ENOUGH MATCHED RECHITS : %i\n",NMatchedHit);
continue;
}
*/
for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
TrajectoryStateOnSurface trajState = measurement_it->updatedState();
if( !trajState.isValid() ) continue;
const TrackingRecHit* hit = (*measurement_it->recHit()).hit();
const SiStripRecHit1D* sistripsimple1dhit = dynamic_cast<const SiStripRecHit1D*>(hit);
const SiStripRecHit2D* sistripsimplehit = dynamic_cast<const SiStripRecHit2D*>(hit);
const SiStripMatchedRecHit2D* sistripmatchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
if(sistripsimplehit){
ComputeChargeOverPath((sistripsimplehit->cluster()).get(), trajState, &iSetup, &track, traj.chiSquared()/ndof);
}else if(sistripmatchedhit){
ComputeChargeOverPath(&sistripmatchedhit->monoCluster(),trajState, &iSetup, &track, traj.chiSquared()/ndof);
ComputeChargeOverPath(&sistripmatchedhit->stereoCluster(),trajState, &iSetup, &track, traj.chiSquared()/ndof);
}else if(sistripsimple1dhit){
ComputeChargeOverPath((sistripsimple1dhit->cluster()).get(), trajState, &iSetup, &track, traj.chiSquared()/ndof);
}else{
}
}
}
}
| void SiStripGainFromData::algoBeginJob | ( | const edm::EventSetup & | iSetup | ) | [private, virtual] |
Reimplemented from ConditionDBWriter< SiStripApvGain >.
Definition at line 338 of file SiStripGainFromData.cc.
References AlgoMode, sistrip::APV, APV_Charge, APV_CumulGain, APV_DetId, APV_Eta, APV_Gain, APV_Id, APV_Momentum, APV_MPV, APV_PathLength, APV_PathLengthM, APV_PrevGain, APV_R, APV_SubDet, APV_Thickness, stAPVGain::APVId, APVsColl, APVsCollOrdered, PV3DBase< T, PVType, FrameType >::basicVector(), DQMStore::book1D(), DQMStore::book2D(), BoundSurface::bounds(), Charge_TEC, Charge_TEC1, Charge_TEC2, Charge_TECM, Charge_TECP, Charge_TIB, Charge_TID, Charge_TIDM, Charge_TIDP, Charge_TOB, Charge_Vs_Alpha, Charge_Vs_Beta, Charge_Vs_PathLength, Charge_Vs_PathLength320, Charge_Vs_PathLength500, Charge_Vs_PathTEC, Charge_Vs_PathTEC1, Charge_Vs_PathTEC2, Charge_Vs_PathTIB, Charge_Vs_PathTID, Charge_Vs_PathTOB, Charge_Vs_TransversAngle, stAPVGain::DetId, cond::rpcobgas::detid, dqmStore_, EEvent, Error_Vs_Entries, Error_Vs_Eta, Error_Vs_MPV, Error_Vs_Phi, ERun, reco::tau::disc::Eta(), stAPVGain::Eta, ETimestamp, cmsRelvalreport::exit, FWHM_Vs_PathLength, FWHM_Vs_PathLength320, FWHM_Vs_PathLength500, stAPVGain::Gain, edm::EventSetup::get(), MonitorElement::getTH1F(), MonitorElement::getTH2F(), HChi2OverNDF, HFirstStrip, HTrackChi2OverNDF, HTrackHits, i, stAPVGain::Index, iSetup_, edm::ESHandleBase::isValid(), j, JobInfo, stAPVGain::MPV, MPV_Vs_Alpha, MPV_Vs_Beta, MPV_Vs_Eta, MPV_Vs_EtaTEC, MPV_Vs_EtaTEC1, MPV_Vs_EtaTEC2, MPV_Vs_EtaTIB, MPV_Vs_EtaTID, MPV_Vs_EtaTOB, MPV_Vs_PathLength, MPV_Vs_PathLength320, MPV_Vs_PathLength500, MPV_Vs_PathTEC, MPV_Vs_PathTEC1, MPV_Vs_PathTEC2, MPV_Vs_PathTIB, MPV_Vs_PathTID, MPV_Vs_PathTOB, MPV_Vs_Phi, MPV_Vs_PhiTEC, MPV_Vs_PhiTEC1, MPV_Vs_PhiTEC2, MPV_Vs_PhiTIB, MPV_Vs_PhiTID, MPV_Vs_PhiTOB, MPV_Vs_R, MPV_Vs_TransversAngle, MPVs, MPVs320, MPVs500, NEvent, NHighStripInCluster, NoMPV_Vs_EtaPhi, NSatStripInCluster, StripTopology::nstrips(), NStrips_Vs_Alpha, NStrips_Vs_Beta, NStrips_Vs_TransversAngle, NumberOfEntriesByAPV, colinearityKinematic::Phi, stAPVGain::Phi, GeomDet::position(), stAPVGain::PreviousGain, dttmaxenums::R, stAPVGain::R, DetId::rawId(), SEvent, stAPVGain::Side, TIDDetId::side(), TECDetId::side(), StripGeomDetUnit::specificTopology(), SRun, STimestamp, stAPVGain::SubDet, DetId::subdetId(), GeomDet::surface(), StripSubdetector::TEC, stAPVGain::Thickness, Bounds::thickness(), StripSubdetector::TIB, StripSubdetector::TID, tmp, StripSubdetector::TOB, Tracks_P_Vs_Eta, and Tracks_Pt_Vs_Eta.
{
iSetup_ = &iSetup;
// TH1::AddDirectory(kTRUE);
tmp = dqmStore_->book1D ("JobInfo" , "JobInfo", 20,0,20); JobInfo = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_DetId" , "APV_DetId" , 72785,0,72784); APV_DetId = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_Id" , "APV_Id" , 72785,0,72784); APV_Id = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_Eta" , "APV_Eta" , 72785,0,72784); APV_Eta = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_R" , "APV_R" , 72785,0,72784); APV_R = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_SubDet" , "APV_SubDet" , 72785,0,72784); APV_SubDet = tmp->getTH1F();
tmp = dqmStore_->book2D ("APV_Momentum" , "APV_Momentum" , 72785,0,72784, 50,0,200); APV_Momentum = tmp->getTH2F();
tmp = dqmStore_->book2D ("APV_Charge" , "APV_Charge" , 72785,0,72784, 1000,0,2000); APV_Charge = tmp->getTH2F();
tmp = dqmStore_->book2D ("APV_PathLength" , "APV_PathLength" , 72785,0,72784, 100,0.2,1.4); APV_PathLength = tmp->getTH2F();
tmp = dqmStore_->book1D ("APV_PathLengthM", "APV_PathLengthM", 72785,0,72784); APV_PathLengthM = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_MPV" , "APV_MPV" , 72785,0,72784); APV_MPV = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_Gain" , "APV_Gain" , 72785,0,72784); APV_Gain = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_PrevGain" , "APV_PrevGain" , 72785,0,72784); APV_PrevGain = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_CumulGain" , "APV_CumulGain" , 72785,0,72784); APV_CumulGain = tmp->getTH1F();
tmp = dqmStore_->book1D ("APV_Thickness" , "APV_Thicknes" , 72785,0,72784); APV_Thickness = tmp->getTH1F();
tmp = dqmStore_->book2D ("Tracks_P_Vs_Eta" , "Tracks_P_Vs_Eta" , 30, 0,3,50,0,200); Tracks_P_Vs_Eta = tmp->getTH2F();
tmp = dqmStore_->book2D ("Tracks_Pt_Vs_Eta" , "Tracks_Pt_Vs_Eta", 30, 0,3,50,0,200); Tracks_Pt_Vs_Eta = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathTIB" , "Charge_Vs_PathTIB" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTIB = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathTID" , "Charge_Vs_PathTID" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTID = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathTOB" , "Charge_Vs_PathTOB" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTOB = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathTEC" , "Charge_Vs_PathTEC" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTEC = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathTEC1", "Charge_Vs_PathTEC1",100,0.2,1.4, 500,0,2000); Charge_Vs_PathTEC1 = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathTEC2", "Charge_Vs_PathTEC2",100,0.2,1.4, 500,0,2000); Charge_Vs_PathTEC2 = tmp->getTH2F();
tmp = dqmStore_->book1D ("Charge_TIB" , "Charge_TIB" ,1000,0,2000); Charge_TIB = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TID" , "Charge_TID" ,1000,0,2000); Charge_TID = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TID+", "Charge_TID+",1000,0,2000); Charge_TIDP = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TID-", "Charge_TID-",1000,0,2000); Charge_TIDM = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TOB" , "Charge_TOB" ,1000,0,2000); Charge_TOB = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TEC" , "Charge_TEC" ,1000,0,2000); Charge_TEC = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TEC1", "Charge_TEC1",1000,0,2000); Charge_TEC1 = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TEC2", "Charge_TEC2",1000,0,2000); Charge_TEC2 = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TEC+", "Charge_TEC+",1000,0,2000); Charge_TECP = tmp->getTH1F();
tmp = dqmStore_->book1D ("Charge_TEC-", "Charge_TEC-",1000,0,2000); Charge_TECM = tmp->getTH1F();
/*
tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS1", "Charge_Vs_PathLength_CS1" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS1 = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS2", "Charge_Vs_PathLength_CS2" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS2 = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS3", "Charge_Vs_PathLength_CS3" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS3 = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS4", "Charge_Vs_PathLength_CS4" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS4 = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS5", "Charge_Vs_PathLength_CS5" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS5 = tmp->getTH2F();
*/
tmp = dqmStore_->book2D ("Charge_Vs_PathLength" , "Charge_Vs_PathLength" , 100,0.2,1.4, 500,0,2000); Charge_Vs_PathLength = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathLength320" , "Charge_Vs_PathLength" , 100,0.2,1.4, 500,0,2000); Charge_Vs_PathLength320 = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_PathLength500" , "Charge_Vs_PathLength" , 100,0.2,1.4, 500,0,2000); Charge_Vs_PathLength500 = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_TransversAngle" , "Charge_Vs_TransversAngle" , 220,-20,200, 500,0,2000); Charge_Vs_TransversAngle = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_Alpha" , "Charge_Vs_Alpha" , 220,-20,200, 500,0,2000); Charge_Vs_Alpha = tmp->getTH2F();
tmp = dqmStore_->book2D ("Charge_Vs_Beta" , "Charge_Vs_Beta" , 220,-20,200, 500,0,2000); Charge_Vs_Beta = tmp->getTH2F();
tmp = dqmStore_->book2D ("NStrips_Vs_TransversAngle", "NStrips_Vs_TransversAngle", 220,-20,200, 10,0,10); NStrips_Vs_TransversAngle = tmp->getTH2F();
tmp = dqmStore_->book2D ("NStrips_Vs_Alpha" , "NStrips_Vs_Alpha" , 220,-20,200, 10,0,10); NStrips_Vs_Alpha = tmp->getTH2F();
tmp = dqmStore_->book2D ("NStrips_Vs_Beta" , "NStrips_Vs_Beta" , 220,-20,200, 10,0,10); NStrips_Vs_Beta = tmp->getTH2F();
tmp = dqmStore_->book1D ("NHighStripInCluster" , "NHighStripInCluster" , 15,0,14); NHighStripInCluster = tmp->getTH1F();
tmp = dqmStore_->book1D ("NSatStripInCluster" , "NSatStripInCluster" , 50,0,50); NSatStripInCluster = tmp->getTH1F();
tmp = dqmStore_->book1D ("TrackChi2OverNDF","TrackChi2OverNDF", 500, 0,10); HTrackChi2OverNDF = tmp->getTH1F();
tmp = dqmStore_->book1D ("TrackHits","TrackHits", 40, 0,40); HTrackHits = tmp->getTH1F();
tmp = dqmStore_->book1D ("FirstStrip","FirstStrip", 800, 0,800); HFirstStrip = tmp->getTH1F();
if( strcmp(AlgoMode.c_str(),"MultiJob")!=0 ){
tmp = dqmStore_->book2D ("MPV_Vs_EtaTIB" , "MPV_Vs_EtaTIB" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTIB = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_EtaTID" , "MPV_Vs_EtaTID" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTID = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_EtaTOB" , "MPV_Vs_EtaTOB" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTOB = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_EtaTEC" , "MPV_Vs_EtaTEC" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTEC = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_EtaTEC1" , "MPV_Vs_EtaTEC1", 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTEC1 = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_EtaTEC2" , "MPV_Vs_EtaTEC2", 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTEC2 = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_PhiTIB" , "MPV_Vs_PhiTIB" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTIB = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_PhiTID" , "MPV_Vs_PhiTID" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTID = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_PhiTOB" , "MPV_Vs_PhiTOB" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTOB = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_PhiTEC" , "MPV_Vs_PhiTEC" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTEC = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_PhiTEC1" , "MPV_Vs_PhiTEC1", 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTEC1 = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_PhiTEC2" , "MPV_Vs_PhiTEC2", 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTEC2 = tmp->getTH2F();
tmp = dqmStore_->book1D ("MPV_Vs_PathTIB" , "MPV_Vs_PathTIB" ,100,0.2,1.4); MPV_Vs_PathTIB = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathTID" , "MPV_Vs_PathTID" ,100,0.2,1.4); MPV_Vs_PathTID = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathTOB" , "MPV_Vs_PathTOB" ,100,0.2,1.4); MPV_Vs_PathTOB = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathTEC" , "MPV_Vs_PathTEC" ,100,0.2,1.4); MPV_Vs_PathTEC = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathTEC1" , "MPV_Vs_PathTEC1" ,100,0.2,1.4); MPV_Vs_PathTEC1 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathTEC2" , "MPV_Vs_PathTEC2" ,100,0.2,1.4); MPV_Vs_PathTEC2 = tmp->getTH1F();
tmp = dqmStore_->book2D ("MPV_Vs_Phi", "MPV_Vs_Phi", 50, -3.2, 3.2 , 300, 0, 600); MPV_Vs_Phi = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_Eta", "MPV_Vs_Eta", 50, -3.0, 3.0 , 300, 0, 600); MPV_Vs_Eta = tmp->getTH2F();
tmp = dqmStore_->book2D ("MPV_Vs_R" , "MPV_Vs_R" , 150, 0.0, 150.0, 300, 0, 600); MPV_Vs_R = tmp->getTH2F();
/*
tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS1" , "MPV_Vs_PathLength_CS1" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS1 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS2" , "MPV_Vs_PathLength_CS2" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS2 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS3" , "MPV_Vs_PathLength_CS3" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS3 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS4" , "MPV_Vs_PathLength_CS4" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS4 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS5" , "MPV_Vs_PathLength_CS5" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS5 = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS1" , "FWHM_Vs_PathLength_CS1", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS1 = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS2" , "FWHM_Vs_PathLength_CS2", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS2 = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS3" , "FWHM_Vs_PathLength_CS3", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS3 = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS4" , "FWHM_Vs_PathLength_CS4", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS4 = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS5" , "FWHM_Vs_PathLength_CS5", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS5 = tmp->getTH1F();
*/
tmp = dqmStore_->book1D ("MPV_Vs_PathLength" , "MPV_Vs_PathLength" , 100, 0.2, 1.4); MPV_Vs_PathLength = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathLength320" , "MPV_Vs_PathLength" , 100, 0.2, 1.4); MPV_Vs_PathLength320 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_PathLength500" , "MPV_Vs_PathLength" , 100, 0.2, 1.4); MPV_Vs_PathLength500 = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength" , "FWHM_Vs_PathLength" , 100, 0.2, 1.4); FWHM_Vs_PathLength = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength320" , "FWHM_Vs_PathLength" , 100, 0.2, 1.4); FWHM_Vs_PathLength320 = tmp->getTH1F();
tmp = dqmStore_->book1D ("FWHM_Vs_PathLength500" , "FWHM_Vs_PathLength" , 100, 0.2, 1.4); FWHM_Vs_PathLength500 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_TransversAngle" , "MPV_Vs_TransversAngle" , 220, -20, 200); MPV_Vs_TransversAngle = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_Alpha" , "MPV_Vs_Alpha" , 220, -20, 200); MPV_Vs_Alpha = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPV_Vs_Beta" , "MPV_Vs_Beta" , 220, -20, 200); MPV_Vs_Beta = tmp->getTH1F();
tmp = dqmStore_->book2D ("Error_Vs_MPV" , "Error_Vs_MPV" ,600,0,600 ,100 ,0 ,50); Error_Vs_MPV = tmp->getTH2F();
tmp = dqmStore_->book2D ("Error_Vs_Entries","Error_Vs_Entries",500,0,10000 ,100 ,0 ,50); Error_Vs_Entries = tmp->getTH2F();
tmp = dqmStore_->book2D ("Error_Vs_Eta" , "Error_Vs_Eta" ,50 ,-3.0,3.0 ,100 ,0 ,50 ); Error_Vs_Eta = tmp->getTH2F();
tmp = dqmStore_->book2D ("Error_Vs_Phi" , "Error_Vs_Phi" ,50 ,-3.2,3.2 ,100 ,0 ,50); Error_Vs_Phi = tmp->getTH2F();
tmp = dqmStore_->book2D ("NoMPV_Vs_EtaPhi" , "NoMPV_Vs_EtaPhi" ,50,-3.0,3.0 ,50 ,-3.2,3.2); NoMPV_Vs_EtaPhi = tmp->getTH2F();
tmp = dqmStore_->book1D ("NumberOfEntriesByAPV" , "NumberOfEntriesByAPV" , 1000, 0,10000); NumberOfEntriesByAPV = tmp->getTH1F();
tmp = dqmStore_->book1D ("Chi2OverNDF","Chi2OverNDF", 500, 0,25); HChi2OverNDF = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPVs", "MPVs", 600,0,600); MPVs = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPVs320", "MPVs320", 600,0,600); MPVs320 = tmp->getTH1F();
tmp = dqmStore_->book1D ("MPVs500", "MPVs500", 600,0,600); MPVs500 = tmp->getTH1F();
// MPV_vs_10RplusEta tmp = dqmStore_->book2D ("MPV_vs_10RplusEta","MPV_vs_10RplusEta", 48000,0,2400, 800,100,500);
}
gROOT->cd();
edm::ESHandle<TrackerGeometry> tkGeom;
iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
vector<GeomDet*> Det = tkGeom->dets();
edm::ESHandle<SiStripGain> gainHandle;
// if(strcmp(AlgoMode.c_str(),"MultiJob")!=0 && !FirstSetOfConstants){
iSetup.get<SiStripGainRcd>().get(gainHandle);
if(!gainHandle.isValid()){printf("\n#####################\n\nERROR --> gainHandle is not valid\n\n#####################\n\n");exit(0);}
// }
unsigned int Id=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;
double Phi = DetUnit->position().basicVector().phi();
double Eta = DetUnit->position().basicVector().eta();
double R = DetUnit->position().basicVector().transverse();
double Thick = DetUnit->surface().bounds().thickness();
for(unsigned int j=0;j<NAPV;j++){
stAPVGain* APV = new stAPVGain;
APV->Index = Id;
APV->DetId = Detid.rawId();
APV->APVId = j;
APV->SubDet = SubDet;
APV->MPV = -1;
APV->Gain = -1;
APV->PreviousGain = 1;
APV->Eta = Eta;
APV->Phi = Phi;
APV->R = R;
APV->Thickness = Thick;
APV->Side = 0;
if(SubDet==StripSubdetector::TID){
TIDDetId detid = TIDDetId(Detid);
APV->Side = detid.side();
}else if(SubDet==StripSubdetector::TEC){
TECDetId detid = TECDetId(Detid);
APV->Side = detid.side();
}
APVsCollOrdered.push_back(APV);
APVsColl[(APV->DetId<<3) | APV->APVId] = APV;
Id++;
APV_DetId ->Fill(Id,APV->DetId);
APV_Id ->Fill(Id,APV->APVId);
APV_Eta ->Fill(Id,APV->Eta);
APV_R ->Fill(Id,APV->R);
APV_SubDet ->Fill(Id,APV->SubDet);
APV_Thickness->Fill(Id,APV->Thickness);
}
}
}
NEvent = 0;
SRun = 0;
SEvent = 0;
STimestamp = 0;
ERun = 0;
EEvent = 0;
ETimestamp = 0;
}
| void SiStripGainFromData::algoBeginRun | ( | const edm::Run & | , |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Reimplemented from ConditionDBWriter< SiStripApvGain >.
Definition at line 971 of file SiStripGainFromData.cc.
References AlgoMode, sistrip::APV, APV_PrevGain, stAPVGain::APVId, APVsCollOrdered, stAPVGain::DetId, cmsRelvalreport::exit, FirstSetOfConstants, edm::EventSetup::get(), stAPVGain::Index, edm::ESHandleBase::isValid(), stAPVGain::PreviousGain, and Validation.
{
edm::ESHandle<SiStripGain> gainHandle;
if((strcmp(AlgoMode.c_str(),"MultiJob")!=0 && !FirstSetOfConstants) || Validation){
iSetup.get<SiStripGainRcd>().get(gainHandle);
if(!gainHandle.isValid()){printf("\n#####################\n\nERROR --> gainHandle is not valid\n\n#####################\n\n");exit(0);}
for(std::vector<stAPVGain*>::iterator it = APVsCollOrdered.begin();it!=APVsCollOrdered.end();it++){
stAPVGain* APV = *it;
if(gainHandle.isValid()){
SiStripApvGain::Range detGainRange = gainHandle->getRange(APV->DetId);
APV->PreviousGain = *(detGainRange.first + APV->APVId);
// APV_PrevGain->Fill(APV->Index,APV->PreviousGain);
APV_PrevGain->SetBinContent(APV_PrevGain->GetXaxis()->FindBin(APV->Index),APV->PreviousGain);
if(APV->PreviousGain<0)APV->PreviousGain = 1;
}else{
printf("GAIN HANDLE IS NOT VALID\n");
}
}
}
}
| void SiStripGainFromData::algoEndJob | ( | ) | [private, virtual] |
Reimplemented from ConditionDBWriter< SiStripApvGain >.
Definition at line 562 of file SiStripGainFromData.cc.
References AlgoMode, sistrip::APV, APV_Charge, APV_CumulGain, APV_Gain, APV_Momentum, APV_MPV, APV_PathLength, APV_PathLengthM, stAPVGain::APVId, APVsColl, APVsCollOrdered, b, pat::BAD, newFWLiteAna::bin, CalibrationLevel, DQMStore::cd(), Charge_TEC, Charge_TEC1, Charge_TEC2, Charge_TECM, Charge_TECP, Charge_TIB, Charge_TID, Charge_TIDM, Charge_TIDP, Charge_TOB, Charge_Vs_Alpha, Charge_Vs_Beta, Charge_Vs_PathLength, Charge_Vs_PathLength320, Charge_Vs_PathLength500, Charge_Vs_PathTEC, Charge_Vs_PathTEC1, Charge_Vs_PathTEC2, Charge_Vs_PathTIB, Charge_Vs_PathTID, Charge_Vs_PathTOB, Charge_Vs_TransversAngle, CheckIfFileExist, stAPVGain::DetId, dqmStore_, EEvent, Error_Vs_Entries, Error_Vs_Eta, Error_Vs_MPV, Error_Vs_Phi, ERun, stAPVGain::Eta, f, mergeVDriftHistosByStation::file, FirstSetOfConstants, FWHM_Vs_PathLength, FWHM_Vs_PathLength320, FWHM_Vs_PathLength500, stAPVGain::Gain, getPeakOfLandau(), pat::GOOD, HChi2OverNDF, HFirstStrip, HTrackChi2OverNDF, HTrackHits, Exhume::I, stAPVGain::Index, j, JobInfo, MaxMPVError, stAPVGain::MPV, MPV_Vs_Alpha, MPV_Vs_Beta, MPV_Vs_Eta, MPV_Vs_EtaTEC, MPV_Vs_EtaTEC1, MPV_Vs_EtaTEC2, MPV_Vs_EtaTIB, MPV_Vs_EtaTID, MPV_Vs_EtaTOB, MPV_Vs_PathLength, MPV_Vs_PathLength320, MPV_Vs_PathLength500, MPV_Vs_PathTEC, MPV_Vs_PathTEC1, MPV_Vs_PathTEC2, MPV_Vs_PathTIB, MPV_Vs_PathTID, MPV_Vs_PathTOB, MPV_Vs_Phi, MPV_Vs_PhiTEC, MPV_Vs_PhiTEC1, MPV_Vs_PhiTEC2, MPV_Vs_PhiTIB, MPV_Vs_PhiTID, MPV_Vs_PhiTOB, MPV_Vs_R, MPV_Vs_TransversAngle, MPVs, MPVs320, MPVs500, NEvent, NHighStripInCluster, NoMPV_Vs_EtaPhi, NSatStripInCluster, NStrips_Vs_Alpha, NStrips_Vs_TransversAngle, NULL, NumberOfEntriesByAPV, OutputGains, OutputHistos, stAPVGain::Phi, stAPVGain::PreviousGain, stAPVGain::R, DQMStore::save(), SEvent, stAPVGain::Side, SRun, stAPVGain::SubDet, StripSubdetector::TEC, stAPVGain::Thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, Tracks_P_Vs_Eta, Tracks_Pt_Vs_Eta, and VInputFiles.
{
unsigned int I=0;
if( strcmp(AlgoMode.c_str(),"WriteOnDB")==0 || strcmp(AlgoMode.c_str(),"Merge")==0){
TH1::AddDirectory(kTRUE);
TFile* file = NULL;
for(unsigned int f=0;f<VInputFiles.size();f++){
printf("Loading New Input File : %s\n", VInputFiles[f].c_str());
if(CheckIfFileExist){
FILE* doesFileExist = fopen( VInputFiles[f].c_str(), "r" );
if(!doesFileExist){
printf("File %s doesn't exist\n",VInputFiles[f].c_str());
continue;
}else{
fclose(doesFileExist);
}
}
file = new TFile( VInputFiles[f].c_str() ); if(!file || file->IsZombie() ){printf("### Bug With File %s\n### File will be skipped \n",VInputFiles[f].c_str()); continue;}
APV_Charge ->Add( (TH1*) file->FindObjectAny("APV_Charge") , 1);
APV_Momentum ->Add( (TH1*) file->FindObjectAny("APV_Momentum") , 1);
APV_PathLength ->Add( (TH1*) file->FindObjectAny("APV_PathLength") , 1);
Tracks_P_Vs_Eta ->Add( (TH1*) file->FindObjectAny("Tracks_P_Vs_Eta") , 1);
Tracks_Pt_Vs_Eta ->Add( (TH1*) file->FindObjectAny("Tracks_Pt_Vs_Eta") , 1);
Charge_Vs_PathTIB ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTIB") , 1);
Charge_Vs_PathTID ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTID") , 1);
Charge_Vs_PathTOB ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTOB") , 1);
Charge_Vs_PathTEC ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTEC") , 1);
Charge_Vs_PathTEC1 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTEC1") , 1);
Charge_Vs_PathTEC2 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTEC2") , 1);
HTrackChi2OverNDF ->Add( (TH1*) file->FindObjectAny("TrackChi2OverNDF") , 1);
HTrackHits ->Add( (TH1*) file->FindObjectAny("TrackHits") , 1);
NHighStripInCluster ->Add( (TH1*) file->FindObjectAny("NHighStripInCluster") , 1);
NSatStripInCluster ->Add( (TH1*) file->FindObjectAny("NSatStripInCluster") , 1);
Charge_Vs_PathLength ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathLength") , 1);
Charge_Vs_PathLength320 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathLength320") , 1);
Charge_Vs_PathLength500 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathLength500") , 1);
Charge_Vs_TransversAngle ->Add( (TH1*) file->FindObjectAny("Charge_Vs_TransversAngle") , 1);
NStrips_Vs_TransversAngle->Add( (TH1*) file->FindObjectAny("NStrips_Vs_TransversAngle"), 1);
Charge_Vs_Alpha ->Add( (TH1*) file->FindObjectAny("Charge_Vs_Alpha") , 1);
NStrips_Vs_Alpha ->Add( (TH1*) file->FindObjectAny("NStrips_Vs_Alpha") , 1);
HFirstStrip ->Add( (TH1*) file->FindObjectAny("FirstStrip") , 1);
TH1F* JobInfo_tmp = (TH1F*) file->FindObjectAny("JobInfo");
NEvent += (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(1));
unsigned int tmp_SRun = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(3));
unsigned int tmp_SEvent = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(4));
unsigned int tmp_ERun = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(6));
unsigned int tmp_EEvent = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(7));
if(SRun==0)SRun = tmp_SRun;
if(tmp_SRun< SRun){SRun=tmp_SRun; SEvent=tmp_SEvent;}
else if(tmp_SRun==SRun && tmp_SEvent<SEvent){SEvent=tmp_SEvent;}
if(tmp_ERun> ERun){ERun=tmp_ERun; EEvent=tmp_EEvent;}
else if(tmp_ERun==ERun && tmp_EEvent>EEvent){EEvent=tmp_EEvent;}
printf("Deleting Current Input File\n");
file->Close();
delete file;
}
}
JobInfo->Fill(1,NEvent);
JobInfo->Fill(3,SRun);
JobInfo->Fill(4,SEvent);
JobInfo->Fill(6,ERun);
JobInfo->Fill(7,EEvent);
if( strcmp(AlgoMode.c_str(),"MultiJob")!=0 ){
TH1D* Proj = NULL;
double* FitResults = new double[5];
I=0;
for(__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual >::iterator it = APVsColl.begin();it!=APVsColl.end();it++){
if( I%3650==0 ) printf("Fitting Histograms \t %6.2f%%\n",(100.0*I)/APVsColl.size());I++;
stAPVGain* APV = it->second;
int bin = APV_Charge->GetXaxis()->FindBin(APV->Index);
Proj = APV_Charge->ProjectionY(" ",bin,bin,"e");
Proj = (TH1D*)Proj->Clone();
if(Proj==NULL)continue;
// ADD PROJECTTIONS COMMING FROM THE SECOND APV IN THE PAIR
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];
int bin2 = APV_Charge->GetXaxis()->FindBin(APV2->Index);
TH1D* Proj2 = APV_Charge->ProjectionY(" ",bin2,bin2,"e");
if(Proj2!=NULL){
Proj->Add(Proj2,1);
}
}else if(CalibrationLevel>1){
// printf("%8i %i--> %4.0f + %4.0f\n",APV->DetId, APV->APVId, 0.0, Proj->GetEntries());
for(__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual >::iterator it2 = APVsColl.begin();it2!=APVsColl.end();it2++){
stAPVGain* APV2 = it2->second;
if(APV2->DetId != APV->DetId)continue;
if(APV2->APVId == APV->APVId)continue;
int bin2 = APV_Charge->GetXaxis()->FindBin(APV2->Index);
TH1D* Proj2 = APV_Charge->ProjectionY(" ",bin2,bin2,"e");
if(Proj2!=NULL){
// printf("%8i %i--> %4.0f + %4.0f\n",APV2->DetId, APV2->APVId, Proj->GetEntries(), Proj2->GetEntries());
Proj->Add(Proj2,1);
}
}
// printf("%8i %i--> %4.0f Full\n",APV->DetId, APV->APVId, Proj->GetEntries());
}
//std::cout << "Proj->GetEntries(): " << Proj->GetEntries() << ", Proj->GetMean(): " << Proj->GetMean() << std::endl;
getPeakOfLandau(Proj,FitResults);
APV->MPV = FitResults[0];
// printf("MPV = %f - %f\n",FitResults[0], FitResults[1]);
if(FitResults[0]!=-0.5 && FitResults[1]<MaxMPVError){
APV_MPV->Fill(APV->Index,APV->MPV);
MPVs ->Fill(APV->MPV);
if(APV->Thickness<0.04) MPVs320->Fill(APV->MPV);
if(APV->Thickness>0.04) MPVs500->Fill(APV->MPV);
MPV_Vs_R->Fill(APV->R,APV->MPV);
MPV_Vs_Eta->Fill(APV->Eta,APV->MPV);
if(APV->SubDet==StripSubdetector::TIB) MPV_Vs_EtaTIB ->Fill(APV->Eta,APV->MPV);
if(APV->SubDet==StripSubdetector::TID) MPV_Vs_EtaTID ->Fill(APV->Eta,APV->MPV);
if(APV->SubDet==StripSubdetector::TOB) MPV_Vs_EtaTOB ->Fill(APV->Eta,APV->MPV);
if(APV->SubDet==StripSubdetector::TEC){ MPV_Vs_EtaTEC ->Fill(APV->Eta,APV->MPV);
if(APV->Thickness<0.04) MPV_Vs_EtaTEC1->Fill(APV->Eta,APV->MPV);
if(APV->Thickness>0.04) MPV_Vs_EtaTEC2->Fill(APV->Eta,APV->MPV);
}
MPV_Vs_Phi->Fill(APV->Phi,APV->MPV);
if(APV->SubDet==StripSubdetector::TIB) MPV_Vs_PhiTIB ->Fill(APV->Phi,APV->MPV);
if(APV->SubDet==StripSubdetector::TID) MPV_Vs_PhiTID ->Fill(APV->Phi,APV->MPV);
if(APV->SubDet==StripSubdetector::TOB) MPV_Vs_PhiTOB ->Fill(APV->Phi,APV->MPV);
if(APV->SubDet==StripSubdetector::TEC){ MPV_Vs_PhiTEC ->Fill(APV->Phi,APV->MPV);
if(APV->Thickness<0.04) MPV_Vs_PhiTEC1->Fill(APV->Phi,APV->MPV);
if(APV->Thickness>0.04) MPV_Vs_PhiTEC2->Fill(APV->Phi,APV->MPV);
}
if(APV->SubDet==StripSubdetector::TIB) Charge_TIB ->Add(Proj,1);
if(APV->SubDet==StripSubdetector::TID){ Charge_TID ->Add(Proj,1);
if(APV->Side==1) Charge_TIDM->Add(Proj,1);
if(APV->Side==2) Charge_TIDP->Add(Proj,1);
}
if(APV->SubDet==StripSubdetector::TOB) Charge_TOB ->Add(Proj,1);
if(APV->SubDet==StripSubdetector::TEC){ Charge_TEC ->Add(Proj,1);
if(APV->Thickness<0.04) Charge_TEC1->Add(Proj,1);
if(APV->Thickness>0.04) Charge_TEC2->Add(Proj,1);
if(APV->Side==1) Charge_TECM->Add(Proj,1);
if(APV->Side==2) Charge_TECP->Add(Proj,1);
}
}
if(APV->SubDet==StripSubdetector::TIB) Charge_TIB ->Add(Proj,1);
if(APV->SubDet==StripSubdetector::TID){ Charge_TID ->Add(Proj,1);
if(APV->Side==1) Charge_TIDM->Add(Proj,1);
if(APV->Side==2) Charge_TIDP->Add(Proj,1);
}
if(APV->SubDet==StripSubdetector::TOB) Charge_TOB ->Add(Proj,1);
if(APV->SubDet==StripSubdetector::TEC){ Charge_TEC ->Add(Proj,1);
if(APV->Thickness<0.04) Charge_TEC1->Add(Proj,1);
if(APV->Thickness>0.04) Charge_TEC2->Add(Proj,1);
if(APV->Side==1) Charge_TECM->Add(Proj,1);
if(APV->Side==2) Charge_TECP->Add(Proj,1);
}
if(FitResults[0]!=-0.5){
HChi2OverNDF->Fill(FitResults[4]);
Error_Vs_MPV->Fill(FitResults[0],FitResults[1]);
Error_Vs_Entries->Fill(Proj->GetEntries(),FitResults[1]);
Error_Vs_Eta->Fill(APV->Eta,FitResults[1]);
Error_Vs_Phi->Fill(APV->Phi,FitResults[1]);
}
NumberOfEntriesByAPV->Fill(Proj->GetEntries());
delete Proj;
Proj = APV_PathLength->ProjectionY(" ",bin,bin,"e");
if(Proj==NULL)continue;
APV_PathLengthM->SetBinContent(APV->Index, Proj->GetMean(1) );
APV_PathLengthM->SetBinError (APV->Index, Proj->GetMeanError(1) );
// delete Proj;
}
unsigned int GOOD = 0;
unsigned int BAD = 0;
double MPVmean = MPVs->GetMean();
MPVmean = 300;
for(__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual >::iterator it = APVsColl.begin();it!=APVsColl.end();it++){
stAPVGain* APV = it->second;
if(APV->MPV>0){
APV->Gain = APV->MPV / MPVmean; // APV->MPV;
GOOD++;
}else{
NoMPV_Vs_EtaPhi->Fill(APV->Eta, APV->Phi);
APV->Gain = 1;
BAD++;
}
if(APV->Gain<=0) APV->Gain = 1;
APV_Gain->Fill(APV->Index,APV->Gain);
if(!FirstSetOfConstants) APV->Gain *= APV->PreviousGain;
APV_CumulGain->Fill(APV->Index,APV->Gain);
}
for(int j=0;j<Charge_Vs_PathLength->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength->SetBinError (j, FitResults[1]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j)) );
FWHM_Vs_PathLength->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j)) );
delete Proj;
}
for(int j=0;j<Charge_Vs_PathLength320->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength320->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength320->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength320->SetBinError (j, FitResults[1]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength320->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j)) );
FWHM_Vs_PathLength320->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j)) );
delete Proj;
}
for(int j=0;j<Charge_Vs_PathLength500->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength500->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength500->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength500->SetBinError (j, FitResults[1]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength500->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j) ));
FWHM_Vs_PathLength500->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j) ));
delete Proj;
}
/*
for(int j=0;j<Charge_Vs_PathLength_CS1->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength_CS1->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength_CS1->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength_CS1->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength_CS1->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j) ));
FWHM_Vs_PathLength_CS1->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j) ));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathLength_CS2->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength_CS2->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength_CS2->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength_CS2->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength_CS2->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j) ));
FWHM_Vs_PathLength_CS2->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j) ));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathLength_CS3->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength_CS3->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength_CS3->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength_CS3->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength_CS3->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j) ));
FWHM_Vs_PathLength_CS3->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j) ));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathLength_CS4->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength_CS4->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength_CS4->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength_CS4->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength_CS4->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j) ));
FWHM_Vs_PathLength_CS4->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j) ));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathLength_CS5->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathLength_CS5->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathLength_CS5->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathLength_CS5->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j));
FWHM_Vs_PathLength_CS5->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j) ));
FWHM_Vs_PathLength_CS5->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j) ));
delete Proj;
}
*/
for(int j=0;j<Charge_Vs_PathTIB->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathTIB->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathTIB->SetBinContent(j, FitResults[0]/Charge_Vs_PathTIB->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathTIB->SetBinError (j, FitResults[1]/Charge_Vs_PathTIB->GetXaxis()->GetBinCenter(j));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathTID->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathTID->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathTID->SetBinContent(j, FitResults[0]/Charge_Vs_PathTID->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathTID->SetBinError (j, FitResults[1]/Charge_Vs_PathTID->GetXaxis()->GetBinCenter(j));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathTOB->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathTOB->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathTOB->SetBinContent(j, FitResults[0]/Charge_Vs_PathTOB->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathTOB->SetBinError (j, FitResults[1]/Charge_Vs_PathTOB->GetXaxis()->GetBinCenter(j));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathTEC->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathTEC->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathTEC->SetBinContent(j, FitResults[0]/Charge_Vs_PathTEC->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathTEC->SetBinError (j, FitResults[1]/Charge_Vs_PathTEC->GetXaxis()->GetBinCenter(j));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathTEC1->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathTEC1->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathTEC1->SetBinContent(j, FitResults[0]/Charge_Vs_PathTEC1->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathTEC1->SetBinError (j, FitResults[1]/Charge_Vs_PathTEC1->GetXaxis()->GetBinCenter(j));
delete Proj;
}
for(int j=0;j<Charge_Vs_PathTEC2->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_PathTEC2->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_PathTEC2->SetBinContent(j, FitResults[0]/Charge_Vs_PathTEC2->GetXaxis()->GetBinCenter(j));
MPV_Vs_PathTEC2->SetBinError (j, FitResults[1]/Charge_Vs_PathTEC2->GetXaxis()->GetBinCenter(j));
delete Proj;
}
for(int j=1;j<Charge_Vs_TransversAngle->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_TransversAngle->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_TransversAngle->SetBinContent(j, FitResults[0]);
MPV_Vs_TransversAngle->SetBinError (j, FitResults[1]);
delete Proj;
}
for(int j=1;j<Charge_Vs_Alpha->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_Alpha->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_Alpha->SetBinContent(j, FitResults[0]);
MPV_Vs_Alpha->SetBinError (j, FitResults[1]);
delete Proj;
}
for(int j=1;j<Charge_Vs_Beta->GetXaxis()->GetNbins();j++){
Proj = Charge_Vs_Beta->ProjectionY(" ",j,j,"e");
getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
MPV_Vs_Beta->SetBinContent(j, FitResults[0]);
MPV_Vs_Beta->SetBinError (j, FitResults[1]);
delete Proj;
}
FILE* Gains = fopen(OutputGains.c_str(),"w");
fprintf(Gains,"NEvents = %i\n",NEvent);
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(std::vector<stAPVGain*>::iterator it = APVsCollOrdered.begin();it!=APVsCollOrdered.end();it++){
stAPVGain* APV = *it;
fprintf(Gains,"%i | %i | PreviousGain = %7.5f NewGain = %7.5f\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain);
}
std::vector<int> DetIdOfBuggedAPV;
fprintf(Gains,"----------------------------------------------------------------------\n");
for(std::vector<stAPVGain*>::iterator it = APVsCollOrdered.begin();it!=APVsCollOrdered.end();it++){
stAPVGain* APV = *it;
if(APV->MPV>0 && APV->MPV<200){
bool tmpBug = false;
for(unsigned int b=0;b<DetIdOfBuggedAPV.size()&&!tmpBug;b++){if(DetIdOfBuggedAPV[b]==APV->DetId)tmpBug=true;}
if(!tmpBug){fprintf(Gains,"%i,\n",APV->DetId);DetIdOfBuggedAPV.push_back(APV->DetId);}
}
}
fclose(Gains);
// delete [] FitResults;
// delete Proj;
}
dqmStore_->cd();
dqmStore_->save(OutputHistos.c_str());
}
| double SiStripGainFromData::ComputeChargeOverPath | ( | const SiStripCluster * | Cluster, |
| TrajectoryStateOnSurface | trajState, | ||
| const edm::EventSetup * | iSetup, | ||
| const Track * | track, | ||
| double | trajChi2OverN | ||
| ) | [private] |
Definition at line 1093 of file SiStripGainFromData.cc.
References a, AllowSaturation, alpha, SiStripCluster::amplitudes(), sistrip::APV, APV_Charge, APV_Momentum, APV_PathLength, APVsColl, beta, Charge_Vs_Alpha, Charge_Vs_Beta, Charge_Vs_PathLength, Charge_Vs_PathLength320, Charge_Vs_PathLength500, Charge_Vs_PathTEC, Charge_Vs_PathTEC1, Charge_Vs_PathTEC2, Charge_Vs_PathTIB, Charge_Vs_PathTID, Charge_Vs_PathTOB, Charge_Vs_TransversAngle, SiStripCluster::firstStrip(), SiStripCluster::geographicalId(), TrajectoryStateOnSurface::globalMomentum(), HFirstStrip, stAPVGain::Index, IsFarFromBorder(), TrajectoryStateOnSurface::localDirection(), PV3DBase< T, PVType, FrameType >::mag(), MaxNrStrips, NHighStripInCluster, NSatStripInCluster, NStrips_Vs_Alpha, NStrips_Vs_Beta, NStrips_Vs_TransversAngle, getHLTPrescaleColumns::path, funct::pow(), stAPVGain::PreviousGain, mathSSE::sqrt(), stAPVGain::SubDet, StripSubdetector::TEC, stAPVGain::Thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, Validation, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by algoAnalyze().
{
LocalVector trackDirection = trajState.localDirection();
double cosine = trackDirection.z()/trackDirection.mag();
// const SiStripCluster* Cluster = (sistripsimplehit->cluster()).get();
// const vector<uint16_t>& Ampls = Cluster->amplitudes();
const vector<uint8_t>& Ampls = Cluster->amplitudes();
uint32_t DetId = Cluster->geographicalId();
int FirstStrip = Cluster->firstStrip();
int APVId = FirstStrip/128;
stAPVGain* APV = APVsColl[(DetId<<3) | APVId];
int Saturation = 0;
bool Overlaping = false;
int Charge = 0;
unsigned int NHighStrip = 0;
if(!IsFarFromBorder(trajState,DetId, iSetup))return -1;
if(FirstStrip==0 )Overlaping=true;
if(FirstStrip==128 )Overlaping=true;
if(FirstStrip==256 )Overlaping=true;
if(FirstStrip==384 )Overlaping=true;
if(FirstStrip==512 )Overlaping=true;
if(FirstStrip==640 )Overlaping=true;
if(FirstStrip<=127 && FirstStrip+Ampls.size()>127)Overlaping=true;
if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlaping=true;
if(FirstStrip<=383 && FirstStrip+Ampls.size()>383)Overlaping=true;
if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlaping=true;
if(FirstStrip<=639 && FirstStrip+Ampls.size()>639)Overlaping=true;
if(FirstStrip+Ampls.size()==127 )Overlaping=true;
if(FirstStrip+Ampls.size()==255 )Overlaping=true;
if(FirstStrip+Ampls.size()==383 )Overlaping=true;
if(FirstStrip+Ampls.size()==511 )Overlaping=true;
if(FirstStrip+Ampls.size()==639 )Overlaping=true;
if(FirstStrip+Ampls.size()==767 )Overlaping=true;
if(Overlaping)return -1;
/*
if(FirstStrip==0 )Overlaping=true;
if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlaping=true;
if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlaping=true;
if(FirstStrip+Ampls.size()==511 )Overlaping=true;
if(FirstStrip+Ampls.size()==767 )Overlaping=true;
if(Overlaping)return -1;
*/
for(unsigned int a=0;a<Ampls.size();a++){Charge+=Ampls[a];if(Ampls[a]>=254)Saturation++;if(Ampls[a]>=20)NHighStrip++;}
double path = (10.0*APV->Thickness)/fabs(cosine);
double ClusterChargeOverPath = (double)Charge / path ;
NSatStripInCluster->Fill(Saturation);
if(Ampls.size()>MaxNrStrips) return -1;
if(Saturation>0 && !AllowSaturation)return -1;
Charge_Vs_PathLength ->Fill(path,Charge);
if(APV->Thickness<0.04) Charge_Vs_PathLength320->Fill(path,Charge);
if(APV->Thickness>0.04) Charge_Vs_PathLength500->Fill(path,Charge);
if(APV->SubDet==StripSubdetector::TIB) Charge_Vs_PathTIB ->Fill(path,Charge);
if(APV->SubDet==StripSubdetector::TID) Charge_Vs_PathTID ->Fill(path,Charge);
if(APV->SubDet==StripSubdetector::TOB) Charge_Vs_PathTOB ->Fill(path,Charge);
if(APV->SubDet==StripSubdetector::TEC){ Charge_Vs_PathTEC ->Fill(path,Charge);
if(APV->Thickness<0.04) Charge_Vs_PathTEC1 ->Fill(path,Charge);
if(APV->Thickness>0.04) Charge_Vs_PathTEC2 ->Fill(path,Charge); }
double trans = atan2(trackDirection.y(),trackDirection.x())*(180/3.14159265);
double alpha = acos(trackDirection.x() / sqrt( pow(trackDirection.x(),2) + pow(trackDirection.z(),2) ) ) * (180/3.14159265);
double beta = acos(trackDirection.y() / sqrt( pow(trackDirection.x(),2) + pow(trackDirection.z(),2) ) ) * (180/3.14159265);
if(path>0.4 && path<0.45){
Charge_Vs_TransversAngle->Fill(trans,Charge/path);
Charge_Vs_Alpha->Fill(alpha,Charge/path);
Charge_Vs_Beta->Fill(beta,Charge/path);
}
NStrips_Vs_TransversAngle->Fill(trans,Ampls.size());
NStrips_Vs_Alpha ->Fill(alpha,Ampls.size());
NStrips_Vs_Beta ->Fill(beta ,Ampls.size());
NHighStripInCluster->Fill(NHighStrip);
// if(NHighStrip==1) Charge_Vs_PathLength_CS1->Fill(path, Charge );
// if(NHighStrip==2) Charge_Vs_PathLength_CS2->Fill(path, Charge );
// if(NHighStrip==3) Charge_Vs_PathLength_CS3->Fill(path, Charge );
// if(NHighStrip==4) Charge_Vs_PathLength_CS4->Fill(path, Charge );
// if(NHighStrip==5) Charge_Vs_PathLength_CS5->Fill(path, Charge );
HFirstStrip ->Fill(FirstStrip);
if(Validation){ClusterChargeOverPath=ClusterChargeOverPath/APV->PreviousGain;}
APV_Charge ->Fill(APV->Index,ClusterChargeOverPath);
APV_Momentum ->Fill(APV->Index,trajState.globalMomentum().mag());
APV_PathLength->Fill(APV->Index,path);
return ClusterChargeOverPath;
}
| SiStripApvGain * SiStripGainFromData::getNewObject | ( | ) | [private, virtual] |
Implements ConditionDBWriter< SiStripApvGain >.
Definition at line 1284 of file SiStripGainFromData.cc.
References a, AlgoMode, sistrip::APV, stAPVGain::APVId, APVsCollOrdered, gather_cfg::cout, stAPVGain::DetId, stAPVGain::Gain, NULL, VarParsing::obj, stAPVGain::PreviousGain, and SiStripApvGain::put().
{
cout << "START getNewObject\n";
// if( !(strcmp(AlgoMode.c_str(),"WriteOnDB")==0 || strcmp(AlgoMode.c_str(),"SingleJob")==0) )return NULL;
if( !(strcmp(AlgoMode.c_str(),"WriteOnDB")==0 || strcmp(AlgoMode.c_str(),"SingleJob")==0) )return new SiStripApvGain();
SiStripApvGain * obj = new SiStripApvGain();
std::vector<float>* theSiStripVector = NULL;
int PreviousDetId = -1;
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;
}
printf("%i | %i | PreviousGain = %7.5f NewGain = %7.5f\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain);
theSiStripVector->push_back(APV->Gain);
// 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);
}
cout << "END getNewObject\n";
return obj;
}
| void SiStripGainFromData::getPeakOfLandau | ( | TH1 * | InputHisto, |
| double * | FitResults, | ||
| double | LowRange = 0, |
||
| double | HighRange = 5400 |
||
| ) | [private] |
Definition at line 1232 of file SiStripGainFromData.cc.
References MaxChi2OverNDF, MinNrEntries, and tablePrinter::width.
Referenced by algoEndJob().
{
double adcs = -0.5;
double adcs_err = 0.;
double width = -0.5;
double width_err = 0;
double chi2overndf = -0.5;
double nr_of_entries = InputHisto->GetEntries();
if( (unsigned int)nr_of_entries < MinNrEntries){
FitResults[0] = adcs;
FitResults[1] = adcs_err;
FitResults[2] = width;
FitResults[3] = width_err;
FitResults[4] = chi2overndf;
return;
}
// perform fit with standard landau
TF1* MyLandau = new TF1("MyLandau","landau",LowRange, HighRange);
MyLandau->SetParameter("MPV",300);
InputHisto->Fit("MyLandau","QR WW");
TF1 * fitfunction = (TF1*) InputHisto->GetListOfFunctions()->First();
// MPV is parameter 1 (0=constant, 1=MPV, 2=Sigma)
adcs = fitfunction->GetParameter("MPV");
adcs_err = fitfunction->GetParError(1);
width = fitfunction->GetParameter(2);
width_err = fitfunction->GetParError(2);
chi2overndf = fitfunction->GetChisquare() / fitfunction->GetNDF();
// if still wrong, give up
if(adcs<2. || chi2overndf>MaxChi2OverNDF){
adcs = -0.5; adcs_err = 0.;
width = -0.5; width_err = 0;
chi2overndf = -0.5;
}
FitResults[0] = adcs;
FitResults[1] = adcs_err;
FitResults[2] = width;
FitResults[3] = width_err;
FitResults[4] = chi2overndf;
delete MyLandau;
}
| bool SiStripGainFromData::IsFarFromBorder | ( | TrajectoryStateOnSurface | trajState, |
| const uint32_t | detid, | ||
| const edm::EventSetup * | iSetup | ||
| ) | [private] |
Definition at line 1194 of file SiStripGainFromData.cc.
References BoundSurface::bounds(), gather_cfg::cout, edm::EventSetup::get(), Bounds::length(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), Parameters::parameters, LocalTrajectoryError::positionError(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().
Referenced by ComputeChargeOverPath().
{
edm::ESHandle<TrackerGeometry> tkGeom; iSetup->get<TrackerDigiGeometryRecord>().get( tkGeom );
LocalPoint HitLocalPos = trajState.localPosition();
LocalError HitLocalError = trajState.localError().positionError() ;
const GeomDetUnit* it = tkGeom->idToDetUnit(DetId(detid));
if (dynamic_cast<const StripGeomDetUnit*>(it)==0 && dynamic_cast<const PixelGeomDetUnit*>(it)==0) {
std::cout << "this detID doesn't seem to belong to the Tracker" << std::endl;
return false;
}
const BoundPlane plane = it->surface();
const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));
double DistFromBorder = 1.0;
//double HalfWidth = it->surface().bounds().width() /2.0;
double HalfLength = it->surface().bounds().length() /2.0;
if(trapezoidalBounds)
{
std::vector<float> const & parameters = (*trapezoidalBounds).parameters();
HalfLength = parameters[3];
//double t = (HalfLength + HitLocalPos.y()) / (2*HalfLength) ;
//HalfWidth = parameters[0] + (parameters[1]-parameters[0]) * t;
}else if(rectangularBounds){
//HalfWidth = it->surface().bounds().width() /2.0;
HalfLength = it->surface().bounds().length() /2.0;
}else{return false;}
// if (fabs(HitLocalPos.x())+HitLocalError.xx() >= (HalfWidth - DistFromBorder) ) return false;//Don't think is really necessary
if (fabs(HitLocalPos.y())+HitLocalError.yy() >= (HalfLength - DistFromBorder) ) return false;
return true;
}
std::string SiStripGainFromData::AlgoMode [private] |
Definition at line 127 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), algoBeginRun(), algoEndJob(), getNewObject(), and SiStripGainFromData().
bool SiStripGainFromData::AllowSaturation [private] |
Definition at line 121 of file SiStripGainFromData.cc.
Referenced by ComputeChargeOverPath(), and SiStripGainFromData().
TH2F* SiStripGainFromData::APV_Charge [private] |
Definition at line 198 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH1F* SiStripGainFromData::APV_CumulGain [private] |
Definition at line 203 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::APV_DetId [private] |
Definition at line 192 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
TH1F* SiStripGainFromData::APV_Eta [private] |
Definition at line 194 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
TH1F* SiStripGainFromData::APV_Gain [private] |
Definition at line 202 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::APV_Id [private] |
Definition at line 193 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
TH2F* SiStripGainFromData::APV_Momentum [private] |
Definition at line 197 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH1F* SiStripGainFromData::APV_MPV [private] |
Definition at line 201 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::APV_PathLength [private] |
Definition at line 199 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH1F* SiStripGainFromData::APV_PathLengthM [private] |
Definition at line 200 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::APV_PrevGain [private] |
Definition at line 204 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoBeginRun().
TH1F* SiStripGainFromData::APV_R [private] |
Definition at line 195 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
TH1F* SiStripGainFromData::APV_SubDet [private] |
Definition at line 196 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
TH1F* SiStripGainFromData::APV_Thickness [private] |
Definition at line 205 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
__gnu_cxx::hash_map<unsigned int, stAPVGain*, __gnu_cxx::hash<unsigned int>, isEqual > SiStripGainFromData::APVsColl [private] |
Definition at line 289 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
std::vector<stAPVGain*> SiStripGainFromData::APVsCollOrdered [private] |
Definition at line 288 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoBeginRun(), algoEndJob(), and getNewObject().
int SiStripGainFromData::CalibrationLevel [private] |
Definition at line 124 of file SiStripGainFromData.cc.
Referenced by algoEndJob(), and SiStripGainFromData().
TH1F* SiStripGainFromData::Charge_TEC [private] |
Definition at line 179 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TEC1 [private] |
Definition at line 180 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TEC2 [private] |
Definition at line 181 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TECM [private] |
Definition at line 183 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TECP [private] |
Definition at line 182 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TIB [private] |
Definition at line 174 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TID [private] |
Definition at line 175 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TIDM [private] |
Definition at line 177 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TIDP [private] |
Definition at line 176 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::Charge_TOB [private] |
Definition at line 178 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::Charge_Vs_Alpha [private] |
Definition at line 252 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_Beta [private] |
Definition at line 256 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathLength [private] |
Definition at line 235 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathLength320 [private] |
Definition at line 236 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathLength500 [private] |
Definition at line 237 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathTEC [private] |
Definition at line 162 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathTEC1 [private] |
Definition at line 163 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathTEC2 [private] |
Definition at line 164 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathTIB [private] |
Definition at line 159 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathTID [private] |
Definition at line 160 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_PathTOB [private] |
Definition at line 161 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::Charge_Vs_TransversAngle [private] |
Definition at line 248 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
bool SiStripGainFromData::CheckIfFileExist [private] |
Definition at line 125 of file SiStripGainFromData.cc.
Referenced by algoEndJob(), and SiStripGainFromData().
bool SiStripGainFromData::CheckLocalAngle [private] |
Definition at line 110 of file SiStripGainFromData.cc.
Referenced by SiStripGainFromData().
DQMStore* SiStripGainFromData::dqmStore_ [private] |
Definition at line 99 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and SiStripGainFromData().
DQMStore* SiStripGainFromData::dqmStore_infile [private] |
Definition at line 100 of file SiStripGainFromData.cc.
Referenced by SiStripGainFromData().
unsigned int SiStripGainFromData::EEvent [private] |
Definition at line 279 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::Error_Vs_Entries [private] |
Definition at line 261 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::Error_Vs_Eta [private] |
Definition at line 262 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::Error_Vs_MPV [private] |
Definition at line 260 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::Error_Vs_Phi [private] |
Definition at line 263 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
unsigned int SiStripGainFromData::ERun [private] |
Definition at line 278 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
TimeValue_t SiStripGainFromData::ETimestamp [private] |
Definition at line 280 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and algoBeginJob().
bool SiStripGainFromData::FirstSetOfConstants [private] |
Definition at line 122 of file SiStripGainFromData.cc.
Referenced by algoBeginRun(), algoEndJob(), and SiStripGainFromData().
TH1F* SiStripGainFromData::FWHM_Vs_PathLength [private] |
Definition at line 243 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::FWHM_Vs_PathLength320 [private] |
Definition at line 244 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::FWHM_Vs_PathLength500 [private] |
Definition at line 245 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::HChi2OverNDF [private] |
Definition at line 141 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::HFirstStrip [private] |
Definition at line 272 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::HitLocalPosition [private] |
Definition at line 267 of file SiStripGainFromData.cc.
TH2F* SiStripGainFromData::HitLocalPositionBefCut [private] |
Definition at line 268 of file SiStripGainFromData.cc.
TH1F* SiStripGainFromData::HTrackChi2OverNDF [private] |
Definition at line 142 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::HTrackHits [private] |
Definition at line 143 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
const edm::Event* SiStripGainFromData::iEvent_ [private] |
Definition at line 108 of file SiStripGainFromData.cc.
Referenced by algoAnalyze().
const edm::EventSetup* SiStripGainFromData::iSetup_ [private] |
Definition at line 107 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
TH1F* SiStripGainFromData::JobInfo [private] |
Definition at line 270 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
double SiStripGainFromData::MaxChi2OverNDF [private] |
Definition at line 113 of file SiStripGainFromData.cc.
Referenced by getPeakOfLandau(), and SiStripGainFromData().
double SiStripGainFromData::MaxMPVError [private] |
Definition at line 112 of file SiStripGainFromData.cc.
Referenced by algoEndJob(), and SiStripGainFromData().
unsigned int SiStripGainFromData::MaxNrStrips [private] |
Definition at line 118 of file SiStripGainFromData.cc.
Referenced by ComputeChargeOverPath(), and SiStripGainFromData().
double SiStripGainFromData::MaxTrackChiOverNdf [private] |
Definition at line 120 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
double SiStripGainFromData::MaxTrackEta [private] |
Definition at line 117 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
double SiStripGainFromData::MaxTrackMomentum [private] |
Definition at line 115 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
unsigned int SiStripGainFromData::MinNrEntries [private] |
Definition at line 111 of file SiStripGainFromData.cc.
Referenced by getPeakOfLandau(), and SiStripGainFromData().
double SiStripGainFromData::MinTrackEta [private] |
Definition at line 116 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
unsigned int SiStripGainFromData::MinTrackHits [private] |
Definition at line 119 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
double SiStripGainFromData::MinTrackMomentum [private] |
Definition at line 114 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
TH1F* SiStripGainFromData::MPV_Vs_Alpha [private] |
Definition at line 253 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_Beta [private] |
Definition at line 257 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_Eta [private] |
Definition at line 186 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_EtaTEC [private] |
Definition at line 148 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_EtaTEC1 [private] |
Definition at line 149 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_EtaTEC2 [private] |
Definition at line 150 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_EtaTIB [private] |
Definition at line 145 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_EtaTID [private] |
Definition at line 146 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_EtaTOB [private] |
Definition at line 147 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathLength [private] |
Definition at line 239 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathLength320 [private] |
Definition at line 240 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathLength500 [private] |
Definition at line 241 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathTEC [private] |
Definition at line 169 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathTEC1 [private] |
Definition at line 170 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathTEC2 [private] |
Definition at line 171 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathTIB [private] |
Definition at line 166 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathTID [private] |
Definition at line 167 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_PathTOB [private] |
Definition at line 168 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_Phi [private] |
Definition at line 185 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_PhiTEC [private] |
Definition at line 155 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_PhiTEC1 [private] |
Definition at line 156 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_PhiTEC2 [private] |
Definition at line 157 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_PhiTIB [private] |
Definition at line 152 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_PhiTID [private] |
Definition at line 153 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_PhiTOB [private] |
Definition at line 154 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::MPV_Vs_R [private] |
Definition at line 187 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPV_Vs_TransversAngle [private] |
Definition at line 249 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPVs [private] |
Definition at line 207 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPVs320 [private] |
Definition at line 208 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::MPVs500 [private] |
Definition at line 209 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
unsigned int SiStripGainFromData::NEvent [private] |
Definition at line 274 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::NHighStripInCluster [private] |
Definition at line 215 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::NoMPV_Vs_EtaPhi [private] |
Definition at line 265 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
TH1F* SiStripGainFromData::NSatStripInCluster [private] |
Definition at line 214 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::NStrips_Vs_Alpha [private] |
Definition at line 254 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::NStrips_Vs_Beta [private] |
Definition at line 258 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and ComputeChargeOverPath().
TH2F* SiStripGainFromData::NStrips_Vs_TransversAngle [private] |
Definition at line 250 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), algoEndJob(), and ComputeChargeOverPath().
TH1F* SiStripGainFromData::NumberOfEntriesByAPV [private] |
Definition at line 140 of file SiStripGainFromData.cc.
Referenced by algoBeginJob(), and algoEndJob().
std::string SiStripGainFromData::OutputGains [private] |
Definition at line 128 of file SiStripGainFromData.cc.
Referenced by algoEndJob(), and SiStripGainFromData().
std::string SiStripGainFromData::OutputHistos [private] |
Definition at line 129 of file SiStripGainFromData.cc.
Referenced by algoEndJob(), and SiStripGainFromData().
unsigned int SiStripGainFromData::SEvent [private] |
Definition at line 276 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
unsigned int SiStripGainFromData::SRun [private] |
Definition at line 275 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
TimeValue_t SiStripGainFromData::STimestamp [private] |
Definition at line 277 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and algoBeginJob().
MonitorElement* SiStripGainFromData::tmp [private] |
Definition at line 135 of file SiStripGainFromData.cc.
Referenced by algoBeginJob().
TH2F* SiStripGainFromData::Tracks_P_Vs_Eta [private] |
Definition at line 137 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
TH2F* SiStripGainFromData::Tracks_Pt_Vs_Eta [private] |
Definition at line 138 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), algoBeginJob(), and algoEndJob().
std::string SiStripGainFromData::TrajToTrackLabel [private] |
Definition at line 131 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
std::string SiStripGainFromData::TrajToTrackProducer [private] |
Definition at line 130 of file SiStripGainFromData.cc.
Referenced by algoAnalyze(), and SiStripGainFromData().
bool SiStripGainFromData::Validation [private] |
Definition at line 123 of file SiStripGainFromData.cc.
Referenced by algoBeginRun(), ComputeChargeOverPath(), and SiStripGainFromData().
vector<string> SiStripGainFromData::VInputFiles [private] |
Definition at line 133 of file SiStripGainFromData.cc.
Referenced by algoEndJob(), and SiStripGainFromData().
1.7.3