Inheritance diagram for dune::michelremoving:

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

virtual	~michelremoving ()

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

Int_t	unbroken_trks

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]

Float_t	trkpitch [kMaxTracks][3][3000]

Float_t	peakT_max [kMaxTracks]

Float_t	peakT_min [kMaxTracks]

Float_t	dist_min [kMaxTracks]

Int_t	adjacent_hits [kMaxTracks]

Int_t	lastwire [kMaxTracks]

Int_t	endtpc [kMaxTracks]

Float_t	lastpeakt [kMaxTracks]

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 84 of file michelremoving_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 143 of file michelremoving_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::michelremoving::~michelremoving ( )

virtual

Definition at line 155 of file michelremoving_module.cc.

155 {

156 }

Member Function Documentation

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

Definition at line 219 of file michelremoving_module.cc.

                                                   {
     reset();  
     std::vector<art::Ptr<recob::Track> > tracklist;
     auto trackListHandle = evt.getHandle< std::vector<recob::Track> >("pandoraTrack");
     if (trackListHandle){
       art::fill_ptr_vector(tracklist, trackListHandle);
     }
     else return;
 
     std::vector<art::Ptr<recob::PFParticle> > pfplist;    
     auto PFPListHandle = evt.getHandle< std::vector<recob::PFParticle> >("pandora"); 
     if (PFPListHandle) art::fill_ptr_vector(pfplist, PFPListHandle);
   
     /******new lines*************************/
     std::vector<art::Ptr<recob::Hit>> hitlist;
     auto hitListHandle = evt.getHandle< std::vector<recob::Hit> >(fHitsModuleLabel) ; // to get information about the hits
     if (hitListHandle) art::fill_ptr_vector(hitlist, hitListHandle);
        
     art::FindManyP<recob::Hit, recob::TrackHitMeta> fmthm(trackListHandle, evt, fTrackModuleLabel); // to associate tracks and hits
 
     art::FindManyP<recob::Hit> fmtht(trackListHandle, evt, fTrackModuleLabel); // to associate tracks and hits
     // art::FindManyP<recob::Track, recob::TrackHitMeta> thass(hitListHandle, evt, fTrackModuleLabel); //to associate hit 
     art::FindManyP<recob::Track> thass(hitListHandle, evt, fTrackModuleLabel); //to associate hit just trying
     
     art::FindManyP<anab::Calorimetry> fmcal(trackListHandle, evt, fCalorimetryModuleLabel);
     if (!fmcal.isValid()){
       std::cout<<"art::Assns<recob::Track,anab::Calorimetry,void> with module label: "<<fCalorimetryModuleLabel<<" not found."<<std::endl;
       return;
     }
 
     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;
     unbroken_trks=0;
     // size_t NHits=hitlist.size();
     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;
       //TVector3 pos, dir_start, dir_end, end;
       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());  
       float startx=pos.X();   
       float starty=pos.Y();
       float startz=pos.Z();
       float endx=end.X();
       float endy=end.Y();
       float endz=end.Z();
 
       float startcosx=dir_start.X();
       float startcosy=dir_start.Y();
       float startcosz=dir_start.Z();
       float endcosx=dir_end.X();
       float endcosy=dir_end.Y();
       float endcosz=dir_end.Z();
       float tracklength=track.Length();
       size_t count = 0;
       /********************************************************************************************************/
 
       /*********************storing hits for each tracks*******************************/
 
       std::vector<art::Ptr<recob::Hit>> allHits=fmtht.at(i); //storing hits for ith track
       std::vector<float> hitpeakT;
       for(size_t h1=0;h1<allHits.size();h1++){
         hitpeakT.push_back(allHits[h1]->PeakTime());
         //      if(allHits[h1]->WireID().Plane==2 && allHits[h1]->PeakTime()>420 && allHits[h1]->PeakTime()<530 && allHits[h1]->WireID().TPC==9 && allHits[h1]->WireID().Wire>60 && allHits[h1]->WireID().Wire<80) cout<<"wire id and peak time from allHits this is from track hit association "<<allHits[h1]->WireID().Wire<<"  "<<allHits[h1]->PeakTime()<<endl;
       }
       float max=-999999;
       float min=-999999;
       min=*min_element(hitpeakT.begin(),hitpeakT.end());
       max=*max_element(hitpeakT.begin(),hitpeakT.end());
       hitpeakT.clear();
 
         float dist0=99999;
         float dist=99999;
         float peaktime=-1;
         int wireno=99999;
         int counter1=99999;
         int tpcno=-1;
 
       if(((std::abs(pos.X())>330 ||  pos.Y()<50 || pos.Y()>550 || pos.Z()<50 || pos.Z()>645) && !(std::abs(end.X())>330 ||  end.Y()<50 || end.Y()>550 || end.Z()<50 || end.Z()>645))||(!(std::abs(pos.X())>330 ||  pos.Y()<50 || pos.Y()>550 || pos.Z()<50 || pos.Z()>645) && (std::abs(end.X())>330 ||  end.Y()<50 || end.Y()>550 || end.Z()<50 || end.Z()>645))){
         stopping_trks++; //these are total stopping tracks, could be broken as well
 
         /**********************************broken tracks removal**************************************************************************************/
         for(size_t k=0;k<NTracks;++k){
           art::Ptr<recob::Track> ptrack_k(trackListHandle, k);
           const recob::Track& track_k = *ptrack_k;
           // TVector3 pos_k, dir_pos_k, dir_end_k, end_k;
           auto pos_k = track_k.Vertex();
           auto dir_pos_k= track_k.VertexDirection();
           auto dir_end_k   = track_k.EndDirection();
           auto end_k = track_k.End();
           if(k==i) continue;
           if((std::abs(((end_k.Y()-pos_k.Y())/(end_k.Z()-pos_k.Z()))*(endz-pos_k.Z())+pos_k.Y()-endy)<30||std::abs(((end_k.Y()-pos_k.Y())/(end_k.Z()-pos_k.Z()))*(startz-pos_k.Z())+pos_k.Y()-starty)<30)&&(std::abs(endcosx*dir_pos_k.X()+endcosy*dir_pos_k.Y()+endcosz*dir_pos_k.Z())>0.97||std::abs(startcosx*dir_pos_k.X()+startcosy*dir_pos_k.Y()+startcosz*dir_pos_k.Z())>0.97||std::abs(endcosx*dir_end_k.X()+endcosy*dir_end_k.Y()+endcosz*dir_end_k.Z())>0.97||std::abs(startcosx*dir_end_k.X()+startcosy*dir_end_k.Y()+startcosz*dir_end_k.Z())>0.97)) break;
           if((std::abs(((end_k.Y()-pos_k.Y())/(end_k.Z()-pos_k.Z()))*(endz-pos_k.Z())+pos_k.Y()-endy)<50||std::abs(((end_k.Y()-pos_k.Y())/(end_k.Z()-pos_k.Z()))*(startz-pos_k.Z())+pos_k.Y()-starty)<50)&&(std::abs(endcosx*dir_pos_k.X()+endcosy*dir_pos_k.Y()+endcosz*dir_pos_k.Z())>0.998||std::abs(startcosx*dir_pos_k.X()+startcosy*dir_pos_k.Y()+startcosz*dir_pos_k.Z())>0.998||std::abs(endcosx*dir_end_k.X()+endcosy*dir_end_k.Y()+endcosz*dir_end_k.Z())>0.998||std::abs(startcosx*dir_end_k.X()+startcosy*dir_end_k.Y()+startcosz*dir_end_k.Z())>0.998)) break;
           count++;
         }        
       
         if(!(count==NTracks-1)|| tracklength<100) continue;
         unbroken_trks++;
         //*******************************new stuff*****************************//
      
         std::vector<int> wirenos;
         std::vector<float> peakts;
 
 
 
         int planenum1=999;
         float xpos=-9999;
         float ypos=-9999;
         float zpos=-9999;
         if(fmthm.isValid()){
           auto vhit=fmthm.at(i);
           auto vmeta=fmthm.data(i);
           for (size_t ii = 0; ii<vhit.size(); ++ii){ //loop over all meta data hit
             bool fBadhit = false;
             if (vmeta[ii]->Index() == static_cast<unsigned int>(std::numeric_limits<int>::max())){
               fBadhit = true;
            
               continue;
             }
             if (vmeta[ii]->Index()>=tracklist[i]->NumberTrajectoryPoints()){
               throw cet::exception("Calorimetry_module.cc") << "Requested track trajectory index "<<vmeta[ii]->Index()<<" exceeds the total number of trajectory points "<<tracklist[i]->NumberTrajectoryPoints()<<" for track index "<<i<<". Something is wrong with the track reconstruction. Please contact tjyang@fnal.gov!!";
             }
             if (!tracklist[i]->HasValidPoint(vmeta[ii]->Index())){
               fBadhit = true;
               continue;
             }
             // TVector3 loc = tracklist[i]->LocationAtPoint(vmeta[ii]->Index());
              auto loc = tracklist[i]->LocationAtPoint(vmeta[ii]->Index());
 
             xpos=loc.X();
             ypos=loc.Y();
             zpos=loc.Z();
             if (fBadhit) continue; //HY::If BAD hit, skip this hit and go next
             if (zpos<-100) continue; //hit not on track
             planenum1=vhit[ii]->WireID().Plane;
             if(planenum1==2){
               if(starty<endy) dist=sqrt(pow(xpos-startx,2)+pow(ypos-starty,2)+pow(zpos-startz,2));
               if(starty>endy) dist=sqrt(pow(xpos-endx,2)+pow(ypos-endy,2)+pow(zpos-endz,2));
               //if(vhit[ii]->WireID().Plane==2 && vhit[ii]->PeakTime()>420 && vhit[ii]->PeakTime()<530 && vhit[ii]->WireID().TPC==9 && vhit[ii]->WireID().Wire>60 && vhit[ii]->WireID().Wire<80){
               //cout<<"wire_no and peak time from metadata"<<vhit[ii]->WireID().Wire<<"  "<<vhit[ii]->PeakTime()<<endl;
 
               wirenos.push_back(vhit[ii]->WireID().Wire);
               peakts.push_back(vhit[ii]->PeakTime());
               // }
 
               if(dist<dist0){
                 dist0=dist;
                 wireno=vhit[ii]->WireID().Wire;
                 peaktime=vhit[ii]->PeakTime();
                 tpcno=vhit[ii]->WireID().TPC;
               } 
             }
           }//loop over vhit
         }//fmthm valid
         /*************************************filling the values****************************/
         counter1=0;
         // bool test=false;
         for(size_t hitl=0;hitl<hitlist.size();hitl++){
           auto & tracks = thass.at(hitlist[hitl].key());
           //auto & vmeta = thass.data(hitlist[hitl].key());
           //if (!tracks.empty()&&tracks[0].key() == ptrack.key()&&vmeta[0]->Index()!=static_cast<unsigned int>(std::numeric_limits<int>::max())) continue;
           //if (!tracks.empty()&&tracks[0].key() == ptrack.key()&&hitlist[hitl]->WireID().Plane==2)
           //if(!(hitlist[hitl]->WireID().Plane==2 && tracks[0].key() != ptrack.key() && hitlist[hitl]->PeakTime()>420 && hitlist[hitl]->PeakTime()<530 && hitlist[hitl]->WireID().TPC==9 && hitlist[hitl]->WireID().Wire>60 && hitlist[hitl]->WireID().Wire<80)) continue;
           //cout<<"wire no and peak time from hit to track association...................."<<hitlist[hitl]->WireID().Wire<<"  "<<hitlist[hitl]->PeakTime()<<endl;
           if (!tracks.empty() && tracks[0].key()!=ptrack.key() && tracklist[tracks[0].key()]->Length()>100) continue;
           bool test=true;
           float peakth1=hitlist[hitl]->PeakTime();
           int wireh1=hitlist[hitl]->WireID().Wire;
           for(size_t m=0;m<wirenos.size();m++){
             if(wireh1==wirenos[m] && peakth1==peakts[m]){
               test=false;
               break;
             }
           }
           if(!test) continue;
 
           int planeid=hitlist[hitl]->WireID().Plane;
           int tpcid=hitlist[hitl]->WireID().TPC;
           if(abs(wireh1-wireno)<6 && abs(peakth1-peaktime)<50 && planeid==2 && tpcid==tpcno){
             counter1++;
           }
         }
         // for(size_t t=0;t<peakts.size();t++) cout<<"wire numbers and peak times in the vector "<<wirenos[t]<<"  "<<peakts[t]<<endl;
 
         wirenos.clear();
         peakts.clear();
      
 
         cout<<"no of hits closeby  "<<counter1<<"   "<<"event "<<event<<" TrkackID "<<track.ID()<<" startx, y, z "<<startx<<" "<<starty<<" "<<startz<<"  wireno, peakt tpcno "<<wireno<<" "<<peaktime<<" "<<tpcno<<" dist "<<dist0<<"min T, max_T"<<min<<" "<<max<<endl; 
       }
 
      
       cross_trks++; // these are good tracks unbroken
       dist_min[cross_trks-1]=dist0;
       adjacent_hits[cross_trks-1]=counter1;
       lastwire[cross_trks-1]=wireno;
       lastpeakt[cross_trks-1]=peaktime;
       endtpc[cross_trks-1]=tpcno;
 
 
      
       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();
       peakT_max[cross_trks-1]=max;
       peakT_min[cross_trks-1]=min;
       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();
           trkpitch[cross_trks-1][planenum][iHit]=(calos[ical]->TrkPitchVec())[iHit];
         } // loop over iHit..
       } // loop over ical 2nd time...
      
     
     } // loop over trks...
   
     fEventTree->Fill();
   } // end of analyze function

void dune::michelremoving::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 160 of file michelremoving_module.cc.

                                {
     std::cout<<"job begin..."<<std::endl;
     art::ServiceHandle<art::TFileService> tfs;
     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("unbroken_trks",&unbroken_trks,"unbroken_trks/I");
     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("peakT_max",peakT_max,"peakT_max[cross_trks]/F");
     fEventTree->Branch("peakT_min",peakT_min,"peakT_min[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");
     fEventTree->Branch("trkpitch",trkpitch,"trkpitch[cross_trks][3][3000]/F");
     fEventTree->Branch("dist_min",dist_min,"dist_min[cross_trks]/F");
     fEventTree->Branch("adjacent_hits",adjacent_hits,"adjacent_hits[cross_trks]/I");
     fEventTree->Branch("lastwire",lastwire,"lastwire[cross_trks]/I");
     fEventTree->Branch("lastpeakt",lastpeakt,"lastpeakt[cross_trks]/F");
     fEventTree->Branch("endtpc",endtpc,"endtpc[cross_trks]/I");
  
   }

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

Definition at line 210 of file michelremoving_module.cc.

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

void dune::michelremoving::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 205 of file michelremoving_module.cc.

                              {     
 
   }

void dune::michelremoving::reset ( )

Definition at line 515 of file michelremoving_module.cc.

                             {
     run = -99999;
     subrun = -99999;
     event = -99999;
     evttime = -99999;
     cross_trks = -99999;
     //all_trks = -99999;
     //unbroken_trks=-99999;
     year_month_date=-99999;
     hour_min_sec=-99999;
     for(int i=0; i<kMaxTracks; i++){
       trackthetaxz[i]=-99999;
       trackthetayz[i]=-99999;
       trkstartx[i]=-99999;
       trkstarty[i]=-99999;
       trkstartz[i]=-99999;
       trkendx[i]=-99999;
       trkendy[i]=-99999;
       trkendz[i]=-99999;
       trklen[i]=-99999;
       TrkID[i]=-99999;
 
       dist_min[i]=-1;
       adjacent_hits[i]=-1;
       lastwire[i]=-1;
       lastpeakt[i]=-1;
       endtpc[i]=-1;
       peakT_max[i]=-99999;
       peakT_min[i]=-99999;
       trkstartcosxyz[i][0]=-99999;
       trkstartcosxyz[i][1]=-99999;
       trkstartcosxyz[i][2]=-99999; 
       trkendcosxyz[i][0]=-99999;
       trkendcosxyz[i][1]=-99999;
       trkendcosxyz[i][2]=-99999;
       ntrkhits[i][0] = -99999;
       ntrkhits[i][1] = -99999;
       ntrkhits[i][2] = -99999;
       for(int j=0; j<3; j++){
         for(int k=0; k<3000; k++){
           trkdqdx[i][j][k]=-99999;
           trkdedx[i][j][k]=-99999;
           trkresrange[i][j][k]=-99999;
           trkhitx[i][j][k]=-99999;
           trkhity[i][j][k]=-99999;
           trkhitz[i][j][k]=-99999;
           trkpitch[i][j][k]=-99999;
         }
       }
     }
   }

Member Data Documentation

Int_t dune::michelremoving::adjacent_hits[kMaxTracks]

private

Definition at line 131 of file michelremoving_module.cc.

Int_t dune::michelremoving::all_trks

private

Definition at line 106 of file michelremoving_module.cc.

Int_t dune::michelremoving::cross_trks

private

Definition at line 104 of file michelremoving_module.cc.

Float_t dune::michelremoving::dist_min[kMaxTracks]

private

Definition at line 130 of file michelremoving_module.cc.

Int_t dune::michelremoving::endtpc[kMaxTracks]

private

Definition at line 133 of file michelremoving_module.cc.

Int_t dune::michelremoving::event

private

Definition at line 100 of file michelremoving_module.cc.

Double_t dune::michelremoving::evttime

private

Definition at line 101 of file michelremoving_module.cc.

std::string dune::michelremoving::fCalorimetryModuleLabel

private

Definition at line 137 of file michelremoving_module.cc.

TTree* dune::michelremoving::fEventTree

private

Definition at line 97 of file michelremoving_module.cc.

std::string dune::michelremoving::fHitsModuleLabel

private

Definition at line 135 of file michelremoving_module.cc.

bool dune::michelremoving::fSaveCaloInfo

private

Definition at line 138 of file michelremoving_module.cc.

bool dune::michelremoving::fSaveTrackInfo

private

Definition at line 139 of file michelremoving_module.cc.

std::string dune::michelremoving::fTrackModuleLabel

private

Definition at line 136 of file michelremoving_module.cc.

Int_t dune::michelremoving::hour_min_sec

private

Definition at line 103 of file michelremoving_module.cc.

Float_t dune::michelremoving::lastpeakt[kMaxTracks]

private

Definition at line 134 of file michelremoving_module.cc.

Int_t dune::michelremoving::lastwire[kMaxTracks]

private

Definition at line 132 of file michelremoving_module.cc.

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

private

Definition at line 120 of file michelremoving_module.cc.

Float_t dune::michelremoving::peakT_max[kMaxTracks]

private

Definition at line 128 of file michelremoving_module.cc.

Float_t dune::michelremoving::peakT_min[kMaxTracks]

private

Definition at line 129 of file michelremoving_module.cc.

Int_t dune::michelremoving::run

private

Definition at line 98 of file michelremoving_module.cc.

Int_t dune::michelremoving::stopping_trks

private

Definition at line 105 of file michelremoving_module.cc.

Int_t dune::michelremoving::subrun

private

Definition at line 99 of file michelremoving_module.cc.

Float_t dune::michelremoving::trackthetaxz[kMaxTracks]

private

Definition at line 108 of file michelremoving_module.cc.

Float_t dune::michelremoving::trackthetayz[kMaxTracks]

private

Definition at line 109 of file michelremoving_module.cc.

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

private

Definition at line 122 of file michelremoving_module.cc.

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

private

Definition at line 121 of file michelremoving_module.cc.

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

private

Definition at line 119 of file michelremoving_module.cc.

Float_t dune::michelremoving::trkendx[kMaxTracks]

private

Definition at line 113 of file michelremoving_module.cc.

Float_t dune::michelremoving::trkendy[kMaxTracks]

private

Definition at line 114 of file michelremoving_module.cc.

Float_t dune::michelremoving::trkendz[kMaxTracks]

private

Definition at line 115 of file michelremoving_module.cc.

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

private

Definition at line 124 of file michelremoving_module.cc.

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

private

Definition at line 125 of file michelremoving_module.cc.

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

private

Definition at line 126 of file michelremoving_module.cc.

Int_t dune::michelremoving::TrkID[kMaxTracks]

private

Definition at line 117 of file michelremoving_module.cc.

Float_t dune::michelremoving::trklen[kMaxTracks]

private

Definition at line 116 of file michelremoving_module.cc.

Float_t dune::michelremoving::trkpitch[kMaxTracks][3][3000]

private

Definition at line 127 of file michelremoving_module.cc.

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

private

Definition at line 123 of file michelremoving_module.cc.

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

private

Definition at line 118 of file michelremoving_module.cc.

Float_t dune::michelremoving::trkstartx[kMaxTracks]

private

Definition at line 110 of file michelremoving_module.cc.

Float_t dune::michelremoving::trkstarty[kMaxTracks]

private

Definition at line 111 of file michelremoving_module.cc.

Float_t dune::michelremoving::trkstartz[kMaxTracks]

private

Definition at line 112 of file michelremoving_module.cc.

Int_t dune::michelremoving::unbroken_trks

private

Definition at line 107 of file michelremoving_module.cc.

Int_t dune::michelremoving::year_month_date

private

Definition at line 102 of file michelremoving_module.cc.

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

protoduneana/protoduneana/singlephase/michelremoving/michelremoving_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation