icarus::crt Namespace Reference

Namespaces
	CRTDataFlags
	icarus::crt::CRTData flag mnemonics

Classes
class	CRTAutoVeto
class	CRTCalibrationAnalysis
class	CRTDataAnalysis
class	AnaProducer
class	CrtCal
class	CrtCalAnalyzer
class	CrtCalTree
class	CRTEventProducer
class	CrtNoiseMonTool
class	CRTPreProcessTree
class	CRTRawTree
class	CRTTiming
class	CRTDetSim
class	CRTGeometryHelper
class	CrtOpHitMatchAnalysis
class	CRTPMTMatchingAna
class	CRTSimAnalysis
	Example analyzer. More...
class	CRTSimHitProducer
class	CRTTrackProducer
struct	CRTavehit
class	CRTTrueHitProducer
class	CRTTruthMatchAnalysis
class	CRTBackTracker
class	CRTCommonUtils
class	CRTDetSimAlg
class	CRTHitRecoAlg
struct	FEB_delay
class	CRTTrackRecoAlg
class	CRTTrueHitRecoAlg
struct	CRTHit
struct	CRTTrack
struct	CRTTzero
struct	CRTData
	Single hit (self trigger) of a CRT board. More...

Typedefs
typedef int	feb_index
typedef std::map< feb_index, FEB_delay >	CRT_delay_map

Functions
void	vmanip (vector< float > v, float *ave, float *rms)
CRTavehit	fillme (uint32_t i, uint16_t j, int32_t k, uint16_t l, float a, float b, float c, float d, float e, float f, float g, int p, string t)
CRTavehit	copyme (sbn::crt::CRTHit myhit)
sbn::crt::CRTTrack	shcut (CRTavehit ppA, CRTavehit ppb, uint32_t time0s, uint16_t terr)
bool	TimeOrderCRTData (std::pair< ChanData, AuxDetIDE > crtdat1, std::pair< ChanData, AuxDetIDE > crtdat2)
bool	sortbytime (const pair< int, ULong64_t > &a, const pair< int, ULong64_t > &b)
ULong64_t	GetMode (std::vector< std::pair< int, ULong64_t >> vector)
CRT_delay_map	LoadFEBMap ()

Variables
struct icarus::crt::CRTavehit	tempah

Typedef Documentation

typedef std::map<feb_index, FEB_delay> icarus::crt::CRT_delay_map

Definition at line 215 of file icaruscode/icaruscode/CRT/CRTUtils/CRTHitRecoAlg.h.

typedef int icarus::crt::feb_index

Definition at line 214 of file icaruscode/icaruscode/CRT/CRTUtils/CRTHitRecoAlg.h.

Function Documentation

CRTavehit icarus::crt::copyme

(

sbn::crt::CRTHit

myhit

)

Definition at line 435 of file icaruscode/icaruscode/CRT/CRTTrackProducer_module.cc.

  {
    CRTavehit h;
    h.ts0_ns     = myhit.ts0_ns;
    h.ts0_ns_err = 0;
    h.ts1_ns     = myhit.ts1_ns;; 
    h.ts1_ns_err = 0;       
    h.x_pos      = myhit.x_pos;
    h.x_err      = myhit.x_err;
    h.y_pos      = myhit.y_pos;
    h.y_err      = myhit.y_err;
    h.z_pos      = myhit.z_pos;
    h.z_err      = myhit.z_err;
    h.pe         = myhit.peshit;
    h.plane      = myhit.plane;
    h.tagger     = myhit.tagger;
    return(h);
  } // copyme()

CRTavehit icarus::crt::fillme	(	uint32_t	i,
		uint16_t	j,
		int32_t	k,
		uint16_t	l,
		float	a,
		float	b,
		float	c,
		float	d,
		float	e,
		float	f,
		float	g,
		int	p,
		string	t
	)

Definition at line 412 of file icaruscode/icaruscode/CRT/CRTTrackProducer_module.cc.

  {
    CRTavehit h;
    h.ts0_ns     = ts0_ns;
    h.ts0_ns_err = ts0_ns_err;
    h.ts1_ns     = ts1_ns; 
    h.ts1_ns_err = ts1_ns_err;                                                        
    
    h.x_pos      = x_pos;
    h.x_err      = x_err;
    h.y_pos      = y_pos;
    h.y_err      = y_err;
    h.z_pos      = z_pos;
    h.z_err      = z_err;
    h.pe         = pe;
    h.plane      = plane;
    h.tagger     = tagger;
    return(h);
  } // fillme()

ULong64_t icarus::crt::GetMode

(

std::vector< std::pair< int, ULong64_t >>

vector

)

Definition at line 1085 of file icaruscode/icaruscode/CRT/CRTUtils/CRTHitRecoAlg.cc.

                                                                      {

        sort(vector.begin(), vector.end(), icarus::crt::sortbytime);

        int modecounter = 0;
        int isnewmodecounter = 0;
        ULong64_t Mode = 0;
        ULong64_t isnewMode = 0;
        bool isFirst = true;
        for (auto i : vector) {
                if (!isFirst) {
                        if (i.second == Mode) modecounter++;
                        else if (i.second !=isnewMode) {
                                isnewMode = i.second;
                                isnewmodecounter = 1;
                        }
                        else if (i.second == isnewMode) {
                                isnewmodecounter++;
                                if (isnewmodecounter > modecounter) {
                                        Mode = isnewMode;
                                        modecounter = isnewmodecounter;
                                }
                        }
                }
                else {
                        isFirst = false;
                        Mode = i.second;
                        modecounter++;
                }
        }
        return Mode;
}

icarus::crt::CRT_delay_map icarus::crt::LoadFEBMap ( )

inline

Definition at line 221 of file icaruscode/icaruscode/CRT/CRTUtils/CRTHitRecoAlg.h.

                                                    {

        CRT_delay_map FEBs;

        FEBs={{271, {81, 198, 271, 283ull, 2000309ull}},
                {270, {119, 197, 270, 298ull, 2000324ull}},
                {269, {87, 196, 269, 313ull, 2000339ull}},
                {268, {92, 195, 268, 329ull, 2000355ull}},
                {267, {180, 194, 267, 344ull, 2000370ull}},
                {266, {97, 193, 266, 359ull, 2000385ull}},
                {265, {174, 192, 265, 374ull, 2000400ull}},
                {251, {238, 178, 251, 390ull, 2000416ull}},
                {237, {234, 164, 237, 405ull, 2000431ull}},
                {297, {189, 224, 297, 420ull, 2000446ull}},
                {296, {190, 223, 296, 436ull, 2000462ull}},
                {295, {80, 222, 295, 451ull, 2000477ull}},
                {294, {162, 221, 294, 466ull, 2000492ull}},
                {293, {64, 220, 293, 482ull, 2000508ull}},
                {255, {172, 182, 255, 298ull, 2000324ull}},
                {254, {114, 181, 254, 313ull, 2000339ull}},
                {253, {100, 180, 253, 328ull, 2000355ull}},
                {252, {150, 179, 252, 344ull, 2000370ull}},
                {238, {176, 165, 238, 359ull, 2000385ull}},
                {224, {67, 151, 224, 374ull, 2000400ull}},
                {223, {138, 150, 223, 390ull, 2000416ull}},
                {209, {170, 136, 209, 405ull, 2000431ull}},
                {195, {101, 122, 195, 420ull, 2000446ull}},
                {181, {142, 108, 181, 435ull, 2000462ull}},
                {279, {139, 206, 279, 451ull, 2000477ull}},
                {280, {185, 207, 280, 466ull, 2000492ull}},
                {182, {6, 109, 182, 481ull, 2000508ull}},
                {196, {177, 123, 196, 497ull, 2000523ull}},
                {210, {61, 137, 210, 512ull, 2000538ull}},
                {256, {123, 183, 256, 298ull, 2000325ull}},
                {242, {116, 169, 242, 314ull, 2000340ull}},
                {241, {104, 168, 241, 329ull, 2000355ull}},
                {240, {91, 167, 240, 344ull, 2000371ull}},
                {239, {88, 166, 239, 360ull, 2000386ull}},
                {225, {120, 152, 225, 375ull, 2000401ull}},
                {211, {132, 138, 211, 390ull, 2000417ull}},
                {197, {95, 124, 197, 405ull, 2000432ull}},
                {183, {232, 110, 183, 421ull, 2000447ull}},
                {281, {165, 208, 281, 436ull, 2000463ull}},
                {282, {148, 209, 282, 451ull, 2000478ull}},
                {184, {237, 111, 184, 467ull, 2000493ull}},
                {198, {102, 125, 198, 482ull, 2000508ull}},
                {212, {94, 139, 212, 497ull, 2000524ull}},
                {226, {130, 153, 226, 513ull, 2000539ull}},
                {257, {181, 184, 257, 284ull, 2000310ull}},
                {243, {124, 170, 243, 299ull, 2000325ull}},
                {229, {152, 156, 229, 314ull, 2000341ull}},
                {228, {98, 155, 228, 329ull, 2000356ull}},
                {227, {173, 154, 227, 345ull, 2000371ull}},
                {213, {169, 140, 213, 360ull, 2000387ull}},
                {199, {144, 126, 199, 375ull, 2000402ull}},
                {185, {239, 112, 185, 391ull, 2000417ull}},
                {283, {147, 210, 283, 306ull, 2000433ull}},
                {284, {105, 211, 284, 421ull, 2000448ull}},
                {186, {231, 114, 186, 437ull, 2000463ull}},
                {200, {117, 127, 200, 452ull, 2000478ull}},
                {214, {126, 141, 214, 467ull, 2000494ull}},
                {215, {90, 142, 215, 482ull, 2000509ull}},
                {201, {183, 128, 201, 498ull, 2000524ull}},
                {187, {241, 114, 187, 513ull, 2000540ull}},
                {285, {113, 212, 285, 528ull, 2000555ull}},
                {258, {233, 185, 258, 283ull, 2000310ull}},
                {244, {164, 171, 244, 299ull, 2000325ull}},
                {230, {161, 157, 230, 314ull, 2000341ull}},
                {231, {203, 158, 231, 329ull, 2000356ull}},
                {232, {122, 159, 232, 345ull, 2000371ull}},
                {218, {2, 145, 218, 360ull, 2000387ull}},
                {204, {112, 131, 204, 375ull, 2000402ull}},
                {190, {62, 117, 190, 391ull, 2000417ull}},
                {288, {133, 215, 288, 406ull, 2000432ull}},
                {287, {168, 214, 287, 421ull, 2000448ull}},
                {189, {182, 116, 189, 436ull, 2000463ull}},
                {203, {107, 130, 203, 452ull, 2000478ull}},
                {217, {252, 144, 217, 467ull, 2000494ull}},
                {216, {141, 143, 216, 482ull, 2000509ull}},
                {202, {160, 129, 202, 498ull, 2000524ull}},
                {188, {137, 115, 188, 513ull, 2000540ull}},
                {286, {179, 213, 286, 528ull, 2000555ull}},
                {259, {66, 186, 259, 298ull, 2000325ull}},
                {245, {247, 172, 245, 314ull, 2000340ull}},
                {246, {198, 173, 246, 329ull, 2000356ull}},
                {247, {243, 174, 247, 344ull, 2000371ull}},
                {248, {72, 175, 248, 360ull, 2000386ull}},
                {234, {250, 161, 234, 375ull, 2000401ull}},
                {220, {249, 147, 220, 390ull, 2000417ull}},
                {206, {248, 133, 206, 405ull, 2000432ull}},
                {192, {60, 119, 192, 421ull, 2000447ull}},
                {290, {145, 217, 290, 436ull, 2000463ull}},
                {289, {110, 216, 289, 451ull, 2000478ull}},
                {191, {59, 118, 191, 467ull, 2000493ull}},
                {205, {202, 132, 205, 482ull, 2000509ull}},
                {219, {135, 146, 219, 497ull, 2000524ull}},
                {233, {246, 160, 233, 513ull, 2000539ull}},
                {260, {253, 187, 260, 342ull, 2000369ull}},
                {261, {245, 188, 261, 358ull, 2000384ull}},
                {262, {65, 189, 262, 373ull, 2000400ull}},
                {263, {57, 190, 263, 388ull, 2000415ull}},
                {249, {63, 176, 249, 404ull, 2000430ull}}, 
                {250, {251, 177, 250, 419ull, 2000445ull}},
                {236, {70, 163, 236, 434ull, 2000461ull}},
                {222, {155, 149, 222, 449ull, 2000476ull}},
                {208, {154, 135, 208, 465ull, 2000491ull}},
                {194, {85, 121, 194, 480ull, 2000507ull}},
                {292, {134, 219, 292, 495ull, 2000522ull}},
                {291, {129, 218, 291, 511ull, 2000537ull}},
                {193, {115, 120, 193, 526ull, 2000553ull}},
                {207, {204, 134, 207, 541ull, 2000568ull}},
                {221, {244, 148, 221, 557ull, 2000583ull}},
                {235, {82, 162, 235, 572ull, 2000598ull}},
                {272, {186, 199, 272, 284ull, 2000310ull}},
                {273, {83, 200, 273, 299ull, 2000326ull}},
                {274, {254, 201, 274, 314ull, 2000341ull}},
                {275, {166, 202, 275, 330ull, 2000356ull}},
                {276, {178, 203, 276, 345ull, 2000371ull}},
                {277, {136, 204, 277, 360ull, 2000387ull}},
                {278, {184, 205, 278, 375ull, 2000402ull}},
                {264, {187, 191, 264, 391ull, 2000417ull}},
                {304, {240, 231, 304, 406ull, 2000433ull}},
                {303, {242, 230, 303, 421ull, 2000448ull}},
                {302, {188, 229, 302, 437ull, 2000463ull}},
                {301, {58, 228, 301, 452ull, 2000479ull}},
                {300, {143, 227, 300, 467ull, 2000494ull}},
                {299, {235, 226, 299, 483ull, 2000509ull}}};
        return FEBs;
}

sbn::crt::CRTTrack icarus::crt::shcut	(	CRTavehit	ppA,
		CRTavehit	ppb,
		uint32_t	time0s,
		uint16_t	terr
	)

Definition at line 455 of file icaruscode/icaruscode/CRT/CRTTrackProducer_module.cc.

  {
    sbn::crt::CRTTrack newtr;
    newtr.ts0_s         = time0s;
    newtr.ts0_s_err     = terr;
    newtr.ts0_ns_h1     = ppA.ts0_ns;
    newtr.ts0_ns_err_h1 = ppA.ts0_ns_err;
    newtr.ts0_ns_h2     = ppB.ts0_ns;
    newtr.ts0_ns_err_h2 = ppB.ts0_ns_err;
    newtr.ts0_ns        = (uint32_t)(0.5*(ppA.ts0_ns+ppB.ts0_ns));
    newtr.ts0_ns_err    = (uint16_t)(0.5*sqrt(ppA.ts0_ns_err*ppA.ts0_ns_err+ppB.ts0_ns_err*ppB.ts0_ns_err));
    newtr.ts1_ns        = (int32_t)(0.5*(ppA.ts1_ns+ppB.ts1_ns));
    newtr.ts1_ns_err    = (uint16_t)(0.5*sqrt(ppA.ts0_ns_err*ppA.ts0_ns_err+ppB.ts0_ns_err*ppB.ts0_ns_err));
    newtr.peshit        = ppA.pe+ppB.pe;
    newtr.x1_pos        = ppA.x_pos;
    newtr.x1_err        = ppA.x_err;
    newtr.y1_pos        = ppA.y_pos;
    newtr.y1_err        = ppA.y_err;
    newtr.z1_pos        = ppA.z_pos;
    newtr.z1_err        = ppA.z_err;
    newtr.x2_pos        = ppB.x_pos;
    newtr.x2_err        = ppB.x_err;
    newtr.y2_pos        = ppB.y_pos;
    newtr.y2_err        = ppB.y_err;
    newtr.z2_pos        = ppB.z_pos;
    newtr.z2_err        = ppB.z_err;
    float deltax        = ppA.x_pos-ppB.x_pos;
    float deltay        = ppA.y_pos-ppB.y_pos;
    float deltaz        = ppA.z_pos-ppB.z_pos;
    newtr.length        = sqrt(deltax*deltax+deltay*deltay+deltaz*deltaz);
    newtr.thetaxy       = atan2(deltax,deltay);
    newtr.phizy         = atan2(deltaz,deltay);
    newtr.plane1        = ppA.plane;
    newtr.plane2        = ppB.plane;

    return(newtr);

  } // shcut()

bool icarus::crt::sortbytime ( const pair< int, ULong64_t > &  a, const pair< int, ULong64_t > &  b  )

inline

Definition at line 68 of file icaruscode/icaruscode/CRT/CRTUtils/CRTHitRecoAlg.h.

  {
    return (a.second < b.second);
  } 

bool icarus::crt::TimeOrderCRTData	(	std::pair< ChanData, AuxDetIDE >	crtdat1,
		std::pair< ChanData, AuxDetIDE >	crtdat2
	)

Definition at line 9 of file CRTDetSimAlg.cc.

                                                                {
        return ( crtdat1.first.ts < crtdat2.first.ts );
    }//TimeOrderCRTData()

void icarus::crt::vmanip	(	std::vector< float >	v,
		float *	ave,
		float *	rms
	)

Definition at line 394 of file icaruscode/icaruscode/CRT/CRTTrackProducer_module.cc.

  {
    *ave=0.0; *rms =0.0;
    if (v.size()>0) {
      //  find the mean and *rms of all the vector elements
      double sum = std::accumulate(v.begin(), v.end(), 0.0);
      double mean = sum / v.size();
      *ave=mean;
      
      if (v.size()>1) {
      double sq_sum = std::inner_product(v.begin(), v.end(), v.begin(), 0.0);
      double stdev = std::sqrt(sq_sum / v.size() - mean * mean);
      *rms=stdev;
      }
    }
  } // vmanip()

Variable Documentation

struct icarus::crt::CRTavehit icarus::crt::tempah