68 #include "TMultiGraph.h"
85 virtual void endJob()
override ;
141 h1BadAPVperEvent_->SetXTitle(
"# Modules with Bad APVs");
142 h1BadAPVperEvent_->SetYTitle(
"Entries");
143 h1BadAPVperEvent_->SetLineWidth(2);
144 h1BadAPVperEvent_->SetLineStyle(2);
147 h1APVCM_->SetXTitle(
"APV CM [adc]");
182 std::vector<uint32_t> detIdV;
185 for(uint32_t
i=0;
i < detIdV.size(); ++
i){
188 pedestals.resize((pedestalsRange.second- pedestalsRange.first)*8/10);
200 for (; itCMDetSetV != moduleCM->end(); ++itCMDetSetV){
202 for(;itCM != itCMDetSetV->
end(); ++itCM)
h1APVCM_->Fill(itCM->adc());
241 uint32_t NBabAPVs = moduleRawDigi->size();
242 std::cout<<
"Number of module with HIP in this event: " << NBabAPVs << std::endl;
245 for (; itRawDigis != moduleRawDigi->end(); ++itRawDigis) {
247 uint32_t detId = itRawDigis->id;
251 if(itDSBaseline->id != detId){
252 std::cout <<
"Collections out of Synch. Something of fishy is going on ;-)" << std::endl;
266 bool restAPV[6] = {0,0,0,0,0,0};
267 int strip =0, totADC=0;
268 int minAPVRes = 7, maxAPVRes = -1;
269 for(;itRaw != itRawDigis->
end(); ++itRaw, ++strip){
270 float adc = itRaw->adc();
278 if(APV>maxAPVRes) maxAPVRes =
APV;
279 if(APV<minAPVRes) minAPVRes =
APV;
285 float minx = -0.5, maxx=767.5;
287 minx = minAPVRes * 128 -0.5;
288 maxx = maxAPVRes * 128 + 127.5;
292 sprintf(detIds,
"%ul", detId);
293 sprintf(evs,
"%llu",
event);
294 sprintf(
runs,
"%u", run);
295 char* dHistoName = Form(
"Id:%s_run:%s_ev:%s",detIds,
runs, evs);
296 h1ProcessedRawDigis_ = sdProcessedRawDigis_.make<TH1F>(dHistoName,dHistoName, bins, minx, maxx);
299 h1Baseline_ = sdBaseline_.make<TH1F>(dHistoName,dHistoName, bins, minx, maxx);
310 h1Clusters_ = sdClusters_.make<TH1F>(dHistoName,dHistoName, bins, minx, maxx);
327 std::vector<int16_t> ProcessedRawDigis(itRawDigis->size());
332 std::vector<int16_t>::const_iterator itProcessedRawDigis;
335 for(itProcessedRawDigis = ProcessedRawDigis.begin();itProcessedRawDigis != ProcessedRawDigis.end(); ++itProcessedRawDigis){
336 if(restAPV[strip/128]){
337 float adc = *itProcessedRawDigis;
350 for ( ; itClusters != clusters->
end(); ++itClusters ){
352 if(itClusters->
id() == detId){
353 int firststrip = clus->firstStrip();
356 for(
auto itAmpl = clus->amplitudes().begin(); itAmpl != clus->amplitudes().end(); ++itAmpl){
int adc(sample_type sample)
get the ADC sample (12 bits)
T getParameter(std::string const &) const
EventNumber_t event() const
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
const_iterator end(bool update=false) const
virtual void endJob() override
#define DEFINE_FWK_MODULE(type)
unsigned long long EventNumber_t
std::vector< TH1F > vBaselineHisto_
T * make(const Args &...args) const
make new ROOT object
data_type const * const_iterator
std::pair< ContainerIterator, ContainerIterator > Range
edm::InputTag srcProcessedRawDigi_
virtual void beginJob() override
id_type id(size_t cell) const
virtual void analyze(const edm::Event &, const edm::EventSetup &) override
~SiStripBaselineAnalyzer()
SiStripBaselineAnalyzer(const edm::ParameterSet &)
std::vector< TH1F > vBaselinePointsHisto_
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
edm::Service< TFileService > fs_
std::vector< TH1F > vClusterHisto_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
iterator end()
Return the off-the-end iterator.
TH1F * h1ProcessedRawDigis_
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
edm::InputTag srcBaseline_
std::vector< int > pedestals
std::vector< TH1F > vProcessedRawDigiHisto_
edm::InputTag srcBaselinePoints_
uint16_t nModuletoDisplay_
iterator begin()
Return an iterator to the first DetSet.
edm::ESHandle< SiStripPedestals > pedestalsHandle
collection_type::const_iterator const_iterator
collection_type::const_iterator const_iterator
std::auto_ptr< SiStripPedestalsSubtractor > subtractorPed_
const_iterator begin(bool update=false) const
static std::auto_ptr< SiStripPedestalsSubtractor > create_SubtractorPed(const edm::ParameterSet &)