Public Member Functions
	Trigger ()

	~Trigger ()

void	AddMuon (simb::MCParticle particle)

int	GetTriggerId ()

double	GetTriggerOffsetX ()

double	GetTriggerOffsetY ()

double	GetTriggerOffsetZ ()

double	GetTriggerOffsetT ()

TLorentzVector	GetTriggerPos ()

double	GetTriggerPosX ()

double	GetTriggerPosY ()

double	GetTriggerPosZ ()

double	GetTriggerPosT ()

TLorentzVector	GetTriggerMom ()

double	GetTriggerMomX ()

double	GetTriggerMomY ()

double	GetTriggerMomZ ()

double	GetTriggerMomT ()

void	MakeTrigger (CLHEP::HepRandomEngine &engine)

double	GetPhi (const double py, const double pz)

void	GetMatrix (double theta, double phi, double(*p_R)[3][3])

void	DoRotation (double urv[], double dir[], double theta, double phi)

bool	Intersect (const double x0[], const double dx[], const double bounds[])

void	ProjectToBoxEdge (const double xyz[], const double indxyz[], const double xlo, const double xhi, const double ylo, const double yhi, const double zlo, const double zhi, double xyzout[])

void	GetTPCSize (double tpc[])

void	GetCryoSize (double cryo[])

void	SetCryoBuffer (std::vector< double > buffer)

void	SetTPCBuffer (std::vector< double > buffer)

Private Attributes
art::ServiceHandle< geo::Geometry >	geom

std::vector< simb::MCParticle >	fMuonList

simb::MCParticle	fTriggerMu

TLorentzVector	fTriggerPos

double	fTriggerID = -999

double	fTriggerPosX

double	fTriggerPosY

double	fTriggerPosZ

double	fTriggerPosT

double	fTriggerOffsetX

double	fTriggerOffsetY

double	fTriggerOffsetZ

double	fTriggerOffsetT

std::vector< double >	fTPCBuffer

std::vector< double >	fCryoBuffer

Detailed Description

Definition at line 124 of file Gen311_module.cc.

Constructor & Destructor Documentation

evgendp::Trigger::Trigger ( )

Definition at line 822 of file Gen311_module.cc.

                        {
   fMuonList.clear();
 }

evgendp::Trigger::~Trigger ( )

Definition at line 826 of file Gen311_module.cc.

826 {}

Member Function Documentation

void evgendp::Trigger::AddMuon ( simb::MCParticle particle )

inline

Definition at line 131 of file Gen311_module.cc.

                                                {
           fMuonList.push_back(particle); return;
         };

void evgendp::Trigger::DoRotation	(	double	urv[],
		double	dir[],
		double	theta,
		double	phi
	)

Definition at line 864 of file Gen311_module.cc.

                                                                                      {
 
   //build the rotation matrix;
   double R[3][3];
   double (*p_R)[3][3]=&R;
   GetMatrix(theta, phi, p_R);
 
   for(int i=0; i<3; i++){
     for(int j=0; j<3; j++){
       dir[i] += (*p_R)[i][j]*urv[j];
       //cout << i << " " << j << " " << urv[j] << " " << (*p_R)[i][j] << " " << dir[i] << endl;
     }
   }
 
   return;
 }

void evgendp::Trigger::GetCryoSize ( double cryo[] )

Definition at line 944 of file Gen311_module.cc.

                                               {
 
   //;
 
   double dummy[6] = {0};
   for(unsigned int c = 0; c < geom->Ncryostats(); ++c){
     geom->CryostatBoundaries(dummy, c);
   }
 
   for(int i=0;i<6;i++){cryo[i] = dummy[i];}
 
   return;
 }

void evgendp::Trigger::GetMatrix	(	double	theta,
		double	phi,
		double(*)	p_R[3][3]
	)

Definition at line 849 of file Gen311_module.cc.

                                                                              {
 
   (*p_R)[0][0]= cos(phi)*cos(theta);
   (*p_R)[0][1]= cos(phi)*sin(theta);
   (*p_R)[0][2]= -sin(phi);
   (*p_R)[1][0]= -sin(theta);
   (*p_R)[1][1]= cos(theta);
   (*p_R)[1][2]= 0;
   (*p_R)[2][0]= sin(phi)*cos(theta);
   (*p_R)[2][1]= sin(phi)*sin(theta);
   (*p_R)[2][2]= cos(phi);
 
   return;
 }

double evgendp::Trigger::GetPhi	(	const double	py,
		const double	pz
	)

Definition at line 828 of file Gen311_module.cc.

                                                                {
 
     if( pz !=0 ){
       if( pz > 0 ){
         return atan( py/pz );
       }
       else if( atan( py/pz ) < 0 ){
         return atan( py/pz )+3.1415926;
       }
       else{
         return atan( py/pz )-3.1415926;
       }
     }
     else if( py > 0 ){
       return 3.1415926/2.;
     }
     else{
       return -3.1415926/2.;
     }
 }

void evgendp::Trigger::GetTPCSize ( double tpc[] )

Definition at line 915 of file Gen311_module.cc.

                                             {
 
   //art::ServiceHandle<geo::Geometry> geom;
 
   double minx=0, miny=0, minz=0;
   double maxx=0, maxy=0, maxz=0;
   for (size_t c = 0; c < geom->Ncryostats(); c++){
          const geo::CryostatGeo& cryostat = geom->Cryostat(c);
          for (size_t t = 0; t < cryostat.NTPC(); t++){
              const geo::TPCGeo& tpcg = cryostat.TPC(t);
              if (tpcg.MinX() < minx) minx = tpcg.MinX();
              if (tpcg.MaxX() > maxx) maxx = tpcg.MaxX();
              if (tpcg.MinY() < miny) miny = tpcg.MinY();
              if (tpcg.MaxY() > maxy) maxy = tpcg.MaxY();
              if (tpcg.MinZ() < minz) minz = tpcg.MinZ();
              if (tpcg.MaxZ() > maxz) maxz = tpcg.MaxZ();
          }
   }
 
   tpc[0] = minx;
   tpc[1] = maxx;
   tpc[2] = miny;
   tpc[3] = maxy;
   tpc[4] = minz;
   tpc[5] = maxz;
 
   return;
 }

int evgendp::Trigger::GetTriggerId ( )

inline

Definition at line 136 of file Gen311_module.cc.

                           {
           return fTriggerID;
         };

TLorentzVector evgendp::Trigger::GetTriggerMom ( )

inline

Definition at line 176 of file Gen311_module.cc.

                                       {
           return fTriggerMu.Momentum();
         };

double evgendp::Trigger::GetTriggerMomT ( )

inline

Definition at line 192 of file Gen311_module.cc.

                                {
           return fTriggerMu.Momentum().T();
         };

double evgendp::Trigger::GetTriggerMomX ( )

inline

Definition at line 180 of file Gen311_module.cc.

                                {
           return fTriggerMu.Momentum().X();
         };

double evgendp::Trigger::GetTriggerMomY ( )

inline

Definition at line 184 of file Gen311_module.cc.

                                {
           return fTriggerMu.Momentum().Y();
         };

double evgendp::Trigger::GetTriggerMomZ ( )

inline

Definition at line 188 of file Gen311_module.cc.

                                {
           return fTriggerMu.Momentum().Z();
         };

double evgendp::Trigger::GetTriggerOffsetT ( )

inline

Definition at line 152 of file Gen311_module.cc.

                                   {
           return fTriggerOffsetT;
         };

double evgendp::Trigger::GetTriggerOffsetX ( )

inline

Definition at line 140 of file Gen311_module.cc.

                                   {
           return fTriggerOffsetX;
         };

double evgendp::Trigger::GetTriggerOffsetY ( )

inline

Definition at line 144 of file Gen311_module.cc.

                                   {
           return fTriggerOffsetY;
         };

double evgendp::Trigger::GetTriggerOffsetZ ( )

inline

Definition at line 148 of file Gen311_module.cc.

                                   {
           return fTriggerOffsetZ;
         };

