Inheritance diagram for dune::XYZcalibration:

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

virtual	~XYZcalibration ()

void	beginJob ()

void	endJob ()

void	beginRun (const art::Run &run)

void	analyze (const art::Event &evt)

void	reset ()

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 (SharedResources const &resources)

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 &)

fhicl::ParameterSetID	selectorConfig () 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
TTree *	fEventTree

Int_t	run

Int_t	subrun

Int_t	event

Double_t	evttime

Int_t	year_month_date

Int_t	hour_min_sec

Int_t	cross_trks

Int_t	stopping_trks

Int_t	all_trks

Float_t	xprojectedlen [kMaxTracks]

Float_t	trackthetaxz [kMaxTracks]

Float_t	trackthetayz [kMaxTracks]

Float_t	trkstartx [kMaxTracks]

Float_t	trkstarty [kMaxTracks]

Float_t	trkstartz [kMaxTracks]

Float_t	trkendx [kMaxTracks]

Float_t	trkendy [kMaxTracks]

Float_t	trkendz [kMaxTracks]

Float_t	trklen [kMaxTracks]

Int_t	TrkID [kMaxTracks]

Float_t	trkstartcosxyz [kMaxTracks][3]

Float_t	trkendcosxyz [kMaxTracks][3]

Int_t	ntrkhits [kMaxTracks][3]

Float_t	trkdqdx [kMaxTracks][3][3000]

Float_t	trkdedx [kMaxTracks][3][3000]

Float_t	trkresrange [kMaxTracks][3][3000]

Float_t	trkhitx [kMaxTracks][3][3000]

Float_t	trkhity [kMaxTracks][3][3000]

Float_t	trkhitz [kMaxTracks][3][3000]

std::string	fHitsModuleLabel

std::string	fTrackModuleLabel

std::string	fCalorimetryModuleLabel

bool	fSaveCaloInfo

bool	fSaveTrackInfo

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

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

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

	Observer (fhicl::ParameterSet const &config)

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

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 77 of file XYZcalibration_module.cc.

Constructor & Destructor Documentation

dune::XYZcalibration::XYZcalibration ( fhicl::ParameterSet const & pset )

explicit

Definition at line 128 of file XYZcalibration_module.cc.

                                                               :
     EDAnalyzer(pset),
     fHitsModuleLabel          (pset.get< std::string >("HitsModuleLabel","")         ),
     fTrackModuleLabel         (pset.get< std::string >("TrackModuleLabel","")        ),
     fCalorimetryModuleLabel   (pset.get< std::string >("CalorimetryModuleLabel","")  ),
     fSaveCaloInfo             (pset.get< bool>("SaveCaloInfo",false)),
     fSaveTrackInfo            (pset.get< bool>("SaveTrackInfo",false))
   {
     if (fSaveTrackInfo == false) fSaveCaloInfo = false;
   }

dune::XYZcalibration::~XYZcalibration ( )

virtual

Definition at line 140 of file XYZcalibration_module.cc.

140 {

141 }

Member Function Documentation

void dune::XYZcalibration::analyze ( const art::Event & evt )

