// // ******************************************************************** // * License and Disclaimer * // * * // * The Geant4 software is copyright of the Copyright Holders of * // * the Geant4 Collaboration. It is provided under the terms and * // * conditions of the Geant4 Software License, included in the file * // * LICENSE and available at http://cern.ch/geant4/license . These * // * include a list of copyright holders. * // * * // * Neither the authors of this software system, nor their employing * // * institutes,nor the agencies providing financial support for this * // * work make any representation or warranty, express or implied, * // * regarding this software system or assume any liability for its * // * use. Please see the license in the file LICENSE and URL above * // * for the full disclaimer and the limitation of liability. * // * * // * This code implementation is the result of the scientific and * // * technical work of the GEANT4 collaboration. * // * By using, copying, modifying or distributing the software (or * // * any work based on the software) you agree to acknowledge its * // * use in resulting scientific publications, and indicate your * // * acceptance of all terms of the Geant4 Software license. * // ******************************************************************** // // // $Id: G4QIsotope.hh,v 1.6 2009/02/23 09:49:24 mkossov Exp $ // GEANT4 tag $Name: hadr-chips-V09-03-08 $ // // // ---------------- G4QIsotope header ---------------- // by Mikhail Kossov, December 2003. // Header of the G4QIsotope class of the CHIPS Simulation Branch in GEANT4 // ---------------------------------------------------------------------------- // Short descriptionIt contains information about natural abundances of stable // and long living isotopes and a NEW "Element" can be initialised for any // isotope set. Randomization of isotopes of the Natural Elements is hardwired and // fast randomization of isotopes of the user defined Elements is a bit slower // CrossSectionWeighted randomisation of isotopes is slow (same for Nat and New) // ------------------------------------------------------------------------------- // **************************************************************************************** // ********* This HEADER is temporary moved from the photolepton_hadron directory ********* // ******* DO NOT MAKE ANY CHANGE! With time it'll move back to photolepton...(M.K.) ****** // **************************************************************************************** // // 1 2 3 4 5 6 7 8 9 //34567890123456789012345678901234567890123456789012345678901234567890123456789012345678901 #ifndef G4QIsotope_hh #define G4QIsotope_hh #include "globals.hh" #include "G4ios.hh" #include "Randomize.hh" #include class G4QIsotope { protected: G4QIsotope(); // ***Singletone*** All natural elements are initialized in the Constructor public: ~G4QIsotope(); // It's public for compilation purposes on Windows, user must not call it! // Create newElement with the Abundancy vector (User must delete elements of the vector) // =================================================================================== // Example of initialization of the new Element (not natural abanduncies of isotopes): // ======= --------------------------------------------------------------------------- //std::vector >*a= new std::vector >; // a->push_back(std::make_pair(n1,abundancy1)); // a->push_back(std::make_pair(n2,abundancy2)); // a->push_back(std::make_pair(n3,1.-abundancy1-abundancy2)); ////Sum of abundancies must be 1 otherwise worning appears & if less LastAbu is increased // G4int Z=1; // G4int ind=1; //// ind>0, if =<0 then the InitElement member function returns the first free index !!! //// For ind>0 if theIndex "ind" already exists, returns the first free index !!! // G4QIsotope::Get()->InitElement(Z, ind, a); // G4QIsotope class is a Singletone // std::for_each(a->begin(), a->end(), void operator()(std::pair >* P) // {delete P;}); //// OR just use the following, which is faser: //// G4int nA=a->size(); //// if(nA) for(G4int i=0; ioperator[](i);} // G4int InitElement(G4int Z, G4int index, std::vector*>* abund); // The highest index defined for Element with Z (Index>0 correspondToUserDefinedElements) // ====================================================================================== G4int GetLastIndex(G4int Z); // Returns the last defined index (if only natural: =0) // Indices can have differen numbers (not 1,2,3,...) & in different sequences (9,3,7,...) // ====================================================================================== G4bool IsDefined(G4int Z, G4int Ind); // Returns true if defined, false, if not defined // A#ofNeutrons in Element with Z & UseDefIndex. Universal for Nat(index=0) & UserDefElem // ====================================================================================== G4int GetNeutrons(G4int Z, G4int index=0);//If theElement doesn't exist, returns negative // #ofProtons in stable isotopes with fixed A=Z+N. Returns length and fils VectOfIsotopes // ====================================================================================== // Example of printing of isotopes with A=152: // ======= ----------------------------------- // G4int A=152; // A can not be more than 269 // std::vector isV(4); // At present A with nIso>4 are not known // G4int nIso= G4QIsotope::Get()->GetProtons(A, isV); // isV is cleaned up before filling // if(nIso)for(G4int i,i& isoV); // Get a pointer to the vector of pairs(N,CrosS), where N is used to calculate CrosS // ================================================================================= // Example of initialization of the Cross Section to randomize weighted isotopes: // ======= ---------------------------------------------------------------------- // std::vector*>* cs= G4QIsotope::Get()->GetCSVector(Z, index); // G4int nIs=cs->size; // A#Of Isotopes in the element // if(nIs) for(G4int i; iat(i)->first; // A#Of neuterons in the isotope // cs->at(i)->second = CalculateCrossSection(particle,Z,N)// Calc particle+A(Z,N) CrosS // } std::vector*>* GetCSVector(G4int Z, G4int index = 0); // Get the abundancy vector for calculation of mean cross sections std::vector*>* GetAbuVector(G4int Z, G4int index = 0); // Get the summed abundancy vector (e.g. for randomization by itself) std::vector*>* GetSumAVector(G4int Z, G4int index = 0); // Calculates the mean Cross Section for the initialized Element(ind=0 Nat,ind>0 UserDef) G4double GetMeanCrossSection(G4int Z, G4int index = 0); // IsoCS's must init, IfNotRet<0 // Randomize A#OfNeutrons in the Isotope weighted by theAbubdancies and theCrossSections G4int GetCSNeutrons(G4int Z, G4int index = 0); // IsoCrosSections must init, IfNotRet<0 static G4QIsotope* Get(); // Get a pointer to the Singletone G4QIsotope private: G4int RandomizeNeutrons(G4int Z); // Gives a#of neutrons in the Random Isotope for the Z private: // Initialized in the constructor static std::vector*>*> natElements; //NaturalElem's static std::vector*>*> natSumAbund; //NatElemSumA's static std::vector*>*> natIsoCrosS; //CSOfNatElem's // It is initialized by user, but it is cleaned up in the destructor static std::vector*>*>*> newElems; static std::vector*>*>*> newSumAb; static std::vector*>*>*> newIsoCS; }; #endif