StRoot: StStrangeMuDstMaker/StDecayMode.cc Source File

00001 /***********************************************************************
00002  *
00003  * StDecayMode: Singleton class for determining decay modes of MC vertices
00004  *
00005  ***********************************************************************/
00006 
00007 #include "StDecayMode.hh"
00008 
00009 StDecayMode* StDecayMode::mInstance = 0;
00010 
00011 ClassImp(StDecayMode)
00012 
00013 StDecayMode::StDecayMode() {}
00014 StDecayMode::~StDecayMode() { mInstance=0; }
00015 
00016 /*
00017     Note: reassign Gamma Id to -1
00018     Decay          Daughter Geant IDs    |sum / product (miss Nu)
00019 11  K+ -> Mu Nu        5 4  (mu+)        |9 / 20 (5 / 5)
00020     K+ -> Pi+ Pi0      8 7               |15 / 56
00021     K+ -> Pi+ Pi+ Pi-  8 8 9             |25 / 576
00022     K+ -> E Nu Pi0     2 4 7 (positron)  |13 / 42 (9 / 14)
00023     K+ -> Mu Nu Pi0    5 4 7             |16 / 140 (12 / 35)
00024     K+ -> Pi+ Pi0 Pi0  8 7 7             |22 / 392
00025 
00026 12  K- -> Mu Nu        6 4   (mu-)       |10 / 24 (6 / 6)
00027     K- -> Pi- Pi0      9 7               |16 / 63
00028     K- -> Pi+ Pi- Pi-  8 9 9             |26 / 648
00029     K- -> E Nu Pi0     3 4 7             |14 / 84 (10 / 21)
00030     K- -> Mu Nu Pi0    6 4 7             |17 / 168 (13 / 42)
00031     K- -> Pi- Pi0 Pi0  9 7 7             |23 / 441
00032 
00033 10  K0l -> Pi+ Pi- Pi0    8 9 7          |24 / 504
00034     K0l -> Pi+ Mu- Nu     8 6 4          |18 / 192 (14 / 48)
00035     K0l -> Pi- Mu+ Nu     9 5 4          |18 / 180 (14 / 45)
00036     K0l -> Pi+ E Nu       8 3 4          |15 / 96 (11 / 24)
00037     K0l -> Pi- E Nu       9 2 4          |15 / 72 (11 / 18)
00038     K0l -> Pi- Pi+        9 8            |17 / 72
00039 
00040 16  K0s -> Pi+ Pi-        8  9           |17 
00041     K0s -> Pi0 Pi0        7  7           |14
00042     K0s -> Pi+ Pi- Gamma  8  9 -1        |16 
00043 
00044  8  Pi+ -> Mu+ Nu         5 4            |9 (5)
00045     Pi+ -> Mu+ Nu Gamma   5 4 -1         |8 (4)
00046     Pi+ -> E Nu           2 4            |6 (2)
00047 
00048  9  Pi- -> Mu- Nu         6 4            |10 (6)
00049     Pi- -> Mu- Nu Gamma   6 4 -1         |9 (5)
00050     Pi- -> E Nu           3 4            |7 (3)
00051 
00052  5  Mu+ -> E Nu Nu        2 4 4          |10 (2,6)
00053     Mu+ -> E Nu Nu Gamma  2 4 4 -1       |9 (1,5)
00054 
00055  6  Mu- -> E Nu Nu        3 4 4          |11 (3,7)
00056     Mu- -> E Nu Nu Gamma  3 4 4 -1       |10 (2,6)
00057 
00058 18  Lambda -> P Pi-        14 9          |23 / 126
00059     Lambda -> N Pi0        13 7          |20 / 91
00060     Lambda -> N Gamma      13 -1         |12 / -13
00061     Lambda -> P Pi- Gamma  14 9 -1       |22 / -126
00062     Lambda -> P E Nu       14 3 4        |21 / 168 (17 / 42)
00063     Lambda -> P Mu Nu      14 6 4        |20 / 336 (24 / 84)
00064 
00065 26  AntiLambda -> AntiP Pi+        15 8       |23 / 120
00066     AntiLambda -> AntiN Pi0        25 7       |32 / 175
00067     AntiLambda -> AntiN Gamma      25 -1      |24 / -25
00068     AntiLambda -> AntiP Pi+ Gamma  15 8 -1    |22 / -120
00069     AntiLambda -> AntiP E Nu       15 2 4     |21 / 120 (17 / 30)
00070     AntiLambda -> AntiP Mu Nu      15 5 4     |20 / 300 (24 / 75)
00071 
00072 23  Cascade     -> Lambda Pi-      18  9      |27  
00073 31  AntiCascade -> AntiLambda Pi+  26  8      |34
00074 
00075 24  Omega       -> Lambda  K-      18 12      |30 / 216
00076     Omega       -> Cascade0 Pi-    22  9      |31 / 198
00077     Omega       -> Cascade- Pi0    23  7      |30 / 161
00078 
00079 32  AntiOmega   -> AntiLambda K+    26 11      |37 / 286
00080     AntiOmega   -> AntiCascade0 Pi+ 30  8      |38 / 240
00081     AntiOmega   -> AntiCascade+ Pi0 31  7      |38 / 217
00082  */
00083 //____________________________________________________________________
00084 Int_t StDecayMode::Process(StMcVertex* mcVertex) 
00085 {
00086   Int_t ID = 0;
00087   Int_t ID2 = 1;
00088   if (!mcVertex->numberOfDaughters()) return kWrongDecay;    // No daughters
00089   if (mcVertex->geantProcess() != 5) return kWrongDecay; // Not a weak decay
00090   const StMcTrack* parent = mcVertex->parent(); 
00091 
00092   if (parent) {
00093     Int_t parentId = parent->geantId();
00094     StSPtrVecMcTrack& Daughters = mcVertex->daughters();
00095 
00096     for (StMcTrackIterator DTrackIt = Daughters.begin();
00097                            DTrackIt != Daughters.end(); DTrackIt++) {
00098       Int_t daughterId = (*DTrackIt)->geantId();
00099       if (daughterId == 1) daughterId = -1;
00100       ID  += daughterId;
00101       ID2 *= daughterId;
00102     }
00103     
00104     switch (parentId) {
00105       case ( 5) : return MuPlusProcess(ID);
00106       case ( 6) : return MuMinusProcess(ID);
00107       case ( 8) : return PiPlusProcess(ID);
00108       case ( 9) : return PiMinusProcess(ID);
00109       case (10) : return KLongProcess(ID,ID2);
00110       case (11) : return KPlusProcess(ID2);
00111       case (12) : return KMinusProcess(ID2);
00112       case (16) : return KShortProcess(ID);
00113       case (18) : return LambdaProcess(ID2);
00114       case (26) : return AntiLambdaProcess(ID,ID2);
00115       case (23) : return XiProcess(ID);
00116       case (31) : return AntiXiProcess(ID);
00117       case (24) : return OmegaProcess(ID2);
00118       case (32) : return AntiOmegaProcess(ID2);
00119     }
00120   }     
00121   return kWrongDecay;
00122 }
00123 //_____________________________________________________________________  
00124 Int_t StDecayMode::KPlusProcess(Int_t ID) 
00125 {
00126   switch (ID) {
00127     case (  5) :
00128     case ( 20) : return kKPlus2MuNu;
00129     case ( 56) : return kKPlus2PiPlusPiZero;
00130     case (576) : return kKPlus2PiPlusPiPlusPiMinus;
00131     case ( 14) :
00132     case ( 42) : return kKPlus2ENuPiZero;
00133     case ( 35) :
00134     case (140) : return kKPlus2MuNuPiZero;
00135     case (392) : return kKPlus2PiPlusPiZeroPiZero;
00136     default    : return kWrongDecay;
00137   }
00138 }
00139 
00140 //_____________________________________________________________________
00141 Int_t StDecayMode::KMinusProcess(Int_t ID)  
00142 {
00143   switch (ID) {
00144     case (  6) :
00145     case ( 24) : return kKMinus2MuNu;
00146     case ( 63) : return kKMinus2PiMinusPiZero;
00147     case (648) : return kKMinus2PiPlusPiMinusPiMinus;
00148     case ( 21) :
00149     case ( 84) : return kKMinus2ENuPiZero;
00150     case ( 42) :
00151     case (168) : return kKMinus2MuNuPiZero;
00152     case (441) : return kKMinus2PiMinusPiZeroPiZero;
00153     default    : return kWrongDecay;
00154   }
00155 }
00156 
00157 //____________________________________________________________________
00158 Int_t StDecayMode::KLongProcess(Int_t ID, Int_t ID2)
00159 {
00160   switch (ID2) {
00161     case (504) : return kKLong2PiPlusPiZeroPiMinus;
00162     case ( 48) :
00163     case (192) : return kKLong2PiPlusMuNu;
00164     case ( 45) :
00165     case (180) : return kKLong2PiMinusMuNu;
00166     case ( 24) :
00167     case ( 96) : return kKLong2PiMinusENu;
00168     case ( 18) :
00169     case ( 72) : return ((ID==17) ? kKLong2PiPlusPiMinus : kKLong2PiPlusENu);
00170     default    : return kWrongDecay;
00171   }
00172 }
00173 //____________________________________________________________________
00174 Int_t StDecayMode::KShortProcess(Int_t ID)
00175 {
00176   switch (ID) {
00177     case (17) : return kKShort2PiPlusPiMinus;
00178     case (14) : return kKShort2PiZeroPiZero;
00179     case (16) : return kKShort2PiPlusPiMinusGamma;
00180     default   : return kWrongDecay;
00181   }
00182 }
00183 //____________________________________________________________________
00184 Int_t StDecayMode::PiPlusProcess(Int_t ID)
00185 {
00186   switch (ID) {
00187     case (5) :
00188     case (9) : return kPiPlus2MuNu;
00189     case (4) :
00190     case (8) : return kPiPlus2MuNuGamma;
00191     case (2) :
00192     case (6) : return kPiPlus2ENu;
00193     default  : return kWrongDecay;
00194   }
00195 }
00196 //____________________________________________________________________
00197 Int_t StDecayMode::PiMinusProcess(Int_t ID)
00198 {
00199   switch (ID) {
00200     case ( 6) :
00201     case (10) : return kPiMinus2MuNu;
00202     case ( 5) :
00203     case ( 9) : return kPiMinus2MuNuGamma;
00204     case ( 3) :
00205     case ( 7) : return kPiMinus2ENu;
00206     default   : return kWrongDecay;
00207   }
00208 }
00209 //____________________________________________________________________
00210 Int_t StDecayMode::MuPlusProcess(Int_t ID)
00211 {
00212   switch (ID) {
00213     case ( 2) :
00214     case ( 6) :
00215     case (10) : return kMuPlus2ENuNu;
00216     case ( 1) :
00217     case ( 5) :
00218     case ( 9) : return kMuPlus2ENuNuGamma;
00219     default   : return kWrongDecay;
00220   }
00221 }
00222 //____________________________________________________________________
00223 Int_t StDecayMode::MuMinusProcess(Int_t ID)
00224 {
00225   switch (ID) {
00226     case ( 3) :
00227     case ( 7) :
00228     case (11) : return kMuMinus2ENuNu;
00229     case ( 2) :
00230     case ( 6) :
00231     case (10) : return kMuMinus2ENuNuGamma;
00232     default   : return kWrongDecay;
00233   }
00234 }
00235 //____________________________________________________________________
00236 Int_t StDecayMode::LambdaProcess(Int_t ID) 
00237 {
00238   switch (ID) {
00239     case ( 126) : return kLambda2ProtonPiMinus;
00240     case (  91) : return kLambda2NeutronPiZero;
00241     case ( -13) : return kLambda2NeutronGamma;
00242     case (-126) : return kLambda2ProtonPiMinusGamma;
00243     case (  42) :
00244     case ( 168) : return kLambda2ProtonENu;
00245     case (  84) :
00246     case ( 336) : return kLambda2ProtonMuNu;
00247     default     : return kWrongDecay;
00248   }
00249 }
00250 //____________________________________________________________________
00251 Int_t StDecayMode::AntiLambdaProcess(Int_t ID, Int_t ID2) 
00252 {
00253   switch (ID) {
00254     case (23) : return kAntiLambda2AntiProtonPiPlus;
00255     case (32) : return kAntiLambda2AntiNeutronPiZero;
00256     case (24) : if (ID2 == -25) return kAntiLambda2AntiNeutronGamma;
00257     case (20) : return kAntiLambda2AntiProtonMuNu;
00258     case (22) : return kAntiLambda2AntiProtonPiPlusGamma;
00259     case (17) :
00260     case (21) : return kAntiLambda2AntiProtonENu;
00261     default   : return kWrongDecay;
00262   }
00263 }
00264 
00265 //____________________________________________________________________
00266 Int_t StDecayMode::XiProcess(Int_t ID) 
00267 {
00268   switch (ID) {
00269     case (27) : return kCascade2LambdaPiMinus;
00270     default   : return kWrongDecay;
00271   }
00272 }
00273 //____________________________________________________________________
00274 Int_t StDecayMode::AntiXiProcess(Int_t ID) 
00275 {
00276   switch (ID) {
00277     case (34) : return kAntiCascade2AntiLambdaPiPlus;
00278     default   : return kWrongDecay;
00279   }
00280 }
00281 //____________________________________________________________________
00282 Int_t StDecayMode::OmegaProcess(Int_t ID) 
00283 {
00284   switch (ID) {
00285     case (216) : return kOmega2LambdaKMinus;
00286     case (198) : return kOmega2CascadePiMinus;
00287     case (161) : return kOmega2CascadePiZero;
00288     default    : return kWrongDecay;
00289   }
00290 }
00291 //____________________________________________________________________
00292 Int_t StDecayMode::AntiOmegaProcess(Int_t ID) 
00293 {
00294   switch (ID) {
00295     case (286) : return kAntiOmega2AntiLambdaKPlus;
00296     case (240) : return kAntiOmega2AntiCascadePiPlus;
00297     case (217) : return kAntiOmega2AntiCascadePiZero;
00298     default    : return kWrongDecay;
00299   }
00300 }
00301 //____________________________________________________________________
00302 Int_t StDecayMode::ParentCharge(Int_t mode) 
00303 {
00304   switch (mode) {
00305     case (kKPlus2MuNu)                   :
00306     case (kKPlus2PiPlusPiZero)           :
00307     case (kKPlus2PiPlusPiPlusPiMinus)    :
00308     case (kKPlus2ENuPiZero)              :
00309     case (kKPlus2MuNuPiZero)             :
00310     case (kKPlus2PiPlusPiZeroPiZero)     :
00311     case (kPiPlus2MuNu)                  :
00312     case (kPiPlus2MuNuGamma)             :
00313     case (kPiPlus2ENu)                   :
00314     case (kMuPlus2ENuNu)                 :
00315     case (kMuPlus2ENuNuGamma)            :
00316     case (kAntiCascade2AntiLambdaPiPlus) :
00317     case (kAntiOmega2AntiLambdaKPlus)    :
00318     case (kAntiOmega2AntiCascadePiPlus)  : return  1;
00319     case (kKMinus2MuNu)                  :
00320     case (kKMinus2PiMinusPiZero)         :
00321     case (kKMinus2PiPlusPiMinusPiMinus)  :
00322     case (kKMinus2ENuPiZero)             :
00323     case (kKMinus2MuNuPiZero)            :
00324     case (kKMinus2PiMinusPiZeroPiZero)   :
00325     case (kPiMinus2MuNu)                 :
00326     case (kPiMinus2MuNuGamma)            :
00327     case (kPiMinus2ENu)                  :
00328     case (kMuMinus2ENuNu)                :
00329     case (kMuMinus2ENuNuGamma)           :
00330     case (kCascade2LambdaPiMinus)        :
00331     case (kOmega2LambdaKMinus)           :
00332     case (kOmega2CascadePiMinus)         : return -1;
00333     default : {}
00334   }
00335   return 0;                                                                     
00336 }
00337 
00338 /***********************************************************************
00339  *
00340  * $Id: StDecayMode.cc,v 3.3 2002/05/29 19:08:44 genevb Exp $
00341  * $Log: StDecayMode.cc,v $
00342  * Revision 3.3  2002/05/29 19:08:44  genevb
00343  * Some additional decay modes, better organization of enumeration
00344  *
00345  * Revision 3.2  2002/04/30 16:02:47  genevb
00346  * Common muDst, improved MC code, better kinks, StrangeCuts now a branch
00347  *
00348  * Revision 3.1  2001/05/04 20:15:13  genevb
00349  * Common interfaces and reorganization of components, add MC event info
00350  *
00351  * Revision 3.0  2000/07/14 12:56:47  genevb
00352  * Revision 3 has event multiplicities and dedx information for vertex tracks
00353  *
00354  * Revision 2.1  2000/06/09 22:17:09  genevb
00355  * Allow MC data to be copied between DSTs, other small improvements
00356  *
00357  * Revision 2.0  2000/06/05 05:19:37  genevb
00358  * New version of Strangeness micro DST package
00359  *
00360  ***********************************************************************/