StRoot: StMuDSTMaker/COMMON/StMuHelix.cxx Source File

00001 /***************************************************************************
00002  *   
00003  * $Id: StMuHelix.cxx,v 1.5 2009/12/01 21:56:35 tone421 Exp $
00004  * Author: Frank Laue, BNL, laue@bnl.gov
00005  *
00006  ***************************************************************************/
00007 
00008 #include "StMuHelix.h"
00009 
00010 StMuHelix::StMuHelix(const StPhysicalHelixD &hh, double field) {
00011   mP = hh.momentum(field*kilogauss);
00012   mOrigin = hh.origin();
00013   mQ = hh.charge( field*kilogauss );
00014   mB = field;
00015 }
00016 
00017 int StMuHelix::bad() const 
00018 {
00019    double qwe;
00020    if(!mOrigin.valid())         return 10;
00021    if(!mP.valid     ())         return 20;
00022 
00023    qwe = fabs(mP.x())+fabs(mP.y());
00024    if (qwe <=1.e-10)            return 21;
00025 
00026    if (abs(mQ)  >    10)        return 30;
00027 
00028    if (!::finite(mB))           return 40;
00029    if (fabs(mB) >  1000)        return 41;
00030    return 0;
00031 }
00032 StPhysicalHelix  StMuHelix::helix() const
00033 {
00034    return StPhysicalHelix(mP,mOrigin,mB*kilogauss,mQ);
00035 }
00036 
00037 ClassImp(StMuHelix)
00038 /***************************************************************************
00039  *
00040  * $Log: StMuHelix.cxx,v $
00041  * Revision 1.5  2009/12/01 21:56:35  tone421
00042  * Implemented changes as per http://www.star.bnl.gov/rt2/Ticket/Display.html?id=1734
00043  *
00044  * Revision 1.4  2009/07/07 19:37:58  perev
00045  * helix() method added
00046  *
00047  * Revision 1.3  2003/10/30 20:08:13  perev
00048  * Check of quality added
00049  *
00050  * Revision 1.2  2003/10/28 18:57:56  perev
00051  * BadData protection added
00052  *
00053  * Revision 1.1  2002/03/08 17:04:18  laue
00054  * initial revision
00055  *
00056  *
00057  **************************************************************************/