forked from JeffersonLab/hd_utilities
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
contains * plugin FCALbadchannels.cc, FCALbadchannels.hh
- Loading branch information
1 parent
c6fcc3e
commit 58733f2
Showing
2 changed files
with
413 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,316 @@ | ||
| // $Id$ | ||
| // | ||
| // File: FCALbadchannels.cc | ||
| // Created: Wed Sep 11 13:39:58 EDT 2019 | ||
| // Creator: cakondi (on Linux ifarm1802 3.10.0-327.el7.x86_64 x86_64) | ||
| // | ||
| #include <iostream> | ||
| #include <iomanip> | ||
| #include <map> | ||
| #include "FCALbadchannels.hh" | ||
| using namespace jana; | ||
| #include <JANA/JEvent.h> | ||
| #include <JANA/JApplication.h> | ||
| #include <JANA/JFactory.h> | ||
| //// DFCAL Hit factory.cc included the following header files | ||
| #include <HDGEOMETRY/DGeometry.h> | ||
| #include "FCAL/DFCALDigiHit.h" | ||
| #include "FCAL/DFCALGeometry.h" | ||
| #include "FCAL/DFCALHit_factory.h" | ||
| #include "DAQ/Df250PulseIntegral.h" | ||
| #include "DAQ/Df250PulsePedestal.h" | ||
| #include "DAQ/Df250Config.h" | ||
| #include "TTAB/DTTabUtilities.h" | ||
| #include <FCAL/DFCALHit.h> | ||
| #include <TRIGGER/DL1Trigger.h> | ||
| #include <DAQ/DL1Info.h> | ||
| #include "DVector2.h" | ||
| #include "units.h" | ||
| #include <DVector3.h> | ||
| #include "FCAL/DFCALCluster.h" | ||
| #include "FCAL/DFCALCluster_factory.h" | ||
| #include "HDGEOMETRY/DGeometry.h" | ||
| #include <FCAL/DFCALShower.h> | ||
| #include <FCAL/DFCALHit.h> | ||
| #include <PID/DEventRFBunch.h> | ||
| #include "TApplication.h" | ||
| #include "TH3.h" | ||
|
|
||
|
|
||
|
|
||
|
|
||
|
|
||
|
|
||
| // Routine used to create our JEventProcessor | ||
| #include <JANA/JApplication.h> | ||
| #include <JANA/JFactory.h> | ||
| extern "C"{ | ||
| void InitPlugin(JApplication *app){ | ||
| InitJANAPlugin(app); | ||
| app->AddProcessor(new FCALbadchannels()); | ||
| } | ||
| } // "C" | ||
|
|
||
|
|
||
| //------------------ | ||
| // FCALbadchannels (Constructor) | ||
| //------------------ | ||
| FCALbadchannels::FCALbadchannels() | ||
| { | ||
|
|
||
| } | ||
|
|
||
| //------------------ | ||
| // ~FCALbadchannels (Destructor) | ||
| //------------------ | ||
| FCALbadchannels::~FCALbadchannels() | ||
| { | ||
|
|
||
| } | ||
|
|
||
| //------------------ | ||
| // init | ||
| //------------------ | ||
| jerror_t FCALbadchannels::init(void) | ||
| { | ||
| char hname[100]; | ||
| const int max_blocks = 2800, max_rows = 59 , max_cols = 59; | ||
| static const int nranges = 5; | ||
|
|
||
| htrigfp = new TH1F("htrigfp", "", 32,0,31); | ||
| hrawtrigbit = new TH1F("hrawtrigbit", "", 32,0,31); | ||
|
|
||
| hAllInte = new TH1F("hAllInte", "", 2000,0,200000); | ||
|
|
||
| hTotalInte = new TH1F("hTotalInte", "", 20000,0.,200000000.); | ||
| hTotalInteCut = new TH1F("hTotalInteCut", "", 20000,0.,200000000.); | ||
| hRanges = new TH1F("hRanges", "", 20000,0.,200000000.); | ||
|
|
||
| hhitcnter = new TH1F("hhitcnter", "", 11,0,10); | ||
| hTotalInteHits = new TH2F("hTotalInteHits", "",20000,0.,200000000., 11,0,10); | ||
|
|
||
| hNhits = new TH1F("hNhits", "", max_blocks,-0.5,max_blocks-0.5); | ||
| hEff = new TH1F("hEff", "", max_blocks,-0.5,max_blocks-0.5); | ||
|
|
||
| hRowColOcc = new TH2F("hRowColOcc", "", max_cols,-0.5,max_cols-0.5, max_rows,-0.5,max_rows-0.5); | ||
| hRowColEff = new TH2F("hRowColEff", "", max_cols,-0.5,max_cols-0.5, max_rows,-0.5,max_rows-0.5); | ||
|
|
||
| for(int ira=0; ira<nranges; ira++){ | ||
| sprintf(hname,"hEffReg_%d", ira); hEffReg_[ira] = new TH1F(hname,"",max_blocks,-0.5,max_blocks-0.5); | ||
| sprintf(hname,"hRowColReg_%d", ira); hRowColReg_[ira] = new TH2F(hname,"",max_cols,-0.5,max_cols-0.5, max_rows,-0.5,max_rows-0.5); | ||
| } | ||
|
|
||
| //----------------- | ||
| //--- for LED time: | ||
| // | ||
| hTotalInteE = new TH1F("hTotalInteE", "", 2000,0.,20000.); | ||
| hHitTime2121= new TH1F("hHitTime2121","",211,-10.5,200.5); | ||
| hTotalInteE_Cut = new TH1F("hTotalInteE_Cut", "", 2000,0.,20000.); | ||
| hRangesE = new TH1F("hRangesE", "", 2000,0.,20000.); | ||
| for(Int_t ibl = 0; ibl < max_blocks; ibl++) { | ||
| sprintf(hname,"hHitTimeRange0_%d",ibl); | ||
| hHitTimeRange0_[ibl]= new TH1F(hname,"",211,-10.5,200.5); | ||
| sprintf(hname,"hHitTimeRange1_%d",ibl); | ||
| hHitTimeRange1_[ibl]= new TH1F(hname,"",211,-10.5,200.5); | ||
| sprintf(hname,"hHitTimeRange2_%d",ibl); | ||
| hHitTimeRange2_[ibl]= new TH1F(hname,"",211,-10.5,200.5); | ||
| sprintf(hname,"hHitTimeRange3_%d",ibl); | ||
| hHitTimeRange3_[ibl]= new TH1F(hname,"",211,-10.5,200.5); | ||
| sprintf(hname,"hHitTimeRange4_%d",ibl); | ||
| hHitTimeRange4_[ibl]= new TH1F(hname,"",211,-10.5,200.5); | ||
| } | ||
|
|
||
|
|
||
| return NOERROR; | ||
| } | ||
|
|
||
| //------------------ | ||
| // brun | ||
| //------------------ | ||
| jerror_t FCALbadchannels::brun(JEventLoop *eventLoop, int32_t runnumber) | ||
| { | ||
| // This is called whenever the run number changes | ||
| return NOERROR; | ||
| } | ||
|
|
||
| //------------------ | ||
| // evnt | ||
| //------------------ | ||
| jerror_t FCALbadchannels::evnt(JEventLoop *loop, uint64_t eventnumber) | ||
| { | ||
| ///// extract the FCAL Geometry information DFCALHits_factory.cc | ||
|
|
||
| vector<const DFCALGeometry*> fcalGeomVect; | ||
| loop->Get( fcalGeomVect ); | ||
| if (fcalGeomVect.size() < 1) | ||
| return OBJECT_NOT_AVAILABLE; | ||
| const DFCALGeometry& fcalGeom = *(fcalGeomVect[0]); | ||
|
|
||
|
|
||
| ///////////////Trigger and Front Panel Trigger ///////////////////////// | ||
|
|
||
| vector<const DL1Trigger *> l1trig; | ||
| loop->Get(l1trig); | ||
|
|
||
| japp->RootFillLock(this); | ||
|
|
||
| int rawtrignumber,trigfpnumber; | ||
| if( l1trig.size() > 0) { | ||
| for(unsigned int trigcounter = 0; trigcounter < 32; trigcounter++){ | ||
| trig_bit[trigcounter+1] = (l1trig[0]->trig_mask & (1 << trigcounter)) ? 1 : 0; | ||
| if (trig_bit[trigcounter+1] == 1) { | ||
| trigbit=trigcounter+1; | ||
| rawtrignumber=trigbit; | ||
| //hrawtrigbit->Fill(trigbit); | ||
| } | ||
| } | ||
| hrawtrigbit->Fill(trigbit); | ||
|
|
||
| for(unsigned int trigfpcounter = 0; trigfpcounter < 32; trigfpcounter++){ | ||
| trig_bit_fp[trigfpcounter+1] = (l1trig[0]->fp_trig_mask & (1 << trigfpcounter)) ? 1 : 0; | ||
| if (trig_bit_fp[trigfpcounter+1] == 1) { | ||
| trigfpbit=trigfpcounter+1; | ||
| trigfpnumber=trigfpbit; | ||
| //htrigfp->Fill(trigfpbit); | ||
| } | ||
| } | ||
| htrigfp->Fill(trigfpbit); | ||
| } | ||
|
|
||
| //----------------- | ||
| if(trigfpnumber!=3) return NOERROR; // trigfpnumber!=3 for FCAL FP trigger | ||
| // if(rawtrignumber!=1) return NOERROR; // trignumber!=1 for rawdata // trignumber !4 for PS | ||
|
|
||
| vector<const DFCALDigiHit*> digihits; | ||
| loop->Get(digihits); | ||
|
|
||
| Double_t RegLo[5] = { 5e+6, 10e+6, 15e+6, 28e+6, 33e+6}; | ||
| // Double_t RegHi[5] = {10e+6, 15e+6, 20e+6, 33e+6, 38e+6}; | ||
| Double_t RegHi[5] = {10e+6, 15e+6, 20e+6, 33e+6, 45e+6};// RunPeriod-2019-11 -- not great | ||
|
|
||
| for(int ira=0; ira<nranges; ira++){ | ||
| hRanges->Fill(RegLo[ira]); | ||
| hRanges->Fill(RegHi[ira]); | ||
| } | ||
|
|
||
| float TotalInte=0; | ||
| int Row[max_blocks], Column[max_blocks]; | ||
| int hitcnter[max_blocks]; | ||
| for(int ibl=0; ibl<max_blocks; ibl++){ | ||
| hitcnter[ibl] = 0; | ||
| } | ||
|
|
||
| for(unsigned int i=0; i < digihits.size(); i++) { | ||
| const DFCALDigiHit *digihit = digihits[i]; | ||
| if (digihit->pulse_time == 0 || digihit->pedestal == 0 || digihit->pulse_peak == 0) continue; | ||
| int ibl = fcalGeom.channel(digihit->row,digihit->column); | ||
| Row[ibl] = digihit->row; | ||
| Column[ibl] = digihit->column; | ||
| hitcnter[ibl]++; | ||
|
|
||
| hNhits->Fill(ibl); | ||
| hRowColOcc->Fill(Column[ibl],Row[ibl]); | ||
|
|
||
| double inte = (digihit->pulse_integral) -(((double)digihit->pedestal/digihit->nsamples_pedestal)*digihit->nsamples_integral); | ||
|
|
||
| hAllInte->Fill(inte); | ||
|
|
||
| TotalInte += inte; | ||
| } | ||
| hTotalInte->Fill(TotalInte); | ||
|
|
||
| if (TotalInte > RegLo[0] && TotalInte <= RegHi[4]){ | ||
| hTotalInteCut->Fill(TotalInte); | ||
|
|
||
| for(int ibl=0; ibl<max_blocks; ibl++){ | ||
| if (hitcnter[ibl] > 0) { | ||
|
|
||
| hEff->Fill(ibl); | ||
| hRowColEff->Fill(Column[ibl],Row[ibl]); | ||
|
|
||
| for(int ira=0; ira<nranges; ira++){ | ||
| if (TotalInte > RegLo[ira] && TotalInte <= RegHi[ira]){ | ||
| hEffReg_[ira]->Fill(ibl); | ||
| hRowColReg_[ira]->Fill(Column[ibl],Row[ibl]); | ||
| } | ||
| } | ||
|
|
||
| int yval; | ||
| if (hitcnter[ibl] <= 10) yval = hitcnter[ibl]; | ||
| if (hitcnter[ibl] > 10) yval = 10; | ||
| hhitcnter->Fill(yval); | ||
| hTotalInteHits->Fill(TotalInte, yval); | ||
|
|
||
| } // hitcnter | ||
| } // ibl | ||
| } // TotalInte threshold | ||
|
|
||
|
|
||
| //----------------- | ||
| //--- for LED time: | ||
| // | ||
| vector<const DFCALHit *> fcalhits; | ||
| loop->Get(fcalhits); | ||
|
|
||
| // RunPeriod-2019-11 | ||
| Double_t RegLoE[5] = {2000, 3200, 4700, 6500, 8000}; | ||
| Double_t RegHiE[5] = {3200, 4700, 6500, 8000, 10000}; | ||
|
|
||
| for(int ira=0; ira<nranges; ira++){ | ||
| hRangesE->Fill(RegLoE[ira]); | ||
| hRangesE->Fill(RegHiE[ira]); | ||
| } | ||
|
|
||
| float TotalInteE=0; | ||
| for (unsigned int i = 0; i < fcalhits.size(); i++){ | ||
| const DFCALHit *fcalhit = fcalhits[i]; | ||
| int ibl = fcalGeom.channel(fcalhit->row, fcalhit->column); | ||
| if(trigfpnumber==3){ | ||
| TotalInteE += fcalhit->E; | ||
| if (ibl == 2121) hHitTime2121->Fill(fcalhit->t); | ||
| } | ||
| } | ||
|
|
||
| hTotalInteE->Fill(TotalInteE); | ||
| for(int ira=0; ira<nranges; ira++){ | ||
| if (TotalInteE > RegLoE[ira] && TotalInteE <= RegHiE[ira]){ | ||
| hTotalInteE_Cut->Fill(TotalInteE); | ||
| for (unsigned int i = 0; i < fcalhits.size(); i++){ | ||
| const DFCALHit *fcalhit = fcalhits[i]; | ||
| int ibl = fcalGeom.channel(fcalhit->row, fcalhit->column); | ||
| if (ira == 0) hHitTimeRange0_[ibl]->Fill(fcalhit->t); | ||
| if (ira == 1) hHitTimeRange1_[ibl]->Fill(fcalhit->t); | ||
| if (ira == 2) hHitTimeRange2_[ibl]->Fill(fcalhit->t); | ||
| if (ira == 3) hHitTimeRange3_[ibl]->Fill(fcalhit->t); | ||
| if (ira == 4) hHitTimeRange4_[ibl]->Fill(fcalhit->t); | ||
| } | ||
| } | ||
| } | ||
|
|
||
|
|
||
|
|
||
| japp->RootFillUnLock(this); | ||
|
|
||
| return NOERROR; | ||
| } | ||
|
|
||
| //------------------ | ||
| // erun | ||
| //------------------ | ||
| jerror_t FCALbadchannels::erun(void) | ||
| { | ||
| // This is called whenever the run number changes, before it is | ||
| // changed to give you a chance to clean up before processing | ||
| // events from the next run number. | ||
| return NOERROR; | ||
| } | ||
|
|
||
| //------------------ | ||
| // fini | ||
| //------------------ | ||
| jerror_t FCALbadchannels::fini(void) | ||
| { | ||
| // Called before program exit after event processing is finished. | ||
| return NOERROR; | ||
| } | ||
|
|
Oops, something went wrong.