Inheritance diagram for simfilter::FilterParticlesActiveVolume:

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

virtual	~FilterParticlesActiveVolume ()

bool	filter (art::Event &)

virtual void	reconfigure (fhicl::ParameterSet const &)

virtual void	beginJob ()

Private Attributes
double	fXmax

double	fXmin

double	fZmax

double	fZmin

double	fYmax

double	fYmin

int	finActive

int	filtpart

Detailed Description

Definition at line 48 of file FilterParticlesActiveVolume_module.cc.

Constructor & Destructor Documentation

simfilter::FilterParticlesActiveVolume::FilterParticlesActiveVolume ( fhicl::ParameterSet const & pset )

explicit

Definition at line 85 of file FilterParticlesActiveVolume_module.cc.

                                                                                         : EDFilter{pset}
   {
     this->reconfigure(pset);
   }

simfilter::FilterParticlesActiveVolume::~FilterParticlesActiveVolume ( )

virtual

Definition at line 92 of file FilterParticlesActiveVolume_module.cc.

93 {

94 }

Member Function Documentation

void simfilter::FilterParticlesActiveVolume::beginJob ( )

virtual

Definition at line 97 of file FilterParticlesActiveVolume_module.cc.

   {
     art::ServiceHandle<geo::Geometry> geo;
   
   }

bool simfilter::FilterParticlesActiveVolume::filter ( art::Event & evt )

