Inheritance diagram for AnalysisExample::RawEVD:

Public Member Functions
	RawEVD (fhicl::ParameterSet const &pset)

virtual	~RawEVD ()

void	beginJob ()

void	beginRun (const art::Run &run)

void	reconfigure (fhicl::ParameterSet const &pset)

void	analyze (const art::Event &evt)

unsigned int	getAPAindex (unsigned int apa)

unsigned int	getTPCindex (unsigned int tpc)

Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDAnalyzer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept

	Analyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	Analyzer (Table< Config > const &config)

void	doBeginJob ()

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

std::string const &	processName () const

bool	wantAllEvents () const

bool	wantEvent (Event const &e)

fhicl::ParameterSetID	selectorConfig () const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
std::string	fRawDigitLabel

unsigned int	fNticks

unsigned int	fUChanMin

unsigned int	fUChanMax

unsigned int	fVChanMin

unsigned int	fVChanMax

unsigned int	fZ0ChanMin

unsigned int	fZ0ChanMax

unsigned int	fZ1ChanMin

unsigned int	fZ1ChanMax

unsigned int	fNofAPA

unsigned int	fChansPerAPA

art::ServiceHandle< geo::Geometry >	fGeom

std::vector< TH2I * >	fTimeChanU

std::vector< TH2I * >	fTimeChanV

std::vector< TH2I * >	fTimeChanZ0

std::vector< TH2I * >	fTimeChanZ1

std::vector< TH2I * >	fTimeChanThumbU

std::vector< TH2I * >	fTimeChanThumbV

std::vector< TH2I * >	fTimeChanThumbZ0

std::vector< TH2I * >	fTimeChanThumbZ1

TH2I *	fChargeSumU

TH2I *	fChargeSumV

TH2I *	fChargeSumZ

unsigned int	UVPlane [120]

unsigned int	ZPlane [240]

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &paths, fhicl::ParameterSet const &config)

detail::ProcessAndEventSelectors &	processAndEventSelectors ()

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 57 of file RawEVD_module.cc.

Constructor & Destructor Documentation

AnalysisExample::RawEVD::RawEVD ( fhicl::ParameterSet const & pset )

explicit

Definition at line 175 of file RawEVD_module.cc.

         : EDAnalyzer(parameterSet)
         {
     this->reconfigure(parameterSet);
   }

AnalysisExample::RawEVD::~RawEVD ( )

virtual

Definition at line 194 of file RawEVD_module.cc.

194 {

195 }

Member Function Documentation

void AnalysisExample::RawEVD::analyze ( const art::Event & evt )

Definition at line 359 of file RawEVD_module.cc.

                                              {
 
         unsigned int tpcid, cryoid;
         std::stringstream  thumbnameZ0, thumbnameZ1;
 
     // get the objects holding all of the raw data information
     art::Handle< std::vector<raw::RawDigit> > Raw;
     event.getByLabel(fRawDigitLabel, Raw);
 
     // put it in a more easily usable form
     std::vector< art::Ptr<raw::RawDigit> >  Digits;
     art::fill_ptr_vector(Digits, Raw);
 
     //loop through all RawDigits (over entire channels)
     for(size_t d = 0; d < Digits.size(); d++){
     art::Ptr<raw::RawDigit> digit;
          digit=Digits.at(d);
 
     // get the channel number for this digit
     uint32_t chan = digit->Channel();
     unsigned int apa = std::floor( chan/fChansPerAPA );
          tpcid=fGeom->ChannelToWire(chan)[0].TPC;
          cryoid=fGeom->ChannelToWire(chan)[0].Cryostat;
 
     std::vector<short> uncompressed(digit->Samples());
     raw::Uncompress(digit->ADCs(), uncompressed, digit->Compression());
 
         if( fGeom->View(chan) == geo::kU ){
                 for(unsigned int l=0;l<uncompressed.size();l++) {
                         if(uncompressed.at(l)!=0){
                                 fTimeChanU[apa]->Fill(chan,l, uncompressed.at(l));
                                 if(uncompressed.at(l)>0) fTimeChanThumbU[apa]->Fill(chan,l, uncompressed.at(l));
                                 fChargeSumU->Fill(1+getAPAindex(apa)%6, 20-(getAPAindex(apa)/6),std::abs(uncompressed.at(l))/2);
                                 }
                         }
                 }
         if( fGeom->View(chan) == geo::kV ){
                 for(unsigned int l=0;l<uncompressed.size();l++) {
                         if(uncompressed.at(l)!=0){
                                 fTimeChanV[apa]->Fill(chan,l, uncompressed.at(l));
                                 if(uncompressed.at(l)>0) fTimeChanThumbV[apa]->Fill(chan,l, uncompressed.at(l));
                                 fChargeSumV->Fill(1+getAPAindex(apa)%6, 20-(getAPAindex(apa)/6),std::abs(uncompressed.at(l))/2);
                                 }
                         }
                 }
         if ( fGeom->View(chan) == geo::kZ && fGeom->ChannelToWire(chan)[0].TPC % 2 == 0 ){
                 for(unsigned int l=0;l<uncompressed.size();l++) {
                         if(uncompressed.at(l)!=0){
                                 fTimeChanZ0[apa]->Fill(chan,l, uncompressed.at(l));
                                 if(uncompressed.at(l)>0) fTimeChanThumbZ0[apa]->Fill(chan,l, uncompressed.at(l));
                                 fChargeSumZ->Fill(1+getTPCindex(cryoid*120+tpcid)%12, 20-(getTPCindex(cryoid*120+tpcid)/12),uncompressed.at(l));
                                 }
                         }
                 }
         if ( fGeom->View(chan) == geo::kZ && fGeom->ChannelToWire(chan)[0].TPC % 2 == 1 ){
                 for(unsigned int l=0;l<uncompressed.size();l++) {
                         if(uncompressed.at(l)!=0){
                                 fTimeChanZ1[apa]->Fill(chan,l, uncompressed.at(l));
                                 if(uncompressed.at(l)>0) fTimeChanThumbZ1[apa]->Fill(chan,l, uncompressed.at(l));
                                 fChargeSumZ->Fill(1+getTPCindex(cryoid*120+tpcid)%12, 20-(getTPCindex(cryoid*120+tpcid)/12),uncompressed.at(l));
                                 }
                         }
                 }
 
     } // end RawDigit loop
 
     return;
   }

void AnalysisExample::RawEVD::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 199 of file RawEVD_module.cc.

                        {
     // Access ART's TFileService, which will handle creating and writing
     // histograms and n-tuples for us. 
     art::ServiceHandle<art::TFileService> tfs;
 
         //Histogram names and titles
         std::stringstream  name, title;
 
     unsigned int UChMin;
     unsigned int UChMax;
     unsigned int VChMin;
     unsigned int VChMax;
     unsigned int Z0ChMin;
     unsigned int Z0ChMax;
     unsigned int Z1ChMin;
     unsigned int Z1ChMax;
     TH2I* TempHisto;
 
     std::ofstream outfile;
         outfile.open("msglog.txt");
 
     // loop through channels in the first APA to find the channel boundaries for each view
     // will adjust for desired APA after
     fUChanMin = 0;
          fNofAPA=fGeom->NTPC()*fGeom->Ncryostats()/2; //NTPC:120 for each Cryostat
     fChansPerAPA = fGeom->Nchannels()/fNofAPA; 
     fZ1ChanMax = fChansPerAPA - 1;
     for ( unsigned int c = fUChanMin + 1; c < fZ1ChanMax; c++ ){
       if ( fGeom->View(c) == geo::kV && fGeom->View(c-1) == geo::kU ){
         fVChanMin = c;
         fUChanMax = c - 1;
       }
       if ( fGeom->View(c) == geo::kZ && fGeom->View(c-1) == geo::kV ){
         fZ0ChanMin = c;
         fVChanMax = c-1;
       }
       if ( fGeom->View(c) == geo::kZ && fGeom->ChannelToWire(c)[0].TPC == fGeom->ChannelToWire(c-1)[0].TPC + 1 ){
         fZ1ChanMin = c;
         fZ0ChanMax = c-1;
       }
     }
 
 outfile<<fChansPerAPA<<"  "<<fGeom->Ncryostats()<<"  "<<fNofAPA<<std::endl;
 
     unsigned int minT = 0;
     unsigned int maxT = 0;
     minT = 0;
     maxT = fNticks;
     unsigned int binT = (maxT-minT);
 
         for(unsigned int i=0;i<fNofAPA;i++){
     UChMin=fUChanMin + i*fChansPerAPA;
     UChMax=fUChanMax + i*fChansPerAPA;
     VChMin=fVChanMin + i*fChansPerAPA;
     VChMax=fVChanMax + i*fChansPerAPA;
     Z0ChMin=fZ0ChanMin + i*fChansPerAPA;
     Z0ChMax=fZ0ChanMax + i*fChansPerAPA;
     Z1ChMin=fZ1ChanMin + i*fChansPerAPA;
     Z1ChMax=fZ1ChanMax + i*fChansPerAPA;
 
     // construct the histograms; TH2 constructors: ("Name", "Title", NxBin, xMin, xMax, NyBin, yMax, yMin)
     name.str("");
     name << "fTimeChanU";
     name <<  i;
     title.str("");
     title << "Time vs Channel(Plane U, APA";
     title << i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), UChMax - UChMin + 1, UChMin, UChMax, binT, minT, maxT);
          fTimeChanU.push_back(TempHisto);
 
     name.str("");
     name << "fTimeChanThumbU";
     name <<  i;
     title.str("");
     title << "Time vs Channel(Plane U, APA";
     title << i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), 32, UChMin, UChMax, 32, minT, maxT);
          fTimeChanThumbU.push_back(TempHisto);
 
     name.str("");
     name << "fTimeChanV";
     name << i;
     title.str("");
     title << "Time vs Channel(Plane V, APA";
     title << i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), VChMax - VChMin + 1, VChMin, VChMax, binT, minT, maxT);
          fTimeChanV.push_back(TempHisto);
 
     name.str("");
     name << "fTimeChanThumbV";
     name << i;
     title.str("");
     title << "Time vs Channel(Plane V, APA";
     title << i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), 32, VChMin, VChMax, 32, minT, maxT);
          fTimeChanThumbV.push_back(TempHisto);
 
     name.str("");
     name << "fTimeChanZ0";
     name << i;
     title.str("");
     title << "Time vs Channel(Plane Z0, APA";
     title <<i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), Z0ChMax - Z0ChMin + 1, Z0ChMin, Z0ChMax, binT, minT, maxT);
          fTimeChanZ0.push_back(TempHisto);
 
     name.str("");
     name << "fTimeChanThumbZ0";
     name << i;
     title.str("");
     title << "Time vs Channel(Plane Z0, APA";
     title <<i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), 32, Z0ChMin, Z0ChMax, 32, minT, maxT);
          fTimeChanThumbZ0.push_back(TempHisto);
 
     name.str("");
     name << "fTimeChanZ1";
     name << i;
     title.str("");
     title << "Time vs Channel(Plane Z1, APA";
     title << i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), Z1ChMax - Z1ChMin + 1, Z1ChMin, Z1ChMax, binT, minT, maxT);
          fTimeChanZ1.push_back(TempHisto);
 
     name.str("");
     name << "fTimeChanThumbZ1";
     name << i;
     title.str("");
     title << "Time vs Channel(Plane Z1, APA";
     title << i<<")";
     TempHisto = tfs->make<TH2I>(name.str().c_str(),title.str().c_str(), 32, Z1ChMin, Z1ChMax, 32, minT, maxT);
          fTimeChanThumbZ1.push_back(TempHisto);
 
     fTimeChanU[i]->SetStats(0);    fTimeChanV[i]->SetStats(0);    
          fTimeChanZ0[i]->SetStats(0);    fTimeChanZ1[i]->SetStats(0);
     fTimeChanThumbU[i]->SetStats(0);    fTimeChanThumbV[i]->SetStats(0);    
          fTimeChanThumbZ0[i]->SetStats(0);    fTimeChanThumbZ1[i]->SetStats(0);
 
     fTimeChanU[i]->GetXaxis()->SetTitle("Channel"); fTimeChanU[i]->GetYaxis()->SetTitle("TDC");
     fTimeChanV[i]->GetXaxis()->SetTitle("Channel"); fTimeChanV[i]->GetYaxis()->SetTitle("TDC");
     fTimeChanZ0[i]->GetXaxis()->SetTitle("Channel"); fTimeChanZ0[i]->GetYaxis()->SetTitle("TDC");
     fTimeChanZ1[i]->GetXaxis()->SetTitle("Channel"); fTimeChanZ1[i]->GetYaxis()->SetTitle("TDC");
         }
 
     fChargeSumU= tfs->make<TH2I>("hChargeSumU","Charge Sum on U-planes", 6,0.5,6.5,20,0.5,20.5 );
     fChargeSumV= tfs->make<TH2I>("hChargeSumV","Charge Sum on V-planes", 6,0.5,6.5,20,0.5,20.5 );
     fChargeSumZ= tfs->make<TH2I>("hChargeSumZ","Charge Sum on Z-planes", 12,0.5,12.5,20,0.5,20.5);
     fChargeSumU->SetStats(0);    fChargeSumV->SetStats(0);    fChargeSumZ->SetStats(0);  
 
   }

void AnalysisExample::RawEVD::beginRun ( const art::Run & run )

Definition at line 352 of file RawEVD_module.cc.

                                       {
 
   }

unsigned int AnalysisExample::RawEVD::getAPAindex ( unsigned int apa )

inline

Definition at line 71 of file RawEVD_module.cc.

                                                    {
          for (int k=0; k<120; k++){if(apa==UVPlane[k]) {return k; }}
          return -1;
         }

unsigned int AnalysisExample::RawEVD::getTPCindex ( unsigned int tpc )

inline

Definition at line 76 of file RawEVD_module.cc.

                                                    {
          for (int k=0; k<240; k++){if(tpc==ZPlane[k]) {return k; }}
          return -1;
         }

void AnalysisExample::RawEVD::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 184 of file RawEVD_module.cc.

                                                     {
     fRawDigitLabel  =  p.get< std::string >("RawDigitLabel");
     auto const *fDetProp = lar::providerFrom<detinfo::DetectorPropertiesService>();
     fNticks         = fDetProp->NumberTimeSamples();
     return;
   }

Member Data Documentation

unsigned int AnalysisExample::RawEVD::fChansPerAPA

private

Definition at line 101 of file RawEVD_module.cc.

TH2I* AnalysisExample::RawEVD::fChargeSumU

