/[PAMELA software]/PamelaDigitizer/Pamelagp2Digits.cxx
ViewVC logotype

Diff of /PamelaDigitizer/Pamelagp2Digits.cxx

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

revision 1.1 by silvio, Thu Sep 13 11:00:53 2007 UTC revision 1.10 by pamelats, Fri Jun 6 15:16:08 2008 UTC
# Line 4  Line 4 
4  //      The conversion is done using a slightly patched version of  //      The conversion is done using a slightly patched version of
5  //      Rene Brun's h2root.  //      Rene Brun's h2root.
6  //  //
7  //      Created : Jan Conrad (conrad@particle.kth.se)  //      Created : 2006 Jan Conrad  (conrad@particle.kth.se)
8  //      Modified: Silvio Orsi (silvio.orsi@roma2.infn.it)  //      Modified: 2007 Silvio Orsi (silvio.orsi@roma2.infn.it)
9  //  //
10  //  //
11  //  HOWTO (RECOMMENDED):  //  HOWTO (RECOMMENDED):
12  //  ./ Pamelagp2Digits filename.his  //  ./ Pamelagp2Digits filename.his
13  //  It generates filename.pam (PAMELA raw data) and filename.gp.root (sim info)  //  It generates filename.pam (PAMELA raw data) and filename.gp.root (sim info)
14  //    //  
15  //  HOWTO (EXPERT, only for testing purposes)  //  NB: more details are present in the file HOW-TO-DIGIT.TXT
16    //
17    //  HOWTO (only for testing purposes)
18  //  Pamelagp2Digits hbook_filename root_filename  raw_filename compress tolower  //  Pamelagp2Digits hbook_filename root_filename  raw_filename compress tolower
19  //  - If the 2nd and 3rd parameters are missing the names are generated from the hbook file name.  //  - If the 2nd and 3rd parameters are missing the names are generated from the hbook file name.
20  //  - The following flags are from h2root and default values are used:  //  - The following flags are from h2root and default values are used:
# Line 24  Line 26 
26  //                convertex to lower case  //                convertex to lower case
27  //  //
28  //  MODIFICATIONS:  //  MODIFICATIONS:
29  //  v. 1.0 (S.Orsi, September 2007)  //  v. 1.5 (S.Orsi, Roma2, 11 October 2007)
30    //  - file HOW-TO-DIGIT.TXT written (mainly by S.Bottai): contains information on how to run the chain from the DB creation to DarthVader through Digitizer, yoda and YodaProfiler
31    //
32    //  v. 1.4 (S.Orsi, Roma2, 11 October 2007)
33    //  - similar to v. 1.3
34    //
35    //  v. 1.3 (S.Orsi, Roma2, 11 October 2007)
36    //  - TOF: + changes and major corrections in ADC and TDC (W.Menn, D.Campana, S.Orsi)
37    //         + flag DEBUG added (default: DEBUG=false; if true, print on screem some messages to debug)
38    //  - AC: + shift register filled;
39    //        + minor bugs fixed
40    //        + counters increment by 1 at each event (they were constant before)
41    //  - S4: + minor bugs fixed (S.Borisov)
42    //        + the routine DigitizeS4() is now inserted in the code. For DV versions earlier than October 2007, the option "-S4" is still needed because S4 calibration is not present
43    //  - corrected severe bug that caused yoda to process only part of large data files (file.pam). S4 data were written to file even if the DigitizeS4() routine was not called. Previous versions of the digitizer should have crashed always, but unexpectedly worked on small files.
44    //  - file HOW-TO-DIGIT.TXT written (mainly by S.Bottai): contains information on how to run the chain from the DB creation to DarthVader through Digitizer, yoda and YodaProfiler. The file is included only from release 1.5
45    //  - digitizer version numbers modified according to CVS repository (v1.00->v1.1, v1.01->v1.2)
46    //  - some flags (LDFLAGS) in the Makefile have been changed; they are needed on some machines (e.g. pamelatov.roma2.infn.it), but give errors on others (e.g. pamelaws02.na.infn.it). See Makefile for details if the compilation fails due to library errors)
47    //
48    //  v. 1.2 (S.Orsi, Roma2, 28 September 2007)
49    //  - TOF TDC inserted (to be tuned)
50    //  - improved tracker dEdx, tracker saturation inserted (S.Bottai)
51    //  - S4 routine inserted in code (not in read-out)
52    //  - bugs with packet length corrected
53    //  - small bugs fixed
54    //  - NB: Run DarthVader with the options: "DarthVader -zerofill -idRun 1 -ALL +RUN +CAL [ --no-crosstalk ] +TRK +TOF +AC +TRG +ND" or  "DarthVader -zerofill -idRun 1 +CAL [ --no-crosstalk ] -S4", i.e. without S4
55    //
56    //  v. 1.1 (S.Orsi, Roma2, 13 September 2007)
57  //  - trigger (dummy), tof (preliminary, adc only), ND (preliminary) inserted  //  - trigger (dummy), tof (preliminary, adc only), ND (preliminary) inserted
58  //  - various changes and bug fixes  //  - various changes and bug fixes
59  //  - padding increased (fPADbuffer=64, fPadding = padbytes+64)  //  - padding increased (fPADbuffer=64, fPadding = padbytes+64)
60  //  //
61  //  v. beta (J.Conrad)  //  v. beta (J.Conrad, KTH, 2006)
62  //  - compiles; includes pscu, calo, trk, ac  //  - compiles; includes pscu, calo, trk, ac; not present on cvs
63  //  //
64  /////////////////////////////////////////////////////////////////////////  /////////////////////////////////////////////////////////////////////////
65    
# Line 49  Line 78 
78  #include "TMath.h"  #include "TMath.h"
79  #include "Digitizer.h"  #include "Digitizer.h"
80    
81    using namespace std;
82    
83  int Error;   //to be removed soon  int Error;   //to be removed soon
84    
85  //  Define the names of the Fortran common blocks for the different OSs  //  Define the names of the Fortran common blocks for the different OSs
# Line 246  Int_t golower  = 1; Line 277  Int_t golower  = 1;
277  Int_t bufsize  = 64000;  Int_t bufsize  = 64000;
278  Int_t optcwn = 1;  Int_t optcwn = 1;
279    
280    char *DigitizerInfo(bool print);
281    
282  int main(int argc, char **argv)  int main(int argc, char **argv)
283    
284  {  {
285      char* cha=DigitizerInfo(1);
286    
287    if (argc < 2) {    if (argc < 2) {
288       printf("Error: Pamelagp2Digits \n");       printf("Error: Pamelagp2Digits \n");
289       printf("Pamelagp2Digits file.hbook file.root file.pam [compress] [tolower] [lrecl] [bufsize] [optcwn] \n");       printf("Pamelagp2Digits file.hbook file.root file.pam [compress] [tolower] [lrecl] [bufsize] [optcwn] \n");
# Line 786  void convert_directory(const char *dir, Line 820  void convert_directory(const char *dir,
820    
821    tree->Write();    tree->Write();
822    
823    std::cout << "Invoking Digitizer" << endl << flush;    cout << "Invoking Digitizer" << endl << flush;
824      char *ndl[8],*q,*q1;
825      q=(char*)malloc(5 *sizeof(char));
826      q1=(char*)malloc(5 *sizeof(char));
827      ndl[0]="nspe";
828      ndl[1]="ntof";
829      ndl[2]="ncat";
830      ndl[3]="ncas";
831      ndl[4]="ncar";
832      ndl[5]="ncal";
833      ndl[6]="nndd";
834      ndl[7]="nstr";
835    
836      int nspe=200,ntof=200,ncat=50,ncas=50,ncar=100,ncal=1000,nnd=200,nstr=1000,par[8],tmp=0,t;
837      par[0]=nspe;
838      par[1]=ntof;
839      par[2]=ncat;
840      par[3]=ncas;
841      par[4]=ncar;
842      par[5]=ncal;
843      par[6]=nnd;
844      par[7]=nstr;
845    
846      ifstream np("vectpar.dat",ios::in);
847      if(!np)printf("ATTENTION: Using default vector legths!\n");
848      while(!np.eof()){
849        np>>q>>tmp;
850        if(np.eof())break;
851        if(tmp<=0){
852          cout<<"ATTENTION: Length of one vector is negative or equal 0!"<<endl<<q<<"="<<tmp<<endl;
853          break;
854        }
855        t=0;
856        while(tmp){
857          q1=ndl[t];
858          if(q[0]==q1[0] && q[1]==q1[1] && q[2]==q1[2] && q[3]==q1[3]){
859            par[t]=tmp;
860            tmp=0;
861            cout<<q<<" "<<par[t]<<"["<<t<<"]"<<endl;
862          }
863          t++;
864        }
865      }
866      nspe=par[0];
867      ntof=par[1];
868      ncat=par[2];
869      ncas=par[3];
870      ncar=par[4];
871      ncal=par[5];
872      nnd=par[6];
873      nstr=par[7];
874    
875    Digitizer* dig = new Digitizer(tree,file_raw);    Digitizer* dig = new Digitizer(tree,file_raw,nspe,ntof,ncat,ncas,ncar,ncal,nnd,nstr);
876    dig->Loop();    dig->Loop();
877    dig->Close();    dig->Close();
878    
879    std::cout << "Finished" << endl << flush;    cout << "Finished" << endl << flush;
880    
881  }  }
882    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.23