CRT::TwoCRTReco Class Reference

Inheritance diagram for CRT::TwoCRTReco:

Classes
struct	recoHits
struct	tempHits
struct	tracksPair

Public Member Functions
	TwoCRTReco (fhicl::ParameterSet const &p)
	TwoCRTReco (TwoCRTReco const &)=delete
	TwoCRTReco (TwoCRTReco &&)=delete
TwoCRTReco &	operator= (TwoCRTReco const &)=delete
TwoCRTReco &	operator= (TwoCRTReco &&)=delete
void	analyze (art::Event const &e) override
bool	moduleMatcher (int module1, int module2)
void	beginJob () override
void	endJob () override
void	createPNG (TH1D *histo)
double	setAngle (double angle)
int	moduletoCTB (int module2, int module1)
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)

Public Attributes
std::string	fTrackModuleLabel = "pandoraTrack"
int	nEvents = 0
int	nHaloMuons = 0
int	nHitsPerEvent =0
int	nMCCMuons =0

Private Attributes
ofstream	logFile
art::InputTag	fCRTLabel
art::InputTag	fCTBLabel
TTree *	fCRTTree
int	run
int	subRun
bool	fMCCSwitch
bool	fCTBTriggerOnly
bool	fSCECorrection
bool	fModuleSwitch
int	fADCThreshold
int	fModuletoModuleTimingCut
int	fFronttoBackTimingCut
int	adcX_F
int	adcX_B
int	adcY_F
int	adcY_B
double	X_F
double	X_B
double	Y_F
double	Y_B
double	Z_F
double	Z_B
int	moduleX_F
int	moduleX_B
int	moduleY_F
int	moduleY_B
int	stripX_F
int	stripX_B
int	stripY_F
int	stripY_B
double	measuredT0
double	CRT_TOF
long long	timeStamp
int	eventNum
std::vector< recoHits >	primaryHits_F
std::vector< recoHits >	primaryHits_B
std::vector< tempHits >	tempHits_F
std::vector< tempHits >	tempHits_B
std::vector< tracksPair >	allTracksPair

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 87 of file TwoCRTReco_module.cc.

Constructor & Destructor Documentation

CRT::TwoCRTReco::TwoCRTReco ( fhicl::ParameterSet const &  p )

explicit

Definition at line 210 of file TwoCRTReco_module.cc.

              :
  EDAnalyzer(p), fCRTLabel(p.get < art::InputTag > ("CRTLabel")),  fCTBLabel(p.get<art::InputTag>("CTBLabel")) {
    consumes < std::vector < CRT::Trigger >> (fCRTLabel);
    consumes < std::vector < art::Assns < sim::AuxDetSimChannel, CRT::Trigger >>> (fCRTLabel); 
  fMCCSwitch=(p.get<bool>("MCC"));
  fCTBTriggerOnly=(p.get<bool>("CTBOnly"));
  fSCECorrection=(p.get<bool>("SCECorrection"));
  }

CRT::TwoCRTReco::TwoCRTReco ( TwoCRTReco const &  )

delete

CRT::TwoCRTReco::TwoCRTReco ( TwoCRTReco &&  )

delete

Member Function Documentation

void CRT::TwoCRTReco::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 291 of file TwoCRTReco_module.cc.

