/[PAMELA software]/PamVMC/trk/include/PamVMCTrkF77GpsSpe.h
ViewVC logotype

Annotation of /PamVMC/trk/include/PamVMCTrkF77GpsSpe.h

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.5 - (hide annotations) (download)
Fri Jun 12 18:39:49 2009 UTC (15 years, 5 months ago) by pam-rm2
Branch: MAIN
CVS Tags: v1r0, HEAD
Changes since 1.1: +0 -0 lines
File MIME type: text/plain 
- Introduced user-defined names of output files and random seeds number.
Users can do it use options of PamVMCApplication constructor:
PamVMCApplication(const char* name,  const char *title, const char*
filename="pamtest", Int_t seed=0).
The Random object that I use is TRandom3 object which has astronomical
large period (in case of default initialization 0). All random generators
in the code use this object by calling of gRandom singleton which keeps
it.

- Corrected TOF digitization routine. No problems with TDC hits due to
hadronic interactions anymore.

- Some small changes was done to compile code under Root 5.23. +
geant4_vmc v. 2.6 without any warnings

- Some classes of PamG4RunConfiguartion was changed for geant4_vmc v.
2.6.Some obsolete classes was deleted as soon as developers implemented
regions.

- Navigation was changed from "geomRootToGeant4" to "geomRoot", because on
VMC web page written that as soon as Geant4 has no option ONLY/MANY
translation of overlapped geometry to Geant4 through VGM could be wrong.
I'd like to stay with Root navigation:
http://root.cern.ch/root/vmc/Geant4VMC.html. This should be default
option.

- New Tracker digitization routine written by Sergio was implemented

- PamVMC again became compatible with geant4_vmc v.2.5 and ROOT 5.20.
 The problem was that ROOT developers introduced in TVirtualMC class a new
method SetMagField and new base class:TVirtualMagField from which
user-defined classes shoukd be derived
1 nikolas 1.1 #ifndef PAMVMCTRK_F77_GPSSPE
2     #define PAMVMCTRK_F77_GPSSPE
3    
4     #include <TObject.h>
5     #include <TClonesArray.h>
6     #include "PamRootManager.h"
7     #include <iostream>
8    
9     using std::cout;
10     using std::endl;
11     // structures for C++ <=> COMMON data exchange
12    
13     #define maxstr 1000
14    
15     struct cGPSSPE {
16     Int_t nstrpx;
17     Int_t nstrpy;
18     Int_t istripx[maxstr];
19     Int_t istripy[maxstr];
20     Float_t qstripx[maxstr];
21     Float_t qstripy[maxstr];
22     Float_t xstripx[maxstr];
23     Float_t ystripy[maxstr];
24     Int_t npstripx[maxstr];
25     Int_t npstripy[maxstr];
26     Int_t ntstripx[maxstr];
27     Int_t ntstripy[maxstr];
28     } ;
29    
30     extern "C" {
31     extern struct cGPSSPE gpsspe_;
32     }
33    
34     class GPSSPEData: public TObject {
35    
36     public:
37    
38     Int_t fistrip;
39     Float_t fqstrip;
40     Float_t fstripc;
41     Int_t fnpstrip;
42     Int_t fntstrip;
43    
44     GPSSPEData(){fistrip=fnpstrip=fntstrip=0; fqstrip=fstripc=0.;}
45    
46     virtual ~GPSSPEData(){; };
47    
48     void Clean(){fistrip=fnpstrip=fntstrip=0; fqstrip=fstripc=0.;}
49    
50     void Clear(Option_t * = "") { this->~GPSSPEData(); }
51    
52     void Print(Option_t * ="" ) const {
53     cout<<" ISTRIP: " <<fistrip<<" (Strip ID in all plane)"<<endl
54     <<" QSTRIP: " <<fqstrip<<" (ADC value in strip)"<<endl
55     <<" CSTRIP: " <<fstripc<<" (Strip coordinate)"<<endl
56     <<" NPSTRIP: " <<fnpstrip<<" (Plane No)"<<endl
57     <<" NTSTRIP: " <<fntstrip<<" (Sensor No)"<<endl;
58     }
59    
60     ClassDef(GPSSPEData,1);
61     };
62    
63     class pGPSSPEHits{
64    
65     private:
66    
67     Int_t fnstrpx;
68     Int_t fnstrpy;
69     TClonesArray * fGPSSPEDataXColl;
70     TClonesArray * fGPSSPEDataYColl;
71    
72     public:
73    
74     pGPSSPEHits():fnstrpx(0), fnstrpy(0){
75     fGPSSPEDataXColl = new TClonesArray("GPSSPEData",1000);
76     fGPSSPEDataYColl = new TClonesArray("GPSSPEData",1000);
77     }
78    
79     virtual ~pGPSSPEHits(){ delete fGPSSPEDataYColl; delete fGPSSPEDataXColl; }
80    
81     void Clean() {
82     fnstrpx=fnstrpy=0;
83     (*fGPSSPEDataXColl).Clear("C");
84     (*fGPSSPEDataYColl).Clear("C");
85     }
86    
87     void Clear( Option_t * = "") { Clean(); }
88    
89     void Compress(){
90     fGPSSPEDataXColl->Compress();
91     fGPSSPEDataYColl->Compress();
92     }
93    
94     GPSSPEData * CreateXHit(Int_t flag){
95     if(!fGPSSPEDataXColl->At(flag)){
96     return (GPSSPEData *)fGPSSPEDataXColl->New(flag);
97     }
98     return (GPSSPEData *)fGPSSPEDataXColl->At(flag);
99     }
100    
101     GPSSPEData * CreateYHit(Int_t flag){
102     if(!fGPSSPEDataYColl->At(flag)){
103     return (GPSSPEData *)fGPSSPEDataYColl->New(flag);
104     }
105     return (GPSSPEData *)fGPSSPEDataYColl->At(flag);
106     }
107    
108     Int_t GetNXHit(){ return fnstrpx; }
109    
110     Int_t GetNYHit(){ return fnstrpy; }
111    
112     GPSSPEData * GetXHit(Int_t flag){
113     if(fGPSSPEDataXColl->At(flag)){
114     GPSSPEData * hit = (GPSSPEData *)fGPSSPEDataXColl->At(flag);
115     return hit;
116     }
117     cout<<"Tracker GPSSPE_X Hit at "<<flag<<" is empry"<<endl;
118     return 0;
119     }
120    
121     GPSSPEData * GetYHit(Int_t flag){
122     if(fGPSSPEDataYColl->At(flag)){
123     GPSSPEData * hit = (GPSSPEData *)fGPSSPEDataYColl->At(flag);
124     return hit;
125     }
126     cout<<"Tracker GPSSPE_Y Hit at "<<flag<<" is empry"<<endl;
127     return 0;
128     }
129    
130     void SetData(cGPSSPE * c){
131    
132     fnstrpx=c->nstrpx;
133     fnstrpy=c->nstrpy;
134     for (Int_t i=0; i<fnstrpx; i++){
135     GPSSPEData *hx = CreateXHit(i);
136     hx->fqstrip=c->qstripx[i];
137     hx->fstripc=c->xstripx[i];
138     hx->fnpstrip=c->npstripx[i];
139     hx->fntstrip=c->ntstripx[i];
140     hx->fistrip=c->istripx[i];
141     }
142     for (Int_t i=0; i<fnstrpy; i++){
143     GPSSPEData *hy = CreateYHit(i);
144     hy->fqstrip=c->qstripy[i];
145     hy->fstripc=c->ystripy[i];
146     hy->fnpstrip=c->npstripy[i];
147     hy->fntstrip=c->ntstripy[i];
148     hy->fistrip=c->istripy[i];
149     }
150    
151     }
152    
153     void Print(Option_t * ="" ) const
154     { cout << "pGPSSPEHitData: "<<endl
155     <<"NSTRPX Number of hits (X View):" << fnstrpx
156     <<" NSTRPY Number of hits (Y View):" << fnstrpy << endl;
157     for (Int_t i=0; i<fnstrpx; i++){
158     cout<<"HIT #" << i << " (X View):" << endl;
159     GPSSPEData * hx = (GPSSPEData *)fGPSSPEDataXColl->At(i);
160     if(hx) hx->Print();
161     }
162     for (Int_t i=0; i<fnstrpy; i++){
163     cout<<"HIT #" << i << " (Y View):" << endl;
164     GPSSPEData * hy = (GPSSPEData *)fGPSSPEDataYColl->At(i);
165     if (hy) hy->Print();
166     }
167     }
168    
169     ClassDef(pGPSSPEHits,1)
170    
171     };
172    
173     #endif
174    
  ViewVC Help
Powered by ViewVC 1.1.23