TLorentzVector evgendp::Trigger::GetTriggerPos ( )

inline

Definition at line 156 of file Gen311_module.cc.

                                       {
           return fTriggerPos;
         };

double evgendp::Trigger::GetTriggerPosT ( )

inline

Definition at line 172 of file Gen311_module.cc.

                                {
           return fTriggerPosT;
         };

double evgendp::Trigger::GetTriggerPosX ( )

inline

Definition at line 160 of file Gen311_module.cc.

                                {
           return fTriggerPosX;
         };

double evgendp::Trigger::GetTriggerPosY ( )

inline

Definition at line 164 of file Gen311_module.cc.

                                {
           return fTriggerPosY;
         };

double evgendp::Trigger::GetTriggerPosZ ( )

inline

Definition at line 168 of file Gen311_module.cc.

                                {
           return fTriggerPosZ;
         };

bool evgendp::Trigger::Intersect	(	const double	x0[],
		const double	dx[],
		const double	bounds[]
	)

Definition at line 881 of file Gen311_module.cc.

                                                                                          {
   //check if the particle intercept the tpc
 
   bool intersects_tpc = false;
       for (int bnd=0; bnd!=6; ++bnd) {
         if (bnd<2) {
           double p2[3] = {bounds[bnd],  x0[1] + (dx[1]/dx[0])*(bounds[bnd] - x0[0]), x0[2] + (dx[2]/dx[0])*(bounds[bnd] - x0[0])};
           if ( p2[1] >= bounds[2] && p2[1] <= bounds[3] &&
                p2[2] >= bounds[4] && p2[2] <= bounds[5] ) {
             intersects_tpc = true;
             break;
           }
         }
         else if (bnd>=2 && bnd<4) {
           double p2[3] = {x0[0] + (dx[0]/dx[1])*(bounds[bnd] - x0[1]), bounds[bnd], x0[2] + (dx[2]/dx[1])*(bounds[bnd] - x0[1])};
           if ( p2[0] >= bounds[0] && p2[0] <= bounds[1] &&
                p2[2] >= bounds[4] && p2[2] <= bounds[5] ) {
             intersects_tpc = true;
       break;
           }
         }
         else if (bnd>=4) {
           double p2[3] = {x0[0] + (dx[0]/dx[2])*(bounds[bnd] - x0[2]), x0[1] + (dx[1]/dx[2])*(bounds[bnd] - x0[2]), bounds[bnd]};
           if ( p2[0] >= bounds[0] && p2[0] <= bounds[1] &&
                p2[1] >= bounds[2] && p2[1] <= bounds[3] ) {
             intersects_tpc = true;
       break;
           }
         }
       }
 
       return intersects_tpc;
 }

void evgendp::Trigger::MakeTrigger ( CLHEP::HepRandomEngine & engine )