{

    nEvents++;  
    eventNum=event.event();
    run=event.run();
    subRun=event.subRun();
  if (fMCCSwitch){
    fModuleSwitch=1;
    fADCThreshold=800;
    fModuletoModuleTimingCut=40;
    fFronttoBackTimingCut=120;
    
}
  else {
    fModuleSwitch=0;
    fADCThreshold=20;
    fModuletoModuleTimingCut=5;
    fFronttoBackTimingCut=8;
    art::ValidHandle<std::vector<raw::RDTimeStamp>> timingHandle = event.getValidHandle<std::vector<raw::RDTimeStamp>>("timingrawdecoder:daq");
    timeStamp=timingHandle->at(0).GetTimeStamp();

} 



   if(!fMCCSwitch){
   //const auto& pdspctbs = event.getValidHandle<std::vector<raw::ctb::pdspctb>>(fCTB_tag);
        art::ValidHandle<std::vector<raw::RDTimeStamp>> timingHandle = event.getValidHandle<std::vector<raw::RDTimeStamp>>("timingrawdecoder:daq");

        const raw::RDTimeStamp& timeStamp = timingHandle->at(0);
        if (fCTBTriggerOnly){
        if(timeStamp.GetFlags()!= 13) return;}
  }
  int nHits = 0;




        //Detector properties service
  primaryHits_F.clear();
  primaryHits_B.clear();
  allTracksPair.clear();
  tempHits_F.clear();
  tempHits_B.clear(); // Arrays to compile hits and move them through


  //Get triggers
  cout << "Getting triggers" << endl;
  const auto & triggers = event.getValidHandle < std::vector < CRT::Trigger >> (fCRTLabel);

  art::FindManyP < sim::AuxDetSimChannel > trigToSim(triggers, event, fCRTLabel);

  //Get a handle to the Geometry service to look up AuxDetGeos from module numbers
  art::ServiceHandle < geo::Geometry > geom;



  //Mapping from channel to trigger
  std::unordered_map < size_t, double > prevTimes;
  int hitID = 0;
  cout << "Looking for hits in Triggers" << endl;

   // Find CTB pixels
        int pixel0 = -1;
        int pixel1 = -1;
        if (!fMCCSwitch)
        {
       const auto& pdspctbs = *event.getValidHandle<std::vector<raw::ctb::pdspctb>>(fCTBLabel);
       std::vector<int> uS, dS;
        const size_t npdspctbs = pdspctbs.size();
        for(size_t j=0;j<npdspctbs;++j)
          {
            const std::vector<raw::ctb::Trigger> HLTriggers = pdspctbs[j].GetHLTriggers();
            const std::vector<raw::ctb::ChStatus> chs = pdspctbs[j].GetChStatusAfterHLTs();
for (size_t k=0; k<HLTriggers.size(); ++k)
              { 
                //cout<<chs[k].timestamp<<endl;
                int num = chs[k].crt;
                //cout<<num<<endl;
        
        const std::string binary = std::bitset<32>(num).to_string();    
        const auto crtmask=chs[k].crt;
         pixel0 = -1;
         pixel1 = -1;
        //cout<<crtmask<<endl;
        for (int i = 0; i<32; ++i){
          if (crtmask & (1<<i)){
            if (i<16){
              pixel0 = i;
            }
            else {
              pixel1 = i;
            }
          }
        }
        if (pixel0!=-1 && pixel1!=-1) {
        //cout<<nEvents<<" TJYang Pixels: "<<pixel0<<","<<pixel1<<endl;
        }
        else if (fCTBTriggerOnly) return;
              }
          }
        }

  int trigID=0;
  for (const auto & trigger: * triggers) {
    const auto & hits = trigger.Hits();
    for (const auto & hit: hits) { // Collect hits on all modules
        //cout<<hits.size()<<','<<hit.ADC()<<endl;
      if (hit.ADC() > fADCThreshold) { // Keep if they are above threshold

        tempHits tHits;
        if (!fMCCSwitch){

        tHits.module = trigger.Channel(); // Values to add to array
        tHits.channel=hit.Channel();
        tHits.adc = hit.ADC();
        tHits.triggerTime=trigger.Timestamp()-timeStamp;
        }
        else{
        tHits.module = trigger.Channel(); // Values to add to array
        tHits.channel=hit.Channel();
        tHits.adc = hit.ADC();
        tHits.triggerTime=trigger.Timestamp();
        }
         //cout<<trigger.Channel()<<','<<hit.Channel()<<','<<hit.ADC()<<endl;
        nHits++;
        tHits.triggerNumber=trigID;
        const auto & trigGeo = geom -> AuxDet(trigger.Channel()); // Get geo  
        const auto & csens = trigGeo.SensitiveVolume(hit.Channel());
        const auto center = csens.GetCenter();
        if (center.Z() < 100) tempHits_F.push_back(tHits); // Sort F/B from Z
        else tempHits_B.push_back(tHits);
        hitID++;
      }
    }
    trigID++;
  }
  nHitsPerEvent=nHits;
  cout << "Hits compiled for event: " << nEvents << endl;
  cout << "Number of Hits above Threshold:  " << hitID << endl;

  for (unsigned int f = 0; f < tempHits_F.size(); f++) {
    for (unsigned int f_test = 0; f_test < tempHits_F.size(); f_test++) {
       if (fabs(tempHits_F[f_test].triggerTime-tempHits_F[f].triggerTime)>fModuletoModuleTimingCut) continue;
      const auto & trigGeo = geom -> AuxDet(tempHits_F[f].module);
      const auto & trigGeo2 = geom -> AuxDet(tempHits_F[f_test].module);

      const auto & hit1Geo = trigGeo.SensitiveVolume(tempHits_F[f].channel);
      const auto hit1Center = hit1Geo.GetCenter();
      // Create 2D hits from geo of the Y and X modules
       const auto & hit2Geo = trigGeo2.SensitiveVolume(tempHits_F[f_test].channel);
      const auto hit2Center = hit2Geo.GetCenter();
      bool moduleMatched;
      moduleMatched=moduleMatcher(tempHits_F[f_test].module, tempHits_F[f].module);
      if (moduleMatched) {
        // Get the center of the hits (CRT_Res=2.5 cm)
        double hitX = hit1Center.X();
        for (unsigned int a = 0; a < tempHits_F.size(); a++)
        {
        if(tempHits_F[a].module==tempHits_F[f].module && (tempHits_F[a].channel-1)==tempHits_F[f].channel) hitX=hit1Center.X()+1.25;
        }
        double hitYPrelim=hit2Center.Y();
        for (unsigned int a = 0; a < tempHits_F.size(); a++)
        {
        if(tempHits_F[a].module==tempHits_F[f_test].module && (tempHits_F[a].channel-1)==tempHits_F[f_test].channel) hitYPrelim=hit2Center.Y()+1.25;
        }
        

        
        double hitY=hitYPrelim;
        double hitZ = (hit1Center.Z() + hit2Center.Z()) / 2.f;

        recoHits rHits;
        rHits.hitPositionX = hitX;
        rHits.hitPositionY = hitY;
        rHits.hitPositionZ = hitZ;
        rHits.geoX=tempHits_F[f].module;
        rHits.geoY=tempHits_F[f_test].module;
        rHits.stripX=tempHits_F[f].channel;
        rHits.stripY=tempHits_F[f_test].channel;

        rHits.adcX=tempHits_F[f].adc;
        rHits.adcY=tempHits_F[f_test].adc;
        rHits.trigNumberX=tempHits_F[f].triggerNumber;
        rHits.trigNumberY=tempHits_F[f_test].triggerNumber;
        rHits.timeAvg = (tempHits_F[f_test].triggerTime+tempHits_F[f].triggerTime)/2.0;
       primaryHits_F.push_back(rHits); // Add array
    }
    }
  }
  for (unsigned int f = 0; f < tempHits_B.size(); f++) {
    for (unsigned int f_test = 0; f_test < tempHits_B.size(); f_test++) { // Same as above but for back CRT
       if (fabs(tempHits_B[f_test].triggerTime-tempHits_B[f].triggerTime)>fModuletoModuleTimingCut) continue;

      const auto & trigGeo = geom -> AuxDet(tempHits_B[f].module);
      const auto & trigGeo2 = geom -> AuxDet(tempHits_B[f_test].module);
      const auto & hit1Geo = trigGeo.SensitiveVolume(tempHits_B[f].channel);
      const auto hit1Center = hit1Geo.GetCenter();

      const auto & hit2Geo = trigGeo2.SensitiveVolume(tempHits_B[f_test].channel);
      const auto hit2Center = hit2Geo.GetCenter();
      bool moduleMatched;
      moduleMatched=moduleMatcher(tempHits_B[f_test].module, tempHits_B[f].module);

      if (moduleMatched) {
        double hitX = hit1Center.X();
        
        
        for (unsigned int a = 0; a < tempHits_B.size(); a++)
        {
        if(tempHits_B[a].module==tempHits_B[f].module && (tempHits_B[a].channel-1)==tempHits_B[f].channel) hitX=hit1Center.X()+1.25;
        }
        
        double hitYPrelim = hit2Center.Y();
        
        for (unsigned int a = 0; a < tempHits_B.size(); a++)
        {
        if(tempHits_B[a].module==tempHits_B[f_test].module && (tempHits_B[a].channel-1)==tempHits_B[f_test].channel) hitYPrelim=hit2Center.Y()+1.25;
        }
        double hitY=hitYPrelim;

        
        double hitZ = (hit1Center.Z() + hit2Center.Z()) / 2.f;

        recoHits rHits;
        rHits.hitPositionX = hitX;
        rHits.hitPositionY = hitY;
        rHits.hitPositionZ = hitZ;
        rHits.geoX=tempHits_B[f].module;
        rHits.geoY=tempHits_B[f_test].module;
        rHits.adcX=tempHits_B[f].adc;
        rHits.adcY=tempHits_B[f_test].adc;
        rHits.stripX=tempHits_B[f].channel;
        rHits.stripY=tempHits_B[f_test].channel;
        rHits.trigNumberX=tempHits_B[f].triggerNumber;
        rHits.trigNumberY=tempHits_B[f_test].triggerNumber;
        rHits.timeAvg = (tempHits_B[f_test].triggerTime+tempHits_B[f].triggerTime)/2.0;
         primaryHits_B.push_back(rHits); 

         }
    }
  }

        std::cout<<"Number of Hits: "<<primaryHits_F.size()<<','<<primaryHits_B.size()<<std::endl;

        std::cout<<"CTB Pixes: "<<pixel0<<','<<pixel1<<std::endl;
    int tempId = 0;
    adcX_F=0;
    adcX_B=0;
    adcY_F=0;
    adcY_B=0;
    int crtPixel0=-1;
    int crtPixel1=-1;
    for (unsigned int f = 0; f < primaryHits_F.size(); f++) {
      if (pixel0==-1 || pixel1==-1) break; 
      for (unsigned int b = 0; b < primaryHits_B.size(); b++) {
      if (fabs(primaryHits_F[f].timeAvg-primaryHits_B[b].timeAvg)>fFronttoBackTimingCut) continue;
      crtPixel0=moduletoCTB(primaryHits_F[f].geoX,primaryHits_F[f].geoY);
      crtPixel1=moduletoCTB(primaryHits_B[b].geoX, primaryHits_B[b].geoY);
      if (crtPixel0!=pixel0 || crtPixel1!=pixel1) continue;     
      //std::cout<<"HEY"<<std::endl;
      if (adcX_F<primaryHits_F[f].adcX && adcY_F<primaryHits_F[f].adcY && adcX_B<primaryHits_B[b].adcX  && adcY_B<primaryHits_B[b].adcY){
      X_F = primaryHits_F[f].hitPositionX;
      Y_F = primaryHits_F[f].hitPositionY;
      Z_F = primaryHits_F[f].hitPositionZ;
      X_B = primaryHits_B[b].hitPositionX;
      Y_B = primaryHits_B[b].hitPositionY;
      Z_B= primaryHits_B[b].hitPositionZ;
      adcX_F=primaryHits_F[f].adcX;
      adcX_B=primaryHits_B[b].adcX;
      adcY_F=primaryHits_F[f].adcY;
      adcY_B=primaryHits_B[b].adcY;
      moduleX_F=primaryHits_F[f].geoX;
      moduleY_F=primaryHits_F[f].geoY;
      moduleX_B=primaryHits_B[b].geoX;
      moduleY_B=primaryHits_B[b].geoY;

      stripX_F=primaryHits_F[f].stripX;
      stripY_F=primaryHits_F[f].stripY;
      stripX_B=primaryHits_B[b].stripX;
      stripY_B=primaryHits_B[b].stripY;
      //std::cout<<"FOUND A COMBO"<<std::endl;
      measuredT0=(primaryHits_F[f].timeAvg+primaryHits_B[b].timeAvg)/2.f;
        CRT_TOF=(primaryHits_F[f].timeAvg-primaryHits_B[b].timeAvg);
        tempId++;

        }       
        }
      }
    if (adcX_F>0) fCRTTree->Fill();
    // Filter return if (adcX_F==0) return false;
    // Filter return
    // return true;
 }

