Ignore:
Timestamp:
Jun 18, 2010, 11:42:07 AM (14 years ago)
Author:
garnier
Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

Location:
trunk/source/processes/electromagnetic/standard/test
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • trunk/source/processes/electromagnetic/standard/test/BremLPMTest.cc

    r1199 r1315  
    2525//
    2626// $Id: BremLPMTest.cc,v 1.2 2008/08/22 08:16:17 schaelic Exp $
    27 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2828//
    2929// ------------------------------------------------------------

  • trunk/source/processes/electromagnetic/standard/test/CrossSectionStd.cc

    r1199 r1315  
    2626//
    2727// $Id: CrossSectionStd.cc,v 1.5 2006/06/29 19:54:07 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2929//
    3030// ------------------------------------------------------------

  • trunk/source/processes/electromagnetic/standard/test/CrossSectionV52.cc

    r1199 r1315  
    2626//
    2727// $Id: CrossSectionV52.cc,v 1.4 2006/06/29 19:54:09 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2929//
    3030// ------------------------------------------------------------

  • trunk/source/processes/electromagnetic/standard/test/G4InitXscPAItest.cc

    r1199 r1315  
    2626//
    2727// $Id: G4InitXscPAItest.cc,v 1.6 2006/06/29 19:54:11 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2929//
    3030//

  • trunk/source/processes/electromagnetic/standard/test/G4PAIdNdxTest.cc

    r1199 r1315  
    2626//
    2727// $Id: G4PAIdNdxTest.cc,v 1.8 2008/04/14 14:53:39 grichine Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2929//
    3030//

  • trunk/source/processes/electromagnetic/standard/test/G4PAIxSectionTest.cc

    r1199 r1315  
    2525//
    2626//
    27 // $Id: G4PAIxSectionTest.cc,v 1.15 2006/06/29 19:54:15 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: G4PAIxSectionTest.cc,v 1.16 2009/12/30 12:57:41 grichine Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2929//
    3030//
    … …  
    432432 outFile <<k<<"\t"<< "  Material : " <<(*theMaterialTable)[k]->GetName() << G4endl ;
    433433  }
    434   G4String testName ;
     434  G4String testName = "Argon" ;
    435435  G4cout<<"Enter material name for test : "<<std::flush ;
    436436  //  G4cin>>testName ;
    … …  
    448448
    449449
    450   for(k=0;k<numOfMaterials;k++)
     450  for( k = 0; k < numOfMaterials; k++ )
    451451  {
    452     //    if((*theMaterialTable)[k]->GetName() != testName) continue ;
     452    if((*theMaterialTable)[k]->GetName() != testName) continue ;
    453453     outFile << "Material : " <<(*theMaterialTable)[k]->GetName() << G4endl ;
    454454     G4cout << "Material : " <<(*theMaterialTable)[k]->GetName() << G4endl ;
    … …  
    543543     //   G4PAIxSection testPAIproton(k,maxEnergyTransfer) ;
    544544
    545      kineticEnergy = 10.0*keV ;  // 110*MeV ;
     545     kineticEnergy = 100.*MeV;  // 10.0*keV ;  // 110*MeV ;
    546546
    547547     //     for(j=1;j<testPAIproton.GetNumberOfGammas();j++)

  • trunk/source/processes/electromagnetic/standard/test/PAIenergyLossTest.cc

    r1199 r1315  
    2626//
    2727// $Id: PAIenergyLossTest.cc,v 1.8 2006/06/29 19:54:17 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2929//
    3030//

  • trunk/source/processes/electromagnetic/standard/test/readBremLPM.py

    r1199 r1315  
    11# $Id: readBremLPM.py,v 1.2 2008/08/22 08:16:17 schaelic Exp $
    2 # GEANT4 tag $Name: geant4-09-03-cand-01 $
     2# GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    33#
    44# works together with BremLPMTest.cc

  • trunk/source/processes/electromagnetic/standard/test/test90Ne10CO2pai.cc

    r1199 r1315  
    2525//
    2626//
    27 // $Id: test90Ne10CO2pai.cc,v 1.6 2008/12/18 13:01:40 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03-cand-01 $
     27// $Id: test90Ne10CO2pai.cc,v 1.8 2009/12/30 12:57:41 grichine Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
    2929//
    3030//
    … …  
    3838// 11.04.08, V. Grichine test of PAI predictions for 90% Ne + 10% CO2
    3939//                       ALICE TPC gas mixture
     40// 8.12.09   V. Grichine update for T2K gas mixture
    4041
    4142#include "G4ios.hh"
    … …  
    150151int main()
    151152{
    152    std::ofstream outFile("90Ne10CO2pai.dat", std::ios::out );
     153  // std::ofstream outFile("90Ne10CO2pai.dat", std::ios::out );
     154   std::ofstream outFile("e5GeVt2kPhilippe.dat", std::ios::out );
    153155   outFile.setf( std::ios::scientific, std::ios::floatfield );
    154156
    … …  
    174176
    175177
     178  a = 1.01*g/mole;
     179  G4Isotope* ih1 = new G4Isotope("Hydrogen",iz=1,n=1,a);
     180
     181  a = 2.01*g/mole;
     182  G4Isotope* ih2 = new G4Isotope("Deuterium",iz=1,n=2,a);
     183
     184  G4Element* elH = new G4Element(name="Hydrogen",symbol="H",2);
     185  elH->AddIsotope(ih1,.999);
     186  elH->AddIsotope(ih2,.001);
     187
    176188  a = 12.01*g/mole;
    177189  G4Element* elC = new G4Element(name="Carbon",symbol="C", ez=6., a);
    … …  
    183195  G4Element* elO = new G4Element(name="Oxygen",symbol="O", ez=8., a);
    184196
     197 
     198  a = 19.00*g/mole;
     199  G4Element* elF  = new G4Element(name="Fluorine", symbol="F", z=9., a);
    185200 
     201  a = 39.948*g/mole;
     202  G4Element* elAr = new G4Element(name="Argon", symbol="Ar", z=18., a);
    186203
    187204  // Neon as detector gas, STP
    … …  
    232249  density *= 0.966/1.01325;
    233250
    234   G4cout<<"density of Ne857CO295N2T292 = "<<density*cm3/mg<<"  mg/cm3"<<G4endl;
     251  // G4cout<<"density of Ne857CO295N2T292 = "<<density*cm3/mg<<"  mg/cm3"<<G4endl;
    235252
    236253  G4Material* Ne857CO295N2T292 = new G4Material(name="Ne857CO295N2T292"  , density,
    … …  
    241258
    242259
     260 
     261  // Ar as detector gas,STP
     262
     263  density = 1.7836*mg/cm3 ;       // STP
     264  G4Material* Argon = new G4Material(name="Argon"  , density, ncomponents=1);
     265  Argon->AddElement(elAr, 1);
     266  /*
     267  // iso-Butane (methylpropane), STP
     268
     269  density = 2.67*mg/cm3 ;
     270  G4Material* isobutane = new G4Material(name="isoC4H10",density,nel=2) ;
     271  isobutane->AddElement(elC,4) ;
     272  isobutane->AddElement(elH,10) ;
     273
     274  // CF4 from ATLAS TRT estimation
     275
     276  G4double TRT_CF4_density = 3.9*mg/cm3;
     277  G4Material* TRT_CF4 = new G4Material(name="TRT_CF4", TRT_CF4_density, nel=2,
     278                                           kStateGas,293.15*kelvin,1.*atmosphere);
     279  */
     280
     281  // Philippe Gros T2K mixture version
     282  // Argon                                                                                                   
     283  /*                                         
     284  density     = 1.66*mg/cm3;
     285  pressure    = 1*atmosphere;
     286  temperature = 288.15*kelvin;
     287  G4Material* Argon = new G4Material(name="Ar", // z=18., a=39.948*g/mole,
     288                                     density, ncomponents=1); // kStateGas,temperature,pressure);
     289  Argon->AddElement(elAr, 1);
     290  */
     291   
     292  // IsoButane   
     293                                                                                                   
     294  density     = 2.51*mg/cm3;
     295  G4Material* Isobu = new G4Material(name="isoC4H10", z=34.,a=58.123*g/mole,
     296                                       density, kStateGas,temperature,pressure);
     297
     298  // Tetrafluoromethane   
     299                                                                                           
     300  density = 3.72*mg/cm3;
     301  G4Material* FlMet = new
     302  G4Material(name="CF4",z=42.,a=88.01*g/mole,density,kStateGas,temperature,pressure);
     303
     304  // Argon + 3% tetrafluoromethane  + 2% iso-butane     
     305                                                           
     306  density = 1.748*mg/cm3;
     307  G4Material* t2kGasMixture = new G4Material(name="t2kGasMixture", density, ncomponents=3);
     308  t2kGasMixture->AddMaterial(Argon, fractionmass = 90.9*perCent);
     309  t2kGasMixture->AddMaterial(FlMet, fractionmass = 6.3*perCent);
     310  t2kGasMixture->AddMaterial(Isobu, fractionmass = 2.8*perCent);
     311
     312
    243313  G4int i, j, jMax, k, numOfMaterials, iSan, nbOfElements, sanIndex, row;
    244314  G4double maxEnergyTransfer, kineticEnergy;
    … …  
    261331
    262332  // G4String testName = "N2";
    263   G4String testName = "Ne10CO2";
     333  // G4String testName = "Ne10CO2";
     334  G4String testName = "t2kGasMixture";
     335  // G4String testName = "Ar";
     336  //  G4String testName = "Argon";
    264337  // G4String testName = "Ne10CO2T293";
    265338  // G4String testName = "Ne857CO295N2T292";
    … …  
    372445
    373446
    374      kineticEnergy = 10.0*keV; // 100.*GeV;    // 10.0*keV;  // 110*MeV; // for proton
     447     // kineticEnergy = 10.0*keV; // 100.*GeV;    // 10.0*keV;  // 110*MeV; // for proton
     448
     449     // kineticEnergy = 5*GeV; // for electrons
     450
     451     kineticEnergy = 5*GeV*proton_mass_c2/electron_mass_c2;
     452
     453     // kineticEnergy = 5*GeV;
    375454
    376455     //     for(j=1;j<testPAIproton.GetNumberOfGammas();j++)
    377456
    378      jMax = 70; // 70;
    379 
    380      outFile<<jMax<<G4endl;
     457     // jMax = 70; // 70;
     458     jMax = 1; // 70;
     459
     460     // outFile<<jMax<<G4endl;
    381461
    382462     for( j = 0; j < jMax; j++ )
    383463     {
    384464       tau      = kineticEnergy/proton_mass_c2;
     465       // tau      = kineticEnergy/electron_mass_c2;
    385466       gamma    = tau +1.0;
    386467       bg2      = tau*(tau + 2.0);
    387468       bg = std::sqrt(bg2);
    388469       beta2    = bg2/(gamma*gamma);
    389        // G4cout<<"bg = "<<bg<<";  b2 = "<<beta2<<G4endl<<G4endl;
     470       G4cout<<"bg = "<<bg<<";  b2 = "<<beta2<<G4endl<<G4endl;
    390471       rateMass = electron_mass_c2/proton_mass_c2;
    391 
    392        Tmax     = 2.0*electron_mass_c2*bg2
    393                    /(1.0+2.0*gamma*rateMass+rateMass*rateMass);
     472       
     473       Tmax = 2.0*electron_mass_c2*bg2/(1.0+2.0*gamma*rateMass+rateMass*rateMass);
     474       // Tmax = 0.5*kineticEnergy;
    394475
    395476       Tkin = maxEnergyTransfer;
    … …  
    404485       }
    405486       G4PAIxSection testPAIproton(k,Tkin,bg2);
    406        
     487       /*       
    407488       G4cout 
    408489         //      << kineticEnergy/keV<<"\t\t"
    … …  
    444525       //             <<testPAIproton.GetPAItable(1,j)*cm<<"\t\t"<<G4endl;
    445526
    446             
    447        /*
     527       */           
     528       
    448529       outFile<<testPAIproton.GetSplineSize()-1<<G4endl;
    449530
    … …  
    451532       {
    452533       outFile
     534               << testPAIproton.GetSplineEnergy(i)/keV       << "\t"
     535         //  << testPAIproton.GetIntegralCerenkov(i)*cm    << "\t"
     536         //     << testPAIproton.GetIntegralMM(i)*cm    << "\t"
     537         //     << testPAIproton.GetIntegralPlasmon(i)*cm     << "\t"
     538         //      << testPAIproton.GetIntegralResonance(i)*cm   << "\t"
     539               << testPAIproton.GetIntegralPAIxSection(i)*cm << "\t"
     540               << G4endl;
     541
     542       G4cout
    453543               << testPAIproton.GetSplineEnergy(i)/keV       << "\t"
    454544               << testPAIproton.GetIntegralCerenkov(i)*cm    << "\t"
    … …  
    458548               << testPAIproton.GetIntegralPAIxSection(i)*cm << "\t"
    459549               << G4endl;
     550
    460551       }
    461552       
    462        */
     553     
    463554
    464555       
Note: See TracChangeset for help on using the changeset viewer.