9 #ifndef ProtoDUNE_opdet_eventdisplay_H 10 #define ProtoDUNE_opdet_eventdisplay_H 15 #include "TLorentzVector.h" 57 #include "art_root_io/TFileService.h" 58 #include "art_root_io/TFileDirectory.h" 61 #include "canvas/Persistency/Common/FindManyP.h" 149 #endif // ProtoDUNE_opdet_eventdisplay_H 170 fTimeEnd = clockData.OpticalClock().FramePeriod();
171 fSampleFreq = clockData.OpticalClock().Frequency();
199 float displace[288]={0};
201 for(
int j=0;j<287;j++) displace[j]=j;
202 for(
int j=0;j<284;j+=4){
207 if((j>=132 && j<=143) || (j>=264 && j<=275)){
225 Float_t binsz[] = {0,10, 29, 48, 66, 84, 102, 120, 138, 156, 174, 192, 211, 230, 240, 257, 278, 296, 314, 332, 350, 368, 386, 404, 422, 441, 460, 470, 489, 508, 526, 544, 562, 580, 598, 616, 634, 650, 671, 690, 700};
226 Int_t binnumz =
sizeof(binsz)/
sizeof(Float_t) - 1;
228 Float_t binsy[]={-600,-540, -480, -420, -360, -300, -240, -180, -120, -60, 0, 60, 120, 180, 240, 300, 360, 420, 480, 540, 600};
229 Int_t binnumy =
sizeof(binsy)/
sizeof(Float_t) - 1;
249 std::vector<TH2D*> DispHist;
253 sprintf(HistName,
"Event %d Disps", evt.
id().
event());
259 for(
unsigned int i = 0; i < FlashHandle->size(); ++i)
270 TH2D * ThisEventDisp =
nullptr;
271 ThisEventDisp = tfs->make<TH2D>(HistName,
"PE from OpHits in OpFlash;z ;DaS RaS (y)",
279 for(
size_t i=0; i!=OpHitHandle->size(); ++i)
282 geom->OpDetGeoFromOpChannel(
fOpChannel).GetCenter(xyz);
287 fPeakTime = OpHitHandle->at(i).PeakTime();
288 fPE = OpHitHandle->at(i).PE();
310 TLine *
line =
new TLine(0,0,600,0);
311 line->SetLineColor(kBlack);
std::vector< int > fFlashFrameVector
std::vector< int > fFlashIDVector
std::vector< float > fAbsTimeVector
Handle< PROD > getHandle(SelectorBase const &) const
std::vector< bool > fInBeamFrameVector
std::string fOpFlashModuleLabel
std::string fOpHitModuleLabel
std::vector< float > fTotalPEVector
EDAnalyzer(fhicl::ParameterSet const &pset)
std::vector< float > fPEsPerFlashPerChannelVector
art framework interface to geometry description
std::vector< float > fYWidthVector
std::vector< float > fZCenterVector
std::vector< float > fZWidthVector
#define DEFINE_ART_MODULE(klass)
T get(std::string const &key) const
std::vector< int > fOnBeamTimeVector
virtual ~ProtoDUNE_opdet_eventdisplay()
std::vector< float > fFlashTimeVector
Encapsulate the geometry of an optical detector.
void analyze(const art::Event &)
void line(double t, double *p, double &x, double &y, double &z)
ProtoDUNE_opdet_eventdisplay(const fhicl::ParameterSet &)
EventNumber_t event() const
Access the description of detector geometry.
std::vector< float > fYCenterVector
pure virtual base interface for detector clocks