void CRT::TwoCRTReco::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 591 of file TwoCRTReco_module.cc.

                             {
        art::ServiceHandle<art::TFileService> fileServiceHandle;
       fCRTTree = fileServiceHandle->make<TTree>("CRTCandidate", "track by track info");

        fCRTTree->Branch("nEvents", &nEvents, "fnEvents/I");
        fCRTTree->Branch("hRun", &run, "run/I");
        fCRTTree->Branch("hEvent",&eventNum,"eventNum/I");


        fCRTTree->Branch("hX_F", &X_F, "X_F/D");
        fCRTTree->Branch("hX_B", &X_B, "X_B/D");
        fCRTTree->Branch("hY_F", &Y_F, "Y_F/D");
        fCRTTree->Branch("hY_B", &Y_B, "Y_B/D");
        fCRTTree->Branch("hZ_F", &Z_F, "Z_F/D");
        fCRTTree->Branch("hZ_B", &Z_B, "Z_B/D");


        fCRTTree->Branch("CRT_TOF", &CRT_TOF, "CRT_TOF/D");

        fCRTTree->Branch("hmoduleX_F", &moduleX_F, "moduleX_F/I");
        fCRTTree->Branch("hmoduleX_B", &moduleX_B, "moduleX_B/I");
        fCRTTree->Branch("hmoduleY_F", &moduleY_F, "moduleY_F/I");
        fCRTTree->Branch("hmoduleY_B", &moduleY_B, "moduleY_B/I");

        fCRTTree->Branch("hstripX_F", &stripX_F, "stripX_F/I");
        fCRTTree->Branch("hstripX_B", &stripX_B, "stripX_B/I");
        fCRTTree->Branch("hstripY_F", &stripY_F, "stripY_F/I");
        fCRTTree->Branch("hstripY_B", &stripY_B, "stripY_B/I");

        fCRTTree->Branch("hadcX_F", &adcX_F, "adcX_F/I");
        fCRTTree->Branch("hadcY_F", &adcY_F, "adcY_F/I");
        fCRTTree->Branch("hadcX_B", &adcX_B, "adcX_B/I");
        fCRTTree->Branch("hadcY_B", &adcY_B, "adcY_B/I");

        fCRTTree->Branch("hmeasuredT0", &measuredT0, "measuredT0/D");





}

