/[PAMELA software]/DarthVader/TrackerLevel2/src/TrkLevel1.cpp
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/TrkLevel1.cpp

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

revision 1.13 by pam-fi, Mon Feb 5 16:01:51 2007 UTC revision 1.16 by pam-fi, Fri Apr 27 10:39:57 2007 UTC
# Line 313  void TrkCluster::Dump(){ Line 313  void TrkCluster::Dump(){
313  /**  /**
314   * Method to fill a level1 struct with only one cluster (done to use F77 p.f.a. routines on a cluster basis).   * Method to fill a level1 struct with only one cluster (done to use F77 p.f.a. routines on a cluster basis).
315   */   */
316  cTrkLevel1* TrkCluster::GetLevel1Struct(){  void TrkCluster::GetLevel1Struct(cTrkLevel1* l1){
317                                    
318  //    cTrkLevel1* l1 = new cTrkLevel1;  //    cTrkLevel1* l1 = new cTrkLevel1;
319    
320      cTrkLevel1* l1 = &level1event_ ;  //    cTrkLevel1* l1 = &level1event_ ;
321                    
322      l1->nclstr1 = 1;      l1->nclstr1 = 1;
323      l1->view[0] = view;      l1->view[0] = view;
# Line 335  cTrkLevel1* TrkCluster::GetLevel1Struct( Line 335  cTrkLevel1* TrkCluster::GetLevel1Struct(
335          l1->clbad[i] = clbad[i];          l1->clbad[i] = clbad[i];
336      };      };
337            
338      return l1;  //    return l1;
339  };  };
340  //--------------------------------------  //--------------------------------------
341  //  //
# Line 397  Float_t TrkCluster::GetETA(Int_t neta, f Line 397  Float_t TrkCluster::GetETA(Int_t neta, f
397  //    cout << "GetETA(neta,angle) "<< neta << " "<< angle;  //    cout << "GetETA(neta,angle) "<< neta << " "<< angle;
398  //      LoadPfaParam();  //      LoadPfaParam();
399    
400        TrkParams::Load(4);
401        if( !TrkParams::IsLoaded(4) ){
402            cout << "int Trajectory::DoTrack2(float* al) --- ERROR --- p.f.a. parameters  not loaded"<<endl;
403            return 0;
404        }
405    
406      float ax = angle;      float ax = angle;
407      int ic = 1;      int ic = 1;
408      GetLevel1Struct();      GetLevel1Struct();
# Line 466  void TrkLevel1::Dump(){ Line 472  void TrkLevel1::Dump(){
472   */   */
473  void TrkLevel1::SetFromLevel1Struct(cTrkLevel1 *l1, Bool_t full){  void TrkLevel1::SetFromLevel1Struct(cTrkLevel1 *l1, Bool_t full){
474    
475    //    cout << "void TrkLevel1::SetFromLevel1Struct(cTrkLevel1 *l1, Bool_t full)"<<endl;
476        
477        Clear();
478      //  ---------------      //  ---------------
479      //  *** CLUSTER ***      //  *** CLUSTER ***
480      //  ---------------      //  ---------------
# Line 522  void TrkLevel1::SetFromLevel1Struct(cTrk Line 531  void TrkLevel1::SetFromLevel1Struct(cTrk
531   * Fills a struct cTrkLevel1 with values from a TrkLevel1 object (to put data into a F77 common).   * Fills a struct cTrkLevel1 with values from a TrkLevel1 object (to put data into a F77 common).
532   */   */
533    
534  cTrkLevel1* TrkLevel1::GetLevel1Struct() {  void TrkLevel1::GetLevel1Struct(cTrkLevel1* l1) {
535            
536      cTrkLevel1 *l1=0;  //    cTrkLevel1* l1 = &level1event_ ;
537      //      
538      for(Int_t i=0; i<12 ; i++){      for(Int_t i=0; i<12 ; i++){
539          l1->good[i] = good[i];          l1->good[i] = good[i];
540          for(Int_t j=0; j<24 ; j++){          for(Int_t j=0; j<24 ; j++){
541              l1->cnev[j][i]    = cn[j][i];              l1->cnev[j][i]    = cn[j][i]  ;
542  //          l1->cnrmsev[j][i] = cnrms[j][i];              l1->cnnev[j][i]   = cnn[j][i] ;
543              l1->cnnev[j][i]   = cnn[j][i];              l1->cnrmsev[j][i] = 0. ;
544          };          };
545            l1->fshower[i] = 0;
546      };      };
547        
548  //  *** CLUSTERS ***      l1->nclstr1=0;
549        l1->totCLlength=0;
550        Int_t index=0;
551      if(Cluster){      if(Cluster){
552          l1->nclstr1 =  Cluster->GetEntries();          Int_t i=0;
553          for(Int_t i=0;i<l1->nclstr1;i++){          for(Int_t ii=0;ii<Cluster->GetEntries();ii++){
554                        TrkCluster *clu = GetCluster(ii);
555              l1->view[i]     = ((TrkCluster *)Cluster->At(i))->view;              // ----------------------------------------
556              l1->maxs[i]     = ((TrkCluster *)Cluster->At(i))->maxs;              // attenzione!!
557              // COMPLETARE //              // se il cluster non e` salvato (view = 0)
558              // COMPLETARE //              // DEVE essere escluso dal common F77
559              // COMPLETARE //              // ----------------------------------------
560              // COMPLETARE //              if(clu->view != 0 ){
561              // COMPLETARE //                  l1->view[i]     = clu->view;
562              // COMPLETARE //                  l1->ladder[i]   = clu->GetLadder();
563                            l1->maxs[i]     = clu->maxs;
564                    l1->mult[i]     = clu->GetMultiplicity();
565                    l1->dedx[i]     = clu->GetSignal();
566                    l1->indstart[i] = index+1;
567                    l1->indmax[i]   = l1->indstart[i] + clu->indmax;
568                    l1->totCLlength += clu->CLlength;
569                    for(Int_t iw=0; iw < clu->CLlength; iw++){
570                        l1->clsignal[index] = clu->clsignal[iw];
571                        l1->clsigma[index]  = clu->clsigma[iw];
572                        l1->cladc[index]    = clu->cladc[iw];
573                        l1->clbad[index]    = clu->clbad[iw];
574                        index++;
575                    }
576                    i++;
577                }
578          }          }
579          // COMPLETARE //          l1->nclstr1 =  i;      
         // COMPLETARE //  
         // COMPLETARE //  
         // COMPLETARE //  
         // COMPLETARE //  
         // COMPLETARE //  
580      }      }
581      return l1;  
582    //    return l1;
583  }  }
584  //--------------------------------------  //--------------------------------------
585  //  //
# Line 608  TrkCluster *TrkLevel1::GetCluster(int is Line 630  TrkCluster *TrkLevel1::GetCluster(int is
630  //  //
631  //  //
632  //--------------------------------------  //--------------------------------------
633  /**  // /**
634   * Load Position-Finding-Algorythm parameters (call the F77 routine).  //  * Load Position-Finding-Algorythm parameters (call the F77 routine).
635   *  //  *
636   */  //  */
637  int TrkLevel1::LoadPfaParam(TString path){  // int TrkLevel1::LoadPfaParam(TString path){
638                    
639      if( strcmp(path_.path,path.Data()) ){  //     if( path.IsNull() ){
640          cout <<"Loading p.f.a. param::eters\n";  //      path = gSystem->Getenv("PAM_CALIB");
641          strcpy(path_.path,path.Data());  //      if(path.IsNull()){
642          path_.pathlen = path.Length();  //          cout << " TrkLevel1::LoadPfaParam() ==> No PAMELA environment variables defined "<<endl;
643          path_.error   = 0;  //          return 0;
644          return readetaparam_();  //      }
645      }    //      path.Append("/trk-param/eta_param-0/");
646      return 0;  //     }
647  }  
648    //     strcpy(path_.path,path.Data());
649    //     path_.pathlen = path.Length();
650    //     path_.error   = 0;
651    //     cout <<"Loading p.f.a. parameters: "<<path<<endl;
652    //     return readetaparam_();
653    // }
654    
655    // /**
656    //  * Load magnetic field parameters (call the F77 routine).
657    //  *
658    //  */
659    // int TrkLevel1::LoadFieldParam(TString path){
660            
661    // //    if( strcmp(path_.path,path.Data()) ){
662    //     if( path.IsNull() ){
663    //      path = gSystem->Getenv("PAM_CALIB");
664    //      if(path.IsNull()){
665    //          cout << " TrkLevel1::LoadFieldParam() ==> No PAMELA environment variables defined "<<endl;
666    //          return 0;
667    //      }
668    //      path.Append("/trk-param/field_param-0/");
669    //     }
670    //     cout <<"Loading magnetic field "<<path<<endl;
671    //     strcpy(path_.path,path.Data());
672    //     path_.pathlen = path.Length();
673    //     path_.error   = 0;
674    //     return readb_();
675    // //    }      
676    // //    return 0;
677    // }
678    // /**
679    //  * Load magnetic field parameters (call the F77 routine).
680    //  *
681    //  */
682    // int TrkLevel1::LoadChargeParam(TString path){
683            
684    // //    if( strcmp(path_.path,path.Data()) ){
685    //     if( path.IsNull() ){
686    //      path = gSystem->Getenv("PAM_CALIB");
687    //      if(path.IsNull()){
688    //          cout << " TrkLevel1::LoadChargeParam() ==> No PAMELA environment variables defined "<<endl;
689    //          return 0;
690    //      }
691    //      path.Append("/trk-param/charge_param-1/");
692    //     }
693    //     cout <<"Loading charge-correlation parameters: "<<path<<endl;
694    //     strcpy(path_.path,path.Data());
695    //     path_.pathlen = path.Length();
696    //     path_.error   = 0;
697    //     return readchargeparam_();
698    // //    }      
699    // //    return 0;
700    // }
701    // /**
702    //  * Load magnetic field parameters (call the F77 routine).
703    //  *
704    //  */
705    // int TrkLevel1::LoadAlignmentParam(TString path){
706            
707    // //    if( strcmp(path_.path,path.Data()) ){
708    //     if( path.IsNull() ){
709    //      path = gSystem->Getenv("PAM_CALIB");
710    //      if(path.IsNull()){
711    //          cout << " TrkLevel1::LoadAlignmentParam() ==> No PAMELA environment variables defined "<<endl;
712    //          return 0;
713    //      }
714    //      path.Append("/trk-param/align_param-0/");
715    //     }
716    //     cout <<"Loading alignment parameters: "<<path<<endl;
717    //     strcpy(path_.path,path.Data());
718    //     path_.pathlen = path.Length();
719    //     path_.error   = 0;
720    //     return readalignparam_();
721    // //    }      
722    // //    return 0;
723    // }
724    // /**
725    //  * Load magnetic field parameters (call the F77 routine).
726    //  *
727    //  */
728    // int TrkLevel1::LoadMipParam(TString path){
729            
730    // //    if( strcmp(path_.path,path.Data()) ){
731    //     if( path.IsNull() ){
732    //      path = gSystem->Getenv("PAM_CALIB");
733    //      if(path.IsNull()){
734    //          cout << " TrkLevel1::LoadMipParam() ==> No PAMELA environment variables defined "<<endl;
735    //          return 0;
736    //      }
737    //      path.Append("/trk-param/mip_param-0/");
738    //     }
739    //     cout <<"Loading ADC-to-MIP conversion parameters: "<<path<<endl;
740    //     strcpy(path_.path,path.Data());
741    //     path_.pathlen = path.Length();
742    //     path_.error   = 0;
743    //     return readmipparam_();
744    // //    }      
745    // //    return 0;
746    // }
747    // /**
748    //  * Load magnetic field parameters (call the F77 routine).
749    //  *
750    //  */
751    // int TrkLevel1::LoadVKMaskParam(TString path){
752            
753    // //    if( strcmp(path_.path,path.Data()) ){
754    //     if( path.IsNull() ){
755    //      path = gSystem->Getenv("PAM_CALIB");
756    //      if(path.IsNull()){
757    //          cout << " TrkLevel1::LoadVKMaskParam() ==> No PAMELA environment variables defined "<<endl;
758    //          return 0;
759    //      }
760    //      path.Append("/trk-param/mask_param-1/");
761    //     }
762    //     cout <<"Loading VK-mask parameters: "<<path<<endl;
763    //     strcpy(path_.path,path.Data());
764    //     path_.pathlen = path.Length();
765    //     path_.error   = 0;
766    //     return readvkmask_();
767    // //    }      
768    // //    return 0;
769    // }
770    
771    // /**
772    //  * Load all (default) parameters. Environment variable must be defined.
773    //  *
774    //  */
775    // int TrkLevel1::LoadParams(){
776    
777    //     int result=0;
778        
779    //     result = result * LoadFieldParam();
780    //     result = result * LoadPfaParam();
781    //     result = result * LoadChargeParam();
782    //     result = result * LoadAlignmentParam();
783    //     result = result * LoadMipParam();
784    //     result = result * LoadVKMaskParam();
785    
786    //     return result;
787    // }
788    
789    
790    
791    int TrkLevel1::GetPfaNbinsAngle(){
792        TrkParams::Load(4);
793        if( !TrkParams::IsLoaded(4) ){
794            cout << "int TrkLevel1::GetPfaNbinsAngle() --- ERROR --- p.f.a. parameters  not loaded"<<endl;
795            return 0;
796        }
797        return pfa_.nangbin;
798    };
799    
800    int TrkLevel1::GetPfaNbinsETA(){
801        TrkParams::Load(4);
802        if( !TrkParams::IsLoaded(4) ){
803            cout << "int TrkLevel1::GetPfaNbinsETA() --- ERROR --- p.f.a. parameters  not loaded"<<endl;
804            return 0;
805        }
806        return pfa_.netaval;
807    };
808    
809  /**  /**
810   *   *
811   *   *
812   */   */
813  float* TrkLevel1::GetPfaCoord(TString pfa, int nview, int nladder, int nang){  float* TrkLevel1::GetPfaCoord(TString pfa, int nview, int nladder, int nang){
814    
815        TrkParams::Load(4);
816        if( !TrkParams::IsLoaded(4) ){
817            cout << "float* TrkLevel1::GetPfaCoord(TString pfa, int nview, int nladder, int nang) --- ERROR --- p.f.a. parameters  not loaded"<<endl;
818            return 0;
819        }
820        
821      int nbins = GetPfaNbinsETA();      int nbins = GetPfaNbinsETA();
822      if(!nbins)return 0;      if(!nbins)return 0;
# Line 655  float* TrkLevel1::GetPfaCoord(TString pf Line 843  float* TrkLevel1::GetPfaCoord(TString pf
843    
844  float* TrkLevel1::GetPfaAbs(TString pfa, int nang){  float* TrkLevel1::GetPfaAbs(TString pfa, int nang){
845        
846        TrkParams::Load(4);
847        if( !TrkParams::IsLoaded(4) ){
848            cout << "float* TrkLevel1::GetPfaAbs(TString pfa, int nang) --- ERROR --- p.f.a. parameters  not loaded"<<endl;
849            return 0;
850        }
851    
852      int nbins = GetPfaNbinsETA();      int nbins = GetPfaNbinsETA();
853      if(!nbins)return 0;      if(!nbins)return 0;
854    
# Line 673  float* TrkLevel1::GetPfaAbs(TString pfa, Line 867  float* TrkLevel1::GetPfaAbs(TString pfa,
867    
868      return fcorr;      return fcorr;
869    
   
870  };  };
871    
872    /**
873     * Method to call the F77 routine that performs level1->level2 processing.
874     * The level2 output is stored in a common block, which can be retrieved
875     * by mean of the method TrkLevel2::SetFromLevel2Struct().
876     * @param pfa Position finding algorythm used to reconstruct the track
877     * Implemented algorythms:
878     * 0  ETA (default)
879     * 1  ---
880     * 2  ETA2
881     * 3  ETA3
882     * 4  ETA4
883     * 10 COG
884     * 11 COG1
885     * 12 COG2
886     * 13 COG3
887     * 14 COG4
888     * NB If the TrkLevel1 object is readout from a tree, and the
889     * TrkLevel1::ProcessEvent(int pfa) is used to reprocess the event, attention
890     * should be payed to the fact that single clusters (clusters not associated
891     * with any track) might not be stored. Full reprocessing should be done from
892     * level0 data.
893     */
894    int TrkLevel1::ProcessEvent(int pfa){
895    
896    //    cout << "int TrkLevel1::ProcessEvent()" << endl;
897        TrkParams::Load( );
898        if( !TrkParams::IsLoaded() )return 0;
899    
900        GetLevel1Struct();
901    
902        analysisflight_(&pfa);
903    
904        return 1;
905    
906    }
907    
908    
909  ClassImp(TrkLevel1);  ClassImp(TrkLevel1);
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.23