AlgDeMuxGolden Class Reference

#include <AlgDeMuxGolden.h>

Inheritance diagram for AlgDeMuxGolden:

AlgBase List of all members.

Public Member Functions

 AlgDeMuxGolden ()
virtual ~AlgDeMuxGolden ()
virtual void Trace (const char *c) const
virtual void RunAlg (AlgConfig &acd, CandHandle &ch, CandContext &cx)

Private Member Functions

void SetPlanesToDeterminedFit (DmxPlaneItr &planeItr, Double_t a, Double_t b, Int_t firstPlane, Int_t lastPlane)
void UseGoldenFit (DmxPlaneItr &leadPlaneItr, DmxPlaneItr &goldenPlaneItr, Int_t vertexPlane, Int_t endPlane)

Private Attributes

Int_t fStrayPlanes
UInt_t fPlanesInSet
Int_t fStrayCut

Detailed Description

Definition at line 21 of file AlgDeMuxGolden.h.

Constructor & Destructor Documentation

AlgDeMuxGolden::AlgDeMuxGolden (  ) 

Definition at line 78 of file AlgDeMuxGolden.cxx.

00078                                :
00079   fStrayPlanes(0),
00080   fPlanesInSet(6),
00081   fStrayCut(6)
00082 {
00083   //default constructor
00084 }

AlgDeMuxGolden::~AlgDeMuxGolden (  )  [virtual]

Definition at line 88 of file AlgDeMuxGolden.cxx.

00089 {
00090 
00091 }

Member Function Documentation

void AlgDeMuxGolden::RunAlg ( AlgConfig &  acd,
CandHandle &  ch,
CandContext &  cx 
) [virtual]

Implements AlgBase.

Definition at line 95 of file AlgDeMuxGolden.cxx.

References DmxUtilities::CheckFit(), DmxUtilities::FillPlaneArray(), DmxUtilities::FindEndPlane(), DmxUtilities::FindVertexPlane(), fPlanesInSet, fStrayPlanes, Registry::GetDouble(), DmxStatus::GetEndPlaneNumber(), DmxStatus::GetEventDeMuxed(), DmxStatus::GetEventNumber(), DmxStatus::GetFigureOfMeritFailure(), Registry::GetInt(), DmxStatus::GetPlaneArray(), DmxStatus::GetVertexPlaneNumber(), DmxStatus::GetVertexZPosition(), DmxUtilities::IsOverlappingMultiple(), Msg::kDebug, KeyOnGoldenU(), KeyOnGoldenV(), KeyOnLeadU(), KeyOnLeadV(), KeyPlane(), DmxPlaneTypes::kShower, Msg::kWarning, MSG, DmxStatus::ResetStatus(), DmxStatus::SetEndPlaneNumber(), DmxStatus::SetFigureOfMerit(), DmxStatus::SetMultipleMuon(), DmxStatus::SetNumberOfPlanes(), DmxStatus::SetUOverlappingMultiple(), DmxStatus::SetUStrayPlanes(), DmxStatus::SetUValidPlanes(), DmxStatus::SetValidPlanesFailure(), DmxStatus::SetVertexPlaneNumber(), DmxStatus::SetVertexPlaneZPosition(), DmxStatus::SetVOverlappingMultiple(), DmxStatus::SetVStrayPlanes(), DmxStatus::SetVValidPlanes(), and UseGoldenFit().

00096 {
00097 
00098   assert( ch.InheritsFrom("CandDeMuxDigitListHandle") );
00099   CandDeMuxDigitListHandle &cdlh = dynamic_cast<CandDeMuxDigitListHandle&>(ch);
00100   //get the AlgConfigDeMux object
00101   fPlanesInSet = acd.GetInt("PlanesInSet");
00102   //MSG("DMX", Msg::kInfo) << "window size = " << fPlanesInSet << "\tcut = " << fStrayCut << endl;
00103 
00104   //get the DmxStatus object
00105   // Find the DmxStatus object or create one for needed scratch space
00106   DmxStatus *status = dynamic_cast<DmxStatus *>(gROOT->GetRootFolder()->FindObject("Loon/DeMux/DmxStatus"));
00107   bool tempStatus = false;
00108   if (status == 0) {
00109     MSG("DMXX", Msg::kDebug) << "//root/Loon/DeMux/DmxStatus not found."
00110                               << " Create a temporary DmxStatus." << endl;
00111     status = new DmxStatus;      // Make a temporary DmxStatus if needed
00112     tempStatus = true;
00113   }
00114   else status->ResetStatus(); 
00115       
00116   //instantiate the DmxUtilities object and then fill DmxStatus with the event information      
00117   DmxUtilities util;
00118   
00119   util.FillPlaneArray(status, cdlh, acd);
00120   const TObjArray *planeArray = status->GetPlaneArray();
00121   //DmxPlane *plane = dynamic_cast<DmxPlane *>(planeArray->First());
00122   //MSG("DMX", Msg::kInfo) <<"Event = " <<  status->GetEventNumber() << endl;
00123   
00124   //create the DmxPlaneItr over planes and program it to sort the planes
00125   DmxPlaneItr planeItr(planeArray);
00126   DmxPlaneKeyFunc *pnKF = planeItr.CreateKeyFunc();
00127   pnKF->SetFun(KeyPlane);
00128   planeItr.GetSet()->AdoptSortKeyFunc(pnKF);
00129   pnKF = 0;
00130   
00131   planeItr.ResetFirst();
00132   status->SetNumberOfPlanes(planeItr.SizeSelect());
00133   status->SetVertexPlaneNumber(util.FindVertexPlane(planeItr));
00134   
00135   //demux the event if there is an identified vertex
00136   if( status->GetVertexPlaneNumber() != -1){
00137     
00138     planeItr.ResetFirst();
00139     planeItr.GetSet()->Slice(status->GetVertexPlaneNumber(), 500);
00140     //set the end plane number since there was a vertex plane
00141     status->SetEndPlaneNumber(util.FindEndPlane(planeItr));
00142     planeItr.GetSet()->Slice();
00143     planeItr.ResetFirst();
00144 
00145     //set the z position of the vertex plane
00146     planeItr.GetSet()->Slice(status->GetVertexPlaneNumber());
00147     planeItr.ResetFirst();
00148     status->SetVertexPlaneZPosition(dynamic_cast<DmxPlane *>(planeItr.Ptr())->GetZPosition());
00149 
00150     planeItr.GetSet()->Slice();
00151     planeItr.ResetFirst();
00152 
00153     //DmxPlane *plane = dynamic_cast<DmxPlane *>(planeArray->First());
00154     //MSG("DMX", Msg::kInfo) <<"Event = " <<  status->GetEventNumber() << endl;
00155     
00156     //create two DmxPlaneItrs
00157     DmxPlaneItr leadPlaneItr(planeArray);
00158     DmxPlaneItr vertexPlaneItr(planeArray);
00159     DmxPlaneItr goldenPlaneItr(planeArray);
00160     
00161     //create a KeyFunc to sort planes by number
00162     DmxPlaneKeyFunc *pnKF = leadPlaneItr.CreateKeyFunc();
00163     DmxPlaneKeyFunc *vertexPlaneKF = vertexPlaneItr.CreateKeyFunc();
00164     DmxPlaneKeyFunc *goldenPlaneKF = goldenPlaneItr.CreateKeyFunc();
00165     
00166     //program the KeyFunc with the sort function
00167     pnKF->SetFun(KeyPlane);
00168     vertexPlaneKF->SetFun(KeyPlane);
00169     goldenPlaneKF->SetFun(KeyPlane);
00170     
00171     //get the NavSet from the iterator and pass the KeyFunc to it
00172     leadPlaneItr.GetSet()->AdoptSortKeyFunc(pnKF);
00173     vertexPlaneItr.GetSet()->AdoptSortKeyFunc(vertexPlaneKF);
00174     goldenPlaneItr.GetSet()->AdoptSortKeyFunc(goldenPlaneKF);
00175     
00176     //clear the KF pointer because we no longer own the KeyFunc
00177     
00178     pnKF = 0;
00179     vertexPlaneKF = 0;
00180     goldenPlaneKF = 0;
00181     leadPlaneItr.ResetFirst();
00182     goldenPlaneItr.ResetFirst();
00183     vertexPlaneItr.ResetFirst();
00184 
00185     //now program a key function to select on plane orientation - U first
00186     DmxPlaneKeyFunc *orientGoldenUKF = goldenPlaneItr.CreateKeyFunc();
00187     DmxPlaneKeyFunc *orientLeadUKF = leadPlaneItr.CreateKeyFunc();
00188     DmxPlaneKeyFunc *orientGoldenVKF = goldenPlaneItr.CreateKeyFunc();
00189     DmxPlaneKeyFunc *orientLeadVKF = leadPlaneItr.CreateKeyFunc();
00190     
00191     Int_t viewCtr = 0;
00192     while( viewCtr < 2 ){
00193      
00194       if( viewCtr == 0){
00195         //program this function with the select function
00196         orientLeadUKF->SetFun(KeyOnLeadU);
00197         orientGoldenUKF->SetFun(KeyOnGoldenU);
00198                 
00199         //adopt it as a selection function
00200         leadPlaneItr.GetSet()->AdoptSelectKeyFunc(orientLeadUKF);
00201         orientLeadUKF = 0;
00202         leadPlaneItr.Reset();
00203         goldenPlaneItr.GetSet()->AdoptSelectKeyFunc(orientGoldenUKF);
00204         orientGoldenUKF = 0;
00205         goldenPlaneItr.Reset();
00206       }
00207       else if(viewCtr == 1){
00208         //program this function with the select function
00209         orientLeadVKF->SetFun(KeyOnLeadV);
00210         orientGoldenVKF->SetFun(KeyOnGoldenV);
00211                 
00212         //adopt it as a selection function
00213         leadPlaneItr.GetSet()->AdoptSelectKeyFunc(orientLeadVKF);
00214         orientLeadVKF = 0;
00215         leadPlaneItr.Reset();
00216         goldenPlaneItr.GetSet()->AdoptSelectKeyFunc(orientGoldenVKF);
00217         orientGoldenVKF = 0;
00218         goldenPlaneItr.Reset();
00219       }
00220 
00221       //slice to those planes from the vertex on
00222       goldenPlaneItr.GetSet()->Slice(status->GetVertexPlaneNumber(), status->GetEndPlaneNumber());
00223 
00224       if(goldenPlaneItr.SizeSelect() >=2 ){
00225 
00226         //MSG("DMX", Msg::kInfo) << "event " << status->GetEventNumber() << endl;
00227         //loop over the golden planes and set them
00228         while(goldenPlaneItr.IsValid() ){
00229 
00230           //MSG("DMX", Msg::kInfo) << dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetPlaneNumber() << endl;
00231           if(dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetPlaneType() == DmxPlaneTypes::kShower){
00232             dynamic_cast<DmxShowerPlane *>(goldenPlaneItr.Ptr())->SetStrips("best");
00233           }
00234           goldenPlaneItr.Next();
00235         }//end setting golden planes
00236 
00237         goldenPlaneItr.Reset();
00238 
00239         UseGoldenFit(leadPlaneItr, goldenPlaneItr, status->GetVertexPlaneNumber(), 
00240                      status->GetEndPlaneNumber());
00241 
00242         //check the fit
00243         fStrayPlanes -= util.CheckFit(leadPlaneItr);
00244 
00245         //what about multiple muons?  look for another vertex after the end plane
00246         leadPlaneItr.GetSet()->Slice();
00247         goldenPlaneItr.GetSet()->Slice();
00248         vertexPlaneItr.GetSet()->Slice(status->GetEndPlaneNumber()+1, 500);
00249         //MSG("DMX1", Msg::kInfo)<< "event =  " << status->GetEventNumber() 
00250         //                   << "\tend plane = " << status->GetEndPlaneNumber() << endl;
00251         
00252         Int_t nextVertex = -1;
00253         if(vertexPlaneItr.SizeSelect() > 3) nextVertex = util.FindVertexPlane(vertexPlaneItr);
00254         
00255         while( nextVertex != -1){
00256           
00257           status->SetMultipleMuon(true); 
00258           vertexPlaneItr.GetSet()->Slice();
00259           vertexPlaneItr.GetSet()->Slice(nextVertex, 500);
00260           
00261           //get the next end plane
00262           Int_t endPlane = util.FindEndPlane(vertexPlaneItr);
00263 
00264           //MSG("DMX1", Msg::kInfo)<< "\tin double muon loop, vertex at plane " << nextVertex 
00265           //             << "\tend plane = " << endPlane << "\tplanes in slice = " 
00266           //             << vertexPlaneItr.SizeSelect() << endl;
00267           
00268           goldenPlaneItr.GetSet()->Slice(nextVertex, endPlane);
00269           
00270           if(goldenPlaneItr.SizeSelect() >= 2 ){
00271             //loop over the golden planes and set them
00272             while(goldenPlaneItr.IsValid() ){
00273               
00274               if(dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetPlaneType() == DmxPlaneTypes::kShower){
00275                 dynamic_cast<DmxShowerPlane *>(goldenPlaneItr.Ptr())->SetStrips("best");
00276               }
00277               goldenPlaneItr.Next();
00278             }
00279             
00280             goldenPlaneItr.Reset();
00281             
00282             UseGoldenFit(leadPlaneItr, goldenPlaneItr, status->GetVertexPlaneNumber(), 
00283                          status->GetEndPlaneNumber());
00284             
00285             //check the fit
00286             fStrayPlanes -= util.CheckFit(leadPlaneItr);
00287             
00288             //CheckFitWithTiming(ch, planeItr);
00289             
00290           }//end if enough golden planes to try
00291           
00292           //look for another vertex after the end plane
00293           leadPlaneItr.GetSet()->Slice();
00294           goldenPlaneItr.GetSet()->Slice();
00295           vertexPlaneItr.GetSet()->Slice();
00296           vertexPlaneItr.GetSet()->Slice(endPlane+1, 500);
00297           //MSG("DMX1", Msg::kInfo)<< "\tlook for new vertex, slice from " << endPlane+1 
00298           //             << "\tto plane 500, planes in slice = " 
00299           //             << vertexPlaneItr.SizeSelect() << endl;
00300           
00301           vertexPlaneItr.Reset();
00302           
00303           if(vertexPlaneItr.SizeSelect() > 0){
00304             nextVertex = util.FindVertexPlane(vertexPlaneItr);
00305           }
00306           else{ nextVertex = -1; }
00307         }//end loop over new verticies
00308       }//end if at least 2 golden planes
00309       else{ 
00310         MSG("DMX", Msg::kWarning)<< "Event " << status->GetEventNumber() 
00311                                  << " not demuxed; not enough golden planes in view" << endl;
00312         status->SetValidPlanesFailure(true); 
00313       } //not enough planes to demux
00314 
00315       status->SetFigureOfMerit(goldenPlaneItr.SizeSelect(), fStrayPlanes);
00316       if(viewCtr == 0 && status->GetEventDeMuxed()){ 
00317         status->SetUStrayPlanes(fStrayPlanes);
00318         status->SetUValidPlanes(goldenPlaneItr.SizeSelect());
00319         if( status->GetFigureOfMeritFailure() )
00320         //see if the event failed the figure of Merit test, if so, see if it is an 
00321         //overlapping muon
00322         status->SetUOverlappingMultiple(util.IsOverlappingMultiple(goldenPlaneItr, status->GetVertexZPosition(), 
00323                                                               acd.GetDouble("HoughInterceptRMSLimit"), 
00324                                                               acd.GetDouble("HoughSlopeRMSLimit"),
00325                                                               acd.GetDouble("HoughPeakLimit")));
00326       }
00327       if(viewCtr == 1 && status->GetEventDeMuxed()){ 
00328         status->SetVStrayPlanes(fStrayPlanes);
00329         status->SetVValidPlanes(goldenPlaneItr.SizeSelect());
00330         //see if the event failed the figure of Merit test, if so, see if it is an 
00331         //overlapping muon
00332         if( status->GetFigureOfMeritFailure() )
00333           status->SetVOverlappingMultiple(util.IsOverlappingMultiple(goldenPlaneItr, status->GetVertexZPosition(),
00334                                                                 acd.GetDouble("HoughInterceptRMSLimit"), 
00335                                                                 acd.GetDouble("HoughSlopeRMSLimit"),
00336                                                                 acd.GetDouble("HoughPeakLimit")));
00337       }
00338 
00339       //done using the first view, so clear the selection function
00340       //MSG("Dmx", Msg::kInfo) <<"Done with View " << viewCtr << endl;
00341       vertexPlaneItr.GetSet()->Slice();
00342       leadPlaneItr.GetSet()->Slice();
00343       leadPlaneItr.GetSet()->AdoptSelectKeyFunc(0);
00344       leadPlaneItr.Reset();
00345       goldenPlaneItr.GetSet()->Slice();
00346       goldenPlaneItr.GetSet()->AdoptSelectKeyFunc(0);
00347       goldenPlaneItr.Reset();
00348 
00349       ++viewCtr;
00350     } //end loop over views
00351     
00352   }//end if event has a vertex
00353   return;
00354 }

void AlgDeMuxGolden::SetPlanesToDeterminedFit ( DmxPlaneItr &  planeItr,
Double_t  a,
Double_t  b,
Int_t  firstPlane,
Int_t  lastPlane 
) [private]

Definition at line 438 of file AlgDeMuxGolden.cxx.

Referenced by UseGoldenFit().

00441 {
00442   planeItr.ResetFirst();
00443  
00444   while( planeItr.IsValid() ){
00445 
00446     DmxPlane *plane = dynamic_cast<DmxPlane *>(planeItr.Ptr());
00447     //MSG("DMX", Msg::kInfo)<<"\tin SetPlanesToDeterminedFit " << plane->GetPlaneNumber() << endl;
00448       
00449     if(plane->GetPlaneNumber()>=firstPlane && plane->GetPlaneNumber()<=lastPlane){
00450       Double_t fitCoG = a + (plane->GetZPosition() * b);
00451       
00452       plane->SetStrips(fitCoG); 
00453     }//end if plane is in the window bounds
00454   
00455     planeItr.Next();
00456   }
00457 
00458   planeItr.Reset();
00459   
00460   return;
00461 }

void AlgDeMuxGolden::Trace ( const char *  c  )  const [virtual]

Reimplemented from AlgBase.

Definition at line 357 of file AlgDeMuxGolden.cxx.

00358 {
00359 }

void AlgDeMuxGolden::UseGoldenFit ( DmxPlaneItr &  leadPlaneItr,
DmxPlaneItr &  goldenPlaneItr,
Int_t  vertexPlane,
Int_t  endPlane 
) [private]

Definition at line 363 of file AlgDeMuxGolden.cxx.

References DmxUtilities::FindPlanesOffFit(), fStrayCut, fStrayPlanes, and SetPlanesToDeterminedFit().

Referenced by RunAlg().

00365 {
00366   
00367   //get the number of loops you have to iterate over.  you have (# golden planes - 1)
00368   //loops to go through them each
00369   Int_t numLoops = goldenPlaneItr.SizeSelect() - 1;
00370   Int_t loopCtr = 0;
00371   DmxUtilities util = DmxUtilities();
00372 
00373   goldenPlaneItr.Reset();
00374   
00375   while(loopCtr < numLoops){
00376     
00377     
00378     Int_t goldCtr = 0;
00379     Float_t goldCoG1 = -1.;
00380     Float_t goldCoG2 = -1.;
00381     Float_t goldZ1 = -1.;
00382     Float_t goldZ2 = -1.;
00383     Int_t firstPlane = -1;
00384     Int_t lastPlane = -1;
00385     //get the cog's for the next two golden planes
00386     while(goldCtr < 2 && goldenPlaneItr.IsValid()){
00387 
00388       //advance the iterator only if on the first plane of the set - that way the first plane
00389       //in the next set is the last in this one
00390       if(goldCtr == 0){
00391         firstPlane = dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetPlaneNumber();
00392         goldCoG1 = dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetCoG();
00393         goldZ1 = dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetZPosition();
00394         goldenPlaneItr.Next();
00395       }
00396       else if(goldCtr == 1){
00397         lastPlane = dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetPlaneNumber();
00398         goldCoG2 = dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetCoG();
00399         goldZ2 = dynamic_cast<DmxPlane*>(goldenPlaneItr.Ptr())->GetZPosition();
00400       }
00401       ++goldCtr;
00402     }//end loop over golden planes
00403 
00404     //find the slope and intercept between the two golden plane centers of gravity
00405     Float_t slope = (goldCoG2 - goldCoG1)/(goldZ2 - goldZ1);
00406     Float_t intercept = goldCoG2 - goldZ2*slope;
00407     
00408     //MSG("DMX", Msg::kInfo) << "first gold plane " << goldCoG1 << " " << goldZ1 
00409     //             << " second " << goldCoG2 << " " << goldZ2 << " " << slope << " " << intercept << endl; 
00410 
00411     //slice the lead plane iterator to be between the golden planes.
00412     //if first loop, slice from vertex to second golden plane
00413     //if last loop, slice from first golden plane to end plane
00414     
00415     if(loopCtr == 0) firstPlane = vertexPlane;
00416     else if(loopCtr == numLoops-1) lastPlane = endPlane;
00417 
00418     leadPlaneItr.GetSet()->Slice(firstPlane, lastPlane);
00419 
00420     SetPlanesToDeterminedFit(leadPlaneItr, intercept, slope, firstPlane, lastPlane);
00421 
00422     fStrayPlanes += util.FindPlanesOffFit(leadPlaneItr, fStrayCut);
00423 
00424     //clear the lead plane slice
00425     leadPlaneItr.GetSet()->Slice();
00426 
00427     ++loopCtr;
00428   }//end loop over sets of goldenPlanes
00429 
00430   leadPlaneItr.Reset();
00431   goldenPlaneItr.Reset();
00432   return;
00433 
00434 }

Member Data Documentation

UInt_t AlgDeMuxGolden::fPlanesInSet [private]

Definition at line 39 of file AlgDeMuxGolden.h.

Referenced by RunAlg().

Int_t AlgDeMuxGolden::fStrayCut [private]

Definition at line 40 of file AlgDeMuxGolden.h.

Referenced by UseGoldenFit().

Int_t AlgDeMuxGolden::fStrayPlanes [private]

Definition at line 38 of file AlgDeMuxGolden.h.

Referenced by RunAlg(), and UseGoldenFit().

The documentation for this class was generated from the following files:
Generated on Wed May 22 22:00:49 2013 for loon by  doxygen 1.4.7