/[PAMELA software]/trieste/pamVMC/run_g4_testTRK.C
ViewVC logotype

Diff of /trieste/pamVMC/run_g4_testTRK.C

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by pamelats, Wed Mar 4 12:51:04 2009 UTC revision 1.2 by pizzolot, Tue Mar 24 14:03:55 2009 UTC
# Line 1  Line 1 
1    
2  // $Id: run_g4.C,v 1.0 2007/06/01  // $Id: run_g4.C,v 1.0 2007/06/01
3  //  //
4  // Macro for running G4PAM_VMC with Geant4.  // Macro for running G4PAM_VMC with Geant4.
# Line 27  void run_g4_testTRK(const TString& confi Line 28  void run_g4_testTRK(const TString& confi
28      PamVMCApplication* appl      PamVMCApplication* appl
29        = new PamVMCApplication("PAMG4_VMC", "PAMELA GEANT4 VMC application");        = new PamVMCApplication("PAMG4_VMC", "PAMELA GEANT4 VMC application");
30            
31    
32      //appl->GetPrimaryGenerator()->SetNofPrimaries(1);        //appl->GetPrimaryGenerator()->SetNofPrimaries(1);  
33          appl->SetVerboseLevel(0);          appl->SetVerboseLevel(0);
34                
# Line 37  void run_g4_testTRK(const TString& confi Line 39  void run_g4_testTRK(const TString& confi
39      // Initialize MC      // Initialize MC
40      appl->InitMC(configMacro);      appl->InitMC(configMacro);
41    
42        // Cecilia: set seed for random generator  
43        TRandom3 *rndgen = new TRandom3(123);
44        ((TGeant4*)gMC)->SetRandom(rndgen);
45    
46    
47      // Customise Geant4 setting after initialization:      // Customise Geant4 setting after initialization:
48      // Physics list      // Physics list
49  //     gMC->SetCut("CUTGAM", 1.0e-6);  //     gMC->SetCut("CUTGAM", 1.0e-6);
# Line 53  void run_g4_testTRK(const TString& confi Line 60  void run_g4_testTRK(const TString& confi
60            
61      ((TGeant4*)gMC)->ProcessGeantMacro("g4config2.in");      ((TGeant4*)gMC)->ProcessGeantMacro("g4config2.in");
62      // Visualization settings      // Visualization settings
63       ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");      ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");
64    
65     ((TGeant4*)gMC)->SetMaxNStep(50000);     ((TGeant4*)gMC)->SetMaxNStep(50000);
66     //appl->GetPrimaryGenerator()->SetParticle(1000060120);     //appl->GetPrimaryGenerator()->SetParticle(1000060120);
67      appl->GetPrimaryGenerator()->SetParticle(kElectron);      appl->GetPrimaryGenerator()->SetParticle(kElectron);
68      //appl->GetPrimaryGenerator()->SetDirection(0.,0.);      //appl->GetPrimaryGenerator()->SetDirection(0.,0.);
69       //appl->GetPrimaryGenerator()->SetKinEnergy(10.);      appl->GetPrimaryGenerator()->SetKinEnergy(10.);
70       //appl->GetPrimaryGenerator()->SetRigidity(1.0); //GV       //appl->GetPrimaryGenerator()->SetRigidity(1.0); //GV
71      //appl->GetPrimaryGenerator()->SetPosition(1.,1.,130.);      //appl->GetPrimaryGenerator()->SetPosition(1.,1.,130.);
72    
73      // Run MC      // Run MC
74      for (Int_t i=0; i<3; i++){      for (Int_t i=0; i<2; i++){
75        cout<<"Event "<<i<<endl;        cout<<"in Loop: Event nr."<<i<<" --<<--<<--<<--<<--<<--<<--"<<endl;
76       appl->GetPrimaryGenerator()->SetKinEnergy(10.);        // ritabrata 04.03.2009
77       if (i==0) appl->GetPrimaryGenerator()->SetKinEnergy(0.01);        //appl->GetPrimaryGenerator()->SetKinEnergy(0.1);
78         appl->GetPrimaryGenerator()->GenSphDist(120.);        //if (i==0) appl->GetPrimaryGenerator()->SetKinEnergy(0.01);
79         //appl->GetPrimaryGenerator()->GenSphDist(120.,0.,90.,0.,360.);        //appl->GetPrimaryGenerator()->GenSphDist(120.);
80        // appl->GetPrimaryGenerator()->GenPosition(-22.4,22.4, -18.5, 18.5, 110.0,110.0);        // endritabrata
81        //appl->GetPrimaryGenerator()->GenDirection(0.,TMath::Pi()/2.,0.,2*TMath::Pi());  
82        //appl->GetPrimaryGenerator()->GenSpe(1.,5.,-2.75);        // cecilia
83          //appl->GetPrimaryGenerator()->GenPosition(-22.4,22.4, -18.5, 18.5, 110.0,110.0);
84          // appl->GetPrimaryGenerator()->GenPosition(-8.,8., -6., 6., 110.0,110.0);
85          appl->GetPrimaryGenerator()->GenSphPhiThe();
86          //appl->GetPrimaryGenerator()->GenDirection(0.,0.1, 0.,2*TMath::Pi()); // 5 gradi
87          // appl->GetPrimaryGenerator()->GenSpe(7.,20.,-2.75);
88        //  appl->GetPrimaryGenerator()->Print();        //  appl->GetPrimaryGenerator()->Print();
89        // if (i==9) ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");        // if (i==9) ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");
90          appl->GetPrimaryGenerator()->Print();
91        appl->RunMC(1);        appl->RunMC(1);
92    
93      }      }
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23