void CRT::TwoCRTReco::createPNG

(

TH1D *

histo

)

Definition at line 269 of file TwoCRTReco_module.cc.

                                          {
  //Save important histograms as PNG
  TCanvas * c = new TCanvas;
  TH1D * h = histo;
  h -> Draw();
  TImage * img = TImage::Create();
  img -> FromPad(c);
  img -> WriteImage((std::string(histo -> GetName()) + ".png").c_str());
  delete h;
  delete c;
  delete img;

}

void CRT::TwoCRTReco::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 633 of file TwoCRTReco_module.cc.

{



}

bool CRT::TwoCRTReco::moduleMatcher	(	int	module1,
		int	module2
	)

Definition at line 224 of file TwoCRTReco_module.cc.

                                                          {
  // Function checking if two hits could reasonably be matched into a 2D hit
  if ((module1 == 6 && (module2 == 10 || module2 == 11)) || (module1 == 14 && (module2 == 10 || module2 == 11)) || (module1 == 19 && (module2 == 26 || module2 == 27)) || (module1 == 31 && (module2 == 26 || module2 == 27)) || (module1 == 7 && (module2 == 12 || module2 == 13)) || (module1 == 15 && (module2 == 12 || module2 == 13)) || (module1 == 18 && (module2 == 24 || module2 == 25)) || (module1 == 30 && (module2 == 24 || module2 == 25)) || (module1 == 1 && (module2 == 4 || module2 == 5)) || (module1 == 9 && (module2 == 4 || module2 == 5)) || (module1 == 16 && (module2 == 20 || module2 == 21)) || (module1 == 28 && (module2 == 20 || module2 == 21)) || (module1 == 0 && (module2 == 2 || module2 == 3)) || (module1 == 8 && (module2 == 2 || module2 == 3)) || (module1 == 17 && (module2 == 22 || module2 == 23)) || (module1 == 29 && (module2 == 22 || module2 == 23))) return 1;
  else return 0;

}

int CRT::TwoCRTReco::moduletoCTB	(	int	module2,
		int	module1
	)

Definition at line 232 of file TwoCRTReco_module.cc.

                                                      {
  if (module1 == 15 && module2 == 12 ) return 15;
  else if (module1 == 15 &&  module2 == 13) return 10;
  else if (module1 == 7 &&  module2 == 12) return 8;
  else if (module1 == 7 &&  module2 == 13) return 9;
  else if (module1 == 19 &&  module2 == 26) return 4;
  else if (module1 == 19 &&  module2 == 27) return 13;
  else if (module1 == 31 &&  module2 == 26) return 3;
  else if (module1 == 31 &&  module2 == 27) return 2;
  else if (module1 == 30 &&  module2 == 24) return 1;
  else if (module1 == 30 &&  module2 == 25) return 0;
  else if (module1 == 18 &&  module2 == 24) return 12;
  else if (module1 == 18 &&  module2 == 25) return 11;
  else if (module1 == 6  &&  module2 == 11) return 7;
  else if (module1 == 6 &&  module2 == 10) return 6;
  else if (module1 == 14  &&  module2 == 11) return 14;
  else if (module1 == 14 &&  module2 == 10) return 5;
  else if (module1 == 0 &&  module2 == 3) return 25;
  else if (module1 == 0 &&  module2 == 2) return 24;
  else if (module1 == 8 &&  module2 == 3) return 26;
  else if (module1 == 8 &&  module2 == 2) return 30;
  else if (module1 == 17 &&  module2 == 23) return 27;
  else if (module1 == 17 &&  module2 == 22) return 28;
  else if (module1 == 29 &&  module2 == 23) return 16;
  else if (module1 == 29 &&  module2 == 22) return 17;
  else if (module1 == 28 &&  module2 == 21) return 18;
  else if (module1 == 28 &&  module2 == 20) return 19;
  else if (module1 == 16 &&  module2 == 21) return 29;
  else if (module1 == 16 &&  module2 == 20) return 20;
  else if (module1 == 9 &&  module2 == 5) return 31;
  else if (module1 == 9 &&  module2 == 4) return 21;
  else if (module1 == 1 &&  module2 == 5) return 23;
  else if (module1 == 1 &&  module2 == 4) return 22;
  else return -1;
}

TwoCRTReco& CRT::TwoCRTReco::operator= ( TwoCRTReco const &  )

delete

TwoCRTReco& CRT::TwoCRTReco::operator= ( TwoCRTReco &&  )

delete

double CRT::TwoCRTReco::setAngle

(

double

angle

)

Definition at line 282 of file TwoCRTReco_module.cc.

                                           {
  if (angle < 0) {
    angle += 3.14159265359;
  }
  angle *= 180.00 / 3.14159265359;
  return angle;
}

Member Data Documentation

int CRT::TwoCRTReco::adcX_B

private

Definition at line 134 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::adcX_F

private

Definition at line 134 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::adcY_B

private

Definition at line 134 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::adcY_F

private

Definition at line 134 of file TwoCRTReco_module.cc.

std::vector< tracksPair > CRT::TwoCRTReco::allTracksPair

private

Definition at line 206 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::CRT_TOF

private

Definition at line 140 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::eventNum

private

Definition at line 142 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::fADCThreshold

private

Definition at line 127 of file TwoCRTReco_module.cc.

art::InputTag CRT::TwoCRTReco::fCRTLabel

private

Definition at line 118 of file TwoCRTReco_module.cc.

TTree* CRT::TwoCRTReco::fCRTTree

private

Definition at line 120 of file TwoCRTReco_module.cc.

art::InputTag CRT::TwoCRTReco::fCTBLabel

private

Definition at line 119 of file TwoCRTReco_module.cc.

bool CRT::TwoCRTReco::fCTBTriggerOnly

private

Definition at line 124 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::fFronttoBackTimingCut

private

Definition at line 129 of file TwoCRTReco_module.cc.

bool CRT::TwoCRTReco::fMCCSwitch

private

Definition at line 123 of file TwoCRTReco_module.cc.

bool CRT::TwoCRTReco::fModuleSwitch

private

Definition at line 126 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::fModuletoModuleTimingCut

private

Definition at line 128 of file TwoCRTReco_module.cc.

bool CRT::TwoCRTReco::fSCECorrection

private

Definition at line 125 of file TwoCRTReco_module.cc.

std::string CRT::TwoCRTReco::fTrackModuleLabel = "pandoraTrack"

Definition at line 91 of file TwoCRTReco_module.cc.

ofstream CRT::TwoCRTReco::logFile

private

Definition at line 115 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::measuredT0

private

Definition at line 139 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::moduleX_B

private

Definition at line 137 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::moduleX_F

private

Definition at line 137 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::moduleY_B

private

Definition at line 137 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::moduleY_F

private

Definition at line 137 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::nEvents = 0

Definition at line 111 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::nHaloMuons = 0

Definition at line 112 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::nHitsPerEvent =0

Definition at line 113 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::nMCCMuons =0

Definition at line 114 of file TwoCRTReco_module.cc.

std::vector< recoHits > CRT::TwoCRTReco::primaryHits_B

private

Definition at line 202 of file TwoCRTReco_module.cc.

std::vector< recoHits > CRT::TwoCRTReco::primaryHits_F

private

Definition at line 201 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::run

private

Definition at line 122 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::stripX_B

private

Definition at line 138 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::stripX_F

private

Definition at line 138 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::stripY_B

private

Definition at line 138 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::stripY_F

private

Definition at line 138 of file TwoCRTReco_module.cc.

int CRT::TwoCRTReco::subRun

private

Definition at line 122 of file TwoCRTReco_module.cc.

std::vector< tempHits > CRT::TwoCRTReco::tempHits_B

private

Definition at line 205 of file TwoCRTReco_module.cc.

std::vector< tempHits > CRT::TwoCRTReco::tempHits_F

private

Definition at line 204 of file TwoCRTReco_module.cc.

long long CRT::TwoCRTReco::timeStamp

private

Definition at line 141 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::X_B

private

Definition at line 135 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::X_F

private

Definition at line 135 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::Y_B

private

Definition at line 135 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::Y_F

private

Definition at line 135 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::Z_B

private

Definition at line 135 of file TwoCRTReco_module.cc.

double CRT::TwoCRTReco::Z_F

private

Definition at line 135 of file TwoCRTReco_module.cc.

---

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

• duneprototypes/duneprototypes/Protodune/singlephase/CRT/TwoCRTReco_module.cc