Definition at line 997 of file Gen311_module.cc.

                                                              {
 
     //get the coordinates of the tpc and an active volume arount it
     double tpc[6] ={0.}; this->GetTPCSize(tpc);
     for(int i=0; i<6; i++){ tpc[i] += fTPCBuffer[i]; }
 
     //get the coordinates of the cryo + a buffer around it
     double cryo[6] ={0.}; this->GetCryoSize(cryo);
     for(int i=0; i<6; i++){ cryo[i] += fCryoBuffer[i]; }
 
     CLHEP::RandFlat flat(engine);
 
     //choose a random muon
     if( fMuonList.size() == 0){
       mf::LogInfo("Gen311") << "Trigger muon not found! Only Background";
       return;
     }
 
     fTriggerMu = fMuonList.at( (int)flat()*fMuonList.size() );
 
     double px=0,py=0,pz=0;
     double p=0, theta=0, phi=0;
 
     px = fTriggerMu.Momentum().X();
     py = fTriggerMu.Momentum().Y();
     pz = fTriggerMu.Momentum().Z();
 
 
     p = sqrt(px*px+py*py+pz*pz);
     theta = acos(px/p);
     phi = this->GetPhi( py, pz );
 
     double xyz[3] ={0.};
     xyz[0]=(tpc[5]-tpc[4])/2-tpc[5];
     xyz[1]=(tpc[1]-tpc[0])/2-tpc[1];
     xyz[2]=(tpc[3]-tpc[2])/2-tpc[2];
 
     //define direction
     double cos_x=cos(theta);
     double cos_y=sin(phi)*sin(theta);
     double cos_z=cos(phi)*sin(theta);
     double dxyz[3] = { cos_z, cos_x, cos_y };
 
     //get cr impact size along u and v
     double length[2]={0.}; //0 along u, 1 along v
     double corner[3]={0.};
     double rdm_start[3] ={0.};
     double xyzo[3]={0.};
 
     //loop over the tpc boundaries and project them on u and v
     for(int i=0; i<2; i++){
       for(int j=2; j<4; j++){
         for(int k=4; k<6; k++){
 
             double dot=0;
             double proj[3]={0.}; //projection on the plane
             double u=0; double v=0;
 
             //coordinates of the tpc corner
             corner[0] = tpc[k];
             corner[1] = tpc[i];
             corner[2] = tpc[j];
 
             //now get the dot product between the corner vector and the plane versor
             for( int p=0; p<3; p++ ){  dot+=(corner[p]-xyz[p])*dxyz[p]; }
             for( int p=0; p<3; p++ ){  proj[p] = corner[p]-(dot*dxyz[p]); }
 
             proj[0] -= (tpc[5]-tpc[4])/2; //NB: not really elegant...makes things work though
 
             //build the rotation matrix;
             double R[3][3];
             double (*p_R)[3][3]=&R;
             this->GetMatrix(theta, phi, p_R);
 
             //project proj along u and pick the longest
             for( int p=0; p<3; p++ ){ u += (*p_R)[p][0]*proj[p]; }
             if(length[0] < u){ length[0] = u; }
 
             //project proj along v and pick the longest projection
             for( int p=0; p<3; p++ ){ v += (*p_R)[p][2]*proj[p]; }
             if(length[1] < v){ length[1] = v; }
 
         }//end for k
       }//end for j
     } //end for i
 
     bool is_inside=false;
     int iteration = 20; //NB: It is harcoded
 
     while(!is_inside && iteration > 0){
 
         //We extract a random start over the plane u, v
         double u = flat()*2*length[0]-length[0];
         double v = flat()*2*length[1]-length[1];;
 
         //now rotate along the r direction
         double urv[3] ={ u,  0,  v };
         double dir[3]={ 0., 0., 0. };
         this->DoRotation( urv, dir, theta, phi );
 
         //parse the correct direction of our reference system
         rdm_start[0] = (double)dir[1];
         rdm_start[1] = (double)dir[2];
         rdm_start[2] = (double)dir[0]+(tpc[5]-tpc[4])/2;
 
         double dx[3]={px,py,pz};
         this->ProjectToBoxEdge(rdm_start, dx, cryo[0], cryo[1], cryo[2], cryo[3], cryo[4], cryo[5], xyzo);
 
         if( this->Intersect(xyzo, dx, tpc) ){
           is_inside=true;
           //mf::LogInfo("Gen311")<< "======> " << iteration;
         }else{
           //mf::LogInfo("Gen311")<< "======> " << iteration;
           iteration--;
         }
     }//end while
 
 
     if(iteration==0){
       mf::LogInfo("Gen311") << "Trigger muon not found in 20 iterations. Only Background";
     }
 
     /*
 
     double rdm_start[3] ={0.};
     rdm_start[0] = 0;
     rdm_start[1] = flat()*(tpc[3]-tpc[2])-(tpc[3]-tpc[2])/2;
     rdm_start[2] = flat()*(tpc[5]-tpc[4])-(tpc[5]-tpc[4])/2;
 
     */
 
     fTriggerID = fTriggerMu.TrackId();
     fTriggerPosX = rdm_start[0];
     fTriggerPosY = rdm_start[1];
     fTriggerPosZ = rdm_start[2];
     fTriggerPosT = 0;
     fTriggerPos.SetXYZT(rdm_start[0],rdm_start[1],rdm_start[2],0);
     fTriggerOffsetX = rdm_start[0] - fTriggerMu.Position().X();
     fTriggerOffsetY = rdm_start[1] - fTriggerMu.Position().Y();
     fTriggerOffsetZ = rdm_start[2] - fTriggerMu.Position().Z();
     fTriggerOffsetT = 0 - fTriggerMu.Position().T();
 
     return;
   }//end make trigger

