Changeset 10060 for trunk/Mars/mcorsika/MCorsikaEvtHeader.h

Timestamp:

11/29/10 14:24:23 (14 years ago)

Author:

rohlfs

Message:

two new command line arguments of readcorsika: -A=arrayNo and -T=telescopeNo. New design of program flow in MCorsikaRead: It is now determined by the order of the data blocks in the input file.

File:

• trunk/Mars/mcorsika/MCorsikaEvtHeader.h (modified) (3 diffs)

trunk/Mars/mcorsika/MCorsikaEvtHeader.h

@@ -21 +21 @@
     UInt_t   fEvtNumber;              // Event number
     UInt_t   fNumReuse;               // Counter of the reuse of the same shower
+    UInt_t   fTotReuse;               //! number of reuse of the same shower     
 //    UInt_t   fParticleID;             // Particle ID (see MMcEvtBasic or CORSIKA manual)
     Float_t  fTotalEnergy;            // [GeV] 
@@ -51 +52 @@
     UInt_t GetEvtNumber() const { return fEvtNumber; }
     UInt_t GetNumReuse() const { return fNumReuse; }
+    UInt_t GetTotReuse() const { return fTotReuse; }
 //    UInt_t GetParticleID() const { return fParticleID; }
 
@@ -67 +69 @@
     Double_t GetImpact() const;
 
+    void GetArrayOffset(Int_t arrayIdx, Float_t & xArrOff, Float_t & yArrOff)
+                     {xArrOff = fTempY[arrayIdx]; yArrOff=-fTempX[arrayIdx]; }
+    void SetTelescopeOffset(Int_t arrayIdx, Float_t xTelOff, Float_t yTelOff)
+                     {fNumReuse = arrayIdx; fX = xTelOff; fY = yTelOff;}
+
     void IncNumReuse() { fNumReuse++; }
-    void ResetNumReuse() { fNumReuse=(UInt_t)-1; }
+    void ResetNumReuse() { fNumReuse=0; }
 
     void InitXY() { fX=fTempY[fNumReuse]; fY=-fTempX[fNumReuse]; }
     void AddXY(Float_t x, Float_t y) { fX+=x; fY+=y; }
 
-    Int_t ReadEvt(MCorsikaFormat *informat);    // read in event header block
+    Int_t ReadEvt(Float_t * f);                 // read in event header block
     Int_t ReadEvtEnd(MCorsikaFormat *informat); // read in event end block