updated STARTLIGHT to r191 (http://starlight.hepforge.org/svn/trunk)
authorcmayer <Christoph.Mayer@cern.ch>
Wed, 12 Nov 2014 11:08:57 +0000 (12:08 +0100)
committercmayer <Christoph.Mayer@cern.ch>
Wed, 12 Nov 2014 11:11:48 +0000 (12:11 +0100)
coverity fixes (warnings due to -Wshadow)

STARLIGHT/starlight/TStarLight/testsl.C
STARLIGHT/starlight/include/inputParameters.h
STARLIGHT/starlight/include/nucleus.h
STARLIGHT/starlight/src/incoherentPhotonNucleusLuminosity.cpp
STARLIGHT/starlight/src/inputParser.cpp
STARLIGHT/starlight/src/photonNucleusCrossSection.cpp
STARLIGHT/starlight/src/wideResonanceCrossSection.cpp
STARLIGHT/starlight/standard/slight.in
STARLIGHT/test/sim.C

index 104ff8f..1e70384 100644 (file)
@@ -10,12 +10,21 @@ void testsl() {
 
   TStarLight* sl = new TStarLight("starlight generator", "title", "");
 
+#if 1
   sl->SetParameter("BEAM_1_Z     =    1    #Z of projectile");
   sl->SetParameter("BEAM_1_A     =    1    #A of projectile");
   sl->SetParameter("BEAM_2_Z     =   82    #Z of target");
   sl->SetParameter("BEAM_2_A     =  208    #A of target");
   sl->SetParameter("BEAM_1_GAMMA = 4264.4  #Gamma of the colliding ions");
   sl->SetParameter("BEAM_2_GAMMA = 1682.4  #Gamma of the colliding ions");
+#else
+  sl->SetParameter("BEAM_1_Z     =   82    #Z of projectile");
+  sl->SetParameter("BEAM_1_A     =  208    #A of projectile");
+  sl->SetParameter("BEAM_2_Z     =   82    #Z of target");
+  sl->SetParameter("BEAM_2_A     =  208    #A of target");
+  sl->SetParameter("BEAM_2_GAMMA = 1470  #Gamma of the colliding ions");
+  sl->SetParameter("BEAM_1_GAMMA = 1470  #Gamma of the colliding ions");
+#endif
   sl->SetParameter("W_MAX        =   15    #Max value of w");
   sl->SetParameter("W_MIN        =    -1   #Min value of w");
   sl->SetParameter("W_N_BINS     =   40    #Bins i w");
@@ -27,7 +36,7 @@ void testsl() {
   sl->SetParameter("CUT_ETA      =    0    #Cut in pseudorapidity? (0 = no, 1 = yes)");
   sl->SetParameter("ETA_MIN      =  -10    #Minimum pseudorapidity");
   sl->SetParameter("ETA_MAX      =   10    #Maximum pseudorapidity");
-  sl->SetParameter("PROD_MODE    =    2    #gg or gP switch (1 = 2-photon, 2 = coherent vector meson (narrow), 3 = coherent vector meson (wide), # 4 = incoherent vector meson, 5 = A+A DPMJet single, 6 = A+A DPMJet double, 7 = p+A DPMJet single, 8 = p+A Pythia single )");
+  sl->SetParameter("PROD_MODE    =    4    #gg or gP switch (1 = 2-photon, 2 = coherent vector meson (narrow), 3 = coherent vector meson (wide), # 4 = incoherent vector meson, 5 = A+A DPMJet single, 6 = A+A DPMJet double, 7 = p+A DPMJet single, 8 = p+A Pythia single )");
   sl->SetParameter("PROD_PID     =  113    #Channel of interest (not relevant for photonuclear processes)");
   sl->SetParameter("RND_SEED     = 34533   #Random number seed");
   sl->SetParameter("BREAKUP_MODE  =   5    #Controls the nuclear breakup");
@@ -40,7 +49,7 @@ void testsl() {
 
   sl->InitStarLight();
   sl->PrintInputs(std::cout);
-  TClonesArray tca("TParticle", 1000);
+  TClonesArray tca("TParticle", 100);
 
   TLorentzVector v[2], vSum;
   TH1* hM  = new TH1D("hM",  "STARLIGHT;M#(){#pi^{+}#pi^{-}}",     100, 0., 2.);
@@ -49,7 +58,7 @@ void testsl() {
 
   std::ofstream ofs("rho0.txt");
   TParticle *p;
-  for (Int_t counter(0); counter<100;) {
+  for (Int_t counter(0); counter<1000;) {
     sl->GenerateEvent();    
     sl->BoostEvent();    
     sl->ImportParticles(&tca, "ALL");
index 6264915..fd11405 100644 (file)
@@ -108,7 +108,12 @@ class parameter : public parameterbase
 public:
 
     // Constructor
-    parameter(const std::string &name, T value, bool required = true) :parameterbase(),_name(name), _value(value), _validate(validate), _required(required) {}
+    parameter(const std::string &name_, T value_, bool required_ = true)
+      : parameterbase()
+      , _name(name_)
+      , _value(value_)
+      , _validate(validate)
+      , _required(required_) {}
 
     virtual ~parameter() {}
 //     T operator()() const {
@@ -128,7 +133,7 @@ public:
     
     void setValue(T v) { _value = v; }
     
-    void setName(std::string name) { _name = name; }
+    void setName(std::string newName) { _name = newName; }
     
     void setRequired(bool r) { _required = r; }
     
index b5edf04..b995b90 100644 (file)
@@ -70,7 +70,10 @@ private:
        double fritiofR0         () const { return _r0 * pow(_A, (1./3.)); }  ///< Fritiof _r0 (rws)/formfactor
 
        double rws(const double r) const
-       { return 1.0 / (1. + exp((r - fritiofR0()) / woodSaxonSkinDepth())); } ///< Wood-Saxon nuclear density
+       { //std::cout << "rws: " << r << " " << fritiofR0() << " " << woodSaxonSkinDepth() << std::endl;
+         const double x=exp(-(r - fritiofR0()) / woodSaxonSkinDepth());
+         return x/(1.+x);
+          return 1.0 / (1. + exp((r - fritiofR0()) / woodSaxonSkinDepth())); } ///< Wood-Saxon nuclear density
 
        int    _Z;                      ///< atomic number of nucleus
        int    _A;                      ///< nucleon number of nucleus
index a4adbe1..d092734 100644 (file)
@@ -84,7 +84,7 @@ void incoherentPhotonNucleusLuminosity::incoherentPhotonNucleusDifferentialLumin
     
   // Write the values of W used in the calculation to slight.txt.  
   wylumfile.open("slight.txt");
-  wylumfile << inputParametersInstance.parameterValueKey() << endl;
+  // wylumfile << inputParametersInstance.parameterValueKey() << endl;
   wylumfile << getbbs().beam1().Z() <<endl;
   wylumfile << getbbs().beam1().A() <<endl;
   wylumfile << getbbs().beam2().Z() <<endl;
@@ -180,6 +180,7 @@ void incoherentPhotonNucleusLuminosity::incoherentPhotonNucleusDifferentialLumin
       }
 
       wylumfile << dndWdY << endl;
+
     }
   }
 
@@ -187,11 +188,8 @@ void incoherentPhotonNucleusLuminosity::incoherentPhotonNucleusDifferentialLumin
   wylumfile << inputParametersInstance.parameterValueKey() << endl;
   wylumfile.close();
   
-//   wylumfile.open("slight.txt",ios::app);
-  cout << "bwnorm: "<< bwnorm <<endl;
-//   wylumfile << bwnorm << endl;
-//   wylumfile << inputParametersInstance.parameterValueKey() << endl;
-//   wylumfile.close();
+  // cout << "bwnorm: "<< bwnorm <<endl;
+
 }
 
 
index ed61607..fbc33df 100644 (file)
@@ -56,7 +56,7 @@ int inputParser::parseFile(std::string filename)
     int lineSize = 256;
     char tmp[lineSize];
     int nParameters = 0;
-    while (!infile.getline(tmp, lineSize).eof())
+    while (infile.getline(tmp, lineSize))
     {
 
         std::string line(tmp);
@@ -69,6 +69,8 @@ int inputParser::parseFile(std::string filename)
 }
 int inputParser::parseString(std::string str)
 {
+  if (str == "")
+    return false;
 
     std::string word;
     std::string name;
index e25fe26..8c281cf 100644 (file)
@@ -20,9 +20,9 @@
 ///////////////////////////////////////////////////////////////////////////
 //
 // File and Version Information:
-// $Rev:: 169                         $: revision of last commit
-// $Author:: jnystrand                $: author of last commit
-// $Date:: 2014-03-28 17:17:37 +0100 #$: date of last commit
+// $Rev:: 189                         $: revision of last commit
+// $Author:: srklein                  $: author of last commit
+// $Date:: 2014-10-27 04:29:14 +0100 #$: date of last commit
 //
 // Description:
 //
@@ -851,7 +851,8 @@ photonNucleusCrossSection::breitWigner(const double W,
        // width depends on energy - Jackson Eq. A.2
        // if below threshold, then return 0.  Added 5/3/2001 SRK
        // 0.5% extra added for safety margin
-       if( _particleType==RHO ||_particleType==RHOZEUS){  
+        // omega added here 10/26/2014 SRK
+       if( _particleType==RHO ||_particleType==RHOZEUS || _particleType==OMEGA){  
                if (W < 2.01*pionChargedMass){
                        nrbw_r=0.;
                        return nrbw_r;
index 3a43ffa..328f320 100644 (file)
@@ -20,9 +20,9 @@
 ///////////////////////////////////////////////////////////////////////////
 //
 // File and Version Information:
-// $Rev:: 164                         $: revision of last commit
-// $Author:: odjuvsla                 $: author of last commit
-// $Date:: 2013-10-06 16:18:08 +0200 #$: date of last commit
+// $Rev:: 191                         $: revision of last commit
+// $Author:: jnystrand                $: author of last commit
+// $Date:: 2014-11-11 11:51:05 +0100 #$: date of last commit
 //
 // Description:
 //
@@ -130,17 +130,9 @@ wideResonanceCrossSection::crossSectionCalculation(const double bwnormsave)
       
                        if(ega1 < Eth) continue;
                        if(ega2 > maxPhotonEnergy()) continue;
-                       // check it !!
-          
-                       if(J == 0){
-                               // >> 1st Point (Calculated only the first time)     =====>>>
-                               //ega1 used.                                                        
-                               csgA1=getcsgA(ega1,W);
-                       }
-                       else{
-                               csgA1 = csgA2;
-                       }
           
+                       csgA1=getcsgA(ega1,W);
+
                        //         >> Middle Point                      =====>>>
                        csgA12=getcsgA(ega12,W);         
 
index 5e743b5..df0bc7c 100644 (file)
@@ -1,32 +1,32 @@
-    BEAM_1_Z     =   1    #Z of projectile
-    BEAM_1_A     =  1    #A of projectile
-    BEAM_2_Z     =   82    #Z of target
-    BEAM_2_A     =  208    #A of target
-    BEAM_1_GAMMA = 4264.4    #Gamma of the colliding ions
-    BEAM_2_GAMMA = 1682.4    #Gamma of the colliding ions
-    W_MAX        =   15.0  #Max value of w
-    W_MIN        =    -1  #Min value of w
-    W_N_BINS     =   50    #Bins i w
-    RAP_MAX      =    1.5   #max y
-    RAP_N_BINS   =   200    #Bins i y
-    CUT_PT       =    0    #Cut in pT? 0 = (no, 1 = yes)
-    PT_MIN       =    1.0  #Minimum pT in GeV
-    PT_MAX       =    3.0  #Maximum pT in GeV
-    CUT_ETA      =    1    #Cut in pseudorapidity? (0 = no, 1 = yes)
-    ETA_MIN      =  -1.46504    #Minimum pseudorapidity
-    ETA_MAX      =   0.53496    #Maximum pseudorapidity
-    PROD_MODE    =   2     #gg or gP switch (1 = 2-photon, 2 = coherent vector meson (narrow), 3 = coherent vector meson (wide), # 4 = incoherent vector meson, 5 = A+A DPMJet single, 6 = A+A DPMJet double, 7 = p+A DPM  Jet single, 8 = p+A Pythia single )  // is N_EVENTS valid
-    N_EVENTS     =   10  #Number of events
-    PROD_PID     =  113    #Channel of interest (not relevant for photonuclear processes)
-    RND_SEED     = 5574533
-    OUTPUT_FORMAT =   2    #Form of the output
-    BREAKUP_MODE  =   5    #Controls the nuclear breakup
-    INTERFERENCE  =   0    #Interference (0 = off, 1 = on)
-    IF_STRENGTH   =   1.   #% of intefernce (0.0 - 0.1)
-   COHERENT      =   1    #Coherent=1,Incoherent=0
-    INCO_FACTOR   =   1.   #percentage of incoherence
-    BFORD         =   9.5  #
-    INT_PT_MAX    =   0.24 #Maximum pt considered, when interference is turned on
-    INT_PT_N_BINS = 120    #Number of pt bins when interference is turned on
-    N_THREADS = 1 # Number of threads used for calculating luminosity (when using the new method)
-    PYTHIA_FULL_EVENTRECORD = 0 # Write full pythia information to output (vertex, parents, daughter etc).
+BEAM_1_Z     =   1    #Z of projectile
+BEAM_1_A     =  1    #A of projectile
+BEAM_2_Z     =   82    #Z of target
+BEAM_2_A     =  208    #A of target
+BEAM_1_GAMMA = 4264.4    #Gamma of the colliding ions
+BEAM_2_GAMMA = 1682.4    #Gamma of the colliding ions
+W_MAX        =   15.0  #Max value of w
+W_MIN        =    -1  #Min value of w
+W_N_BINS     =   50    #Bins i w
+RAP_MAX      =    1.5   #max y
+RAP_N_BINS   =   200    #Bins i y
+CUT_PT       =    0    #Cut in pT? 0 = (no, 1 = yes)
+PT_MIN       =    1.0  #Minimum pT in GeV
+PT_MAX       =    3.0  #Maximum pT in GeV
+CUT_ETA      =    0    #Cut in pseudorapidity? (0 = no, 1 = yes)
+ETA_MIN      =  -1.46504    #Minimum pseudorapidity
+ETA_MAX      =   0.53496    #Maximum pseudorapidity
+PROD_MODE    =   4     #gg or gP switch (1 = 2-photon, 2 = coherent vector meson (narrow), 3 = coherent vector meson (wide), # 4 = incoherent vector meson, 5 = A+A DPMJet single, 6 = A+A DPMJet double, 7 = p+A DPM  Jet single, 8 = p+A Pythia single )
+N_EVENTS     =   10  #Number of events
+PROD_PID     =  113    #Channel of interest (not relevant for photonuclear processes)
+RND_SEED     = 5574533
+OUTPUT_FORMAT =   2    #Form of the output
+BREAKUP_MODE  =   5    #Controls the nuclear breakup
+INTERFERENCE  =   0    #Interference (0 = off, 1 = on)
+IF_STRENGTH   =   1.   #% of intefernce (0.0 - 0.1)
+COHERENT      =   1    #Coherent=1,Incoherent=0
+INCO_FACTOR   =   1.   #percentage of incoherence
+BFORD         =   9.5  #
+INT_PT_MAX    =   0.24 #Maximum pt considered, when interference is turned on
+INT_PT_N_BINS = 120    #Number of pt bins when interference is turned on
+N_THREADS = 1 # Number of threads used for calculating luminosity (when using the new method)
+PYTHIA_FULL_EVENTRECORD = 0 # Write full pythia information to output (vertex, parents, daughter etc).
index 185b968..352118c 100644 (file)
@@ -7,7 +7,7 @@ void sim(Int_t nev=100) {
   simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON ZDC PMD T0 VZERO FMD");
   simulator.SetMakeDigitsFromHits("ITS TPC");
 
-  simulator.SetDefaultStorage("alien://Folder=/alice/data/2010/OCDB");
+  simulator.SetDefaultStorage("alien://Folder=/alice/data/2013/OCDB");
 
 //   simulator.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
 //    simulator.SetSpecificStorage("GRP/GRP/Data",