TrigT1TGC : Implemented L1 muon endcap thresholds for run -3 (ATR-19291)

Trigger/TrigT1/TrigT1TGC/TrigT1TGC/TGCRPhiCoincidenceMap.hh

@@ -19,8 +19,12 @@ namespace LVL1TGCTrigger {
 class TGCRPhiCoincidenceMap {
 public:
   bool test(int octantId, int moduleId, int subsector, 
-	    int type, int pt, 
-	    int dr, int dphi) const;
+            int type, int pt, 
+            int dr, int dphi) const;// for Run2
+
+  int test_Run3(int octantId, int moduleId, int subsector, 
+                int type, int dr, int dphi) const; //for Run3
+
 
   int   getMapType(int hlwire, int hlstrip) const;
   const std::string& getVersion() const;
@@ -30,9 +34,9 @@ public:
   void  setFullCW( bool val);
 
   TGCRPhiCoincidenceMap(TGCArguments*,
-			const SG::ReadCondHandleKey<TGCTriggerData>& readCondKey,
+                        const SG::ReadCondHandleKey<TGCTriggerData>& readCondKey,
                         const std::string& version,
-			int   sideId=0, int octantId=0);
+                        int   sideId=0, int octantId=0);
 
   virtual ~TGCRPhiCoincidenceMap();
 
@@ -40,7 +44,8 @@ public:
   TGCRPhiCoincidenceMap(const TGCRPhiCoincidenceMap& right);
   TGCRPhiCoincidenceMap& operator=(const TGCRPhiCoincidenceMap& right);
 
-  bool readMap();  
+  bool readMap(); // for Run2
+  bool readMap_Run3();// for Run3
 
   TGCArguments* tgcArgs() const;
 
@@ -58,10 +63,21 @@ protected:
  
   enum {TMap_HH=0, TMap_HL, TMap_LH, TMap_LL, N_TMap};
   enum {N_PT_THRESH=6};
+  enum {NumberOfCoincidenceType=4};
+  enum {NumberOfModuleType=12};
   enum {DR_offset=-15, DPhi_offset=-7};
 
+ //converter from char to pTthre Number. the sign mean muon sign.
+  std::map<char, int> m_pTdef={{'X',0},
+              {'A',1},{'B',2},{'C',3},{'D',4},{'E',5},{'F',6},{'G',7},{'H',8},{'I',9},{'J',10},{'K',11},{'L',12},{'M',13},{'N',14},{'O',15},
+              {'a',-1},{'b',-2},{'c',-3},{'d',-4},{'e',-5},{'f',-6},{'g',-7},{'h',-8},{'i',-9},{'j',-10},{'k',-11},{'l',-12},{'m',-13},{'n',-14},{'o',-15} };
+
+
 private:
-  std::map<int, std::map<int,int> > m_mapDB[N_PT_THRESH];
+  std::map<int, std::map<int,int> > m_mapDB[N_PT_THRESH];//for Run2 [ptLevel]< RoI(&type),<RNumber,RWindow> >
+  std::map<int, std::map<int, std::map<int, char> > > m_mapDB_Run3;//for Run3 <RoI(&type),<R,<Phi,pT(char)> > >
+
+
 
   int m_numberOfDR,m_numberOfDPhi;
   std::string m_verName;

Trigger/TrigT1/TrigT1TGC/TrigT1TGC/TGCRPhiCoincidenceOut.hh

@@ -28,11 +28,13 @@ public:
   bool getInnerCoincidenceFlag() const {return m_innerCoincidenceFlag;};
 
   bool getHit(int ptIn) const { return m_hit[ptIn];};
+  int  getpT() const {return m_pT;}; // for Run3 Algo
   void setHit(int ptIn){ m_hit[ptIn]=true;};
   void clearHit(int ptIn) { m_hit[ptIn] = false;};
-  int  getPtLevel() const;
+  int  getPtLevel() const; // for Run2 Algo
 
   void setIdSSC(int idSSCIn){ m_idSSC = idSSCIn;};
+  void setpT(int pTIn){ m_pT=pTIn;};
   void setR(int rIn){ m_r=rIn;};
   void setPhi(int phiIn){ m_phi=phiIn;};
   void setDR(int drIn) { m_dR = drIn; };
@@ -56,7 +58,8 @@ public:
 
 private:
   int m_idSSC;
-  bool m_hit[NumberOfPtLevel+1];
+  bool m_hit[NumberOfPtLevel+1]; // for Run2 Algo
+  int m_pT; // for Run3 Algo
   int m_phi;
   int m_r;
   int m_dR;

Trigger/TrigT1/TrigT1TGC/src/TGCRPhiCoincidenceMap.cxx

@@ -60,6 +60,40 @@ bool TGCRPhiCoincidenceMap::test(int octantId, int moduleId, int subsector,
   else return false;
 }
 
+
+
+
+int TGCRPhiCoincidenceMap::test_Run3(int octantId, int moduleId, int subsector, 
+                                     int type, int dr, int dphi) const
+{  
+  // check pt range
+  if (type<0 || type>=N_TMap ) return false; 
+
+  int sector=(moduleId-2)/3+octantId*3;
+  int phimod2 = (moduleId==2||moduleId==5||moduleId==8)&&(sector%2==1) ? 1 : 0;
+  int addr=SUBSECTORADD(subsector, moduleId, phimod2,type);
+
+  std::map<int, std::map<int, std::map<int, char> > >::const_iterator it=m_mapDB_Run3.find(addr);
+  if (it==m_mapDB_Run3.end()) return false;
+
+  std::map<int, std::map<int, char> > mapR = it->second;
+  std::map<int, std::map<int, char> >::const_iterator itR=mapR.find(dr);
+  if (itR==mapR.end()) return false;
+
+  std::map<int, char> mapPhi = itR->second;
+  std::map<int, char>::const_iterator itPhi=mapPhi.find(dphi);
+  if (itPhi==mapPhi.end()) return false;
+
+  char pT_char=itPhi->second;
+  int pT_int=m_pTdef.find(pT_char)->second;
+
+  return  pT_int;
+}
+
+
+
+
+
 TGCRPhiCoincidenceMap::TGCRPhiCoincidenceMap(TGCArguments* tgcargs,
 					     const SG::ReadCondHandleKey<TGCTriggerData>& readCondKey,
                                              const std::string& version,
@@ -73,16 +107,15 @@ TGCRPhiCoincidenceMap::TGCRPhiCoincidenceMap(TGCArguments* tgcargs,
    m_readCondKey(readCondKey)
 {
   if (!tgcArgs()->USE_CONDDB()) {
-    if (!checkVersion()){
+    if(!tgcArgs()->useRun3Config()){
+      if (!checkVersion()){
       m_verName = "NA";
       return;
-    }
-
-    if(!tgcArgs()->useRun3Config()){
+      }
       this->readMap();  // read Coincidence Map for Run2 (6 thresholds)
     }
     else{
-      //this -> readMap_Run3(); This function will be implemented. // read Coincidence Map for Run3 (15 thresholds)
+      this -> readMap_Run3();// read Coincidence Map for Run3 (15 thresholds)
     }
 
   } 
@@ -170,7 +203,7 @@ TGCRPhiCoincidenceMap::TGCRPhiCoincidenceMap(const TGCRPhiCoincidenceMap& right)
   m_octant=right.m_octant;
   m_fullCW=right.m_fullCW;
   if(m_verName!="NA" && !tgcArgs()->useRun3Config()) this->readMap();  // read Coincidence Map for Run2 (6 thresholds)
-  if(m_verName!="NA" &&  tgcArgs()->useRun3Config()){/*this -> readMap_Run3(); This function will be implemented. */}// read Coincidence Map for Run3 (15 thresholds)
+  if(m_verName!="NA" &&  tgcArgs()->useRun3Config()){this -> readMap_Run3();}// read Coincidence Map for Run3 (15 thresholds)
 }
 
 
@@ -184,7 +217,7 @@ TGCRPhiCoincidenceMap& TGCRPhiCoincidenceMap::operator=(const TGCRPhiCoincidence
     m_octant=right.m_octant;
     m_fullCW=right.m_fullCW;
     if(m_verName!="NA" && !tgcArgs()->useRun3Config()) this->readMap();  // read Coincidence Map for Run2 (6 thresholds)
-    if(m_verName!="NA" &&  tgcArgs()->useRun3Config()){ /*this -> readMap_Run3(); This function will be implemented.*/} // read Coincidence Map for Run3 (15 thresholds)
+    if(m_verName!="NA" &&  tgcArgs()->useRun3Config()){ this -> readMap_Run3();} // read Coincidence Map for Run3 (15 thresholds)
    }
 
   return *this;
@@ -295,5 +328,122 @@ bool TGCRPhiCoincidenceMap::readMap()
 }
 
 
+bool TGCRPhiCoincidenceMap::readMap_Run3() 
+{
+  const int moduleNumber[NumberOfModuleType]  =
+    {  0,  1,   2,   2,  3,  4,   5,   5,  6,  7,   8,  8 };
+  const std::string moduleName[NumberOfModuleType]=
+    {"0a","1a","2a","2b","3a","4a","5a","5b","6a","7a","8a","8b"};
+  const std::string sideName[NumberOfSide] = {"a","c"};
+  const std::string octantName[NumberOfOctant] =
+    {  "0", "1", "2", "3", "4", "5", "6", "7"};
+  const std::string coincidenceTypeName[NumberOfCoincidenceType] = {"HH","HL","LH","LL"};
+
+
+  IMessageSvc* msgSvc = 0;
+  ISvcLocator* svcLocator = Gaudi::svcLocator();
+  if (svcLocator->service("MessageSvc", msgSvc) == StatusCode::FAILURE) {
+    return false;
+  }
+  MsgStream log(msgSvc, "TGCRPhiCoincidenceMap::TGCRPhiCoincidenceMap");
+
+
+  // initialize
+  std::string buf;
+  std::string fn, fullName, tag;
+  int ssId;
+  char delimiter = '\n';
+
+  // loop over all files...
+  for(int iModule=0; iModule<NumberOfModuleType; iModule+=1) {
+    for(int iCoinType=0; iCoinType!=NumberOfCoincidenceType; iCoinType++){
+      int phimod2=moduleName[iModule].find("b")!=std::string::npos ? 1 : 0;
+      std::string fn;
+      if (m_fullCW) {
+        fn = "/BW/cm_"+ sideName[m_side]+octantName[m_octant]+moduleName[iModule]+coincidenceTypeName[iCoinType]+"_v0001.db";}
+      else{
+        fn = "/BW/cm_"+ sideName[m_side]+moduleName[iModule]+coincidenceTypeName[iCoinType]+"_Octant_v0001.db";
+      }
+      bool Forward_type1=(moduleName[iModule]=="2b"||moduleName[iModule]=="5a"||moduleName[iModule]=="8b");
+      bool Forward_type2=(moduleName[iModule]=="2a"||moduleName[iModule]=="5b"||moduleName[iModule]=="8a");
+      if(m_octant%2==0 && Forward_type1){continue;}
+      if(m_octant%2==1 && Forward_type2){continue;}
+      int type = -1;
+      int lDR, hDR, lDPhi, hDPhi;
+
+      fullName = PathResolver::FindCalibDirectory("TrigT1TGC_CW")+fn;
+      if( fullName.length() == 0 ) { 
+        log << MSG::ERROR
+            << " Could not found "
+            << fn.c_str()
+            <<endmsg;
+        continue;  
+      } 
+
+      std::ifstream file(fullName.c_str(),std::ios::in);   
+      if(!file){
+        log << MSG::ERROR
+            << " Could not found "
+            << fullName.c_str()
+            <<endmsg;
+        continue; 
+      } 
+      while(getline(file,buf,delimiter)){
+        std::istringstream header(buf);
+        header>>tag;
+        if(tag=="#") { // read header part.     
+          header>>ssId>>lDR>>hDR>>lDPhi>>hDPhi;
+          type = getTYPE( lDR, hDR, lDPhi, hDPhi );
+          // check moduleNumber and ptLevel
+          if( type<0 ) {
+            log << MSG::WARNING
+                << " illegal parameter in database header : "
+                <<  header.str()
+                <<" in file "
+                << fn
+                <<endmsg;
+            break;
+          }
+
+          // get window data
+          std::map<int, std::map<int, char> >  bWindow;//<R,<~>>
+          char pT;
+          for(int ir=0; ir<=hDR-lDR; ir++) {
+            getline(file,buf,delimiter);
+            std::map<int, char> aWindow;//<Phi,pT>
+            for(int iphi=0; iphi<=hDPhi-lDPhi; iphi++){
+              pT = buf[iphi];
+              if (pT=='X'){continue;} // none of window is opened in this dR
+              aWindow[iphi+DPhi_offset] = pT;
+            }
+            // Warning : no window 
+            if (aWindow.size()==0) {
+              log << MSG::DEBUG
+                  << " No window is opened for (ptLevel,ssId,mod) = ("
+                  <<  ssId
+                  <<" ) "
+                  <<endmsg;
+            } else {
+              bWindow[ir+DR_offset]=aWindow;
+            }
+          }
+
+          int addr = SUBSECTORADD(ssId,moduleNumber[iModule],phimod2,type);
+          if (m_mapDB_Run3.find(addr)!=m_mapDB_Run3.end()) {
+            log << MSG::DEBUG
+                << " This subsector was already reserved."
+                <<endmsg;
+          } else {
+            m_mapDB_Run3[addr]=bWindow;
+          }
+        }
+      }
+      
+    }
+  }
+  return true;
+
+}
+
 
 } //end of namespace bracket

Trigger/TrigT1/TrigT1TGC/src/TGCRPhiCoincidenceMatrix.cxx

@@ -86,12 +86,12 @@ TGCRPhiCoincidenceOut* TGCRPhiCoincidenceMatrix::doCoincidence()
     // calculate pT of muon candidate
     if(tgcArgs()->useRun3Config()){
       //Algorithm for Run3
-      /*int pt=map->test_Run3(sectorLogic->getOctantID(),sectorLogic->getModuleID(),
-	subsector,type,dR,dPhi[j]); // this function will be implemented. 
-	ptOut = std::abs(pt)-1;
-	chargeOut = pt<0 ? 0:1;
-	//isgoodMFOut : will be set.
-      */
+      int pt=m_map->test_Run3(m_sectorLogic->getOctantID(),m_sectorLogic->getModuleID(),
+                            subsector,type,m_dR,m_dPhi[j]); // this function will be implemented. 
+      ptOut = std::abs(pt);
+      chargeOut = pt<0 ? 0:1;
+      //isgoodMFOut : will be set.
+      
       CoincidenceTypeOut=(type==0);
     }
     else{
@@ -110,7 +110,8 @@ TGCRPhiCoincidenceOut* TGCRPhiCoincidenceMatrix::doCoincidence()
       ptMax = ptOut;
       out->clear();    
       out->setIdSSC(m_SSCId);
-      out->setHit(ptMax+1);   
+      if(!tgcArgs()->useRun3Config()){out->setHit(ptMax+1);}// for Run2 Algo
+      else{out->setpT(ptMax);}// for Run3 Algo
       out->setR(m_r);
       out->setPhi(m_phi[j]);
       out->setDR(m_dR);

Trigger/TrigT1/TrigT1TGC/src/TGCRPhiCoincidenceOut.cxx

@@ -12,6 +12,7 @@ namespace LVL1TGCTrigger {
 
 TGCRPhiCoincidenceOut::TGCRPhiCoincidenceOut():
   m_idSSC(-1),
+  m_pT(0),
   m_phi(-1),
   m_r(-1),
   m_dR(0),

Trigger/TrigT1/TrigT1TGC/src/TGCSectorLogic.cxx

@@ -187,7 +187,10 @@ void TGCSectorLogic::clockIn(const SG::ReadCondHandleKey<TGCTriggerData> readCon
     }
     ////////////////////////////////////////////
     // do coincidence with Inner Tracklet of EIFI and/or TileMu
-    if (m_useEIFI) doInnerCoincidence(readCondKey, SSCid, coincidenceOut);
+    if (m_useEIFI){ 
+      if(!tgcArgs()->useRun3Config()){doInnerCoincidence(readCondKey, SSCid, coincidenceOut);}
+      else{/*InnerCoincidence Algorithm for Run3 will be implemented;*/}
+    }
 
     if(coincidenceOut){
       if(tgcArgs()->useRun3Config()){m_trackSelector.input(coincidenceOut);}// TrackSelector for Run3

Trigger/TrigT1/TrigT1TGC/src/TGCTrackSelector.cxx

@@ -26,7 +26,7 @@ namespace LVL1TGCTrigger {
 
   void TGCTrackSelector::input(TGCRPhiCoincidenceOut* rPhiOut){
     if(rPhiOut!=0){
-      if(rPhiOut->hasHit()){
+      if(rPhiOut->getpT()!=0){
         m_coincidenceIn[m_numberOfCandidate_In].reset(rPhiOut);
         m_numberOfCandidate_In++;
       }
@@ -46,7 +46,7 @@ namespace LVL1TGCTrigger {
       for(int track2=track1+1;track2!=m_numberOfCandidate_In;track2++){
         bool compare_result=compare(m_coincidenceIn[track1].get(),m_coincidenceIn[track2].get()); //1>2 true:1 1<2 false:0
 
-	//TrackPriorityRank : Smaller value is high priority.
+       //TrackPriorityRank : Smaller value is high priority.
         m_trackPriorityRank[track1]+=(int)!compare_result;
         m_trackPriorityRank[track2]+=(int)compare_result;
       }
@@ -61,7 +61,7 @@ namespace LVL1TGCTrigger {
         int R= 2*m_coincidenceIn[track]->getIdSSC()+m_coincidenceIn[track]->getR() - (m_sectorLogic->getRegion()==Endcap ? 1 : 0);
         TrackcandidateOut->setR(m_numberOfCandidate_Out,R);
         TrackcandidateOut->setPhi(m_numberOfCandidate_Out,m_coincidenceIn[track]->getPhi());
-        TrackcandidateOut->setPtLevel(m_numberOfCandidate_Out,m_coincidenceIn[track]->getPtLevel());
+        TrackcandidateOut->setPtLevel(m_numberOfCandidate_Out,m_coincidenceIn[track]->getpT());
         TrackcandidateOut->setDR(m_numberOfCandidate_Out,2*m_coincidenceIn[track]->getDR());
         TrackcandidateOut->setDPhi(m_numberOfCandidate_Out,m_coincidenceIn[track]->getDPhi());
         TrackcandidateOut->setInnerVeto(m_numberOfCandidate_Out,m_coincidenceIn[track]->getInnerVeto());
@@ -85,8 +85,8 @@ namespace LVL1TGCTrigger {
 
 
     // The definition of priority is not fixed. This function wiil be updated.
-    if(track1->getPtLevel()>track2->getPtLevel()){return true;}
-    if(track1->getPtLevel()<track2->getPtLevel()){return false;}
+    if(track1->getpT() > track2->getpT()){return true;}
+    if(track1->getpT() < track2->getpT()){return false;}
     else{
       // when the pt level of 2 tracks are same, track with large R is selected.
       if(track1->getIdSSC()<track2->getIdSSC()){return true;}