private

Definition at line 121 of file RawEVD_module.cc.

TH2I* AnalysisExample::RawEVD::fChargeSumV

private

Definition at line 122 of file RawEVD_module.cc.

TH2I* AnalysisExample::RawEVD::fChargeSumZ

private

Definition at line 123 of file RawEVD_module.cc.

art::ServiceHandle<geo::Geometry> AnalysisExample::RawEVD::fGeom

private

Definition at line 109 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fNofAPA

private

Definition at line 100 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fNticks

private

Definition at line 88 of file RawEVD_module.cc.

std::string AnalysisExample::RawEVD::fRawDigitLabel

private

Definition at line 84 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanThumbU

private

Definition at line 116 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanThumbV

private

Definition at line 117 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanThumbZ0

private

Definition at line 118 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanThumbZ1

private

Definition at line 119 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanU

private

Definition at line 111 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanV

private

Definition at line 112 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanZ0

private

Definition at line 113 of file RawEVD_module.cc.

std::vector<TH2I*> AnalysisExample::RawEVD::fTimeChanZ1

private

Definition at line 114 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fUChanMax

private

Definition at line 92 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fUChanMin

private

Definition at line 91 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fVChanMax

private

Definition at line 94 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fVChanMin

private

Definition at line 93 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fZ0ChanMax

private

Definition at line 96 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fZ0ChanMin

private

Definition at line 95 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fZ1ChanMax

private

Definition at line 98 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::fZ1ChanMin

private

Definition at line 97 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::UVPlane[120]

private

Initial value:

={
        119,118,117,59,58,57,
        113,112,111,53,52,51,
        107,106,105,47,46,45,
        101,100,99,41,40,39,
        95,94,93,35,34,33,
        89,88,87,29,28,27,
        83,82,81,23,22,21,
        77,76,75,17,16,15,
        71,70,69,11,10,9,
        65,64,63,5,4,3,
        116,115,114,56,55,54,
        110,109,108,50,49,48,
        104,103,102,44,43,42,
        98,97,96,38,37,36,
        92,91,90,32,31,30,
        86,85,84,26,25,24,
        80,79,78,20,19,18,
        74,73,72,14,13,12,
        68,67,66,8,7,6,
        62,61,60,2,1,0}

Definition at line 125 of file RawEVD_module.cc.

unsigned int AnalysisExample::RawEVD::ZPlane[240]

private

Initial value:

={
        239,238,237,236,235,234,119,118,117,116,115,114,
        227,226,225,224,223,222,107,106,105,104,103,102,
        215,214,213,212,211,210,95,94,93,92,91,90,
        203,202,201,200,199,198,83,82,81,80,79,78,
        191,190,189,188,187,186,71,70,69,68,67,66,
        179,178,177,176,175,174,59,58,57,56,55,54,
        167,166,165,164,163,162,47,46,45,44,43,42,
        155,154,153,152,151,150,35,34,33,32,31,30,
        143,142,141,140,139,138,23,22,21,20,19,18,
        131,130,129,128,127,126,11,10,9,8,7,6,
        233,232,231,230,229,228,113,112,111,110,109,108,
        221,220,219,218,217,216,101,100,99,98,97,96,
        209,208,207,206,205,204,89,88,87,86,85,84,
        197,196,195,194,193,192,77,76,75,74,73,72,
        185,184,183,182,181,180,65,64,63,62,61,60,
        173,172,171,170,169,168,53,52,51,50,49,48,
        161,160,159,158,157,156,41,40,39,38,37,36,
        149,148,147,146,145,144,29,28,27,26,25,24,
        137,136,135,134,133,132,17,16,15,14,13,12,
        125,124,123,122,121,120,5,4,3,2,1,0
        }

Definition at line 147 of file RawEVD_module.cc.

The documentation for this class was generated from the following file:

dunetpc/dune/RawdataDisplay/RawEVD_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation