AltDeMuxDisplay Class Reference

#include <AltDeMuxDisplay.h>

List of all members.

Public Member Functions
	AltDeMuxDisplay ()
void	Ana (const MomNavigator *mom)
bool	Draw (int ipad)
int	MaxPlane (int irun)
void	FindNearestEdge (float x, float y, int &idege, float &d)
Private Attributes
bool	amPSing
bool	amCutting
bool	amVetoing
int	_maxPlane
int	_lastRun
int	xMaxPlane
int	NdemuxedHitsU
float	demuxedHitPlaneU [20000]
float	demuxedHitStripU [20000]
float	demuxedHitQU [20000]
float	demuxedHitTU [20000]
float	QSPlane [520]
float	QPlane [520]
int	NdemuxedHitsV
float	demuxedHitPlaneV [20000]
float	demuxedHitStripV [20000]
float	demuxedHitQV [20000]
float	demuxedHitTV [20000]
int	NdemuxedPairsV
int	NdemuxedPairsU
float	demuxedPairPlaneV [20000]
float	demuxedPairPlaneU [20000]
float	demuxedPairTimeV [20000]
float	demuxedPairTimeU [20000]
int	NVeto
float	xVeto [1000]
float	yVeto [1000]
float	zVeto [1000]
std::vector< PlexSEIdAltL * >	PlanesAltListsE [500]
std::vector< PlexSEIdAltL * >	PlanesAltListsW [500]
std::vector< DeMuxedPair >	_demuxedPairs
TCanvas *	fCanvas
TPad *	fUpad
TPad *	fVpad
TPad *	fXYpad
TPad *	fZQpad
TPad *	fZTpad
TH2F *	fUview
TH2F *	fVview
TH2F *	fXYview
TH1F *	fZQview
TH2F *	fZTview
TList *	fList
TPolyLine *	pZTFit
TPolyLine *	pCut1
TPolyLine *	pCut2
TPolyLine *	pSM1
TPolyLine *	pSM2
TPolyLine *	pBoxSM1
TPolyLine *	pBoxSM2
TPolyLine *	pOctagon
TPolyLine *	pOctagonC
TArc *	pCircle

---

Constructor & Destructor Documentation

AltDeMuxDisplay::AltDeMuxDisplay (   )

Definition at line 34 of file AltDeMuxDisplay.cxx. References _lastRun, _maxPlane, amCutting, amPSing, amVetoing, DISPLAY_TYPE, fCanvas, fList, fUpad, fVpad, fXYpad, fZQpad, fZTpad, pBoxSM1, pBoxSM2, pCircle, pOctagon, pOctagonC, and xMaxPlane. { amVetoing = false; amPSing = false; amCutting = false; xMaxPlane = 0; _maxPlane = 0; _lastRun = 0; fList = new TList(); fCanvas = new TCanvas("AltDeMuxDisplay","MyAnalysis",750,900); // Create the pads; Pad takes xlo, ylo, xhi, yhi in NC fUpad = new TPad("upad","U pad",0.010,0.670,0.990,0.990); fVpad = new TPad("vpad","V pad",0.010,0.340,0.990,0.660); fXYpad = new TPad("xypad","XY pad",0.010,0.010,0.620,0.330); fZQpad = new TPad("zqpad","ZQ pad",0.640,0.010,0.990,0.165); fZTpad = new TPad("zqpad","Zt pad",0.640,0.175,0.990,0.330); // white background fUpad->SetFillColor(10); fVpad->SetFillColor(10); fXYpad->SetFillColor(10); fZQpad->SetFillColor(10); fZTpad->SetFillColor(10); // tweak the margins to maximize histogram use of pads fUpad->SetLeftMargin(0.05); fUpad->SetRightMargin(0.0001); fUpad->SetTopMargin(0.0001); fUpad->SetBottomMargin(0.055); fVpad->SetLeftMargin(0.05); fVpad->SetRightMargin(0.0001); fVpad->SetTopMargin(0.0001); fVpad->SetBottomMargin(0.055); fXYpad->SetLeftMargin(0.05); fXYpad->SetRightMargin(0.0001); fXYpad->SetTopMargin(0.0001); fXYpad->SetBottomMargin(0.055); fZQpad->SetLeftMargin(0.05); fZQpad->SetRightMargin(0.0001); fZQpad->SetTopMargin(0.0001); fZQpad->SetBottomMargin(0.055); fZTpad->SetLeftMargin(0.05); fZTpad->SetRightMargin(0.0001); fZTpad->SetTopMargin(0.0001); fZTpad->SetBottomMargin(0.055); // draw the pads and the canvas if(DISPLAY_TYPE==1){ fUpad->Draw(); fVpad->Draw(); fXYpad->Draw(); fZQpad->Draw(); fZTpad->Draw(); } fCanvas->Draw(); float x[9]; float y[9]; x[0] = 0.; x[1] = 249.; x[2] = 249.; x[3] = 0.; x[4] = 0.; y[0] = 0.; y[1] = 0.; y[2] = 192.; y[3] = 192.; y[4] = 0.; pBoxSM1 = new TPolyLine(5,x,y); pBoxSM2 = 0; x[0] = -1.65; x[1] = +1.64; x[2] = +4; x[3] = +4; x[4] = +1.65; x[5] = -1.65; x[6] = -4; x[7] = -4; x[8] = -1.65; y[0] = -4; y[1] = -4; y[2] = -1.65; y[3] = +1.65; y[4] = +4; y[5] = +4; y[6] = +1.65; y[7] = -1.65; y[8] = -4; pOctagon = new TPolyLine(9,x,y); x[0] = -1.44375; x[1] = +1.44375; x[2] = +3.5; x[3] = +3.5; x[4] = +1.44375; x[5] = -1.44375; x[6] = -3.5; x[7] = -3.5; x[8] = -1.44375; y[0] = -3.5; y[1] = -3.5; y[2] = -1.44375; y[3] = +1.44375; y[4] = +3.5; y[5] = +3.5; y[6] = +1.44375; y[7] = -1.44375; y[8] = -3.5; pOctagonC = new TPolyLine(9,x,y); pCircle = new TArc(0.,0.,3.5, 0., 360.); }

---

Member Function Documentation

void AltDeMuxDisplay::Ana (  const MomNavigator *  mom  )

Definition at line 163 of file AltDeMuxDisplay.cxx. References _demuxedPairs, amVetoing, MomNavigator::At(), UgliStripHandle::ClearFiber(), demuxedHitPlaneU, demuxedHitPlaneV, demuxedHitQU, demuxedHitQV, demuxedHitStripU, demuxedHitStripV, demuxedPairPlaneU, demuxedPairPlaneV, demuxedPairTimeU, demuxedPairTimeV, Draw(), fCanvas, fList, fUpad, fUview, fVpad, fVview, fXYview, fZQview, fZTview, PlexSEIdAltL::GetBestItem(), PlexSEIdAltL::GetBestSEId(), PlexSEIdAltL::GetBestWeight(), PlexSEIdAltL::GetCurrentSEId(), CandHandle::GetDaughterIterator(), PlexSEIdAltL::GetEnd(), MomNavigator::GetFragment(), UgliStripHandle::GetHalfLength(), PlexSEIdAltLItem::GetPE(), PlexSEIdAltL::GetPlane(), PlexSEIdAltL::GetPlaneView(), CandDigitHandle::GetPlexSEIdAltLWritable(), RawRecord::GetRawHeader(), RawDaqHeader::GetRun(), PlexSEIdAltLItem::GetSigCorr(), RawDaqSnarlHeader::GetSnarl(), PlexStripEndId::GetStrip(), UgliGeomHandle::GetStripHandle(), PlexSEIdAltLItem::GetTime(), RecMinos::GetVldContext(), UgliStripHandle::GlobalPos(), UgliStripHandle::IsValid(), PlexSEIdAltL::IsValid(), PlexSEIdAltL::IsVetoShield(), MaxPlane(), MSG, NdemuxedHitsU, NdemuxedHitsV, NdemuxedPairsU, NdemuxedPairsV, NVeto, PlanesAltListsE, PlanesAltListsW, QPlane, QSPlane, PlexSEIdAltL::SetFirst(), UgliStripHandle::WlsPigtail(), xMaxPlane, xVeto, yVeto, and zVeto. { cout << " *** AltDeMuxDisplay::Ana *** " << endl; NdemuxedHitsU=0; NdemuxedHitsV=0; NdemuxedPairsV=0; NdemuxedPairsU=0; NVeto = 0; // Find PrimaryCandidateRecord fragment in MOM. CandRecord* candrecord = dynamic_cast<CandRecord*>(mom->GetFragment("CandRecord","PrimaryCandidateRecord")); if (candrecord == 0) { MSG("MyAna", Msg::kWarning) << "No PrimaryCandidateRecord in MOM." << endl; return; } int run; int snarl; TObject* obj; for (int i=0; (obj=mom->At(i)); ++i) { const RawRecord* rr = dynamic_cast<RawRecord*>(obj); if (rr) { const RawDaqSnarlHeader* rdsh = dynamic_cast<const RawDaqSnarlHeader*>(rr->GetRawHeader()); if (rdsh) { run = rdsh->GetRun(); snarl = rdsh->GetSnarl(); } } } this->MaxPlane(run); UgliStripHandle ushE; UgliStripHandle ushW; UgliStripHandle ushV; for(int i=0;i<500;i++){ QPlane[i]=0.0; QSPlane[i]=0.0; } _demuxedPairs.erase(_demuxedPairs.begin(),_demuxedPairs.end()); for(int i=0; i<500;i++){ PlanesAltListsE[i].erase(PlanesAltListsE[i].begin(),PlanesAltListsE[i].end()); PlanesAltListsW[i].erase(PlanesAltListsW[i].begin(),PlanesAltListsW[i].end()); } // clear the primitives from the pads fUpad->Clear(""); fVpad->Clear(""); // remove old histograms if they exists if ( fUview ) { delete fUview; fUview=0; } if ( fVview ) { delete fVview; fVview=0; } if ( fXYview ) { delete fXYview; fXYview=0; } if ( fZQview ) { delete fZQview; fZQview=0; } if ( fZTview ) { delete fZTview; fZTview=0; } // create histograms that we'll use for World Coordinates fUview = new TH2F("Uview","U view",1,0.,520.,1,0.,193.); fVview = new TH2F("Vview","V view",1,0.,520.,1,0.,193.); fXYview = new TH2F("XYview","XY view",1,-7.5,7.5,1,-5.,5.); fZQview = new TH1F("ZQview","ZQ view",520,0.,520.); fZTview = new TH2F("ZTview","ZT view",1,0.,520.,1,0.,50.); // don't draw the Stats box fUview->SetStats(kFALSE); fVview->SetStats(kFALSE); fXYview->SetStats(kFALSE); fZQview->SetStats(kFALSE); fZTview->SetStats(kFALSE); // take away any association with a "directory" for these histograms fUview->SetDirectory(0); fVview->SetDirectory(0); fXYview->SetDirectory(0); fZQview->SetDirectory(0); fZTview->SetDirectory(0); PlexSEIdAltL altlist; PlexSEIdAltL* paltlist; // Find PrimaryCandidateRecord fragment in MOM. const CandDigitListHandle *canddigit = DataUtil::GetCandidate<CandDigitListHandle>(mom, "CandDigitListHandle", "canddigitlist"); if (canddigit == 0) { MSG("MyAna",Msg::kWarning) << "Failed to get CandDigitListHandle!\n"; return; } // get handle to UgliGeometry const VldContext* vldc = candrecord->GetVldContext(); UgliGeomHandle ugh(*vldc); CandDigitHandleItr cdhItr(canddigit->GetDaughterIterator()); int nplanesU=0; int nplanesV=0; while ( CandDigitHandle *cdh = cdhItr() ) { // for each digit paltlist = &(cdh->GetPlexSEIdAltLWritable()); float plane = paltlist->GetPlane(); if(plane>499)cout << " HITS in Planes : " << plane << endl; if(paltlist->IsVetoShield()){ cout << "VETO SHIELD " << endl; paltlist->SetFirst(); if(paltlist->IsValid()){ cout << " VETO : IS valid " << plane << endl; ushV = ugh.GetStripHandle(paltlist->GetCurrentSEId()); if(ushV.IsValid()){ xVeto[NVeto] = (ushV.GlobalPos(0).X()/Munits::cm)/100.; yVeto[NVeto] = (ushV.GlobalPos(0).Y()/Munits::cm)/100.; zVeto[NVeto] = (ushV.GlobalPos(0).Z()/Munits::cm)/100.; NVeto++; }else{ cout << " VETO : invalid " << plane << endl; } }else{ cout << " VETO : paltlisrt invalid " << plane << endl; } } if(amVetoing&&NVeto>0)return; if(!paltlist->IsVetoShield()){ if(paltlist->GetBestWeight()>0.99){ int iplane = static_cast<int>(plane); if(iplane>xMaxPlane)xMaxPlane = iplane; if(paltlist->GetEnd()==StripEnd::kEast){ PlanesAltListsE[iplane].push_back(paltlist); } if(paltlist->GetEnd()==StripEnd::kWest){ PlanesAltListsW[iplane].push_back(paltlist); } if(plane<500){ PlaneView::PlaneView_t kView = paltlist->GetPlaneView(); PlexStripEndId seid = paltlist->GetBestSEId(); float strip = seid.GetStrip(); float QPE = paltlist->GetBestItem().GetPE(); float Q = paltlist->GetBestItem().GetSigCorr()/60.; if(Q/QPE>5. || QPE/Q > 5.)Q=QPE; if(plane>248)plane+=20; if(iplane>248)iplane+=20; if(kView==PlaneView::kU){ demuxedHitPlaneU[NdemuxedHitsU] = plane; demuxedHitStripU[NdemuxedHitsU] = strip; demuxedHitQU[NdemuxedHitsU] = Q; NdemuxedHitsU++; if(QPlane[iplane]<0.001)nplanesU++; } if(kView==PlaneView::kV){ demuxedHitPlaneV[NdemuxedHitsV] = plane; demuxedHitStripV[NdemuxedHitsV] = strip; demuxedHitQV[NdemuxedHitsV] = Q; NdemuxedHitsV++; if(QPlane[iplane]<0.001)nplanesV++; } QSPlane[iplane] += Q*strip; QPlane[iplane] += Q; } } } } float stmin=99999.; for(int iplane = 0; iplane <500; iplane++){ if(PlanesAltListsE[iplane].size()&&PlanesAltListsW[iplane].size()){ for(unsigned int iE=0; iE<PlanesAltListsE[iplane].size();iE++ ){ for(unsigned int iW=0; iW<PlanesAltListsW[iplane].size();iW++ ){ PlexSEIdAltL* pAltE = PlanesAltListsE[iplane][iE]; PlexSEIdAltL* pAltW = PlanesAltListsW[iplane][iW]; if(pAltE->GetBestSEId().GetStrip() == pAltW->GetBestSEId().GetStrip()){ PlaneView::PlaneView_t kView = pAltE->GetPlaneView(); Double_t tE= pAltE->GetBestItem().GetTime(); Double_t tW= pAltW->GetBestItem().GetTime(); ushE = ugh.GetStripHandle(pAltE->GetBestSEId()); float clearE = ushE.ClearFiber(StripEnd::kEast); float wlsE = ushE.WlsPigtail(StripEnd::kEast); ushW = ugh.GetStripHandle(pAltW->GetBestSEId()); float clearW = ushW.ClearFiber(StripEnd::kWest); float wlsW = ushW.WlsPigtail(StripEnd::kWest); float halflength = ushW.GetHalfLength(); // float dt = (tW - tE)*1E9-(clearW-clearE)*6.0-(wlsW-wlsE)*4.6; float st = (tW + tE)*1E9-(clearW+clearE)*6.0-(wlsW+wlsE)*4.6; st = st/2.-halflength*6.0; if(st<stmin)stmin=st; int jplane = iplane; if(iplane>248)jplane+=20; if(kView==PlaneView::kU){ demuxedPairPlaneU[NdemuxedPairsU] = jplane; demuxedPairTimeU[NdemuxedPairsU] = st; NdemuxedPairsU++; } if(kView==PlaneView::kV){ demuxedPairPlaneV[NdemuxedPairsV] = jplane; demuxedPairTimeV[NdemuxedPairsV] = st; NdemuxedPairsV++; } } } } } } for(int iU=0;iU<NdemuxedPairsU;iU++)demuxedPairTimeU[iU]-=stmin; for(int iV=0;iV<NdemuxedPairsV;iV++)demuxedPairTimeV[iV]-=stmin; bool ok=false; if(nplanesU>2&&nplanesV>2){ ok = true; ok = this->Draw(1); if(ok){ cout << "tagged as interesting " << endl; } } if(amPSing)fCanvas->SaveAs("MyAna.eps"); fList->Delete(); return; }

bool AltDeMuxDisplay::Draw (  int  ipad  )

Definition at line 461 of file AltDeMuxDisplay.cxx. References _maxPlane, demuxedHitPlaneU, demuxedHitPlaneV, demuxedHitQU, demuxedHitQV, demuxedHitStripU, demuxedHitStripV, demuxedHitTU, demuxedHitTV, demuxedPairPlaneU, demuxedPairPlaneV, demuxedPairTimeU, demuxedPairTimeV, DISPLAY_TYPE, fCanvas, FindNearestEdge(), fList, fUpad, fUview, fVpad, fVview, fXYpad, fXYview, fZQpad, fZQview, fZTpad, fZTview, NdemuxedPairsU, NdemuxedPairsV, NVeto, pBoxSM1, pBoxSM2, pCircle, pCut1, pOctagon, pOctagonC, QPlane, QSPlane, xVeto, and yVeto. Referenced by Ana(). { cout << "AltDeMuxDisplay::Draw()" << endl; TPolyMarker* pTPU1; TPolyMarker* pTPV1; TPolyMarker* pTPU2; TPolyMarker* pTPV2; TPolyMarker* pTPU3; TPolyMarker* pTPV3; TPolyMarker* pTPU4; TPolyMarker* pTPV4; TPolyMarker* pTPU5; TPolyMarker* pTPV5; TPolyMarker* pTPXY1; TPolyMarker* pTPXY2; TPolyMarker* pTPXY3; TPolyMarker* pTPXY4; TPolyMarker* pTPXY5; TPolyMarker* pTPZTU; TPolyMarker* pTPZTV; TPolyMarker* pTPVeto1; int fidPlane[500]; float xU1[1000]; float yU1[1000]; float xU2[1000]; float yU2[1000]; float xU3[1000]; float yU3[1000]; float xU4[1000]; float yU4[1000]; float xU5[1000]; float yU5[1000]; int nU1=0; int nU2=0; int nU3=0; int nU4=0; int nU5=0; float xV1[1000]; float yV1[1000]; float xV2[1000]; float yV2[1000]; float xV3[1000]; float yV3[1000]; float xV4[1000]; float yV4[1000]; float xV5[1000]; float yV5[1000]; int nV1=0; int nV2=0; int nV3=0; int nV4=0; int nV5=0; float x1[1000]; float y1[1000]; float x2[1000]; float y2[1000]; float x3[1000]; float y3[1000]; float x4[1000]; float y4[1000]; float x5[1000]; float y5[1000]; int n1=0; int n2=0; int n3=0; int n4=0; int n5=0; float zU[5000]; float zV[5000]; float tU[5000]; float tV[5000]; int nU=0; int nV=0; float QFid=0.; float d; int iedge; float QFiducial[12]; float QFiducialPC[12]; for(int i=0;i<12;i++)QFiducial[i]=0.0; for(int i=0;i<500;i++)fidPlane[i]=0; for(int i=0;i<12;i++)QFiducialPC[i]=0.0; for(int i=0;i<NdemuxedHitsU;i++){ for(int j=i+1;j<NdemuxedHitsU;j++){ if(demuxedHitPlaneU[i]==demuxedHitPlaneU[j]){ if(demuxedHitStripU[i]==demuxedHitStripU[j]){ demuxedHitQU[i] +=demuxedHitQU[j]; demuxedHitQU[j] = 0; } } } } float planeMin =1000; float planeMax =0; for(int i=0;i<NdemuxedHitsU;i++){ int iU = static_cast<int>(demuxedHitPlaneU[i]); fZQview->Fill(demuxedHitPlaneU[i],demuxedHitQU[i]); zU[nU] = demuxedHitPlaneU[i]; tU[nU] = demuxedHitTU[i]; nU++; bool ok = false; float V; float U = demuxedHitStripU[i]; float x; float y; if(iU==0){ }else{ if(QPlane[iU-1]>0 && QPlane[iU+1]>0){ float sm = QSPlane[iU-1]/QPlane[iU-1]; float sp = QSPlane[iU+1]/QPlane[iU+1]; V = (sm+sp)/2.0; x = ((U-V)*0.0417)/1.4142; y = ((U+V)*0.0417-8.0)/1.4142; ok = true; } if(!ok){ if(QPlane[iU-3]>0 && QPlane[iU-1]>0){ float smm = QSPlane[iU-3]/QPlane[iU-3]; float sm = QSPlane[iU-1]/QPlane[iU-1]; V = (3*sm-smm)/2.0; x = ((U-V)*0.0417)/1.4142; y = ((U+V)*0.0417-8.0)/1.4142; ok = true; } } if(!ok){ if(QPlane[iU+3]>0 && QPlane[iU+1]>0){ float spp = QSPlane[iU+3]/QPlane[iU+3]; float sp = QSPlane[iU+1]/QPlane[iU+1]; V = (3*sp-spp)/2.0; x = ((U-V)*0.0417)/1.4142; y = ((U+V)*0.0417-8.0)/1.4142; ok = true; } } } if(ok){ this->FindNearestEdge(x,y,iedge,d); if(demuxedHitPlaneU[i]<=5)iedge=0; if(amCutting){ if(demuxedHitPlaneU[i]>=119 &&demuxedHitPlaneU[i]<=124 )iedge=9; if(demuxedHitPlaneU[i]>=119)iedge=9; }else{ if(demuxedHitPlaneU[i]>=243&&demuxedHitPlaneU[i]<=249)iedge=9; if(demuxedHitPlaneU[i]>=269&&demuxedHitPlaneU[i]<=275)iedge=10; if(demuxedHitPlaneU[i]>=_maxPlane+15)iedge=11; } if(d<0.5||iedge==0||iedge>=9)QFiducial[iedge]+=demuxedHitQU[i]; if(d>=0.5&&iedge!=0&&iedge<9)fidPlane[static_cast<int>(demuxedHitPlaneU[i])]++; if(demuxedHitPlaneU[i]>=planeMin && demuxedHitPlaneU[i]<=planeMax){ if(d<0.5||iedge==0||iedge>=9)QFiducialPC[iedge]+=demuxedHitQU[i]; } if(demuxedHitQU[i]>=0.1 && demuxedHitQU[i]<=2.0){ x1[n1] = x; y1[n1] = y; n1++; } if(demuxedHitQU[i]>=2 && demuxedHitQU[i]<=5){ x2[n2] = x; y2[n2] = y; n2++; } if(demuxedHitQU[i]>=5 && demuxedHitQU[i]<=25){ x3[n3] = x; y3[n3] = y; n3++; } if(demuxedHitQU[i]>=25 && demuxedHitQU[i]<=100){ x4[n4] = x; y4[n4] = y; n4++; } if(demuxedHitQU[i]>=100){ x5[n5] = x; y5[n5] = y; n5++; } } if(demuxedHitQU[i]>=1 && demuxedHitQU[i]<=3){ xU1[nU1] = demuxedHitPlaneU[i]; yU1[nU1] = demuxedHitStripU[i]; nU1++; } if(demuxedHitQU[i]>=3 && demuxedHitQU[i]<=10){ xU2[nU2] = demuxedHitPlaneU[i]; yU2[nU2] = demuxedHitStripU[i]; nU2++; } if(demuxedHitQU[i]>=10 && demuxedHitQU[i]<=50){ xU3[nU3] = demuxedHitPlaneU[i]; yU3[nU3] = demuxedHitStripU[i]; nU3++; } if(demuxedHitQU[i]>=50 && demuxedHitQU[i]<=200){ xU4[nU4] = demuxedHitPlaneU[i]; yU4[nU4] = demuxedHitStripU[i]; nU4++; } if(demuxedHitQU[i]>=200){ xU5[nU5] = demuxedHitPlaneU[i]; yU5[nU5] = demuxedHitStripU[i]; nU5++; } } for(int i=0;i<NdemuxedHitsV;i++){ for(int j=i+1;j<NdemuxedHitsV;j++){ if(demuxedHitPlaneV[i]==demuxedHitPlaneV[j]){ if(demuxedHitStripV[i]==demuxedHitStripV[j]){ demuxedHitQV[i] +=demuxedHitQV[j]; demuxedHitQV[j] = 0; } } } } for(int i=0;i<NdemuxedHitsV;i++){ int iV = static_cast<int>(demuxedHitPlaneV[i]); fZQview->Fill(demuxedHitPlaneV[i],demuxedHitQV[i]); zV[nU] = demuxedHitPlaneV[i]; tV[nU] = demuxedHitTV[i]; nV++; bool ok = false; float U; float V = demuxedHitStripV[i]; float x; float y; if(iV==0){ }else{ if(QPlane[iV-1]>0 && QPlane[iV+1]>0){ float sm = QSPlane[iV-1]/QPlane[iV-1]; float sp = QSPlane[iV+1]/QPlane[iV+1]; U = (sm+sp)/2.0; x = ((U-V)*0.0417)/1.4142; y = ((U+V)*0.0417-8.0)/1.4142; ok = true; } if(!ok){ if(QPlane[iV-3]>0 && QPlane[iV-1]>0){ float smm = QSPlane[iV-3]/QPlane[iV-3]; float sm = QSPlane[iV-1]/QPlane[iV-1]; U = (3*sm-smm)/2.0; x = ((U-V)*0.0417)/1.4142; y = ((U+V)*0.0417-8.0)/1.4142; ok = true; } } if(!ok){ if(QPlane[iV+3]>0 && QPlane[iV+1]>0){ float spp = QSPlane[iV+3]/QPlane[iV+3]; float sp = QSPlane[iV+1]/QPlane[iV+1]; U = (3*sp-spp)/2.0; x = ((U-V)*0.0417)/1.4142; y = ((U+V)*0.0417-8.0)/1.4142; ok = true; } } } if(ok){ this->FindNearestEdge(x,y,iedge,d); if(demuxedHitPlaneV[i]<=5)iedge=0; if(amCutting){ if(demuxedHitPlaneV[i]>=119 &&demuxedHitPlaneV[i]<=124 )iedge=9; if(demuxedHitPlaneV[i]>=119)iedge=9; }else{ if(demuxedHitPlaneV[i]>=243&&demuxedHitPlaneV[i]<=249)iedge=9; if(demuxedHitPlaneV[i]>=269&&demuxedHitPlaneV[i]<=275)iedge=10; if(demuxedHitPlaneV[i]>=_maxPlane+15)iedge=11; } if(d<0.5||iedge==0||iedge>=9)QFiducial[iedge]+=demuxedHitQV[i]; if(d>=0.5&&iedge!=0&&iedge<9)fidPlane[static_cast<int>(demuxedHitPlaneV[i])]++; if(demuxedHitPlaneV[i]>=planeMin && demuxedHitPlaneV[i]<=planeMax){ if(d<0.5||iedge==0||iedge>=9)QFiducialPC[iedge]+=demuxedHitQV[i]; } if(demuxedHitQV[i]>=0.1 && demuxedHitQV[i]<=2){ x1[n1] = x; y1[n1] = y; n1++; } if(demuxedHitQV[i]>=2 && demuxedHitQV[i]<=5){ x2[n2] = x; y2[n2] = y; n2++; } if(demuxedHitQV[i]>=5 && demuxedHitQV[i]<=25){ x3[n3] = x; y3[n3] = y; n3++; } if(demuxedHitQV[i]>=25 && demuxedHitQV[i]<=100){ x4[n4] = x; y4[n4] = y; n4++; } if(demuxedHitQV[i]>=100){ x5[n5] = x; y5[n5] = y; n5++; } } if(demuxedHitQV[i]>=0.1 && demuxedHitQV[i]<=2){ xV1[nV1] = demuxedHitPlaneV[i]; yV1[nV1] = demuxedHitStripV[i]; nV1++; } if(demuxedHitQV[i]>=2 && demuxedHitQV[i]<=5){ xV2[nV2] = demuxedHitPlaneV[i]; yV2[nV2] = demuxedHitStripV[i]; nV2++; } if(demuxedHitQV[i]>=5 && demuxedHitQV[i]<=25){ xV3[nV3] = demuxedHitPlaneV[i]; yV3[nV3] = demuxedHitStripV[i]; nV3++; } if(demuxedHitQV[i]>=25 && demuxedHitQV[i]<=100){ xV4[nV4] = demuxedHitPlaneV[i]; yV4[nV4] = demuxedHitStripV[i]; nV4++; } if(demuxedHitQV[i]>=100){ xV5[nV5] = demuxedHitPlaneV[i]; yV5[nV5] = demuxedHitStripV[i]; nV5++; } } pTPVeto1 = new TPolyMarker(NVeto,xVeto,yVeto); pTPVeto1->SetMarkerStyle(20); pTPVeto1->SetMarkerColor(kRed); pTPVeto1->SetMarkerSize(1.0); pTPZTU = new TPolyMarker(NdemuxedPairsU,demuxedPairPlaneU,demuxedPairTimeU); pTPZTU->SetMarkerStyle(20); pTPZTU->SetMarkerColor(kBlue); pTPZTU->SetMarkerSize(0.5); pTPZTV = new TPolyMarker(NdemuxedPairsV,demuxedPairPlaneV,demuxedPairTimeV); pTPZTV->SetMarkerStyle(20); pTPZTV->SetMarkerColor(kRed); pTPZTV->SetMarkerSize(0.5); pTPU1 = new TPolyMarker(nU1,xU1,yU1); pTPU1->SetMarkerStyle(20); pTPU1->SetMarkerColor(kYellow); pTPU1->SetMarkerSize(0.5); pTPU2 = new TPolyMarker(nU2,xU2,yU2); pTPU2->SetMarkerStyle(20); pTPU2->SetMarkerColor(kCyan); pTPU2->SetMarkerSize(0.5); pTPU3 = new TPolyMarker(nU3,xU3,yU3); pTPU3->SetMarkerStyle(20); pTPU3->SetMarkerColor(kBlue); pTPU3->SetMarkerSize(0.5); pTPU4 = new TPolyMarker(nU4,xU4,yU4); pTPU4->SetMarkerStyle(20); pTPU4->SetMarkerColor(kBlack); pTPU4->SetMarkerSize(0.5); pTPU5 = new TPolyMarker(nU5,xU5,yU5); pTPU5->SetMarkerStyle(20); pTPU5->SetMarkerColor(kRed); pTPU5->SetMarkerSize(0.5); pTPV1 = new TPolyMarker(nV1,xV1,yV1); pTPV1->SetMarkerStyle(20); pTPV1->SetMarkerColor(kYellow); pTPV1->SetMarkerSize(0.5); pTPV2 = new TPolyMarker(nV2,xV2,yV2); pTPV2->SetMarkerStyle(20); pTPV2->SetMarkerColor(kCyan); pTPV2->SetMarkerSize(0.5); pTPV3 = new TPolyMarker(nV3,xV3,yV3); pTPV3->SetMarkerStyle(20); pTPV3->SetMarkerColor(kBlue); pTPV3->SetMarkerSize(0.5); pTPV4 = new TPolyMarker(nV4,xV4,yV4); pTPV4->SetMarkerStyle(20); pTPV4->SetMarkerColor(kBlack); pTPV4->SetMarkerSize(0.5); pTPV5 = new TPolyMarker(nV5,xV5,yV5); pTPV5->SetMarkerStyle(20); pTPV5->SetMarkerColor(kRed); pTPV5->SetMarkerSize(0.5); pTPXY1 = new TPolyMarker(n1,x1,y1); pTPXY1->SetMarkerStyle(20); pTPXY1->SetMarkerColor(kYellow); pTPXY1->SetMarkerSize(0.5); pTPXY2 = new TPolyMarker(n2,x2,y2); pTPXY2->SetMarkerStyle(20); pTPXY2->SetMarkerColor(kCyan); pTPXY2->SetMarkerSize(0.5); pTPXY3 = new TPolyMarker(n3,x3,y3); pTPXY3->SetMarkerStyle(20); pTPXY3->SetMarkerColor(kBlue); pTPXY3->SetMarkerSize(0.5); pTPXY4 = new TPolyMarker(n4,x4,y4); pTPXY4->SetMarkerStyle(20); pTPXY4->SetMarkerColor(kBlack); pTPXY4->SetMarkerSize(0.5); pTPXY5 = new TPolyMarker(n5,x5,y5); pTPXY5->SetMarkerStyle(20); pTPXY5->SetMarkerColor(kRed); pTPXY5->SetMarkerSize(0.5); // Keep track of the TPolyMarkers - delete them later fList->Add(pTPV1); fList->Add(pTPV2); fList->Add(pTPV3); fList->Add(pTPV4); fList->Add(pTPV5); fList->Add(pTPU1); fList->Add(pTPU2); fList->Add(pTPU3); fList->Add(pTPU4); fList->Add(pTPU5); fList->Add(pTPVeto1); fList->Add(pTPXY1); fList->Add(pTPXY2); fList->Add(pTPXY3); fList->Add(pTPXY4); fList->Add(pTPXY5); fList->Add(pTPZTU); fList->Add(pTPZTV); if(DISPLAY_TYPE==1)fUpad->cd(); if (fUview) fUview->Draw(); pTPU1->Draw(); pTPU2->Draw(); pTPU3->Draw(); pTPU4->Draw(); pTPU5->Draw(); pBoxSM1->Draw(); pBoxSM2->Draw(); if(amCutting)pCut1->Draw(); if(DISPLAY_TYPE==1)fVpad->cd(); if (fVview) fVview->Draw(); pTPV1->Draw(); pTPV2->Draw(); pTPV3->Draw(); pTPV4->Draw(); pTPV5->Draw(); pBoxSM1->Draw(); pBoxSM2->Draw(); if(amCutting)pCut1->Draw(); if(DISPLAY_TYPE==1)fXYpad->cd(); if (fXYview) fXYview->Draw(); pTPXY1->Draw(); pTPXY2->Draw(); pTPXY3->Draw(); pTPXY4->Draw(); pTPXY5->Draw(); pTPVeto1->Draw(); pOctagon->Draw(); pOctagonC->Draw(); pCircle->Draw(); fZTpad->cd(); if(fZTview)fZTview->Draw(); pTPZTU->Draw(); pTPZTV->Draw(); fZQpad->cd(); if(fZQview)fZQview->Draw(); fCanvas->Update(); cout << " PEs outside FIDUCIAL cut : " << QFid << endl; int icount=0; int jcount=0; for(int i=0;i<12;i++){ cout << " Edge " << i << " : " << QFiducial[i] << " : " << QFiducialPC[i] << endl; if(QFiducial[i]>1.0)icount++; if(QFiducialPC[i]>1.0)jcount++; } int nInFid=0; for(int i=0;i<500;i++){ if(fidPlane[i]>0)nInFid++; } cout << "Hit Planes in Fiducial region : " << nInFid << endl; if(nInFid<6)return false; if(icount==0)cout << " CONTAINED " << endl; if(icount==1&&jcount==0)cout << " PARTIALLY CONTAINED " << endl; if(icount==1&&jcount==1)cout << " STOPPER " << endl; if(icount>=2&&jcount==1)cout << " THROUGH-GOING " << endl; if(icount>0&&jcount>0)return false; getchar(); return true; }

void AltDeMuxDisplay::FindNearestEdge (  float  x, float  y, int &  idege, float &  d )

Definition at line 402 of file AltDeMuxDisplay.cxx. Referenced by Draw(). { float d1 = fabs(x); float d2 = fabs(y); float d3 = 0.70711*fabs(x-y); float d4 = 0.70711*fabs(x+y); if(d1>d2){ if(d3>d4){ if(d1>d3){ d = 4.0-d1; if(x>0)iedge=3; if(x<=0)iedge=7; }else{ d = 4.0-d3; if(x-y>0)iedge=4; if(x-y<=0)iedge=8; } }else{ if(d1>d4){ d = 4.0-d1; if(x>0)iedge=3; if(x<=0)iedge=7; }else{ d = 4.0-d4; if(x+y>0)iedge=2; if(x+y<=0)iedge=6; } } }else{ if(d3>d4){ if(d2>d3){ d = 4.0-d2; if(y>0)iedge=1; if(y<=0)iedge=5; }else{ d = 4.0-d3; if(x-y>0)iedge=4; if(x-y<=0)iedge=8; } }else{ if(d2>d4){ d = 4.0-d2; if(y>0)iedge=1; if(y<=0)iedge=5; }else{ d = 4.0-d4; if(x+y>0)iedge=2; if(x+y<=0)iedge=6; } } } return; }

int AltDeMuxDisplay::MaxPlane (  int  irun  )

Definition at line 1041 of file AltDeMuxDisplay.cxx. References _lastRun, _maxPlane, pBoxSM2, and pCut1. Referenced by Ana(). { if(irun==_lastRun)return _maxPlane; if(irun<= 1000)_maxPlane = 497; if(irun>= 8000 && irun < 8504 )_maxPlane = 293; // 8039 8493 : 293 if(irun>= 8504 && irun < 8803 )_maxPlane = 305; // 8580 8802 : 305 if(irun>= 8803 && irun < 9279 )_maxPlane = 313; if(irun>= 9279 && irun < 9279 )_maxPlane = 317; if(irun>= 9279 && irun < 9530 )_maxPlane = 321; if(irun>= 9530 && irun < 9700 )_maxPlane = 325; if(irun>= 9700 && irun < 9988 )_maxPlane = 333; if(irun>= 9988 && irun <10002 )_maxPlane = 341; if(irun>=10002 && irun <10299 )_maxPlane = 345; if(irun>=10299 && irun <10700 )_maxPlane = 349; if(irun>=10700 && irun <11400 )_maxPlane = 365; if(irun>=11400 && irun <12433 )_maxPlane = 377; if(irun>=12433 && irun <99999 )_maxPlane = 498; _lastRun = irun; float x[5]; float y[5]; x[0] = 269.; x[1] = _maxPlane+20.; x[2] = _maxPlane+20.; x[3] = 269.; x[4] = 269.; y[0] = 0.; y[1] = 0.; y[2] = 192.; y[3] = 192.; y[4] = 0.; if(pBoxSM2)delete pBoxSM2; pBoxSM2 = new TPolyLine(5,x,y); x[0] = 124.; y[0] = 0.; x[1] = 124.; y[1] = 192.; pCut1 = new TPolyLine(2,x,y); return _maxPlane; }

---

Member Data Documentation

std::vector<DeMuxedPair> AltDeMuxDisplay::_demuxedPairs [private]

Definition at line 79 of file AltDeMuxDisplay.h. Referenced by Ana().

int AltDeMuxDisplay::_lastRun [private]

Definition at line 45 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and MaxPlane().

int AltDeMuxDisplay::_maxPlane [private]

Definition at line 44 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), Draw(), and MaxPlane().

bool AltDeMuxDisplay::amCutting [private]

Definition at line 41 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay().

bool AltDeMuxDisplay::amPSing [private]

Definition at line 40 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay().

bool AltDeMuxDisplay::amVetoing [private]

Definition at line 42 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Ana().

float AltDeMuxDisplay::demuxedHitPlaneU[20000] [private]

Definition at line 50 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedHitPlaneV[20000] [private]

Definition at line 57 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedHitQU[20000] [private]

Definition at line 52 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedHitQV[20000] [private]

Definition at line 59 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedHitStripU[20000] [private]

Definition at line 51 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedHitStripV[20000] [private]

Definition at line 58 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedHitTU[20000] [private]

Definition at line 53 of file AltDeMuxDisplay.h. Referenced by Draw().

float AltDeMuxDisplay::demuxedHitTV[20000] [private]

Definition at line 60 of file AltDeMuxDisplay.h. Referenced by Draw().

float AltDeMuxDisplay::demuxedPairPlaneU[20000] [private]

Definition at line 65 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedPairPlaneV[20000] [private]

Definition at line 64 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedPairTimeU[20000] [private]

Definition at line 67 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::demuxedPairTimeV[20000] [private]

Definition at line 66 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

TCanvas* AltDeMuxDisplay::fCanvas [private]

Definition at line 81 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), Ana(), and Draw().

TList* AltDeMuxDisplay::fList [private]

Definition at line 93 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), Ana(), and Draw().

TPad* AltDeMuxDisplay::fUpad [private]

Definition at line 82 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), Ana(), and Draw().

TH2F* AltDeMuxDisplay::fUview [private]

Definition at line 87 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

TPad* AltDeMuxDisplay::fVpad [private]

Definition at line 83 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), Ana(), and Draw().

TH2F* AltDeMuxDisplay::fVview [private]

Definition at line 88 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

TPad* AltDeMuxDisplay::fXYpad [private]

Definition at line 84 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Draw().

TH2F* AltDeMuxDisplay::fXYview [private]

Definition at line 89 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

TPad* AltDeMuxDisplay::fZQpad [private]

Definition at line 85 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Draw().

TH1F* AltDeMuxDisplay::fZQview [private]

Definition at line 90 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

TPad* AltDeMuxDisplay::fZTpad [private]

Definition at line 86 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Draw().

TH2F* AltDeMuxDisplay::fZTview [private]

Definition at line 91 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

int AltDeMuxDisplay::NdemuxedHitsU [private]

Definition at line 48 of file AltDeMuxDisplay.h. Referenced by Ana().

int AltDeMuxDisplay::NdemuxedHitsV [private]

Definition at line 56 of file AltDeMuxDisplay.h. Referenced by Ana().

int AltDeMuxDisplay::NdemuxedPairsU [private]

Definition at line 63 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

int AltDeMuxDisplay::NdemuxedPairsV [private]

Definition at line 62 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

int AltDeMuxDisplay::NVeto [private]

Definition at line 69 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

TPolyLine* AltDeMuxDisplay::pBoxSM1 [private]

Definition at line 99 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Draw().

TPolyLine* AltDeMuxDisplay::pBoxSM2 [private]

Definition at line 100 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), Draw(), and MaxPlane().

TArc* AltDeMuxDisplay::pCircle [private]

Definition at line 103 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Draw().

TPolyLine* AltDeMuxDisplay::pCut1 [private]

Definition at line 95 of file AltDeMuxDisplay.h. Referenced by Draw(), and MaxPlane().

TPolyLine* AltDeMuxDisplay::pCut2 [private]

Definition at line 96 of file AltDeMuxDisplay.h.

std::vector<PlexSEIdAltL*> AltDeMuxDisplay::PlanesAltListsE[500] [private]

Definition at line 77 of file AltDeMuxDisplay.h. Referenced by Ana().

std::vector<PlexSEIdAltL*> AltDeMuxDisplay::PlanesAltListsW[500] [private]

Definition at line 78 of file AltDeMuxDisplay.h. Referenced by Ana().

TPolyLine* AltDeMuxDisplay::pOctagon [private]

Definition at line 101 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Draw().

TPolyLine* AltDeMuxDisplay::pOctagonC [private]

Definition at line 102 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Draw().

TPolyLine* AltDeMuxDisplay::pSM1 [private]

Definition at line 97 of file AltDeMuxDisplay.h.

TPolyLine* AltDeMuxDisplay::pSM2 [private]

Definition at line 98 of file AltDeMuxDisplay.h.

TPolyLine* AltDeMuxDisplay::pZTFit [private]

Definition at line 94 of file AltDeMuxDisplay.h.

float AltDeMuxDisplay::QPlane[520] [private]

Definition at line 55 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::QSPlane[520] [private]

Definition at line 54 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

int AltDeMuxDisplay::xMaxPlane [private]

Definition at line 46 of file AltDeMuxDisplay.h. Referenced by AltDeMuxDisplay(), and Ana().

float AltDeMuxDisplay::xVeto[1000] [private]

Definition at line 70 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::yVeto[1000] [private]

Definition at line 71 of file AltDeMuxDisplay.h. Referenced by Ana(), and Draw().

float AltDeMuxDisplay::zVeto[1000] [private]

Definition at line 72 of file AltDeMuxDisplay.h. Referenced by Ana().

---

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

• AltDeMuxDisplay.h
• AltDeMuxDisplay.cxx

---