Definition at line 129 of file FilterParticlesActiveVolume_module.cc.

   {
     bool interactionDesired(false);
     //get the list of particles from this event
     art::ServiceHandle<geo::Geometry> geom;
 
       
     // * MC truth information
 
       //std::vector< art::Handle< std::vector<simb::MCTruth> > > allmclists;
       //evt.getManyByType(allmclists);
       auto allmclists = evt.getMany< std::vector<simb::MCTruth> >();
 
 
       std::cout << fXmin << " " << fXmax << " " << fYmin << " " <<fYmax << " " << fZmin << " " << fZmax << std::endl;
  
       geo::CryostatGeo const& cryo0 = geom->Cryostat(0);
       geo::CryostatGeo const& cryo1 = geom->Cryostat(1);
 
       geo::TPCGeo const& tpc00 = cryo0.TPC(0);
       TVector3 xyzcenter00 = tpc00.GetActiveVolumeCenter();
       std::cout << xyzcenter00[0] << " " << xyzcenter00[1] << " " << xyzcenter00[2] << std::endl;
 
       geo::TPCGeo const& tpc01 = cryo0.TPC(1);
       TVector3 xyzcenter01 = tpc01.GetActiveVolumeCenter();
       std::cout << xyzcenter01[0] << " " << xyzcenter01[1] << " " << xyzcenter01[2] << std::endl;
 
       geo::TPCGeo const& tpc10 = cryo1.TPC(0);
       TVector3 xyzcenter10 = tpc10.GetActiveVolumeCenter();
       std::cout << xyzcenter10[0] << " " << xyzcenter10[1] << " " << xyzcenter10[2] << std::endl;
 
       geo::TPCGeo const& tpc11 = cryo1.TPC(1);
       TVector3 xyzcenter11 = tpc11.GetActiveVolumeCenter();
       std::cout << xyzcenter11[0] << " " << xyzcenter11[1] << " " << xyzcenter11[2] << std::endl;
 
       double h00=tpc00.ActiveHalfHeight();
       double h01=tpc01.ActiveHalfHeight();
       double h10=tpc10.ActiveHalfHeight();
       double h11=tpc11.ActiveHalfHeight();
       double w00=tpc00.ActiveHalfWidth();
       double w01=tpc01.ActiveHalfWidth();
       double w10=tpc10.ActiveHalfWidth();
       double w11=tpc11.ActiveHalfWidth();
       double l00=tpc00.ActiveLength();
       double l01=tpc01.ActiveLength();
       double l10=tpc10.ActiveLength();
       double l11=tpc11.ActiveLength();
       
       std::cout << h00 << " " << w00 << " " << l00 << std::endl;
       std::cout << h01 << " " << w01 << " " << l01 << std::endl;
       std::cout << h10 << " " << w10 << " " << l10 << std::endl;
       std::cout << h11 << " " << w11 << " " << l11 << std::endl;
       
       for(size_t mcl = 0; mcl < allmclists.size(); ++mcl){
           art::Handle< std::vector<simb::MCTruth> > mclistHandle = allmclists[mcl];
           for(size_t m = 0; m < mclistHandle->size(); ++m){
               art::Ptr<simb::MCTruth> mct(mclistHandle, m);
               for(int ipart=0;ipart<mct->NParticles();ipart++){
                   int pdg=mct->GetParticle(ipart).PdgCode();
                   double xx=mct->GetParticle(ipart).Vx();
                   double yy=mct->GetParticle(ipart).Vy();
                   double zz=mct->GetParticle(ipart).Vz();
       
                   if (finActive==1 && pdg==filtpart)
                     {
                       if (xx>(xyzcenter00[0]-w00) && xx<(xyzcenter00[0]+w00) && yy>(xyzcenter00[1]-h00) && yy<(xyzcenter00[1]+h00) && zz>(xyzcenter00[2]-l00/2) && zz<(xyzcenter00[2]+l00/2))
                         {
                           interactionDesired = true;
                         }
                       if (xx>(xyzcenter01[0]-w01) && xx<(xyzcenter01[0]+w01) && yy>(xyzcenter01[1]-h01) && yy<(xyzcenter01[1]+h01) && zz>(xyzcenter01[2]-l01/2) && zz<(xyzcenter01[2]+l01/2))
                         {
                           interactionDesired = true;
                         }
                       if (xx>(xyzcenter10[0]-w10) && xx<(xyzcenter10[0]+w10) && yy>(xyzcenter10[1]-h10) && yy<(xyzcenter10[1]+h10) && zz>(xyzcenter10[2]-l10/2) && zz<(xyzcenter10[2]+l10/2))
                         {
                           interactionDesired = true;
                         }
                       if (xx>(xyzcenter11[0]-w11) && xx<(xyzcenter11[0]+w11) && yy>(xyzcenter11[1]-h11) && yy<(xyzcenter11[1]+h11) && zz>(xyzcenter11[2]-l11/2) && zz<(xyzcenter11[2]+l11/2))
                         {
                           interactionDesired = true;
                         }
                         
                       if (finActive==0 && pdg==filtpart)
                         {
                           if (xx>fXmin && xx<fXmax && yy>fYmin && yy<fYmax && zz>fZmin && zz<fZmax)
                             {
                               interactionDesired = true;
                             }
                         }
         
                     }
           //    std::cout << "FilterNoDirtParticles: i is " << i << std::endl ;
           // Now walk through trajectory and see if it enters the TPC
            // trajectory loop
          // end Genie particle
               }
           }
       }// loop on MCPHandle
 
   return interactionDesired;
     
   } // end FilterNoDirtParticles()function

void simfilter::FilterParticlesActiveVolume::reconfigure ( fhicl::ParameterSet const & p )

virtual

Definition at line 104 of file FilterParticlesActiveVolume_module.cc.

   {
       finActive        = p.get< int >("inActive");
       if (finActive==0) {
       fXmax            = p.get< double >("Xmax");
       fYmax            = p.get< double >("Ymax");
       fZmax            = p.get< double >("Zmax");
       fXmin            = p.get< double >("Xmin");
       fYmin            = p.get< double >("Ymin");
       fZmin            = p.get< double >("Zmin");
       }
       if (finActive==1) {
           fXmax=-999;
           fYmax=-999;
           fZmax=-999;
           fXmin=-999;
           fYmin=-999;
           fZmin=-999;
       }
       filtpart       =  p.get< int >("filterpart");
       
     return;
   }