void evgendp::Trigger::ProjectToBoxEdge	(	const double	xyz[],
		const double	indxyz[],
		const double	xlo,
		const double	xhi,
		const double	ylo,
		const double	yhi,
		const double	zlo,
		const double	zhi,
		double	xyzout[]
	)

Definition at line 958 of file Gen311_module.cc.

                                                           {
 
 
   //we want to project backwards, so take mirror of momentum
   const double dxyz[3]={-indxyz[0],-indxyz[1],-indxyz[2]};
 
   // Compute the distances to the x/y/z walls
   double dx = 99.E99;
   double dy = 99.E99;
   double dz = 99.E99;
   if      (dxyz[0] > 0.0) { dx = (xhi-xyz[0])/dxyz[0]; }
   else if (dxyz[0] < 0.0) { dx = (xlo-xyz[0])/dxyz[0]; }
   if      (dxyz[1] > 0.0) { dy = (yhi-xyz[1])/dxyz[1]; }
   else if (dxyz[1] < 0.0) { dy = (ylo-xyz[1])/dxyz[1]; }
   if      (dxyz[2] > 0.0) { dz = (zhi-xyz[2])/dxyz[2]; }
   else if (dxyz[2] < 0.0) { dz = (zlo-xyz[2])/dxyz[2]; }
 
 
   // Choose the shortest distance
   double d = 0.0;
   if      (dx < dy && dx < dz) d = dx;
   else if (dy < dz && dy < dx) d = dy;
   else if (dz < dx && dz < dy) d = dz;
 
   // Make the step
   for (int i = 0; i < 3; ++i) {
     xyzout[i] = xyz[i] + dxyz[i]*d;
   }
 
 }

void evgendp::Trigger::SetCryoBuffer ( std::vector< double > buffer )

inline

Definition at line 216 of file Gen311_module.cc.

                                                       {
           for(int i=0; i<6; i++){ fCryoBuffer.assign( buffer.begin(), buffer.end() ); }
         }

void evgendp::Trigger::SetTPCBuffer ( std::vector< double > buffer )

inline

Definition at line 219 of file Gen311_module.cc.

                                                      {
           for(int i=0; i<6; i++){ fTPCBuffer.assign( buffer.begin(), buffer.end() ); }
         }

Member Data Documentation

std::vector<double> evgendp::Trigger::fCryoBuffer

private

Definition at line 241 of file Gen311_module.cc.

std::vector<simb::MCParticle> evgendp::Trigger::fMuonList

private

Definition at line 228 of file Gen311_module.cc.

std::vector<double> evgendp::Trigger::fTPCBuffer

private

Definition at line 240 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerID = -999

private

Definition at line 231 of file Gen311_module.cc.

simb::MCParticle evgendp::Trigger::fTriggerMu

private

Definition at line 229 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerOffsetT

private

Definition at line 239 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerOffsetX

private

Definition at line 236 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerOffsetY

private

Definition at line 237 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerOffsetZ

private

Definition at line 238 of file Gen311_module.cc.

TLorentzVector evgendp::Trigger::fTriggerPos

private

Definition at line 230 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerPosT

private

Definition at line 235 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerPosX

private

Definition at line 232 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerPosY

private

Definition at line 233 of file Gen311_module.cc.

double evgendp::Trigger::fTriggerPosZ

private

Definition at line 234 of file Gen311_module.cc.

art::ServiceHandle<geo::Geometry> evgendp::Trigger::geom

private

Definition at line 226 of file Gen311_module.cc.

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

dunesim/dunesim/EventGenerator/CORSIKAprotodunedp/Gen311_module.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation