/[PAMELA software]/DarthVader/CalorimeterLevel2/src/CaloLevel1.cpp
ViewVC logotype

Diff of /DarthVader/CalorimeterLevel2/src/CaloLevel1.cpp

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

revision 1.6 by mocchiut, Fri Mar 30 11:17:16 2007 UTC revision 1.9 by mocchiut, Fri Nov 9 10:38:25 2007 UTC
# Line 8  Line 8 
8  ClassImp(CaloStrip);  ClassImp(CaloStrip);
9  ClassImp(CaloLevel1);  ClassImp(CaloLevel1);
10    
11    Float_t CaloStrip::UXal = CTX;
12    Float_t CaloStrip::UYal = CTY;
13    Float_t CaloStrip::UZal = CTZ;
14    Bool_t CaloStrip::paramload  = false;
15    
16  /**  /**
17   * CaloStrip default constructor   * CaloStrip default constructor
18  **/  **/
# Line 24  CaloStrip::CaloStrip(Bool_t mechalig) { Line 29  CaloStrip::CaloStrip(Bool_t mechalig) {
29    this->Clear();    this->Clear();
30    if ( mechalig ){    if ( mechalig ){
31      ismech = true;      ismech = true;
32        paramload = true;
33      UXal = MECHCTX;      UXal = MECHCTX;
34      UYal = MECHCTY;      UYal = MECHCTY;
35      UZal = MECHCTZ;      UZal = MECHCTZ;
# Line 51  CaloStrip::CaloStrip(CaloLevel1 *calo, B Line 57  CaloStrip::CaloStrip(CaloLevel1 *calo, B
57    this->Clear();    this->Clear();
58    if ( mechalig ){    if ( mechalig ){
59      ismech = true;      ismech = true;
60        paramload = true;
61      UXal = MECHCTX;      UXal = MECHCTX;
62      UYal = MECHCTY;      UYal = MECHCTY;
63      UZal = MECHCTZ;      UZal = MECHCTZ;
# Line 78  void CaloStrip::Clear() {   Line 85  void CaloStrip::Clear() {  
85   **/   **/
86  void CaloStrip::UseStandardAlig(){  void CaloStrip::UseStandardAlig(){
87    //    //
88    ismech = false;    if ( !paramload ){
   //  
   stringstream aligfile;  
   Int_t error = 0;  
   FILE *f = 0;  
   ifstream badfile;  
   GL_PARAM *glparam = new GL_PARAM();  
   //  
   TString host = "mysql://localhost/pamelaprod";  
   TString user = "anonymous";  
   TString psw = "";  
   const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");  
   const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");  
   const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");  
   if ( !pamdbhost ) pamdbhost = "";  
   if ( !pamdbuser ) pamdbuser = "";  
   if ( !pamdbpsw ) pamdbpsw = "";  
   if ( strcmp(pamdbhost,"") ) host = pamdbhost;  
   if ( strcmp(pamdbuser,"") ) user = pamdbuser;  
   if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;  
   TSQLServer *dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());  
   //  
   UXal = 0.;  
   UYal = 0.;  
   UZal = 0.;  
   //  
   if ( dbc ){  
     //  
     // determine where I can find calorimeter ADC to MIP conversion file    
     //  
     printf(" Querying DB for calorimeter parameters files...\n");  
89      //      //
90        paramload = true;
91        ismech = false;
92      //      //
93        stringstream aligfile;
94        Int_t error = 0;
95        FILE *f = 0;
96        ifstream badfile;
97        GL_PARAM *glparam = new GL_PARAM();
98        //
99        TString host = "mysql://localhost/pamelaprod";
100        TString user = "anonymous";
101        TString psw = "";
102        const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
103        const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
104        const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
105        if ( !pamdbhost ) pamdbhost = "";
106        if ( !pamdbuser ) pamdbuser = "";
107        if ( !pamdbpsw ) pamdbpsw = "";
108        if ( strcmp(pamdbhost,"") ) host = pamdbhost;
109        if ( strcmp(pamdbuser,"") ) user = pamdbuser;
110        if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
111        TSQLServer *dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
112        //
113        UXal = 0.;
114        UYal = 0.;
115        UZal = 0.;
116      //      //
117      error = 0;      if ( dbc ){
118      error = glparam->Query_GL_PARAM(1000,102,dbc);        //
119      if ( error >= 0 ){        // determine where I can find calorimeter ADC to MIP conversion file  
120        //        //
121        aligfile.str("");        printf(" Querying DB for calorimeter parameters files...\n");
       aligfile << glparam->PATH.Data() << "/";  
       aligfile << glparam->NAME.Data();  
122        //        //
123        printf("\n Using parameter file: \n %s \n",aligfile.str().c_str());        //
124        f = fopen(aligfile.str().c_str(),"rb");        //
125        if ( f ){        error = 0;
126          error = glparam->Query_GL_PARAM(1000,102,dbc);
127          if ( error >= 0 ){
128          //          //
129          fread(&UXal,sizeof(UXal),1,f);          aligfile.str("");
130          fread(&UYal,sizeof(UYal),1,f);          aligfile << glparam->PATH.Data() << "/";
131          fread(&UZal,sizeof(UZal),1,f);          aligfile << glparam->NAME.Data();
132            //
133            printf("\n Using parameter file: \n %s \n",aligfile.str().c_str());
134            f = fopen(aligfile.str().c_str(),"rb");
135            if ( f ){
136              //
137              fread(&UXal,sizeof(UXal),1,f);
138              fread(&UYal,sizeof(UYal),1,f);
139              fread(&UZal,sizeof(UZal),1,f);
140              //
141              fclose(f);
142            };
143          //          //
         fclose(f);  
144        };        };
145          dbc->Close();
146        };
147        if ( !UXal ){
148        //        //
149          printf(" No able to query DB for calorimeter parameters files\n Using hard-coded parameters \n");
150          UXal = CTX;
151          UYal = CTY;
152          UZal = CTZ;
153      };      };
154      dbc->Close();      //  
   };  
   if ( !UXal ){  
     //  
     printf(" No able to query DB for calorimeter parameters files\n Using hard-coded parameters \n");  
     UXal = CTX;  
     UYal = CTY;  
     UZal = CTZ;  
155    };    };
156    //    //
157  };  };
158    
   
159  /**  /**
160   * Given a strip returns its position in the PAMELA reference system   * Given a strip returns its position in the PAMELA reference system
161  **/  **/
# Line 176  void CaloStrip::Set(Int_t view, Int_t pl Line 187  void CaloStrip::Set(Int_t view, Int_t pl
187    fView = view + 1;    fView = view + 1;
188    fPlane = plane + 1;    fPlane = plane + 1;
189    fStrip = strip + 1;    fStrip = strip + 1;
190      //
191    if ( fPlane%2 ){    if ( fPlane%2 ){
192      lShift = +0.5;      lShift = +0.5;
193    } else {    } else {
# Line 352  Float_t CaloLevel1::GetEstrip(Int_t svie Line 364  Float_t CaloLevel1::GetEstrip(Int_t svie
364   * Given estrip entry returns energy plus view, plane and strip numbers   * Given estrip entry returns energy plus view, plane and strip numbers
365  **/  **/
366  Float_t CaloLevel1::DecodeEstrip(Int_t entry, Int_t &view, Int_t &plane, Int_t &strip){  Float_t CaloLevel1::DecodeEstrip(Int_t entry, Int_t &view, Int_t &plane, Int_t &strip){
367      Bool_t sat = false;
368      Float_t mip=this->DecodeEstrip(entry,view,plane,strip,sat);
369      return(mip);
370    };
371    
372    /**
373     * Given estrip entry returns energy plus view, plane, strip numbers and saturation info
374    **/
375    Float_t CaloLevel1::DecodeEstrip(Int_t entry, Int_t &view, Int_t &plane, Int_t &strip, Bool_t &saturated){
376    //    //
377    if ( entry>istrip ) return(0.);    if ( entry>istrip ){
378        //
379        printf(" ERROR: problems decoding entry %i, it seems that number of entries is %i\n",entry,istrip);
380        //
381        return(0.);
382      };
383    //    //
384    //  printf(" num lim %f \n",std::numeric_limits<Float_t>::max());    //  printf(" num lim %f \n",std::numeric_limits<Float_t>::max());
385    //  printf(" estrip.At(%i) = %i \n",entry,estrip.At(entry));    //  printf(" estrip.At(%i) = %i \n",entry,estrip.At(entry));
# Line 397  Float_t CaloLevel1::DecodeEstrip(Int_t e Line 423  Float_t CaloLevel1::DecodeEstrip(Int_t e
423    //    //
424    Float_t mip = ((Float_t)(eval - fbi*1000000000 -plom*10000000 -strip*100000))/tim;    Float_t mip = ((Float_t)(eval - fbi*1000000000 -plom*10000000 -strip*100000))/tim;
425    //    //
426   //  if ( view == 0 ){    saturated = false;
427  //     if ( plane%2 ){    if ( mip > 5000. ){
428  //       plane -= 1;      mip -= 5000.;
429  //     } else {      saturated = true;
430  //       plane += 1;    };
 //     };  
 //   };  
   //  
   //  
431    if ( mip > 0. && mip < 99999. ) return(mip);    if ( mip > 0. && mip < 99999. ) return(mip);
432    //    //
433    printf(" WARNING: problems decoding value %i at entry %i \n",estrip.At(entry),entry);    printf(" ERROR: problems decoding value %i at entry %i \n",estrip.At(entry),entry);
434    //    //
435    view = -1;    view = -1;
436    plane = -1;    plane = -1;
# Line 420  Float_t CaloLevel1::DecodeEstrip(Int_t e Line 442  Float_t CaloLevel1::DecodeEstrip(Int_t e
442   * Returns energy released on plane nplane (where 0<= nplane <= 43, 0 = 1Y, 1 = 1X, 2 = 2Y, 3 = 2X, etc. etc.).   * Returns energy released on plane nplane (where 0<= nplane <= 43, 0 = 1Y, 1 = 1X, 2 = 2Y, 3 = 2X, etc. etc.).
443   */   */
444  Float_t CaloLevel1::qtotpl(Int_t nplane){  Float_t CaloLevel1::qtotpl(Int_t nplane){
445      Bool_t sat = false;
446      Float_t mip = this->qtotpl(nplane,sat);
447      return(mip);
448    };
449    
450    /*
451     * Returns energy released on plane nplane (where 0<= nplane <= 43, 0 = 1Y, 1 = 1X, 2 = 2Y, 3 = 2X, etc. etc.).
452     */
453    Float_t CaloLevel1::qtotpl(Int_t nplane, Bool_t &sat){
454    //    //
455      sat = false;
456    Int_t sview = 1;    Int_t sview = 1;
457    if ( nplane%2 ) sview = 0;    if ( nplane%2 ) sview = 0;
458    //    //
459    Int_t splane = nplane-(sview+1)/2;    Int_t splane = nplane-(sview+1)/2;
460    //    //
461  //   if ( sview == 0 ){    Float_t totmip = qtotpl(sview,splane,sat);
 //     if ( splane%2 ){  
 //       splane += 1;  
 //     } else {  
 //       splane -= 1;  
 //     };  
 //   };  
   //  
   Float_t totmip = qtotpl(sview,splane);  
462    //    //
463    return(totmip);    return(totmip);
464    //    //
# Line 444  Float_t CaloLevel1::qtotpl(Int_t nplane) Line 468  Float_t CaloLevel1::qtotpl(Int_t nplane)
468   * Returns energy released on view "view" (0 = X, 1 = Y) and plane "plane" ( 0 <= plane <= 21 ).   * Returns energy released on view "view" (0 = X, 1 = Y) and plane "plane" ( 0 <= plane <= 21 ).
469   */   */
470  Float_t CaloLevel1::qtotpl(Int_t sview, Int_t splane){  Float_t CaloLevel1::qtotpl(Int_t sview, Int_t splane){
471      Bool_t sat = false;
472      Float_t mip = this->qtotpl(sview,splane,sat);
473      return(mip);
474    };
475    
476    /*
477     * Returns energy released on view "view" (0 = X, 1 = Y) and plane "plane" ( 0 <= plane <= 21 ).
478     */
479    Float_t CaloLevel1::qtotpl(Int_t sview, Int_t splane, Bool_t &sat){
480    //    //
481    Int_t view = -1;    Int_t view = -1;
482    Int_t plane = -1;    Int_t plane = -1;
483    Int_t strip = -1;    Int_t strip = -1;
484      Bool_t lsat = false;
485      sat = false;
486    //    //
487    Float_t mip = 0.;    Float_t mip = 0.;
488    Float_t totmip = 0.;    Float_t totmip = 0.;
# Line 456  Float_t CaloLevel1::qtotpl(Int_t sview, Line 491  Float_t CaloLevel1::qtotpl(Int_t sview,
491    //    //
492    for (Int_t i = 0; i<istrip; i++ ){    for (Int_t i = 0; i<istrip; i++ ){
493      //      //
494      mip = DecodeEstrip(i,view,plane,strip);      mip = DecodeEstrip(i,view,plane,strip,lsat);
495      //      //
496      if ( view == sview && splane == plane ) totmip += mip;      if ( view == sview && splane == plane ){
497          if ( lsat ) sat = true;
498          totmip += mip;
499        };
500      //      //
501      // entry are ordered by strip, plane and view number. Go out if you pass the input strip      // entry are ordered by strip, plane and view number. Go out if you pass the input strip
502      //      //
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23