Definition at line 197 of file XYZcalibration_module.cc.

                                                   {
     reset();  
      
    
     std::vector<art::Ptr<recob::Track> > tracklist;
     std::vector<art::Ptr<recob::PFParticle> > pfplist;
   
     auto trackListHandle = evt.getHandle< std::vector<recob::Track> >(fTrackModuleLabel);
     if (trackListHandle) art::fill_ptr_vector(tracklist, trackListHandle);
 
     auto PFPListHandle = evt.getHandle< std::vector<recob::PFParticle> >("pandora");
     if (PFPListHandle) art::fill_ptr_vector(pfplist, PFPListHandle);
   
     art::FindManyP<anab::Calorimetry> fmcal(trackListHandle, evt, fCalorimetryModuleLabel);
     art::FindManyP<anab::T0> trk_t0_assn_v(PFPListHandle, evt ,"pandora");
     art::FindManyP<recob::PFParticle> pfp_trk_assn(trackListHandle,evt,"pandoraTrack");
     art::FindManyP<anab::T0> fmT0(trackListHandle, evt ,"pmtrack");
 
     run = evt.run();
     subrun = evt.subRun();
     event = evt.id().event();
     art::Timestamp ts = evt.time();
     TTimeStamp tts(ts.timeHigh(), ts.timeLow());
     evttime=tts.AsDouble();
      
     UInt_t year=0;
     UInt_t month=0;
     UInt_t day=0;
      
     year_month_date=tts.GetDate(kTRUE,0,&year,&month,&day);
      
     UInt_t hour=0;
     UInt_t min=0;
     UInt_t sec=0;
      
     hour_min_sec=tts.GetTime(kTRUE,0,&hour,&min,&sec);
   
     cross_trks=0;
     stopping_trks=0;
     all_trks=0;
   
     size_t NTracks = tracklist.size();
     for(size_t i=0; i<NTracks;++i){
       art::Ptr<recob::Track> ptrack(trackListHandle, i);
       if(fTrackModuleLabel=="pandoraTrack"){
         std::vector<art::Ptr<recob::PFParticle>> pfps=pfp_trk_assn.at(i);
         if(!pfps.size()) continue;
         std::vector<art::Ptr<anab::T0>> t0s=trk_t0_assn_v.at(pfps[0].key());
         if(!t0s.size()) continue;
         //auto t0 = t0s.at(0);
         // double t_zero=t0->Time();
       }
       if(fTrackModuleLabel=="pmtrack"){
         std::vector<art::Ptr<anab::T0>> T0s=fmT0.at(i);
         if(T0s.size()==0)
           continue;
       }
       all_trks++;
       std::vector<art::Ptr<anab::Calorimetry>> calos=fmcal.at(i);     
 
       const recob::Track& track = *ptrack;
       auto pos = track.Vertex();
       auto dir_start = track.VertexDirection();
       auto dir_end   = track.EndDirection();
       auto end = track.End();
       double theta_xz = std::atan2(dir_start.X(), dir_start.Z());
       double theta_yz = std::atan2(dir_start.Y(), dir_start.Z());
       //deltaX->Fill(X);
       cross_trks++;
       if(std::abs(pos.X())<340 && std::abs(end.X())<340 &&  pos.Y()>20 && pos.Y()<580 && end.Y()>20 && end.Y()<580 &&  pos.Z()>20 && pos.Z()<675 && end.Z()>20 && end.Z()<675) stopping_trks++;//this infact gives both ends confined, we have not counted the stopping muon tracks in this module
       trackthetaxz[cross_trks-1]=theta_xz;
       trackthetayz[cross_trks-1]=theta_yz;
       trkstartx[cross_trks-1]=pos.X();
       trkstarty[cross_trks-1]=pos.Y();
       trkstartz[cross_trks-1]=pos.Z();
       trkendx[cross_trks-1]=end.X();
       trkendy[cross_trks-1]=end.Y();
       trkendz[cross_trks-1]=end.Z();
       trklen[cross_trks-1]=track.Length();
       TrkID[cross_trks-1]=track.ID();
       trkstartcosxyz[cross_trks-1][0]=dir_start.X();
       trkstartcosxyz[cross_trks-1][1]=dir_start.Y();
       trkstartcosxyz[cross_trks-1][2]=dir_start.Z();
       trkendcosxyz[cross_trks-1][0]=dir_end.X();
       trkendcosxyz[cross_trks-1][1]=dir_end.Y();
       trkendcosxyz[cross_trks-1][2]=dir_end.Z();
                                    
       for(size_t ical = 0; ical<calos.size(); ++ical){
         if(!calos[ical]) continue;
         if(!calos[ical]->PlaneID().isValid) continue;
         int planenum = calos[ical]->PlaneID().Plane;
         if(planenum<0||planenum>2) continue;
         const size_t NHits = calos[ical] -> dEdx().size();
         ntrkhits[cross_trks-1][planenum]=int(NHits);
         for(size_t iHit = 0; iHit < NHits; ++iHit){
           const auto& TrkPos = (calos[ical] -> XYZ())[iHit];
           trkdqdx[cross_trks-1][planenum][iHit]=(calos[ical] -> dQdx())[iHit];
           trkdedx[cross_trks-1][planenum][iHit]=(calos[ical] -> dEdx())[iHit];
           trkresrange[cross_trks-1][planenum][iHit]=(calos[ical]->ResidualRange())[iHit];
           trkhitx[cross_trks-1][planenum][iHit]=TrkPos.X();
           trkhity[cross_trks-1][planenum][iHit]=TrkPos.Y();
           trkhitz[cross_trks-1][planenum][iHit]=TrkPos.Z();
         } // loop over iHit..
       } // loop over ical 2 nd time...
     } // loop over trks...
     fEventTree->Fill();
   } // end of analyze function

void dune::XYZcalibration::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 145 of file XYZcalibration_module.cc.

                                {
     std::cout<<"job begin..."<<std::endl;
     art::ServiceHandle<art::TFileService> tfs;
     // deltaX = tfs->make<TH1D>("deltaX","Plot of deltaX",400,0,400);
     // deltaX->Sumw2();
     fEventTree = tfs->make<TTree>("Event", "Event Tree from Reco");
     fEventTree->Branch("event", &event,"event/I");
     fEventTree->Branch("evttime",&evttime,"evttime/D");
     fEventTree->Branch("run", &run,"run/I");
     fEventTree->Branch("subrun", &subrun,"surbrun/I");
     fEventTree->Branch("year_month_date", &year_month_date,"year_month_date/I");
     fEventTree->Branch("hour_min_sec", &hour_min_sec,"hour_min_sec/I");
     fEventTree->Branch("cross_trks",&cross_trks,"cross_trks/I");
     fEventTree->Branch("stopping_trks",&stopping_trks,"stopping_trks/I");
     fEventTree->Branch("all_trks",&all_trks,"all_trks/I");
     fEventTree->Branch("xprojectedlen",xprojectedlen,"xprojectedlen[all_trks]/F");
     fEventTree->Branch("trackthetaxz",trackthetaxz,"trackthetaxz[cross_trks]/F");
     fEventTree->Branch("trackthetayz",trackthetayz,"trackthetayz[cross_trks]/F");
     fEventTree->Branch("trkstartx",trkstartx,"trkstartx[cross_trks]/F");
     fEventTree->Branch("trkstarty",trkstarty,"trkstarty[cross_trks]/F");
     fEventTree->Branch("trkstartz",trkstartz,"trkstartz[cross_trks]/F");
     fEventTree->Branch("trkendx",trkendx,"trkendx[cross_trks]/F");
     fEventTree->Branch("trkendy",trkendy,"trkendy[cross_trks]/F");
     fEventTree->Branch("trkendz",trkendz,"trkendz[cross_trks]/F");
     fEventTree->Branch("trklen",trklen,"trklen[cross_trks]/F");
     fEventTree->Branch("TrkID",TrkID,"TrkID[cross_trks]/I");
     fEventTree->Branch("trkstartcosxyz",trkstartcosxyz,"trkstartcosxyz[cross_trks][3]/F");
     fEventTree->Branch("trkendcosxyz",trkendcosxyz,"trkendcosxyz[cross_trks][3]/F");
     fEventTree->Branch("ntrkhits",ntrkhits,"ntrkhits[cross_trks][3]/I");
     fEventTree->Branch("trkdqdx",trkdqdx,"trkdqdx[cross_trks][3][3000]/F");
     fEventTree->Branch("trkdedx",trkdedx,"trkdedx[cross_trks][3][3000]/F");
     fEventTree->Branch("trkresrange",trkresrange,"trkresrange[cross_trks][3][3000]/F");
     fEventTree->Branch("trkhitx",trkhitx,"trkhitx[cross_trks][3][3000]/F");
     fEventTree->Branch("trkhity",trkhity,"trkhity[cross_trks][3][3000]/F");
     fEventTree->Branch("trkhitz",trkhitz,"trkhitz[cross_trks][3][3000]/F");
   }

void dune::XYZcalibration::beginRun ( const art::Run & run )

Definition at line 188 of file XYZcalibration_module.cc.

                                             {
     mf::LogInfo("XYZcalibration")<<"begin run..."<<std::endl;
   }

void dune::XYZcalibration::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 183 of file XYZcalibration_module.cc.

                              {     
 
   }

void dune::XYZcalibration::reset ( )

Definition at line 306 of file XYZcalibration_module.cc.

                             {
     run = -9999;
     subrun = -9999;
     event = -9999;
     evttime = -9999;
     cross_trks = -9999;
     all_trks = -9999;
     year_month_date=-9999;
     hour_min_sec=-9999;
     for(int i=0; i<kMaxTracks; i++){
       trackthetaxz[i]=-9999;
       trackthetayz[i]=-9999;
       trkstartx[i]=-9999;
       trkstarty[i]=-9999;
       trkstartz[i]=-9999;
       trkendx[i]=-9999;
       trkendy[i]=-9999;
       trkendz[i]=-9999;
       trklen[i]=-9999;
       TrkID[i]=-9999;
       xprojectedlen[i]=-9999;
       trkstartcosxyz[i][0]=-9999;
       trkstartcosxyz[i][1]=-9999;
       trkstartcosxyz[i][2]=-9999; 
       trkendcosxyz[i][0]=-9999;
       trkendcosxyz[i][1]=-9999;
       trkendcosxyz[i][2]=-9999;
       ntrkhits[i][0] = -9999;
       ntrkhits[i][1] = -9999;
       ntrkhits[i][2] = -9999;
       for(int j=0; j<3; j++){
         for(int k=0; k<3000; k++){
           trkdqdx[i][j][k]=-9999;
           trkdedx[i][j][k]=-9999;
           trkresrange[i][j][k]=-9999;
           trkhitx[i][j][k]=-9999;
           trkhity[i][j][k]=-9999;
           trkhitz[i][j][k]=-9999;
         }
       }
     }
   }

Member Data Documentation

Int_t dune::XYZcalibration::all_trks

private

Definition at line 99 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::cross_trks

private

Definition at line 97 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::event

private

Definition at line 93 of file XYZcalibration_module.cc.

Double_t dune::XYZcalibration::evttime

private

Definition at line 94 of file XYZcalibration_module.cc.

std::string dune::XYZcalibration::fCalorimetryModuleLabel

private

Definition at line 122 of file XYZcalibration_module.cc.

TTree* dune::XYZcalibration::fEventTree

private

Definition at line 90 of file XYZcalibration_module.cc.

std::string dune::XYZcalibration::fHitsModuleLabel

private

Definition at line 120 of file XYZcalibration_module.cc.

bool dune::XYZcalibration::fSaveCaloInfo

private

Definition at line 123 of file XYZcalibration_module.cc.

bool dune::XYZcalibration::fSaveTrackInfo

private

Definition at line 124 of file XYZcalibration_module.cc.

std::string dune::XYZcalibration::fTrackModuleLabel

private

Definition at line 121 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::hour_min_sec

private

Definition at line 96 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::ntrkhits[kMaxTracks][3]

private

Definition at line 113 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::run

private

Definition at line 91 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::stopping_trks

private

Definition at line 98 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::subrun

private

Definition at line 92 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trackthetaxz[kMaxTracks]

private

Definition at line 101 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trackthetayz[kMaxTracks]

private

Definition at line 102 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkdedx[kMaxTracks][3][3000]

private

Definition at line 115 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkdqdx[kMaxTracks][3][3000]

private

Definition at line 114 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkendcosxyz[kMaxTracks][3]

private

Definition at line 112 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkendx[kMaxTracks]

private

Definition at line 106 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkendy[kMaxTracks]

private

Definition at line 107 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkendz[kMaxTracks]

private

Definition at line 108 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkhitx[kMaxTracks][3][3000]

private

Definition at line 117 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkhity[kMaxTracks][3][3000]

private

Definition at line 118 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkhitz[kMaxTracks][3][3000]

private

Definition at line 119 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::TrkID[kMaxTracks]

private

Definition at line 110 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trklen[kMaxTracks]

private

Definition at line 109 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkresrange[kMaxTracks][3][3000]

private

Definition at line 116 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkstartcosxyz[kMaxTracks][3]

private

Definition at line 111 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkstartx[kMaxTracks]

private

Definition at line 103 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkstarty[kMaxTracks]

private

Definition at line 104 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::trkstartz[kMaxTracks]

private

Definition at line 105 of file XYZcalibration_module.cc.

Float_t dune::XYZcalibration::xprojectedlen[kMaxTracks]

private

Definition at line 100 of file XYZcalibration_module.cc.

Int_t dune::XYZcalibration::year_month_date

private

Definition at line 95 of file XYZcalibration_module.cc.

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

protoduneana/protoduneana/singlephase/XYZcalibration/XYZcalibration_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation